地震波動の3次元数値シミュレーションのための微分演算子の最適化に関する研究

1994 Fiscal Year Annual Research Report

• Project/Area Number: 06740351
• Research Institution: Hokkaido University of Education
• Principal Investigator: 古村 孝志 北海道教育大学, 教育学部・岩見沢校, 助教授 (80241404)
• Keywords: 地震 / 数値シミュレーション / 並列計算 / Pseudospectral法 / 数値微分

Research Abstract

Pseudospectral法を用いて地震波動計算を行う際に発生する高周波ノイズを軽減するための,新しい数値微分(Fourier微分)法を開発した。従来のFourier微分法ではNyquist波数までの全ての波数を厳密に扱うため,計算を進めるにつれて高波数ノイズが蓄積し,計算が発散するという問題があった。このような計算ノイズの累積は計算領域の次元が高いほど顕著に現れるため,これまでのPseudpspectral法による波動計算では3次元波動計算ができなかった。本研究では,まず,Fourier微分でノイズが発生するメカニズムを,計算過程の可視化ビデオ画像をもとに調査した。そして,このノイズ軽減のために,Lanchosフイルタの利用と計算の過程に減衰の効果を組み込む工夫の二つについて検討した。また、Pseusoepectral法で一般に用いられている,不連続で安定なFourier微分法(対称法,古村・竹中,1992)が高波数ノイズを地表付近にトラップさせるという問題も明らかにした。次に,3次元波動計算を高速に行うための並列計算の可能性について調査した。現有の計算機環境では実用的な3次元計算が困難であるため,計算は8台のワークステーションクラスタを用いた疎結合の並列システム上で行った。プロセッサ間のデータ通信のようすをビデオ装置を用いて可視化し,並列化の効果が最大となる最適な分散形態を求めた。Pseudospectral法ではFourier微分が独立した計算の基本単位となるが,この演算を複数のワークステーションで分散処理させた場合にはデータ転送と回収の時間が全計算のボトルネックとなる。このため,分割数が4つ,転送データサイズは4KByte程度(Ethernetの場合)に計算の効率が最大となることがわかった。

Research Products

• [Publications] 古村孝志,竹中博士: "Pseudospectral法による理論地震記象-高周波ノイズ軽減のための微分演算子の最適化-" 物理探査学会第91回学術講演会論文集. 91. 277-280 (1994)
• [Publications] Furumura.T.and H.Takenaka: "A wroparound elimination technigue for the pseudospectral wave synthesis using on antiperiodic extension of the wavefield." Geophysics. 59. 1-9 (1995)
• [Publications] 古村孝志: "Pseudospectral法による地震波動場の3次元計算-現状と課題" 東京大学地震研究所一般共同研究「不均質な地下構造における地震伝達のシミュレーション」報告書. (印刷中). (1995)