ICESat-2植被冠层高度和地表高程数据产品用于森林高度提取的效果评价

由美国宇航局(NASA)研制的新一代冰、云和陆地高程卫星(ICESat-2)于2018-09-15发射成功,其搭载的先进地形激光测高系统(ATLAS)采用微脉冲多波束光子计数激光雷达技术,可用于全球高程线采样数据的获取.目前公开发布了9种数据产品,其中包括植被冠层高度和地表高程数据产品(ATL08),为全球森林结构参数...     

一次新型液态CO2播撒效果的数值模拟

利用三维积层混合云人工增雨数值模式对2002年7月11日的一次天气过程进行了由播撒液态CO2引起的微物理量变化及云动力影响的数值模拟。结果表明:播撒后,云中最大上升气流速度增大,由未播撒时的0.37 m/s增大到播撒后的0.54 m/s,播撒使云中出现最大上升速度W的时间比未播撒提前了4 min,表明播撒液态CO2影响了云的动力过程。同时与未播撒相比:云中雨水含量最大值由1.04 g/kg增加到1.40 g/kg;冰粒子含水量的出现提前了88 min,最大值的出现提前了76min;冰粒子浓度的出现提前了72 min,最大值的出现提前了72 min;雪粒子含水量的出现提前了72 min,最大值的出现提前了128 min;云水含量最大值由1.21 g/kg减小到0.87 g/kg。证明了播撒液态CO2后可影响云的微物理过程,从而导致地面降水的增加。     

基于测井频谱分析的松科二井登娄库组地层沉积速率研究

松辽盆地科学钻探工程松科二井获取了白垩纪陆相沉积资料,其中连续、多参数、高分辨率的地球物理测井资料为研究松辽盆地地层的沉积速率变化规律提供了机会。本文利用GR、Th、K、Th/K多种测井数据,采用天文旋回的方法来计算松科二井登娄库组地层的沉积速率,并探讨影响地层沉积的主控因素。多种测井数据的频谱分析结果表明了登娄库组地层记录着米兰科维奇旋回信息,这说明登娄库组沉积过程受天文轨道驱动力影响。Th测井数据对天文旋回信息的敏感性强,综合考虑多种测井数据的测试结果获得了连续、相对准确的地层沉积速率。登娄库组地层沉积速率从下往上整体呈现为由高到低的趋势。登二段至登三段时期地层沉积主要受盆地断陷活动控制,地层整体具有高沉积速率特征,最高达到16.2 cm/ka。随后至登四段时期由于断陷活动逐渐减弱,盆地向区域坳陷的构造格局转化,登四段整体具有沉积速率较低的特征,最低为5.9cm/ka。该研究建立了松科二井登娄库组地层连续的沉积速率剖面,为揭示松辽盆地由断陷活动过渡到坳陷活动的地质规律提供了测井证据。     

Limit states and failure mechanisms of viscous dampers and the implications for large earthquakes

Fluid viscous dampers are used to control story drifts and member forces in structures during earthquake events. These elements provide satisfactory performance at the design‐level or maximum considered earthquake. However, buildings using fluid viscous dampers have not been subjected to very large earthquakes with intensities greater than the design and maximum considered events. Furthermore, an extensive database of viscous damper performance during large seismic events does not exist. To address these issues, a comprehensive analytical and experimental investigation was conducted to determine the performance of damped structures subjected to large earthquakes. A critical component of this research was the development and verification of a detailed viscous damper mathematical model that incorporates limit states. The development of this model and the laboratory and simulation results conclude good correlation with the new model and the damper limit states and provide superior results compared with the typical damper model when considering near collapse evaluation of structures. Copyright © 2010 John Wiley & Sons, Ltd.     

A coupled discrete element model for the simulation of soil and water flow through an orifice

Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd.     

Numerical analysis of the effect of natural microcracks on the supercritical CO2 fracturing crack network of shale rock based on bonded particle models

The problem of predicting the geometric structure of induced fractures is highly complex and significant in the fracturing stimulation of rock reservoirs. In the traditional continuous fracturing models, the mechanical properties of reservoir rock are input as macroscopic quantities. These models neglect the microcracks and discontinuous characteristics of rock, which are important factors influencing the geometric structure of the induced fractures. In this paper, we simulate supercritical CO₂ fracturing based on the bonded particle model to investigate the effect of original natural microcracks on the induced‐fracture network distribution. The microcracks are simulated explicitly as broken bonds that form and coalesce into macroscopic fractures in the supercritical CO₂ fracturing process. A calculation method for the distribution uniformity index (DUI) is proposed. The influence of the total number and DUI of initial microcracks on the mechanical properties of the rock sample is studied. The DUI of the induced fractures of supercritical CO₂ fracturing and hydraulic fracturing for different DUIs of initial microcracks are compared, holding other conditions constant. The sensitivity of the DUI of the induced fractures to that of initial natural microcracks under different horizontal stress ratios is also probed. The numerical results indicate that the distribution of induced fractures of supercritical CO₂ fracturing is more uniform than that of common hydraulic fracturing when the horizontal stress ratio is small.     

