泥质砂岩中接触胶结泥质定量估算及对砂岩弹性的影响

在泥质砂岩的岩石物理建模中,明确泥质砂岩中泥质胶结物的接触类型及其含量对正确认识泥质的胶结作用对泥质砂岩声速的影响以及合理地建立岩石物理模型至关重要.现阶段,尚未有实验室定量估算胶结泥质的方法,导致应用胶结砂岩理论模型预测胶结砂岩地层的声速时往往由于胶结物含量被高估从而导致预测声速结果偏高.本文通过观察铸体薄片中泥质与颗粒之间的接触关系和相对分布提出了一种区分胶结泥质和分散泥质的方法:与两个或两个以上颗粒接触的连续分布的泥质为胶结泥质;与一个颗粒接触或者不与颗粒接触的泥质为分散泥质.基于这一准则,本文基于像素拾取法估算了人造泥质砂岩的胶结泥质含量,并将胶结泥质含量作为胶结砂岩模型的输入参数优化CCT模型.对比原始模型,本文方法声速误差下降了20%,预测准确度显著提高.本文方法适用于弱胶结地层的岩石物理建模,能够准确的预测声速以结合地震和测井资料识别有利储层,定量评价储层参数.     

Numerical modeling of formation damage by two-phase particulate transport processes during CO2 injection in deep heterogeneous porous media

Prediction of CO₂ injection performance in deep subsurface porous media relies on the ability of the well to maintain high flow rates of carbon dioxide during several decades typically without fracturing the host formation or damaging the well. Dynamics of solid particulate suspensions in permeable media are recognized as one major factor leading to injection well plugging in sandstones. The invading supercritical liquid-like fluid can contain exogenous fine suspensions or endogenous particles generated in situ by physical and chemical interactions or hydrodynamic release mechanisms. Suspended solids can plug the pores possibly leading to formation damage and permeability reduction in the vicinity of the injector. In this study we developed a finite volume simulator to predict the injectivity decline near CO₂ injection wells and also for production wells in the context of enhanced oil recovery. The numerical model solves a system of two coupled sets of finite volume equations corresponding to the pressure-saturation two-phase flow, and a second subsystem of solute and particle convection-diffusion equations. Particle transport equations are subject to mechanistic rate laws of colloidal, hydrodynamic release from pore surfaces, blocking in pore bodies and pore throats, and interphase particle transfer. The model was validated against available laboratory experiments at the core scale. Example results reveal that lower CO₂ residual saturation and formation porosity enhance CO₂-wet particle mobility and clogging around sinks and production wells. We conclude from more realistic simulations with heterogeneous permeability spanning several orders of magnitude that the control mode of mobilization, capture of particles, and permeability reduction processes strongly depends on the type of permeability distribution and connectivity between injection and production wells.     

横观各向同性页岩岩石物理模型建立——以龙马溪组页岩为例

在龙马溪页岩微观物性特征分析的基础上,综合利用测井解释、微观测试分析资料,建立了一种适用于龙马溪页岩的横观各向同性岩石物理模型,该模型建模过程：将各向异性SCA和DEM模型联合模拟得到的黏土和干酪根混合物作为背景介质;采用SCA模型对脆性矿物混合物进行模拟,利用各向异性DEM将脆性矿物混合物添加到背景介质;进一步将空孔隙添加到页岩基质,并利用Brown-Korringa模型进行各向异性条件下的流体替换,从而得到横观各向同性页岩岩石物理模型.通过对四川盆地A井龙马溪页岩进行岩石物理建模分析,计算了孔隙纵横比、纵横波速、各向异性系数和弹性参数,检验了模型的准确性.研究结果表明：矿物颗粒和孔隙纵横比是影响模型精度的关键参数,黏土和干酪根颗粒纵横比为0.05,图像识别获得的脆性矿物颗粒纵横比主要分布于0.45~1.0（集中分布于0.5~0.85）,横波波速反演获得的孔隙纵横比主要分布于0.1~0.3（平均值约为0.22）;模型预测和实测纵波波速之间误差为-2.40%~2.21%（平均绝对误差仅1.20%）,预测和实测横波波速之间误差为-1.93%~1.42%（平均绝对误差仅0.64%）,证实了本文模型的准确性和精度.本文模型能够准确计算页岩5个独立的刚度系数,为页岩弹性参数、声波波速、各向异性和脆性分析提供了有效手段,也为后续地球物理和工程地质参数分析提供了重要依据.     

