The motion of sediment-water mixtures during intense bedload transport: computer simulations

A new model, which couples fluid and particle dynamics, has been developed to study the motion of the sediment-water mixture during intense bedload transport, including the velocity profiles of both sediment and water, the roughness length of an upper plane bed and the thickness of moving sediment layers. Standard mixing length theory is used to model the motion of water above the boundary between the overlying water and the sediment-water mixture. The turbulent flow within the moving sediment layers is described by a shear stress model, in which the effective viscosity of the flowing water is proportional to the velocity difference between the fluid and the sediment. The particle dynamics method, in which the equations of motion of each of many particles are solved directly, is applied to model the movement of sediment particles. The particle-fluid interaction is expressed by a velocity-squared fluid drag force exerted on each sediment particle. Both computer simulation results and theoretical analysis have shown that the velocities of both sediment and fluid during intense sediment transport decrease exponentially with depth in the top layers of a fast-moving sediment—water mixture. The thickness of the moving sediment layers, obtained from the computer simulation results, is proportional to the shear stress, which agrees with previous experimental observations.     

Numerical simulation of turbidity current flow and sedimentation: I. Theory

A computer-based numerical model of turbidity current flow and sedimentation is presented that integrates geological observations with basic equations for fluid and sediment motion. The model quantifies those aspects of turbidity currents that make them different from better-understood fluvial processes, including water mixing across the upper flow boundary and the interactions between the suspended-sediment concentration and the flow dynamics and sedimentation. The model includes three numerical components: (1) a layer-averaged three-equation flow model for tracing downslope flow evolution using continuity and momentum equations, (2) a sedimentation/fluidization model for tracing sediment-size fractionation in sedimenting multicomponent suspensions and (3) a concentration-viscosity model for quantifying the changes in resistance of such suspensions toward fluid and sediment motion. The model traces the evolution of a model turbidity current in terms the layer-averaged flow velocity, flow thickness, sediment concentration distribution, and the rate of sedimentation and sediment size fractionation. It generates synthetic turbidites with downslope variations in thickness and grain-size structuring at each point along the flow path. This study represents an effort to evaluate quantitatively the effects of basin geometry, sediment supply and sediment properties on the mechanics of turbidity current flow and sedimentation and on the geometry and grain size characteristics of the resulting deposits.     

扇三角洲砂体几何形态沉积数值模拟方法研究

扇三角洲是陆相碎屑岩重要的沉积类型之一.利用砂体沉积过程数值模拟方法可以预测扇三角洲砂体的几何形态.在利用泥沙冲淤动力学模式建立扇三角洲沉积过程的数学方程后,根据现代三角洲沉积特征,设计了其沉积过程的模拟条件,计算域长100 km、宽50 km,流量按50年一遇洪水设计,模拟过程到2 000年时,扇三角洲沉积过程基本达到平衡状态,此时前缘复合砂体最大厚度约39 m.通过对扇三角洲前缘砂体几何形态的模拟研究,总结出扇三角洲沉积体储层建筑结构的模拟与预测方法.计算过程中可识别出4种沉积砂体,包括水下分流河道(砂体的平均长宽比为2.68,平均宽厚比为79.8)、河口砂坝(平均长宽比为2.02,平均宽厚比为68.2)、远砂坝(平均长宽比为1.65,平均宽厚比为58.3)、水下溢岸沉积.主要沉积单元砂体几何参数之间的相关关系较好.实验结果表明,利用泥沙冲淤动力模式可以较好地揭示扇三角洲发育过程,进而可以预测扇三角洲砂体的几何形态.     

Fully coupled generalised hybrid-mixed finite element approximation of two-phase two-component flow in porous media. Part II: numerical scheme and numerical results

