1.2 震源モデル(2)yagi-y/text/2017geodynamics04.pdf · p軸は σ 1方向 t軸はσ 3方向...
TRANSCRIPT
放射パターン(震源メカニズム)の記述 3Dの震源メカニズムを2Dで表現する。
断層面の記述 断層面の記述方法を理解する。
断層パラメターと巨視的な断層モデル ハスケルモデルや円形クラックモデルの概要を理解する。
断層運動による波動 波動が観測されるのか単純化して理解する。
第4回 地球変動科学
1.2 震源モデル(2)
参考文献: ① Aki & Richards, Quantitative Seismology ② 松浦,地球連続体力学(5章),岩波講座地球惑星科学
1.2 震源モデル
放射パターンの表示方法(震源メカニズム解)その1震源球に押し引き分布を書き込んでも3D情報なので、紙に書きにくい。 → 等積投影法を用いて下の震源球を2Dで表記する。
1.2 震源モデル
放射パターンの表示方法(震源メカニズム解)その2典型的な断層運動と震源メカニズム解 震源メカニズムと応力場
の関係 第0近似で、
P軸は σ1方向 T軸はσ3方向 B軸はσ2方向 に近いはず。
クーロンの破壊規準を理解していれば、この近似はバイアスがかかる近似である事が分かる。(調べてみよう)
北
正断層
逆断層
横ずれ断層
T軸 P軸
T軸P軸
T軸P軸
東落ちの正断層西落ちの正断層
東落ちの逆断層西落ちの逆断層
北
北
北
北
左横ずれ断層
北
右横ずれ断層
正断層
逆断層
横ずれ断層
T軸
P軸
B軸2
3 3
2
T軸P軸
B軸
1
3
1
3
T軸
P軸B軸
2
1 1
2
1.2 震源モデル
断層面の記述 その1
断層の幾何学
断層面は、走向と傾斜 すべり方向は、すべり角 三つのパラメターで断層運動を記述できる。
走向:北から断層面までの角度で時計回りに計測(0~360度) 傾斜角:水平面と断層面がなす角で時計回りに計測(0~90度) すべり角:走向ベクトルとすべりベクトルのなす角度で、反時計回りに計測(-180~180度)
走向ベクトルの方向はどっち?
答え:Dipの角度が0~90度に収まる方
n
Strike direc
tion
Slip directionRake
Strike
Dip
North
Depth
1.2 震源モデル
断層面の記述 その2
法線ベクトルと走向( )・傾斜( )
地表面に断層を投影し、投影された法線ベクトルと北方向を指すベクトルとの角度に着目すると、
φs δ
s
n1
n2
n
North
East
tan 180 − φs + 90( )⎡⎣ ⎤⎦ = − n2n1
∴φs = arctan − n1n2
⎛⎝⎜
⎞⎠⎟
断層を真横からみると、
sin 90 −δ( ) = sin δ − 90( ) = cosδ = − n3n
∴δ = arccos −n3( )
n
Depth
cross section
90n3
法線ベクトルから走向と傾斜が求まった。
1.2 震源モデル
断層面の記述 その3
すべり角度( )は?すべりベクトル を、断層の断面図に投影すると、
λ
− ν3sinδ
走向・傾斜を軸にして見てみると、
sin 180 − λ( ) = −sin λ −180( ) =−v3 sinδd
∴λ = arcsin − v3sinδ
⎛⎝⎜
⎞⎠⎟
すべり角度が求まる。
となる。
ν
Depth
cross section
v3sinv3
Strike
Dip
v3sin
1.2 震源モデル
断層面の記述 その4逆に、断層面パラメターから法線ベクトルは、
n =−sinδ sinφssinδ cosφs−cosδ
⎛
⎝
⎜⎜⎜
⎞
⎠
⎟⎟⎟
スリップベクトルは、
ν =cosλ cosφs + cosδ sinλ sinφscosλ sinφs − cosδ sinλ cosφs
−sinδ sinλ
⎛
⎝
⎜⎜⎜
⎞
⎠
⎟⎟⎟
となる。
レポート①:P軸ベクトル、T軸ベクトル、B軸ベクトルを求めよ。
n
Strike direc
tion
Slip directionRake
Strike
Dip
North
Depth
1.2 震源モデル
断層パラメター(断層運動に関わるもの) その1
断層運動により地震波が発生する。 断層運動のプロセスを総称して「震源過程」と呼ぶ。
M6クラスでは、3~5秒かけて断層を破壊、断層差し渡りの長さ10km程度
M8クラスでは、30秒~1分かけて断層を破壊、断層差し渡りの長さ100km程度
このようなものが、断層運動を表記する、断層パラメター
静的震源パラメター
最終的な断層の長さ(L)幅(W)食い違い量(D)
1.2 震源モデル
巨視的な断層モデル その1ハスケルモデル巨大地震に主に適用される巨視的断層モデル
長さ(L)幅(W)の断層の端から走向方向に破壊が破壊伝播速度(Vr)で進行。 断層上の各点のすべり継続時間(τ)は一定 つまり、ある瞬間にすべっている幅は一定。
円形クラックモデル中小地震に主に適用される巨視的断層モデル
半径(a)の円形の断層モデル 破壊は中心から開始し、外側に向かって破壊伝播速度(Vr)で進行。
Tim
e
Origintime
Fault area
Slipping area
Rupture direction
Vr
Vr
Vr
Tim
e
Origintime
Fault area
Slipping area
Rupture directionVr
Vr
Vr
1.2 震源モデル
巨視的な断層モデル その2
大地震の断層の幅は地震が発生しうる領域によって上限が定まる。 → ハスケルモデル
中地震は、上限値とは無関係に破壊が広がる。 → 円形クラックモデル
とは言えフラクタル的な性質により、 「大きな地震でも、小さな地震でもすべり分布は複雑」 つまり、巨視的な断層モデルには限界がある。
地表付近の横ずれ断層の場合は、地震発生領域の幅が狭いため、第0近似として
Plastsphere
SchizosphereHypocenter
Centroid
Large Earthquake
Surface
1.2 震源モデル断層運動による波動 その1震源時間関数とモーメント速度関数
断層面の各点から震源までの距離は変化する 面倒だから、遠方近似をする。
震源時間関数Scは、断層面のすべり量分布を距離変化による時間遅れ効果を加えて面積分したものになる。遠方近似を用いると、
Sc t,γ( ) = µ !D ξ ,t + ξ iγ c( ) dS∫∫
仮想的に断層面直上では なので ξiγ = 0
!M 0 t( ) = µ !D ξ,t( )dS∫∫これをモーメント速度関数と呼ぶ。
だけ波の到達時間が変化する。
断層面上のすべり速度関数
0
・
t = ξ iγ c
遠方近似
十分にはなれている観測点では、波の波源が変化した場合
1.2 震源モデル断層運動による波動 その2シンプルな場合を考える 断層すべりが場所に依存しない場合
すべり速度の関数は、
と書けるので、モーメント速度関数は、
!M 0 t( ) = µ !D t − s( ) !S s( )ds∫ = µ !D∗ !S
点ξに破壊が到達する時間 !D ξ ,t( ) = !D t − s ξ( ) ( )
モーメント速度関数はすべり速度関数と面積の拡大速度のコンボリューションとなる。
地震モーメントと震源時間関数
Sc t,γ( )dt0
∞
∫ = dt0
∞
∫ µ !D ξ,t − s( )dS∫∫= µD ξ,∞( )dS∫∫ = M 0
よって、観測点の位置にかかわらず、Scを時間積分したものはM0となる。
やや深発地震のモーメント速度関数
補足説明コンボリューション(畳み込み積分)その1
線形現象の場合、足し合わせる事によって、その影響を考慮する事ができる。
A t( ) = B t( )*C t( )
= B τ( )C t −τ( )dτ−∞
∞
∫デルタ関数の場合は、
A t( ) = δ τ( )C t −τ( )dτ−∞
∞
∫ = C t( )
A t( ) = δ τ − t1( )+δ τ − t2( )⎡⎣ ⎤⎦C t −τ( )dτ−∞
∞
∫ = C t − t1( )+C t − t2( )
デルタ関数が二つの場合は、
* =
t t t
t
convolution
t t
t
convolution
t
0 0 0 0
0 0 0 0 0
* =
t = 2 t = 4t = 3t =0 t = 1
t = 4 t = 6t = 0 t = 2
補足説明コンボリューション(畳み込み積分)その2
A t( ) = B t( )*C t( )
= B τ( )C t −τ( )dτ−∞
∞
∫
フーリエ変換すると
A ω( ) = B ω( )× C ω( )
単純なかけ算となる。
レポート②:コンボリューションとフーリエ変換を調べて、それぞれ1枚(A4)のレポートにまとめよ。 また、長さが違うボックスカー関数をコンボリューションすると台形になることを確認せよ。
コンボリューションの式は、 フーリエ変換(復習)
A ω( ) = A t( )e− iωt dt−∞
∞
∫逆変換
時間or空間変数の関数を周波数変数の関数に変換で、可逆変換。
A t( ) = 12π A ω( )eiωt dω−∞
∞
∫
相関関数を思い出すと、フーリエ変換は、ある周波数をもった関数との相関値を表す関数を作り出すことだと分かる(物理のがきしっぽ)