An investigation into the effects of particle texture,water content and parallel plates' diameters on rheological behavior of fine sediment

Siltation,a phenomenon resulted from the presence of fine particles in an aqueous environment,dominated by silt and clay,is a known and common environmental issue worldwide.The accumulation of fine sediments engenders murky water with low oxygen levels,which leads to the death of aquatic life.Thus,investigating the physical and mechanical properties of fine sediment by rheological methods has expanded.Rheology is the science of deformation and flow of matter in stress.This survey investigates the rheological behavior of six samples of soil as the fine particles structure(D 63 μm) from different regions of Malaysia by using a rotational rheometer with a parallel-plate measuring(using two sizes:25 mm and 50 mm) device to explore the flow and viscoelastic properties of fine particles.The samples were examined in two rheological curve and amplitude sweep test methods to investigate the effect of water content ratio,texture,and structure of patticles on rheological properties.It was found that the content of fine sand,clay,and silt had an effect on the stiffness,structural stability,and shear behavior.Thus,the pseudoplastic and viscoelastic behaviot are respectively shown.Moreover,the amount of fine sediments present in water i.e.the concentration of these particles,has a direct effect on the rheological curve.A reduction in viscosity of samples with higher concentrations of water has been observed.As a consequence,a considerable quantity of fine sediments are distributed within the water body and remain suspended over the time.As a result,the sedimentation rate slows down.It needs to be asserted that the stotage modulus G′,loss modulus G″,and yield point can vary depending on particle type.The G′ and G″were instigated for samples(70%and 45%concentrations) that demonstrated viscoelastic characteristics using the same rotational rheometer with a parallel-plate measuring device.     

Characterization of particle size and settling velocity of cohesive sediments affected by a neutral exopolymer

In natural waters,exopolymers or extracellular polymeric substances(EPS) exuded by microorganisms interact with clay particles,resulting in the flocculation of clays and hence alteration to the properties of suspended cohesive sediments.To investigate and further understand how neutral EPS affect cohesive sediment transport and the final sediment yield,an experimental study was conducted on laboratory-prepared clay and guar gum(used as an analog for neutral EPS) suspensions to characterize EPS-induced flocculation and the settling velocity of resultant floes.Four different clays consisting of kaolinite,illite,Ca-montmorillonite,and Na-montmorillonite were studied to examine the influence of different layer charges on clay flocculation induced by neutral EPS.Floc size was determined by a laser particle size analyzer,and settling velocity estimated by analyzing the time-series floc settling images captured by an optical microscope.Results indicate that neutral EPS promote clay-EPS flocculation for all four clays with the particle/floc size significantly increased from～0.1-60μm to as large as～600μm.Clays’ layer charge has a profound influence on the clay-EPS flocculation.With the same floc size,the settling velocity of clay-EPS flocs is typically smaller than that of pure clay flocs,which is attributed to the reduced density of flocs caused by the EPS. However,for flocs of the same composition(e.g.pure clay or hybrid clay-EPS mixture),the settling velocity increases with size.The fractal dimension of these clay-EPS flocs estimated from settling velocity ranges from 1.39 to 1.47,which are smaller than that of pure clay flocs,indicating that these flocs are less compacted than the pure clay flocs.     

Nonlinear analysis of masonry structures using fiber-section line elements

