44
1 Non-DC složky Non-DC složky vícenásobných vícenásobných zemětřesení. zemětřesení. J. Zahradník J. Zahradník
1
Non-DC složky vícenásobných Non-DC složky vícenásobných zemětřesení.zemětřesení.
J. ZahradníkJ. Zahradník
2
3
4 x M5, 20070410
4
Rozklad momentového tenzoru
M=M_iso + M_deviaM_devia=M_dc + M_clvd
Vl. čísla M_devia: 1, -f, f-10<=f<=0.5(1, -f, f-1) = (1-2f) . (1,0,-1) + 2f . (1,-1/2,-1/2)
(1-2f).100% … DC%
(pozor SED !)
5
Skalární seismický moment ?
(1, -f, f-1) = (1-2f) . (1,0,-1) + 2f . (1,-1/2,-1/2)
Moment jako max. vl. čísloMoment jako |Eukl. norma / sqrt(2)Moment jako aritm. průměr 1 a |f-1| ‘best DC (Harvard)’
6
Zakynthos, April 2006, 3x M~5.5
3Events DATE ORIGIN Event20060411 0002 41.25 #120060411 1729 28.64 #220060412 1652 01.43 #3
7
agenturní řešení
8
93 24 97
61
74 30 96
DC%
9
83% 21% 58%
naše řešení
10
redukce variance ~ 0.8
11
12
13
14
15
Závěr 1. části: DC% je různé, nejmenší pro jev č. 2
Následuje 2. část:interpretace non-DC
(porovnáním devia a DC-constr.)
16
17
18
19Sub1_DEV=Sub1_DC + Sub2_DC
20Sub1_DC =DC část celkového zdroje
21
Automatický výsledek iterat. dekon. je blízký jinému formálnímu
rozkladu momentového tenzoru:
M_devia=M_dc + M_clvdM_devia=M_dc1 + M_dc2
(1, -f, f-1) = 1. (1,0,-1) + f . (0,-1,1) P,T,N N,P,T
22
Závěr 2. části:
ekvivalence mezi modelem jednoduchého non-DC zdroje a
modelem dvojice DC zdrojů
23
Část 3:
Lze věřit jen součtu dvou DC zdrojů nebo každému zvlášť ?
(syntetické testy bez šumu)
24
25
10-sec apart
two identical
single
~1se
c ap
art
26
Závěr 3. části:
obecně nelze zaručit, že každý ze dvou DC zdrojů je správný
27
Část 4
Jak moc vadí malý počet stanic ?Proč se DC% silně mění s časem?Proč se tolik liší agenturní DC%?
(syntetický test)
28
MOMENT STRIKE DIP RAKE .183458E+18 211. 80. 122. .713815E+17 290. 76. -88. mutual time shift: 1 secfreq range: 0.02 0.025 0.08 0.10
tensor sum (no time shift) Mo: .182E+18DC: 32%strike, dip, rake 210.1 86.01 125.6
29
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
2,3,4 6,7,8
all
30
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
2,3,4 6,7,8
all0 0.2 0.4 0.6 0.8 1
CO RRELATIO N
0
20
40
60
80
100
DC
%
all VLF
31
Závěr 4. části
rychlá změna DC% s časem není důsledkem malého počtu stanic
změna DC% se zpomalí při nižší frekvenci, pak se ale DC% změní
32
Zásadně jiná je situace, kdy „studované frekvence vidí“ i prostorový posun zdrojů.
DC% násobných jevů je pak veličinou
azimutálně závislou,
podobně jako zdánlivá časová funkce.
330 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
2,3,4 6,7,8
6,7,8 VLF
0 0.2 0.4 0.6 0.8 1
CO RRELATIO N
0
20
40
60
80
100
DC
%
0 0.2 0.4 0.6 0.8 1
CORRELATIO N
0
20
40
60
80
100
DC
%
3480 90 100 110 120
tim e (se c)
0
2E+014
4E+014
6E+014
8E+014
1E+015
mo
me
nt (
Nm
)
4925
2194
77
54
37
12
51
55
36
18
53
74
36
16
21
9384
5026
0
20
40
60
80
100
DC
%DC% =
3580 90 100 110 120tim e (se c)
0
2E+014
4E+014
6E+014
8E+014
1E+015
mo
me
nt (
Nm
)
4925
2194
77
54
37
12
51
55
36
18
53
74
36
16
21
9384
5026
0
100
200
300
400
stri
ke (
de
g)
DC% =
36
Amfilochia20021231
Mw 4.9(P. Adamová)
40%
30%
24%
DC%
DC=24% získánoaž inverzí synt.datodpovídajících dvojitému zdrojipro f < 0.025 Hz !
37
Stejný frekv. obor jako pro Zakynthos,
ale menší M.Přesto však
„menší překryv“ ??Prodleva 4 sec !!!
Neplatí sub1_devia=
sub1_DC+sub2_DC
65%61%
38
Dvojitý model zdroje vysvětlil agenturní DC%, protože umožnil
syntézu nízkých frekvencí.
39
40
Závěry• Na f<1/T (překryv) se vícenásobný zdroj jeví jako non-DC, DC%
vykazuje variace desítek procent i pro optimalizovanou polohu zdroje a stabilní orientaci MT.
• DC% závisí na frekv. oboru; stabilní DC% pro nízké frekvence f<<1/T.
• K odhadu DC% pro f<<1/T lze použít synt. model násobného zdroje odvozený na f<1/T.
• Agenturní nízké hodnoty DC% pro záp. Řecko velmi běžné a použitelné jen s velkou opatrností.
41
42
43
44
