Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Liczenie MT oparte jest na programie hybridmt-v1.2.2 Grzegorz Kwiatka (Wiejacz 1992, Kwiatek i Martinez-Garzon 2016). Poniższa tabela pokazuje odpowiedniki składowych tensora momentu sejsmicznego:

...

Mechanizm jest liczony z amplitud pierwszych wejść fal P lub pierwszych wejść fal P i S. Zakłada się zbliżony do prostokątnego proces rozrywu  co pozwala na pomińcie określania pełnej funkcji Greena i przebiegu czasowego procesu rozrywu. Według Aki i Richardsa (2002), Fitch i in. (1980) oraz De Natale i in. (1987), zarejestrowane przemieszczenie dla fali P wynosi

...

gdzie ρ to średnia gęstość medium, r to odległość źródło-odbiornik, αβ to średnie prędkości fal P i S u źródła, TP i TS to czasy podróży fal P i S , M jest tensorem momentu sejsmicznego, l jest kierunkiem fali w odbiorniku, γ jest kierunkiem fali w źródle, a ̂Pγ i ̂PI są kierunkami poprzecznymi w źródle i odbiorniku. MT uzyskuje się przez rozwiązanie układu N powyższych równań, gdzie N jest liczbą faz sejsmicznych na stacjach rejestrujących zdarzenie. Tensor momentu otrzymujemy rozwiązując zbiór N równań na na U, gdzie N jest liczbą stacji, które zarejestrowały zdarzenie. Sześć niezależnych składowych tensora momentu wymaga minimum sześć równań, ale im więcej, tym lepiej. System równań jest traktowany wtedy jako nadmiarowy i rozwiązany za pomocą metody najmniejszych kwadratów (norma L2), gdzie funkcja z kosztu jest  sumą kwadratów reszt. Niepewność oszacowania MT można jest liczona za pomocą znormalizowanego błędu średniej kwadratowej (RMS) między teoretyczną a szacowaną amplitudą (Stierle i in., 2014a, 2014b)

W pierwszej kolejności liczony jest pełny tensor. Następnie warunek śladu zerowego (trace null - TN) jest nakładany na rozwiązanie i minimalizuje się funkcję kosztu powtórnie. Daje to nam rozwiązanie mechanizmu wstrząsu bez zmiany objętości w źródle. Kiedy warunki zerowy ślad i wyznacznika zerowy są ustawione, rozwiązanie jest ograniczone do podwójnej pary (double couple -DC), czyli ścinania prostego.. Cały moment roztwór można również rozłożyć na składowe: izotropową (ISO), skompensowany liniowy dipol wektorowy (CLVD) i części mechanizmu ścinający (DC).

...

Sposób stary oparty jest bezpośrednio na pracy Wiejacza *(1992) i programie hybridmt-v1.2.2. Liczy on MT tylko ze składowych pionowych wyznaczając amplitudę fali P według wzoru

...