断層すべり方向と継続時間を考慮した鹿児島県北西部地震の地震動の再現
SYNTHESIS OF STRONG GROUND MOTION FOR THE 1997 KAGOSHIMAKEN-HOKUSEIBU
EARTHQUAKE

大角恒雄*・原田隆典**
Tsuneo OHSUMI and Takanori HARADAIn this paper,the motions from small earthquakes are produced by a stochastic model
where the spectrum amplitude is given from a stochastic source model with ω− 2 property
and the duration of motion depends on the observer relative to the fault. By summing thesimulated motions and making use of the scaling relations between large and small
earthquakes, the ground motions from an extended rupture too large to be treated as apoint source model are simulated. As a verif ication of the proposed method, simulated
ground motions are compared with the K-Net recorded ground motions during theKagoshimaken-hokuseibu earthquake（M6.3,1997）.
Key Words : Green's functions, stochastic simulation, empirical simulation, Kagoshimaken-hokuseibu earthquake

１．まえがき点に加えて、グリーン関数として用いることのできる適切観測された地震動記録には、震源の影響と、震源から観な観測波形がない建設地点においても震源特性を考慮した測点までに地震波が伝播する影響が反映されている。した

強震動の波形合成が可能となる。このような本方法の考えがって、大地震による地震動を予測する場合に、大地震の方は、すでにIrikura1）、Joyner and Boore（1986）2）、釜江・震源近くで発生した同じような震源特性を持つ小地震の観他（1991）3）、佐藤・他（1994）4）、原田・他（1995）5）が示してい測波形を用いると、この観測波形にはすでに地震波の伝播るが、上述したように地震波の伝播特性が近似的に評価さ特性が反映されているので、あとはこの小地震と予測の対れるという弱点があり、また、小地震によるグリーン関数

象である大地震の震源特性の違いのみを修正すればよいこに断層破壊方向と継続時間の関係が考慮されていないため、とになる。したがって、このような経験的グリーン関数法

これらの適用に当たっては、用いたモデル・パラメータとの長所は、小地震による観測波形をグリーン関数として用の関係を明確にしながら、実地震波形との比較を通じて問いるところにある。しかし、本文では、確率論的震源モデ題点や有効性を検討しておかなければならない。そこで、ルによる波形をグリーン関数として用い、大地震の断層破

本文では、提案モデルによる合成波形と実地震による観測壊伝播、断層の大きさと観測点の幾何学的関係等を考慮し波形を比較してモデルの適用性を検討した。て大地震による強震動波形を合成する。経験的グリーン関

そこで本文では、小地震によるグリーン関数において断数法は、現実地盤における地震波の伝播特性を評価してい層すべり方向と継続時間を考慮し、鹿児島県北西部地震るという長所を有するが、本方法では、この伝播特性が近（Ｍ＝6 . 3、1997.３.26）を再現した。比較検討では、余震似的に評価されるが、従来の方法のように大地震の断層破（Ｍ＝3.8、1997.３.27）による観測波形を用いる経験的グリ

壊伝播特性や観測点との幾何学的関係が考慮されるというーン関数法並びに、一様確率分布位相角を有する正弦波の級数和から計算するShinozuka・他6）の手法から計算した波形をグリーン関数として用いる確率論的グリーン関数法7）＊
中央研究所 開発研究部 主任研究員＊＊ 宮崎大学教授 工学部土木環境工学科 工博により本震の記録を再現し、防災科学技術研究所の地震観1

測ネット（K-Net）で得られた観測波形と比較検討するこ

このような断層破壊の不均質性の影響が混在した結果としととした。て地震動が観測されるため、記録の解析から一義的に、高振動数地震波に影響するこれらの種々の不均質性を分離す２．モデルの概要ることは不可能である。したがって、これらの物理現象を2.1

決定論的に仮定して、地震動予測を行なうことは困難であグリーン関数を用いた強震動の波形合成法

小地震の地震動u0（χ,ω）から近似的に大地震の地震動る。そこで、現状では経験的に予測できるような簡易なモu（χ,ω）
を求める計算式は次式で与えられる（Figure.１）。

デルが工学的に意義があるので、これらの断層破壊の不均質性を考慮する因子として、κによって震源スペクトルの高振動数成分が調整可能な式sのTmn（ω）の右辺の第２因子のようなモデルを提案している。ここでは日向灘地震の検この式やそのパラメータの詳細は原田・他（1995）に譲るが、討結果7)からκ＝３を用いた。

ここでは、後述の観測波形との比較において重要と思われ2.2る高振動数地震波に関する本モデルの特徴を記述する。

確率論的震源モデルによるグリーン関数

確率論的グリーン関数として用いる小地震の地震動波形
a（t）は非定常パワースペクトルSaa（t、ω）を用いて次式に式aを時間領域に変換したものは、Irikura（1986）と一致するが、伝達関数T mn（ω）の評価が異なる。伝達関数T mnより計算する。

