多孔介质多重介质模型的表征范文

《工程力学杂志》2015年第S1期

摘要：

裂缝性岩体及缝洞型碳酸盐岩储层属于复杂多孔介质，岩体内部空隙为从纳米或微米级的微孔隙到数十厘米以上的大裂缝和溶蚀孔洞，空隙尺度跨越的范围很大。对于很多实际问题，由于单重介质模型计算不准确，且又不可能对复杂空隙空间结构在微观水平上进行精确描述，因此只能应用多重连续介质方法研究复杂多孔介质中的输运问题。该文给出了多重介质模型表征单元体的定义。在所研究问题的尺度范围内，这一多重介质模型表征单元体存在，复杂多孔介质的多重连续介质理论才能够成立。还以孔隙度和渗透率作为相关状态变量进行分析，研究了单重介质模型表征单元体和多重介质模型表征单元体间的联系，给出了复杂多孔介质多重介质模型的建立方法。

关键词：

多孔介质；三重介质模型；多重介质模型；表征单元体；连续介质方法；缝洞型油藏

多孔介质是固相单元(球形、片状、针状等颗粒)的集合体，各固相单元之间的空隙形成孔隙。孔隙内充满流体，是典型的固-流两相介质。在外荷载的作用下，由于固相骨架同孔隙中流相的相互作用，固-流两相介质的力学特性和数学描述完全不同于纯固相介质。由于多孔介质微观尺度上数学形态和固体-流体相互作用的复杂性，对于实际工程中土体和岩石这类地质材料，往往不可能以任何精确的数学方法在微观水平上进行描述，只能在宏观水平上应用连续介质的分析方法。基于达西定律的渗流力学方法和岩石力学中的弹塑性分析方法，都是连续介质理论在多孔介质力学方面的应用。

采用连续介质理论研究多孔介质，即是用一个假想的连续介质替代实际的、遍布空隙(微裂纹、微孔洞、孔隙、裂隙等)的多孔介质；将多孔介质物理参数等状态变量在某一体积元上的平均值作为该假想连续体参数的值。同单个孔隙空间相比，这一体积元的尺寸要足够大(超过某一临界尺寸)，以使平均出的物理参数具有代表性；这一临界值即是该参数对应的表征单元体(RepresentativeElementaryVolume，REV)。Bear等[1―2]和Prandtl等[3]指出，所研究问题的尺度远大于表征单元体的尺度时，则表征单元体存在，多孔介质的连续介质理论才能够成立。对于土体来说，砂粒颗粒的粒径一般在2mm以下，粘粒颗粒的粒径在0.005mm以下。因此，在大部分土力学工程问题涉及的尺度范围内，土体的表征单元体通常是存在的，因而连续介质理论的分析方法广为适用。对于岩体这类复杂多孔介质来说，情况则有很大的不同，裂隙岩体中遍布微裂纹、微孔洞、裂隙和裂缝等，不连续性、非均匀性、非弹性和各向异性的特点十分突出。在渗流力学方面，早期的研究都将岩体假设成为一个连续介质进行研究，即单重介质理论。

1960年Barenblatt等[4]提出了双孔隙度模型，即假定周期性分布的裂缝将岩石基质块分割开来，形成空间上重叠的两个各向同性的连续体：裂缝系统和基质系统；两个系统之间有流体输运。Warren和Root[5]则对Barenblatt的双孔隙度模型做了进一步的发展，除了不再忽略裂缝系统的储积能力外，还假设岩体中发育有均质的、正交的、互相连通的裂隙系统。为了更准确的描述裂缝系统同岩石基质系统的流体交换，Kazemi[6]提出了不仅依赖压力还依赖于压力梯度的交换函数。同样基于Barenblatt的双孔隙度模型，Hill和Thomas[7]提出了双渗透率模型，不仅考虑基质系统的储积能力，还充分考虑了基质块之间的流体输运。这些研究工作，形成了较为完整的双重介质理论。基于双重介质概念，很多学者提出了不同的三重介质(多重介质)模型[8―10]。Pruess和Narasimhan[11]针对裂缝岩石提出了多重相关连续介质模型，Wu等[12―13]归纳了的微孔隙、裂缝、溶孔不同组合模式，根据相应的概念模型，提出了适用于碳酸盐岩缝洞型多孔介质的三重介质模型。由于岩体内空隙从纳米或微米级的微裂缝和微孔隙到尺度数十厘米以上的大裂缝和溶蚀孔洞，空隙尺度跨越的范围很大，因此自20世纪90年代以来，这类复杂多孔介质的表征单元体的研究越来越引起关注[14―23]。

