Localization with multicomponent seismic array

International audienceSeismo-volcano source localization is essential to improve our understanding of volcano systems. The lack of clear seismic wave phases prohibits the use of classical location methods. Seismic antennas composed of one-component (1C) seismometers provide a good estimate of the back-azimuth of the wavefield. The depth estimation, on the other hand, is difficult or impossible to determine. In order to determine the source location parameters (back-azimuth and depth), we extend the 1C seismic antenna approach to 3Cs. This communication discusses a high resolution location method using a 3C array survey (3C-MUSIC algorithm) with data from two seismic antennas installed on an andesitic volcano in Peru (Ubinas volcano). After introducing the 3C MUSIC processing, we evaluate the robustness of the location method on a full wavefield 3D synthetic dataset generated using a digital elevation model of Ubinas volcano and an homogeneous velocity model. Results show that the back-azimuth determined using the 3C array has a smaller error than a 1C array. Only the 3C method allows the recovery of the source depths. Finally, we applied the 3C-MUSIC to two seismic events recorded in 2009. Therefore, extending 1C arrays to 3C arrays in volcano monitoring allows a more accurate determination of the source epicenter and now an estimate for the depth

Seismo-volcano source localization with triaxial broad-band seismic array

International audienceSeismo-volcano source localization is essential to improve our understanding of eruptive dynamics and of magmatic systems. The lack of clear seismic wave phases prohibits the use of classical location methods. Seismic antennas composed of one-component (1C) seismometers provide a good estimate of the backazimuth of the wavefield. The depth estimation, on the other hand, is difficult or impossible to determine. As in classical seismology, the use of three-component (3C) seismometers is now common in volcano studies. To determine the source location parameters (backazimuth and depth), we extend the 1C seismic antenna approach to 3Cs. This paper discusses a high-resolution location method using a 3C array survey (3CMUSIC algorithm) with data from two seismic antennas installed on an andesitic volcano in Peru (Ubinas volcano). One of the main scientific questions related to the eruptive process of Ubinas volcano is the relationship between the magmatic explosions and long-period (LP) swarms. After introducing the 3C array theory, we evaluate the robustness of the location method on a full wavefield 3-D synthetic data set generated using a digital elevation model of Ubinas volcano and an homogeneous velocity model. Results show that the backazimuth determined using the 3C array has a smaller error than a 1C array. Only the 3C method allows the recovery of the source depths. Finally, we applied the 3C approach to two seismic events recorded in 2009. Crossing the estimated backazimuth and incidence angles, we find sources located 1000 ± 660 m and 3000 ± 730 m below the bottom of the active crater for the explosion and the LP event, respectively. Therefore, extending 1C arrays to 3C arrays in volcano monitoring allows a more accurate determination of the source epicentre and now an estimate for the depth

Modulated Oscillations in Three Dimensions

The analysis of the fully three-dimensional and time-varying polarization characteristics of a modulated trivariate, or three-component, oscillation is addressed. The use of the analytic operator enables the instantaneous three-dimensional polarization state of any square-integrable trivariate signal to be uniquely defined. Straightforward expressions are given which permit the ellipse parameters to be recovered from data. The notions of instantaneous frequency and instantaneous bandwidth, generalized to the trivariate case, are related to variations in the ellipse properties. Rates of change of the ellipse parameters are found to be intimately linked to the first few moments of the signal's spectrum, averaged over the three signal components. In particular, the trivariate instantaneous bandwidth---a measure of the instantaneous departure of the signal from a single pure sinusoidal oscillation---is found to contain five contributions: three essentially two-dimensional effects due to the motion of the ellipse within a fixed plane, and two effects due to the motion of the plane containing the ellipse. The resulting analysis method is an informative means of describing nonstationary trivariate signals, as is illustrated with an application to a seismic record.Comment: IEEE Transactions on Signal Processing, 201

Geophysical methods to detect tunnelling at a geological repository site : Applicability in safeguards