For seismic analysis of unreinforced masonry (URM) buildings characterized by a box-like behavior, a widely accepted model is based on the equivalent frame idealization of walls. The equivalent frame model uses 1D elements to represent the vertical piers and horizontal spandrels which are connected by rigid nodes. The mechanical characterization of the elements is one of the crucial aspects to predict reasonably the building seismic behavior. Through the comparison with pseudo-static and dynamic experimental tests performed on two-story full-scale buildings, this paper validates the frame modeling in the OpenSees framework, which includes a fiber-section force-based beam element for the axial-flexural behavior, coupled with a cyclic shear-deformation phenomenological law.     

A three-scale model for ionic solute transport in swelling clays incorporating ion–ion correlation effects

A new three-scale model is proposed to describe the movement of ionic species of different valences in swelling clays characterized by three separate length scales (nano, micro, and macro) and two levels of porosity (nano- and micropores). At the finest (nano) scale the medium is treated as charged clay particles saturated by aqueous electrolyte solution containing monovalent and divalent ions forming the electrical double layer. A new constitutive law is constructed for the disjoining pressure based on the numerical resolution of non-local problem at the nanoscale which, in contrast to the Poisson–Boltzmann theory for point charge ions, is capable of capturing the short-range interactions between the ions due to their finite size. At the intermediate scale (microscale), the two-phase homogenized particle/electrolyte solution system is represented by swollen clay clusters (or aggregates) with the nanoscale disjoining pressure incorporated in a modified form of Terzaghi’s effective principle. At the macroscale, the electro-chemical–mechanical couplings within clay clusters is homogenized with the ion transport in the bulk fluid lying in the micro pores. The resultant macroscopic picture is governed by a three-scale model wherein ion transport takes place in the bulk solution strongly coupled with the mechanics of the clay clusters which play the role of sources/sinks of mass to the bulk fluid associated with ion adsorption/desorption in the electrical double layer at the nanoscale. Within the context of the quasi-steady version of the multiscale model, wherein the electrolyte solution in the nanopores is assumed at instantaneous thermodynamic equilibrium with the bulk fluid in the micropores, we build-up numerically the ion-adsorption isotherms along with the constitutive law of the retardation coefficients of monovalent and divalent ions. In addition, the constitutive law for the macroscopic swelling pressure is reconstructed numerically showing patterns of attractive forces between particles for bivalent ions for particular ranges of bulk concentrations. The three-scale model is applied to numerically simulate ion diffusion in a compacted clay liner underneath a sanitary landfill. Owing to the distinct constitutive behavior of the swelling pressure and partition coefficient for each ionic species, different compaction regimes and diffusion/adsorption patterns, with totally different characteristic time scales, are observed for sodium and calcium migration in the clay liner.     

Effect of particle density and inflow concentration of suspended sediment on bedload transport in rill flow

Laboratory flume experiments were carried out to evaluate the effect of particle density on bedload transport of sand‐sized particles and the effect of a suspended load of clay particles (kaolinite) on bedload transport of sand‐sized particles in rill flow conditions. Three materials in the range 400–600 µm were selected to simulate bedload transport of primary particles and aggregates: sand (2650 kg/m³), crushed brick (2450 kg/m³) and anthracite (1300–1700 kg/m³). In the two first experiments, two different methods were applied to determine bedload transport capacity of coarse particles for various conditions of flow discharge (from 2 to 15 L/min) and slope (2.2, 3 and 4%). In the third experiment, clear water was replaced with kaolinite–water mixture and bedload transport capacity of crushed brick particles was determined for a 4% slope and different concentrations of kaolinite (0, 7, 41 and 84 g/L). The results showed that bedload transport increased significantly with the decrease in particle density but the effect of particle density on transport rates was much less important than flow discharge. Velocity measurements of clear flow, flow mixed with coarse particles and coarse particles confirmed the existence of a differentiation between suspended load and bedload. In these experimental conditions, suspended load of kaolinite did not affect bedload rates of crushed brick particles. Three transport capacity formulae were tested against observed bedload rates. A calibration of the Foster formula revealed that the shear stress exponent should be greater than 1.5. The Low and the Govers unit stream power (USP) equations were then evaluated. The Low equation was preferred for the prediction of bedload rates of primary particles but it was not recommended in the case of aggregates of low density because of the limited experimental conditions applied to derive this equation. Copyright © 2008 John Wiley & Sons, Ltd.     

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

