泾阳南塬蒋刘4#滑坡特征及成因机制

2016年3月6日,泾阳南塬蒋刘村发生了约2×104 m3的小型黄土滑坡。滑坡发生后,在对滑坡现场进行详细的地质调查基础上,结合无人机航测和三维激光扫描综合手段,对滑坡的基本特征及成因机理进行了初步分析。结果表明:长期引水灌溉是滑坡发生的基本条件;坡脚滞水软化作用、与灌溉渠较短的渗透距离、坡脚冻结黄土春季解冻是导致滑体发生的直接原因;底部饱和泥化黄土的铺垫及导流作用,是导致其滑程较远的根本原因。该类滑坡近几年来发生频率较高,为泾阳南塬滑坡演化后期的主要发生类型,具有一定的特殊性,值得深入研究。     

大兴安岭呼中钼多金属矿床成矿流体特征及成矿机制

大兴安岭呼中镇钼多金属矿床产于早白垩世碎裂花岗岩中,矿体形态似层状、脉状或透镜状。成矿阶段分为辉钼矿阶段、铜铅锌阶段、方解石阶段等3个阶段。本矿床中的包裹体主要有气液H_2O包裹体、气相H_2O包裹体、含子矿物H_2O包裹体,CO_2包裹体、CO_2-H_2O包裹体等五大类。辉钼矿阶段均一温度为209.2℃413.4℃,盐度ω(NaCleq)为0.53%16.53%,平均值7.3%,属中低盐度,成矿压力在53.3 68.5 MPa,相应的成矿深度为2 052 2 807 m;铜铅锌阶段为191.5℃351.0℃,盐度ω(NaCleq)为0.88%34.24%,平均值3.96%,属低盐度;方解石阶段207.2℃233.5℃,盐度ω(NaCleq)为3.55%4.18%,平均值3.77%,为低盐度。从辉钼矿阶段到铜铅锌阶段再到方解石阶段,温度及盐度逐渐降低。辉钼矿成矿压力在53.3 72.9 MPa,成矿深度为2 052 2 807 m。成矿过程中发生了CO_2-H_2O不混溶作用,这是导致辉钼矿沉淀的主要机制。矿床类型属中低温浅成热液型矿床。     

从成岩-蚀变特征探讨塔深1井白云岩储层成因

塔深1井钻遇的寒武系地层属于中寒武统阿瓦塔格组和上寒武统丘里塔格下亚群。阿瓦塔格组以灰色-深灰色粉晶-细晶藻白云岩为主,丘里塔格下亚群以灰色-浅灰色晶粒结构白云岩与残余亮晶结构白云岩为主。藻白云岩在阴极发光下以明亮的中等强度砖红色-桔黄色为主要特点,具有相对偏重的δ18O值(-6.3‰~-4.5‰,平均为-5.3‰,VPDB)、相对较高的稀土总量以及与正常海相碳酸盐岩相似的稀土配分样式。晶粒白云岩与残余亮晶结构白云岩在阴极发光下以偏暗的中等强度蓝紫色-砖红色为主要特点,具有相对偏轻的δ18O值(-9.2‰~-5.7‰,平均为-7.5‰,VPDB)、相对偏低的稀土总量、Eu正异常与重稀土明显亏损为特点的稀土配分模式。阿瓦塔格组和丘里塔格下亚群白云岩在上述岩石学特征和地球化学特征上的差异表明储层发育的成岩-蚀变过程存在差异。研究认为,埋藏较深的阿瓦塔格组优质储层的形成受多种因素控制,其中准同生期大气淡水作用下石膏结核的溶解及非稳定碳酸盐矿物的溶解奠定了基础,后期成岩改造及流体作用进行了强化,局部遭充填破坏。     

华北克拉通北缘东段“开原岩群”的解体及重新认识

根据近年来的同位素年代学资料,对分布于辽北地区的一套构造变质岩系（"开原岩群"或"清河镇岩群"）进行了重新认识和划分,并将其重新定义为清河构造混杂岩,由中太古代、新太古代、中元古代、新元古代、二叠纪等不同时代、不同构造环境的地质体组成.其中,中太古代、新太古代、中元古代岩石来自华北克拉通,而新元古代、二叠纪岩石则来自兴蒙造山带.将原沈家堡子岩组重新划分为新太古代变质表壳岩和中-新太古代变质深成岩,原板石沟岩组（或芦家堡子岩组）厘定为中元古代石门岩组,原照北山岩组重新划分为新元古代南平岩片和晚二叠世照北山岩组,将原佟家屯岩组和尖山子火山岩合并为晚二叠世佟家屯岩组.最后,提出了清河构造混杂岩带为华北克拉通与兴蒙造山带的界线以及古亚洲洋在华北板块北缘东段的最终闭合时间为中三叠世早期（约245Ma）的认识.     

Investigation of gas content of organic-rich shale:A case study from Lower Permian shale in southern North China Basin,central China

