岩体裂隙结构模型及其渗流规律研究

通过对岩体裂隙结构特征的分析,应用随机原理模拟生成岩体裂隙的结构及形态。最后结合渗流过程分析,得出工程岩体中裂隙赋存状态对渗流过程的影响。     

三维应力作用下岩体单个裂隙的渗流特性分析

天然岩体中存在着大量的孔隙和裂隙,这些缺陷不仅改变了岩体的力学性质,也严重影响了岩体的渗流特性。在对现有裂隙岩体渗流特性研究成果进行分析的基础上,讨论了岩体单个裂隙的力学性质和渗流对单个裂隙岩体产生的力学作用,研究了岩体单裂隙渗透系数与岩体三维应力的关系,考虑了裂隙粗糙度对渗流的影响以及不同方向应力对渗透性影响的差异,分析了单个裂隙岩体在三向应力作用下的渗流特征,得出了裂隙所受三维应力与渗透系数关系式,认为垂直于裂隙面的应力对岩体渗透性起主导作用,岩体渗透系数随垂直裂隙面应力的增加而迅速减小。通过与渗流规律试验结果对比分析,证明了所得单个裂隙岩体渗透系数表达式的正确性。     

裂隙岩体渗流与应力耦合的试验研究

通过对较大尺寸的裂隙岩体试块进行不同侧压力和加载条件下的渗流试验结果研究，分析了裂隙岩体渗流与应力的耦合机理，得出了不同应力条件下裂隙岩体渗流量与应力成四次方的关系。并且得出并非压应力都引起裂隙岩体的渗流量减小，当裂隙岩体受平行于裂隙面方向的单向压应力时，渗流量随着压应力的增加而增加。     

低温裂隙岩体水-热耦合模型研究及数值分析

裂隙渗流会引起裂隙周围岩体中的温度场变化,在低温岩体中其影响更为明显;此外,裂隙水与周围低温岩石介质发生热交换会引起裂隙中的水冰相变过程发生,而裂隙水冻结将阻碍裂隙渗流,引起裂隙渗流场的变化。因此,低温下的裂隙岩体水-热相互作用是一个强耦合过程。考虑裂隙中的水冰相变过程和渗流作用,建立了低温冻结条件下裂隙岩体水-热耦合模型;以冻结法施工为例,考察了低温冻结过程中裂隙水渗流对裂隙冻结交圈的影响。研究结果表明：由于裂隙渗流的存在,距裂隙较远处岩石先冻结,裂隙冻结所需时间远大于周围岩石;裂隙宽度和裂隙水压力差都会影响冻结交圈时间,裂隙越宽、水压力差越大,裂隙冻结需要时间越长;随着冻结时间的推进,裂隙水渗流速度逐渐降低,当裂隙冻结后裂隙渗流停止。最后通过构建随机裂隙网络模型,利用所建立的水-热耦合模型考察了裂隙网络渗流对冻结交圈的影响,说明了在冻结法施工中考虑裂隙的重要性。     

降雨作用下路基裂隙渗流分析

为探讨路基裂隙在降雨作用下渗流特征及其对路基含水率影响作用,分析了路基裂隙产生的原因及土质裂隙特征,研究了土质裂隙渗流特点。根据降雨IEM及FDEM模型,推导了降雨作用下路基裂隙渗流公式。利用Hydrus渗流软件,分别模拟了路基中的裂隙及路堑边坡裂隙渗流下含水率分布状况特征。模拟结果表明,降雨过程中路基裂隙渗流能够较显著的改变局部路基含水率,在工程中值得重视。     

渗透系数的各向异性对裂隙水渗流潜水面的影响

天然岩体是多相的不连续介质，在其中分布有大量的节理裂隙，这些裂隙的渗透能力要比岩石高几个数量级，构成了地下水流的主要通道，并控制着地下水渗流的运动规律．由于裂隙分布的方向性，导致裂隙岩体的渗流具有明显的各向异性．本文便在规则裂隙岩体模型的基础上，改变每组裂隙的渗透系数，利用有限元计算出裂隙水渗流的潜水面，分析了裂隙岩体渗透系数的各向异性对裂隙水渗流潜水面的影响．     

压剪作用下裂隙扩展过程渗流与应力耦合模型研究

