はじめに

地震ハザード評価のための広帯域地震動計算（概ね0.05〜20秒程度）
長周期（低周波数：LF）側と短周期（高周波数：HF）側で別々に計算してから、接続周波数で足し合わせる手法がとられることが多い。
＝ハイブリッド合成法

長周期地震動評価の広帯域化に向けた検討¾ 低周波数側 (LF)：理論的・決定論的手法 FDM, FEM, etc.
•

震源過程と3次元的な波動伝播過程（サイト効果含む）を考慮できる•

震源・地下構造モデルの精度、計算コストэ

せいぜい1 Hz以下

¾ 高周波数側 (HF)：半経験的/統計的手法

岩城麻子（防災科学技術研究所）•

統計的な地震動特性をもつ波形ができる•

特徴的な地震動の特性が反映されず、LF地震動と整合的でないことがある

1

はじめに

2

本研究の目的¾ 理論計算による低周波数側（LF）の地震動評価の高度化

約1秒程度まで低周波数地震動が計算されていることを前提とし、それを使って高周波数地震動を合成する方法（岩城・藤原, 2013）• 地下構造モデル：J-SHIS、全国一次地下構造モデル、浅部を含む高精度化• 震源モデル：知見の蓄積、不確実性の考慮

LF ( f < ~ 1 Hz ) FDM• 計算機環境：多数のシナリオ計算に基づくハザード評価が現実的になりつつある

HF ( f > ~ 1 Hz )

100
10-1
地下構造モデル・震源モデル
長周期地震動予測地図作成等支援事業(2012,2013)101

Freq. (Hz)•

観測記録を用いるので観測と整合的な地震動が得られる•

短波長の地盤モデル、震源モデルを特には必要としない•（HFの計算には）計算コストかからないэ新しい広帯域地震動予測手法として提案することを目指す
34

茨城県沖MW 5.9
K-NET つくば ǻaNP

加速度エンベロープに着目

提案手法の概略ɶᙹ೉‫ע‬ᩗƷᚇยᚡ᥵Ɣǒኺ᬴ႎƳ᧙̞ࡸਁЈ¾ 異なる周波数帯間の加速度エンベ
ロープに着目すると，S波立ち上がり部分までは互いに似た形状をしている

Obs.
~1 Hz

Obs.
2-4 HzObs.
1-2 HzObs.
4-8 HzObs.
8-16 HzLF¾ 形状の違いは主として周波数帯間伝播特性とサイト特性の違いにの伝起因していると考えられる差分法等で計算したターゲット地震の低周波数エンベロープ~1.0 Hz

LF
周波数帯間の伝播・サイト特性の経験的な関係性を地点ごとに評価

エンベロープ比関数

HF

エンベロープ比関数ǿȸǲȃȈ‫ע‬ᩗ

LF

FDM
~1 Hz

HF
低周波数(LF)地震動から

Ǩȳșȭȸȗൔ᧙ૠƱˮႻǛɨƑƯӳ঺Syn.
2-4 HzSyn.
1-2 HzSyn.
4-8 HzSyn.
8-16 Hz

LF

高周波数(HF)地震動を合成

広帯域予測地震動

HF

65

エンベロープ比関数 Hn,n-1
S波立ち上がり部以降のエンベロープを強強震部と減衰部に分けて以下の様にモデル化できると仮定

ࢍᩗᢿan-1bn-1

大地震のHF加速度波形の合成˯ԗඬૠ中小地震観測記録から得られたパラメータ

強震部：振幅一定減衰部：‫ି ܴ ݂ ܩ ݂ ܵ = )ݐ ;݂(ܧ‬ఊ ‫)ݐ ;݂(ܣ‬En-1ᩗเǵǤȈ࠹˴ถᘛϋᢿถᘛan⇛⇊⇮‫ࠢف‬ཎࣱ‒ᩗเ⇟∀⇕⇮∑ൔ⁅⁠‒‡⁅⁠‟‣‒ࢍᩗዒዓ଺᧓↝ࠀ‒˯ԗඬૠ↝ࢍᩗዒዓ଺᧓‒ถᘛཎࣱ‒ߨ݂‫ݐ‬‫ܣ = )ݐ ;݂(ܣ‬଴ exp െ
ܳ(݂)

bn

ターゲット地震（大地震）の情報⇎∙⇿∓∞⇽ൔ᧙ૠ›⁠„⁠‟‣‒

᭗ԗඬૠ

En

周波数帯域n(n=0, 1, …, 4)
強震継続時間bnとする‫ܧ‬௡ ‫= ݐ‬
n=0 : 0.5 – 1 Hz
n=1 : 1 – 2 Hz
n=2 : 2 – 4 Hz
n=3 : 4 – 8 Hz
n=4 : 8 – 16 Hzܽ௡ (‫ ݐ‬൑ ܾ௡ )ܽ௡ exp െܿ௡ ‫ ݐ‬െ ܾ௡ )(‫ܾ > ݐ‬௡ )位相スペクトル
ターゲット地震の低周波数加速度エンベロープ

