Integrating new and traditional approaches for the estimate of slip-rates of active faults: examples from the Mw 6.3, 2009 L’Aquila earthquake area, Central Italy

This thesis is based on the integration of traditional and innovative approaches aimed at improving the normal faults seimogenic identification and characterization, focusing mainly on slip-rate estimate as a measure of the fault activity. The L’Aquila Mw 6.3 April 6, 2009 earthquake causative fault, namely the Paganica - San Demetrio fault system (PSDFS), was used as a test site. We developed a multidisciplinary and scale‐based strategy consisting of paleoseismological investigations, detailed geomorphological and geological field studies, as well as shallow geophysical imaging and an innovative application of physical properties measurements. We produced a detailed geomorphological and geological map of the PSDFS, defining its tectonic style, arrangement, kinematics, extent, geometry and internal complexities. The PSDFS is a 19 km-long tectonic structure, characterized by a complex structural setting and arranged in two main sectors: the Paganica sector to the NW, characterized by a narrow deformation zone, and the San Demetrio sector to SE, where the strain is accommodated by several tectonic structures, exhuming and dissecting a wide Quaternary basin, suggesting the occurrence of strain migration through time. The integration of all the fault displacement data and age constraints (radiocarbon dating, optically stimulated luminescence (OSL) and tephrochronology) helped in calculating an average Quaternary slip-rate representative for the PSDFS of 0.27 - 0.48 mm/yr. On the basis of its length (ca. 20 km) and slip per event (up to 0.8 m) we also estimated a max expected Magnitude of 6.3-6.8 for this fault. All these topics have a significant implication in terms of surface faulting hazard in the area and may contribute also to the understanding of the PSDFS seismic behavior and of the local seismic hazard.

ATP modeling of internal transformer faults for relay performance testing

Transformers are very important elements of any power system. Unfortunately, they are subjected to through-faults and abnormal operating conditions which can affect not only the transformer itself but also other equipment connected to the transformer. Thus, it is essential to provide sufficient protection for transformers as well as the best possible selectivity and sensitivity of the protection. Nowadays microprocessor-based relays are widely used to protect power equipment. Current differential and voltage protection strategies are used in transformer protection applications and provide fast and sensitive multi-level protection and monitoring. The elements responsible for detecting turn-to-turn and turn-to-ground faults are the negative-sequence percentage differential element and restricted earth-fault (REF) element, respectively. During severe internal faults current transformers can saturate and slow down the speed of relay operation which affects the degree of equipment damage. The scope of this work is to develop a modeling methodology to perform simulations and laboratory tests for internal faults such as turn-to-turn and turn-to-ground for two step-down power transformers with capacity ratings of 11.2 MVA and 290 MVA. The simulated current waveforms are injected to a microprocessor relay to check its sensitivity for these internal faults. Saturation of current transformers is also studied in this work. All simulations are performed with the Alternative Transients Program (ATP) utilizing the internal fault model for three-phase two-winding transformers. The tested microprocessor relay is the SEL-487E current differential and voltage protection relay. The results showed that the ATP internal fault model can be used for testing microprocessor relays for any percentage of turns involved in an internal fault. An interesting observation from the experiments was that the SEL-487E relay is more sensitive to turn-to-turn faults than advertized for the transformers studied. The sensitivity of the restricted earth-fault element was confirmed. CT saturation cases showed that low accuracy CTs can be saturated with a high percentage of turn-to-turn faults, where the CT burden will affect the extent of saturation. Recommendations for future work include more accurate simulation of internal faults, transformer energization inrush, and other scenarios involving core saturation, using the newest version of the internal fault model. The SEL-487E relay or other microprocessor relays should again be tested for performance. Also, application of a grounding bank to the delta-connected side of a transformer will increase the zone of protection and relay performance can be tested for internal ground faults on both sides of a transformer.

Transformer modeling ATP : internal faults & high-frequency discretization

Power transformers are key components of the power grid and are also one of the most subjected to a variety of power system transients. The failure of a large transformer can cause severe monetary losses to a utility, thus adequate protection schemes are of great importance to avoid transformer damage and maximize the continuity of service. Computer modeling can be used as an efficient tool to improve the reliability of a transformer protective relay application. Unfortunately, transformer models presently available in commercial software lack completeness in the representation of several aspects such as internal winding faults, which is a common cause of transformer failure. It is also important to adequately represent the transformer at frequencies higher than the power frequency for a more accurate simulation of switching transients since these are a well known cause for the unwanted tripping of protective relays. This work develops new capabilities for the Hybrid Transformer Model (XFMR) implemented in ATPDraw to allow the representation of internal winding faults and slow-front transients up to 10 kHz. The new model can be developed using any of two sources of information: 1) test report data and 2) design data. When only test-report data is available, a higher-order leakage inductance matrix is created from standard measurements. If design information is available, a Finite Element Model is created to calculate the leakage parameters for the higher-order model. An analytical model is also implemented as an alternative to FEM modeling. Measurements on 15-kVA 240?/208Y V and 500-kVA 11430Y/235Y V distribution transformers were performed to validate the model. A transformer model that is valid for simulations for frequencies above the power frequency was developed after continuing the division of windings into multiple sections and including a higher-order capacitance matrix. Frequency-scan laboratory measurements were used to benchmark the simulations. Finally, a stability analysis of the higher-order model was made by analyzing the trapezoidal rule for numerical integration as used in ATP. Numerical damping was also added to suppress oscillations locally when discontinuities occurred in the solution. A maximum error magnitude of 7.84% was encountered in the simulated currents for different turn-to-ground and turn-to-turn faults. The FEM approach provided the most accurate means to determine the leakage parameters for the ATP model. The higher-order model was found to reproduce the short-circuit impedance acceptably up to about 10 kHz and the behavior at the first anti-resonant frequency was better matched with the measurements.

