waterPMLコマンドにおける地下構造の設定ファイルの書式(layer)

(A structure file format (layer) for waterPML command)


◆概要 (Overview)


◆書式 (Format)


◆例 (Examples)

以下の例において[TAB]はタブを表す。 またP波速度を\(V_p\)、S波速度を\(V_s\)、密度を\(\rho\)と表す。
In the following examples, [TAB] represents a tab. The P-wave velocity, S-wave velocity, and density are expressed as \(V_p\), \(V_s\), and \(\rho\), respectively.

(1)均質媒質 (A homogeneous medium)

物理領域の\(z\)座標範囲を\(-2000\leq z\leq 1000\) [m]とし、 その全域において \(V_p=3000\) [m/s], \(V_s=V_p/\sqrt{3}=1732.1\) [m/s], \(\rho=2300\) [kg/m\(^3\)]の均質媒質とする。 この地下構造は以下のexample1.iniによって実現できる。 ここで地下構造の定義域を物理領域よりも少し広く \(-2010\leq z\leq 1010\) [m]の範囲で取った。
Let the \(z\)-coordinate range of the physical volume be \(-2000\leq z\leq 1000\) [m], and let the medium be homogeneous with \(V_p=3000\) [m/s], \(V_s=V_p/\sqrt{3}=1732.1\) [m/s], and \(\rho=2300\) [kg/m\(^3\)]. This structure is realized by example1.ini below, where the definition range of the subsurface structure was taken as \(-2010\leq z\leq 1010\) [m], slightly wider than the physical volume.

(example1.ini)
1010.0[TAB]3000.0[TAB]1732.1[TAB]2300.0
-2010.0[TAB]3000.0[TAB]1732.1[TAB]2300.0


(2)標高の区分線形関数で与えられる構造 (A structure given as piecewise linear functions of the altitude)

物理領域の\(z\)座標範囲を\(-10000 \leq z\leq 0\) [m]とし、 この範囲でのP波速度、S波速度、密度が標高\(z\)の区分線形関数として 表1のように与えられる場合を考える。 この地下構造は以下のexample2.iniによって実現できる。 ここで地下構造の定義域を物理領域よりも少し広く \(-10000.1 \leq z\leq 0.1\) [m]の範囲で取った。
Let the \(z\)-coordinate range of the physical volume be \(-10000\leq z\leq 0\) [m], and let the P-wave velocity, S-wave velocity, and density be piecewise linear functions of the altitude \(z\) given as Table 1. This structure is realized by example2.ini below, where the definition range of the subsurface structure was taken as \(-10000.1 \leq z\leq 0.1\) [m], slightly wider than the physical volume.

表1(Table 1)
\(z\) [m] \(V_p\) [m/s] \(V_s\) [m/s] \(\rho\) [kg/m\(^3\)]
0 2000 1000 2100
-500 2500 1500 2200
-2000 3500 2000 2400
-5000 5000 3000 2600
-10000 7000 4000 2800

(example2.ini)
0.1[TAB]2000.0[TAB]1000.0[TAB]2100.0
-500.0[TAB]2500.0[TAB]1500.0[TAB]2200.0
-2000.0[TAB]3500.0[TAB]2000.0[TAB]2400.0
-5000.0[TAB]5000.0[TAB]3000.0[TAB]2600.0
-10000.1[TAB]7000.0[TAB]4000.0[TAB]2800.0


(3)不連続な水平成層構造 (A discrete horizontally stratified medium)

物理領域の\(z\)座標範囲を\(-10000 \leq z\leq 0\) [m]とし、 この範囲でのP波速度、S波速度、密度が 表2のように与えられる場合を考える。 この地下構造は以下のexample3.iniによって実現できる。 この例では速度・密度不連続を非常に薄い層を挟むことで表現している。
Let the \(z\)-coordinate range of the physical volume be \(-10000\leq z\leq 0\) [m], and let the P-wave velocity, S-wave velocity, and density be given as Table 2. This structure is realized by example3.ini below, where the velocity and density discontinuities are expressed by extremely thin layers.

表2(Table 2)
標高の範囲[m]
(Altitude range [m])
\(V_p\) [m/s] \(V_s\) [m/s] \(\rho\) [kg/m\(^3\)]
\(-1000\leq z\leq 0\) 2000 1000 2100
\(-2000\leq z\leq -1000\) 2500 1500 2200
\(-4000\leq z\leq -2000\) 3500 2000 2400
\(-10000\leq z\leq -4000\) 5000 3000 2600

(example3.ini)
0.1[TAB]2000.0[TAB]1000.0[TAB]2100.0
-1000.0[TAB]2000.0[TAB]1000.0[TAB]2100.0
-1000.1[TAB]2500.0[TAB]1500.0[TAB]2200.0
-2000.0[TAB]2500.0[TAB]1500.0[TAB]2200.0
-2000.1[TAB]3500.0[TAB]2000.0[TAB]2400.0
-4000.0[TAB]3500.0[TAB]2000.0[TAB]2400.0
-4000.1[TAB]5000.0[TAB]3000.0[TAB]2600.0
-10000.1[TAB]5000.0[TAB]3000.0[TAB]2600.0


(4)地表面に平行な境界面を含む構造 (A structure with layer boundaries parallel to the ground surface)

地表面の標高が水平位置\((x,y)\)の関数として\(h(x,y)\)で与えられる場合を考える。 物理領域の下端を\(z=-10000\) [m]とし、 地下構造が表3のように与えられるものとする。 この地下構造は以下のexample4.iniによって実現できる。
Let \(h(x,y)\) be the ground surface altitude that is a function of horizontal locations \((x,y)\). Let the lower limit of the physical volume be \(z=-10000\) [m], and let the subsurface structure be given by Table 3. This structure is realized by example4.ini below.

表3(Table 3)
標高の範囲[m]
(Altitude range [m])
\(V_p\) [m/s] \(V_s\) [m/s] \(\rho\) [kg/m\(^3\)]
\(h(x,y)-500\leq z\leq h(x,y)\) 2500 1500 2200
\(0\leq z\leq h(x,y)-500\) 3500 2000 2400
\(-10000\leq z\leq 0\) 5000 3000 2600

(example4.ini)
surface+0.1[TAB]2500.0[TAB]1500.0[TAB]2200.0
surface-500.0[TAB]2500.0[TAB]1500.0[TAB]2200.0
surface-500.1[TAB]3500.0[TAB]2000.0[TAB]2400.0
0.0[TAB]3500.0[TAB]2000.0[TAB]2400.0
-0.1[TAB]5000.0[TAB]3000.0[TAB]2600.0
-10000.1[TAB]5000.0[TAB]3000.0[TAB]2600.0