Polyaxial stress-dependent permeability of a three-dimensional fractured rock layer

A study about the influence of polyaxial (true-triaxial) stresses on the permeability of a three-dimensional (3D) fractured rock layer is presented. The 3D fracture system is constructed by extruding a two-dimensional (2D) outcrop pattern of a limestone bed that exhibits a ladder structure consisting of a “through-going” joint set abutted by later-stage short fractures. Geomechanical behaviour of the 3D fractured rock in response to in-situ stresses is modelled by the finite-discrete element method, which can capture the deformation of matrix blocks, variation of stress fields, reactivation of pre-existing rough fractures and propagation of new cracks. A series of numerical simulations is designed to load the fractured rock using various polyaxial in-situ stresses and the stress-dependent flow properties are further calculated. The fractured layer tends to exhibit stronger flow localisation and higher equivalent permeability as the far-field stress ratio is increased and the stress field is rotated such that fractures are preferentially oriented for shearing. The shear dilation of pre-existing fractures has dominant effects on flow localisation in the system, while the propagation of new fractures has minor impacts. The role of the overburden stress suggests that the conventional 2D analysis that neglects the effect of the out-of-plane stress (perpendicular to the bedding interface) may provide indicative approximations but not fully capture the polyaxial stress-dependent fracture network behaviour. The results of this study have important implications for understanding the heterogeneous flow of geological fluids (e.g. groundwater, petroleum) in subsurface and upscaling permeability for large-scale assessments.     

Silt- and Bioclastic-Rich Flocs and Their Relationship to Sedimentary Structures: Modern Observations from the Petitcodiac River Estuary

Fine-grained sediment deposition is often conceived to happen in still water conditions and to represent slow sedimentation rates. However, it has been increasingly noted that mud sedimentation commonly occurs under high-energy conditions. The macrotidal Petitcodiac River estuary is used as an example to investigate the significance of flocculation as an important process in the rapid removal of large amounts of sediments from suspension under turbulent flow conditions. A range of physical sedimentary structures was observed at the Petitcodiac River estuary intertidal flats and bars, including low-angle and horizontal planar lamination, current and climbing ripples, surficial fluid mud, soft-sediment deformation, microfaults, and mud rip-up clasts. Fluid mud within the study area contains considerable proportions of clay (21–67%) and contributes to the formation of creeping fluid-mud sheets and streams. Detailed examination of the naturally occurring clay flocs shows that they contain up to 77% of the entangled silt- and sand-sized grains. SEM and microscopic imaging of fluid mud reveal a substantial amount of bioclastic material within the flocs and dispersed among the sediments. These observations show that physicochemical and biological processes influence silt, plankton, and clay aggregation. Water samples and observations from the Petitcodiac River estuary confirm that flocs form in the water column and then settle to the tidal-channel floor and flanking intertidal flats. Laboratory experiments, using Petitcodiac sediment, show how the clay flocs have the potential to sweep the water of silt-sized suspended grains. Additionally, diatoms and their associated extracellular polymeric substance (EPS) sheaths might play a role in mineral aggregation and increased sediment cohesiveness.     

Sedimentologic and stratigraphic constraints on emplacement of the Star Kimberlite, east–central Saskatchewan

Diamond-bearing kimberlites in the Fort à la Corne region, east–central Saskatchewan, consist primarily of extra-crater pyroclastic deposits which are interstratified with Lower Cretaceous (Albian and Cenomanian) marine, marginal marine and continental sediments. Approximately 70 individual kimberlite occurrences have been documented. The Star Kimberlite, occurring at the southeastern end of the main Fort à la Corne trend, has been identified as being of economic interest, and is characterized by an excellent drill core database. Integration of multi-disciplinary data-sets has helped to refine and resolve models for emplacement of the Star Kimberlite. Detailed core logging has provided the foundation for sedimentological and volcanological studies and for construction of a regionally consistent stratigraphic and architectural framework for the kimberlite complex. Micropaleontologic and biostratigraphic analysis of selected sedimentary rocks, and U–Pb perovskite geochronology on kimberlite samples have been integrated to define periods of kimberlite emplacement. Radiometric age determination and micropaleontologic evidence support the hypothesis that multiple kimberlite eruptive phases occurred at Star. The oldest kimberlite in the Star body erupted during deposition of the predominantly continental strata of the lower Mannville Group (Cantuar Formation). Kimberlites within the Cantuar Formation include terrestrial airfall deposits as well as fluvially transported kimberlitic sandstone and conglomerate. Successive eruptive events occurred contemporaneous with deposition of the marginal marine upper Mannville Group (Pense Formation). Kimberlites within the Pense Formation consist primarily of terrestrial airfall deposits. Fine- to medium-grained cross-stratified kimberlitic (olivine-dominated) sandstone in this interval reflects reworking of airfall deposits during a regional marine transgression. The location of the source feeder vents of the Cantuar and Pense kimberlite deposits has not been identified. The youngest and volumetrically most significant eruptive events associated with the Star Kimberlite occur within the predominantly marine Lower Colorado Group (Joli Fou and Viking Formations). Kimberlite beds, which occur at several horizons within these units, consist of subaerial and marine fall deposits, the latter commonly exhibiting evidence of wave-reworking. Black shale-encased resedimented kimberlite beds, likely deposited as subaqueous debris flows and turbidites, are particularly common in the Lower Colorado Group. During its multi-eruptive history, the Star Kimberlite body is interpreted to have evolved from a feeder vent and overlying positive-relief tephra ring, into a tephra cone. Initial early Joli Fou volcanism resulted in formation of a feeder vent (200 m diameter) and tephra ring. Subsequent eruptions, dominated by subaerial deposits, partly infilled the crater and constructed a tephra cone. A late Joli Fou eruption formed a small (70 m diameter) feeder pipe slightly offset to the NW of the early Joli Fou feeder vent. Deposits from this event further infilled the crater, and were deposited on top of early Joli Fou kimberlite (proximal to the vent) and sediments of the Joli Fou Formation (distal to the vent). The shape of the tephra cone was modified during multiple marine transgression and regression cycles coeval with deposition of the Lower Colorado Group, resulting in wave-reworked kimberlite sand along the fringes of the cone and kimberlitic event deposits (tempestites, turbidites, debris flows) in more distal settings.     