ABSTRACT Generating power with nuclear energy accumulates radioactive spent nuclear fuel, anticipated not to be diversified into any unknown purposes. Nuclear safeguards include bookkeeping of nuclear fuel inventories, frequent checking, and monitoring to confirm nuclear non-proliferation. Permanent isolation of radionuclides from biosphere by disposal challenges established practices, as opportunities for monitoring of individual fuel assemblies ceases. Different concepts for treatment and geological disposal of spent nuclear fuel exist. Spent nuclear fuel disposal facility is under construction in Olkiluoto in Southwest Finland. Posiva Oy has carried out multidisciplinary bedrock characterization of crystalline bedrock for siting and design of the facility. Site description involved compilation of geological models from investigations at surface level, from drillholes and from underground rock characterization facility ONKALO. Research focused on long term safety case (performance) of engineered and natural barriers in purpose to minimize risks of radionuclide release. Nuclear safeguards include several concepts. Containment and surveillance (C/S) are tracking presence of nuclear fuel through manufacturing, energy generation, cooling, transfer, and encapsulation. Continuity of knowledge (CoK) ensures traceability and non-diversion. Design information provided by the operator to the state and European Commission (Euratom), and further to IAEA describes spent nuclear fuel handling in the facility. Design information verification (DIV) using timely or unannounced inspections, provide credible assurance on absence of any ongoing undeclared activities within the disposal facility. Safeguards by design provide information applicable for the planning of safeguards measures, e.g., surveillance during operation of disposal facility. Probability of detection of an attempt to any undeclared intrusion into the repository containment needs to be high. Detection of such preparations after site closure would require long term monitoring or repeated geophysical measurements within or at proximity of the repository. Bedrock imaging (remote sensing, geophysical surveys) would serve for verifying declarations where applicable, or for characterization of surrounding rock mass to detect undeclared activities. ASTOR working group has considered ground penetrating radar (GPR) for DIV in underground constructed premises during operation. Seismic reflection survey and electrical or electromagnetic imaging may also apply. This report summarizes geophysical methods used in Olkiluoto, and some recent development, from which findings could be applied also for nuclear safeguards. In this report the geophysical source fields, involved physical properties, range of detection, resolution, survey geometries, and timing of measurements are reviewed for different survey methods. Useful interpretation of geophysical data may rely on comparison of results to declared repository layout, since independent understanding of the results may not be successful. Monitoring provided by an operator may enable alarm and localization of an undeclared activity in a cost-effective manner until closure of the site. Direct detection of constructed spaces, though possible, might require repeated effort, have difficulties to provide spatial coverage, and involve false positive alarms still requiring further inspection

On Some DOA Algorithms for Tri-axial Geophone

In this paper a short study of some basic methods of DOA of a seismic wave using so called tri-axial geophone has been presented. The proposed methods exploit the properties of Rayleigh surface plane wave to find DOA of an incoming seismic wave using inner products of appropriately filtered signals recorded by geophones. The advantage of the proposed method is its simplicity and ease of implementation in small DSP or application processors still retaining pretty good accuracy. A number of example results for real data have been given

The study of local sources of ULF geoelectric signals with steep fronts

The method for recognition and monitoring of ULF electrotelluric signals with steep fronts is presented. The method is based on the multicomponent measurements of the geoelectric field and recognition of the clusters of signals using the data of spatial-polarization filtering. This detects the local sources of the signals associated with the variation of mode of deformation of the medium. The method was applied to the study of geoelectric sources associated with aftershocks of the Racha-Djava earthquake. The field observations show with a fair degree of confidence that these signals are generated by a small number of local zones in the vicinity of the observer (at a range of tens-hundreds meters). The spatial distribution and the dynamic properties of the local sources are described. The analysis of observations leads to the conclusion that the dynamics of the local sources of ULF signals with steep fronts may be considered possible electrotelluric precursors of the earthquakes. Consideration is given to the feasibility of the well-known VAN method for earthquake prediction. The overall methodological aspects of the current studies of the electrotelluric precursors of the earthquakes are discussed

Joint detection and classification of rockfalls in a microseismic monitoring network

