Hydraulic fracture treatments induce microseismicity (i.e., discrete events plus noise) which can be recorded only by seismometers located in or near the treatment well bore. Seismograms recorded in the treatment well bore are composed of interacting phases which to date defy discrete identification and preclude standard inversion techniques to determine source characteristics, propagation path characteristics, or both. To understand the recorded data, synthetic seismograms have been generated using a finite difference approximation for a two-dimensional hydraulic fracture model which includes a treatment-induced, low-velocity zone (LVZ). The LVZ is produced by pore dilatancy in response to formation fluid pressurization caused by pressure diffusion from the treatment fluid and fracture opening. Borehole, seismometer coupling, fluid transport, and dynamic fracture effects are not considered. The models examined showed that a fixed-aperture fracture plus surrounding LVZ acted as a leaky waveguide confining the seismic energy and altering propagation characteristics which resulted in unconventional seismograms. The seismograms were composed of multiple transient phases caused by phase interactions and conversions at the waveguide boundaries and stationary phases corresponding to the eigenfrequency structure of the fracture-LVZ waveguide.
BACKGROUNDHydraulic fracturing for enhancing the recovery of hydrocarbons has become a common reservoir engineering practice since it was first introduced at the Hugoton gas field in western Kansas in 1947 [Howard and Fast, 1970]. Over the intervening years, treatment procedures have evolved in an attempt to create fractures having specific or determinable characteristics (i.e., orientation, length, height, and aperture). For example, treatment parameters (e.g., injection rate, injection pressure, proppant concentration, and fracture fluid biochemistry and rheology) are varied in ways thought to maximize fracture fluid efficiency by decreasing fluid loss to the formation, to prevent fracture growth into formations other than the desired formation, and to produce maximum propped fracture length and permeability. During the actual treatment, designed procedures frequently are adjusted to compensate for the effects of unexpected formation heterogeneities, spatially varying stress and equipment malfunctions. Because of the complexity of the treatment and the setting upon which the treatment is applied, a need exists for determining the final effect of the treatment.A number of authors have published models which define the final fracture geometry by using the treatment parameters [e.g., Advanti et al., 1985; Ahmed, 1984; Geertsma and De-Klerk, 1969; Nolte, 1983; Perkins and Kern, 1961; Settari and Cleary, 1984]. These models are strongly dependent upon assumed formation characteristics and their interactions with actual treatment parameters. As a result, techniques for delineating the final fracture or fractures which are independent of formation and treatment assumptions or procedures...