Οι σεισμικές ακολουθίες χαρακτηρίζονται συχνά από διακριτά μοτίβα, όπως χωρικές και χρονικές συστάδες, μονόπλευρη, αμφίπλευρη ή ακτινική μετανάστευση υποκέντρων και εξάρσεις στο ρυθμό σεισμικότητας. Αυτά εξαρτώνται από τη γεωμετρία του δικτύου ρηγμάτων, το καθεστώς τάσεων περιφερειακής κλίμακας με τοπικές ανωμαλίες εξαιτίας ανακατανομής τάσεων ύστερα από ισχυρά γεγονότα ή ασεισμικούς παράγοντες όπως ερπυσμό και διάχυση ρευστών υπό πίεση. Είναι σύνηθες στη χωρικά συσταδοποιημένη σεισμικότητα να εμφανίζονται σεισμοί με όμοιες εστιακές παραμέτρους (σμηνοσεισμοί, ή multiplets), οι οποίοι παράγουν όμοιες σεισμικές αναγραφές. Στην παρούσα διδακτορική διατριβή αξιοποιείται το φαινόμενο της ομοιότητας κυματομορφών με σκοπό τη βελτίωση της χωροχρονικής κατανομής των σεισμών ενός καταλόγου, τόσο σε ποσότητα δεδομένων, μέσω μεθόδων ημιαυτόματης ανίχνευσης και εντοπισμού, όσο και σε ποιότητα, εφαρμόζοντας επαναπροσδιορισμό υψηλής ευκρίνειας με τη μέθοδο των διπλών διαφορών που ελαχιστοποιεί τις αβεβαιότητες των σχετικών υποκεντρικών θέσεων. Η προτιμώμενη μέθοδος για τη μέτρηση της ομοιότητας κυματομορφών είναι το μέγιστο της συνάρτησης ετεροσυσχέτισης στο πεδίο του χρόνου. Τα χαρακτηριστικά της εξετάζονται ως προς διάφορες παραμέτρους, όπως τα συχνοτικά φίλτρα που εφαρμόζονται, το μήκος παραθύρου και η διάρκεια της σεισμικής πηγής, καθώς και τη βέλτιστη τιμή κατωφλίου για συσταδοποίηση με σύνδεση του «πλησιέστερου γείτονα». Αναπτύχθηκε μια ημιαυτόματη μέθοδος προσδιορισμού χρόνων άφιξης, βασισμένη στην τεχνική του «ανιχνευτή συσχέτισης» που αξιοποιεί την ομοιότητα κυματομορφών σε πολλαπλούς σταθμούς. Η μέθοδος αυτή έχει τροποποιηθεί σε έναν υβριδικό αλγόριθμο που συνδυάζει την αυξημένη ανιχνευσιμότητα ενός μεμονωμένου σταθμού με έναν υψηλής ευαισθησίας αυτόματο αλγόριθμο προσδιορισμού του χρόνου άφιξης σεισμικών φάσεων, ο οποίος βασίζεται στο Κριτήριο Πληροφορίας του Akaike και σε στατιστική ανώτερης τάξεως (κύρτωση και λοξότητα). Περιγράφεται, επίσης, μία μέθοδος προσδιορισμού του μεγέθους σεισμικής ροπής μέσα από φασματική ανάλυση, με έμφαση σε τροποποιήσεις που επιτρέπουν τον υπολογισμό του σε τοπικούς μικροσεισμούς με στενό συχνοτικό εύρος και χαμηλό λόγο σήματος προς θόρυβο. Οι αλγόριθμοι που αναπτύχθηκαν επιτρέπουν τον δεκαπλασιασμό της ποσότητας σεισμών ενός σεισμικού καταλόγου και τη βελτίωση της πληρότητάς του κατά περίπου μία τάξη μεγέθους.Παρουσιάζονται διάφορες μελέτες περιπτώσεων που περιλαμβάνουν πρότυπα μετασεισμικών ακολουθιών (όπως για τους σεισμούς της Ανδραβίδας το 2008, της Κεφαλονιάς το 2014 και της Λευκάδας το 2015) καθώς και σεισμικών σμηνών (όπως στη Λίμνη Τριχωνίδα το 2007, την Οιχαλία το 2011 και την Ελίκη το 2013). Η κύρια περιοχή μελέτης είναι ο δυτικός Κορινθιακός Κόλπος, με ειδική έμφαση σε σεισμούς που καταγράφηκαν την περίοδο 2000-2007 από το τοπικό δίκτυο Corinth Rift Laboratory (CRL). Οι επαναπροσδιορισμένοι κατάλογοι, ύστερα από χωρική συσταδοποίηση και διαχωρισμό σε διαδοχικές χρονικές περιόδους, διερευνώνται μέσω πολλαπλών εγκάρσιων τομών και χωροχρονικών προβολών. Η εξέλιξη του ιστορικού των σμηνοσεισμών δίνει μία εικόνα για τον τρόπο με τον οποίο εξαπλώνεται η σεισμικότητα σε νέα τμήματα ρηγμάτων που δεν είχαν διαρρηχθεί, την προϋπάρχουσα δραστηριότητα σε περιοχές που αργότερα έδωσαν ισχυρές σεισμικές συστάδες και την επαναλαμβανόμενη σεισμικότητα σε τμήματα ρηξιγενών ζωνών. Όπου τα δεδομένα το επιτρέπουν, κατασκευάζονται σύνθετοι μηχανισμοί γένεσης για χωρικές συστάδες ή ομάδες σμηνοσεισμών που βοηθούν συμπληρωματικά στην ερμηνεία της χωρικής ανάλυσης και τον προσδιορισμό των ενεργοποιημένων ρηξιγενών δομών. Αναζητούνται μοτίβα χωροχρονικής σεισμικής μετανάστευσης που αποκαλύπτουν την πιθανή επίδραση ρευστών κατά τη διάρκεια σεισμικών σμηνών και προσδιορίζονται παράμετροι όπως η υδραυλική διαχυτικότητα και η ταχύτητα μετανάστευσης. Τέλος, εξετάζεται η συνεισφορά ασεισμικών παραγόντων μέσα από μοντέλα ETAS (Μετασεισμική Ακολουθία Σποραδικού Τύπου) που παρέχουν, επιπροσθέτως, πληροφορία για την παραγωγικότητα μετασεισμών και τον ρυθμό απόσβεσης.
3D Seismic Tomography to Image the Subsurface Structure of “IY” Geothermal Field Using Double-Difference Method and Waveform Cross-Correlation: Preliminary Results