Discussion of “Geology and diamond distribution of the 140/141 kimberlite, Fort à la Corne, central Saskatchewan, Canada”, by A. Berryman, B.H. Scott-Smith and B.C. Jellicoe (Lithos v. 76, p. 99–114)

A wide variety of geological data and geological observations by numerous geoscientists do not support a two-stage crater excavation and in-fill model, or a champagne glass-shaped geometry for the 169 or 140/141 kimberlite bodies in the Fort à la Corne kimberlite field, Saskatchewan as described by Berryman, A., Scott Smith, B.H., Jellicoe, B., (2004). Rather, these kimberlite bodies are best described as polygenetic kimberlite tephra cones and tuff rings with associated feeder vents of variable geometry as shown by previous workers for the 169 kimberlite, the 140/141 kimberlite and the Star kimberlite. The domal tephra cone geometry is preserved due to burial by conformable Cretaceous marine mudstones and siltstones and is not an artifact of Quaternary glacial processes.     

Prediction of in-situ block size distributions with reference to armourstone for breakwaters

Understanding a quarry in terms of its potential for breakwater construction materials presents a special challenge for the engineering geologist. Unlike blasting in aggregates and mining operations, optimisation of the extraction process has a focus on the potential for production of large blocks for armourstone. These blocks weighing many tonnes are used for cover layers to resist wave action. The quarry-run is used for breakwater core. If the quarry has been developed as a source of materials dedicated to a breakwater construction project, the success of the project depends greatly on the blasting and production of rock sizes that are required and the avoidance of leaving a massive quantity of unused materials behind in the quarry after project completion. Prediction of in-situ block sizes such as from joint spacing data, provides the most critical input for the prediction of the blast pile block size distribution (BBSD), which in turn is a vital early design input if the constructed breakwater is to be economical as well as effective.This paper is part of a series of papers that introduces the coastal engineering motivation for this work on engineering geology, giving reasons why the prediction of the fragmentation curve of the blast products in a dedicated quarry is of such economic importance for breakwater projects. The first step towards blasted block size distribution (BBSD) prediction is the prediction of the in-situ block sized distribution (IBSD), the main subject of this paper. Drawing together research methods from the 1990s and the rock mechanics principles of discontinuity analysis, a practical step by step methodology for IBSD assessment that includes approaches that are not reliant on specialised computer software is presented. Continuing on the practical theme, a new extension of the volumetric joint count approach is suggested for IBSD prediction for the case when sparse borehole data is all that is available. A case study of IBSD assessment and the associated BBSD and blast assessment is presented from a Carboniferous limestone quarry. For clarity, details of blast design and yield curve prediction that are recommended for use in the context of armourstone production, have been presented in a companion paper. The Rosin-Rammler equation is used as an example form for the BBSD prediction of a dedicated quarry and the potential for breakwater project optimisation is illustrated. The final section sets out a method for directly comparing yield curves together with the demand for materials set by the breakwater design. On the same plot, sizes where there is a relative shortfall in production can be identified. The dependence of effective breakwater design on accurate quarry yield prediction and quarry blasting performance is discussed.     

Modelling of fluid–solid interactions using an adaptive mesh fluid model coupled with a combined finite–discrete element model

