12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 1 Stochastic Ground Motion Simulation for Crustal Earthquakes in Japan Tatsuya Itoi Associate Professor, Graduate School of Engineering, the University of Tokyo, Tokyo, Japan ABSTRACT: In this study, an empirical simulation model of stochastic ground motion for crustal earthquakes in Japan is proposed based on the ground motion records observed by K-NET in Japan from 1997 to 2011. The proposed model is developed based on that by Rezaeian & Der Kiureghian (2010). The characteristics of ground motion depends on various parameters such as magnitude, distance from source to site and local site conditions (shallow soil as well as deep subsurface structure), which is included in the parameters in the proposed model. The proposed model is considered useful when conducting probabilistic risk analysis of structures such as base-isolated or vibration-controlled buildings. 1. INTRODUCTION For probabilistic seismic risk analysis of structure such as base-isolated or vibration-controlled structures, time history analysis is employed. Ground motion time history is needed for the analysis. For conventional seismic risk analysis, artificially-generated ground motion time history where amplitude is adjusted so that it fits ground motion intensity with some exceedance probability is used. Exceedance probability of the ground motion intensity measure for amplitude adjustment is obtained based on probabilistic seismic hazard analysis for the intensity measure. The ground motion parameter such as its duration, however, is conservatively assumed frequently, which prevents realistic estimate by overestimating or underestimating the risk of the structure. On the other hand, ground motion simulation considering the fault rupture process such as the empirical Greenβs function method is widely used. It is, however, not appropriate to be utilized in probabilistic seismic hazard analysis. This is because it requires much time and resource as well as because there is no established method for variability in simulated ground motion. Therefore, in this study, an empirical simulation model for ground motion time history is developed based on that proposed by Rezaeian & Der Kiureghian (2008, 2010), and an empirical attenuation relation of ground motion time history is proposed, based on the observed records by K-NET in Japan. 2. MODELLING OF GROUND MOTION 2.1. Overview In this study, a model for velocity time history of ground motion is developed based on Rezaeian & Der Kiureghian (2008, 2010). Velocity time history π£(π‘) is assumed as follows: π£(π‘) β π(π‘, πΆ)π’0(π‘) (1) where, π(π‘, πΆ) is a time modulating function and π’0(π‘) is the stochastic process where variance is constant in time. The gamma distribution is employed for the time modulating function π(π‘, πΆ) as follows: π(π‘, πΆ) β π‘πΌ1exp(βπΌ2π‘) (2) where, πΌ1 and πΌ2 are the constants. These constants are determined by determining π‘π, the time where the amplitude of envelope reaches maximum, and π‘π, the time where the amplitude 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 2 of envelope becomes 0.1 times the maximum amplitude. The amplitude of velocity time history π£(π‘) is defined by integrated squared velocity πΌπ as follows: πΌπ = β« π£2(π‘)ππ‘β0 (3) The stochastic process π’0(π‘) in Equation (1) is modelled by the weighted average of two stochastic processes, π’1(π‘) and π’2(π‘), as follows: π’0(π‘) = βπ(π‘)π’1(π‘) + β(1 β π(π‘))π’2(π‘) (4) where, the direct S waves are modeled using π’1(π‘), while the later phases are modeled using π’2(π‘). The stochastic process, π’1(π‘) and π’2(π‘), respectively are modeled as the response of 1 DOF system with different natural frequency ππ and damping coefficient ππ which is excited by Gaussian process π€π(π) as follows: π’π(π‘) β β« βπ(π‘ β π|ππ, ππ)π€π(π)πππ‘ββ (5) The weight π(π‘) in Equation (4) is assumed proportional to time t as follows: π(π‘) = {1 βπ‘π‘π (0 < π‘ β€ π‘π)0 (π‘ > π‘π) (6) The model proposed here needs eight parameters, πΌποΌπ1οΌπ2οΌπ1οΌπ2οΌπ‘cοΌπ‘ποΌπ‘π , to simulate velocity time history of ground motion. 2.2. Ground motion parameters Ground motion parameters, πΌποΌπ1οΌπ2οΌπ1οΌπ2οΌπ‘cοΌπ‘ποΌπ‘π, are estimated for each recorded ground motion time history. The procedure for parameter estimation is almost same as that in Rezaeian & Der Kiureghian (2008, 2010). 3. GROUND MOTION DATABASE Database for ground motion time histories recorded at K-NET stations at fault distances less than 100 km is compiled from crustal earthquakes of moment magnitude (Mw) larger than 5, which occurred from 1997 to 2011. Figure 1 shows the location of earthquake hypocenter, while the relation between magnitude and fault distance for each record is shown in Figure 2. The range of site condition for recorded motions is 200m/s β€ππ30 β€ 700m/s and π1500 β€ 2000m. Ground motion parameters, πΌποΌπ1οΌπ2οΌπ1οΌπ2 οΌ π‘c οΌ π‘π οΌ π‘π , are estimated for each component record. e Figure 1: Location of earthquake hypocenter Figure 2: Relation between magnitude and fault distance for ground motion database 4. PREDICTION EQUATION 4.1. Transformation of probability distribution for ground motion parameters First, statistics for each ground motion parameters are analyzed. Table 1 shows the probability 0 20 40 60 80 10055.566.57Fault distance (km)MW Strike-slipNormalReverseReverse(oblique)12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 3 distribution and its parameters for each ground motion parameter. Table 1 shows the parameters for probability distribution for each ground motion parameter estimated by the maximum likelihood method. Table 1: Parameter values. Then, each ground motion parameter is transformed so that it fits the standard normal distribution as follows: π£π = π·β1 (πΉπ©π(ππ)) (7) where, π·β1( ) is the inverse of the standard normal distribution. The probability density function for the gamma distribution is as follows: ππ(π₯) =1π½πΌπ€(πΌ)π₯πΌβ1πβπ₯π½ (8) where, the πΌ and π½ are distribution parameters, and π€( ) is the gamma function. The probability density function for the beta distribution is as follows: ππ(π₯) =1π΅(π, π)π₯πβ1(1 β π₯)πβ1 (9) where, and π and π are the distribution parameters, and π΅( ) is the beta function. The probability density function for the lognormal distribution is as follows: ππ(π₯) =1π₯πlnπβ2πππ₯π {β(lnπ₯ β πlnπ)22πlnπ2} (10) where, πlnπ and πlnπ are the distribution parameters. Mean ground motion parameters π£AMπ , which are the average of those for each component, π£NSπ and π£EWπ , is calculated as follows: π£AMπ =12(π£NSπ + π£EWπ) (11) Standard deviation of π£NSπ β π£AMπ is estimated as follows: πcompπ = βVar(π£NSπ β π£AMπ) (12) Table 2 shows the estimated standard deviation for each ground motion parameter. Table 2: Standard deviation of π£NSπ β π£AMπ 4.2. Prediction equation In this section, a prediction equation for each ground motion parameter is obtained. The explanatory variables for the equation are those for source, propagation path and site amplification characteristics as follows: π£AMπππ = π(πππ, π·π, π ππ , ππ30π, π1500π)+ πππ (13) where, πππ , π·π , π ππ , ππ30π and π1500π are moment magnitude, hypocenter depth, shortest fault distance, 30 meter average shear wave velocity and depth to the layer with shear wave velocity of 1500m/s. The regression equation for each ground motion parameter is as follows: π£AM1= π10 + π11ππ6.0+ π12π·10km+ π13log10 (π + π1610π17ππ40km)+ π14log10 (ππ30400m/s)+ π15min(π1500, ππ9)100m+ π1 (14) π£AMπ= ππ0 + ππ1ππ6.0+ ππ2π·10km+ ππ3π 40km+ ππ4log10 (min(ππ30, ππ8)400m/s)+ ππ5log10 (π1500100m) + ππ (π = 2,3) (15) (unit) Distibution BoundI V (v 1) (m2/s) Lognorm l (0, β) mlnEv -8.308 slnEv 2.777f 1 (v 2) (Hz) Gamma (0, β) af1 4.549 bf1 0.8055f 2 (v 3) (Hz) Gamma (0, β) af2 1.597 bf2 1.328z 1 (v 4) - Beta (0, 1) qz1 1.0147 rz1 5.484z 2 (v 5) - Beta (0, 1) qz2 0.8117 rz2 2.553t c (v 6) (s) Gamma (0, β) atc 3.707 btc 5.143t p (v 7) (s) Gamma (0, β) atp 1.415 btp 3.226t d - t p (v 8) (s) Lognormal (0, β) mln(td-tp) 3.488 sln(td-tp) 0.8019Distribution parametersParamet rv 1 v 2 v 3 v 4 v 5 v 6 v 7 v 80.083 0.385 0.194 0.569 0.389 0.627 0.311 0.23312th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 4 π£AMπ = ππ0 + ππ1ππ6.0+ ππ2π·10km+ ππ3π 40km+ ππ4log10 (ππ30400m/s)+ ππ5π1500100m+ ππ (π = 4,5,7) (16) π£AMπ = ππ0 + ππ1ππ6.0+ ππ2π·10km+ ππ3min(π , ππ6)40km+ ππ4log10 (ππ30400m/s)+ ππ5min(π1500, ππ9)100m+ ππ(π = 6,8) (17) The regression results and the standard deviation of residuals are shown in Table 3, while the correlation coefficients between residuals are shown in Table 4. 4.3. Simulation of ground motion by proposed equations Ground motion time history of north-south component is simulated using the proposed equations. Figure 3 shows some samples for velocity and acceleration time histories for simulated ground motions. Ground motions are simulated for the source with πποΌ6.5, π·=15km and the site for π =10km, π1500 =1000m, ππ30=500m/s. Figure 3 shows the comparison with respect to acceleration response spectra between the simulated ground motions and existing attenuation equations. Simulated acceleration response spectra are higher in the short period range. Effects of deep subsurface structure are simulated by the proposed model, while existing attenuation equations do not. Figure 4 shows the amplification of response spectrum (T=2s) due to deep subsurface structure, which is harmonic with the existing empirical relationship proposed by Itoi & Takada (2012). Table 3: Regression coefficients and standard deviation of prediction error for Equations (14) β (18) Table 4: Correlation coefficients between prediction errors for Equations (14) β (18). i 0ai 1ai 2ai 3ai 4ai 5ai 6ai 7ai 8ai 9seiIVvGM1-8.046 8.400 0.254 -2.193 -1.498 0.099 31.65 0.394 - 444.1 0.386f1vGM23.022 -3.248 0.163 0.207 0.914 -0.330 - - 393.9 - 0.863f2vGM33.558 -3.212 0.301 -0.121 2.371 -0.673 - - 253.6 - 0.718z1 GM4-1.460 2.015 -0.133 -0.161 1.499 -0.009 - - - - 0.851z2vGM5-0.640 0.807 -0.144 0.096 2.099 -0.021 - - - - 0.803tcvGM6-3.726 2.215 -0.054 0.730 -0.712 0.118 65.77 - - 664.0 0.934tpvGM7-4.885 4.706 -0.315 0.303 -0.527 0.051 - - - - 0.690td-tpvGM8-3.865 2.116 -0.263 1.432 -0.648 0.185 48.81 - - 640.0 0.723GMPvGM1vGM2vGM3vGM4vGM5vGM6vGM7vGM8vGM11 -0.229 .061 - .075 -0.311 -0.162 0.008 -0.203vGM21 0.332 -0.338 0.287 -0.027 -0.006 -0.072vGM31 -0.134 -0.221 -0.176 -0.227 -0.537vGM41 0.157 -0.181 0.068 0.103vGM51 0.010 -0.139 -0.039vGM6sym. 1 0.109 0.444vGM71 0.330vGM8112th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 5 Figure 3: Samples of simulated ground motions (north-south component) (πποΌ6.5, π·=15km and the site for π =10km, π1500=1000m, ππ30=500m/s) 5. CONCLUSIONS In this study, an attenuation relation for velocity ground motion time history is proposed for crustal earthquake in Japan. The proposed model is harmonic with existing attenuation relation for acceleration response spectra. The proposed relation can be applied to crustal earthquake with 5.1 β€ ππ β€ 6.9, and at the site for the shortest distance π β€ 100km and the soil condition 200m/s β€ ππ30 β€ 700m/s and π1500 β€ 2000m. The proposed model, however, is not applicable to the near-field ground motion at a soft soil site, because the non-linearity of soil amplification is not included in the proposed model, which will be discussed in the future study. 6. ACKNOWLEDGEMENT Ground motion data are provided by NIED (K-NET). (a) π =10km π1500=20m (b) π =10km π1500=1000m Figure 4: Comparison between simulated ground motions and conventional attenuation relation for response spectra (πποΌ6.5, π·=15km, ππ30=500m/s) Figure 5: Amplification for Sa(T=2s) due to deep subsurface structure (πποΌ6.5, π·=15km, ππ30=500m/s) 0 20 40-200002000cm/s20 20 40-1000100cm/s0 20 40-200002000cm/s20 20 40-1000100cm/s0 20 40-200002000cm/s20 20 40-1000100cm/s0 20 40-200002000cm/s20 20 40-1000100cm/s0 20 40-200002000cm/s2time(s)0 20 40-1000100cm/stime(s)10-110010-1100101102103Period(s)Pseudo velocity(h=0.05)(cm/s) Kanno et al.Uchiyama & MidorikawaSatohSimulated (median)Simulated (+-sigma)10-110010-1100101102103Period(s)Pseudo velocity(h=0.05)(cm/s) Kanno et al.Uchiyama & MidorikawaSatohSimulated (median)Simulated (+-sigma)0 500 1000 1500 200011.522.5Z1500(m)Amplification due to deep subsurface structure Simulated (median)Itoi & Takada (2012)12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 6 7. REFERENCES Itoi, T. & Takada, T. (2012). βEmpirical correction of ground motion prediction equations of response spectra at rock sites for near-field earthquakes considering amplification effect with deep subsurface structureβ, Journal of Japan Association for Earthquake Engineering, 12(1), 1_43-1_61. (in Japanese) Kanno, T. et al. (2006). βA New Attenuation Relation for Strong Ground Motion in Japan Based on Recorded Dataβ Bulletin of the Seismological Society of America, 96(3), 879β897. Rezaeian & Der Kiureghian (2008). βA stochastic ground motion model with separable temporal and spectral nonstationaritiesβ, Earthquake Engineering and Structural Dynamics, 37(13), 1565-1584. Rezaeian & Der Kiureghian (2010). βSimulation of synthetic ground motions for specified earthquake and site characteristicsβ, Earthquake Engineering and Structural Dynamics, 39(10), 1155-1180. Satoh, T. (2008). βAttenuation relations of horizontal and vertical ground motions for P wave, S wave, and all duration of crustal earthquakesβ Journal of structural and construction engineering 73(632), 1745-1754 . (in Japanese) Uchiyama, Y. & Midorikawa, S. (2006). βAttenuation relationship for response spectra on engineering bedrock considering effects of focal depthβ Journal of structural and construction engineering (606), 81-88. (in Japanese)
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015) /
- Stochastic gound motion simulation for crustal earthquakes...
Open Collections
International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015)
Stochastic gound motion simulation for crustal earthquakes in Japan Itoi, Tatsuya Jul 31, 2015
pdf
Page Metadata
Item Metadata
Title | Stochastic gound motion simulation for crustal earthquakes in Japan |
Creator |
Itoi, Tatsuya |
Contributor | International Conference on Applications of Statistics and Probability (12th : 2015 : Vancouver, B.C.) |
Date Issued | 2015-07 |
Description | In this study, an empirical simulation model of stochastic ground motion for crustal earthquakes in Japan is proposed based on the ground motion records observed by K-NET in Japan from 1997 to 2011. The proposed model is developed based on that by Rezaeian & Der Kiureghian (2010). The characteristics of ground motion depends on various parameters such as magnitude, distance from source to site and local site conditions (shallow soil as well as deep subsurface structure), which is included in the parameters in the proposed model. The proposed model is considered useful when conducting probabilistic risk analysis of structures such as base-isolated or vibration-controlled buildings. |
Genre |
Conference Paper |
Type |
Text |
Language | eng |
Notes | This collection contains the proceedings of ICASP12, the 12th International Conference on Applications of Statistics and Probability in Civil Engineering held in Vancouver, Canada on July 12-15, 2015. Abstracts were peer-reviewed and authors of accepted abstracts were invited to submit full papers. Also full papers were peer reviewed. The editor for this collection is Professor Terje Haukaas, Department of Civil Engineering, UBC Vancouver. |
Date Available | 2015-05-15 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivs 2.5 Canada |
DOI | 10.14288/1.0076066 |
URI | http://hdl.handle.net/2429/53204 |
Affiliation |
Non UBC |
Citation | Haukaas, T. (Ed.) (2015). Proceedings of the 12th International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP12), Vancouver, Canada, July 12-15. |
Peer Review Status | Unreviewed |
Scholarly Level | Faculty |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/2.5/ca/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 53032-Paper_185_Itoi.pdf [ 578.33kB ]
- Metadata
- JSON: 53032-1.0076066.json
- JSON-LD: 53032-1.0076066-ld.json
- RDF/XML (Pretty): 53032-1.0076066-rdf.xml
- RDF/JSON: 53032-1.0076066-rdf.json
- Turtle: 53032-1.0076066-turtle.txt
- N-Triples: 53032-1.0076066-rdf-ntriples.txt
- Original Record: 53032-1.0076066-source.json
- Full Text
- 53032-1.0076066-fulltext.txt
- Citation
- 53032-1.0076066.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.53032.1-0076066/manifest