We consider the modeling and simulation of compositional two-phase flow in a porous medium, where one phase is allowed to vanish or appear. The modeling of Marchand et al. (in review) leads to a nonlinear system of two conservation equations. Each conservation equation contains several nonlinear diffusion terms, which in general cannot be written as a function of the gradients of the two principal unknowns. Also the diffusion coefficients are not necessarily explicit local functions of them. For the generalised mixed finite elements approximation, Lagrange multipliers associated to each principal unknown are introduced, the sum of the diffusive fluxes of each component is explicitly eliminated and the static condensation leads to a “global” nonlinear system of equations only in the Lagrange multipliers also including complementarity conditions to cope with vanishing or appearing phases. After time discretisation, this system can be solved at each time step using a semi-smooth Newton method. The static condensation involves “local” nonlinear systems of equations associated to each element, solved also by a semismooth Newton method. The algorithm is successfully applied to 1D and 2D examples of water–hydrogen flow involving gas phase appearance and disappearance.     

A parallel, implicit, cell‐centered method for two‐phase flow with a preconditioned Newton–Krylov solver

A new parallel solution technique is developed for the fully implicit three‐dimensional two‐phase flow model. An expandedcell‐centered finite difference scheme which allows for a full permeability tensor is employed for the spatial discretization, and backwardEuler is used for the time discretization. The discrete systems are solved using a novel inexact Newton method that reuses the Krylov information generated by the GMRES linear iterative solver. Fast nonlinear convergence can be achieved by composing inexact Newton steps with quasi‐Newton steps restricted to the underlying Krylov subspace. Furthermore, robustness and efficiency are achieved with a line‐search backtracking globalization strategy for the nonlinear systems and a preconditioner for each coupled linear system to be solved. This inexact Newton method also makes use of forcing terms suggested by Eisenstat and Walker which prevent oversolving of the Jacobian systems. The preconditioner is a new two‐stage method which involves a decoupling strategy plus the separate solutions of both nonwetting‐phase pressure and saturation equations. Numerical results show that these nonlinear and linear solvers are very effective.     

Numerical fixed domain mapping solution of free-surface flows coupled with an evolving interior field

The evolution of a gravity-driven free-surface flow of varying horizontal extent which couples with a field evolving within the flow is solved using a finite difference discretization of a mapping of the problem onto the unit square. Since the size of the solution domain may show several orders of magnitude of variation, while the normalized geometry of the domain and the internal field may not vary significantly, this procedure avoids excessively fine or coarse discretizations, as well as interpolations at the boundary. The parabolic and hyperbolic evolution equations for the internal field are considered. The evolution of the coupled system is solved by an implicit marching scheme. The discretizations in space and in time are accurate to second order. Multipoint upwinding is used to avoid an instability arising advective terms are large. The evolution equations are nonlinear, and are solved using a nested Newton–Raphson procedure. The nesting is achieved by using successively better approximations to the ture evolution equations. The matrix equation that arises is solved by a conjugate-gradient-like (ORTHOMIN) iteration procedure with an incomplete Cholesky factorization preconditioning. The method has a wide variety of potential applications in the earth sciences, with the ability to describe glacier flow, lava flow, avalanching and landslides. Some calculations of the thermomechanical evolution of ice-sheets are given as illustrations, and the possible existence of thermally induced instabilities is considered.     

湖北省晚二叠世龙潭组沉积期岩相古地理格局与成煤环境

晚二叠世龙潭组沉积期是湖北省境内一次重要的成煤时期,由于煤层的不稳定性,其煤层层数、厚度、煤质等存在明显差异。通过岩性、煤层及古生物研究对比,查明龙潭组由下向上共出现5个煤层,煤层的形成与沉积特征明显受岩相古地理条件所制约。在龙潭组沉积期,湖北省境内总体为"两陆夹一盆"的基本格局,形成潮坪-泻湖-障壁岛沉积体系,该沉积体系是华南海北部沉积区的一个组成部分,由秦淮古陆、海湾泻湖及江南古陆组成,沉积区主要为泻湖相和潮坪相单元,其中沼泽亚相为煤层的最佳沉积环境。     

Parallel inexact constraint preconditioning for ill-conditioned consolidation problems

