Ninth World Conference on Earthquake Engineering August 2-9, 1988, Tokyo-Kyoto, Japan

ANALYSIS OF STRONG GROUND MOTION OF PLAIN AND BASIN BEING COMPOSED OF SOFT SOIL BY FAULT MODEL AND BOUNDARY ELEMENT METHOD
Yoshiaki HISADA, Shunroku YAMAMOTO, and Sukenobu TANI
Tani Laboratory (51-18F-04), School of Science and Engineering, Waseda Univ., Ohkubo 3-4-1, Shinjuku-ku, Tokyo 160, Japan


SUMMARY
SH and Love waves in Kanto sedimentary basin are investigated by 2-D BEM including seismic source fault. Near sources are treated, and following conclusions are obtained: 1) Waveforms in the basin strongly depend on source location to the basin; 2) In particular, when a source is located shallowly and laterally away from the basin, large Love waves are generated at the basin edge; 3) The Love waves predominate only in comparatively long period. Thus, they probably influence long-period structures on the basin, if the natural period of the structures is nearly equal to the dominant period of the Love waves.

INTRODUCTION
Recently, earthquake damage of long-period structures on alluvial plain or sedimentary basin is on the increase. For example, 1) oil overflow of storage tanks in Niigata city, Japan during 1983 Nihonkai-chubu earthquake, 2) severance of elevator tail cords of high-rise buildings in Shinjuku, Tokyo during 1984 Naganoken-seibu earthquake (Ref.1), 3) collapse of multistory buildings in Mexco city during 1985 Michoacan earthquake. Many large cities are situated on soft sediment basins, and a lot of long-period structures are located on them. Thus, to investigate seismic waves of longer period in sediment basin is very important for earthquake engineering.
In order to simulate longer period (1-10 sec) seismic waves in Kanto sedimentary basin (in and around Tokyo, Japan), it is necessary to deal with a huge-scale basin model including a source fault. Because, seismic waves observed on the basin indicate that the longer period waves strongly depend on following three factors: 1) physical property of both bedrock and thick sediments (2-3km depth), 2) three dimensional irregular shape of the basin, 3) source mechanism and its location to the basin (e.g. Refs.2,3).
Various types of methods including a source have been proposed to simulate seismic waves in a basin (e.g. Refs.4,5,6). Recently, boundary element method (BEM) is widely used to compute wave propagation and scattering in a basin (e.g. Refs.7,8). Because it naturally contains the radiation condition of waves in a infinite domain, and it needs much fewer discrete elements than domain element methods, such as finite difference method (FDM) and finite element method (FEM). Therefore, we proposed 2-D BEM including a source fault in 1987 (Ref.9).
The purpose of this paper is to study longer period SH and Love waves in Kanto sedimentary basin with the BEM including a source, and to discuss influence of the waves on long-period structures in the sedimentary basin.

BOUNDARY ELEMENT METHOD INCLUDING SEISMIC SOURCE FAULT
Fundamental Equation of BEM Including Source Fault (Ref.9) Anti-plane waves in 2-D medium are treated. The medium consists of a bedrock and a basin sediment, as illustrated in Fig.1. The bedrock is assumed to be elastic, isotropic and homogeneous material. Its boundaries are composed of free surface, infinite boundary, internal boundary and irregular shaped boundary between bedrock and basin sediment. The free surface and infinite boundary are not treated, when we formulate the integral equation of this model with the fundamental solution of semi-infinite medium. The internal boundary is formulated as the dislocation model of the crack theory or seismology. In this study, the dislocation is modeled after Haskell's source model (Ref.10) as illustrated in Fig.2. Under the foregoing assumptions, the boundary integral equation of the bedrock in frequency domain is given by
...............(1)
...............(2)
...............(3)
...............(4)
...............(5)
...............(6)
where, V:displacement, P:stress on boundary, X,Y:point on smooth boundary, Xf:point on source fault, V':the fundamental solution (assuming exp(i*ƒÖ*t)), P':stress of V' on boundary, i:imaginary number, K(=ƒÖ/Vs):wave number, HN:Hankell function of the Nth order of the second kind, and W,ƒÂ,D,ƒÑ,Vr:fault parameters (width, dip angle, fault displacement, rise time of ramp function, rupture velocity).
The irregular shaped boundary between the bedrock and the basin sediment is treated by the BEM. In this study, the basin sediment is also formulated by the BEM.
Quality Factor (Internal Damping) The quality factor (Q) is introduced into the model by multiplying R1 and R2 of the fundamental solation (eq.(3)) by the following term.
...............(7)

