Search USGSSearch

USGS · 70012497

On the mechanical interaction between a fluid-filled fracture and the earth's surface

Abstract

The mechanical interaction between a fluid-filled fracture (e.g., hydraulic fracture joint, or igneous dike) and the earth's surface is analyzed using a two-dimensional elastic solution for a slit of arbitrary inclination buried beneath a horizontal free surface and subjected to an arbitrary pressure distribution. The solution is obtained by iteratively superimposing two fundamental sets of analytical solutions. For uniform internal pressure the slit behaves essentially as if it were in an infinite region if the depth-to-center is three times greater than the half-length. For shallower slits interaction with the free surface is pronounced: stresses and displacements near the slit differ by more than 10% from values for the deeply buried slit. The following changes are noted as the depth-to-center decreases: 1. (1) the mode I stress intensity factor increases for both ends of the slit, but more rapidly at the upper end; 2. (2) the mode II stress-intensity factor is significantly different from zero (except for vertical slits) suggesting propagation out of the original plane of the slit; 3. (3) displacements of the slit wall are asymmetric such that the slit gaps open more widely near the upper end. Similar changes are noted if fluid density creates a linear pressure gradient that is smaller than the lithostatic gradient. Under such conditions natural fractures should propagate preferentially upward toward the earth's surface requiring less pressure as they grow in length. If deformation near the surface is of interest, the model should account explicitly for the free surface. Stresses and displacements at the free surface are not approximated very well by values calculated along a line in an infinite region, even when the slit is far from the line. As depth-to-center of a shallow pressurized slit decreases, the following changes are noted: 1. (1) displacements of the free surface increase to the same order of magnitude as the displacements of the slit walls, 2. (2) tensile stresses of magnitude greater than the pressure in the slit are concentrated along the free surface. The relative surface displacements over a shallow vertical slit are downward over the slit and upward to both sides of this area. The tensile stress acting parallel to the free surface over a shallow vertical slit is concentrated in two maxima adjacent to a point of very low stress immediately over the slit. The solution is used to estimate the length-to-depth ratio at which igneous sills have gained sufficient leverage on overlying strata to bend these strata upward and form a laccolith. The pronounced mode II stress intensity associated with shallow horizontal slits explains the tendency for some sills to climb to higher stratigraphie horizons as they grow in length. The bimodal tensile stress concentration over shallow vertical slits correlates qualitatively with the distribution of cracks and normal faults which flank fissure eruptions on volcanoes. The solution may be used to analyze surface displacements and tilts over massive hydraulic fractures in oil fields and to understand the behavior of hydraulic fractures in granite quarries.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

David D. Pollard, Gary Holzhausen. 1979. On the mechanical interaction between a fluid-filled fracture and the earth's surface. https://doi.org/10.1016/0040-1951(79)90353-6

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics

Fairweather transform boundary Oligocene to present orogenesis: Fairweather Range vertical extrusion and rotation of the Yakutat microplate at ca. 3 Ma

Oblique-slip along transform fault boundaries is often partitioned between a strike-slip system and thrust faults that accommodate contraction. However, topography along the Yakutat-North American transform (Fairweather fault), is asymmetric with low-terrain above active thrusts on the western, Yakutat side of the transform and high topography on the continental side with peaks >4500 m (Mount Fairweather: 4671 m) to the west of the Border Ranges fault, limited recorded earthquakes >M4, and no apparent reverse faults to generate the highest terrain. In this study we compile, for the first time, published U-Pb zircon, 40 Ar/ 39 Ar and K-Ar (hornblende, muscovite, and biotite) and U-Th/He and fission-track (zircon and apatite) bedrock ages (109) from 75 samples to investigate the exhumation history of the Fairweather Range region, complemented by a published detrital sample (ZFT and AFT) and 13 new 40 Ar/ 39 Ar (hornblende, biotite, and K-feldspar) ages on 9 bedrock samples from both sides of the Fairweather fault. Additionally, we examined published seismicity and geodetic data of the Fairweather region and assessed if plate paleo-vectors correlate with the cooling history of the Fairweather Range. Cooling age, seismic, and block-motion patterns indicate the Fairweather Range has been vertically extruded between the Fairweather and the Border Ranges faults as a coherent block since ca. 25 Ma. The pre-6 Ma Pacific plate motion (N30°W) aligns with the N33°W strike of the Fairweather Fault whereas a hypothetical pre-6 Ma Yakutat microplate paleo-vector of (N39°W) does not: indicating a post-6 Ma timing for Yakutat microplate counter-clockwise rotation (9°). We infer that rotation and impingement of the Yakutat microplate along the Fairweather fault at ca. 3 Ma led to the development of the Fairweather restraining bend and increased cooling rates. The resultant thickened Fairweather welt and the ∼30 km thick southeast end of the Yakutat microplate compounded double-indenter tectonics into Alaska's southeast convergent corner

Tectonophysics

A scaling relationship for the width of secondary deformation around strike-slip faults

Simple mechanical arguments suggest that slip along interlocked, rough faults, damages surrounding rocks. The same arguments require that the scale of secondary damage is proportional to the size of geometric irregularities along the main fault. This relationship could apply at all scales, but has, so far, been difficult to observe at the 10s to 100 s of km scales of large, natural faults, often because large-scale deformation is distributed across wide, complex plate-boundary fault systems, like the San Andreas Fault. The geometry and geology of another large-scale plate-boundary strike slip fault—the Queen Charlotte Fault (QCF)—is, in contrast, especially simple. Here, we show that observations of secondary deformation are well-aligned with predictions of stress variations caused by geometric irregularities along the QCF, suggesting a geometric relationship between primary fault geometry and secondary deformation. The analytic stress solution reveals that the highest stresses and highest likelihood of failure are confined to a zone of influence (ZOI) with a width quantified by ZOI = λ / 2 π "> ZOI=λ/2π , where λ is the wavelength of geometric variations along the main fault. This simple model is consistent with ∼100-km-scale observations along the QCF and can theoretically be used to predict the width of secondary deformation at all scales.

Queen Charlotte Fault