Constraint preconditioners have proved very efficient for the solution of ill-conditioned finite element (FE) coupled consolidation problems in a sequential computing environment. Their implementation on parallel computers, however, is not straightforward because of their inherent sequentiality. The present paper describes a novel parallel inexact constraint preconditioner (ParICP) for the efficient solution of linear algebraic systems arising from the FE discretization of the coupled poro-elasticity equations. The ParICP implementation is based on the use of the block factorized sparse approximate inverse incomplete Cholesky preconditioner, which is a very recent and effective development for the parallel preconditioning of symmetric positive definite matrices. The ParICP performance is experimented with in real 3D coupled consolidation problems, proving a scalable and efficient implementation of the constraint preconditioning for high-performance computing. ParICP appears to be a very robust algorithm for solving ill-conditioned large-size coupled models in a parallel computing environment.     

中国沙漠物源研究:回顾与展望

沙漠物源研究不仅在风沙地貌学上有重大的理论和实践意义,而且对联结大气粉尘排放、黄土堆积、气候系统和海洋生物地球化学循环也有重要价值。在简要回顾中国沙漠研究的基础上,梳理了中国沙漠物源研究的理论、方法和主要成果,结合国际上沉积物物源分析的趋向,指出目前世界沙漠物源的研究都是基于沉积物组份属性统计的反演模型,这种研究范式在数据获取、处理与解释方面存在缺陷,如沉积物的取样设计与测试分析、未消除"粒级依赖"对沉积物组份影响、数据未进行对数变换、忽视Dickinson图解应用的前提条件等。提出今后中国沙漠物源研究的方向:①采用正确统一的取样设计和分析方法对各沙漠的沉积物组份属性进行系统研究,建立中国沙漠沉积物组份的属性数据库;②选择若干典型沙漠,利用其周边山地详尽的地质构造、母岩和气候数据,定量构建沉积物生成的正演模型,模拟源区生成沉积物的数量、成分和结构,并用建成的沉积物组份属性数据验证和校正;③定量评估河流冲积物、冲积-湖积物、洪积-冲积物和基岩风化的残积、坡积物对中国各沙漠物源的贡献率与迁移路径,研究中国各沙漠中细颗粒物质的形成机制,对比中国沙漠与低纬度沙漠物源的形成机制;④定量研究历史和地质时间尺度沙漠-黄土-深海沉积物物源的内在联系及其驱动因素,建立陆地-大气-海洋物质循环的机理模型。     

线性回归分析在邯峰矿区2~#煤层瓦斯含量预测中的应用

以邯峰矿区53个钻孔中2#煤层资料为基础数据进行线性回归分析,首先初步确定原煤水分、挥发分、灰分、煤厚、煤层倾角、埋深、顶板(岩性、厚度)和底板(岩性、厚度)等因素与煤层瓦斯含量的相关程度,并将其定量化,经一元线性拟合筛选,删除对2#煤层瓦斯含量影响较小的灰分和底板因素,对余下因素再进行多元线性拟合,建立2#煤层瓦斯含量与2#煤的水分、挥发分、埋深、煤厚、煤层倾角和顶板(岩性、厚度)的多元线性回归方程。经检验该方程较好地表达了瓦斯含量与相关因素的定量关系,预测值与实测值接近,可作为邯峰矿区预测其它钻孔2#煤层瓦斯含量的数学模型。     

A hydrologic regression sediment-yield model for two ungaged watershed outlet stations in Africa

A hydrologic regression sediment-yield model was established to determine the relationship between water discharge and suspended sediment discharge at the Blue Nile and the Atbara River outlet stations during the flood season. The model consisted of two main submodels: (1) a suspended sediment discharge model, which was used to determine suspended sediment discharge for each basin outlet; and (2) a sediment rating model, which related water discharge and suspended sediment discharge for each outlet station. Due to the absence of suspended sediment concentration measurements at or near the outlet stations, a minimum norm solution, which is based on the minimization of the unknowns rather than the residuals, was used to determine the suspended sediment discharges at the stations. In addition, the sediment rating submodel was regressed by using an observation equations procedure. Verification analyses on the model were carried out and the mean percentage errors were found to be +12.59 and –12.39, respectively, for the Blue Nile and Atbara. The hydrologic regression model was found to be most sensitive to the relative weight matrix, moderately sensitive to the mean water discharge ratio, and slightly sensitive to the concentration variation along the River Nile's course.     

