现代油气渗流力学体系及其发展趋势

油气渗流是油气藏开采的科学核心问题,以连续介质假设和达西方程为基础的传统渗流力学在常规油气资源开发中发挥了重要作用.近年来,非常规油气藏成为油气行业勘探开发的主要阵地,其开发理论和技术尤其是渗流问题也成为学术界研究的热点问题.非常规油气藏其岩石多孔介质具有明显的多尺度特征,尺度差异达到6个数量级之多,而且采用大规模的水平井分段压裂开发模式,储层应力强烈作用.因此,传统的油气渗流理论已无法准确描述非常规油气藏的流动特征.实质上,非常规油气资源的开发过程是一个典型的多场作用下的多流动模式的多相流体(油气水)在多尺度多孔介质的流动过程.为此,本文提出了"多场作用下的多流动模式的多相流体在多尺度多孔介质中流动动力学体系"的现代渗流力学体系的概念,并从纳微尺度油气流动模拟、流动模拟的多尺度升级、非常规油气藏的宏观流动模拟、大尺度缝洞碳酸盐岩油藏的流动模拟以及油气渗流物理模拟等方面系统阐述了其研究现状及发展趋势.     

考虑微裂缝的数字岩心多点统计学构建方法

页岩的微观孔隙结构决定了其宏观性质和流动规律,其微观孔隙结构呈现出典型多尺度特点,孔隙尺寸从纳米到微米不等,同时还伴有天然发育的微裂缝.物理实验方法受限于分辨率和样品尺寸之间的矛盾,无法获取能够同时含有微纳米孔隙信息和微裂缝特点的表征单元体,因此数值重建算法显得尤为重要.传统的重建算法只能构建单一孔隙类型的数字岩心,针对页岩多尺度特点,本文提出了一种基于模式的多点统计学方法,结合训练图像和数据模板捕获模式库和条件概率函数,利用模式法的插值方式对待模拟空间进行构建.通过自相关函数、孔隙形态拓扑结构及渗透率方面,对本文提出的数值重建算法进行评价.在不同类型模拟算例中,本文算法呈现较好的可靠性.与传统重建算法相比,本文算法能够构建页岩多相多尺度数字岩心.     

离散位错动力学算法及其在材料塑性行为模拟中的应用

晶体材料的塑性变形由位错的运动演化而引起.离散位错动力学(discrete dislocation dynamics, DDD)通过直接模拟大量位错的演化而研究材料的塑性变形,因此能够揭示材料微结构-位错微结构-塑性力学行为之间内在的物理关联,并能够自然而然地捕捉塑性变形微米/亚微米特征尺度下本征的尺度效应.它所能模拟的尺度介于微观分子动力学模拟和宏观有限元模拟之间,在多尺度算法中起到承上启下的作用.本文首先系统地发展、完善和丰富了离散位错动力学-有限元(finite element method, FEM)叠加算法、DDD-FEM直接耦合算法(discrete-continuous method, DCM)以及离散位错动力学-扩展有限元(extended finite element method, XFEM)耦合算法等框架体系.在此基础上,利用这些方法对单晶镍基高温合金的塑性变形机理、晶体材料的断裂和损伤变形行为以及塑性行为的微尺度和微结构效应3个方面开展了系统的研究.所得模拟结果指导了基于微结构和位错机制的单晶镍基高温合金晶体塑性本构模型的建立,丰富和加深了人们对材料强化、循环塑性、断裂、损伤、尺度效应和微结构效应的认识.此外,离散位错动力学可进一步应用于诸如高温、高压、高应变率、化学腐蚀环境、高辐照等极端条件下晶体材料塑性行为的研究,是材料力学行为多尺度模拟研究中的重要一环.     

旋转交错网格有限差分及其完全匹配层吸收边界条件

