Estimation of S-wave Velocity Structures by Using Microtremor Array ...

4 downloads 0 Views 2MB Size Report
Abstract. Jakarta is located on a thick sedimentary layer that potentially has a very high seismic wave amplification. However, the available information.
J. Math. Fund. Sci., Vol. 46, No. 3, 2014, 313-327

313

Estimation of S-wave Velocity Structures by Using Microtremor Array Measurements for Subsurface Modeling in Jakarta Mohamad Ridwan1*, Afnimar2, Sri Widiyantoro2,3, Masyhur Irsyam3,4 & Hiroaki Yamanaka5 1

Graduate Research in Earthquakes and Active Tectonics, Faculty of Earth Science and Technology, Institut Teknologi Bandung 2 Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung 3 Research Center for Disaster Mitigation, Institut Teknologi Bandung 4 Faculty of Civil and Environmental Engineering, Institut Teknologi Bandung Jalan Ganesa No. 10, Bandung 40132, Indonesia 5 Department of Environmental Science and Technology, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8550, Japan *Email: [email protected] Abstract. Jakarta is located on a thick sedimentary layer that potentially has a very high seismic wave amplification. However, the available information concerning the subsurface model and bedrock depth is insufficient for a seismic hazard analysis. In this study, a microtremor array method was applied to estimate the geometry and S-wave velocity of the sedimentary layer. The spatial autocorrelation (SPAC) method was applied to estimate the dispersion curve, while the S-wave velocity was estimated using a genetic algorithm approach. The analysis of the 1D and 2D S-wave velocity profiles shows that along a north-south line, the sedimentary layer is thicker towards the north. It has a positive correlation with a geological cross section derived from a borehole down to a depth of about 300 m. The SPT data from the BMKG site were used to verify the 1D S-wave velocity profile. They show a good agreement. The microtremor analysis reached the engineering bedrock in a range from 359 to 608 m as depicted by a cross section in the north-south direction. The site class was also estimated at each site, based on the average S-wave velocity until 30 m depth. The sites UI to ISTN belong to class D (medium soil), while BMKG and ANCL belong to class E (soft soil). Keywords: engineering bedrock; microtremor array; S-wave velocity; site class; Jakarta.

1

Introduction

Jakarta is one of the major cities of Indonesia for which the government has prioritized the preparation of seismic microzonation maps. This was done under consideration of various aspects, such as seismicity, infrastructure, population, etc., in which Jakarta has a high risk. Site response analysis is part of the Received December 20th, 2013, Revised October 9th, 2014, Accepted for publication November 3rd, 2014. Copyright © 2014 Published by ITB Journal Publisher, ISSN: 2337-5760, DOI: 10.5614/j.math.fund.sci.2014.46.3.9

314

Mohamad Ridwan, et al.

microzonation studies that consider the effects of ground shaking. This is to be calculated by using the results of a complete geotechnical investigation. Subsurface information, such as bedrock depths and soil dynamic properties, are important parameters in the dynamic response analysis of the ground related to a hazard analysis for the surface. Previous studies concerning the subsurface geology of Jakarta have been conducted by Padmosukismo and Yahya [1] and showed the configuration of the Northwest Java basin and the Ciputat sub-basin in which Jakarta is located. Turkandi, et al. [2] have conducted a regional geological mapping and Fachri, et al. [3] studied the hydrostratigraphy of the groundwater basins. The results of these studies provide very important information, but further studies for earthquake engineering purposes are still required to obtain more physical parameters and the position of bedrock. In this study, a microtremor array method was applied to investigate the subsurface underneath Jakarta and to construct a stratigraphic model based on S-wave velocity parameters. The advantage of this method is that it simplifies the field operation because it does not require active vibration sources or boreholes. Microtremor array observation was first developed by Aki [4] to obtain the dispersion curve by using a circular array employing the SPAC method. Then Okada [5 and 6] expanded this method for subsurface exploration and modified the array configuration into a triangular shape using only 4 sensors. Furthermore, Morikawa, et al. [7] developed 2sSPAC by employing only 2 sensors in an effort to reduce difficulties in the field work. Microtremor array measurements were conducted in Jakarta at 10 sites along a north-south line, from ANCL to UI (Figure 1), in order to determine bedrock depth and site class. To achieve the objective above, we used small and large arrays for shallow and deep surveys, respectively. Borehole data to a depth of 150 m at the BMKG site were available, so they could be used to verify the 1D S-wave velocity model derived from the analysis of the microtremor data along with the standard penetration test (SPT) data.

