Η παρούσα διπλωματική εργασία ασχολείται με την ανάλυση της τυρβώδους ροής σε ορθογωνικό, ανοικτό αγωγό τριών διαστάσεων, στον οποίον ενυπάρχουν εμπόδια της μορφής προβόλων (groynes). Τέτοιου είδους κατασκευές χρησιμοποιούνται σε πάρα πολλούς ποταμούς ώστε να διατηρείται η επιθυμητή διατομή και το επιθυμητό βάθος ροής αλλά και να αποφεύγεται η διάβρωση των οχθών σε έντονα πλημμυρικά φαινόμενα. Για την επίλυση του προβλήματος αξιοποιήθηκαν οι εξισώσεις RANS, ενώ για το κλείσιμο της τύρβης χρησιμοποιήθηκε το μοντέλο δύο εξισώσεων και το μοντέλο . Η διαχείριση της ελεύθερης επιφάνειας έγινε με την μέθοδο Volume of Fluid (VOF), ενώ η αριθμητική επίλυση βασίστηκε στην μέθοδο των πεπερασμένων όγκων και πραγματοποιήθηκε με τον εμπορικό κώδικα FLUENT 6.1.2. Για την ροή στον υπό εξέταση αγωγό θεωρήθηκε αριθμός Reynolds , ύψος τραχύτητας τοιχωμάτων και κλίση πυθμένα . Για λόγους ελέγχου της ακρίβειας της αριθμητικής μεθόδου που χρησιμοποιήθηκε, αρχικά επιλύθηκε η περίπτωση τρισδιάστατου καναλιού ορθογωνικής διατομής χωρίς προβόλους και τα αποτελέσματα που προέκυψαν συγκρίθηκαν με αναλυτικά αποτελέσματα δισδιάστατης ροής (κατακόρυφο επίπεδο) υπεράνω επίπεδου πυθμένα. Τα αποτελέσματα βρέθηκαν σε καλή συμφωνία μεταξύ τους, γεγονός που επιβεβαίωσε την καταλληλότητα της μεθόδου. Για το τρισδιάστατο πρόβλημα με τους εγκάρσιους προβόλους, η ανάλυση έδειξε ότι το διάμηκες προφίλ της ελεύθερης επιφάνειας παρουσιάζει τοπική ταπείνωση στην θέση των προβόλων, ενώ η κατανομή των διατμητικών τάσεων παρουσιάζει μια σημαντική αύξηση στο μέτωπο των προβόλων και μια μείωση κατάντη των προβόλων στον πυθμένα και στο πλαϊνό τοίχωμα συγκριτικά με τις τιμές της διατμητικής τάσης στο κανάλι χωρίς τα εμπόδια. Επιπροσθέτως, παρατηρήθηκε η δημιουργία μιας περιοχής ανακυκλοφορίας της ροής κατάντη κάθε προβόλου. Η επανακόλληση της ροής κατάντη των προβόλων υπολογίστηκε ότι συμβαίνει σε μια απόσταση 11-19 φορές το μήκος του προβόλου πλησίον του πυθμένα ενώ πλησίον της ελεύθερης επιφάνειας αυτή η απόσταση ισούται με 6-7 φορές το μήκος του προβόλου ανάλογα με τον αριθμό και την απόσταση που απέχουν οι πρόβολοι μεταξύ τους. Αξίζει να σημειωθεί επίσης ότι οι τάσεις στο πλαϊνό τοίχωμα και στον πυθμένα στην περιοχή κατάντη των προβόλων παρουσίασαν μειωμένη κατά απόλυτο τιμή συγκριτικά με την περίπτωση του καναλιού χωρίς τα εμπόδια, γεγονός που συνηγορεί στη χρησιμοποίηση των προβόλων για προστασία διάβρωσης του ακτής. / This thesis deals with the analysis of turbulent flow in rectangular, three-dimensional open channel, in which inherent cantilever barriers (groynes). Such structures are used in many rivers to maintain the desired cross section and the desired depth of flow and to avoid corrosion of the banks to intense flooding.

