Roughness control on hydraulic conductivity in fractured rocks

The influence of joint roughness on the typologies of fluid flow inside fractures is well known and, thanks to experiences in the field of hydraulics, it has been studied from both a physical and mathematical point of view. Nevertheless, the formulations adopted by traditional hydraulic models are hardly applicable in the geological field, because of the difficulty encountered in the roughness parameter estimation. Normally this parameter can be estimated using the joint roughness coefficient (JRC), which considers both the asperity height and its regularity and directional trend. The main advantage in using the JRC arises from the fact that it can easily be obtained from geological-technical surveys and from comparison with the standard Barton profiles. Some relationships have been built up that allow for the estimation of the hydraulic conductivity tensor (an essential parameter for understanding water flow in fractured rock masses), not only as a function of traditional parameters like aperture, spacing, dip and dip direction, etc., but also of joint roughness, precisely expressed in terms of the roughness coefficient. These relationships have been studied initially from a theoretical point of view and then practically, through laboratory investigations.

---

Resumen  Se conoce muy bien la influencia de la rugosidad de las grietas en las tipologías del flujo de fluidos a lo interior de las fracturas y gracias a las experiencias en el campo de hidráulica ha sido posible estudiarla desde puntos de vista matemáticos y físicos. Sin embargo, las formulaciones adoptadas por los modelos hidráulicos tradicionales tienen poca aplicabilidad en el campo geológico debido a la dificultad relacionada con la estimación del parámetro de rugosidad. Normalmente este parámetro puede estimarse usando el coeficiente de rugosidad de grieta (JRC) el cual considera tanto la altura de la aspereza como su regularidad y tendencia direccional. La principal ventaja de utilizar el JRC se deriva del hecho que puede obtenerse fácilmente de levantamientos técnico-geológicos y de la comparación con los perfiles Standard Barton. Se han construido algunas relaciones que permiten la estimación del tensor de conductividad hidráulica (un parámetro esencial para el entendimiento del flujo de agua en masas de roca fracturadas), no solo en función de parámetros tradicionales como apertura, espaciado, buzamiento y dirección de buzamiento, etc., sino también en función de la rugosidad de la grieta estimada con precisión en términos del coeficiente de rugosidad. Estas relaciones se han estudiado inicialmente desde un punto de vista teórico y luego de modo práctico a través de investigaciones de laboratorio.

---

Résumé  L’influence de la rugosité des joints sur les types d’écoulement de fluide dans les fractures est bien connue et a été étudiée aussi bien du point de vue physique que mathématique grace à des expériences menées dans le domaine de l’hydraulique. Cependant les formulations adoptées dans les modèles hydrauliques traditionnels sont difficilement applicables dans le domaine de la géologie à cause de la difficulté rencontrée pour estimer la rugosité. Ce paramètre peut normalement être apprécié grace au coefficient de rugosité du joint (JRC), lequel prend en compte à la fois la hauteur de l’aspérité ainsi que sa régularité et sa direction. Le principal avantage dans l’utilisation du JRC réside dans le fait qu’il peut facilement être obtenu à partir d’études techniques-géologiques et par comparaison avec la classification de Barton. Des relations qui permettent une estimation du tenseur de conductivité hydraulique (un paramètre essentiel pour comprendre l’écoulement de l’eau dans les masses rocheuses fracturées) ont été élaborées, pas seulement en fonction de paramètres traditionnels tels que l’ouverture, l’espacement, l’inclinaison et la direction d’inclinaison, etc , mais aussi en prenant en compte la rugosité des joints à travers le coefficient de rugosité. Ces relations ont initialement été étudiées d’un point de vue théorique puis expérimentalement à travers des recherches en laboratoire.

     

On the quantitative evaluation of metamorphic fluid composition: Verification of the results of physicochemical simulations for water-mineral-rock reactions

This paper reports the results of a comparison of the qualitative physicochemical simulations (by the Winsel program complex) of the composition of the reacting fluid with experimental data on the water-electrolyte (NaCl, HCl, NaOH, and KOH)-mineral (quartz, corundum, microcline, and plagioclase) system and the water-electrolyte-rock (granite and pelite) system at 400–800°C and 1–10 kbar. Constraints are proposed for the temperature, pressure, and the composition of the electrolyte at which the simulation results are consistent with the experimental data.     

Mineral Components, Texture, and Forming Conditions of Hydrothermal Chimney on the East Pacific Rise at 9°-10°N