SH AND LOVE WAVES ANALYSIS IN KANTO SEDIMENTARY BASIN
SH and Love Wave Analysis in the Basin by the BEM According to observation of seismic waves on Kanto basin, it is certain that the seismic waves in the basin are greatly effected by source location to the basin, as mentioned in INTRODUCTION. In particular, when the sources are located under the basin (source type A), the waves seem to consist of only body waves and the duration is short. On the contrary, when the sources are located shallowly and laterally away from the basin (type B), large surface waves are observed inside the basin and the duration becomes much longer than the one of the type A (e.g. Refs.2,3).
We simulate SH and Love waves in Kanto basin by the BEM and explain the foregoing phenomena. Kanto basin model and two type sources A and B are illustrated in Fig.3. The model consist of a bedrock and a basin sediment. The right edge of the sediment is not discretised in order to transmit almost all waves in right direction. In the basin, dominant period of SH wave is 10 sec. Type A source locates under the basin, and type B locates shallowly and laterally away from the basin. A simple pulse wave is generated at those sources by using fault parameters listed in Fig.3. Fourier amplitude spectrum of the source becomes zero at circular frequency (ƒÖ)=ƒÎ, and ƒÖ<ƒÎ (period>2 sec) is treated. However, static term (ƒÖ=0) is neglected, because displacement amplitude diverges due to the 2-D source. Thus, amplitude of calculated displacement waves does not start from 0.
Fig.4 shows displacement waves and Fourier amplitude spectra of velocity; they are simulated on the basin. Dotted lines are computed by the BEM, and solid lines are by 1-D multiple reflection method (only vertical SH waves in the sediment are dealt with). The waves simulated at every 10km are vertically arranged. The top waves are obtained at the nearest point to the basin edge. The four left waves are obtained by type A source. The waves almost consist of SH waves vertically traveling from bedrock, and the duration of the strong part is not long. Waves and spectra of the BEM are almost equal to ones of the 1-D method. The four right waves are obtained by the type B source. In this case, the duration becomes longer, as the observational points are located away from the basin edge. This is because large waves are generated at the edge and they are slowly propagated through the sediment. In this case, the 1-D method is only applicable to SH wave part, and it underestimates Fourier amplitude spectra of the BEM waves. These results agree with the ones of the foregoing observation on Kanto basin.
Identification of Love Wave On the analogy of several researchers' studies (e.g. Refs.11,12), the large waves following SH waves are expected to be Love waves; they are generated at the basin edge. In this study, the waves are identified as mainly fundamental Love mode by a theoretical method. In the method, infinitely lateral sediment is assumed and a point force is given on free surface to generate Love waves. The force point corresponds to the basin edge point on free surface. Fourier spectrum of the point force is obtained by multiplying the source spectrum by a correction factor of amplitude. The equation is given by,
...............(8)
where, X:lateral length from the force point (basin edge point) to observational point, V1:the theoretical Love wave displacement, Vf:the source spectrum of displacement, Rf:length from the source to the basin edge, T:transmission coefficient from the bedrock to the sediment, A1:Harkrider's Love wave amplitude response (Ref.13), C1:the amplitude correction factor (4*ƒÊ is used in this study, where ƒÊ is shear modules of the sediment), C:phase velocity of Love wave.
Fig.5 shows comparison of the BEM waves and the theoretical Love waves, in case of the source type B. The large BEM waves following SH waves nearly equal to the theoretical Love waves. Thus, it is confirmed that the latter waves are mainly the Love waves of fundamental mode and they are generated at basin edge.
Influence of Internal Damping (Q) on the Waves Fig.6 shows comparisons of Fourier spectrum and wave with and without Q in case of the source type B. Solid line is calculated with Q and doted line is without Q. The Q is assumed as 50 in sediment and 100 in bedrock. The wave comparison indicates that SH waves are little attenuated, because they shortly travel in low-Q sediment. On the contrary, the Love waves, especially short period components, are largely attenuated, because they travel long through the sediment. However, it is confirmed by the spectrum that long-period components are little attenuated. The similar result was obtained by Horike using the Aki-Larner method in 1987 (Ref.14).
Characteristics of Love Waves in a Sediment Basin and Influence of the Waves on Structures on the Basin As Drake had already pointed out with a FEM alluvial valley model (Ref.12) and we also confirmed in this study, the Love waves mainly propagate longer period components. Because they consist of principally fundamental mode (its dominant period is about 9 sec in this model) and their short period components are largely attenuated in the low-Q sediment. Therefore, short-period structures on the basin are not probably influenced by the waves. However, long-period structures are likely influenced, if the natural period of the structures is nearly equal to the dominant period of the Love waves. Although we investigate only one case in this study, it is analogized that the long-period structures' damage, as mentioned in INTRODUCTION, is partially caused by the local surface waves in the basins.