2

Geological Setting

Jakarta is located in the Ciputat sub-basin of the Northwest Java basin. The Ciputat sub-basin is filled by Tertiary and Quarter sediments and controlled by geological structures with a dominant north-south direction that form heights and depression zones [1]. Referring to the regional geological maps [2 and 3], the stratigraphy sequence from old to young is as follows: first, the Rengganis formation (Early Miocene),

Estimation of S-wave Velocity Structures in Jakarta

315

then the Bojongmanik formation in the western part, which is inter-fingered with the Jatiluhur formations in the eastern part lying on top of it (Middle Miocene), while above the two formations there are the Klapanunggal formation (Late Miocene), the Genteng formation (Early Pliocene), and the Serpong formation (Late Pliocene). The Quarter sedimentary rocks covering the whole Jakarta area are: Banten tuff (Plio-Pleistocene), young volcanic rocks (Pleistocene), alluvium fan, coastal sediment and alluvium.

3

Data and Method

3.1

Data Acquisition

6°6'S

Microtremor array measurements were performed at 10 sites in Jakarta along the ANCL-UI line (Figure 1) using triangular arrays with 4 sensors, where one sensor was placed at the center and 3 sensors were deployed on a circular array. (See Figures 2a). Figure 2b shows an example of the array configuration for microtremor measurements carried out at the BMKG site.

6°8'S

A # ANCL BMKG

6°10'S

ISTN

6°12'S

#

6°16'S

6°14'S

#MRIA # KNGN # DRN3 CPT

#

6°18'S 6°20'S 6°22'S

#!

#KBGS Legend

#ZPUR

!

Deep Borehole

#

Microtremors Array Observation Sites Jakarta City

106°42'E

106°44'E

106°46'E

-

A’ #UI

106°48'E

106°50'E

106°52'E

106°54'E

106°56'E

106°58'E

107°0'E

Figure 1 Locations of microtremor array measurements along the UI-ANCL line and deep borehole.

316

Mohamad Ridwan, et al.

1

4 2

3

a)

b)

Source= 0.0km

Channel number

0

10

Time (s) 20

30

40

50

60

70

80

90

100

110

120

130

140

110030

150

160

170 0.70801V

110032

0.10577V

110035

0.20287V

110038

0.15001V

E:\MT_Array Jak-14\Pd.Ranggon\McSEISMT_NEO_0508200420_0016.sg2 Source= 0.0km Time (s)

Channel number

0

4

8

12

110030

16

20 0.70801V

110032

0.10577V

110035

0.20287V

110038

0.15001V

E:\MT_Array Jak-14\Pd.Ranggon\McSEISMT_NEO_0508200420_0016.sg2

c) Figure 2 a) Array shape and location of 4 sensors, b) example of microtremor array measurements at the BMKG site, and c) example of microtremor data record.

The array size was designed in consideration of the target layer depth. Array sizes of 600 m, 300 m and 150 m were used for deep layer surveys, while for shallow layers the array sizes were 30 m and 15 m. Microtremor data recording was conducted in 1-1.5 hours for the large arrays, while for the small arrays it was approximately 30-45 minutes [5 and 6]. The data records were digitized with a sampling interval of 0.01 seconds and then divided into several blocks consisting of 16,384 samples for each block to be used in the data processing. Figure 2c shows an example of the data recorded at the BMKG site for a time block with a duration of 163.84 seconds. The instruments used in the field work were: four units of microtremor equipment (McSEIS Neo-MT) produced by OYO Corporation completed with GPS for time synchronization. For the data processing, SeisImager and Yamanaka software was used.