裂隙的剪胀特性及扩展演化规律对岩体的渗流特性具有重要影响。为揭示裂隙剪胀及扩展演化对岩体渗流的影响,基于残余强度提出了一种能较好描述岩石硬化-软化特性的全剪切本构关系;结合剪切变形与裂隙开度的关系,利用最小势能原理和立方定律,建立了压剪作用下考虑裂隙剪胀特性的渗流应力耦合模型;假定压剪作用下裂隙发生Ⅰ型扩展,提出了伴有翼型裂纹的渗流模型,该模型不仅考虑了岩石的剪胀特性,更反映了裂隙扩展过程渗流的演化规律。对不同裂隙粗糙度的剪切应力-位移曲线进行分析,全剪切本构模型表现出较高的拟合精度。在剪切应力-位移关系基础上,通过剪切渗流试验数据对压剪作用下渗流模型进行验证,结果表明,该模型能较好地描述岩体剪胀阶段渗透性演化规律。利用等效裂隙简化裂隙网络,并通过试验数据进行验证,证实了裂隙扩展过程渗流-应力耦合模型的准确性与适用性。     

三峡围堰渗流状态的边界元计算分析

本文对三峡围堰的渗流稳定性进行了边界元计算分析，共计算了１８个方案，分析了防渗墙插入基岩深度、风化砂密实性、防渗墙产生裂隙以及基岩渗透系数对堰体渗透性能的影响，得出了若干对围堰工程设计和施工十分有意义的结论。     

三峡船闸高边坡降雨入渗的三维数值仿真

基于反映裂隙岩体渗透张量各向异性特性的广义达西定律与相对透水率等裂隙岩体渗流基本理论,应用三维非线性有限元方法,对长江三峡工程永久船闸高陡边坡在连续降雨条件下裂隙岩体饱和-非饱和、非恒定渗流扩散过程进行了三维数值仿真分析,得出了连续降雨条件下船闸高陡边坡裂隙岩体渗流场的时空分布规律及典型降雨时段闸室混凝土衬砌上的最不利外水压力荷载,对暴雨入渗情况下闸室边墙稳定性与安全性进行了评价。分析结果表明,闸室混凝土衬砌具有较高的安全度;船闸高陡边坡内排水系统效果显著,即使遭遇长期强降雨,闸室衬砌上的外水压力荷载受降雨的影响也不大。     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

应用随裂缝倾角变化的孔隙模型改进储层饱和度计算

复杂的孔隙结构及多变的孔隙类型使得裂缝性储层成为测井评价的难点,通过假定不同孔隙类型间的导电方式来建立合理的孔隙模型,进而结合阿尔奇方程实现饱和度计算是目前裂缝性储层定量评价的主要方法和途径。在经典的双孔隙模型中求取岩石系统的总孔隙指数m时,未考虑裂缝与电流线的夹角变化带来的影响。通过在裂缝性储层中同时考虑裂缝曲折度、裂缝与电流线夹角双重因素,构建了改进双重孔隙模型,并结合电成像资料提供的倾角信息进行了实例处理。结果表明,改进模型的计算值与取心分析有着更好的一致性;若不考虑裂缝与电流线夹角等因素,将带来321%的相对误差。     

XFEM-based cohesive zone approach for modeling near-wellbore hydraulic fracture complexity

Near-wellbore fracture tortuosity has important impacts on the productivity of fractured oil and gas wells and the injectivity of CO₂ or solids disposal injectors. Previous models for simulating near-wellbore fracture tortuosity usually assume fracture growth in linear-elastic media, without considering the effects of porous features of the rock. In this paper, a 2D fully coupled model is developed to simulate near-wellbore fracturing using the XFEM-based cohesive segment method. The model takes into account a variety of crucial physical aspects, including fracture extension and turning, fluid flow in the fracture, fluid leak-off through wellbore wall and fracture surfaces, pore fluid flow, and rock deformation. The proposed model was verified against two sets of published experimental results. Numerical examples were carried out to investigate the effects of various parameters on near-wellbore fracture trajectory, injection pressure, and fracture width. Results show that near-wellbore fracture behaviors are not only dependent on rock elastic properties and field stresses, but also greatly influenced by porous properties of the rock, such as permeability and leak-off coefficient. Some field implications were provided based on the simulation results. By overcoming some limitations of the previous models, the proposed model predicts more realistic fracture evolution in the near-wellbore region and provides an attractive tool for design and evaluation of many field operations, for which near-wellbore fracture behaviors play an important role on their successes.

     

Modified governing equation and numerical simulation of seepage flow in a single fracture with three-dimensional roughness