周创兵等[17―19]、陈卫忠等[20]、Kibok等[21]利用离散裂缝网络模型描述裂隙岩体，研究了裂隙岩体等效弹性参数的解析表达式，并提出估算裂隙岩体等效弹性模量REV的方法。Wang等[22]、Li等[23]、Kjetil等[24]对三维裂隙网络单相流体流动进行了渗流分析，应用渗流量等效原则、忽略基质和裂隙间的渗流耦合效应，计算了裂隙岩体的等效渗透率张量并估计了渗透率张量REV的大小。周创兵等[25―26]介绍了岩体表征单元体方面的研究进展，归纳、总结并提出了确定岩体REV的能量叠加法、地质统计法、数值模拟法3种方法。然而，上述这些研究本质上都是基于单重介质概念的，即将复杂多孔介质等效成一个连续体、并试图通过数值、解析或统计的方法研究其表征单元体的存在性和大小。由于岩体空隙空间的复杂性，三重介质理论能够更好地描述复杂多孔介质中的输运过程。各种多重介质理论从本质上来看，都是将复杂多孔介质中的多相流动问题，等效成为若干个连续介质中的多相流动问题，这若干个连续介质储集能力和渗透能力不同，空间上重叠在一起，它们之间有流体交换。

与单重介质模型将复杂多孔介质等效成一个连续体不同的是，三重介质模型中将岩体空隙空间划分成基质系统、裂缝系统和溶蚀孔洞系统，分别等效为在空间上重叠在一起的3个连续体(基质介质、裂缝介质和溶蚀孔洞介质)，这3个连续体上物理参数和力学变量等状态变量等效所涉及的表征单元体问题，相关研究在国内外还较少。本文以孔隙度作为代表性状态变量，首先研究了单重介质模型表征单元体，然后引出了多重介质模型表征单元体的概念，分析了复杂多孔介质的单重介质模型表征单元体和多重介质模型表征单元体间的联系；最后研究了复杂多孔介质渗透率对应的多重介质模型表征单元体，从理论上说明了多重介质模型的必要性和建立准则。

1单重介质模型的表征单元体

单重介质模型，即是把复杂多孔介质等效为单一的连续体，这个假想连续介质中的任一点上的状态变量是多孔介质中以该点为中心的某一体积内相应状态变量的统计平均(如图1所示)。如果这一体积大小的选择使得统计平均得到的状态变量是其中心点空间坐标和时间的连续函数，则该体积即是相应状态变量的表征单元体。

2多重介质模型的表征单元体

现以三重介质模型为例来说明多重介质模型REV的定义。三重介质模型将复杂多孔介质域()划分成基质系统()M、裂缝系统()F和溶蚀孔洞系统()V，并分别等效为在空间上重叠在一起的3个连续体：基质介质、裂缝介质和溶蚀孔洞介质。

3多重介质模型REV同单重介质模型REV的关系

复杂多孔介质在多重连续体上的每个内涵量都有一个对应REV。现以孔隙度为例，研究多孔介质三重介质模型的REV同其单重介质模型REV的关系。根据式(19)和式(20)，可分析多重介质模型REV同单重介质模型REV之间的关系。对于一复杂多孔介质来说，如果其单重介质模型REV为0U，则式(6)和式(7)成立，由式(19)和式(20)可知，式(16)和式(17)必成立，亦即多重介质模型REV必定存在且为0U。如果多重介质模型的REV存在且为0U，则式(16)和式(17)成立，由式(19)可知式(6)必成立，但由式(20)可知式(7)不一定成立，亦即在0U尺度上单重介质模型REV不一定存在。可见，采用连续介质理论描述多孔介质时，同单重介质模型相比，多重介质模型更能精细地描述多孔介质的输运过程。在尺度0U上，单重介质模型REV不存在而多重介质模型REV存在应满足的条件为。

4渗透率的REV