Variations of Strength,Resistivity and Thermal Parameters of Clay after High Temperature Treatment

This paper reports the variations of strength, resistivity and thermal parameters of clay after high-temperature heating. Experiments were carried out to test the physical properties of clay heated at temperatures ranging from room temperature to 800°C in a furnace. The experiment results show that below 400°C the uniaxial compressive strength and resistivity change very little. However, above 400°C, both increase rapidly. At a temperature under 400°C, the thermal conductivity and specific heat capacity decrease significantly. The thermogravimetric analysis (TG) and differential scanning calorimeter (DSC) test indicate that a series of changes occur in kaolinite at temperatures from 400 to 600°C, which is considered the primary cause of the variation of physical and mechanical properties of clay under high temperatures.     

Computational modeling of flow and sediment transport and deposition in meandering rivers

A computational modeling analysis of the flow and sediment transport, and deposition in meandering-river models was performed. The Reynolds stress transport model of the FLUENT^TM code was used for evaluating the river flow characteristics, including the mean velocity field and the Reynolds stress components. The simulation results were compared with the available experimental data of the river model and discussed. The Lagrangian tracking of individual particles was performed, and the transport and deposition of particles of various sizes in the meandering river were analyzed. Particular attention was given to the sedimentation patterns of different size particles in the river-bend model. The flow patterns in a physical river were also studied. A Froude number based scale ratio of 1:100 was used, and the flow patterns in the physical and river models are compared. The result shows that the mean-flow quantities exhibit dynamic similarity, but the turbulence parameters of the physical river are different from the model. More strikingly, the particle sedimentation features in the physical and river models do not obey the expected similarity scaling.     

Mechanism of Suspended Kaolinite Particle Clogging in Porous Media During Managed Aquifer Recharge

Managed aquifer recharge is an effective strategy for urban stormwater management. Chemical ions are normally retained in stormwater and groundwater and may accelerate clogging during the recharge process. However, the effect of water chemistry on physical clogging has not previously been investigated. In this study, we investigated the hydrogeochemical mechanism of saturated porous media clogging in a series of column experiments. The column was packed with river sand and added suspensions of kaolinite particles. Calcium chloride and sodium chloride are used as representative ions to study chemical effects. We found that an increase in ionic strength resulted in retention of kaolinite solids in the column, with a breakthrough peak of C/C₀ value of 1 to 0.2. The corresponding hydraulic conductivity decreased with increased solids clogging. Divalent cations were also found to have a greater influence on kaolinite particle clogging than monovalent cations. The enhanced hydrochemical-related clogging was caused by kaolinite solids flocculating and increasing the deposition rate coefficient by 1 to 2 times in high ionic strength conditions. Three clogging mechanisms of kaolinite solids are proposed: surface filtration, inner blocking, and attachment. This study further deepens the understanding of the mechanisms of solids clogging during aquifer recharge and demonstrates the significance of ionic strength on recharge clogging risk assessments.     

骨架导电的混合泥质砂岩通用孔隙结合电阻率模型研究