LFとHFのエンベロープの比‫ܪ‬௡,௡ିଵ ‫ܧ ؠ‬௡ /‫ܧ‬௡ିଵ
782003年十勝沖地震(MW8.3)への適用

3次元差分法による低周波数(LF)地震動

1. 震源モデルと地下構造モデル
Aoi et al. (2008, JGR) に基づき、３次元差分法で低低周波数(LF)地震動計算
2. 20個の余震(5.5 < MW < 6.6)の地表強震観測記録を用いてエンベロープ比関数を作成•

計算は3次元差分法による地震動シミュ
レータGMS（青井・他, 2004）による•

Aoi et al. (2008)の震源・地下構造モデル
は0.3 Hz以下を対象としているが、ここでは0.5 Hzまでの計算波形を使用•

観測速度波形(0.05 - 0.5 Hz)を概ね良く再現できている•

一部高周波数側で振幅が足りない観測点もある

3. 高周波数(HF)地震動の合成• すべりの大きい領域（アスペリティA1, A2, A3）とその他の領域（背景領域）に分けてLF地震動を計算• それぞれを別々に広帯域化し、最後に波形を足し合わせる

9高周波数(HF)地震動の合成

10

広帯域地震動（観測・合成）

広帯域(0.05-16 Hz)速度波形・スペクトルレベルは観測を概ねよく再現している¾ 余震観測記録からエンベロープ比関数を作成• S波到着時刻から120秒間を対象• メカニズムの補正はしていない• 0.5-1, 1-2, 2-4, 4-8, 8-16 Hz‾‸⇎∙⇿∓∞⇽‒• ターゲット地震（本震）の各アス
ペリティ・背景領域ごとに、円形クラック断層におけるȦ-2のコーナー周波数を仮定
ĺ震源スペクトル比 ܵ௡ /ܵ௡ିଵ

¾ HF波
波形の合成
• LF加速度エンベロープにエンベロープ比関数を順次掛け合わせる• 位相：0.5-0.6Hzの間でなめらかにつながる乱数を選択

1112

広帯域地震動（観測・合成）

統計的グリーン関数法(SGM)
によるハイブリッド法との比較

 低周波数(LF)地震動で再現性の良くない観測点

эHFでも良くない LFでの再現性は悪くなかったが 、HFの波を再現できていない⇒ サイト特性や伝播経路特性の影響は線形の範囲では経験式で考慮されているはずなので、震源特性の影響と考えられる。

13

まとめと今後の課題¾ 周波数帯域間の加速度エンベロープ特性の関係性を利用した高周波数地震動合成手法を、M8クラスの地震に適用する方法を検討した¾ 3次元差分法による低周波数地震動シミュレーション（f<0.5Hz）と組み合わせて、2003年十勝沖地震の広帯域地震動シミュレーションを行った。
•

差分法で0.5Hzまでよく再現されている観測点では、広帯域(<16Hz)でも波形を概ねよく再現できた。

•

一部観測点では低周波数帯域では見られないパルスが高周波数帯域では観測されており、それを再現することができなかった。
⇒ 周波数依存性を持つ震源モデルには対応できていない¾ 現時点での手法の適用範囲•震源距離：数十km以上〜200km程度以内•

地震規模：M7, M8クラス•

ターゲット地震の近くで中小地震記録が得られている¾ 今後の課題•

予測問題におけるLF地震動の不確実性•

震源距離の短い観測点での地震動•（周波数依存性を持つ震源モデル）（＋パラメータのバラツキ）

1514