对于多孔介质不同的内涵量，其对应REV的大小和计算方法并不相同。孔隙度反映的是孔隙的空间展布，空隙空间具备数学上的线性叠加特性。渗透率是孔隙结构流体通过能力的平均统计特性，应以二阶张量表示，不能直接使用线性叠加原理。设k(x)表示复杂多孔介质域()内以点x为中心的孔隙管道的渗透率值。对以点0x为中心的空间区域(U)内若干这样的孔隙管道进行统计平均，可得到区域(U)上的等效渗透率张量。设0(,)iKxU是该等效渗透率张量的3个主值、()ikx是k(x)在等效渗透率张量3个主轴上的分量。根据式(32)可知，对于单重介质模型REV不存在仅是因为等效渗透率随空间位置变化不够平滑的情况(即式(24)成立而式(25)不成立)，如果能够通过对空隙空间合理的划分(见式(8))，使得式(27)成立，则多重介质模型REV存在，多重介质模型成立，复杂介质区域()内的多相流动问题可以等效成基质系统、裂缝系统和溶蚀孔洞系统这3个连续介质中的多相流动问题。将式(25)、式(27)和式(32)结合起来看，单重介质模型REV不存在而多重介质模型REV存在的原因在于渗透率随空间位置的变化过大(即式(32)右侧第2项大于0)。式(32)从机理上解释了双重介质模型或多重介质模型比单重介质模型能更准确地描述复杂孔隙介质中流体的流动。

5碳酸盐岩储层多重介质模型的建立

世界范围内，碳酸盐岩仅占沉积岩的20%，但却占油气探明储量的50%以上。近年来，碳酸盐岩油气藏已成为国内外勘探开发的主战场之一。缝洞型碳酸盐岩储层属于复杂多孔介质，岩体内空隙从纳米或微米级的微裂缝和微孔隙到尺度数10cm以上的大裂缝和溶蚀孔洞，空隙尺度跨越的范围很大。对于实际碳酸盐岩油藏的数值模拟，往往很难基于精确的数学方法对复杂多孔介质空隙空间内多相流体的流动进行计算，只能在宏观水平上应用连续介质的分析方法(单重介质模型或者多重介质模型)来模拟。假设某一缝洞型碳酸盐岩储层()的空隙空间结构和分布已知，对这一储层进行油藏数值模拟时的研究尺度(网格尺度)为L。下面选择孔隙度和渗透率作为状态变量，根据本文的多重介质模型REV理论来建立该碳酸盐岩储层的三重连续介质模型。

首先，尝试建立单重介质模型。由于碳酸盐岩储层空隙域(p)已知，对()内的每一点可由式(6)、式(7)、式(24)和式(25)确定孔隙度和渗透率对应的单重介质模型表征单元体U0。如果U0远比研究区域的尺度小，即满足下式：30UL(33)则在研究区域内单重介质模型REV存在，可以将该复杂多孔介质等效为单一的连续体，即储层的单重介质模型成立。如果表征单元体U0不满足式(33)，则表明在所研究问题的尺度范围内表征单元体不存在，不能将该储层等效为单一的连续体，单重介质模型不成立。由前述分析可知，单重介质模型REV不存在而多重介质模型REV存在满足式(21)的条件。因此，在单重介质模型REV不存在的情况下，可尝试建立多重介质模型。建立多重介质模型时，主要问题是空隙域的划分。按照式(8)，将储层()划分成基质系统()M、裂缝系统()F和溶蚀孔洞系统()V；与此划分对应，可在()内的每一点上根据式(16)、式(17)、式(26)和式(27)确定多重介质模型的表征单元体0U。如果存在这样一种划分，它对应的表征单元体0U满足式(33)，则表明在研究区域内三重介质模型REV存在，可以将该储层的多相流动问题等效成为3个连续介质中的多相流动问题，即三重介质模型成立；相应的划分就是三重介质模型的建立方法。如果这样的划分不存在(通常都是由缝洞型碳酸盐岩储层中大裂缝和大的溶蚀孔洞引起的)，则需要把这些大尺度的裂缝和溶洞从()中剔出出来，作为独立的连续区域单独处理；从而使得改变后的()内，多重介质模型REV存在。

6结论

