Compaction control of topography and fault network structure along strike-slip faults in sedimentary basins

Strike-slip faults commonly display structurally complex areas of positive or negative topography. Understanding the development of such areas has important implications for earthquake studies and hydrocarbon exploration. Previous workers identified the key factors controlling the occurrence of both topographic modes and the related structural styles. Kinematic and stress boundary conditions are of first-order relevance. Surface mass transport and material properties affect fault network structure. Experiments demonstrate that dilatancy can generate positive topography even under simple-shear boundary conditions. Here, we use physical models with sand to show that the degree of compaction of the deformed rocks alone can determine the type of topography and related surface fault network structure in simple-shear settings. In our experiments, volume changes of ∼5% are sufficient to generate localized uplift or subsidence. We discuss scalability of model volume changes and fault network structure and show that our model fault zones satisfy geometrical similarity with natural flower structures. Our results imply that compaction may be an important factor in the development of topography and fault network structure along strike-slip faults in sedimentary basins.

Constraint effects on multiple void interaction in pressure sensitive plastic solids

In this work, effects of pressure sensitive yielding and plastic dilatancy on void growth and void interaction mechanisms in fracture specimens displaying high and low constraint levels are investigated. To this end, large deformation finite element simulations are carried out with discrete voids ahead of the notch. It is observed that multiple void interaction mechanism which is favored by high initial porosity is further accelerated by pressure sensitive yielding, but is retarded by loss of constraint. The resistance curves predicted based on a simple void coalescence criterion show enhancement in fracture resistance when constraint level is low and when pressure sensitivity is suppressed.

A dilation-driven vortex flow in sheared granular materials explains a rheometric anomaly

Granular flows occur widely in nature and industry, yet a continuum description that captures their important features is yet not at hand. Recent experiments on granular materials sheared in a cylindrical Couette device revealed a puzzling anomaly, wherein all components of the stress rise nearly exponentially with depth. Here we show, using particle dynamics simulations and imaging experiments, that the stress anomaly arises from a remarkable vortex flow. For the entire range of fill heights explored, we observe a single toroidal vortex that spans the entire Couette cell and whose sense is opposite to the uppermost Taylor vortex in a fluid. We show that the vortex is driven by a combination of shear-induced dilation, a phenomenon that has no analogue in fluids, and gravity flow. Dilatancy is an important feature of granular mechanics, but not adequately incorporated in existing models.

An Approximate Method for Evaluating the Shear Band Thickness in Saturated Soils

formula for the thickness of a shear band formed in saturated soils under a simple shear or a combined stress state has been proposed. It is shown that the shear band thickness is dependent on the pore pressure properties of the material and the dilatancy rate, but is independent of the details of the combined stress state. This is in accordance with some separate experimental observations.

Formation of shear bands in plane sheet

The formation of shear bands in plane sheet is studied, both analytically and experimentally, to enhance the fundamental understanding of this phenomenon and to develop a capability for predicting material failure. The evolution of voids is measured and its interaction with the process of shear banding is examined. The evolving dilatancy in plasticity is shown to have a vital role in analysing the shear-band type of bifurcation, and tremendously reduces the theoretical value of critical stresses. The analyses, referring to both localized and diffuse modes of bifurcation, fairly explain the corresponding observations obtained through testing a dual-phase steer sheet and provide a justification of the constitutive model used.

Plastic stress-strain fields and blunting effects during large deformations

Plastic stress-strain fields of two types of steel specimens loaded to large deformations are studied. Computational results demonstrate that, owing to the fact that the hardening exponent of the material varies as strain enlarges and the blunting of the crack tip, the well known HRR stress field in the plane strain model can only hold for the stage of a small plastic strain. Plastic dilatancy is shown to have substantial effects on strain distributions and blunting. To justify the constitutive equations used for analysis and to check the precision of computations, the load-deflection of a three-point bend beam and the load-elongation of an axisymmetric bar notched by a V-shaped cut were tested and recorded. The computed curves are in good accordance with experimental data.

Slope Stability Analysis Using Anisotropic and Nonlinear Failure Criteria