Estimation of S-wave Velocity Structures in Jakarta

3.2

317

SPAC Method

Following, Okada [6] and Morikawa, et al. [7], phase velocity calculations were performed using the SPAC method for the vertical component of microtremors. For the triangular array configuration, the coherence function between two stations within a distance r was calculated. SPAC coefficient ρ(ω,r) is the average of the coherence function for all pairs of stations on the circular array, which can be directly calculated by the following equation:

 (; r ) 

1 2

Re  SCX (; r , )

2

 S 0

C

(;0,0)   S X (; r , )

d

(1)

where Sc (ω,0,0) and Sx (ω,r,θ) are the power spectrum of the microtremors at stations C and X, while Scx (ω,r,θ) is the cross spectrum between stations C and X. Since the array configuration is circular, the SPAC coefficient can be expressed as a Bessel function with variables r and k, as follows: ρ (ω, r) = J0(r,k), where k = ω/c(ω) ; c(ω) = phase velocity, therefore:  r    c( ) 

(2)

 (; r )  J 0 

Jo = Bessel function of the first kind zero order, c(ω) = phase velocity as a function of frequency ω, r = distance between two stations.

3.3

Inversion of Dispersion Curves Using Genetic Algorithm

The S-wave velocity structures were estimated from the dispersion curves through an inversion technique using a genetic algorithm. This technique was introduced in seismology by Yamanaka and Ishida [8] as a global optimizing method. The inversion process is conducted on the searched area as the population by minimizing the misfit function (ϕ), which can be defined by the difference between observed phase velocity Uo(Ti) and calculated phase velocity Uc(Ti,), and can be expressed as follows: ∑

(

( ) ( )

( )

)

(3)

where N is the number of observed data and σ(Ti) is the standard deviation of the observed phase velocity over a period of Ti.

318

Mohamad Ridwan, et al.

A genetic algorithm was applied to the phase velocities of the fundamental modes of the Rayleigh waves, assuming the layer models are horizontal, isotropic, and homogenous. In the inversion, the model parameters of S-wave velocities and thicknesses are the unknown parameters to be estimated, while the density is fixed. In this study, the inversion was conducted by performing 20 populations with 10 iterations of the genetic algorithm and the parameters of genetic operation following Yamanaka and Ishida [8]. The inversion was conducted on the initial population by all genetic operations, where the model that had the smallest misfit survived in the next generation, while bad models were replaced.

4

Results

4.1

Dispersion Curves

The results of the estimated dispersion curves using the SPAC method for deep layers is shown in Figure 3a. The dispersion curves of microtremors were obtained in a frequency range of 0.2-5.0 Hz and have various phase velocities. Figure 3b shows the dispersion curves for the shallow layers, which have a higher frequency range (3-16 Hz). The variation of phase velocities is due to local geological conditions. 500

1200

BMKG

800

ISTN MRIA KNGN

600

DRN3 CPT

400

KBGS ZPUR

200

ANCL

400

ANCL

Phase Velocity (m/s)

Phase Velocity (m/s)

Dispersion Curve

Dispersion Curve

1000

BMKG ISTN

300

MRIA KNGN

DRN3

200

CPT KBGS ZPUR

100

UI

a)

UI

b)

0

0 0

1

2

3

Frequency (Hz)

4

5

0

5

10

15

20

25

Frequency (Hz)

Figure 3 Dispersion curves of microtremors for: a) deep layers, and b) shallow layers.

The inversion technique was performed for specified initial models, considering geological conditions. This study assumed 5 layered models for deep layers and 4 layered models for shallow layers. These initial models were

Estimation of S-wave Velocity Structures in Jakarta

319

inverted by using a genetic algorithm to generate close fits of dispersion curves between observed and calculated ones, as shown in Figures 4 and 5.

0 0 1000

1

2

3

4

800 600 400 200 0

1

2

3

0

1

2

3

4

400 200

0

Frequency (Hz)

1

2

3

1

2

