Sedimentation–consolidation of a double porosity material

This paper studies the sedimentation–consolidation of a double porosity material, such as lumpy clay. Large displacements and finite strains are accounted for in a multidimensional setting. Fundamental equations are derived using a phenomenological approach and non-equilibrium thermodynamics, as set out by Coussy [Coussy, Poromechanics, Wiley, Chichester, 2004]. These equations particularise to three non-linear partial differential equations in one dimensional context. Numerical implementation in a finite element code is currently being undertaken.     

Importance of borehole deviation surveys for monitoring of hydraulic fracturing treatments

During seismic monitoring of hydraulic fracturing treatment, it is very common to ignore the deviations of the monitoring or treatment wells from their assumed positions. For example, a well is assumed to be perfectly vertical, but in fact, it deviates from verticality. This can lead to significant errors in the observed azimuth and other parameters of the monitored fracture‐system geometry derived from microseismic event locations. For common hydraulic fracturing geometries, a 2° deviation uncertainty on the positions of the monitoring or treatment well survey can cause a more than 20° uncertainty of the inverted fracture azimuths. Furthermore, if the positions of both the injection point and the receiver array are not known accurately and the velocity model is adjusted to locate perforations on the assumed positions, several‐millisecond discrepancies between measured and modeled SH‐P traveltime differences may appear along the receiver array. These traveltime discrepancies may then be misinterpreted as an effect of anisotropy, and the use of such anisotropic model may lead to the mislocation of the detected fracture system. The uncertainty of the relative positions between the monitoring and treatment wells can have a cumulative, nonlinear effect on inverted fracture parameters. We show that incorporation of borehole deviation surveys allows reasonably accurate positioning of the microseismic events. In this study, we concentrate on the effects of horizontal uncertainties of receiver and perforation positions. Understanding them is sufficient for treatment of vertical wells, and also necessary for horizontal wells.     

Geochemistry and Fluid Inclusions Analysis of Vein Quartz in the Multiple Hydrothermal Systems of Mankayan Mineral District,Philippines

Several high‐sulfidation epithermal gold orebodies in the Mankayan Mineral District were formed in an environment that has been already affected by earlier porphyry‐type mineralization. This study reports the geologic and geochemical characteristics of the Carmen and Florence epithermal orebodies, which are located in the south of the Lepanto main enargite–gold orebody. The gold‐bearing epithermal quartz veins in the Carmen and Florence areas are of two types: (i) the enargite‐rich veins and (ii) the quartz–pyrite–gold (QPG) veins. The two types of veins are mainly hosted by the Cretaceous Lepanto Metavolcanics basement rocks, with minor veins cutting the Pleistocene Imbanguila Dacite Pyroclastics. The mineral assemblages and homogenization temperatures of fluid inclusions indicate that the Carmen and Florence orebodies were deposited by fluids varying from high to very high sulfidation state. The enargite and QPG epithermal veins of Carmen and Florence cut porphyry‐type quartz veinlet stockworks and veins that host polyphase hypersaline fluid inclusions that did not homogenize at or below 400°C. These high‐temperature quartz exhibits distinctly different mineral chemistry from the quartz of the QPG and enargite‐rich epithermal veins. In particular, the Ti content of quartz of the porphyry‐type veinlet stockwork is elevated (>100 ppm), whereas the Ti concentration of the epithermal vein quartz crystals are below detection limits. The Fe concentration of quartz is high in epithermal vein quartz (>300 ppm), whereas nearly undetected in the porphyry‐type stockwork veinlet quartz. Multiple generations of quartz with different mineral chemistry, fluid inclusions morphology, temperature, salinity and bulk gas compositions, and stable isotopic ratios indicate the variable hydrothermal conditions throughout the mineralization history of the Mankayan District. The temperature, pH, sulfidation state, oxidation state, and fluid composition vary among the orebodies in Carmen and Florence areas. Furthermore, the characteristics of earlier alteration affected the apparent characteristics of subsequent mineralization.     

Slab rollback and microcontinent subduction in the evolution of the Zambales Ophiolite Complex(Philippines):A review

New radiolarian ages show that the island arc-related Acoje block of the Zambales Ophiolite Complex is possibly of Late Jurassic to Early Cretaceous age.Radiometric dating of its plutonic and volcanichypabyssal rocks yielded middle Eocene ages.On the other hand,the paleontological dating of the sedimentary carapace of the transitional mid-ocean ridge-island arc affiliated Coto block of the ophiolite complex,together with isotopic age datings of its dikes and mafic cumulate rocks,also yielded Eocene ages.This offers the possibility that the Zambales Ophiolite Complex could have:(1)evolved from a Mesozoic arc(Acoje block)that split to form a Cenozoic back-arc basin(Coto block),(2)through faulting,structurally juxtaposed a Mesozoic oceanic crust with a younger Cenozoic lithospheric fragment or(3)through the interplay of slab rollback,slab break-off and,at a later time,collision with a microcontinent fragment,caused the formation of an island arc-related ophiolite block(Acoje)that migrated trench-ward resulting into the generation of a back-arc basin(Coto block)with a limited subduction signature.This Meso-Cenozoic ophiolite complex is compared with the other oceanic lithosphere fragments along the western seaboard of the Philippines in the context of their evolution in terms of their recognized environments of generation.     