Fluid–structure interactions are modelled by coupling the finite element fluid/ocean model ‘Fluidity-ICOM’ with a combined finite–discrete element solid model ‘Y3D’. Because separate meshes are used for the fluids and solids, the present method is flexible in terms of discretisation schemes used for each material. Also, it can tackle multiple solids impacting on one another, without having ill-posed problems in the resolution of the fluid’s equations. Importantly, the proposed approach ensures that Newton’s third law is satisfied at the discrete level. This is done by first computing the action–reaction force on a supermesh, i.e. a function superspace of the fluid and solid meshes, and then projecting it to both meshes to use it as a source term in the fluid and solid equations. This paper demonstrates the properties of spatial conservation and accuracy of the method for a sphere immersed in a fluid, with prescribed fluid and solid velocities. While spatial conservation is shown to be independent of the mesh resolutions, accuracy requires fine resolutions in both fluid and solid meshes. It is further highlighted that unstructured meshes adapted to the solid concentration field reduce the numerical errors, in comparison with uniformly structured meshes with the same number of elements. The method is verified on flow past a falling sphere. Its potential for ocean applications is further shown through the simulation of vortex-induced vibrations of two cylinders and the flow past two flexible fibres.

     

Modelling of fluid–solid interactions using an adaptive mesh fluid model coupled with a combined finite–discrete element model

Fluid–structure interactions are modelled by coupling the finite element fluid/ocean model ‘Fluidity-ICOM’ with a combined finite–discrete element solid model ‘Y3D’. Because separate meshes are used for the fluids and solids, the present method is flexible in terms of discretisation schemes used for each material. Also, it can tackle multiple solids impacting on one another, without having ill-posed problems in the resolution of the fluid’s equations. Importantly, the proposed approach ensures that Newton’s third law is satisfied at the discrete level. This is done by first computing the action–reaction force on a supermesh, i.e. a function superspace of the fluid and solid meshes, and then projecting it to both meshes to use it as a source term in the fluid and solid equations. This paper demonstrates the properties of spatial conservation and accuracy of the method for a sphere immersed in a fluid, with prescribed fluid and solid velocities. While spatial conservation is shown to be independent of the mesh resolutions, accuracy requires fine resolutions in both fluid and solid meshes. It is further highlighted that unstructured meshes adapted to the solid concentration field reduce the numerical errors, in comparison with uniformly structured meshes with the same number of elements. The method is verified on flow past a falling sphere. Its potential for ocean applications is further shown through the simulation of vortex-induced vibrations of two cylinders and the flow past two flexible fibres.     

A Blind Comparison Between Results of Four Image Analysis Systems Using a Photo-Library of Piles of Sieved Fragments

Four image analysis systems for measuring rock fragmentation: FragScan, PowerSieve®, Split and WipFrag, have been compared under conditions necessary to provide an objective though limited assessment of their capabilities. The analysis of results is based on a sample of ten photographs taken from a series of photographs of controlled artificial muckpiles. These were created from dumping a blended mixture of sieved samples of limestone aggregate, in order to create a range of near perfect Rosin-Rammler sieve size distributions. Results from the various systems are compared with sieved results using both histogram and cumulative forms, with and without fines corrections in the case of Split and Wipfrag. Statistical indicators are evaluated to examine the match between system prediction values and sieving values. Commentaries on the results by the inventors of each system have been incorporated. All four systems were found to perform both well in some cases and poorly in others. From a detailed examination of the results, some insight into the strengths and weaknesses of the various systems is presented.     

The effects of a large destructive local earthquake on earthquake preparedness as assessed by an earthquake preparedness scale

A multi-act scale for measuring the earthquake preparedness of individuals and small businesses was developed and used to assess the earthquake preparedness and the perceived difficulty of becoming prepared for earthquakes of 291 University of Southern California undergraduate students approximately 3 weeks prior to the 5.9 magnitude Whittier Narrows earthquake on 1 October 1987. These data were then compared with similar information collected over a 2 1/2 month period following the earthquake from randomly selected samples of subjects that had participated in the original survey. In one case, levels of preparedness of a single group were measured at approximately 2 week intervals (the repeated measures study). In a second case, this information was collected at different points in time following the earthquake from different groups of subjects (the between groups study). Results indicated an initial increase in earthquake preparedness which was significant for subjects in the repeated measures group and which approached significance for subjects in the between groups. This increase in preparedness was maintained for subjects in the repeated measure study but gradually declined to pre-earthquake levels for subjects in the between subjects study. Subjects in the repeated measures study also perceived earthquake preparedness as a significantly less difficult task following the earthquake. Level of perceived difficulty continued to decline over the 2 1/2 month study period. Declines in perceived difficulty for the between subjects study were more erratic, and were only approximately 1/3 of that for the repeated measures group at the end of 2 1/2 months.