（ω）とは、大地震と小地震のすべり速度時間関数の違いを表わす伝達関数であり、本論文では、次式のものを用いる。

ここに、ωuは振動数の上限値で、これ以上の振動数では非定常スペクトルの振幅が数学的または物理的理由から十分に小さく、零とみなされる振動数である。式dにおいて、
φjは（0, 2π）の間で一様分布を持つ互いに独立なランダム
ここに、VRは断層の破壊速度、CSはS波の伝播速度を表わ位相角である。非定常パワースペクトルは、次のように与す。式sの提案モデルは、指数関数形のすべり時間関数をえられるものとする。

採用し、震源スペクトルの高振動数領域の値をκによってコントロールできるモデルである。地震断層の破壊による高振動数地震波は、断層のすべりがスムーズでないこと、断層面上においてすべり量や応力降下量が一様でないこと、断層破壊速度が一定でないこと等によって発生することがここに、W（t,ω）は非定常特性を表現する修正関数を表わ考えられるが、これらの断層破壊の不均質性を取り入れた

すが、本文では、大崎（1994）8）の提案式のように時間のみ種々の物理モデルが提案されている。ここで重要な点は、の関数W（t,ω）＝W（t）を用いる。

Teは強震動の継続時間（有効継続時間：ほぼ定常部分の継続時間）で、次式から計算する（Ohsumi・他（1997））9）。

ここに、L：地表面への断層投影長さ（km）、VR：断層破壊速度（km/s）、CS：Ｓ波速度（km/s）。Td：強震動振幅が90パーセント以上の継続時間を表わす。式kの継続時間と地震マグニチュードMの関係を久田の経験式10）（logTd＝0.312

M−0.774）ともに示す（Figure.３）。式kでは、logL＝0.5
M−1.88、VR/CS＝0.7を用いた。断層線と観測点の角度θ
によって継続時間に違いが生じるが、経験式は理論式のほぼ中間値を与えている。
式hのC、A（ω）
、A（ω）
、A（ω）
は、それぞれ、振幅SDA
定数、震源スペクトル、伝播経路特性、およびサイト特性を表す係数である。振幅定数と震源スペクトルは次式で与えられる。

高振動数の断振動数に関しては、釜江・他の経験式3）を用いて以下のような値を決定した。
ここに、は、震源からの地震波の放射パターンに対する平均補正係数（0.63）、Fは自由表面の影響による波動の増地盤の材料減衰特性係数Q値は、Kiyono（1992）12）の式よ幅補正係数（2.0）、Vは２方向の水平成分のエネルギー分り以下のように定めた。

配を説明する低減係数（0.5）である。ρは媒質の密度、CSは媒質のせん断波速度である。
なお、コーナー振動数ωcは次式から定めた11）。

本文では、地震基盤から工学基盤面までの地震波の増幅特性AAをBoore（1983）11）により式¡3左辺の第１因子のように評価し、工学的基盤から地表までを第２因子として、金ここに、ωCの単位は（rad/sec）、m0の単位は（dyn-cm）、井・田治見スペクトルで評価した。ωgは表層地盤の卓越振

Δσは応力の単位（bars）で地震の規模によらず一定であると動数、nは増幅率を制御するパラメーターである。硬質地盤仮定する。

では、ωg＝15.6（rad/sec）、hg＝0.6を用いた。

減衰特性係数と表層地盤による地震波の増幅特性は次式

３．観測記録による検証例で与える。

3.1

経験的グリーン関数による波形合成

本方法の検証のためにK-Netを利用し、鹿児島県北西部地震（Ｍ＝6.3、1997.３.26）と、その余震（Ｍ＝3.8、1997.
３.27）記録を用いて検討した。
解析に用いた断層モデルは田中・他13）の地殻変動から求めた断層モデルを用いた（Figure.４、Table.１）。また、震源からの破壊伝播は、震源（Ｎ35.68°、Ｅ136.55°）を破壊開始点とし、角田・他13）によるとbilateralに広がったとされているが、ここでは、Vr＝2.5km/sで放射状に破壊伝播するものとした。
K-Netで得られた観測波形は、震源域の中心から８km以内を取り囲む方向に分布した６ヶ所（阿久根、出水、大口、宮之城、横川、川内：Figure.5）：KGS002〜007（Table.２）とした。経験的手法に用いた余震は、本震がＭ6.3とあまり規模の大きくない地震であるので、サブフォルトの分割数が６以上となるＭ＝3.8の余震を採用した。
3.2

確率論的グリーン関数による波形合成

次に、確率論的震源モデルから計算で求めたＭ3.8の地震波形を余震の波形とみなして、本震（Ｍ6.3）による加速度波形の再現を試みる。計算にあたり、本震の震源パラメータおよび断層モデルはTable.１のものを用いた。確率論的震源モデルによる地震波形の作成に当たっては、震源位置と3地震のマグニチュードは余震と同じとする。