Roughness and tortuosity influence groundwater flow through a fracture. Steady flow through a single fracture can be described primitively by the well-known Cubic Law and Reynolds equation with the assumption that the fracture is made of smooth parallel plates. However, ignoring the roughness and tortuosity of the fracture will lead to inaccurate estimations of the flow rate. To obtain a more accurate flow rate through a rough fracture, this paper has derived a modified governing equation, taking into account the three-dimensional effect of the roughness. The equation modifies the Reynolds equation by adding correction coefficients to the terms of the flow rates, which are relative to the roughness angles in both the longitudinal and transverse directions. Experiments of steady seepage flow through sawtooth fractures were conducted. The accuracy of the modified equation has been verified by comparing the experimental data and the theoretical computational data. Furthermore, three-dimensional numerical models were established to simulate the steady flow in rough fractures with the triangular, sinusoidal surfaces and the typical joint roughness coefficient (JRC) profiles. The simulation results were compared with the calculation results of the modified equation and the current equations. The comparison indicates that the flow rate calculated by the modified equation is the closest to the numerical result.     

基于连续损伤的水平井射孔-近井筒三维破裂模拟

射孔作为井筒与储层之间的液流通道,是水力压裂过程中的重要可控性参数。为研究水平井射孔-近井筒破裂机制,采用岩层变形-流体渗流方程描述应力状态变化,应用连续损伤破裂单元表征三维破裂位置与形态演化,并开发有限元求解程序模拟分析了射孔对水平井初始破裂压力、破裂位置及近井筒裂缝复杂性的调控作用。通过与解析模型及射孔压裂物理模型试验结果对比,验证了模型及有限元程序的有效性;水平井破裂压力数值分析结果与现场测试数据吻合较好。研究表明：射孔可调控水平井破裂压力与初始破裂位置,同时对近井筒区域裂缝扩展形态影响显著。通过优化射孔参数可以引导初始破裂向最优破裂面扩展、有效降低破裂压力,减小由于螺旋射孔空间排布引起的水平井近井筒裂缝迂曲与复杂程度,提高致密油气藏压裂改造效果。     

考虑非达西效应的酸蚀裂缝流场数值模拟

酸压改造后的酸蚀裂缝是储层流体的主要流动通道,研究流体在酸蚀裂缝中的流动规律是十分必要的。本文介绍了目前裂隙流动研究中的4种常用控制方程及各自的适用条件,确定了Navier-Stokes方程为酸蚀裂缝流场数值模拟的控制方程。为了分析酸蚀裂缝流场分布规律以及非达西效应对裂缝导流能力的影响,利用逆向工程技术对两种不同刻蚀形态的酸蚀岩样壁面进行了实体重构,并采用有限元数值模拟方法进行了不同流量下的酸蚀裂缝流动试验。结果表明：酸刻蚀形态和裂缝接触关系对流场分布的影响较大。隙宽分布平滑的均匀刻蚀裂缝,其流态稳定、曲折度低,但由于裂缝开度较窄且比表面较高,边界层效应明显,流体流动时产生的黏滞阻力较高;沟槽裂缝曲面粗糙,隙宽分布复杂,流体流动时流态不稳定,曲折度高,在高流量下会产生明显的涡流流动,增大惯性阻力。裂缝通道存在缩颈现象会引起流体绕流和多次加减速,产生额外的压力损耗。随着模拟流量的增大,压降与流量会逐渐偏离线性关系而呈现出非达西流动现象。酸蚀裂缝壁面越粗糙,产生非达西效应时的临界流量和临界雷诺数就越小,在相同流量下非达西效应就越强烈,导流能力下降速率更快。     

Effect of macropore tortuosity and morphology on preferential flow through saturated soil: A Lattice Boltzmann study

Unlike micropores where water moves upward or downward based on hydraulic gradient, in macropores, water flows predominantly downward due to the gravity. Therefore, models based on capillary flow are not capable of simulating macropore flow. There are attempts to model the macropore flow using two domains, one for capillary flow and another one for macropores. These models use Richard’s equation for capillary flow and Poiseuille’s law for macropores in which the macropore is approximated to be cylindrical or planar. This study quantifies the magnitudes of the errors induced by this assumption. Influence of macropore shapes and tortuosity was quantified by using a 3D Lattice Boltzmann model, which is capable of simulating fluid flow in micropores as well as macropores of cracked clays. Artificial macropores of constant sectional area and volume, but different shapes were generated in 3D and the influence of macropore shapes, shape related parameters, and tortuosity were systematically investigated. Macropore flow rate decreases with different shapes compared to cylindrical macropores and increase in aspect ratio of sectional shape leads to decrease in macropore flow rate. The maximum effect of bends/turnings along the tortuous macropore was about 25% on overall decrease of flow rate due to tortuosity. However, more detailed study is required on the influence of bends on macropore flow rate. The macropore flow rate reduces by about 70% for tortuosity of 1.41. A prediction equation is verified to predict the flow rate of different shapes and tortuous macropores based on straight cylindrical macropore using aspect ratio and tortuosity factor.     