To characterize the hydrothermal processes of East Pacific rise at 9°-10°N, sulfide mineral compositions, textural, and geochemical features of chimney ores were studied using ore microscope,scanning electron microscope, X-ray diffraction analysis, and electron microprobe techniques. Results show that there are three mineral assemblages for the hydrothermal chimney ores, namely：（i） anhydrite ＋ marcasite ＋ pyrite, （ii） pyrite ＋ sphalerite ＋ chalcopyrite, and （iii） chalcopyrite ＋ bornite ＋ digenite ＋ covellite. Mineral assemblages, zonational features, and geochemical characteristics of the ore minerals indicate that ore fluid temperature changed from low to high then to low with a maximum temperature up to 400 ℃. The chimney is a typical black smoker. The initial structure of the chimney was formed by the precipitation of anhydrites, and later the sulfides began to precipitate in the inner wall.     

Geology and Metallogenesis of the Sawayaerdun Gold Deposit in the Southwestern Tianshan Mountains, Xinjiang, China

The Sawayaerdun gold deposit, located in Wuqia County, Southwest Tianshan, China, occurs in Upper Silurian and Lower Devonian low‐grade metamorphic carbonaceous turbidites. The orebodies are controlled by a series of NE‐NNE‐trending, brittle–ductile shear zones. Twenty‐four gold mineralized zones have been recognized in the Sawayaerdun ore deposit. Among these, the up to 4‐km‐long and 200‐m wide No. IV mineralized zone is economically the most important. The average gold grade is 1–6 g/t. Gold reserves of the Sawayaerdun deposit have been identified at approximately 37 tonnes and an inferred resource of 123 tonnes. Hydrothermal alteration is characterized by silicification, pyritization, arsenopyritization, sericitization, carbonatization and chloritization. On the basis of field evidence and petrographic analysis, five stages of vein emplacement and hydrothermal mineralization can be distinguished: stage 1, early quartz stage, characterized by the occurrence of quartz veins; stage 2, arsenopyrite–pyrite–quartz stage, characterized by the formation of auriferous quartz veinlets and stockworks; stage 3, polymetallic sulfide quartz stage, characterized by the presence of auriferous polymetallic sulfide quartz veinlets and stockworks; stage 4, antimony–quartz stage, characterized by the formation of stibnite–jamesonite quartz veins; and stage 5, quartz–carbonate vein stage. Stages 2 and 3 represent the main gold mineralization, with stage 4 representing a major antimony mineralization episode in the Sawayaerdun deposit. Two types of fluid inclusion, namely H₂O–NaCl and H₂O–CO₂–NaCl types, have been recognized in quartz and calcite. Aqueous inclusions show a wide range of homogenization temperatures from 125 to 340°C, and can be correlated with the mineralization stage during which the inclusions formed. Similarly, salinities and densities of these fluids range for each stage of mineralization from 2.57 to 22 equivalent wt% NaCl and 0.76 to 1.05 g/cm³, respectively. The ore‐forming fluids thus are representative of a medium‐ to low‐temperature, low‐ to medium‐salinity H₂O–NaCl–CO₂–CH₄–N₂ system. The δ³⁴S_CDT values of sulfides associated with mineralization fall into a narrow range of ?3.0 to +2.6‰ with a mean of +0.1‰. The δ¹³C_PDB values of dolomite and siderite from the Sawayaerdun gold deposit range from ?5.4 to ?0.6‰, possibly reflecting derivation of the carbonate carbon from a mixed magmatic/sedimentary source. Changes in physico‐chemical conditions and composition of the hydrothermal fluids, water–rock exchange and immiscibility of hydrothermal fluids are inferred to have played important roles in the ore‐forming process of the Sawayaerdun gold–antimony deposit.     

On the flow structure and turbulence during sweep and ejection events in a wind-tunnel model canopy

Particle image velocimetry (PIV) data obtained in a wind-tunnel model of a canopy boundary layer is used to examine the characteristics of mean flow and turbulence. The vector spacing varies between 1.7 and 2.5 times the Kolmogorov scales. Conditional sampling based on quadrants, i.e. based on the signs of velocity fluctuations, reveals fundamental differences in flow structure, especially between sweep and ejection events, which dominate the flow. During sweeps, the downward flow generates a narrow, highly turbulent, shear layer containing multiple small-scale vortices just below canopy height. During ejections, the upward flow expands this shear layer and the associated small-scale flow structures to a broad region located above the canopy. Consequently, during sweeps the turbulent kinetic energy (TKE), Reynolds stresses, as well as production and dissipation rates, have distinct narrow peaks just below canopy height, whereas during ejections these variables have broad maxima well above the canopy. Three methods to estimate the dissipation rate are compared, including spectral fits, measured subgrid-scale (SGS) energy fluxes at different scales, and direct measurements of slightly underresolved instantaneous velocity gradients. The SGS energy flux is 40–60% of the gradient-based (direct) estimates for filter sizes inside the inertial range, while decreasing with scale, as expected, within the dissipation range. The spectral fits are within 5–30% of the direct estimates. The spectral fits exceed the direct estimates near canopy height, but are lower well above and below canopy height. The dissipation rate below canopy height increases with velocity magnitude, i.e. it has the highest values during sweep and quadrant 1 events, and is significantly lower during ejection and quadrant 3 events. Well above the canopy, ejections are the most dissipative. Turbulent transport during sweep events acts as a source below the narrow shear layer within the canopy and as a sink above it. Transport during ejection events is a source only well above the canopy. The residual term in the TKE transport equation, representing mostly the effect of pressure–velocity correlations, is substantial only within the canopy, and is dominated by sweeps.     