Αντικείμενο της παρούσας διδακτορικής διατριβής αποτελεί η μελέτη και ανάπτυξη νέων μεθοδολογιών αυτόματης επεξεργασίας σεισμολογικών δεδομένων, µε σκοπό την επίλυση σημαντικών προβλημάτων που συναντώνται στα πεδία των επιστημών της σεισμολογίας και της γεωφυσικής όπως: 1) η ανίχνευση μικροσεισμικών γεγονότων από µία καταγραφή, µε άλλα λόγια ο διαχωρισμός της καταγραφής σε τμήματα που αποτελούνται από εδαφικό θόρυβο και σε τμήματα που περιέχουν την χρήσιμη πληροφορία (σεισμικά γεγονότα) για τους γεωεπιστήμονες και 2) η εκτίμηση των χρόνων άφιξης των διαμήκων (P-) καθώς και των εγκαρσίων (S-) σεισμικών φάσεων. Πιο αναλυτικά, η διατριβή είναι δομημένη ως εξής: Το πρώτο κεφάλαιο αποτελεί την εισαγωγή της διατριβής. Στο δεύτερο κεφάλαιο συγκεντρώνονται και κατηγοριοποιούνται όλες οι υπάρχουσες τεχνικές που έχουν αναπτυχθεί για την επίλυση του προβλήματος της αυτόματης ανίχνευσης σεισμικών γεγονότων καθώς και τον αυτόματο προσδιορισμό του χρόνου άφιξης των P και S σεισμικών φάσεων. Συγκεκριμένα γίνεται κατηγοριοποίηση αυτών σε τεχνικές που στηρίζονται στην ανάλυση και επεξεργασία των σεισμικών καταγραφών στα πεδία του χρόνου και της συχνότητας, στη χρήση νευρωνικών δικτύων, στην ανάλυση χρονικών σειρών και αυτοπαλινδρόμησης, στην ανάλυση της πόλωσης των κυμάτων, στις στατιστικές υψηλότερης τάξης, μεθόδους ασαφούς λογικής, κυματιδιακές μεθόδους κτλ. Στο τρίτο κεφάλαιο, αναπτύσσεται νέα τεχνική για την επίλυση του προβλήματος της αυτόματης ανίχνευσης σεισμικών γεγονότων από μία καταγραφή, η οποία βασίζεται σε μία μη αυστηρή διαδικασία ελέγχου υποθέσεων. Η προτεινόμενη τεχνική πραγματοποιείται σε δύο στάδια. Κατά το πρώτο στάδιο εκτιμώνται οι εμπειρικές συναρτήσεις πυκνότητας πιθανότητας που προκύπτουν τόσο από τον εδαφικό θόρυβο όσο και από τα υπόλοιπα που προέκυψαν από την λεύκανση αυτού. Κατά το δεύτερο στάδιο προτείνεται στατιστικό τεστ τύπου κατωφλίωσης για την αυτόματη ανίχνευση μικροσεισμικών γεγονότων. Η προτεινόμενη τεχνική εφαρμόζεται σε συνθετικά και πραγματικά δεδομένα και συγκρίνεται με τον γνωστό αλγόριθμο του λόγου βραχυπρόθεσμου προς μακροπρόθεσμο μέσο (STA/LTA). Στο τέταρτο κεφάλαιο παρουσιάζεται μέθοδος για την επίλυση του προβλήματος του αυτόματου προσδιορισμό του χρόνου άφιξης της P φάσης κάνοντας χρήση στατιστικών ανώτερης τάξης. Συγκεκριμένα, γίνεται χρήση των ποσοτήτων της λοξότητας, της κύρτωσης και μίας εκτίμησης της αντιεντροπίας ως γραμμικός συνδυασμός των παραπάνω. Επιπλέον παρουσιάζονται τα αποτελέσματα από την εφαρμογή της συγκεκριμένης τεχνικής σε συνθετικά αλλά και πραγματικά δεδομένα μικροσεισμικού δικτύου, κατάλληλα για χρήση στην παθητική σεισμική τομογραφία υψηλής ευκρίνειας. Τα αποτελέσματα αυτά συγκρίνονται με γνωστές ενεργειακές μεθόδους. Στο πέμπτο κεφάλαιο, αναπτύσσεται νέα τεχνική για την επίλυση του προβλήματος της αυτόματης εκτίμησης του χρόνου άφιξης της S φάσης. Η προτεινόμενη τεχνική βασίζεται στην στατιστική επεξεργασία συγκεκριμένης χαρακτηριστικής συνάρτησης, η οποία προκύπτει από τις ιδιότητες πόλωσης των σεισμικών κυμάτων που έχουν καταγραφεί. Επιπλέον, για να ελαττωθεί η εξάρτηση του προτεινόμενου αλγορίθμου από το χρησιμοποιούμενο παράθυρο, ακολουθείται μια πολυ-παραθυρική προσέγγιση του προβλήματος σε συνδυασμό με χρήση συναρτήσεων βαρών οι οποίες εκτιμώνται αυτόματα και βασίζονται στις μεταβολές της ενέργειας του σήματος κατά τη S άφιξη. Τέλος, παρουσιάζονται τα αποτελέσματα της εφαρμογής της μεθόδου σε πραγματικά δεδομένα καθώς και η αξιολόγησή τους σε περιβάλλον θορύβου. Στο έκτο κεφάλαιο, παρουσιάζονται τα αποτελέσματα της εφαρμογής των προτεινόμενων τεχνικών σε δεδομένα μικροσεισμικού δικτύου και συγκεκριμένα σε δεδομένα που προέκυψαν από πειράματα παθητικής σεισμικής τομογραφίας και τεχνητής υδραυλικής διάρρηξης που έλαβαν χώρα στην περιοχή Δέλβινο της ΝΔ Αλβανίας. Επιπλέον, γίνεται ανάλυση των αποτελεσμάτων βάσεις των δεικτών αβεβαιότητας που επέλεξαν οι αναλυτές στις εκτιμήσεις τους, καθώς και βάσει των λόγων σήματος θορύβου των καταγραφών. Στο έβδομο κεφάλαιο παρατίθενται τα συμπεράσματα της παρούσας διδακτορικής διατριβής, καθώς και πιθανές μελλοντικές προεκτάσεις. / The problems of seismic event detection and P- and S-phase arrival time estimation constitute important and vital tasks for the geoscientists. The solution of the aforementioned problems provides with important geophysical and seismological information, that can be used in a number of problems such as the structure of the earth’s interior, geotectonic settings, hypocentric and epicentric coordinates of an earthquake, the seismicity of an area and seismic hazard assessment. Traditionally, human experts have carried out this task. Nevertheless, during the last three decades due to the progress in computer technology, several methods have been developed for the automatic seismic event detection and P- and S- phase identification. After the introduction of the first chapter, in the second chapter the majority of the existing methods that have been developed and applied up to now, are gathered and categorized. These methods involve energy criteria, the seismic wave polarity assumption, artificial neural networks, higher order statistics, maximum likelihood methods, fuzzy logic methods etc. In the third chapter, a new thresholding type technique is proposed, tailored to fit real world situations where our knowledge on the statistical characteristics of the background noise process are unknown and a strict hypothesis testing framework can not be followed. In such cases the replacement of the unknown probability density function under the null hypothesis by its empirical counterpart, constitutes a possibility. In this work, a two stage procedure is proposed. The first one concerns the estimation of the empirical functions of the noise process itself as well as its whitened counterpart. In the second stage, using the above empirical functions, a thresholding scheme is proposed in order to solve the problem of the detection of seismic events in a non strict hypothesis testing framework. The performance of the proposed technique is confirmed by its application in a series of experiments both in synthetic and real seismic datasets. In the fourth chapter, the problem of automatic P-phase identification is solved using higher order statistics. The first- and second-order statistics (such as mean value, variance, autocorrelation, and power spectrum) are extensively used in signal processing to describe linear and Gaussian processes. In practice, many processes deviate from linearity and Gaussianity. Higher order statistics can be used for the study of such processes. The P-phase arrival time is estimated using these HOS parameters and additionally, an estimation of the negentropy defined as a linear combination of skewness and kurtosis. According to the implemented algorithm a moving window “slides” on the recorded signal, estimating skewness, kurtosis, and negentropy. Skewness can be considered as a measure of symmetry of the distribution, while kurtosis is a measure of heaviness of the tails, so they are suitable for detecting parts of the signal that do not follow the amplitude distribution of ambient noise. Seismic events have higher amplitudes in comparison to the seismic noise, and these higher values occupy the tails of the distribution (high degree of asymmetry of distribution). In the case of seismic events, skewness and kurtosis obtain high values, presenting maxima in the transition from ambient noise to the seismic events (P-arrival). The proposed algorithms are applied on synthetic as well as real seismic data and compared to well known energy based methods. Algorithms that deal with the automatic S-onset time identification problem, is a topic of ongoing research. Modern dense seismic networks used for earthquake location, seismic tomography investigations, source studies, early warning etc., demand accurate automatic S-wave picking. Most of the techniques that have been proposed up to now are mainly based on the polarization features of the seismic waves. In the fifth chapter, a new time domain method for the automatic determination of the S-phase arrival onsets is proposed and its implementation on local earthquake data is presented. Eigevalue analysis is taking place over small time intervals, and the maximum eigenvalue which is obtained on each step is retained for further processing. In this way a time series of maximum eigenvalues is formed, which serves as a characteristic function. A first S-phase arrival time estimation is obtained by applying the kurtosis criterion on the derived characteristic function. Furthermore, a multi-window approach combined with an energy-based weighting scheme is also applied, in order to reduce the algorithm’s dependence on the moving window’s length and provide a weighted S phase onset. Automatic picks are compared against manual reference picks and moreover the proposed technique is subjected to a noise robustness test. In the sixth chapter, the results of the implementation of the proposed techniques on microseismic data are presented. Specifically, the proposed methods are applied on two real sets of data. One dataset was been recorded during a Passive Seismic Tomography (PST) experiment, while the second one during the seismic monitoring of fracking operations. Both experiments took place in a hydrocarbon field in Delvina, SW Albania. These results are also analyzed, based on the arrival times and their uncertainty as they were evaluated by human analysts as well as the corresponding signal to noise ratio of the seismic records. Finally, the seventh chapter concludes this work and possible future extensions are discussed.