3

600 400

4

200 0 5

600 400 200

4

0

Frequency (Hz)

1

2

3

0

5 1000

1

2

3

600 400 200

4

600 400

Obs 200 Cal 0 4 5 0

5

Frequency (Hz)

2

3

0

1000

1

2

3

4

5

Frequency (Hz) ANCL

800 600 400 200

Obs Cal

1

Obs Cal

0

Frequency (Hz) ZPUR

800

KNGN

800

Obs Cal

0

Frequency (Hz) KBGS

800

Obs Cal

0

0

MRIA

800

Obs Cal

1000

1000

Phase Velocity (m/s)

200

5

Frequency (Hz) CPT

600

5

400

0

800

4

600

Obs Cal

0 5

Obs Cal

0

200

1000

Frequency (Hz) DRN3

Phase Velocity (m/s)

Phase Velocity (m/s)

400

Obs Cal

Phase Velocity (m/s)

200

600

ISTN

800

Phase Velocity (m/s)

400

1000

BMKG

800

Phase Velocity (m/s)

600

1000

Phase Velocity (m/s)

Phase Velocity (m/s)

Phase Velocity (m/s)

ANCL

800

Phase Velocity (m/s)

1000

1000

Obs Cal

0

4

5

0

Frequency (Hz)

1

2

3

4

5

Frequency (Hz)

Figure 4 Comparison between observed (circles) and calculated (solid lines) dispersion curves for deep layers.

0

3

6

9

12

Phase Velocity (m/s)

400 300 200 100 0 0

3

6

9

12

Frequency (Hz)

0

3

6

9

12

400 300 200 100

15

100

0 0

3

6

9

12

Frequency (Hz)

0 500

Obs Cal 15

300 200 100

Obs Cal

0

15

Frequency (Hz) CPT

500

Obs Cal

200

Obs Cal

0

15

Frequency (Hz) DRN3

500

100

300

3

6

9

12

15

300 200 100

0

3

6

9

12

Frequency (Hz)

3

6

9

300 200 100

15

300 200

12

0

3

6

9

12

Frequency (Hz)

Obs Cal

0 0

3

6

9

12

15

Frequency (Hz) U I

500 400 300 200

Obs Cal

0

100

15

Frequency (Hz) ZPUR

400

Obs Cal

0

0

500

KNGN

400

Obs Cal

0

Frequency (Hz) KBGS

400

Phase Velocity (m/s)

Obs Cal

0

Phase Velocity (m/s)

200

500

MRIA

400

Phase Velocity (m/s)

100

300

Phase Velocity (m/s)

200

500

ISTN

400

Phase Velocity (m/s)

300

500

BMKG

400

Phase Velocity (m/s)

ANCL

400

Phase Velocity (m/s)

500

Phase Velocity (m/s)

Phase Velocity (m/s)

500

100

15

Obs Cal

0 0

3

6

9

12

15

Frequency (Hz)

Figure 5 Comparison between observed (circles) and calculated (solid lines) dispersion curves for shallow layers.

4.2

One-dimensional S-Wave Velocity Profiles

The 1D S-wave velocity profiles resulting from the individual inversion of each dispersion curve are illustrated in Figure 6 for deep layers and in Figure 7 for shallow depths. The bedrock depth with S-wave velocity value > 750 ms-1 [9] can be estimated from each 1D S-wave velocity profile. Generally, the bedrock depth increases toward northern Jakarta. In southern Jakarta (UI) the bedrock depth is about 378 m for Vs = 875 ms-1 and in the north (BMKG) the bedrock depth is about 651 m for Vs = 965 ms-1.

Mohamad Ridwan, et al.

S-wave Velocity (m/s)

750 1000 1250

0

250

500

-750

-500

-750

-1000

-1000

0

250

500

-500

250

500

-500

-1000

S-wave Velocity (m/s) 250

500

S-wave Velocity (m/s)

750 1000 1250

0

0

250

500

750 1000 1250

0

KBGS

-500

750 1000 1250

-750

0

ZPUR