A rockfall (RF) is a ubiquitous geohazard that is difficult to monitor or predict and poses a significant risk for people and transportation in several hilly and mountainous environments. The seismic signal generated by RF carries abundant physical and mechanical information. Thus, signals can be used by researchers to reconstruct the event location, onset time, volume and trajectory, and develop an efficient early warning system. Therefore, the precise automatic detection and classification of RF events are important objectives for scientists, especially in seismic monitoring arrays. An algorithm called DESTRO (DEtection and STorage of ROckfalls) aimed at combining seismic event automatic detection and classification was implemented ad hoc within the MATLAB environment. In event detection, the STA/LTA (short-time-average through long-time-average) method combined with other parameters, such as the minimum duration of an RF and the minimum interval time between two continuous seismic events is used. Furthermore, nine significant features based on the frequency, amplitude, seismic waveform, duration and multiple station attributes are newly proposed to classify seismic events in a RF environment. In particular, a three-step classification method is proposed for the discrimination of five different source types: RFs, earthquakes (EQs), tremors, multispike events (MSs) and subordinate MS events. Each component (vertical, east–west and north–south) at each station within the monitoring network is analysed, and a three-step classification is performed. At a given time, the event series detected from each component are integrated and reclassified component by component and station by station into a final event-type series as an output result. By this algorithm, a case study of the seven-month-long seismic monitoring of a former quarry in Central Italy was investigated by means of four triaxial velocimeters with continuous acquisition at a sampling rate of 200 Hz. During this monitoring period, a human-induced RF simulation was performed, releasing 95 blocks (in which 90 blocks validated) of different sizes from the benches of the quarry. Consequently, 64.9 per cent of EQs within 100 km were confirmed in a one-month monitoring period, 88 blocks in the RF simulation were classified correctly as RF events and 2 blocks were classified as MSs given their small energy. Finally, an ad hoc section of the algorithm was designed specifically for RF classification combined with EQ recognition. The algorithm could be applied in slope seismic monitoring to monitor the dynamic states of rock masses, as well as in slope instability forecasting and risk evaluation in EQ-prone areas

Understanding magmatic processes and seismo-volcano source localization with multicomponent seismic arrays

Dans cette thèse, nous étudions le problème de la localisation de sources sismo-volcanique, à partir des données enregistrées par des réseaux de capteurs composés de nouveaux sismomètres à trois composantes (3C). Nous nous concentrerons sur le volcan Ubinas, l'un des plus actifs au Pérou. Nous développons une nouvelle approche (MUSIC-3C) basée sur la méthode MUSIC permetant de retourner les 3 paramètres utiles (lenteur, azimut et incidence). Pour valider notre méthodologie, nous analysons des sources synthétiques propagées en tenant compte de la topographie du volcan Ubinas. Dans cette expérience, les données synthétiques ont été générées pour plusieurs sources situées à différentes profondeurs sous le cratère Ubinas. Nous utilisons l'algorithme MUSIC-3C pour les relocaliser. Nous traitons également des données réelles provenant d'une expérience de terrain menée sur le volcan Ubinas (Pérou) en 2009 par les équipes de recherche de l'IRD-France (Institut de Recherche pour le Déveleppment), UCD l'Irlande (projet VOLUME) et l'Institut de Géophysique du Pérou (IGP). Nous utilisons l'algorithme MUSIC-3C pour localiser les événements explosifs (type vulcanien), ce qui nous permet d'identifier et d'analyser les processus physiques de ces événements, à la suite de cette analyse, nous avons trouvé deux sources pour chaque explosion situées à 300 m et 1100 m en dessous du fond du cratère actif. Basé sur les mécanismes éruptifs proposés pour d'autres volcans du même type, nous interprétons la position de ces sources ainsi que les limites du conduit éruptif impliqué dans le processus de fragmentation.In this thesis, we study the seismo-volcanic source localization using data recorded by new sensor arrays composed of three-component (3C) seismometers deployed on Ubinas stratovolcano (Peru). We develop a new framework (MUSIC-3C) of source localization method based on the well-known MUSIC algorithm. To investigate the performance of the MUSIC-3C method, we use synthetic datasets designed from eight broadband isotropic seismic sources located beneath the crater floor at different depths. The fundamental scheme of the MUSIC-3C method exploits the fact of the cross-spectral matrix of 3C array data, corresponding to the first seismic signal arrivals, provides of useful vector components (slowness, back-azimuth and incidence angle) from the seismic source. Application of the MUSIC-3C method on synthetic datasets shows the recovery of source positions. Real data used in this study was collected during seismic measurements with two seismic antennas deployed at Ubinas volcano in 2009, whose experiment conduced by volcanic teams of IRD-France (l'Institute de Recherche pour le Déveleppment), Geophysics group University College Dublin Ireland and Geophysical Institute of Peru (IGP). We apply the MUSIC-3C algorithm to investigate wave fields associated with the magmatic activity of Ubinas volcano. These analysis evidence a complex mechanism of vulcanian eruptions in which their seismic sources are found at two separated sources located at depths of 300 m and 1100 m beneath the crater floor. This implies the reproduction of similar mechanisms into the conduit. Based on the eruptive mechanisms proposed for other volcanoes of the same type, we interpret the position of this sources as the limits of the conduit portion that was involved in the fragmentation process.SAVOIE-SCD - Bib.électronique (730659901) / SudocGRENOBLE1/INP-Bib.électronique (384210012) / SudocGRENOBLE2/3-Bib.électronique (384219901) / SudocSudocFranceF