将完全匹配层吸收边界条件引入到旋转网格有限差分中, 用以解决在非均匀弹性和孔隙弹性介质情况下数值模拟中的吸收边界问题, 同时, 将旋转交错网格有限差分法用于数值求解等效弹性介质、孔隙弹性介质和各向异性弹性介质的波动方程. 与普通的交错网格有限差分方法相比, 旋转交错网格有限差分的好处在于不同的物理量只位于两个不同的位置: 应力和应变(或质点速度和位移)位于离散单元的中心, 而质点速度或位移(或应力和应变)位于单元的顶点. 经过这样的处理后, 弹性常数就不再需要进行平均(模拟非均匀介质时)或内插(模拟各向异性介质时), 因为此时所有的弹性模量都位于相同的位置, 且与应力或应变的位置相对应. 为了验证新算法, 对相同的模型采用了不同的算法计算和比 较. 研究结果表明, 旋转交错网格有限差分算法和普通交错网格有限差分算法的结果吻合很好. 不仅如此, 新算法能很方便地处理声阻抗差别较大的非均匀介质, 特别是在模拟充有液体或气体的裂缝介质时更具优势. 另一方面, 应用完全匹配层吸收边界条件可以大大减少人工界面产生的反射波. 当边界层的厚度超过半个波长时, 在人工界面几乎无反射波产生. 理论和数值模拟结果表明, 旋转网格中的完全匹配层吸收条件与普通交错网格中的吸收效果和处理方法几乎相同. 此外, 建立了在等效弹性、孔隙弹性和各向异性弹性介质中的完全匹配层吸收条件的速度-应力差分方程系统.     

有限容积法与格子Boltzmann方法耦合模拟传热流动问题

自然界和工程领域中的许多物理现象的发生通常涵盖几个数量级的几何空间及时间范围, 我们将其统称为多尺度物理现象. 在模拟多尺度问题时, 如果仅采用宏观方法, 则会存在一些不足, 如无法预知微小部分的细节以及引入复杂的经验关联式; 如果仅采用介观/微观方法, 则需要消耗大量的计算资源. 构造宏观-介观、宏观-微观、宏观-介观-微观等多种层次上方法的耦合体系, 可以在很大程度上克服这些不足. 构造了宏观有限容积法(FVM)与介观格子Boltzmann方法(LBM)的耦合模型(CFVLBM), 给出了由宏观物理量重构密度分布函数和温度分布函数的两个重构算子, 解决了LBM与宏观方法耦合的关键难题. 选取二维、三维典型传热流动问题对耦合模型进行了考核, 计算结果同基准解符合得很好. 最后将CFVLBM应用于计算多孔介质内的复杂流动问题. 研究表明, 基于文中重构算子的CFVLBM可以准确有效地应用于模拟传热流动问题.     

页岩水力压裂的关键力学问题

页岩气的开采成为我国绿色能源开发的新领域.尽管北美页岩气革命取得了成功,水力压裂是成功的开采方式,目前采收率仅为5%~15%,问题出在哪里呢?由此给力学家提出了巨大的挑战和机遇.本文针对页岩水力压裂的关键力学问题,阐述理论、计算和实验的研究进展和技术难点,主要内容有:页岩人工裂缝扩展的大型物理实验模拟平台;考虑时间相关性的各向异性本构模型;页岩起裂、分叉及多裂缝相互作用的断裂力学准则和模拟方法;裂缝簇稳定性扩展的力学条件和创造缝网的多尺度有限元模型;耦合断裂力学和流场压力的裂缝网扩展数值模拟方法.     

泥石流动力学模型与数值模拟

泥石流是介于崩塌、滑坡等块体运动和山洪水流之间的一种物理过程,既具有土体的结构性,又具有水体的流动性。从物质组成上看,泥石流是一种由水、岩土体和气体组成的多相介质,具有多种内部结构。其中的固相组分颗粒形状极其不规则,尺度跨越范围较大。浆体与固相颗粒以及颗粒间的相互作用非常复杂。因此,完全考虑各种因素建立全描述的泥石流动力学模型比较困难。从描述组成物质和运动的观点来看,目前泥石流的动力学模型可划分为连续介质、离散介质和混合介质模型。基于不同动力学模型的泥石流运动数值模拟,可广泛应用于流量过程反演、危险范围预测、风险评估、防治工程评估等方面;但泥石流本身启动和产汇流机理涉及多门学科,是需要进一步研究的重要课题。同时,泥石流形成过程和运动过程的数值模拟存在着时空尺度上的差异,如何实现两者的耦合求解仍需深入探讨。     

磁场-渗流场耦合作用下的铁磁流体多孔介质流动数值模拟

铁磁流体作为一种新型功能材料,由于同时具备磁性及流体性质被广泛应用于各个技术领域.本文从数值模拟理论上研究了利用磁场控制的铁磁流体驱油问题:针对外磁场作用下的铁磁流体多孔介质流动物理过程,将磁力项引入铁磁流体相运动方程从而耦合磁场-渗流场给出其流动方程,并采用全隐式有限体积法形成相应数值离散格式;通过二维均质填砂平板流动实验验证了模型和算法的正确性,在此基础上,对比了平面及纵向非均质多孔介质中注水及注铁磁流体的驱油效果,计算结果表明通过控制铁磁流体的驱替路径可以提高驱替波及范围继而提高采收率,为提高原油采收率提供了新的思路和方法.     