-250

Depth (m)

Depth (m)

-500

500

KNGN

-500

750 1000 1250

CPT

250

-250

-1000

0

-250

0

-750

0

DRN3

S-wave Velocity (m/s)

750 1000 1250 0

S-wave Velocity (m/s)

750 1000 1250

0

-250

500

MRIA

-1000

0

250

-250

S-wave Velocity (m/s)

750 1000 1250

0 0

-750

S-wave Velocity (m/s) 500

S-wave Velocity (m/s)

750 1000 1250

-250

Depth (m)

Depth (m)

-500

250

500

ISTN

-250

0

250

BMKG

-250

Depth (m)

0 0

ANCL

Depth (m)

S-wave Velocity (m/s)

750 1000 1250

0

Depth (m)

500

UI

-250

Depth (m)

250

Depth (m)

S-wave Velocity (m/s) 0 0

-500

-250

Depth (m)

320

-500

-500

-750

-750

-750

-750

-750

-1000

-1000

-1000

-1000

-1000

Figure 6 1D S-wave velocity profiles at the 10 surveyed sites for deep layers.

The estimated 1D S-wave velocity profiles for shallow depths were also modeled into four layers down to 40 m depth. The site class determination was based on Vs30, where sites UI to ISTN belong to class D (medium soil), while BMKG and ANCL belong to class E (soft soil). 400

S-wave Velocity (m/s)

500

0

100

200

300

-20

400

0

100

200

300

400

100

200

300

400

500

-20

Vs30 = 268 ms-1

S-wave Velocity (m/s)

500

0

100

200

300

400

S-wave Velocity (m/s)

500

0

100

200

300

400

0

ZPUR

UI -10

-20

-30

Vs30 = 235 ms-1 -40

300

-40

-10

-30

200

-30

0

-20

Vs30 = 300 ms-1 -40

400

100

KNGN

-40

0

0

Vs30 = 250 ms-1

-10

-30

Vs30 = 200 ms-1

S-wave Velocity (m/s)

500

-20

KBGS

Depth (m)

Depth (m)

-30

400

-30

0

-20

300

-10

S-wave Velocity (m/s)

500

-10

-20

200

MRIA

-20

CPT

-10

100

0

Vs30 = 227 ms-1

0

DRN3

0

-10

S-wave Velocity (m/s)

500

S-wave Velocity (m/s)

500

-40

S-wave Velocity (m/s) 0

400

-30

-40

300

300

0

Vs30 = 172 ms-1

-40

200

200

ISTN

-30

Vs30 = 156 ms-1

100

100

-10

Depth (m)

Depth (m)

Depth (m)

-20

0

0

BMKG -10

-30

Depth (m)

S-wave Velocity (m/s)

500 0

ANCL -10

-40

400

0

Depth (m)

300

Depth (m)

200

Depth (m)

100

Depth (m)

S-wave Velocity (m/s) 0 0

-20

-30

Vs30 = 270 ms-1

Vs30 = 282 ms-1 -40

Figure 7 Same as Figure 6, but for shallow depths.

-40

500

Estimation of S-wave Velocity Structures in Jakarta

4.3

321

Comparison with Borehole Data

S-wave velocity parameters can also be calculated from standard penetration test (SPT) data using an empirical formula, as shown in Table 1. All these empirical formulas can be used for Vs calculation by taking the average of the results. At the BMKG site, SPT was conducted down to a depth of 100 m in a borehole near the microtremor measurement site. Comparison of the S-wave velocity profiles derived from SPT and microtremor array measurements demonstrates a strong similarity (Figure 8) and has a positive correlation with the lithological condition. Empirical correlation between N-SPT and S-wave velocity [10]. Vs (m/sec) Correlation Vs = 96.9 N0.314 Vs = 91 N0.337 Vs = 85.3 N0.341 Vs = 101 N0.29

References Imai and Tonouchi [11] Imai [12] Ohta and Goto [13] Sykora and Stokoe [14] Shear Wave Velocity (m/s) 0

250

500

750

0

