Στατιστική των σεισμών Κεφ.13 Θ.Σώκος Εργαστήριο Σεισμολογίας Τμήμα Γεωλογίας Η στατιστική ή των σεισμών ώ ασχολείται λ ί με τη μελέτη της κατανομής των σεισμών λαμβάνοντας υπ’ όψη σαν κύρια παράμετρο το σεισμικό μ μ μέγεθος γ ς Η ανάλυση γίνεται με καθαρά πιθανολογικά κριτήρια χωρίς να λαμβάνεται υπ’ όψη το σεισμοτεκτονικό ό καθεστώς θ ώ Βασικές έννοιες από τη θεωρία των πιθανοτήτων Η έννοια της Πιθανότητας – Έστω ότι από ένα σύνολο n παρατηρήσεων μιας παραμέτρου Υ βρίσκουμε την τιμή Υ=Υi , ni φορές, η ποσότητα ni/n ονομάζεται σχετική συχνότητα εμφάνισης και εκφράζει την στατιστική πιθανότητα P(Y) της παραμέτρου Υ. ni P(Y ) = n Βασικές έννοιες από τη θεωρία των πιθανοτήτων Παράδειγμα Σε μια περιοχή έχουν γίνει n σεισμοί από τους οποίους ί οι k είναι ί μεγαλύτεροι λύ από ό6 6. P((> > 6) = k / n n−k k F= = 1− = P ⇒ F + P = 1 n n Βασικές έννοιες από τη θεωρία των πιθανοτήτων Συνάρτηση κατανομής πιθανότητας Έστω Μ τυχαίο σεισμικό μέγεθος και m ένας πραγματικός αριθμός που μπορεί να λάβει τιμές m1 ≤ m ≤ mn. Έστω Ρ η πιθανότητα το τυχαίο μέγεθος Μ να είναι μικρότερο του m FM(m)=P(M<m) Στην περίπτωση που ο αριθμός m παίρνει όλες τις δυνατές τιμές στο διάστημα (m1,mn) τότε η πιθανότητα FM(m) περιγράφει μια συνάρτηση που ονομάζεται συνάρτηση κατανομής πιθανότητας της τυχαίας μεταβλητής Μ. Βασικές έννοιες από τη θεωρία των πιθανοτήτων FM((m)) Fm(mn)=1 1 Fm(m1)=0 0 m1 Γενική μορφή της συνάρτησης κατανομής πιθανότητας mn m Η κατανομή Poisson Έστω λ έναςς θετικόςς αριθμός ρ μ ς και έστω ότι η τυχαία μεταβλητή Ν παίρνει τιμές μ ς 0,1,2,3… , , , Αν η πιθανότητα η P(Ν=k (Ν=k) παρέχεται από τη σχέση −λ e λ P( N = k ) = k! k Τότε η αντίστοιχη συνάρτηση κατανομής πιθανότητας ονομάζεται κατανομή Poisson Η πιθανότητα να συμβούν k «γεγονότα» αν ο ρυθμός εμφάνισης είναι λ. Στατιστική των σεισμών Σ Στατιστικά ά το απλούστερο λ ύ μοντέλο έλ για την περιγραφή ή της σεισμικότητας με το χρόνο είναι αυτό της κατανομής Poisson. Η κατανομή αυτή υποθέτει ότι οι σεισμοί είναι μεμονωμένα γεγονότα και η γένεση έ του ενός ό δεν δ επηρεάζει άζ την γένεση έ των υπολοίπων. λ ί Εάν α είναι ο μέσος αριθμός εμφάνισης σεισμών σε ένα χρονικό διάστημα t , η πιθανότητα n σεισμοί να γίνουν μέσα σ’ αυτό το χρονικό διάστημα είναι (at ) n e − ( at ) P( N = n, t ) = n! Wanner, 1937 Στατιστική των σεισμών Π άδ Παράδειγμα … Αν σε μία περιοχή έχουμε 5 σεισμούς κάθε αιώνα με μέγεθος > 6 τότε η πιθανότητα να έχουμε ένα σεισμό με μέγεθος Μ>6 σε 10 χρόνια θα είναι (α=5/100 α=5/100): ): ((5 / 100) *10)1 e − ( 5 /100 )*10 = 0.3 P(n = 1, t = 10) = 1! Η πιθανότητα η να μην μη συμβεί μβ κανέναςς σεισμός μ ς Μ>6 P(N=0)=e-(5/100)*10=0.61 Κατανομή της συχνότητας μεγεθών Σχέση Gutenberg – Richter Ένας από τους σημαντικότερους στατιστικούς νόμους στη Σεισμολογία χρησιμοποιείται για να περιγράψει την κατανομή των σεισμικών μεγεθών σε μια ή Προτάθηκε Π άθ ό τους G t b περιοχή. από Gutenberg – Richter το 1941 μετά από μελέτη της σεισμικότητας στην Καλιφόρνια. Κατανομή της συχνότητας μεγεθών Σχέση Gutenberg – Richter Οι Gutenberg – Richter βρήκαν ότι ο αριθμός ρ μ ς των σεισμών μ n((Μ)) μεγέθους μ γ ς Μ ΔΜ σε μια περιοχή και για χρονικό διάστημα k ετών, ετών εκφράζεται από την εμπειρική σχέση Logn(M)=a Logn(M) Logn(M)=aa-bM Κατανομή της συχνότητας μεγεθών Σχέση Gutenberg – Richter O Utsu το 1961 έδειξε ότι είναι προτιμότερο ρ μ ρ να χρη χρησιμοποιείται μ ο αριθμός ρ μ ς των σεισμών με μέγεθος ίσο ή μεγαλύτερο από Μ (δηλ (δηλ. η συσσωρευτική συχνότητα) logNk(M) (M)=a ak-bM bM Ν(Μ)=Σn(M)dm συσσωρευτική ρ ή συχνότητα χ η Όσο αυξάνει ξ το μ μέγεθος γ ς μειώνεται ο αριθμός των σεισμών σε μία περιοχή Κατανομή της συχνότητας μεγεθών Σχέση Gutenberg – Richter Gutenberg-Richter: logN=a-bM •N αριθμός σεισμών μεγαλύτερων από Μ lo ogN •a, b σταθερές 10a b 0 Magnitude a, b υπολογίζονται από τους σεισμικούς καταλόγους Κατανομή της συχνότητας μεγεθών Σχέση Gutenberg – Richter Ο παράμετροι Οι ά a, b είναι ί πολύ λύ σημαντικοίί στη Σεισμολογία, ιδιαίτερα η b, – η παράμετρος a εκφράζει τον λογάριθμο του αριθμού των σεισμών με μέγεθος 0 και μεγαλύτερο λύ και η τιμής ή της μεταβάλλεται βάλλ από ό περιοχή σε περιοχή Εξ ά από ό: – Εξαρτάται – Τη χρονική περίοδο των παρατηρήσεων – Την έκταση της περιοχής έρευνας – Την η σεισμικότητα μ η Κατανομή της συχνότητας μεγεθών Η παράμετρος b Η παράμετρος ά b εκφράζει άζ τον ρυθμό θ ό αύξησης του αριθμού των σεισμών καθώς μειώνεται το μέγεθός τους. Οι τιμές μ ς της ης παραμέτρου ρ μ ρ b κυμαίνονται μ από 0.4 έως 1.4 Η παράμετρος b εξαρτάται από τις τάσεις και από τις μηχανικές ιδιότητες του υλικού της εστιακής περιοχής περιοχής, άρα και από την τεκτονική. Κατανομή της συχνότητας μεγεθών Η παράμετρος b Επίσης το b μεταβάλλεται με το βάθος (καθώς το βάθος αυξάνει, το b μειώνεται) καθώςς και με μ τη η γεωλογική γ γ ή ηλικία η της ης περιοχής, υψηλές τιμές έχουν βρεθεί για τις Αλπικές ζώνες (1 (1.0 (1.00-1.8) 01 8) και ενδιάμεσες για τις ηπειρωτικές πλατφόρμες (0.6 (0.6--0.7) Κατανομή του b στην Ελλάδα Κατανομή της συχνότητας μεγεθών Ποσοτικά μέτρα Από την σχέση GutenbergGutenberg-Richter προκύπτουν διάφορα ποσοτικά μέτρα για τον στατιστικό καθορισμό της σεισμικότητας π.χ. Ο ετήσιοςς αριθμόςς σεισμών με μέγεθοςς Μ ή μεγαλύτερο a 10 N M = bM 10 Κατανομή της συχνότητας μεγεθών Ποσοτικά μέτρα Η μέση περίοδος επανάληψης των σεισμών ΤΜ bM 10 TM = = 1/ NM a 10 Κατανομή της συχνότητας μεγεθών Ποσοτικά μέτρα Τ πιθανότερο Το θ ό μέγιστο έ μέγεθος έ θ Μ*, Μ* σε k χρόνια a + log k M* = b Κατανομή της συχνότητας μεγεθών Ποσοτικά μέτρα To συχνότερα παρατηρούμενο μέγιστο ετήσιο μέγεθος Μ1=a/b H πιθανότητα P να συμβεί σεισμός μεγέθους έθ Μ ή μεγαλύτερου λύ σε χρόνο ό t P = 1 − exp(−10 a−bM t) Σχέση Gutenberg – Richter για την Ελλάδα Μέση περίοδος επανάληψης Μέση περίοδος επανάληψης Υπολογισμός των παραμέτρων a a,b b Μέθοδος της μέσης τιμής – Σεισμικός μ ς κατάλογος γ ς ((πληρότητα, ηρ η , κλίμακα μ μεγέθους, εξαρτημένα γεγονότα) – Η συχνότητα των σεισμών κάθε μεγέθους ανάγεται σε χρονικό διάστημα ενός έτους και στην συνέχεια υπολογίζεται η αθροιστική συχνότητα Μέθοδος των ακραίων τιμών Χρησιμοποιούμε μόνο τις τιμές των μέγιστων μεγεθών των σεισμών Υπολογισμός των παραμέτρων a,b ab Καταγράφονται οι σεισμοί ανάλογα με το μέγεθός τους σε διάφορες τάξεις μεγεθών Μ1,Μ Μ2…Μ Μn με ορισμένο βήμα ταξινόμησης ΔΜ, έτσι για κάθε τιμή μεγέθους Μi αντιστοιχεί ορισμένος αριθμός σεισμών n1,n2…nn . Στη συνέχεια υπολογίζεται η συσσωρευτική κατανομή των μεγεθών και δημιουργείται το γράφημα LogN – M. Οι τιμές των ρ μ ρ a,, b παραμέτρων υπολογίζονται με την μέθοδο των ελαχίστων τετραγώνων. Ο σεισμικός κύκλος Σεισμικές ακολουθίες Σεισμική Σ ή ακολουθία λ θί ή σεισμική ή σειρά ά ονομάζεται το σύνολο των σεισμικών δονήσεων που συμβαίνουν σε ένα μικρό χρονικό διάστημα Ο μεγαλύτερος σεισμός της ακολουθίας καλείται σεισμοί α ε α κύριος ύρ ος σεισμός σε σμός εενώ ώ οι ο σε σμο που προηγούνται και ακολουθούν τον κύριο σεισμό ονομάζονται, προσεισμοί και μετασεισμοί αντίστοιχα. Σεισμικές ακολουθίες - Προσεισμοί Μ<7 7<Μ<7.8 Δεν παρατηρείται ελάττωση του χρόνου αύξησης της προσεισμικής δράσης με το μέγεθος έ θ του κύριου ύ σεισμού. Ρυθμός εμφάνισης προσεισμών με τον χρόνο n=at-c Μ>7.8 Όπου α,c σταθερές με την c να προσεγγίζει το 1 Σεισμικές ακολουθίες - Προσεισμοί Από τη μελέτη των προσεισμών σεισμών με μέγεθος > 6 προκύπτει η παρακάτω σχέση για την Ελλάδα g logP=2.13-0.57Μ for όπου Mfor είναι το μέγεθος του προσεισμού προσεισμού. Από την παραπάνω σχέση προκύπτει ότι η πιθανότητα να συμβεί προσεισμός μεγέθους >4.4 είναι 42% , για κύριους σεισμούς > 6.0 Σεισμικές ακολουθίες - Μετασεισμοί Η συχνότητα εμφάνισης μετασεισμών είναι αντιστρόφως ρ φ ς ανάλογη γη του εστιακού βάθους και της ηλικίας των πετρωμάτων στην εστιακή περιοχή Σεισμικές ακολουθίες Η σχέση Gutenberg – Richter ισχύει και για την περίπτωση των σεισμικών ακολουθιών. ακολουθιών Έχει δειχθεί ότι η παράμετρος b είναι μικρότερη για τους προσεισμούς. προσεισμούς Μέσες τιμές είναι 0.67 για προσεισμούς και 0.92 για μετασεισμούς Σεισμικές ακολουθίες – Σμηνοσειρές (Swarm) Στην περίπτωση που όλοι οι σεισμοί μιας μ ής ακολουθίαςς συμβαίνουν μβ σε μια μ σεισμικής πολύ περιορισμένη περιοχή και κανένας σεισμός δεν έχει μέγεθος αρκετά μεγαλύτερο από του υπόλοιπους η ακολουθία αυτή ονομάζεται σμηνοσειρά (swarm) και οι σεισμοί , σμηνοσεισμοί σμηνοσεισμοί.. Σεισμικές ακολουθίες – Σμηνοσειρές (Swarm) Οι σμηνοσεισμοί συμβαίνουν συνήθως σε ηφ ηφαιστειακέςς περιοχές ρ χ ς και έχουν χ μ μικρό ρ μέγεθος κατά κανόνα Πολλές φορές παρουσιάζουν περιοδικότητα Η παράμετρος b παρουσιάζει μεγάλη απόκλιση από την τιμή 1 (συνήθως 2,2.5) για τις σμηνοσειρές Παράδειγμα ρ γμ σμηνοσειράς μη ρ ς Είδη σεισμικών ακολουθιών Σχέση μεταξύ μεγέθους κύριου σεισμού και μετασεισμού Συνήθως η διαφορά των μεγεθών είναι της ξης του 1.2 και αυξάνεται ξ γγια σεισμούς μ ς τάξης με βάθος μεγαλύτερο των 50km. 50km. Ο χρόνος που μεσολαβεί ανάμεσα στον κύριο σεισμό και τον ισχυρότερο ό του δεν είναι ί θ ό και μετασεισμό σταθερός ξ ρ από το μ μέγεθος γ ς το β βάθοςς του εξαρτάται σεισμού καθώς και από τα γεωτεκτονικά χαρακτηριστικά της περιοχής. Σχέση μεταξύ μεγέθους κύριου σεισμού και μετασεισμού Η πιθανότητα να συμβεί ο μεγαλύτερος μετασεισμός σε χρόνο μεγαλύτερο από Τ1 μετά τη γένεση του κυρίου σεισμού δίνεται από: ( 1))=0.53 )=0.53--0.2Τ1 Ρ(Τ Από την παραπάνω σχέση η πιθανότητα να συμβεί ο μεγαλύτερος μετασεισμός μετά την πρώτη μέρα είναι 53% και μετά την πρώτη εβδομάδα 34%. Χρονική κατανομή της συχνότητας των μετασεισμών Η συχνότητα του αριθμού των μετασεισμών μετά από κάθε κύριο σεισμό ελαττώνεται με την πάροδο του χρόνου χρόνου. Η ελάττωση αυτή ακολουθεί ένα νόμο που είναι γνωστός ως ο νόμος του Omori και περιγράφεται από την σχέση: c n= p (k + t ) Όπου n η συχνότητα των μετασεισμών k, k c, c και p είναι μετασεισμών, σταθερές που εξαρτώνται από το μέγεθος του σεισμού , το p παίρνει τιμές από 1.0-1.4. 1014 Χωρική κατανομή των μετασεισμών και τρόπος έκλυσης της μετασεισμικής ενέργειας Η χωρική ή κατανομή ή των μετασεισμών, ώ η συχνότητα και το μέγεθός τους καθορίζονται από την κατανομή των κλείθρων κλείθρων, τις φυσικές ιδιότητες της ρηξιγενούς επιφάνειας και από το μέγεθος της αποθηκευμένης δυναμικής ενέργειας ενέργειας. Εμπειρική σχέση ανάμεσα σε «μετασεισμικό όγκο,,εμβαδόν» και μέγεθος κυρίου σεισμού όγκο logV=9.58+1.47M logV=9 58+1 47M logA=1.02M+6.0 logA=1 02M+6 0 Έχει σημασία για τον καθορισμό της ρηξιγενούς επιφάνειας να χρησιμοποιούνται ύ τα δ δεδομένα δ έ των πρώτων ώ 1 1-2 2 ημερών ώ Χωρική και χρονική κατανομή των μετασεισμών του σεισμού μ της ης Sumatra ((26/12/2004 - 28/03/2005))
© Copyright 2024 Paperzz