Measuring gas content is an essential step in estimating the commerciality of gas reserves. In this study,eight shale core samples from the Mouye-1 well were measured using a homemade patented gas desorption apparatus to determine their gas contents. Due to the air contamination that is introduced into the desorption canister, a mathematical method was devised to correct the gas quantity and quality.Compared to the chemical compositions of desorbed gas, the chemical compositions of residual gas are somewhat different. In residual gas, carbon dioxide and nitrogen record a slight increase, and propane is first observed. This phenomenon may be related to the exposure time during the transportation of shale samples from the drilling site to the laboratory, as well as the differences in the mass, size and adsorptivity of different gas molecules. In addition to a series of conventional methods, including the USBM direct method and the Amoco Curve Fit(ACF) method, which were used here for lost gas content estimation, a Modified Curve Fit(MCF) method, based on the 'bidisperse' diffusion model, was established to estimate lost gas content. By fitting the ACF and MCF models to gas desorption data, we determined that the MCF method could reasonably describe the gas desorption data over the entire time period, whereas the ACF method failed. The failure of the ACF method to describe the gas desorption process may be related to its restrictive assumption of a single pore size within shale samples. In comparison to the indirect method, this study demonstrates that none of the three methods studied in this investigation(USBM, ACF and MCF) could individually estimate the lost gas contents of all shale samples and that the proportion of free gas relative to total gas has a significant effect on the estimation accuracy of the selected method. When the ratio of free gas to total gas is lower than 45%, the USBM method is the best for estimating the lost gas content, whereas when the ratio ranges from 45% to 75% or is more than 75%, the ACF and MCF methods, are the best options respectively.     

Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.     

Effects of Wave–Current Interaction on Salt Intrusion During a Typhoon Event in a Highly Stratified Estuary

Wave–current interaction (WCI) is important in modulating hydrodynamics and water mixing in estuaries, and thereby the transport of water-borne materials. However, the effects of WCI on salt transport and salt intrusion in estuaries during storm events have been rarely examined. In the present study, we use a coupled atmosphere–ocean–wave–sediment transport (COAWST) modeling system to investigate the effects of WCI on salt intrusion in the highly stratified Modaomen Estuary during Typhoon Hagupit (2008). The model is validated by the measured wave, water elevation, and surface salinity data, and several diagnostic model experiments are conducted. WCI increases the storm surge by 0.8 m at the peak surge (25% of the total surge height). The wave-breaking-induced momentum flux and the Stokes drift increase the magnitude of the landward flow by 0.3 m s^?1 (30% of the total landward flow). In addition, the waves increase water mixing by 2–4 times compared with that without waves. Hence, WCI significantly increases the landward advective salt transport and decreases the steady shear transport. The net effect of the WCI is a significant increase of salt import and salt intrusion during the typhoon event. However, in the aftermath of the storm, the imported salt water is rapidly flushed out by the increased river discharge, and the estuary regains its stratification within one day.     

Numerical analysis of effect of baffle configuration on impact force exerted from rock avalanches

In mountainous areas, channelized rock avalanches swarm downslope leading to large impact forces on building structures in residential areas. Arrays of rock avalanche baffles are usually installed in front of rigid barriers to attenuate the flow energy of rock avalanches. However, previous studies have not sufficiently addressed the mechanisms of interaction between the rock avalanches and baffles. In addition, empirical design approaches such as debris flow (Tang et al., Quat Int 250:63–73, 2012), rockfall (Spang and Rautenstrauch, 1237–1243, 1988), snow avalanches (Favier et al., 14:3–15, 2012), and rock avalanches (Manzella and Labiouse, Landslides 10:23–36, 2013), which are applied in natural geo-disasters mitigation cannot met construction requirements. This study presents details of numerical modeling using the discrete element method (DEM) to investigate the effect of the configuration of baffles (number and spacing of baffle columns and rows) on the impact force that rock avalanches exert on baffles. The numerical modeling is firstly conducted to provide insights into the flow interaction between rock avalanches and an array of baffles. Then, a modeling analysis is made to investigate the change pattern of the impact force with respect to baffle configurations. The results demonstrate that three crucial influencing factors (baffle row numbers, baffle column spacing, and baffle row spacing) have close relationship with energy dissipation of baffles. Interestingly, it is found that capacity of energy dissipation of baffles increases with increasing baffle row numbers and baffle row spacing, while it decreases with increasing baffle column spacing. The results obtained from this study are useful for facilitating design of baffles against rock avalanches.     

Application of incomplete similarity theory to the estimation of the mean velocity of debris flows

The mean velocity of debris flow is one of the most important parameters in the design of mitigation structures and in quantitative risk analysis. This study develops a model to predict the mean debris flow velocity observed in the field by applying the incomplete similarity argument. An equation for estimating the Darcy-Weisbach resistance coefficient for debris flows with a volumetric sediment concentration larger than 0.19 is accordingly derived using 128 sets of observation data from nine Chinese gullies, in which both the effect of the volumetric sediment concentration and channel slope on resistance are considered. The derived equation is then verified and compared against five previously published equations by using 61 sets of published observation data from six gullies located in four countries. The applications of the proposed equation are discussed, and the improvements made using the proposed equation are clearly very significant when compared with the previously published equations.