Variational inequalities for modeling flow in heterogeneous porous media with entry pressure

One of the driving forces in porous media flow is the capillary pressure. In standard models, it is given depending on the saturation. However, recent experiments have shown disagreement between measurements and numerical solutions using such simple models. Hence, we consider in this paper two extensions to standard capillary pressure relationships. Firstly, to correct the nonphysical behavior, we use a recently established saturation-dependent retardation term. Secondly, in the case of heterogeneous porous media, we apply a model with a capillary threshold pressure that controls the penetration process. Mathematically, we rewrite this model as inequality constraint at the interfaces, which allows discontinuities in the saturation and pressure. For the standard model, often finite-volume schemes resulting in a nonlinear system for the saturation are applied. To handle the enhanced model at the interfaces correctly, we apply a mortar discretization method on nonmatching meshes. Introducing the flux as a new variable allows us to solve the inequality constraint efficiently. This method can be applied to both the standard and the enhanced capillary model. As nonlinear solver, we use an active set strategy combined with a Newton method. Several numerical examples demonstrate the efficiency and flexibility of the new algorithm in 2D and 3D and show the influence of the retardation term. This work was supported in part by IRTG NUPUS.     

Muglad盆地Sufyan凹陷下白垩统AG组2段沉积特征与成因模式

下白垩统AG组2段(后简称AG2段)是苏丹Muglad盆地Sufyan凹陷的主力烃源岩和最重要的含油层系。本文基于岩石学、测井相、地震相和地震属性研究,对Sufyan凹陷下白垩统AG2段沉积体系进行精细研究,结合其沉积期构造特征、古地貌特征、古气候环境及内物源和外物源发育特征,建立了其断陷期沉积模式。Sufyan凹陷AG组断陷初期为双断模式,南北两侧高陡断层发育规模相当,随着盆地的演化和湖平面的变化,到AG2段进入高水位期,凹陷北部逐渐演变为缓坡沉积环境,南部仍然为陡坡沉积环境。凹陷北部为Babanusa凸起和中非走滑带,为Muglad盆地的盆外物源,凹陷南部以Tomat凸起与Nugara凹陷分割,发育内物源。从钻井和岩心资料来看,凹陷北部以分选较好、细粒沉积为主,在断层下降盘发现中粗粒重力流沉积,反映了北部物源长距离搬运和长期供应的特点,据此建立了断陷后期AG2段北部曲流河三角洲沉积模式,北部物源为主力物源区;南部临近凹陷断裂带钻井揭示其发育较高砂地比的细粒沉积,显示物源供给的周期性和近源性,研究认为发育辨状河三角洲,沉积模式类似于扇三角洲,总体物源供给能力较弱,具有近源水下搬运特征。南部物源和西部物源一起,构成凹陷的次物源。本文综合编制了Sufyan凹陷AG2段沉积期的两期沉积相图,建立了立体沉积模式,认为盆地优质储层应以凹陷中部近烃源岩的区带为目标进行寻找。     

北海盆地维京地堑新近纪下切谷发育主控因素及其内部充填特征*

下切谷是陆地上一种常见的侵蚀地貌,古今均十分发育,但是不同盆地、不同时代发育的下切谷,无论是在沉积充填特征还是平面展布形态等多个方面都大相径庭。为此,作者依据曼宁公式进行理论推导,再结合分析北海盆地维京地堑发育的典型下切谷实例以及前人进行的相关水槽实验结果总结认为: 地形坡度、基准面下降速率及幅度与下切谷所侵蚀地层的岩性是影响下切谷发育模式的最主要因素。这些因素共同控制了: (1)下切谷的弯度指数及宽深比,(2)下切谷平面上发育密度及其规模,(3)与下切谷伴生的陆棚三角洲、陆棚边缘三角洲沉积厚度及展布范围。海侵过程中下切谷的沉积充填类型及岩性,受平均基准面上升速率及沉积物供应量的共同控制,根据沉积充填类型及岩性的不同可将其划分为富砂型陆棚三角洲充填及富泥型河口湾充填两类。     