小地震の方向性を考慮することはともかく、確率論的グリ3.3ーン関数法においても断層すべり方向と継続時間の関係を加速度時刻歴の計算例考慮することにより、観測波形の特徴をより精度よく再現

Figure.６に地震波形の比較を示す。経験的グリーン関数できるものと思われる。

法によるものは、最大加速度値、水平１方向のみの計測震
度：Ｉ、ならびに継続時間において、観測記録とよく対応４．おわりにしている。特に出水地点では継続時間が他の地点と比較して短いが、経験的グリーン関数法においても継続時間の短

本文では、サイト特性として金井・多治見スペクトルにい波形が得られている。経験的グリーン関数法に比べ、確よる工学的基盤面（VS＝250km/sec程度）での地震波を考慮

率論的グリーン関数法の継続時間はやや長めであるが、出したが、特に軟弱層による影響が大きい阿久根地点では、水地点の継続時間の特性がよく表現できている。

サイト特性を加味した増幅地震動が必要となる。別途解析

阿久根は、港湾地域での液状化現象が確認されており、等を実施することも考えられるが、地域防災計画や地震計阿久根地点の本震記録にみられる長周期成分の原因がそのの少ない地域の直後予測システム等では、当該再現地震波液状化のためと考えられがちであるが、余震記録（Acc19ga1）と一次元波動論等と組み合わせる程度で、従来の加速度のでも長周期成分が卓越していることを考えると（Figure.７）、経験式（アッテネーション式）に比べ精度よく算定が可能

観測地点の表層地盤が軟弱であることが推察される。よっであると思われる。また、経験的手法において余震データて、経験的グリーン関数法および確率論的グリーン関数法を用いる場合も、サイト（求める地震動位置）と断層の位とも観測波の比較および計測震度において、継続時間を考置関係を考慮し、方向が一致する位置での余震観測記録を慮することで再現性が向上すると思われる。

用いることが精度の高い予測地震動を得るためには重要とFigure.８に55の観測地点における最大加速度値と確率論思われる。

的グリーン関数法から求めた地震動の最大値を示す。断層［参考文献］
１ ） Irikura, K., Prediction of Strong Accelerations Motions using
Empirical Green's Function, Proc. of 7th Japan Earthquake
Engineering Symposium, pp.151-156, 1986.
２）Joyner, W. B. and Boore, D.M.,On Simulating Large Earthquakes by Green's Function Addition of Smaller Earthquakes,
Earthquake Source Mechanics, Maurice Ewing Series 6, edited

すべり方向を図中に示す。断層の上盤に相当する北側地域での加速度が大きいことが示され、断層の傾斜角が断層すべり方向と継続時間の関係13）を考慮することが重要であると思われる。また、経験的グリーン関数法においても用いる

4

統計的特性と地盤の非線形性を考慮した強震動予測、日本建築学会構造系論文報告集、第463号、pp.27-37、1994。
５）原田隆典、田中 剛、田村好文：震源特性を考慮した地震動のシミュレーション手法、土木学会論文報告集、第507号/I-30, pp.209-217、1995。
６）Shinozuka, M. and Sato, Y., Simu1ation of nonstationary
random process, Journal of the Engineering Mechanics Division, ASCE, Vol.93, No. EM1, pp.11-40, 1967.
７）大角恒雄、原田隆典、泉谷恭男：確率論的グリーン関数を用いた強震動の波形合成法とその検証例による震源パラメータの考察、
土木学会・応用力学論文集、Vol.1、1998.8
８）新・地震動のスペクトル解析入門、大崎順彦、鹿島出版会
９）Ohsumi, T., Harada, T., and Darama, H.（1997）, Engineering simulation of ground motions using a seismologicalmodel, Proc. of the 7th international Conference on StructuralSafety and Reliability, Kyoto, November 24-28.
10）Hisada, T. and Ando, H.：Relation between Duration ofEarthquake Motion and the Magnitude, Kajima Institute ofTechnology, June 1976.
11） Boore,D.M.,Stochastic Simulation of High Frequency
Ground Motion Based on Seismological Model of RadiatedSpectra, Bull. of Seism. Soc. of Am., Vo1.73, pp.1865-1894,1983.
12）Kiyono, J., Identification and Synthesis of Seismic Ground
Motion in Structural Response Analysis, PhD thesis, Dept. ofCivil Engineering, Kyoto University, 1992.
13）平成９年度建設技術講演会「鹿児島県北西部地震について」7鹿
児島県建設技術センター（1997.11.17）
14）大角恒雄、奥倉英世、原田隆典：地震断層を考慮した強振動の波
形合成、土木学会第52回年次学術講演会（1997.9）

by Das, S. and Boatwright, J., pp.269-273, 1986.
３）釜江克宏、入倉孝次郎、福知保長：地震のスケーリング則に基づいた大地震時の強震動予測、日本建築学会構造系論文報告集、第430号、pp.1-9、1991。
４）佐藤智美、川瀬 博、佐藤俊明：観測記録から同定した地震動の5

6