Nb–Ta fractionation induced by fluid‐rock interaction in subduction‐zones: constraints from UHP eclogite‐ and vein‐hosted rutile from the Dabie orogen,Central‐Eastern China

Niobium and Ta concentrations in ultrahigh‐pressure (UHP) eclogites and rutile from these eclogites and associated high pressure (HP) veins were used to study the behaviour of Nb–Ta during dehydration and fluid‐rock interaction. Samples were collected through a ～2 km profile at the Bixiling complex in the Dabie orogenic belt, Central‐Eastern China. All but one eclogite away from veins (EAVs) display nearly constant Nb/Ta ratios ranging from 16.1 to 19.2, with an average of 16.9 ± 0.8 (2 SE), similar to that of their gabbroic protolith from the Yangtze Block. Nb/Ta ratios of rutile from the EAVs range from 12.7 to 25.3 among different individual grains, with the average values close to those of the corresponding bulk rocks. These observations show that Nb and Ta were not significantly fractionated by prograde metamorphism up to eclogite facies when no significant fluid‐rock interaction occurs. In contrast, Nb/Ta ratios of rutile from eclogites close to veins (ECVs) are highly variable from 17.8 to 49.8, which are systematically higher (by up to 17) than those of rutile from the veins. These observations demonstrate that Nb and Ta were mobilized and fractionated during localized fluid flow and intensive fluid‐rock interaction. This is strongly supported by Nb/Ta zoning patterns in single rutile grains revealed by in situ LA‐ICP‐MS analysis. Ratios of Nb/Ta in the ECV‐hosted rutile decrease gradually from cores towards rims, whereas those in the EAV‐hosted rutile are nearly invariable. Furthermore, the vein rutile shows Nb/Ta zoning patterns that are complementary to those in rutile from their immediate hosts (ECVs), suggesting an internal origin for the vein‐forming fluids. The Nb/Ta ratios of such fluids evolved from low values at the early stage of subduction to higher values at later supercritical conditions with increased temperature and pressure. Quantitative modelling was conducted to constrain the compositional evolution of metamorphic fluids during dehydration and fluid‐rock interaction focusing on Nb–Ta distribution. The modelling results based on our proposed multistage fluid phase evolution path can essentially reproduce the natural observations reported in the present study.     

Geodynamic modelling of the Century deposit,Mt Isa Province,Queensland

This paper is concerned with an example of quantitative modelling of orebody formation as a guide to reducing the risk for future mineral exploration. Specifically, the paper presents a detailed 3–D numerical model for the formation of the Century zinc deposit in northern Queensland. The model couples fluid flow with deformation, thermal transport and chemical reactions. The emphasis of the study is a systems approach where the holistic mineralising system is considered rather than concentrating solely on the mineral deposit. In so doing the complete plumbing system for mineralisation is considered with a view to specifying the critical conditions responsible for the ore deposit occurring where it does and having the size and metal grades that are observed. The numerical model is based on detailed geological, tectonic, isotopic and mineralogical data collected over the past 20 years. The conclusions are that the Century zinc deposit is located where it is because of the following factors: (i) a thermal anomaly is associated with the Termite Range Fault due to advection of heat from depth by fluid flow up the Termite Range Fault; (ii) bedding‐plane fissility in the shale rocks hosting the Century zinc deposit has controlled the wavelength and nature of D₁ folding in the vicinity of the deposit and has also controlled increases in permeability due to hydrofracture of the shales; such hydrofracture is also associated with the production of hydrocarbons as these shales passed through the ‘oil‐window’; (iii) Pb–Zn leached from crustal rocks in the stratigraphic column migrated up along faults normal to the Termite Range Fault driven by topographic relief associated with inversion at the end of the Isan Orogeny; these fluids mixed with H₂S derived at depth moving up the Termite Range Fault to mix with the crustal fluids to precipitate Pb–Zn in a plume downstream from the point of mixing. Critical factors to be used as exploration guides are high temperatures, carbonaceous fissile shales now folded into relatively tight D₁ folds, fault‐controlled plumbing systems that enable fluid mixing, depletion of metals upstream of the deposit and,in particular,a very wide Fe‐depletion halo upstream of the deposit.     