有限容积法离散方程截断误差本质及其来源

本文从有限容积法的原理出发,指出有限容积法离散方程截断误差的本质为离散表达式与守恒型积分方程的近似程度.基于该思想,推导了有限容积法离散方程截断误差表达式,该表达式与有限差分法截差表达式有显著的不同;并指出有限容积法中涉及的截断误差主要包含3个层次,分别为由界面上值的近似、界面通量的近似和离散方程的近似带来的截断误差.有限容积法界面通量截断误差包含法向和切向两部分,有限差分法不存在界面未知值近似和切向截断误差,因而即使对于常物性二维/三维问题,两者截差表达式也不相同.与有限差分法相比,有限容积法误差来源较多,单纯提高待求变量的离散精度并不能有效地提高离散方程的整体精度.此外,通过理论推导和数值试验证明了对于规则区域正交非均分网格,计算区域采用内节点法离散的计算精度总体上优于外节点法.     

格子Boltzmann方法的工程热物理应用

格子Boltzmann方法在过去的20年里已经发展成为一种有效的数值模拟方法, 它是介于微观分子动力学方法和基于连续介质假设的宏观方法之间的一种介观方法. 该方法与传统的流体模拟方法不同, 它基于分子动理论, 通过跟踪粒子分布函数的输运而后对分布函数求矩来获得宏观平均特性. 基于本课题成员的一些工作, 本文对格子Boltzmann方法在工程热物理领域中研究进展做了简要回顾, 包括: (1) 格子Boltzmann方法的发展简介; (2) 沿着Boltzmann方程 —— Maxwell分布 —— Boltzmann-BGK方程 —— 格子Boltzmann-BGK方程这一主线对格子Boltzmann方法的基础理论和基本模型做了简要介绍; (3) 在格子Boltzmann方法的模型发展方面, 介绍了用于理想气体可压缩流动与传热的耦合双分布函数模型和用于非平衡态气体流动模拟的格子Boltzmann模型. 在格子Boltzmann方法的边界处理方面, 介绍了反弹与镜面漫反射格式及质量修正格式. 在格子Boltzmann方法的计算格式方面, 介绍了显式-隐式有限差分格式, 该格式能够有效地提高计算效率; (4) 在格子Boltzmann方法的应用方面, 着重对其在交变流动、可压缩流动、多孔介质流动、微尺度气体流动及热声中的应用做了详细介绍; 最后指出了格子Boltzmann方法在工程热物理领域中需要进一步研究的内容.     

磁性液体数值模拟研究进展

磁性液体(magneticfluid)是一种力学性能受磁场调控的磁流变智能材料,在医学与工业生产中有广泛应用.近年来,随着计算机性能的发展,数值模拟日益成为研究磁性液体力学行为细观机理的重要方法.本文综述并评价了磁性液体理论与数值模拟领域的最新研究进展,介绍了磁性液体在剪切、挤压和阀模式3种工作模式下的力学模型,对磁性液体的现有模拟方法进行了总结,讨论了磁性液体数值模拟的研究现状,并展望了数值模拟的发展前景.磁性液体的数值模拟亟待开展以下三方面研究:(1)建立涵盖多种微观相互作用的精细理论模型,研究颗粒表面包覆、添加剂等非磁性成分对磁流变效应的影响;(2)将不同数值模拟方法相结合,建立磁性液体的多尺度模型,进一步提高模拟精度,利用机器学习协调不同尺度的数值模拟,压缩计算量,已成为一种可行思路;(3)将力学模型与电、磁学模型相结合,发展多尺度、多物理场耦合的数值模拟方法,模拟磁性液体其他物理性能.最终为高性能磁性液体的研制及其应用研究提供技术和理论支撑.     

全球有限区数值预报模式动力框架的试验验证

