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

地震波三维数值模拟微分算子优化研究

• 批准号: 06740351
• 负责人: 古村 孝志
• 金额: $ 0.77万
• 依托单位: Hokkaido University of Education
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Encouragement of Young Scientists (A)
• 财政年份: 1994
• 资助国家: 日本
• 起止时间: 1994 至 无数据
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-06740351/
• 关键词: 地震; 数値シミュレーション; 並列計算; Pseudospectral法; 数値微分

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

Pseudospectral method is used to calculate the seismic ratio, and the new Fourier differential method is developed. The Fourier differential method is used to calculate the Nyquist wave number and the total wave number. The calculation of the cumulative ratio of the calculation field is based on the Pseudpspectral method. In this study, the Fourier derivative was investigated by using the visualization of the computational process. However, during this period of time, the effectiveness of Lanchos 'utilization and calculation process in reducing degradation has been carefully combined. The Fourier differential method (Symmetry method, Furumura·Takenaka, 1992) is generally used to solve the problem of high wave numbers. Second, the possibility of calculating the three-dimensional ratio at high speed and parallel calculation is investigated. In the existing computer environment, it is difficult to calculate the three-dimensional space, and it is difficult to calculate the three-dimensional space. The best way to achieve this goal is to visualize and optimize the communication between servers. Pseudospectral method: Fourier differentiation is independent of calculation, and calculation is complex and decentralized. The number of divisions is 4, and the transmission rate is 4 KByte (Ethernet). The calculation efficiency is maximum.

项目成果

古村孝志: "Pseudospectral法による地震波動場の3次元計算-現状と課題" 東京大学地震研究所一般共同研究「不均質な地下構造における地震伝達のシミュレーション」報告書. (印刷中). (1995)

Takashi Furumura：《利用伪谱法进行地震波场的三维计算——现状与挑战》东京大学地震研究所联合研究《非均质地下结构中地震传播的模拟》报告（正在出版）。 (1995)