Search

Abstract The Irpinia Fault, also known as the Monte Marzano Fault System, located in the Southern Apennines (Italy), is one of the most seismically active structures in the Mediterranean. It is the source of the 1980, Ms 6.9, multi-segment rupture earthquake that caused significant damage and nearly 3,000 casualties. Paleoseismological surveys indicate that this structure has generated at least four Mw ~ 7 surface-rupturing earthquakes in the past 2 ka. This paper presents a comprehensive, high-resolution geophysical investigation focused on the southernmost fault segment of the Monte Marzano Fault System, i.e., the Pantano-Ripa Rossa Fault, outcropping within the Pantano di San Gregorio Magno intramontane basin. The project, named TEst Site IRpinia fAult (TESIRA), was supported by the University of Napoli Federico II to study the near-surface structure of this intra-basin fault splay that repeatedly ruptured co-seismically in the past thousands of years. Our imaging approach included 2D and 3D electrical and seismic surveys, gravimetry, 3D FullWaver electrical tomography, drone-borne GPR and magnetic surveys, and CO2 soil flux assessment across the surface rupture. This multidisciplinary investigation improved our understanding of the basin shallow structure, providing an image of a rather complex subsurface fault and basin geometry. Seismic data suggest that fault activity at the Pantano segment of MMFS is characterized by a near-surface cumulative displacement greater than previous estimations, calling into question earlier assumptions about the timing of its activation. Despite some challenges with our drone-mounted survey equipment, the integrated dataset provides a comprehensive and reliable image of the subsurface structure. This work demonstrates the utility of developing an integrated approach at high-resolution geophysical imaging and interpretation of fault zones with weak morphological expressions.

The textural and geometrical properties of the pore networks (i.e., such as pore size distribution, pore shape, connectivity, and tortuosity) provides a primary control on the fluid storage and migration of geofluids within porous carbonate reservoirs. These properties are highly variable because of primary depositional conditions, diagenetic processes and deformation. This issue represents an important challenge for the characterization and exploitation plan in this type of reservoirs. In this study, the complementary properties of neutrons and X-ray experiments are carried out to better understand the effects of pore network properties on the hydraulic behavior of porous carbonates. Neutrons have unique properties and are particularly suitable for this study due to the sensitivity of neutrons to hydrogen-based fluids. The used methodology combines dynamic neutron radiography (NR), integrated X-ray and neutron tomography (XCT, NCT), and computational fluid dynamics simulations (lattice-Boltzmann method) of porous carbonate reservoir analogs from central and southern Italy. Dynamic 2D NR images provide information regarding the fluid transport and the wetting front dynamics related to the effect of heterogeneities (e.g., fractures and deformation bands) at the microscale. The combination of NCT (dry and wet samples) and XCT (dry), generates more information regarding the effective pore space contribution to fluid flow. The fluid flow simulations generate information about the connected pore network and the permeability evaluated rock sample at saturated condition.

Geofluid reservoirs located in heterolithic successions (e.g., turbidites) can be affected by vertical and lateral compartmentalization due to interbedded fine-grained facies (i.e., shale, siltstones) and the presence of faults, respectively. A fault can behave as a conduit or barrier to fluid flow depending on its architecture and the individual hydraulic behavior of its components (i.e., fault core, damage zone). The fault core, normally composed by fault rock or smeared clay material, commonly acts as a flow inhibitor across the fault. Fault-related fractures (macro- and microscopic) in the damage zone generally increase the permeability parallel to the fault, except when they are cemented or filled with gouge material. Although macrofractures (which define the fracture porosity) dominate fluid flow, the matrix porosity (including microfractures) begins to have a more important role in fluid flow as the aperture of macrofractures is occluded, particularly at greater depth. This study investigates the variation in matrix permeability in fault zones hosted in heterolithic successions due to fault architecture and stratigraphy of host rock (i.e., sand-rich turbidites). Two key areas of well-exposed, faulted Miocene turbidites located in central and southern Italy were selected. For this study, six separate fault zones of varying offset were chosen. Each impacts heterolithic successions that formed under similar tectonic conditions and burial depths. Across the selected fault zones, an extensive petrophysical analysis was done in the field and laboratory, through air permeameter measurements, thin section, and synchrotron analysis in both host rock, damage zone, and fault core. Results suggest that the amount and distribution of clay layers in a heterolithic sequence affects fluid flow across the fault, regardless of fault offset.

Fluid flow through a single fracture is traditionally described by the cubic law, which is derived from the Navier-Stokes equation for the flow of an incompressible fluid between two smooth-parallel plates. Thus, the permeability of a single fracture depends only on the so-called hydraulic aperture which differs from the mechanical aperture (separation between the two fracture wall surfaces). This difference is mainly related to the roughness of the fracture walls, which has been evaluated in previous works by including a friction factor in the permeability equation or directly deriving the hydraulic aperture. However, these methodologies may lack adequate precision to provide valid results. This work presents a complete protocol for fracture surface mapping, roughness evaluation, fracture modeling, fluid flow simulation, and permeability estimation of individual fracture (open or sheared joint/pressure solution seam). The methodology includes laboratory-based high-resolution structure from motion (SfM) photogrammetry of fracture surfaces, power spectral density (PSD) surface evaluation, synthetic fracture modeling, and fluid flow simulation using the Lattice-Boltzmann method. This work evaluates the respective controls on permeability exerted by the fracture displacement (perpendicular and parallel to the fracture walls), surface roughness, and surface pair mismatch. The results may contribute to defining a more accurate equation of hydraulic aperture and permeability of single fractures, which represents a pillar for the modeling and upscaling of the hydraulic properties of a geofluid reservoir.