BMKG

400

Shear Wave Velocity (m/s) 600

200

40

300

60

500 600

80

120 140

800

160

900

180

1000

200

From Microtremor

a)

0

From SPT

400

600

BMKG 10

20

30

_ _ _ _ _ _ _ _ _ _ _ _ _ _

Clay

. . . . . . . . . . . . . . . . . . . . . . .

Sand

From Microtremor

b)

200

0

_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ . . . . . . . . . _ _ _ _ _ _ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . _ _ _ _ _ _ . . . . . . .. . . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ . . . . . . . . . . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ . . . . . . . . . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _

100

700

From SPT

BMKG

20

Depth (m)

Depth (m)

200

0

400

Clay and sand (Japan) Clay and sand (Japan) Clay and sand (Japan) Clay (USA)

Shear Wave Velocity (m/s)

1000 1250

0 100

Soil Type

Depth (m)

Table 1

40 From SPT

c)

From Microtremor

d)

Figure 8 Comparison of 1D S-wave velocity profiles derived from microtremor array measurements (dash lines) and SPT (solid lines) at BMKG. a) and b) are for deep layers, c) lithology, and d) for shallow layers.

4.4

Two-dimensional S-Wave Velocity Profiles

The results of the individual inversions of the dispersion curves for each site show a lateral variation in the S-wave velocity. The 1D S-wave velocity profiles in central and northern Jakarta that resulted from the first inversion indicate a

322

Mohamad Ridwan, et al.

four-layer model over the engineering bedrock, while in the southern part there is a three-layer model. Basically, it is difficult to conduct a stratigraphic correlation based on such S-wave velocity structures. Hence, the second inversion was applied to the dispersion curve with a constant velocity as the initial model (cf. Yamanaka and Yamada [15]). According to the 1D S-wave velocity profiles (Figure 6), the four-layer model was assumed for the initial model with the S-wave velocity of the second, the third and the fourth layers being 500, 700 and 900 ms-1, respectively, determined by averaging the S-wave velocity of each layer. This initial model was inverted again for each dispersion curve to obtain the depth of each fixed velocity. The results are shown in Figure 9. S-wave Velocity (m/s) 0

250

500

-500

-750

-1000

-1000

0

250

500

-500

DRN3

Depth (m)

-500

500

-500

500

-500

-750

-1000

0

250

500

S-wave Velocity (m/s)

750 1000 1250

0

0

KBGS

250

500

750 1000 1250

0

ZPUR

UI

-250

-500

750 1000 1250

KNGN

S-wave Velocity (m/s)

750 1000 1250

-250

-500

250

-250

-1000

250

0

CPT

0

-750

0

-250

S-wave Velocity (m/s)

750 1000 1250 0

S-wave Velocity (m/s)

750 1000 1250

0

-250

500

MRIA

-1000

0

250

-250

S-wave Velocity (m/s)

750 1000 1250

Depth (m)

500

0 0

-750

S-wave Velocity (m/s) 250

S-wave Velocity (m/s)

750 1000 1250

-250

Depth (m)

Depth (m)

Depth (m)

-750

0

500

ISTN

-250

-500

250

BMKG

-250

Depth (m)

0 0

ANCL

-500

S-wave Velocity (m/s)

750 1000 1250

0

Depth (m)

750 1000 1250

Depth (m)

500

Depth (m)

250

-250

Depth (m)

S-wave Velocity (m/s) 0 0

-500

-750

-750

-750

-750

-750

-1000

-1000

-1000

-1000

-1000

Figure 9 1D S-wave velocity profiles of the 10 surveyed sites derived by second inversion using the initial model, in which the S-wave velocity for each layer is fixed.

The 1D S-wave velocity profiles resulting from the second inversion were then utilized to perform a stratigraphic correlation. The result of S-wave velocity correlation along the ANCL-UI line forms a cross section depicting a subsurface model for Jakarta (Figure 10). This model consists of four layers with the S-wave velocity from the second to the fourth layers being 500, 700, and 900 ms-1, respectively, while the S-wave velocity varies laterally in the first layer. A subsurface model for a similar north-south line derived from borehole data down to 350 m depth was conducted by Fachri, et al. [3]. The comparison of the stratigraphic model from the microtremor measurements and their geological cross section derived from the borehole data [3] indicates a good correlation (see Figure 11).