在所研究问题的尺度范围内表征单元体存在，多孔介质的连续介质理论才能够成立。单重介质模型成立的充要条件是单重介质模型REV存在。多重介质模型，即是将复杂多孔介质中的多相流动问题，等效成为若干个连续介质中的多相流动问题，这若干个连续介质储集能力和渗透能力不同，空间上重叠在一起，它们之间有流体交换。多重介质模型成立的充要条件是多重介质模型REV存在。本文以孔隙度和渗透率为例，给出了多重介质模型REV的定义，分析了复杂多孔介质的单重介质模型REV和多重介质模型REV间的联系。在所研究问题的尺度范围内，如果单重介质模型REV存在，则多重介质模型REV必定存在；如果多重介质模型REV存在，而单重介质模型REV则不一定存在。并给出了单重介质模型REV不存在而多重介质模型REV存在应满足的条件。如果单重介质模型REV不存在是因为任一点在不同邻域上的等效渗透率波动较大造成的，则相应的多重介质模型REV也不存在，此时无论是单重介质模型还是多重介质模型都不能用来描述多孔介质中的流动。对于单重介质模型REV不存在仅仅是因为等效渗透率随空间位置变化不够平滑的情况，则有可能通过对空隙空间合理的划分，使得多重介质模型REV存在，多重介质模型成立。本文通过数学公式，从机理上解释了多重介质模型比单重介质模型能更精确地描述复杂孔隙介质中流体的流动，同时以缝洞型碳酸盐岩储层为例，给出了多重介质模型的建立方法。

参考文献：

[1]BearJ.DynamicsofFluidsinPorousMedia[M].NewYork:AmericanElsevier,1972:16―21.

[2]BachmatY,BearJ.Macroscopicmodelingoftransportphenomenainporousmedia:1.Thecontinuumapproach[J].TransportinPorousMedia,1986,1(3):213―240.

[3]PrandtlL,TietjensOG.FundamentalsofHydroandAeromechanics[M].NewYork:McGraw-Hill,1934.

[4]BarenblattGI,ZheltovIP,KochinaIN.Basicconceptsinthetheoryofseepageofhomogeneousliquidsinfissuredrocks[J].JournalofAppliedMathematicsandMechanics,1960,24(5):852―864.

[5]WarrenJE,RootPJ.Thebehaviorofnaturallyfracturedreservoirs[J].SocietyofPetroleumEngineersJournal,1963,228:245―255.

[6]KazemiH.Pressuretransientanalysisofnaturallyfracturedreservoirswithuniformfracturedistribution[J].SocietyofPetroleumEngineersJournal,AIME,1969,246:451―462.

[7]HillAC,ThomasGW.ANewApproachforSimulatingComplexFracturedReservoirs[C].SPE13537-MSMiddleEastOilTechnicalConferenceandExhibition,Bahrain,1985.

[8]ClosemannPJ.Theaquifermodelforfissuredfracturedreservoir[J].Soc.Pet.Eng.J.,1975,385―398.

[9]AbdassahD,ErshaghisI.Triple-porositysystemforrepresentingnaturallyfracturedreservoirs[J].SPEFormationEvaluation,1986,1:113―127.

[10]BaiM,ElsworthD,RoegiersJC.Multiporosity/multipermeabilityapproachtothesimulationofnaturallyfracturedreservoirs[J].WaterResourcesResearch,1993,29:1621―1633.[11]PruessK,NarasimhanTN.Apracticalmethodformodelingfluidandheatflowinfracturedporousmedia[J].SocietyofPetroleumEngineersJournal,1985,25:14―26.[12]WuYS,LiuHH,BodvarssonGS.Atriple-continuumapproachformodelingflowandtransportprocessesinfracturedrock[J].JournalofContaminantHydrology,2004,73():145―179.

[13]Yu-ShuWu,YuanDi,ZhijiangKang,PeraponFakcharoenphol.Amultiple-continuummodelforsimulatingsingle-phaseandmultiphaseflowinnaturallyfracturedvuggyreservoirs[J].JournalofPetroleumScienceandEngineering,2011,78(1):13―22.

[14]PintoA,DaCH.ScaleEfectsinRockMasses[M].Rotterdam:BalkemaAAPublishers,1993.

[15]卢波,葛修润,朱冬林,陈剑平.节理岩体表征单元体的分形几何研究[J].岩石力学与工程学报,2005,24(8):1355―1361.LuBo,GeXiurun,ZhuDonglin,ChenJianping.Fractalstudyontherepresentaltiveelementaryvolumeofjointedrockmasses[J].ChineseJournalofRockMechanicsandEngineering,2005,24(8):1355―1361.(inChinese)