Pore scale study of permeability and tortuosity for flow through particulate media using Lattice Boltzmann method

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic permeability of granular media. The obtained results show that particles' size diameter has a two‐fold effect on the coefficient of permeability: one is through specific surface and the other is by tortuosity factor. Numerical study also reveals that tortuosity of granular soils decreases almost linearly with increasing the porosity. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing flow paths in a karst aquifer based on multiple dye tracing tests using stochastic simulation and the MODFLOW-CFP code

Karst systems show high spatial variability of hydraulic parameters over small distances and this makes their modeling a difficult task with several uncertainties. Interconnections of fractures have a major role on the transport of groundwater, but many of the stochastic methods in use do not have the capability to reproduce these complex structures. A methodology is presented for the quantification of tortuosity using the single normal equation simulation (SNESIM) algorithm and a groundwater flow model. A training image was produced based on the statistical parameters of fractures and then used in the simulation process. The SNESIM algorithm was used to generate 75 realizations of the four classes of fractures in a karst aquifer in Iran. The results from six dye tracing tests were used to assign hydraulic conductivity values to each class of fractures. In the next step, the MODFLOW-CFP and MODPATH codes were consecutively implemented to compute the groundwater flow paths. The 9,000 flow paths obtained from the MODPATH code were further analyzed to calculate the tortuosity factor. Finally, the hydraulic conductivity values calculated from the dye tracing experiments were refined using the actual flow paths of groundwater. The key outcomes of this research are: (1) a methodology for the quantification of tortuosity; (2) hydraulic conductivities, that are incorrectly estimated (biased low) with empirical equations that assume Darcian (laminar) flow with parallel rather than tortuous streamlines; and (3) an understanding of the scale-dependence and non-normal distributions of tortuosity.     

考虑接触面积影响的粗糙节理渗流分析

节理中接触面积的存在使得渗流曲折效应更加明显,将对渗流产生显著影响。为了研究节理接触面积对节理渗流的影响,建立节理渗流概念模型,推导得到了节理渗流的接触面积影响系数;然后在节理渗流计算公式中引入该影响系数,得到了考虑接触面积影响的渗流计算公式。同时,由于接触面积的存在使得渗流曲折因子的计算更为复杂,在获得节理隙宽分布的基础上提出了计算节理渗流曲折因子的五点比较法。然后以人工大理岩节理试件为研究对象,在获得其三维表面形貌的基础上,一方面对其进行不同接触状态下5级法向应力作用下的渗流试验;另一方面,在根据节理试件的三维表面形貌数据计算节理的渗流平均曲折因子和接触面积比的基础上,分别采用考虑接触面积影响的渗流计算公式和Zimmerman计算公式对节理试件在不同法向应力作用下的渗流体积流量进行计算,并将计算结果与实测值进行比较。结果表明：考虑接触面积影响的渗流计算公式计算结果与实测值吻合较好,从而验证了考虑接触面积影响的渗流计算公式的正确性,而Zimmerman计算公式高估了节理的渗流流量。     

土壤干缩裂缝网络几何特征

通过室内试验研究了土壤干缩裂缝发育规律,采用数字图像处理技术与AutoCAD测量功能,分析了不同含水率下土壤裂缝网络几何形态特征.结果表明:土壤含水率达到43%时出现裂缝,裂缝骨架长度、周长和面积分别在含水率达到23%、23%、7%时趋于稳定,裂缝平均宽度的发育经历了波动增长、线性增长和稳定3个阶段,并在含水率达到7%时保持不变;裂缝条数与节点数之比随含水率降低而减小,并最终保持在1.67,此结果与收缩块区分散度和裂缝交叉角度分布规律共同表明,裂缝网络形态呈方形正交网络分布趋势,且处于方型网络与“T”型正交网络之间;裂缝平均弯曲度随含水率减小,由1.06减小至1.02并保持恒定,揭示了裂缝整体弯曲程度与不同时期裂缝弯曲的变化规律.