A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

Reef exploration in the East Sengkang Basin,Sulawesi, Indonesia

Reconnaissance seismic shot in 1971/72 showed a number of well defined seismic anomalies within the East Sengkang Basin which were interpreted as buried reefs. Subsequent fieldwork revealed that Upper Miocene reefs outcropped along the southern margin of the basin. A drilling programme in 1975 and 1976 proved the presence of shallow, gas-bearing, Upper Miocene reefs in the northern part of the basin. Seismic acquisition and drilling during 1981 confirmed the economic significance of these discoveries, with four separate accumulations containing about 750 × 10⁹ cubic feet of dry gas in place at an average depth of 700 m. Kampung Baru is the largest field and contains over half the total, both reservoir quality and gas deliverability are excellent. Deposition in the East Sengkang Basin probably started during the Early Miocene. A sequence of Lower Miocene mudstones and limestones unconformably overlies acoustic basement which consists of Eocene volcanics. During the tectonically active Middle Miocene, deposition was interrupted by two periods of deformation and erosion. Carbonate deposition became established in the Late Miocene with widespread development of platform limestones throughout the East Sengkang Basin. Thick pinnacle reef complexes developed in the areas where reef growth could keep pace with the relative rise in sea level. Most reef growth ceased at the end of the Miocene and subsequent renewed clastic sedimentation covered the irregular limestone surface. Late Pliocene regression culminated in the Holocene with erosion. The Walanae fault zone, part of a major regional sinistral strike-slip system, separates the East and West Sengkang Basins. Both normal and reverse faulting are inferred from seismic data and post Late Pliocene reverse faulting is seen in outcrop.     

Microstructures of co-axially folded vein segments and crenulation cleavage： Evidence for dissolution phenomena in the Chamba Thrust Sheet, western Himalayas

The metasediments in the Chamba region experienced three phases of deformation: DF1, DF2 and DF3.Folded quartz veins are co-folded with the F2 crenulation folds. Their geometric and tectonic significance is studied by microstructures and shortening adjacent to the discrete crenulation cleavage, S2. Microstructures of folded vein segments, their geometric changes and truncation to cleavage (S2) are mainly due to pressure-solution phenomena and the estimated volume loss from reconstructed vein segments range from 16 to 25.5%,which is closely related to volume decrease (26%) calculated from the polydeformed slates of North Wales areas.     

Comparison of marine gas hydrates in sediments of an active and passive continental margin

Two sites of the Deep Sea Drilling Project in contrasting geologic settings provide a basis for comparison of the geochemical conditions associated with marine gas hydrates in continental margin sediments. Site 533 is located at 3191 m water depth on a spit-like extension of the continental rise on a passive margin in the Atlantic Ocean. Site 568, at 2031 m water depth, is in upper slope sediment of an active accretionary margin in the Pacific Ocean. Both sites are characterized by high rates of sedimentation, and the organic carbon contents of these sediments generally exceed 0.5%. Anomalous seismic reflections that transgress sedimentary structures and parallel the seafloor, suggested the presence of gas hydrates at both sites, and, during coring, small samples of gas hydrate were recovered at subbottom depths of 238m (Site 533) and 404 m (Site 568). The principal gaseous components of the gas hydrates wer methane, ethane, and CO₂. Residual methane in sediments at both sites usually exceeded 10 mll^?1 of wet sediment. Carbon isotopic compositions of methane, CO₂, and ΣCO₂ followed parallel trends with depth, suggesting that methane formed mainly as a result of biological reduction of oxidized carbon. Salinity of pore waters decreased with depth, a likely result of gas hydrate formation. These geochemical characteristics define some of the conditions associated with the occurrence of gas hydrates formed by in situ processes in continental margin sediments.     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

Plane of elastic non-resisting tension material under foundation structures

In the frame of 2D-static problems one approaches the problem of elastic-NRT (not-resisting tension) semi-plane loaded on its limit line. This problem is intended to model the stress situation induced in the soil by a foundation structure. The solution, in terms of activated stress field, is searched for in the class of stress fields satisfying equilibrium and admissibility conditions, by applying an energy approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Fluid production and isograd reactions at contacts of carbonate-rich and carbonate-poor layers during progressive metamorphism

ABSTRACT With increasing temperature during prograde metamorphism reactions will occur first at the lithological contacts of mixed pelite and calcsilicate terranes. At these interfaces, a fluid of lower chemical potential of H₂O and CO₂ than that required to produce a fluid in either layer can be produced whether reaction is caused by fluid infiltration or is initially fluid absent. If the interface region does not allow fluid transport then as temperature increases, a fluid pressure greater than lithostatic can develop. At some degree of over-pressure relative to rock pressure, the fluid hydraulically fractures the rock and a gradient in fluid composition away from the contact can be produced. These phenomena occur at the compositional interfaces whenever univariant reactions in the differing layers cross on a temperature vs. mole fraction of CO₂ diagram with slopes of opposite sign. The first occurrence of these reaction products at lithological contacts delineates an isograd that defines temperature as well as the mole fraction of CO₂ at constant pressure in systems open to fluid transport. These isograds can be contrasted with fluid-producing isograds in closed systems. As an illustration of possible effects, the reactions quartz + clinozoisite + muscovite = anorthite + K-feldspar + H₂O and phlogopite + quartz + calcite = tremolite + K-feldspar + H₂O + CO₂ at 4 kbar are analysed and equations for fluid production and transport are developed.