09:15 〜 09:30
[S01-11] 人工震源装置「アクロス」を用いた振幅変化検出手法の開発と実データへの適用
我々は、伝播振幅の変化を人工震源装置「アクロス」を用いて取得した伝達関数から検出する手法を開発し、その手法を先行研究Ikuta et al. (2002)で取得された伝達関数に適用した。その結果、地震時の強振動や注水実験に伴う振幅変化等を含めて、ノイズの影響を考慮して連続的に振幅変化を推定できることを確認した。
減衰の時間変化は、間隙流体の移動やクラックの開閉といった地下構造の変化を反映すると考えられる。従来は自然地震のスペクトルの傾きを比較する手法で主に調べられ、地震の発生に伴う減衰の増加(Chun et al. 2004; Kelly et al. 2013)や、火山活動に伴う減衰の時間変化(Titzschkau et al. 2010)が報告されている。
本研究は、このような減衰の時間変化を、人工震源装置「アクロス」の特長を利用して推定する新手法を開発した。本手法の特徴は、ノイズの変化に影響されずに減衰変化を推定できることにある。周波数領域ではアクロスの信号は飛び飛びの周波数成分で表現できるため、隣接する周波数成分の分散はその周波数の伝達関数に含まれているノイズの分散と同一であると仮定することができる。信号の周波数成分の２乗からノイズ分散を差し引いたものが、求める信号の２乗の推定値となる。これにより、エンベロープの比を計算する手法や、単純にエネルギーの比を計算する手法に比べてノイズの影響を受けずに推定することが可能になった。
新たに開発した手法が実データに適用できることを確認するためIkuta et al. (2002)で取得され、Ikuta and Yamaoka (2004)で走時変化に基づいて発生した現象の詳細な検討が行われている伝達関数に適用した。この伝達関数は、2000年1月から2001年4月にかけて、野島断層の分岐断層近傍に設置されたテストサイトで得られた。この期間には、比較的規模が大きな地動を伴う地震が4回発生しており、また、断層に水を圧入する注水実験が合計4回行われた。この伝達関数は、走時変化に関して詳細に解析がなされており、地震時の変化に関しては、Ikuta et al. (2002)で走時変化の解析が行われ、2000年鳥取県西部地震(MW=6.6, 震央距離165km)と、2001年芸予地震(MW=6.4, 震央距離215km)の2回の地震に関して走時変化が起こったことが、注水実験に関しては、山岡・生田(2001)で4回行われた実験のうち、2回目と4回目の注水に伴う走時変化が報告されている。さらに、地震時の変化の原因については、Ikuta and Yamaoka (2004)で異方性に注目した詳細な検討が行われ、間隙流体圧の上昇によるクラックの選択的な開口が原因であったと推測されている。
本研究の手法を適用して得られた振幅変化の結果には、地震時の振幅変化と注水実験に伴う変化、および降雨に伴う変化や長期的な変化が見られた。振幅変化は、鉛直成分のP波と、水平2成分のS波に着目して計算を行い、同時に、先行研究の手法で計算した走時変化も計算した。求めた振幅変化と走時変化を図に示す。地震時の振幅変化は、2000年鳥取県西部地震に関しては、最大で5%にも及ぶ振幅の急激な減少と、その後約1週間程度で振幅がもとの水準に回復する様子が見られたが、芸予地震に関してはわずかなステップ上の振幅変化が見られたのみであった。注水実験に伴う変化は、最大で5%程度で、4回行われた注水実験のうち、2回目と4回目のみに顕著な変化が観測された。その変化は、どちらの回でも一度振幅が上昇し、その後元に戻るような変化であった。
Ikuta and Yamaoka (2004)で指摘されているように、地震時の変化がクラックへの流体の流入によって起こったならば、粘性損失によって振幅の減少が期待される。我々の解析でも、地震時には振幅が減少し、彼らの考えと矛盾しない。一方、注水実験に伴う変化は、振幅の増大と、それに引き続く振幅減少とみなすことができる。振幅の増加は注水による応力の増大によって周辺岩石中のクラックが閉じることが原因、振幅の減少は注水された水が拡散することで粘性損失の範囲が拡大することが原因と解釈できる。
減衰の時間変化は、間隙流体の移動やクラックの開閉といった地下構造の変化を反映すると考えられる。従来は自然地震のスペクトルの傾きを比較する手法で主に調べられ、地震の発生に伴う減衰の増加(Chun et al. 2004; Kelly et al. 2013)や、火山活動に伴う減衰の時間変化(Titzschkau et al. 2010)が報告されている。
本研究は、このような減衰の時間変化を、人工震源装置「アクロス」の特長を利用して推定する新手法を開発した。本手法の特徴は、ノイズの変化に影響されずに減衰変化を推定できることにある。周波数領域ではアクロスの信号は飛び飛びの周波数成分で表現できるため、隣接する周波数成分の分散はその周波数の伝達関数に含まれているノイズの分散と同一であると仮定することができる。信号の周波数成分の２乗からノイズ分散を差し引いたものが、求める信号の２乗の推定値となる。これにより、エンベロープの比を計算する手法や、単純にエネルギーの比を計算する手法に比べてノイズの影響を受けずに推定することが可能になった。
新たに開発した手法が実データに適用できることを確認するためIkuta et al. (2002)で取得され、Ikuta and Yamaoka (2004)で走時変化に基づいて発生した現象の詳細な検討が行われている伝達関数に適用した。この伝達関数は、2000年1月から2001年4月にかけて、野島断層の分岐断層近傍に設置されたテストサイトで得られた。この期間には、比較的規模が大きな地動を伴う地震が4回発生しており、また、断層に水を圧入する注水実験が合計4回行われた。この伝達関数は、走時変化に関して詳細に解析がなされており、地震時の変化に関しては、Ikuta et al. (2002)で走時変化の解析が行われ、2000年鳥取県西部地震(MW=6.6, 震央距離165km)と、2001年芸予地震(MW=6.4, 震央距離215km)の2回の地震に関して走時変化が起こったことが、注水実験に関しては、山岡・生田(2001)で4回行われた実験のうち、2回目と4回目の注水に伴う走時変化が報告されている。さらに、地震時の変化の原因については、Ikuta and Yamaoka (2004)で異方性に注目した詳細な検討が行われ、間隙流体圧の上昇によるクラックの選択的な開口が原因であったと推測されている。
本研究の手法を適用して得られた振幅変化の結果には、地震時の振幅変化と注水実験に伴う変化、および降雨に伴う変化や長期的な変化が見られた。振幅変化は、鉛直成分のP波と、水平2成分のS波に着目して計算を行い、同時に、先行研究の手法で計算した走時変化も計算した。求めた振幅変化と走時変化を図に示す。地震時の振幅変化は、2000年鳥取県西部地震に関しては、最大で5%にも及ぶ振幅の急激な減少と、その後約1週間程度で振幅がもとの水準に回復する様子が見られたが、芸予地震に関してはわずかなステップ上の振幅変化が見られたのみであった。注水実験に伴う変化は、最大で5%程度で、4回行われた注水実験のうち、2回目と4回目のみに顕著な変化が観測された。その変化は、どちらの回でも一度振幅が上昇し、その後元に戻るような変化であった。
Ikuta and Yamaoka (2004)で指摘されているように、地震時の変化がクラックへの流体の流入によって起こったならば、粘性損失によって振幅の減少が期待される。我々の解析でも、地震時には振幅が減少し、彼らの考えと矛盾しない。一方、注水実験に伴う変化は、振幅の増大と、それに引き続く振幅減少とみなすことができる。振幅の増加は注水による応力の増大によって周辺岩石中のクラックが閉じることが原因、振幅の減少は注水された水が拡散することで粘性損失の範囲が拡大することが原因と解釈できる。