Silurian and Devonian systems

The Ediacara mineral field is situated 30 km W of Beltana on the western margins of the Flinders Ranges, South Australia, and consists of silver‐lead and copper deposits in lower Cambrian carbonate rocks that contain anomalous base‐metal contents throughout the Adelaide Geosyncline. The lower Cambrian rocks, which consist of the basal Parachilna Formation and overlying Ajax Limestone, rest disconformably on the Precambrian, and at Ediacara occupy a shallow N‐S elongate syncline near the hinge zone of the Adelaide Geosyncline. The main primary ore minerals of the silver‐lead mineralization are galena and pyrite, with very minor chalcopyrite and sphalerite, and rare tetrahedrite and pearceite. The gangue consists mainly of silica (both chalcedony and quartz), with minor dolomite and rare barite. The mineralization is stratabound and occurs in conformable zones, the lowest of which commences about 30–50 m above the base of the Cambrian sequence. The host to the silver‐lead mineralization, the Ajax Limestone, can be subdivided into three units which represent a set of lithologies, structures and organic traces indicative of a shallow near‐shore carbonate environment. The silver‐lead mineralization is mainly present in sandy and laminated dolomites which were deposited in an environment ranging from sub‐tidal to bar and channel and tidal flat, respectively. Four types of mineralization have been recognized; disseminated sulphides of syngenetic and/or diagenetic origin and epigenetic concentrations along stylolites, in veins and as breccia fillings. Post‐depositional solution activity has affected a large proportion of the carbonate sequence. The effects of this activity range from stylolites through stylobreccias to solution collapse breccias. The epigenetic concentrations of mineralizations have apparently been formed by the remobilization of the disseminated sulphides during solution activity. The ore and gangue minerals of the epigenetic mineralization display both euhedral forms and distinct colloform banding, and framboidal textures have also been observed in both pyrite and galena. There is evidence of repeated episodic precipitation and no simple paragenetic sequence can be recognized. Fluid inclusions in silica and dolomite associated with the epigenetic mineralization have homogenization temperatures of 159 to 199°C and freezing temperatures that indicate the fluids to be saline brines containing NaCl with CaCl₂ and/or MgCl₂. Sulphur isotope analyses show a range of 8³⁴S values from ‐12.5 to +8.6 per mil, with no evidence of significant differences between the four types of mineralization. The data suggest deposition of the disseminated sulphides as a result of biological reduction of seawater sulphate in a system partially open with respect to sulphate supply. Subsequent remobilization of sulphides apparently involved little or no sulphur isotope fractionation. The Ediacara silver‐lead deposits have many features in common with Mississippi Valley‐type lead‐zinc deposits and appear to have similarities in terms of genesis, in that the epigenetic mineralization has been formed as a result of post‐depositional solution activity during diagenesis in a sedimentary basin. The scale of transport of the metals deposited as the epigenetic mineralization at Ediacara appears, however, to have been very much less than that of the metals in other Mississippi Valley‐type deposits.     

The role of eddies in solute transport and recovery in rock fractures: Implication for groundwater remediation

A better understanding of solute transport and retention mechanism in rock fractures has been challenging due to difficulty in their direct observations in microscale rough‐walled fractures. Six representative troughs in a rough‐walled fracture were selected for microscale observations of eddy formation with increasing flow velocity and its effect on spatiotemporal changes of solute concentration. This experimental study was enabled by a microscale visualization technique of micro particle image velocimetry. With increasing flow velocity (Re ≤ 2.86), no eddies were generated, and solutes along the main streamlines transported rapidly, whereas those near the wall moved slowly. A larger amount of solutes remained trapped at all troughs at Re = 2.86 than Re < 1. For Re = 8.57, weak eddies started to be developed at the troughs on the lee side, which little contributed to overall solute flushing in the fracture. Accordingly, a large of amount of water was needed for solute flushing. The flow condition of 1 < Re < 10, before a full development of eddies, was least favourable in terms of time and amount of remediation fluid required to reach a target concentration. After large eddies were fully developed at troughs on the lee side for Re = 17.13, solutes were substantially reduced by eddies with less amount of water. Fully developed eddies were found to enhance solute transport and recovery, as opposed to a general consensus that eddies trap and delay solutes. Direct inflow into troughs on the stoss side also made a great contribution to solute flushing out of the troughs. This study indicates that fully developed eddies or strong inflows at troughs are highly possible to form for Re > 10 and this flow range could be favourable for efficient remediation.