Penny-shaped fractures revisited

All methods of seismic characterization of fractured reservoirs are based on effective media theories that relate geometrical and material properties of fractures and surrounding rock to the effective stiffnesses. In exploration seismology, the first-order theory of Hudson is the most popular. It describes the effective model caused by the presence of a single set of thin, aligned vertical fractures in otherwise isotropic rock. This model is known to be transversely isotropic with a horizontal symmetry axis (HTI). Following the theory, one can invert the effective anisotropy for the crack density and type of fluid infill of fractures, the quantities of great importance for reservoir appraisal and management.Here I compute effective media numerically using the finite element method. I deliberately construct models that contain a single set of vertical, ellipsoidal, non-intersecting and non-interconnected fractures to check validity of the first-order Hudson’s theory and establish the limits of its applicability. Contrary to conventional wisdom that Hudson’s results are accurate up to crack density e ≈ 0.1, I show that they consistently overestimate the magnitudes of all effective anisotropic coefficients ε^(V), δ^(V), and γ^(V). Accuracy of theoretically derived anisotropy depends on the type of fluid infill and typically deteriorates as e grows. While the theory gives | ε^(V)|, |δ^(V)|, |γ^(V)| and close to the upper bound of the corresponding numerically obtained values for randomly distributed liquid-filled fractures, theoretical predictions of ε^(V), δ^(V) are not supported by numerical computations when the cracks are dry. This happens primarily because the first-order Hudson’s theory makes no attempt to account for fracture interaction which contributes to the final result much stronger for gas- than for liquid-filled cracks. I find that Mori-Tanaka’s theory is superior to Hudson’s for all examined crack densities and both types of fluid infill.The paper was presented at the 11th International Workshop on Seismic Anisotropy (11IWSA) held in St. John’s, Canada in 2004.     

Laser Ablation ICP-MS Developments and Trends for 2003

This annual review of laser-ablation inductively coupled plasma-mass spectrometry (LA-ICP-MS) covers the year 2003. Significant advances were made in understanding laser-sample interactions. In particular, research defined the distribution of particle sizes produced by the interplay of laser wavelength, laser pulse width and the gas environment of ablation. A link between particle sizes and elemental and isotopic fractionation at both the ablation site and in the ICP was established. Experimental 15 7 nm and femtosecond laser systems were tested with promising results. The number of applications of LA-ICP-MS in geology and environmental Earth science continued to grow with particular interest in element concentration and isotope ratio profiling of materials, linking composition to time scales. In situ isotopic ratio measurements were increasingly made using multicollector magnetic sector ICP-MS instruments. Other applications of wide interest included bulk sampling of rocks and ores prepared as lithium borate glasses; low level analysis of platinum-group elements, rhenium and gold in sulfides, metal and silicates; in situ uranium-lead zircon geochronology; and melt and fluid inclusion analysis.     

钻孔灌注桩后压浆作用机理与应用

通过对桩侧和桩底后压浆,可大幅度提高钻孔灌注桩单桩承载力.应用桩-土理论分析了后压浆中浆液对土体及桩体的作用,并从分析影响桩基承载力的因素出发,研究了后压浆的作用机理.在渗透率小的地层中后压浆以压密注浆为主,在渗透率大的地层中后压浆以渗透注浆为主.为避免过分压裂土层,在压浆过程中应控制灌浆压力.压浆量是影响后压浆质量和提高单桩承载力的关键因素之一.     

Investigation into the effect of fault properties on wave transmission

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.     

A smoothed particle hydrodynamics modelling of soil–water mixing and resulting changes in average strength

Soil–water interaction is a pivotal process in many underwater geohazards such as underwater landslides where soil sediments gradually evolve into turbidity currents after interactions with ambient water. Due to the large deformations, multiphase interactions and phase changes this involves, investigations from numerical modelling of the transition process have been limited so far. This study explores a simple numerical replication of such soil–water mixing with respect to changes in average strength using smoothed particle hydrodynamics (SPH). A uniform viscoplastic model is used for both the solid-like and fluid-like SPH particles. The proposed numerical solution scheme is verified by single-phase dam break tests and multiphase simple shear tests. SPH combinations of solid-like and fluid-like particles can replicate the clay–water mixture as long as the liquidity index of the solid-like particles is larger than unity. The proposed numerical scheme is shown to capture key features of an underwater landslide such as hydroplaning, water entrainment and wave generation and thus shows promise as a tool to simulate the whole process of subaquatic geohazards involving solid–fluid transition during mass transport.     

From soil aggregates to riverine flocs: a laboratory experiment assessing the respective effects of soil type and flow shear stress on particles characteristics

Particles eroded from hillslopes and exported to rivers are recognized to be composite particles of high internal complexity. Their architecture and composition are known to influence their transport behaviour within the water column relative to discrete particles. To‐date, hillslope erosion studies consider aggregates to be stable once they are detached from the soil matrix. However, lowland rivers and estuaries studies often suggest that particle structure and dynamics are controlled by flocculation within the water column. In order to improve the understanding of particle dynamics along the continuum from hillslopes to the lowland river environment, soil particle behaviour was tested under controlled laboratory conditions. Seven flume erosion and deposition experiments, designed to simulate a natural erosive event, and five shear cell experiments were performed using three contrasting materials: two of them were poorly developed and as such can not be considered as soils, whilst the third one was a calcareous brown soil. These experiments revealed that soil aggregates were prone to disaggregation within the water column and that flocculation may affect their size distribution during transport. Large differences in effective particle size were found between soil types during the rising limb of the bed shear stress sequence. Indeed, at the maximum applied bed shear stress, the aggregated particles median diameter was found to be three times larger for the well‐developed soil than for the two others. Differences were smaller in the falling limb, suggesting that soil aggregates underwent structural changes. However, characterization of particles strength parameters showed that these changes did not fully turn soil aggregates into flocs, but rather into hybrid soil aggregate–floc particles. Copyright © 2013 John Wiley & Sons, Ltd.     

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling the flow of self‐compacting concrete

A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd.     

单片卫星影像在航测测图中的利用

在边境测图过程中,经常遇到较大面积云层问题,使航测测图无法顺利进行。能充分利用的影像资源是第一次全国地理国情普查的影像,有国产高分2号、SPOT-6、资源3号等单片卫星影像,无立体像对影像。通过对国情普查的影像分析选取,选取一试验区域进行试验,用国产高分2号卫星影像与SPOT-6卫星影像组成立体像对测制地物地貌,用Worldview卫星立体像对影像测制同一区域的地物地貌,将两种影像源组成立体模型测制的地物地貌进行比较及精度统计,说明单片卫星影像在航测测图中达到的精度。