Microfabric memory of vein quartz for strain localization in detachment faults: A case study on the Simplon fault zone

This manuscript deals with the adaptation of quartz-microfabrics to changing physical deformation conditions, and discusses their preservation potential during subsequent retrograde deformation. Using microstructural analysis, a sequence of recrystallization processes in quartz, ranging from Grain-Boundary Migration Recrystallization (GBM) over Subgrain-Rotation Recrystallization (SGR) to Bulging Nucleation (BLG) is detected for the Simplon fault zone (SFZ) from the low strain rim towards the internal high strain part of the large-scale shear zone. Based on: (i) the retrograde cooling path; (ii) estimates of deformation temperatures; and (iii) spatial variation of dynamic recrystallization processes and different microstructural characteristics, continuous strain localization with decreasing temperature is inferred. In contrast to the recrystallization microstructures, crystallographic preferred orientations (CPO) have a longer memory. CPO patterns indicative of prism and rhomb glide systems in mylonitic quartz veins, overprinted at low temperatures (�400 �C), suggest inheritance of a high-temperature deformation. In this way, microstructural, textural and geochemical analyses provide information for several million years of the deformation history. The reasons for such incomplete resetting of the rock texture is that strain localization is caused by change in effective viscosity contrasts related to temporal large- and small-scale temperature changes during the evolution of such a long-lived shear zone. The spatially resolved, quantitative investigation of quartz microfabrics and associated recrystallization processes therefore provide great potential for an improved understanding of the geodynamics of large-scale shear zones.

Long-term fluid circulation in extensional faults in the central Catalan Coastal Ranges: P-T constraints from neoformed chlorite and K-white mica

The neoformation of chlorite and K-white mica in fault rocks from two main faults of the central Catalan Coastal Ranges, the Vallès and the Hospital faults, has allowed us to constrain the P–T conditions during fault evolution using thermodynamic modeling. Crystallization of M1 and M2 muscovite and microcline occured as result of deuteric alteration during the exhumation of the pluton (290 °C > T > 370 °C) in the Permian. After that, three tectonic events have been distinguished. The first tectonic event, attributed to the Mesozoic rifting, is characterized by precipitation of M3 and M4 phengite together with chlorite and calcite C1 at temperatures between 190 and 310 °C. The second tectonic event attributed to the Paleogene compression has only been identified in the Hospital fault with precipitation of low-temperature calcite C2. The shortcut produced during inversion of the Vallès fault was probably the responsible for the lack of neoformed minerals within this fault. Finally, the third tectonic event, which is related to the Neogene extension, is characterized in the Vallès fault by a new generation of chlorite, associated with calcite C4 and laumontite, formed at temperatures between 125 and 190 °C in the absence of K-white mica. Differently, the Hospital fault is characterized by the precipitation of calcite C3 during the syn-rift stage at temperatures around 150 °C and by low-temperature fluids precipitating calcites C5, C6 and PC1 during the post-rift stage. During the two extensional events (Mesozoic and Neogene), faults acted as conduits for hot fluids producing anomalous high geothermal gradients (50 °C/km minimum).

Strike-slip faults mediate the rise of crustal-derived fluids and mud volcanism in the deep sea

We report on newly discovered mud volcanoes located at about 4500 m water depth 90 km west of the deformation front of the accretionary wedge of the Gulf of Cadiz, and thus outside of their typical geotectonic environment. Seismic data suggest that fluid flow is mediated by a >400-km-long strike-slip fault marking the transcurrent plate boundary between Africa and Eurasia. Geochemical data (Cl, B, Sr, 87Sr/86Sr, Delta18O, DeltaD) reveal that fluids originate in oceanic crust older than 140 Ma. On their rise to the surface, these fluids receive strong geochemical signals from recrystallization of Upper Jurassic carbonates and clay-mineral dehydration in younger terrigeneous units. At present, reports of mud volcanoes in similar deep-sea settings are rare, but given that the large area of transform-type plate boundaries has been barely investigated, such pathways of fluid discharge may provide an important, yet unappreciated link between the deeply buried oceanic crust and the deep ocean.

Detection and Isolation of Simultaneous Additive and Parametric Faults in Nonlinear Stochastic Dynamical Systems

This paper presents a new fault detection and isolation scheme for dealing with simultaneous additive and parametric faults. The new design integrates a system for additive fault detection based on Castillo and Zufiria, 2009 and a new parametric fault detection and isolation scheme inspired in Munz and Zufiria, 2008 . It is shown that the so far existing schemes do not behave correctly when both additive and parametric faults occur simultaneously; to solve the problem a new integrated scheme is proposed. Computer simulation results are presented to confirm the theoretical studies.

New selective earth faults only current directional method for isolated neutral systems

In order to achieve total selectivity at electrical distribution networks it is of great importance to analyze the defect currents at ungrounded power systems. This information will help to grant selectivity at electrical distribution networks ensuring that only the defect line or feeder is removed from service. In the present work a new selective and directional protection method for ungrounded power systems is evaluated. The new method measures only defect currents to detect earth faults and works with a directional criterion to determine the line under faulty conditions. The main contribution of this new technique is that it can detect earth faults in outgoing lines at any type of substation avoiding the possible mismatch of traditional directional earth fault relays. This detection technique is based on the comparison of the direction of a reference current to the direction of all earth fault capacitive currents at all the feeders connected to the same bus bars. This new method has been validated through computer simulations. The results for the different cases studied are remarkable, proving total validity and usefulness of the new method.