Amoroso, Ortensia

Sub surface operations for energy production such as gas storage, fluid injection or hydraulic fracking modify the physical properties of the crust, in particular seismic velocity and anelastic attenuation. Continuously measuring these properties may be crucial to monitor the status of the reservoir. Here we propose a not usual use of the empirical ground-motion prediction equations (GMPEs) to monitor large-scale medium properties variations in a reservoir during fluid injection experiments. In practice, peak-ground velocities recorded during field operations are used to update the coefficients of a reference GMPE whose variation can be physically interpreted in terms of anelastic attenuation and seismic velocity. We apply the technique to earthquakes recorded at The Geysers geothermal field in Southern California and events occurred in the St. Gallen (Switzerland) geothermal field. Our results suggest that the GMPEs can be effectively used as a proxy for some reservoir properties variation by using induced earthquakes recorded at relatively dense networks.

We present a detailed seismic imaging of the harnessed Nesjavellir geothermal area, SW-Iceland, which is one of several geothermal fields on the flanks of the Hengill volcano. We map the vP , vS , and vPvS ratio using seismic data recorded in 2016–2020 and compare them with both a resistivity model of the same area and the rock temperature as measured in boreholes. The results show that the shallower crust (depth less than 1 km) is characterized by low vP and vS , and high vPvS ratio (around 1.9). Shallow low resistivity values at similar depths in the same area have been interpreted as the smectite clay cap of the system. At depths between 1 and 3 km the observed low vPvS ratio of 1.64–1.70 correlates with high resistivity values. In this area, characterized by temperatures larger than 240°C, the seismicity appears to be sparse and located close to the production wells. This seismicity has been interpreted as induced by the production in combination with naturally occurring earthquakes. At depths greater than 4 km, high vPvS ratio of 1.9 correlates well with low resistivity values. In the valley of Nesjavellir, a deep-seated conductive body, domes up at about 4.500 m b.sl. and coincides spatially with a significant high vPvS ratio anomaly (>1.9). Above these anomalies an elevated temperature is registered according to borehole temperature data. This is proposed here to be caused by hot 600°C–900°C cooling intrusives, close to the brittle ductile transition—probably the heat source(s) of the geothermal field above. These anomalies are at the same location as the last fissure eruption in Hengill almost 2,000 years ago. The NNE-SSW trending, deeper seismic cluster at 3–6 km depth is located at the edge of this high vPvS anomaly. The heat source of the Nesjavellir geothermal field is most likely connected to this most recent volcanism as reflected by the deep-seated low resistivity body and high vPvS ratio, located beneath the deep fault that connects the flow path of the high temperature geothermal fluid, resulting in an actively producing reservoir.

A tomographic analysis of Mt. Pollino area (Italy) has been performed using earthquakes recorded in the area during an intense seismic sequence that occurred between 2010 and 2014. 870 local earthquakes with magnitude ranging from 1.8 to 5.0 were selected considering the number of recording stations, the signal quality, and the hypocenter distribution. P- and S-wave arrival times were manually picked and used to compute 3D velocity models through tomographic seismic inversion. The resulting 3D distributions of VP and VS are characterized by high resolution in the central part of the investigated area and from surface to about 10 km below sea level. The aim of the work is to obtain high- quality tomographic images to correlate with the main lithological units that characterize the study area. The results will be important to enhance the seismic hazard assessment of this complex tectonic region. These images show the ductile Apennine platform (VP = 5.3 km/s) overlaying the brittle Apulian platform (VP=6.0 km/s) at depth of around 5 km. The central sector of the area shows a clear fold and thrust interface. Along this structure,most of the seismicity occurred, including the strongest event of the sequence (M W 5.0). High V P (>6.8 km/s) and high V P /VS (>1.9) patterns, intersecting the southern edge of this western seismogenic volume, have been interpreted as water saturated rocks, in agreement with similar geological context in the Apennines. These fluids could have played a role in nucleation and development of the seismic sequence. A recent study revealed the occurrence of clusters of earthquakes with similar waveforms along the same seismogenic volume. The hypocenters of these cluster events have been compared with the events re-located in this work. Jointly, they depict a 10 km × 4 km fault plane, NW-SE oriented, deepening towards SW with a dip angle of 40–45° . Instead, the volume of seismicity responsible for the M L 4.3 earthquake developed as a mainshock-aftershock sequence, occurring entirely within the average-to-low VP /VS Apennine platform. Our results agree with other independent geophysical analyses carried out in this area, and they could significantly improve the actual knowledge of the main lithologic units of this complex tectonic area.

Geo-resources are widely exploited in our society, with huge benefits for both economy and communities. Nevertheless, with benefits come risks and impacts. Understanding how such risks and impacts are intrinsically borne in a given project is of critical importance for both industry and society. In particular, it is crucial to distinguish between the specific impacts related to exploiting a given energy resource and those shared with the exploitation of other energy resources. A variety of different approaches can be used to identify and assess such risks and impacts. In particular, Life Cycle Assessment (LCA) and risk assessments (RAs) are the most commonly adopted. Although both are widely used to support decision making in environmental management, they are rarely used in combination perhaps because they have been developed by largely different groups of specialists. By analyzing the structure and the ratio of the two tools, we have developed an approach for combining and harmonizing LCA and MRA; the resulting protocol envisages building MRA upon LCA both qualitatively and quantitatively. We demonstrate the approach in a case study using a virtual site (based on a real one) for geothermal energy production

Focal mechanisms of selected earthquakes, recorded in the Mount Pollino region (southern Italy) from 2010 through 2014, are used to infer the pore fluid pressure at hypocenter depths. The 3-D excess pore pressure field provides evidence that the sequence occurs in a fluid-filled volume with values reaching 35 MPa. The mechanisms underlying this swarm-like sequence and the triggering of earthquakes are investigated by computing the cumulative static Coulomb stress change at hypocenter depths and analyzing the pore-pressure diffusion mechanism. The results indicate that static Coulomb stress change was lower than 0.01 MPa, which is the value generally assumed as threshold for the triggering, and seismicity distribution was actually driven by pore-pressure diffusion with relatively low diffusivity value. This latter mechanism could also explain the delayed triggering of the two larger events ML 4.3 and ML 5.0, respectively, that occurred about 150 days apart.

Reliable seismic hazard analyses are crucial to mitigate seismic risk. When dealing with induced seismicity the standard Probabilistic Seismic Hazard Analysis (PSHA) has to be modified because of the peculiar characteristics of the induced events. In particular, the relative shallow depths, small magnitude, a correlation with field operations, and eventually non-Poisson recurrence time. In addition to the well-known problem of estimating the maximum expected magnitude, it is important to take into account how the industrial field operations affect the temporal and spatial distribution of the earthquakes. In fact, during specific stages of the project the seismicity may be hard to be modelled as a Poisson process—as usually done in the standard PSHA—and can cluster near the well or migrate toward hazardous known or—even worse—not known faults. Here we present a technique in which we modify the standard PSHA to compute time-dependent seismic hazard. The technique allows using non-Poisson models (BPT, Weibull, gamma and ETAS) whose parameters are fitted using the seismicity record during distinct stages of the field operations. As a test case, the procedure has been implemented by using data recorded at St. Gallen deep geothermal field, Switzerland, during fluid injection. The results suggest that seismic hazard analyses, using appropriate recurrence model, ground motion predictive equations, and maximum magnitude allow the expected ground-motion to be reliably predicted in the study area. The predictions can support site managers to decide how to proceed with the project avoiding adverse consequences.