全球有限区通用数值预报模式是21世纪数值预报模式的发展趋势, 为此“十五”期间中国气象局发展了一个GRAPES(global/regional assimilation and prediction system)全球有限区多尺度非静力数值预报模式. 为了验证GRAPES模式方案设计的科学性、程序编写的正确性, 针对模式的特点设计了一套完整的理想试验方案来对GRAPES动力模式的性能进行检验, 包括用于检验中尺度模式的密度流和三维地形波试验; 用于检验近极地半拉格朗日上游点和极区离散方案正确性的越极地气流试验等. 密度流试验显示模式框架对精细尺度的非线性流及其瞬变特征有较强的模拟能力; 三维地形波试验说明GRAPES模式动力框架能够较好地模拟气流过山时重力内波的水平与垂直传播. 而越极地气流试验验证了模式对极地半拉格朗日上游点及近极地的离散处理是正确的. 实际个例试验表明, GRAPES全球模式对影响中国夏季的大尺度环流形势有较好的预报能力.     

利用LBM-FVM-CA耦合方法模拟管表面上的颗粒沉积与脱离过程

发展了一个包含积灰形状演化、颗粒沉积和脱离机理的数值模型和计算方法,模拟了单排管表面的积灰过程.首先建立了流动的格子Boltzmann方法-有限容积方法(LBM-FVM)耦合的计算模型,给出了通过宏观参数构建多松弛LBM分布函数的重构算子;并结合元胞自动机模型、能量平衡模型以及力和力矩分析,模拟了颗粒的运动、沉积和碰撞过程;其次,针对模拟的时间步长相对于实际积灰时间较短的困难,提出了用于将模拟时间换算到实际时间的比例;然后,具体模拟和分析了不同直径颗粒在不同入口速度下的积灰过程.结果显示,积灰面积随时间呈指数增长并趋近于一个平衡值;在颗粒质量分数一定时,存在一个积灰速率较高的速度范围;颗粒脱离在迎风面和管子正后方较为严重,但在背风面侧面相对缓和;积灰层在整个背风面生长,在迎风面形成锥形,改变了流动并阻止了颗粒的进一步沉积.     

基于格子Boltzmann方法的复杂多孔介质内双扩散效应的对流传热传质机理研究

基于TD2G9不可压格子Boltzmnann模型, 通过引入第3个分布函数表征浓度场的演化, 并在标准演化方程后附加源项, 构造了用于模拟多孔介质内在多物理场(浓度场, 温度场)下交叉耦合效应的自然对流传热传质格子Boltzmann模型. 基于非平衡态不可逆热力学的基本原理, 引入Boussinesq假设, 在考虑了耦合扩散效应的基础上建立了可用于描述多物理场耦合效应下的自然对流传热传质的控制方程. 采用提出的格子Boltzmann模型结合多孔介质构造算法从孔隙尺度对规则以及随机多孔介质内双扩散效应的自然对流传热传质过程进行了模拟, 研究了不同瑞利数Ra, 不同孔隙率下的多孔介质内传热传质特征, 考察了温度梯度等因素对多孔介质内传质过程的影响, 创新地从孔隙尺度对多孔介质内的耦合对流扩散过程的传热传质机理进行了研究.     

一种抑制漂移-扩散方程数值解振荡的理论与方法

一、引言 半导体器件经典模型是漂移-扩散模型,在漂移项相对扩散项占优时,对该模型方程应用经典的有限元、有限差分或有限箱法离散,解会出现数值振荡现象,这种振荡本质上是由于数值离散造成的“负扩散”引起的。为抑制振荡和提高解的精度,人们提出了不少处理方法,如Scharfetter-Gummel方法(简称S-G方法)、SUPG方法.在一维情况下,SUPG和     

页岩微裂缝内气-水两相流动规律

气-水两相流动广泛存在于页岩气藏压裂液返排及页岩气生产阶段,直接决定页岩气藏压裂及开发效果.页岩微裂缝是储层流体运移的主要通道.然而由于物理实验尺度及精度的限制,准确模拟页岩微裂缝中的气-水两相流动存在巨大挑战.为探究页岩微裂缝中气-水两相流动规律,本文系统考虑了气-水两相流体多重微尺度效应,包括:(1)气相滑移;(2)水相滑移及边界层黏度变化;(3)气驱水条件下,存在沿壁面流动的水膜.建立气-水两相流体微尺度流动模型,采用侵入逾渗判断两相分布,求解真实页岩微裂缝内气-水相对渗透率,并验证模型的正确性.研究结果表明,微裂缝开度小于3μm时,微尺度效应影响不可忽略;微裂缝开度及水膜对气-水相对渗透率的影响均取决于各相流体所占流动空间的相对大小;气藏开采过程中压力降低,导致气相表观渗透率升高及微裂缝开度减小,均会造成气相相对渗透率减小;水膜的存在对气相流动能力造成影响的微裂缝开度界限为0.65μm.该研究揭示了页岩微裂缝内气-水两相流动规律,为评价页岩气井压裂及生产效果提供了理论基础.     

