データ同化を用いた巨大地震後の余効すべりと粘弾性変形の同時推定

利用数据同化同时估计大地震后的后效滑移和粘弹性变形

• 批准号: 20K14574
• 负责人: 大谷 真紀子
• 金额: $ 2.66万
• 依托单位: Kyoto University (2022)The University of Tokyo (2020-2021)
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Early-Career Scientists
• 财政年份: 2020
• 资助国家: 日本
• 起止时间: 2020-04-01 至 2024-03-31
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/en/grant/KAKENHI-PROJECT-20K14574/
• 关键词: データ同化; 断層すべり; アジョイント法; EnKF; 余効変動; 粘弾性変形; 余効すべり; 地殻変動; 粘弾性

本研究では、巨大地震発生後に観測される余効変動データから、断層面上で起こる余効すべりと地下深部の粘弾性領域で起こる粘弾性変動の両者を切り分け、変動の推移及び物性パラメタを推定することを目的とした手法開発を行う。推定手法にはデータと数値モデルの統合を行うデータ同化手法を用いる。データ同化手法には大きく分けて、逐次データ同化と変分法の二つの手法が存在する。断層すべりのような非線形モデルに用いられる逐次データ同化手法としてEnKFが挙げられるが、これは四次元変分法に比べて導入が容易であるため、本研究ではこれまでに巨大地震の発生による応力擾乱によって駆動される地下の粘弾性変形問題において逐次データ同化手法であるEnKFの適用可能性を調べ、その有効性を数値実験で確認した。一方、四次元変分法であるアジョイント法にはEnKFよりも少ない計算量で行えるなどの利点がある。非線形モデルにおいては、これらEnKF, アジョイント法と用いる手法が違えば、推定のされ方や推定結果も異なってくる。そこで、これら手法による推定結果の違いを調べるために、昨年度は断層すべり問題として既にEnKFを適用した先行研究が存在する、繰り返すゆっくりすべり(SSE)問題に対してアジョイント法の適用を行い、数値実験を実施した。一周期分のSSE地殻変動データに対して、断層の初期摩擦強度分布を更新しながらアジョイント法による摩擦パラメタ推定を繰り返し行う二段階反復法を新たに考案し、断層強度分布及びパラメタを推定することに成功した。昨年度は周期的SSE１周期分の人工データを用いた結果を示したが、本年度SSE半周期分のデータを用いた際にも推定が上手くいくことを示した。また次のSSEの予測性を数値実験によって調べたところ、本手法がSSEの繰り返し間隔を良く推定し、次のSSE発生時期の予測性能はEnKFを上回ることが分かった。

This study で は, huge earthquake 発 postnatal に 観 measuring さ れ る yu sharper - move デ ー タ か ら, fault surface で こ る yu unseen す べ り と deep underground の stick で 弾 sex field case こ る sticky 弾 sex - move の struck the を り け, cutting - move の passes and び property パ ラ メ タ を presumption す る こ と を purpose と し た technique open 発 う を line. Presumption of gimmick に は デ ー タ と the numerical モ デ ル の integration を line う デ ー タ assimilation technique を with い る. デ ー タ assimilation technique に は big き く points け て, successive デ ー タ assimilation と の two methods - つ の gimmick が exist す る. Fault す べ り の よ う な nonlinear モ デ ル に with い ら れ る successive デ ー タ assimilation technique と し て EnKF が 挙 げ ら れ る が, こ れ は に four yuan - points method than べ て import が easy で あ る た め, this study で は こ れ ま で に huge earthquake の 発 raw に よ る 応 force disrupt に よ っ て 駆 dynamic さ れ る underground の sticky 弾 sex - shaped に お い て successive デ ー The タ assimilation method であるEnKF <s:1> application possibility を the べ and そ are effective を numerical experiments で confirm the た. The one-dimensional and four-dimensional variation method であるアジョ the であるアジョ ト the ト ト method に the <s:1> EnKFよ the <s:1> less な the な computational load で rows えるな <s:1> the <s:1> illy point がある. Nonlinear モ デ ル に お い て は, こ れ ら EnKF, ア ジ ョ イ ン と ト method with い る gimmick が violations え ば presumption, の さ れ party や constructive results も different な っ て く る. そ こ で, こ れ ら gimmick に よ る constructive results の violations い を adjustable べ る た め に, yesterday's annual は fault す べ り problem と し て both に EnKF を applicable し た leading research が exist す る, Qiao り return す ゆ っ く り す べ り (SSE) problem に し seaborne て ア ジ ョ イ ン ト の applicable を い, the numerical be 験 を be applied し た. A periodic points の SSE crust - move デ ー タ に し seaborne て intensity distribution, fault early の friction を update し な が ら ア ジ ョ イ ン ト method に よ る friction パ ラ メ タ presumption を Qiao り return line し う two Duan Jie anadiplosis を new た し に test case, fault intensity distribution and び パ ラ メ タ を presumption す る こ と に successful し た. Yesterday annual は cycle SSE1 cycle points の artificial デ ー タ を with い を た results shown し た が, SSE this year and a half cycle の デ ー タ を with い た interstate に も presumption が overhand く い く こ と を shown し た. ま た times の SSE sex を の to test the numerical be 験 に よ っ て adjustable べ た と こ ろ, this technique が SSE の Qiao り presumption to return good し interval を く し, の SSE 発 の born period to measure performance は EnKF を last る こ と が points か っ た.

项目成果

巨大地震発生機構の理解と予測可能性に関する地震発生サイクルシミュレーション

地震发生周期模拟，用于了解和预测特大地震发生机制

• 发表时间: 2021
• 作者: Baba;S.;Obara;K.;Takemura;S.;Takeo;A.;& Abers;G. A.;三好雅也・新宮竜太朗;大谷真紀子
• 通讯作者: 大谷真紀子

豊後SSEの断層摩擦特性・すべり推定の4DVAR数値実験

用于估算Bungo SSE断层摩擦特性和滑移的4DVAR数值实验

• 发表时间: 2022
• 作者: Escobar Ma. Teresa;Takahata Naoto;Obata Hajime;Sano Yuji;大谷真紀子・亀伸樹・加納将行
• 通讯作者: 大谷真紀子・亀伸樹・加納将行