The shear strength of soils or rocks developed in a landslide usually exhibits anisotropic and nonlinear behavior. The process of sedimentation and subsequent consolidation can cause anisotropy of sedimentary soils or rocks, for instance. Nonlinearity of failure envelope could be attributed to "interlocking" or "dilatancy" of the material, which is generally dependent upon the stress level. An analytical method considering both anisotropy and nonlinearity of the failure envelops of soil and rocks is presented in the paper. The nonlinearfailure envelopes can be determined from routine triaxial tests. A spreadsheet program, which uses the Janbu's Generalized Procedure of Slice and incorporates anisotropic, illustrates the implementation of the approach and nonlinearfailure envelops. In the analysis, an equivalent Mohr-Coulomb linear failure criterion is obtained by drawing a tangent to the nonlinear envelope of an anisotropic soil at an appropriate stress level. An illustrative example is presented to show the feasibility and numerical efficiency of the method.

Thermo-hydro-geomechanical simulations of methane gas production from deep sea methane hydrate formations

A new constitutive model called Methane Hydrate Critical State (MHCS) model was conducted to investigate the geomechanical response of the gas-hydrate-bearing sediments at the Nankai Trough during the wellbore construction process. The strength and dilatancy of gas-hydrate-bearing soil would gradually disappear when the bonds are destroyed because of excessively shearing, which are often observed in dense soils and also in bonded soils such as cemented soil and unsaturated soil. In this study, the MHCS model, which presents such softening features, would be incorporated into a staged-finite-element model in ABAQUS, which mainly considered the loading history of soils and the interaction between cement-casing-formation. This model shows the influence of gas-hydrate-bearing soil to the deformation and stability of a wellbore and the surrounding sediments during wellbore construction. At the same time, the conventional Mohr-Coulomb model was used in the model to show the advantages of MHCS model by comparing the results of the two models.

Initiation of immersed granular avalanches.

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.

裂隙/孔隙型弹性介质中地震波场模拟与叠前逆时偏移

With the development of oil/gas seismic exploration, seismic survey for fracture/porosity type reservoir is becoming more and more important. As for China, since it has over 60% store of low porosity and low permeability oil/gas reservoir, it’s more urgent to validly describe fracture/porosity type oil/gas trap and proposing the related, developed seismic technique. To achieve mapping fracture/porosity region and its development status, it demands profound understanding of seismic wave propagation discipline in complex fractured/pored media. Meanwhile, it has profound scientific significance and applied worth to study forward modeling of fracture/porosity type media and pre-stacked reverse time migration. Especially, pre-stacked reverse-time migration is the lead edge technique in the field of seismology and seismic exploration. In this paper, the author has summarized the meaning, history and the present state of numerical simulation of seismic propagation in fractured/pored media and seismic exploration of fractured/pored reservoirs. Extensive Dilatancy Anisotropy (EDA) model is selected as media object in this work. As to forward modeling, due to local limitation of solving spatial partial derivative when using finite-difference and finite-element method, the author turns to pseudo-spectral method (PSM), which is based on the global characteristic of Fourier transform to simulate three-component elastic wave-field. Artifact boundary effect reduction and simulation algorithm stability are also discussed in the work. The author has completed successfully forward modeling coding of elastic wave-field and numerical simulation of two-dimensional and three-dimensional EDA models with different symmetric axis. Seismic dynamic and kinematical properties of EDA media are analyzed from time slices and seismic records of wave propagation. As to pre-stacked reverse-time migration for elastic wave-field in fractured/pored media, based on the successful experience in forward modeling results with PSM, the author has studied pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field in two dimensional EDA media induced by preferred fracture/pore distribution. At the same time, different image conditions will bring up what kind of migration result is detailed in this paper. The author has worded out software for pre-stacked reverse-time depth-domain migration of elastic wave-field in EDA media. After migration processing of a series of seismic shot gathers, influences to migration from different isotropic and anisotropy models are described in the paper. In summary, following creative research achievements are obtained:  Realizing two-dimensional and three-dimensional elastic wave-field modeling for fractured/pored media and related software has been completed.  Proposed pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field.  Through analysis of the seismic dynamic and kinematical properties of EDA media, the author made a conclusion that collection of multi-component seismic data can provide important data basis for locating and describing the fracture/pore regions and their magnitudes and the preferred directions.  Pre-stacked reverse-time depth-domain migration technique has the ability to reconstruct complex geological object with steep formations and tilt fracture distribution. Neglecting seismic anisotropy induced by the preferred fracture/pore distribution, will lead to the disastrous imaging results.