基于定量古地理的BQART模型深时古地势重建方法

源—汇系统分析下的古地势重建是当前沉积学领域的一个热点方向。古地势重建对于理解区域构造演化、气候变化、地表风化、物源供给及其相互关系具有重要作用。BQART沉积通量模型提供了入海流域中沉积通量与地质背景B、水流量Q、流域面积A、最大地势高度R、年均温度T之间的经验关系,逐渐发展为深时尺度下古流域地势重建的重要工具。本文在综合层序—古地理、旋回地层学、古水文比例关系和古地貌比例关系等方法的基础上,提出基于定量古地理的BQART模型深时古地势重建方法,结合沉积体积回填法可进一步获得构造抬升幅度与速率。通过该方法获得峨眉山大火成岩省内带古地势在晚二叠世期间持续缓慢上升,至晚二叠世末达200余米,累计构造抬升幅度约500~750 m。该深时古地势方法适用于气候温暖湿润、构造相对稳定条件下的中大型外流盆地,对能源盆地的勘探开发具有重要意义。     

利用瞬变电磁测量数据进行岩性、岩相划分及沉积体系识别

沙枣泉盆地内铀矿化与扇前缘亚相关系密切。根据测区岩性视电阻率的差异,开展盆地内瞬变电磁测量及相应的数据处理,建立测区视电阻率与不同岩性的对应关系,划分了不同岩性、岩相在平面上的分布范围,识别了冲积扇沉积体系,解释出扇前缘亚相的空间位置。从后期钻探验证结果看,效果非常明显。     

Efficient reservoir history matching using subspace vectors

In this paper, we describe a method of history matching in which changes to the reservoir model are constructed from a limited set of basis vectors. The purpose of this reparameterization is to reduce the cost of a Newton iteration, without altering the final estimate of model parameters and without substantially slowing the rate of convergence. The utility of a subspace method depends on several factors, including the choice and number of the subspace vectors to be used. Computational gains in efficiency result partly from a reduction in the size of the matrix system that must be solved in a Newton iteration. More important contributions, however, result from a reduction in the number of sensitivity coefficients that must be computed, reduction in the dimensions of the matrices that must be multiplied, and elimination of matrix products involving the inverse of the prior model covariance matrix. These factors affect the efficiency of each Newton iteration. Although computation of the optimal set of subspace vectors may be expensive, we show that the rate of convergence and the final results are somewhat insensitive to the choice of subspace vectors. We also show that it is desirable to start with a small number of subspace vectors and gradually increase the number at each Newton iteration until an acceptable level of data mismatch is obtained.     

河网非恒定水沙数学模型研究

建立了能够适用于多种情况的河网非恒定流水沙数学模型。推导提出了河网泥沙方程组的分级解法,提出了适用于河网计算的汊点分沙模式,分析和对比了多种分沙模式的适用性,对河网区域中可蓄水汊点也提出了相应的处理办法。采用汉江杜家台分蓄洪区的资料作了模型的验证和应用。     

The application of HVSR method in detecting sediment thickness in karst collapse area of Pearl River Delta,China

The microtremor horizontal-to-vertical spectral ratio (HVSR, or H/V spectral ratio) method is an effective tool for detecting sediment thickness. Firstly, single-station microtremor measurements were taken beside 52 boreholes located in Pearl River Delta, China. Sediment thicknesses revealed by those boreholes range from 7.9 to 39.6 m. Then, those microtremor data were analyzed using HVSR method, and peak frequency of each site was extracted. According to those peak frequencies of HVSR and corresponding sediment thicknesses, the frequency-to-thickness fitting (f–Z) equation which is suitable for Pearl River Delta area was established. And this equation was compared with several fitting equations derived from different regions by other researchers. Finally, this fitting equation was applied in a collapse area in Guangzhou city to estimate site sediment thickness. Compared with data from five boreholes in this collapse area, the consistency between estimated sediment thickness and real value was relatively good. Karst sinkholes are generally located in regions where thickness of sediments is shallow. Therefore, the method described in this paper can provide support for delineating potential collapse areas.