Non-DC složky vícenásobných zemětřesení.

Post on 18-Jan-2016

36 views 0 download

description

Non-DC složky vícenásobných zemětřesení. J. Zahradník. 4 x M5, 20070410. Rozklad momentov é ho tenzoru. M=M_iso + M_devia M_devia=M_dc + M_clvd Vl. čísla M _ devia : 1, -f, f-1 0

transcript

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