在水介质中顺序添加分散粘土颗粒、油珠、导电骨架颗粒、层状泥质,并对每一种成分进行连续积分,建立了一种适用骨架导电及含有分散粘土和层状泥质的泥质砂岩通用电阻率模型.通过对该模型的影响因素分析,发现泥质分布形式对模型计算的含水饱和度有很大影响;对应两个不同粘土颗粒电阻率或骨架颗粒电阻率的地层电导率之差,几乎与总含水饱和度无关,而对应两个不同层状泥质电阻率的地层电导率之差,随总含水饱和度增大而增大;骨架胶结指数变化对地层电导率与总含水饱和度关系曲线的影响最大,而粘土胶结指数变化对地层电导率与总含水饱和度关系曲线的影响最小;饱和度指数对地层电导率与总含水饱和度关系曲线的影响随总含水饱和度的增大而减小.通过一组骨架导电的人造岩样的试验,表明当地层水电阻率.小于颗粒电阻率时,该模型可以用于不含粘土的骨架导电的岩石.通过两组分散泥质砂岩岩样实验测量数据和一组层状泥质砂岩测井资料及实际测井资料的计算,表明本文给出的电阻率模型既适用于分散泥质砂岩地层解释又适用于层状泥质砂岩地层解释,同时,还适用于含有分散粘土和层状泥质的混合泥质砂岩地层解释.     

软弱场地地铁车站结构地震反应的数值模拟与模型试验对比分析

根据软弱场地土上地铁车站结构大型振动台模型试验结果,以软件ABAQU S为平台,采用记忆型嵌套面黏塑性动力本构模型和动塑性损伤模型,分别模拟土体和车站结构混凝土的动力特性,建立了土-地铁车站结构非线性动力相互作用二维和三维有限元分析模型,对各种试验工况下地基土-地铁车站结构体系的地震反应进行了数值模拟,并与试验结果进行了对比。结果表明：二维、三维数值模拟与振动台模型试验结果基本一致,三维模型可更好地模拟软弱场地与地铁车站结构的动力相互作用及模型结构的动力反应。数值模拟结果和振动台试验结果可相互验证其可靠性。     

Numerical study of pyroclastic surges

The dynamics of pyroclastic surges accompanied by co-ignimbrite plumes is investigated numerically. The numerical simulations are performed with a newly developed numerical model, which is based on the Navier–Stokes equations for time-dependent flows of a compressible fluid in two-dimensional Cartesian coordinates. We regard pyroclastic surges as dilute turbulent suspensions in which hot gases and fine solid particles are homogeneously mixed owing to vigorous turbulence. In other words, the gas–particle mixture is treated as a single-phase fluid whose bulk density is represented by averaging the density of each component in the numerical model. We focus on the effect of buoyancy forces generated by the thermal expansion of the air mixed into pyroclastic surges from the calm surroundings. For our purpose, the numerical model is designed to simulate relatively simple flows spreading over a horizontal flat surface. Topographic irregularity and the sedimentation process of solid particles are neglected in the present simulations. The motion of pyroclastic surges is generated by the instantaneous release of a gas–particle mixture whose density is initially larger than the ambient air density and changes nonlinearly with the temperature and concentration of suspended solid particles. Turbulent mixing is evaluated by adopting the Smagorinsky model. By employing cubic interpolated pseudo-particle (CIP) method and C-CUP method, we obtain the fine structure of flows. The behavior of calculated flows agrees fairly well with observed pyroclastic surges in nature. The current head, which remains hot and dense, keeps spreading over a horizontal surface at a speed of about 20 m s⁻¹. The spreading speed is of the order of the speed of a gravity current that excludes the influence of thermal expansion. Besides, turbulent mixing between the basal dense layer and the ambient air is enhanced by the successive development of an interfacial less-dense layer. This results in the formation of a number of buoyant plumes rising above a horizontally spreading current. Consequently, the tails of the current thickens as time progresses. A parametric study shows that the initial temperature of a gas–particle mixture should be higher than about 600 K when buoyant plumes occur owing to the thermal expansion of mixed air. The result is quantitatively interpreted by introducing a diagram that describes the relationship among the bulk density, temperature and concentration of solid particles suspended in pyroclastic surges.     

Impact of small amounts of swelling clays on the physical properties of debris‐flow‐like granular materials. Implications for the study of alpine debris flow

