11:30 〜 11:45
[S22-08] 震央方位を考慮した機械学習による地点固有の地震動評価モデルの検討
新たな観点の地震動評価による新たな知見獲得を目指し、堆積平野内や周辺で得られた水平地震動観測記録を教師データとする機械学習により、地点固有の地震動評価モデルの作成を試みた。各地点では地盤特性が共通なので、震源・伝播特性を記述する各パラメータを入力する「特徴量」とし、観測記録から求められた地震動指標を「目的変数」とすることにより、両者を関連付けて機械学習を実施した。
機械学習には、勾配ブースティングと決定木を組み合わせた手法「勾配ブースティング木」を用いた。「目的変数」として、地震動の最大加速度PGA [cm/s2]および複数の周期T [s]での減衰定数h＝5%の擬似速度応答スペクトルPSV [cm/s]と速度応答継続時間スペクトルTSV [s]（パラメータp1＝0.03，p2＝0.95）を検討した。大振幅になるとデータ数が減るPGAとPSVは、モデル精度を高めるために機械学習に適した解析上のデータ分布となる常用対数（log10PGA，log10 PSV）に変換して入力した。「特徴量」としては、地震のモーメントマグニチュードMW、震源深さH [km]、震源距離X [km]、震央方位Λ [degree]を与えた。Λは真北を0°として時計回りに定めたが、真北を境に不連続量となるので、解析上のデータにはsinΛとcosΛのペアを与えた。一特徴量のデータ列のみシャッフルして置き換え機械学習により再評価した時に評価精度の悪化した度合いにより、モデルへのその特徴量の「影響度」とした。
K-NET観測点TKY028での「地震動評価モデルT」の検討例を図1～4に示す。Λの影響度は特にTSVで大きく、周期によってはMW・X・Hと同等以上になった。この例では、評価値／観測値の比の平均は約1.1、その常用対数標準偏差は振幅では0.2強、継続時間では0.1強で、評価値の大半は観測値の倍半分以内に収まった。従来の距離減衰式等で扱われなかった震央方位や応答継続時間を考慮し、地点・周期毎に地震動特性が変化する様子の定性的・定量的な分析・検討も可能で、地震波の三次元伝播特性の違いを簡便な評価モデルに反映出来る可能性がある。
AIとBig Dataを生かす上で良質・大量のデータと情報の得られる地点毎に地震動評価モデルを作成することには大きな利点がある。今後、データの粗密に起因するモデル精度のアンバランスを改良するための工夫や、データ数が少なく外挿領域に当たる大地震・大振幅地震動・長時間地震動の評価上の工夫も必要となろう。
機械学習には、勾配ブースティングと決定木を組み合わせた手法「勾配ブースティング木」を用いた。「目的変数」として、地震動の最大加速度PGA [cm/s2]および複数の周期T [s]での減衰定数h＝5%の擬似速度応答スペクトルPSV [cm/s]と速度応答継続時間スペクトルTSV [s]（パラメータp1＝0.03，p2＝0.95）を検討した。大振幅になるとデータ数が減るPGAとPSVは、モデル精度を高めるために機械学習に適した解析上のデータ分布となる常用対数（log10PGA，log10 PSV）に変換して入力した。「特徴量」としては、地震のモーメントマグニチュードMW、震源深さH [km]、震源距離X [km]、震央方位Λ [degree]を与えた。Λは真北を0°として時計回りに定めたが、真北を境に不連続量となるので、解析上のデータにはsinΛとcosΛのペアを与えた。一特徴量のデータ列のみシャッフルして置き換え機械学習により再評価した時に評価精度の悪化した度合いにより、モデルへのその特徴量の「影響度」とした。
K-NET観測点TKY028での「地震動評価モデルT」の検討例を図1～4に示す。Λの影響度は特にTSVで大きく、周期によってはMW・X・Hと同等以上になった。この例では、評価値／観測値の比の平均は約1.1、その常用対数標準偏差は振幅では0.2強、継続時間では0.1強で、評価値の大半は観測値の倍半分以内に収まった。従来の距離減衰式等で扱われなかった震央方位や応答継続時間を考慮し、地点・周期毎に地震動特性が変化する様子の定性的・定量的な分析・検討も可能で、地震波の三次元伝播特性の違いを簡便な評価モデルに反映出来る可能性がある。
AIとBig Dataを生かす上で良質・大量のデータと情報の得られる地点毎に地震動評価モデルを作成することには大きな利点がある。今後、データの粗密に起因するモデル精度のアンバランスを改良するための工夫や、データ数が少なく外挿領域に当たる大地震・大振幅地震動・長時間地震動の評価上の工夫も必要となろう。