Στην περιοχή του ΒΔ Λασιθίου απαντώνται οι γεωλογικοί σχηματισμοί των πλακώδων ασβεστόλιθων και των ασβεστόλιθων της ζώνης Τρίπολης. οι σχηματισμοί αυτοί είναι έντονα καρστικοποιημένοι και αποτελούν τους κύριους υδροφόρους της περιοχής έρευνας. Μεταξύ των δύο υδροφόρων, στις περιοχές Νεάπολης και Ποτάμων παρεμβάλλονται τα αδιαπέρατα στρώματα των Φυλλιτών-Χαλαζιτών. Ο υδροφόρος που αναπτύσσεται στους ανθρακικούς σχηματισμούς της ζώνης Τρίπολης κατηγοριοποιείται σε τρεις ζώνες σύμφωνα με τον υδροχημικό χαρακτήρα του υπόγειου νερού. η πρώτη ζώνη αποτελείται από τις περιοχές τροφοδοσίας του υδροφόρου (Δράσι, Άγιος Κωνσταντίνος και Ποτάμοι) με εξαίρετη ποιότητα νερού, η δεύτερη ζώνη συνιστάται από τις περιοχές των Λακωνίων και Κριτσάς, όπου το νερό του υδροφόρου έχει επηρεαστεί από τη διείσδυση του θαλασσινού νερού. Τέλος, η τρίτη ζώνη αποτελείται από τις παράκτιες περιοχές του Αγίου Νικολάου και Αλμυρού με έντονη επιρροή από την διαδικασία της υφαλμύρινσης Ο υδροφόρος που αναπτύσσεται στους Πλακώδεις ασβεστόλιθους είναι ανοικτός στην θάλασσα και επηρεάζεται στο σύνολο του από την διείσδυση του θαλασσινού νερού. Με τη χρήση των κανονικοποιημένων διαγραμμάτων των σπάνιων γαιών, τόσο σε δείγματα πετρωμάτων όσο και σε δείγματα νερού προέκυψε ότι οι δύο καρστικοί υδροφόροι επικοινωνούν υδραυλικά στην περιοχή του Δάμακα Β-ΒΔ του Αγίου Νικολάου. / The region of NW Lasithi is dominated by plattenkalk and Tripoli’s zone limestones. These geological formations are characterized by karst and host the main aquifers of study area. Impermeable layers of phyllite-quartzite formation occur between these limestone formations, in the areas of Neapolis and Potami. The aquifer hosted by the carbonate formations of Tripoli’s zone is divided in three zones based on the hydrochemical characteristics of the underground water. The first zone includes the areas Drasi, Agios Konstantinos and Potami that are the sources of the aquifer, with excellent water quality. The second zone consists of the areas Lakonia and Kritsa, in which the underground water have been infiltrated by sea water. Finally the third consists of the coastal zones of Agios Nikolaos and Almyros, which have been affected by the process of brackish water. The aquifer hosted by the plattenkalk limestones is open to the sea and is affected by the intrusion of sea water. It can be shown, using normalized rare earth element diagrams of analyses of both rocks and underground waters, that the two karstified aquifers communicate hydraulically in the area of Damakas, N-NW of city Agios Nikolaos.