Synthetic Seismic Reflection Modelling in a Supercritical Geothermal System: An Image of the K-Horizon in the Larderello Field (Italy)

• Published in: Geofluids Vol. 2019; no. 2019; pp. 1 - 21
• Main Authors: Manzella, Adele, Montegrossi, G., Caielli, G., Scrocca, Davide, Petracchini, L., de Franco, R., Santilano, Alessandro
• Format: Journal Article
• Language: English
• Published: Cairo, Egypt Hindawi Publishing Corporation 01.01.2019; Hindawi; John Wiley & Sons, Inc; Hindawi Limited; Hindawi-Wiley
• Subjects: Calibration; Computational fluid dynamics; Cores; Crops; Discontinuity; Fault lines; Fluids; Geological processes; Geology; Geophysical data; Geophysics; Geothermal fields; Heat flow; Heat transmission; Horizon; Hypotheses; Lithology; Modelling; Outcrops; Phase transitions; Plate tectonics; Reflectance; Reflection; Seismic activity; Seismic profiles; Seismic reflection method; Seismic response; Seismic surveys; Seismic velocities; Seismic wave velocities; Temperature; Thickness; Three dimensional models; Two dimensional models; Velocity; Velocity distribution; Wave propagation; Wave velocity; Apennines; Italy
• ISSN: 1468-8115; 1468-8123
• DOI: 10.1155/2019/8492453

We tested synthetic seismic reflection modelling along a seismic line (CROP-18A) in the geothermal field at Larderello (Italy). This seismic line is characterized by a discontinuous but locally very bright seismic marker, named K-horizon, which has been associated with various geological processes, including the presence of fluids at supercritical conditions. Geological and geophysical data were integrated in order to develop a 3D subsurface model of a portion of the Larderello field, where extremely high heat flow values have been recorded. In the study area, the K-horizon is particularly shallow and supercritical deep conditions were accessed at depth. The 2D model of the main geological units up to the K-horizon was extracted from the 3D model along the CROP-18A and used to generate the synthetic TWT stacked seismic sections which were then compared with the observed stacked CROP-18A seismic section. To build the synthetic sections, generated through the exploding reflector approach, a 2D velocity model was created assigning to each pixel of the model a constant P-wave velocity corresponding to the related geological unit. The geophysical parameters and the geological model reconstructions used in the modelling process derive from a multidisciplinary integration process including geological outcrop analogues, core samples, and geophysical and laboratory information. Two geophysical models were used to test the seismic response of the K-horizon, which is associated with (1) a lithological discontinuity or (2) a physically perturbed layer, represented by a randomized velocity distribution in a thin layer. For the latter geophysical model (i.e., the physically perturbed layer), we have tested three different scenarios changing the shape and the thickness of the modelled layer. Despite the reliable calibration implied by the use of homogeneous units, the seismic modelling clearly shows that the physically perturbed layer provides a better explanation of the reflectivity features associated with the K-horizon.