The effect of the small fraction of clays on the rheological behaviour of alpine debris flow is poorly understood. This is partly due to the complexity of the debris flow mineralogy and the broad particle size distribution. This study has investigated this issue by simulating an alpine debris flow with a mixture of well characterized fractions and then varying the clay fraction composition. Four samples were tested, ranging from a clay fraction made up of only kaolinite (1:1 type clay) to samples where 80 per cent of the kaolinite is replaced by bedeillite (a 2:1 type clay similar to smectite). Changing the content of 2:1 type clay has a strong influence on the behaviour of the whole material, despite its low weight fraction of around 2 per cent. The tests carried out on these reconstituted materials were compared with the results obtained for natural debris flow materials and showed some common trends: in particular, the rheological parameters for materials with and without 2:1 clays with respect to yield stress as a function of solid content. Copyright © 2006 John Wiley & Sons, Ltd.     

Sedimentation model and analysis for differential settling of two - particle - size suspensions in the Stokes region

A model is developed for predicting the settling velocity in suspensions of particles of two different sizes based on experimental data for the settling rate of two-size suspensions in various liquids using particles of equal density. In these experiments, the retarding effect of the smaller particles on the settling velocities of the larger ones is taken into account. The model considers Steinour’s fundamental equation and assumes a fixed arrangement of particles and constant velocity in a single-size susp...     

Response of suspended sediment particle size distributions to changes in water chemistry at an Andean mountain stream confluence receiving arsenic rich acid drainage

Acid drainage is an important water quality issue in Andean watersheds, affecting the sustainability of urban, agricultural, and industrial activities. Mixing zones receiving acid drainage are critical sites where changes in pH and chemical environment promote the formation and dissolution of iron and aluminum oxy/hydroxides. These particles can significantly change the speciation of toxic metals and metalloids throughout drainage networks via sorption, desorption, and settling processes. However, little is known about the behavior of particle size distributions (PSDs) in streams affected by acid drainage and their relationship to metal speciation. This work studied: (a) the PSDs for a wide range of mixing ratios found at a fluvial confluence affected by acid drainage, and (b) the response of PSDs and arsenic speciation to environmental changes found when the particles approach complete mixing conditions. The confluence between the Azufre River (pH ~ 2, high concentration of dissolved metals) and Caracarani River (pH = 8.6, low concentration of dissolved metals) was used as a representative model for study. Field measurements show a bimodal PSD with modal diameters of ~50 and 300 μm. At shorter distances from the junction, the smaller modes with smaller particle volumes were dominant across the stream cross‐sections. A systematic shift towards larger particle sizes and larger particle volumes occurred downstream. The analysis of laboratory PSDs for Azufre/Caracarani mixing ratios between 0.01 and 0.5 (pH from 6.2 to 2.3) showed a bimodal trend with ~15 and 50 μm characteristic diameters; larger particles formed at pH>4. When particle suspensions were transferred in laboratory experiments from very low pH to full mixing conditions (pH ~ 2.8 and mixing ratio ~ 0.25) particle sizes varied, and the dissolved arsenic concentration decreased. The observed reaction kinetics were slow compared to the time scale of advective transport, creating opportunities for engineered controls for arsenic. This work contributes to a better understanding of the chemical‐hydrodynamic interactions in watersheds affected by mining, and identifying opportunities to improve water quality at points of use.     

Modeling oil spill trajectory in coastal waters based on fractional Brownian motion

This paper proposes a numerical method to simulate oil spill trajectories, which are affected by the combination of advection, turbulent diffusion and mechanical spreading process, based on a particle tracking algorithm. Recent studies have shown that the trajectories of drifters on the ocean surface have a fractal structure that is far from being described using ordinary Brownian motion. Thus, in modeling the diffusion process, a discrete method has been employed for the generation of fractional Brownian motion (fBm) to illustrate superdiffusive transport. The algorithm is implemented to predict oil slick trajectories following the “Arteaga” oil spill accident that occurred near the Dalian coastal region in 2005. When compared with the observed data and the results of traditional diffusion modeling, the numerical results based on the fBm model are encouraging.