The celestial mechanics approach: application to data of the GRACE mission

The celestial mechanics approach (CMA) has its roots in the Bernese GPS software and was extensively used for determining the orbits of high-orbiting satellites. The CMA was extended to determine the orbits of Low Earth Orbiting satellites (LEOs) equipped with GPS receivers and of constellations of LEOs equipped in addition with inter-satellite links. In recent years the CMA was further developed and used for gravity field determination. The CMA was developed by the Astronomical Institute of the University of Bern (AIUB). The CMA is presented from the theoretical perspective in (Beutler et al. 2010). The key elements of the CMA are illustrated here using data from 50 days of GPS, K-Band, and accelerometer observations gathered by the Gravity Recovery And Climate Experiment (GRACE) mission in 2007. We study in particular the impact of (1) analyzing different observables [Global Positioning System (GPS) observations only, inter-satellite measurements only], (2) analyzing a combination of observations of different types on the level of the normal equation systems (NEQs), (3) using accelerometer data, (4) different orbit parametrizations (short-arc, reduced-dynamic) by imposing different constraints on the stochastic orbit parameters, and (5) using either the inter-satellite ranges or their time derivatives. The so-called GRACE baseline, i.e., the achievable accuracy of the GRACE gravity field for a particular solution strategy, is established for the CMA.     

Effect of channel deepening on tidal flow and sediment transport—part II: muddy channels

Natural tidal channels often need deepening for navigation purposes (larger vessels). The depth increase may lead to tidal amplification, salt intrusion over longer distances, and increasing sand and mud import. Increasing fine sediment import, in turn, may start a process in which the sediment concentration progressively increases until the river becomes hyper-turbid, which may lead to increased dredging volumes and to decreased ecological values. These effects can be modeled and studied using detailed 3D models. Reliable simplified models for a first quick engineering evaluation are however lacking. In this paper, we apply both simplified and detailed 3D models to analyze the effects of channel deepening in prismatic and weakly converging tidal channels with saturated mud flow. The objective is to gain quantitative understanding of the effects of channel deepening on mud transport. We developed a simplified tidal mud model describing most relevant processes and effects in saturated mud flows with only minor horizontal transport gradients (quasi uniform conditions). The simplified model is not valid for non-saturated mud flow conditions. This model can either be used in standalone mode or in post-processing mode with computed near-bed velocities from a 3D hydrodynamic model as an input. The standalone model has been compared to various field data sets. Mud transport processes in the mouth region of muddy tidal channels can be realistically represented by the simplified model, if sufficient salinity and sediment data are available for calibration. The simulation of tidal mud transport and the behavior of an estuarine turbidity maximum (ETM) in saturated and non-saturated mud flow conditions cannot be represented by the simplified model and requires the application of a detailed 3D model.     

Experimental verification of effective anisotropic crack theories in variable crack aspect ratio medium

Physical modelling of cracked/fractured media using downscaled laboratory experiments has been used with great success as a useful alternative for understanding the effect of anisotropy in the hydrocarbon reservoir characterization and in the crustal and mantle seismology. The main goal of this work was to experimentally verify the predictions of effective elastic parameters in anisotropic cracked media by Hudson and Eshelby–Cheng's effective medium models. For this purpose, we carried out ultrasonic measurements on synthetic anisotropic samples with low crack densities and different aspect ratios. Twelve samples were prepared with two different crack densities, 5% and 8%. Three samples for each crack density presented cracks with only one crack aspect ratio, whereas other three samples for each crack density presented cracks with three different aspect ratios in their composition. It results in samples with aspect ratio values varying from 0.13 to 0.26. All the cracked samples were simulated by penny‐shaped rubber inclusions in a homogeneous isotropic matrix made with epoxy resin. Moreover, an isotropic sample for reference was constructed with epoxy resin only. Regarding velocity predictions performed by the theoretical models, Eshelby–Cheng shows a better fit when compared with the experimental results for samples with single and mix crack aspect ratio (for both crack densities). From velocity values, our comparisons were also performed in terms of the ε, γ, and δ parameters (Thomsen parameters). The results show that Eshelby–Cheng effective medium model fits better with the measurements of ε and γ parameters for crack samples with only one type of crack aspect ratio.