低渗透非均质砂砾岩的三维重构与水压致裂模拟

运用?CT成像、X射线衍射等方法,构建了考虑砂砾矿物成分与分布特征的非均质砂砾岩的三维重构模型,采用基于连续介质的离散元程序(CDEM)模拟分析了不同水平地应力下砂砾岩水力压裂裂缝的起裂与扩展行为.通过与均质砂岩计算结果的对比,阐述了水平应力比变化和砾石颗粒对砂砾岩压裂裂缝起裂、扩展和空间展布的影响.同时,完成了非均质砂砾岩和均质砂岩物理模型的水力压裂试验,并与数值模拟结果进行了对比.研究表明,砂砾岩的非均质性显著影响压裂裂缝的起裂和扩展行为.相同应力条件下砂砾岩压裂裂缝的起裂压力小于均质砂岩起裂压力,随水平应力比增大两者起裂压力的差异逐渐减小,非均质性强的砂砾岩的起裂压力受水平应力比变化的影响较小.无论非均质砂砾岩或均质砂岩,水平应力比变化均影响压裂裂缝的空间展布和形态.     

耗散粒子动力学方法在生物学领域的应用与研究进展:从蛋白质结构到细胞力学

耗散粒子动力学(dissipative particle dynamics, DPD)是近年发展起来的一种介观尺度的数值模拟方法,是研究软物质和复杂流体动力学行为的一种重要手段.这种新型介观模拟方法采用粗粒化粒子模型描述具有关联性的原子团或物质团,并通过简单的软排斥作用力描述粗粒化粒子间的相互作用,从而实现更大时间和空间尺度的复杂系统模拟计算,如油/水/表面活性剂体系、聚合物和胶体溶液的化学形态、微观形貌、相分离以及复杂流体流变特性的模拟等.本文首先介绍了DPD方法的理论框架,继而详细综述了DPD方法在生物系统中的应用.具体地,在分子尺度,我们重点介绍了该方法在蛋白质结构及其相互作用、两亲性脂质分子膜的结构与动力学、脂质膜与蛋白分子相互作用、纳米颗粒与脂质膜相互作用等方面的研究现状和研究热点.在细胞尺度,我们归纳了DPD方法在模拟血液微循环系统中血细胞的流动和血液流变学行为等方面的应用进展,包括红细胞的变形及流动,白细胞边聚及黏附行为,血小板边聚、黏附及聚集行为,健康与疾病状态下血液流变学特征,循环肿瘤细胞迁移、黏附及分选富集等.此外,我们总结了用于模拟血细胞变形及血液流动的其他数值模...     

一种优化的交错变网格有限差分法及其在井间声波中的应用

提出了一种优化交错变网格有限差分算法, 并在二维速度-应力关系的弹性波方程中实现. 利用频散关系守恒准则构造了四阶精度的差分算子, 该算法属于连续变网格方法, 不需要在精细网格和粗糙网格之间进行插值. 将优化算法的数值结果与解析解及八阶规则交错网格差分算法进行了比较, 验证了该算法的精度. 与基于Taylor展开的变网格有限差分算法比较可知, 优化算法的频散特性较好, 在数值模拟中可使用更粗糙的网格. 将提出的优化算法应用于复杂的井间声波模型. 该数值实例表明, 优化算法可以节省大量的计算内存和计算时间, 同时具有优良的稳定性.     

粒介质在类固-液相变过程中的时空参数特性

颗粒介质的类固-液相变过程除受材料性质影响外, 还与其运动速率和密集度有密切关系, 而接触时间数和配位数是表征该相变过程中颗粒间相互作用的重要时间和空间参量. 本文采用三维离散元方法对不同切变速率和密集度下颗粒介质的动力学行为进行了数值模拟, 确定了颗粒材料在类固-液相变过程中接触时间数和配位数的参数特性和演化过程, 并结合宏观应力的分布特性, 进一步确定了颗粒介质在液态和固态相互转化中的动力学机理, 讨论了颗粒介质在由快速流动向慢速、准静态转化的相变过程. 通过对颗粒单元相互作用的细观数值模拟, 获取了颗粒介质在宏观上的动力学行为, 为研究其在不同相态下的本构模型提供了依据.