[16]宁宇,徐卫亚,郑文棠,等.柱状节理岩体随机模拟及其表征单元体尺度研究[J].岩石力学与工程学报,2008,27(6):1202―1208.NingYu,XuWeiya,ZhengWentang,etal.Studyofrandomsimulationofcolumnarjointedrockmassesanditsrepresentativevolumescale[J].ChineseJournalofRockMechanicsandEngineering,2008,27(6):1202―1208.(inChinese)

[17]周创兵,熊文林.双场耦合条件下裂隙岩体的渗透张量[J].岩石力学与工程学报,1996,15(4):338―344.ZhouChuangbing,XiongWenlin.Permeabilitytensorforjointedrockmassesincoupledseepageandstressfield[J].ChineseJournalofRockMechanicsandEngineering,1996,15(4):338―344.(inChinese)

[18]周创兵,於三大.论岩体表征单元体积REV岩体力学参数取值的一个基本问题[J].工程地质学报,1999,7(4):332―336.ZhouChuangbing,YuSanda.RepresentativeelementaryvolumeREV—Afundamentalproblemforselectingthemechanicalparametersofjointedrockmass[J].JournalofEngineeringGeology,1999,7(4):332―336.(inChinese)

[19]荣冠,周创兵,王恩志.裂隙岩体渗透张量计算及其表征单元体积初步研究[J].岩石力学与工程学报,2007,26(4):740―746.RongGuan,ZhouChuangbing,WangEnzhi.Preliminarystudyonpermeabilitytensorcalculationoffracturedrockmassanditsrepresentativeelementaryvolume[J].ChineseJournalofRockMechanicsandEngineering,2007,26(4):740―746.(inChinese)

[20]陈卫忠,杨建平,邹喜德,周春宏.裂隙岩体宏观力学参数研究[J].岩土力学与工程学报,2008,27(8):1569―1575.ChenWeizhong,YangJianping,ZouXide,ZhouChunhong.Researchonmacromechanicalparametersoffracturedrockmasses[J].ChineseJournalofRockMechanicsandEngineering,2008,27(8):1569―1575.(inChinese)

[21]KibokM,JingLR.Numericaldeterminationoftheequivalentelasticcompliancetensorforfracturedrockmassesusingthedistinctelementmethod[J].InternationalJournalofRockMechanicsandMiningSciences,2003,40(6):795―8l6.

[22]WangM,KulatilakeaPHSW,UreaJ,etal.EstimationofREVsizeandthree-dimensionalhydraulicconductivitytensorforafracturedrockmassthroughasinglewellpackertestanddiscretefracturefluidflowmodeling[J].InternationalJournalofRockMechanicsandMiningSciences,2002,39:887―904.

[23]LiJH,ZhangLM,WangY,FredlundDG.Permeabilitytensorandrepresentativeelementaryvolumeofsaturatedcrackedsoil[J].CanadianGeotechnicalJournal,2009,46(8):928―942.

[24]KjetilN,PhilipSR.Identifyingtherepresentativeelementaryvolumeforpermeabilityinheterolithicdepositsusingnumericalrockmodels[J].MathematicalGeosciences,2008,40(7):753―771.

[25]周创兵,陈益峰,姜清辉.岩体表征单元体与岩体力学参数[J].岩土工程学报,2007,29(8):1135―1142.ZhouChuangbing,ChenYifeng,JiangQinghui.Representativeelementaryvolumeandmechanicalparametersoffracturedrockmasses[J].ChineseJournalofGeotechnicalEngineering,2007,29(8):1135―1142.(inChinese)

[26]向文飞,周创兵.裂隙岩体表征单元体研究进展[J].岩石力学与工程学报,2005,24(Suppl2):5686―5692.XiangWenfei,ZhouChuangbing.Theadvancesininvestigationofrepresentativeelementaryvolumeforfracturedrockmass[J].ChineseJournalofRockMechanicsandEngineering,2005,24(Suppl2):5686―5692.(inChinese)

作者：邸元 康志江 代亚非 WU Yu-Shu 单位：北京大学工学院 中国石油化工股份有限公司石油勘探开发研究院 科罗拉多矿业大学石油工程系