CONCLUSIONS
SH and Love waves in Kanto sedimentary basin are studied by the 2-D BEM including source fault. The basin is modeled as a single layer (dominant period of SH wave is 10 second). Near sources are treated and the waves of 2 sec and over period are simulated. With the BEM simulation and the theoretical Love wave analysis, following conclusions are obtained;
1) Waveform in the basin are strongly effected by source location. When a source is located under the basin, the waves in the basin consist of mainly SH waves; the waves are vertically traveling from the underlying bedrock. In this case, it is reasonable to use the 1-D multiple reflection method of SH waves. On the contrary, when a source is located shallowly and laterally away from the basin, the waves consist of not only the SH waves, but also Love waves. The Love waves are generated at basin edge and slowly propagated through the basin sediment. In this case, the 1-D method is possible to simulate only SH waves part.
2) The Love waves in the sediment basin predominant only in comparatively long period, because they consist of mainly fundamental mode and their short period components are largely attenuated in the low-Q (high-damped) sediment. Consequently, the Love waves probably influence long-period structures on the basin, if the natural period of the structure is nearly equal to the dominant period of the Love waves.

ACKNOWLEDGMENTS
The authors would like to thank Prof. Kanto Seo and Mr.Makoto Yamada for their helpful information on the seismic waves in Kanto basin, Dr.Keiichi Kasahara and Prof. Shigeru Noda for useful suggestions on this study, and Profs. Shouichi Kobayashi and Michio Iguchi for valuable advice on the BEM technique.

REFERENCES
  1. Muramatsu,I., "Shindo to Higai to Jishindo (Sono 1)", Abstracts of SSJ, No.1, pp.280, 1988, (in Japanese)
  2. Yamada,M., Nasu,N., Takeuchi,M. and Morioka,T., "Some characteristics of the Earthquake Ground Motions Observed in Tokyo, Especially those of Just-under-seated Earthquake", Proc. of the 6th Japan Earthq. Eng. Symp. 121-128, 1982
  3. Seo,K. and Kobayashi,H., "The Influences of GroundStructure on Earthquake Ground Motions", 3rd International Earthq. Microzonation Conf., 637-648, 1982
  4. Bouchon,M. and Aki,K., "Near-Field of a Seismic Source in a Layered Medium with Irregular Interfaces", Geophs.J.R.astr.Soc, Vol.50, 669-684, 1977
  5. Kohketsu,K., "2-D Reflectivity Method and Synthetic Seismograms for Irregularly Layered Structures -1. Sh-Wave Generation", Geophys.J..R.astr.Soc., Vol.89, 821-838, 1987
  6. Vidale,J.E. and Helmberger,D.V., "Elastic Finite Difference Modeling of the 1971 San Fernando, California Earthquake", Bull.SSA, Vol.78, No.1, 122-141, 1988
  7. Sanchez-Sesma,F.J. and Eaquivel,J.A., "Ground Motion on Alluvial Valleys under Incident SH Waves", Bull.SSA, Vol.69, No.4, 1107-1120, 1979
  8. Dravinski.M., "Influence of Interface Depth upon Strong Ground Motion", Bull.SSA, Vol.72, 597-614, 1982
  9. Hisada Y., Yamamoto S. and Tani S., "A study of Seismic Waves on a Sedimentary Plain by Boundary Element Method and Fault Model", Proc.Symp. on Compu. Methods in Structural Eng. and Related Fields, 123-128, 1987, (in Japanese)
  10. Haskell,N.A., "Total Energy and Energy Spectral Density of Elastic Wave Radiation from Propagating Faults", Bull.SSA, Vol.54, No.6, 1811-1841, 1964
  11. Bard,P.-Y. and Bouchon,M., "The Seismic Response of Sediment-Filled Valleys. Part 1 : The Case of Incident SH Waves", Bull.SSA, Vol.70, No.4, 1263-1286, 1980
  12. Drake,L.A., "Love and Rayleigh Waves in an Irregular Soil Layer", Bull.SSA, Vol.70, No.2, 571-582, 1980
  13. Harkrider,D.G., "Surface Waves in Multilayered Elastic Media. Part 2. Higher Mode Spectra and Spectral Ratios from Point Sources in Plane Layered Earth Models", Bull.SSA, Vol.60, No.6, 1937-1987, 1970
  14. Horike,M., "Extension of the Aki and Larner Method to Absorbing Media with Plural Curved Interfaces and Several Characteristics of a Seismic Response on a Sedimentary Bas in", J.SSJ, Vol.40, No.2, 247-259, 1987, (in Japanese)