Estimation of S-wave Velocity Structures in Jakarta

NORTH

323

SOUTH

Meter

Meter

100 ANCL

BMKG

MRIA

ISTN

KNGN

0 218

264

263

282

283

DRN3 296

CPT 360

ZPUR

KBGS

398

341

UI

100

443

0 -100

-100

-200

-200

500 m/s -300

-300

-400

-400

700 m/s -500

-500

-600

-600

900 m/s -700

-700

-800

-800

Figure 10 Cross section of the S-wave velocity structure along the ANCL-UI line. The solid line depicts the estimated bedrock depth. NORTH

SOUTH

Meter

Meter

100 ANCL

0

BMKG

218 264

ISTN 263

MRIA 282

283

DRN3 296

CPT 360

KBGS 341

500 m/s

-300

UI

100

443 398

1

0

2

3 4

-100

-200

KNGN

ZPUR

-100

-200

-300

-400

-400 700 m/s

-500

-500

-600

-600 900 m/s

-700

-700

-800

-800

Figure 11 Comparison of microtremor analysis (solid and dash lines) and geological cross section derived from borehole data (modified from Fachri, et al. [3]).

324

Mohamad Ridwan, et al.

Referring to SNI-1726-2012 [9], in which for engineering purposes the bedrock is defined by Vs > 750 ms-1, the top of the fourth layer on the cross section, which has an S-wave velocity of about 900 ms-1, represents the engineering bedrock depth. In most cases, the subsurface model along the north-south line shows a bedrock morphology that is slanting towards the northern part of Jakarta with a depth range of 359-608 m.

5

Discussion

The results of the microtremor array analysis of each site show various dispersion curves that have been influenced by the geological conditions, as expressed by the variation in S-wave velocities. The section along the northsouth line across Jakarta (Figure 10) depicts the subsurface model consisting of four layers, in which the first layer has S-wave velocities in the range of 218443 ms-1, associated with Quarter sediments. The second and third layers have S-wave velocities of about 500 and 700 ms-1, respectively, associated with Tertiary sediments. The fourth layer, which has an S-wave velocity of about 900 ms-1, is estimated as the bedrock for earthquake engineering purposes and has a depth range from 359-608 m. Comparison with the previous geological cross section from Fachri, et al. [3], depicted in Figure 11, shows that the results of the microtremors analysis are in good agreement, whereby the sediment thickness increases northward. Nevertheless, the presence of faults around the Babakan site at a depth of 100300 m, as suggested by Fachri, et al. [3], cannot be identified well by the microtremor data. Therefore, future work is needed to add more measurement data from sites around the envisaged fault area. In brief, the available limited data derived from other methods cause difficulties in verifying the results of our microtremor data analysis for each site. However, the borehole data down to a depth of 150 m at the BMKG site, which has SPT data until 100 m depth, could be used to verify the microtremor data analysis results. The comparison of the S-wave velocity profiles derived from both microtremor and borehole data shows that these profiles agree well (see Figure 8). All in all we note that more information about subsurface conditions in Jakarta is essential for further studies.

6

Conclusions

This study obtained a subsurface model for Jakarta by using microtremor array measurements. Based on the results of measurements at 10 surveyed sites, the stratigraphic model for Jakarta consists of 4 layers, with S-wave velocity values ranging from 218 ms-1 to 443 ms-1 in the uppermost layer, and increasing from

Estimation of S-wave Velocity Structures in Jakarta

325