Analysis of the impact of conservation measures on stream flow regime in catchments of the Loess Plateau,China

Catchments in the Loess Plateau have been under the influence of human activities for centuries. In the last four decades, soil conservation measures have accelerated and intensified. These measures were designed to reduce soil erosion, improve agricultural productivity, and enhance environmental quality. It is important to evaluate the effects of these measures on hydrology in order to develop sustainable catchment management plans in the region. This study evaluated changes in stream flow data for four selected catchments in the Loess Plateau following large‐scale soil conservation measures. The non‐parametric Mann–Kendall test was used to identify trends in annual stream flow and the results showed significant downward trends in three of the four catchments. The Pettitt test indicated that a change point occurred in 1978 in these three catchments. Annual precipitation in all the catchments showed no significant trend during the period of record. Comparison of daily flow duration curves for two 20‐year periods (1957–1978) and (1979–2003) showed significant changes in stream flow regime. Reduction in most percentile flows varied between 20 and 45%, and the reduction in low flows was greatest. Overall, the reductions in daily flow were increasing with time, with significant changes occurring in the 1990s. However, it is not clear whether these catchments have seen the full effects of the soil conservation measures, so the results of this study might underestimate the final impact of soil conservation on stream flow regime. Copyright © 2007 John Wiley & Sons, Ltd.     

一种基于图像金字塔光流的特征跟踪方法

推导并实现了一种基于图像金字塔光流的角点特征跟踪方法。实验结果表明,该方法在不同运动幅度和运动方式下的检测跟踪性能较好,能够有效地应用于长序列图像的特征跟踪。     

Incipient sediment motion based on turbulent fluctuations

The current study modifies the representation of the Shields parameter using turbulent strength,i.e.the root-mean-square(rms) fluid velocity.Experiments were done under a steady,uniform flow using eight sediment sizes with particle Reynolds numbers(Rep) ranging between 1.0(fine sediment) and 183.4(coarse sediment).Utilising the peak rms horizontal(u) values,the critical Shields parameter,θ_c,was calculated and a trend similar to the well-established Shields curve was developed.The analysis was extended to the Shields curve obtained based on the critical shear velocity,Reynolds shear stress,and data extracted fro m the oscillating grid-turbulence experiments.Results show that turbulent fluctuations are crucial for the incipient sediment motion and are essentially better predictors than the commonly used critical shear velocity.A quadrant analysis to identify the role of turbulent bursting events in incipient sediment motion also was done where sweeps and ejections are dominant for finer and coarser sediment sizes,respectively.     

Mass flows and river response in rapid uplifting regions – A case of lower Yarlung Tsangpo basin,southeast Tibet,China

The fluvial geomorphology in tectonically active (particularly rapid uplift) regions often undergoes continuous change. The rapid uplift is coincident with high erosion rates; consequently, incised valleys are formed. Mass flows (for example, avalanches, landslides, and debris flows) in incised valleys can markedly influence fluvial processes and even reshape valley geomorphology. However, these processes and long-term evolution corresponding to mass flows require further clarification. Field campaigns were carried out in the region near the Yigong Tsangpo and Palong Tsangpo Rivers (hereafter the Yigong and Palong Rivers), the two largest tributaries of the lower Yarlung Tsangpo River, to examine the feedback between fluvial processes and mass flows. Remote sensing images from recent decades were used to compare the channel morphology before and after typical mass flows (particularly catastrophic ones). The morphology of the lower Yigong River has evidently been impacted by landslides, while that of the Palong River has mainly been shaped by glacial processes and debris flows. At present, the morphology of the latter consists of alternating sections of gorges and wide valleys, with a staircase-like longitudinal profile. The gorge sections exhibit single and deeply incised channels with a high-gradient channel bed and terraces. In contrast, the wide valley sections consist of lakes, braided or anabranching channels, gentle bed gradients, and thick alluvial deposits. Debris flows occur more frequently in gullies in the reaches of the gorge sections and rarely in gullies along the wide valley sections. The occurrence of mass flow events has resulted in an imbalance of the previous (quasi-)equilibrium in the river morphology; however, this has triggered negative feedback that is driving the transient river morphology to a new state of (quasi-)equilibrium.