Home About us Contact
|
Home About us Contact | ||
|
|
||
Elastic Medium (elastic + medium)
Selected AbstractsThree-dimensional models of elastostatic deformation in heterogeneous media, with applications to the Eastern California Shear ZoneGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2009Sylvain Barbot SUMMARY We present a semi-analytic iterative procedure for evaluating the 3-D deformation due to faults in an arbitrarily heterogeneous elastic half-space. Spatially variable elastic properties are modelled with equivalent body forces and equivalent surface traction in a ,homogenized' elastic medium. The displacement field is obtained in the Fourier domain using a semi-analytic Green function. We apply this model to investigate the response of 3-D compliant zones (CZ) around major crustal faults to coseismic stressing by nearby earthquakes. We constrain the two elastic moduli, as well as the geometry of the fault zones by comparing the model predictions to Synthetic Aperture Radar inferferometric (InSAR) data. Our results confirm that the CZ models for the Rodman, Calico and Pinto Mountain faults in the Eastern California Shear Zone (ECSZ) can explain the coseismic InSAR data from both the Landers and the Hector Mine earthquakes. For the Pinto Mountain fault zone, InSAR data suggest a 50 per cent reduction in effective shear modulus and no significant change in Poisson's ratio compared to the ambient crust. The large wavelength of coseismic line-of-sight displacements around the Pinto Mountain fault requires a fairly wide (,1.9 km) CZ extending to a depth of at least 9 km. Best fit for the Calico CZ, north of Galway Dry Lake, is obtained for a 4 km deep structure, with a 60 per cent reduction in shear modulus, with no change in Poisson's ratio. We find that the required effective rigidity of the Calico fault zone south of Galway Dry Lake is not as low as that of the northern segment, suggesting along-strike variations of effective elastic moduli within the same fault zone. The ECSZ InSAR data is best explained by CZ models with reduction in both shear and bulk moduli. These observations suggest pervasive and widespread damage around active crustal faults. [source] The role of friction and secondary flaws on deflection and re-initiation of hydraulic fractures at orthogonal pre-existing fracturesGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2006Xi Zhang SUMMARY In this study, we explore the nature of plane-strain hydraulic fracture growth in the presence of pre-existing fractures such as joints without or with secondary flaws. The 2-D plane-strain fracture studied can be taken as a cross-section through the short dimensions of an elongated 3-D fracture or as an approximate representation of the leading edge of a 3-D fracture where the edge curvature is negligible. The fluid-driven fracture intersects a pre-existing fracture to which it is initially perpendicular and is assumed not to immediately cross, but is rather deflected into the pre-existing fracture. The intersection results in branching of the fracture and associated fluid flow into the pre-existing fracture. Further growth results in opening and frictional sliding along the pre-existing fracture. Fracture propagation in an impermeable homogeneous elastic medium and fluid invasion into a pre-existing fracture are both driven by an incompressible, Newtonian fluid injected at a constant rate. The frictional stress on the surfaces of pre-existing fractures is assumed to obey the Coulomb law. The governing equations for quasi-static fluid-driven fracture growth are given and a scaling is introduced to help identify important parameters. The displacement discontinuity method and the finite difference method are employed to deal with this coupling mechanism of rock fracture and fluid flow. In order to account for fluid lag, a method for separately tracking the crack tip and the fluid front is included in the numerical model. Numerical results are obtained for internal pressure, frictional contact stresses, opening and shear displacements, and fluid lag size, as well as for fracture re-initiation from secondary flaws. After fracture intersection, the hydraulic fracture growth mode changes from tensile to shearing. This contributes to increased injection pressure and to a reduction in fracture width. In the presence of pre-existing fractures, the fluid-driven cracks can be arrested or retarded in growth rate as a result of diversion of fluid flow into and frictional sliding along the pre-existing fractures. Frictional behaviour significantly affects the ability of the fluid to enter or penetrate the pre-existing fracture only for those situations where the fluid front is within a certain distance from the intersecting point. Importantly, fluid penetration requires higher injection pressure for frictionally weak pre-existing fractures. Fracture re-initiation from secondary flaws can reduce the injection pressure, but re-initiation is suppressed by large sliding on pre-existing fractures that are frictionally weak. [source] On strike-slip faulting in layered mediaGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2002Maurizio Bonafede Summary We study the effects of structural inhomogeneities on the stress and displacement fields induced by strike-slip faults in layered media. An elastic medium is considered, made up of an upper layer bounded by a free surface and welded to a lower half-space characterized by different elastic parameters. Shear cracks with assigned stress drop are employed as mathematical models of strike-slip faults, which are assumed to be vertical and planar. If the crack is entirely embedded within the lower medium (case A), a Cauchy-kernel integral equation is obtained, which is solved by employing an expansion of the dislocation density in Chebyshev polynomials. If the crack is within the lower medium but it terminates at the interface (case B), a generalized Cauchy singularity appears in the integral kernel. This singularity affects the singular behaviour of the dislocation density at the crack tip touching the interface. Finally, the case of a crack crossing the interface is considered (case C). The crack is split into two interacting sections, each placed in a homogeneous medium and both open at the interface. Two coupled generalized Cauchy equations are obtained and solved for the dislocation density distribution of each crack section. An asymptotic study near the intersection between the crack and the interface shows that the dislocation densities for each crack section are bounded at the interface, where a jump discontinuity is present. As a corollary, the stress drop must be discontinuous at the interface, with a jump proportional to the rigidity contrast between the adjoining media. This finding is shown to have important implications for the development of geometrical complexities within transform fault zones: planar strike-slip faults cutting across layer discontinuities with arbitrary stress drop values are shown to be admissible only if the interface between different layers becomes unwelded during the earthquake at the crack/interface junction. Planar strike-slip faulting may take place only in mature transform zones, where a repetitive earthquake cycle has already developed, if the rheology is perfectly elastic. Otherwise, the fault cannot be planar: we infer that strike-slip faulting at depth is plausibly accompanied by en-echelon surface breaks in a shallow sedimentary layer (where the stress drop is lower than prescribed by the discontinuity condition), while ductile deformation (or steady sliding) at depth may be accommodated by multiple fault branching or by antithetic faulting in the upper brittle layer (endowed with lower rigidity but higher stress). [source] Vertically fractured transversely isotropic media: dimensionality and deconstructionGEOPHYSICAL PROSPECTING, Issue 2 2009Michael A. Schoenberg ABSTRACT A vertically fractured transversely isotropic (VFTI) elastic medium is one in which any number of sets of vertical aligned fractures (each set has its normal lying in the horizontal x1, x2 -plane) pervade the medium and the sets of aligned fractures are the only features of the medium disturbing the axi-symmetry about the x3 -axis implying that in the absence of fractures, the background medium is transversely isotropic (TI). Under the assumptions of long wavelength equivalent medium theory, the compliance matrix of a fractured medium is the sum of the background medium's compliance matrix and a fracture compliance matrix. For sets of parallel rotationally symmetric fractures (on average), the fracture compliance matrix is dependent on 3 parameters , its normal and tangential compliance and its strike direction. When one fracture set is present, the medium is orthorhombic and the analysis is straightforward. When two (non-orthogonal) or more sets are present, the overall medium is in general elastically monoclinic; its compliance tensor components are subject to two equalities yielding an 11 parameter monoclinic medium. Constructing a monoclinic VFTI medium with n embedded vertical fracture sets, requires 5 TI parameters plus 3×n fracture set parameters. A deconstruction of such an 11 parameter monoclinic medium involves using its compliance tensor to find a background transversely isotropic medium and several sets of vertical fractures which, in the long wavelength limit, will behave exactly as the original 11 parameter monoclinic medium. A minimal deconstruction, would be to determine, from the 11 independent components, the transversely isotropic background (5 parameters) and two fracture sets (specified by 2 × 3 = 6 parameters). Two of the background TI medium's compliance matrix components are known immediately by inspection, leaving nine monoclinic components to be used in the minimal deconstruction of the VFTI medium. The use of the properties of a TI medium, which are linear relations on its compliance components, allows the deconstruction to be reduced to solving a pair of non-linear equations on the orientations of two fracture sets. A single root yielding a physically meaningful minimum deconstruction yields a unique minimal representation of the monoclinic medium as a VFTI medium. When no such root exists, deconstruction requires an additional fracture set and uniqueness is lost. The boundary between those monoclinic media that have a unique minimal representation and those that do not is yet to be determined. [source] The effect of inertial coupling on seismic reflection amplitudesGEOPHYSICAL PROSPECTING, Issue 5 2008Ashish Arora ABSTRACT A problem of reflection and transmission of elastic waves at a plane interface between a uniform elastic solid half-space and a porous elastic half-space containing two immiscible fluids is investigated. The theory developed by Lo, Sposito and Majer for porous media containing two immiscible fluids is employed to find out the reflection and transmission coefficients. The incident wave is assumed to propagate through the uniform elastic half-space and two cases are considered. In the first case, a beam of plane longitudinal wave is assumed to be incident and in the second case, a beam of transverse wave is assumed to be incident at the interface. By taking granite as impervious elastic medium and columbia fine sandy loam containing air-water mixture as porous medium, reflection and transmission coefficients are obtained. By neglecting the inertial coupling coefficients, these coefficients are reduced to those obtained by Tomar and Arora using the theory of Tuncay and Corapcioglu. It is found that the inertial coupling parameters significantly affect the phase speeds and the amplitude ratios of the transmitted waves. [source] Uniform asymptotic Green's functions for efficient modeling of cracks in elastic layers with relative shear deformation controlled by linear springsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 3 2009Anthony P. Peirce Abstract We present a uniform asymptotic solution (UAS) for a displacement discontinuity (DD) that lies within the middle layer of a three-layer elastic medium in which relative shear deformation between parallel interfaces is controlled by linear springs. The DD is assumed to be normal to the two interfaces between the elastic media. Using the Fourier transform method we construct a leading term in the asymptotic expansion for the spectral coefficient functions for a DD in a three-layer-spring medium. Although a closed-form solution will require a solution in terms of an infinite series, we demonstrate how this UAS can be used to construct highly efficient and accurate solutions even in the case in which the DD actually touches the interface. We compare the results using the Green's function UAS solution for a crack crossing a soft interface with results obtained using a multi-layer boundary element method. We also present results from an implementation of the UAS Green's function approach in a pseudo-3D hydraulic fracturing simulator to analyze the effect of interface shear deformation on the fracture propagation process. These results are compared with field measurements. Copyright © 2008 John Wiley & Sons, Ltd. [source] Analytical solution of the harmonic waves diffraction by a cylindrical lined cavity in poroelastic saturated mediumINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 5 2007Y. S. Karinski Abstract This paper presents a model for the analysis of plane waves diffraction at a cavity in an infinite homogeneous poroelastic saturated medium, lined by a lining composed of four equal segments. An elastic boundary layer is placed between the cavity lining and the infinite porous medium. The boundary layer is simulated by ,elastic boundary conditions' in which the bulk matrix stress is proportional to the relative displacement between the lining and the surrounding medium matrix boundary. In addition, fluid impermeability through the intermediate layer is assumed. For the frequencies, that differ from the pseudoresonanse frequencies, the problem was reduced to the problem of an ideal elastic medium. A closed-form analytical solution of the problem was obtained using Fourier,Bessel series, the convergence of which was proven. It was shown that the number of series terms required to obtain a desired level of accuracy can be determined in advance. The influence of the medium porosity on the medium dynamic stress concentration was studied. Copyright © 2006 John Wiley & Sons, Ltd. [source] Self-similar solution of a plane-strain fracture driven by a power-law fluidINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 6 2002J. I. Adachi Abstract This paper analyses the problem of a hydraulically driven fracture, propagating in an impermeable, linear elastic medium. The fracture is driven by injection of an incompressible, viscous fluid with power-law rheology and behaviour index n,0. The opening of the fracture and the internal fluid pressure are related through the elastic singular integral equation, and the flow of fluid inside the crack is modelled using the lubrication theory. Under the additional assumptions of negligible toughness and no lag between the fluid front and the crack tip, the problem is reduced to self-similar form. A solution that describes the crack length evolution, the fracture opening, the net fluid pressure and the fluid flow rate inside the crack is presented. This self-similar solution is obtained by expanding the fracture opening in a series of Gegenbauer polynomials, with the series coefficients calculated using a numerical minimization procedure. The influence of the fluid index n in the crack propagation is also analysed. Copyright © 2002 John Wiley & Sons, Ltd. [source] State-space time integration with energy control and fourth-order accuracy for linear dynamic systemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2006Steen Krenk Abstract A fourth-order accurate time integration algorithm with exact energy conservation for linear structural dynamics is presented. It is derived by integrating the phase-space representation and evaluating the resulting displacement and velocity integrals via integration by parts, substituting the time derivatives from the original differential equations. The resulting algorithm has an exact energy equation, in which the change of energy is equal to the work of the external forces minus a quadratic form of the damping matrix. This implies unconditional stability of the algorithm, and the relative phase error is of fourth-order. An optional high-frequency algorithmic damping is constructed by optimal combination of three different damping matrices, each proportional to either the mass or the stiffness matrix. This leads to a modified form of the undamped algorithm with scalar weights on some of the matrices introducing damping of fourth-order in the frequency. Thus, the low-frequency response is virtually undamped, and the algorithm remains third-order accurate even when algorithmic damping is included. The accuracy of the algorithm is illustrated by an application to pulse propagation in an elastic medium, where the algorithmic damping is used to reduce dispersion due to the spatial discretization, leading to a smooth solution with a clearly defined wave front. Copyright © 2005 John Wiley & Sons, Ltd. [source] A dual mesh multigrid preconditioner for the efficient solution of hydraulically driven fracture problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 13 2005A. P. Peirce Abstract We present a novel multigrid (MG) procedure for the efficient solution of the large non-symmetric system of algebraic equations used to model the evolution of a hydraulically driven fracture in a multi-layered elastic medium. The governing equations involve a highly non-linear coupled system of integro-partial differential equations along with the fracture front free boundary problem. The conditioning of the algebraic equations typically degrades as O(N3). A number of characteristics of this problem present significant new challenges for designing an effective MG strategy. Large changes in the coefficients of the PDE are dealt with by taking the appropriate harmonic averages of the discrete coefficients. Coarse level Green's functions for multiple elastic layers are constructed using a single dual mesh and superposition. Coarse grids that are sub-sets of the finest grid are used to treat mixed variable problems associated with ,pinch points.' Localized approximations to the Jacobian at each MG level are used to devise efficient Gauss,Seidel smoothers and preferential line iterations are used to eliminate grid anisotropy caused by large aspect ratio elements. The performance of the MG preconditioner is demonstrated in a number of numerical experiments. Copyright © 2005 John Wiley & Sons, Ltd. [source] The elastic echo problemMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 2 2003Mongi Mabrouk Abstract We consider a homogeneous isotropic unbounded linear elastic medium ,,,3, having a free boundary ,. A forcing f(t,x) creates an incident displacement field u0(t,x). This primary field is scattered by , giving rise to a secondary field or echo, for which we determine the asymptotic behaviour in time. These results are obtained via the use of an tension of the time-dependent scattering theory of C. Wilcox. Copyright © 2003 John Wiley & Sons, Ltd. [source] Existence of solution in elastic wave scattering by unbounded rough surfacesMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 6 2002T. Arens We consider the two-dimensional problem of the scattering of a time-harmonic wave, propagating in an homogeneous, isotropic elastic medium, by a rough surface on which the displacement is assumed to vanish. This surface is assumed to be given as the graph of a function ,,C1,1(,). Following up on earlier work establishing uniqueness of solution to this problem, existence of solution is studied via the boundary integral equation method. This requires a novel approach to the study of solvability of integral equations on the real line. The paper establishes the existence of a unique solution to the boundary integral equation formulation in the space of bounded and continuous functions as well as in all Lp spaces, p,[1, ,] and hence existence of solution to the elastic wave scattering problem. Copyright © 2002 John Wiley & Sons, Ltd. [source] The effect of a penalty term involving higher order derivatives on the distribution of phases in an elastic medium with a two-well elastic potentialMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 4 2002M. Bildhauer Abstract We consider the problem of minimizing 0 The effect of a surface energy term on the distribution of phases in an elastic medium with a two-well elastic potentialMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 2 2002Michael Bildhauer Abstract We consider the problem of minimizing among functions u:,d,,,,d, u,,,=0, and measurable subsets E of ,. Here fh+, f, denote quadratic potentials defined on ,¯×{symmetric d×d matrices}, h is the minimum energy of fh+ and ,(u) is the symmetric gradient of the displacement field u. An equilibrium state û, Ê of J(u,E) is called one-phase if E=, or E=,, two-phase otherwise. For two-phase states, ,,,E,,, measures the effect of the separating surface, and we investigate the way in which the distribution of phases is affected by the choice of the parameters h,,, ,>0. Additional results concern the smoothness of two-phase equilibrium states and the behaviour of inf J(u,E) in the limit ,,0. Moreover, we discuss the case of additional volume force potentials, and extend the previous results to non-zero boundary values. Copyright © 2002 John Wiley & Sons, Ltd. [source]
| | ||||||||||||