the second to the forth layers, i.e. 500 ms-1, 700 ms-1, and 900 ms-1, respectively. In general, the thickness of each layer increases northward, particularly the second layer. This result is in good agreement with the geological conditions in Jakarta, where sediments are very thick in the northern part. From the cross section (Figure 10), this study identifies the bedrock position in the fourth layer. The depth of the bedrock is diverse, within a range of 359-608 m. The deepest position of the bedrock is in the northern Jakarta area (ANCL), which reaches 608 m. Therefore it is expected that seismic wave amplification will be higher in the north than in the south, due to the very thick sediments.

Acknowledgments The authors gratefully acknowledge the Research Institute for Human Settlement (RIHS), Ministry of Public Works, Republic of Indonesia and the team of microzonation projects for supporting data used in this study. This research is part of the microzonation study for Jakarta that is supported by RIHS and ITB through the “Riset Unggulan” program.

References [1] [2]

[3] [4] [5] [6] [7]

[8]

Padmosukismo, S. & Yahya, I., The Basement Configuration of the North West Java Area, IPA Proceedings of the 3rd Annual Convention, Jakarta, 1, pp. 129-152, 1974. Turkandi, T., Sidarto, Agustiyanto, D.A. & Hadiwodjoyo, M.M.P., Geological Map of Jakarta and Kepulauan Seribu Quadrangles Jawa, Geological Research and Development Centre, Ministry of Energy and Mineral Resources, 1992. Fachri, M., Djuhaeni, Hutasoit, L.M. & Ramdhan, A.M., Stratigraphy and Hydrostratigraphy of Ground Water Basin for Jakarta, Bulletin Geology, 34(3), pp. 169-189, 2002. (Text in Indonesian) Aki, K., Space and Time Spectra of Stationary Stochastic Waves, with Special Reference to Microtremors, Bulletin of the Earthquake Research Institute, University of Tokyo, 35, pp. 415-456, 1957. Okada, H., Microtremors as an Exploration Method, Geoexploration Handbook, 2, Society of the Exploration Geophysicists of Japan, 1998. Okada, H., The Microtremors Survey Method, Geophysical Monograph Series, 12, 2003. Morikawa, H., Sawada, S. & Akamatsu, J., A Method to Estimate Phase Velocities of Rayleigh Wave Using Microseisms Simultaneously Observed at Two Sites, Bulletin of the Seismological Society of America, 94(3), pp. 961-976, 2004. Yamanaka, H. & Ishida, H., Application of Genetic Algorithms to an Inversion of Surface Wave Dispersion Data, Bulletin of the

326

[9] [10] [11] [12] [13] [14]

[15]

Mohamad Ridwan, et al.

Seismological Society of America, 86(2), pp. 436-444, 1996. SNI3-03-1726-2012, Guidelines of Earthquake Resistant Design for Building Structures, The National Standardization Agency of Indonesia, 2012. (Text in Indonesian) Irsyam, M., Sahadewa, A. & Darjanto, H., Soil Dynamics and Machine Foundation, Institut Teknologi Bandung, ISBN 978-979-1344-24-1, 2008. (Text in Indonesian) Imai, T. & Tonouchi, K., Correlation of N-value with S-Wave Velocity and Shear Modulus, Proceedings of the 2nd European Symposium on Penetrating Testing, Amsterdam, Netherlands, pp. 67-72, 1982. Imai, T., P- and S-Wave Velocities of the Ground in Japan, Proc. 9th, Int. Conf. Soil Mechanics and Foundation Eng., Tokyo, pp. 257-260, 1997. Ohta, Y. & Goto, N., Empirical Shear Wave Velocity Equation in Terms of Characteristic Soil Indexes, Earthquake Engineering and Structural Dynamics, 6, pp. 167-187, 1978. Sykora, D.W. & Stokoe, K.H., Correlation of In Situ Measurement in Sands of Shear Wave Velocity, Soil Characteristics, and Site Conditions, Geotechnical Engineering Report GR83-33, The University of Texas, Austin, 1983. Yamanaka, H. & Yamada, N., Estimation of 3D S-wave Velocity Model of Kanto Basin, Japan, Using Rayleigh Wave Phase Velocity, Bulletin of the Earthquake Research Institute, University of Tokyo, 81, pp. 295-301, 2006.