WIHMコマンドで用いている計算式

(Formula used in WIHM command)


1. 共通の記号 (Common notations)



2. シングルフォースが作る変位波形 (Displacement waveforms generated by a single force)

Aki and Richards (2002)の(4.27)式を使用する。 シングルフォースの時間関数を \(\myvector{F}(t)=(F_1(t),F_2(t),F_3(t))^T\) とおくと、変位波形は以下の式で与えられる。
Eq. (4.27) of Aki and Richards (2002) is used. Let \(\myvector{F}(t)=(F_1(t),F_2(t),F_3(t))^T\) be the time function of the single force. Then the displacement waveform is given by the formula below.

\[\begin{equation} u_n(t)=u_n^{N}(t)+u_n^{FP}(t)+u_n^{FS}(t) \hspace{1em} (n=1,2,3) \label{eq.F} \end{equation}\]


3. モーメントテンソルが作る変位波形 (Displacement waveforms generated by a moment tensor)

Aki and Richards (2002)の(4.29)式を使用する。 モーメントテンソルの時間関数を \[\begin{equation} \myvector{M}(t)= \begin{pmatrix} M_{11}(t) & M_{12}(t) & M_{13}(t) \\ M_{21}(t) & M_{22}(t) & M_{23}(t) \\ M_{31}(t) & M_{32}(t) & M_{33}(t) \end{pmatrix} \label{eq.moment} \end{equation}\] とおくと、変位波形は以下の式で与えられる。
Eq. (4.29) of Aki and Richards (2002) is used. Let \(\myvector{M}(t)\) (eq. \ref{eq.moment}) be the time function of the moment tensor. Then the displacement waveform is given by the formula below.

\[\begin{equation} u_n(t)=u_n^{N}(t)+u_n^{IP}(t)+u_n^{IS}(t)+u_n^{FP}(t)+u_n^{FS}(t) \hspace{1em} (n=1,2,3) \label{eq.M} \end{equation}\]


4. 計算方法 (Computation method)

WIHMコマンドでは\(F(t)\)や\(M(t)\)の時間関数を解析式で与える。 また、それらの時間関数の導関数についても解析式が利用できる。 したがって中間項や遠地項はすべて解析式を用いて 時刻\(t\)での値を計算できる。
Time functions of \(F(t)\) and \(M(t)\) are given by analytical formulas in WIHM command. Analytical formulas are also available for the derivatives of these time functions. Therefore, the values of all intermediate- and far-field terms at time \(t\) can be computed analytically.

WIHMコマンドでは数値積分による誤差を避けるため、 近地項についても解析式を導出してそれを用いている。 \(G(t)\)を2階微分可能な関数とし、 その1階導関数を\(g_1(t)\equiv\dot{G}(t)\)、 2階導関数を\(g_2(t)\equiv\dot{g_1}(t)=\ddot{G}(t)\)とする。 \(a\), \(b\)を定数として \[\begin{equation} I(g_2;a,b;t) \equiv \int_a^b \tau g_2(t-\tau) d\tau \label{eq.I.definition} \end{equation}\] を考える。 この積分は次のように変形できる。 \[\begin{eqnarray} I(g_2;a,b;t) &=& \int_a^b \tau g_2(t-\tau) d\tau \nonumber \\ &=& \left[-\tau g_1(t-\tau)\right]_a^b -\int_a^b -g_1(t-\tau) d\tau \nonumber \\ &=& \left[\left\{-b g_1(t-b)\right\}-\left\{-a g_1(t-a)\right\}\right] +\int_a^b g_1(t-\tau) d\tau \nonumber \\ &=& -\left[b g_1(t-b) - a g_1(t-a)\right] +\left[-G(t-\tau)\right]_a^b \nonumber \\ &=& -\left[b g_1(t-b) - a g_1(t-a)\right] -\left[G(t-b)-G(t-a)\right] \label{eq.I.analytical} \end{eqnarray}\] (\ref{eq.F.N})式は\(I(F;r/\alpha,r/\beta;t)\)、 (\ref{eq.M.N})式は\(I(M;r/\alpha,r/\beta;t)\) である。 したがって震源時間関数(\(F(t)\)や\(M(t)\))を (\ref{eq.I.analytical})式の\(g_2(t)\)と見なすことができ、 \(g_1(t)\)は震源時間関数の1階積分、 \(G(t)\)は震源時間関数の2階積分である。 ymaeda_opentoolsでは時間関数の1階積分や2階積分の解析式も利用できるので、 (\ref{eq.I.analytical})式を用いれば近地項も解析的に計算できる。
To avoid errors caused by numerical integration, the near-field term is also composed analytically in WIHM command. Let \(G(t)\) be a function that can be differentiated twice, and \(g_1(t)\equiv\dot{G}(t)\) and \(g_2(t)\equiv\dot{g_1}(t)=\ddot{G}(t)\) be the 1st- and 2nd-order derivatives of \(G(t)\), respectively. Let \(a\) and \(b\) be constants, and an integral \(I(g_2;a,b;t)\) be defined by eq. (\ref{eq.I.definition}). This integral can be arranged as (\ref{eq.I.analytical}). Note that eqs. (\ref{eq.F.N}) and (\ref{eq.M.N}) are equivalent to \(I(F;r/\alpha,r/\beta;t)\) and \(I(M;r/\alpha,r/\beta;t)\), respectively. Therefore the source time functions (\(F(t)\) and \(M(t)\)) can be regarded as \(g_2(t)\) in eq. (\ref{eq.I.analytical}); \(g_1(t)\) and \(G(t)\) are 1st- and 2nd-order integrals of the source time function, respectively. Analytical formulas for the 1st- and 2nd-order integrals of time functions are available in ymaeda_opentools, meaning that the near-field term can be computed analytically using eq. (\ref{eq.I.analytical}).