储油（气）岩石的渗透率.ppt

第三章 储油（气）岩石的渗透率,重点解决以下几个问题,孔隙介质中流体渗流特征的刻画 （达西定律及其三种不同的表达方式：折算压力表达式、一般表达式、微分形式） 2. 达西定律的含义、适用范围、及其校正 3. 渗透率的获取方法,概念,在一定的压差下，岩石允许流体通过的性质称为岩石的渗透性。从数量上度量岩石渗透性的参数就叫岩石的渗透率。渗透率就是岩石允许流体通过的能力。,基本知识,孔隙度-度量岩石储存能力的参数，它是一个没有方向性的标量。 渗透率-度量岩石渗透能力的参数，是一个具有方向性的向量。 渗透率又可分为：绝对渗透率、相渗透率与相对渗透率。,本章着重讨论绝对渗透率，相渗透率与相对渗透率将在第三篇中介绍。,第三章 储油（气）岩石的渗透率,1 达西定律及其表达式,一、达西定律及其表达式 达西定律是1856年法国亨利达西在解决城市供水问题时，用直立均质未胶结砂柱做水流渗滤试验，得出的一个经验公式，后人为纪念他，把这一公式命名为达西公式或称达西定律。,在砂柱中，顶底分别用渗 透性铁丝网封住，紧靠砂柱顶 底分别与测压管相连接，当水 流通过砂柱时，水在测压管内 分别上升到相对于任一基准面 以上h1和h2的高度， 实验中发现，无论砂柱中 砂层类型如何改变，流量总是 与测压管水柱高差、及砂柱横 截面积成正比，而与砂柱的长 度成反比。,式中； Q 总流量； A 截面积； v渗流速度，可以理解为单位时间内单位截面积的注入量（cm/s）； h相对于某个基准面压力计的液面高差（cm）； k比例常数，也叫介质的渗流系数（cm2）。,1）渗流的液体是均质的、不可压缩的水，水的粘度不变， 因此没有考虑粘度对渗流规律的影响； 2）均质砂柱由极细小的细砂组成，具微小的连通孔隙通道， (达西改变砂子类型，实际上仅改变了k的大小)； 3）渗流速度较小，且变化不大； 4）试验装置始终保持在垂直条件下；,在该项实验中，其边界条件如下：,之后，曾有他人在改变边界条件4（即将实验装置摆放成各种角度的倾斜位置）重复进行达西实验，结果发现不管装置倾斜程度如何，只要测验管水头差（h1-h2）相同，则流量相同,另外，人们通过改变边界条件2，用实际岩心代替砂柱进行实验，证明达西定律是成立的，但介质特性（k）对流量有影响； 当在改变边界条件1时，即用各种液体而不仅仅是水作实验时达定律仍成立，但发现流体粘度对流量有影响； 因此达西公式进一步表示为：,上述实验表明，不管如何改变边界条件，达西定律是成 立的。改变不同介质与流体所导致的对流量的影响主要是因 为渗流系数发生了改变。 因此原始达西公式中的k只代表了某种特定流体在特定介质条件下的渗流能力。 由此可看出，不同的流体、不同的介质条件，其渗流系 数是不同的。,我们可将水头高h1、h2分别换算成液面h高度时的压力 Pr1和Pr2（称为折算压力），即： Pr1=gh1 Pr2=gh2 将上述折算压力代入达西公式，即：,达西公式中的h1和h2代表了渗流液体液面相对于某一基准面水柱的高度,注：Pr的大小与选用的基准面有关，称为基准压力或折算压力该公式实际上是以压力形式表示油层中各点液体所具有的总能量）,h1=Pr1/g h2=Pr2/g,达西公式,设k=Kg K=k/g，则,此公式即为达西公式的折算压力表达式,由于总水压头（总能量） PrghP（压力计压能）+gZ（势能） 故 Pr1gh1P1gZ1 Pr2gh2P2gZ2,代入达西折算压力公式：,当Z=0时，即流体为水平流动时,达西定律的 一般表达式,变换上式，得：,K即为岩石的渗透率（cm2） 当流体性质不变情况下，岩石渗透率仅仅是与多孔介质（岩石性质）有关的参数。 上述达西公式是均质孔隙介质中单相流体在作直线稳定渗流情况下推导出来的一个平均关系式。如不满足上述条件，岩石的渗透率就会降低。实际上，孔隙介质是不均匀的，流体在孔隙介质中的渗流也常常表现为非稳定的线性渗流。经大量实验证明，很多渗流是符合达西定律的。但对于高速流动的液体，以及速度极低或极高的气体，达西定律就不适用了。,对于实际中不均匀的孔隙介质，加上不均质的流体（即多相）流体同时渗流时，常作非平面、非稳定的线性渗流。大量实验证明，达西定律也是适用的。 达西公式的一般表达式为：,当岩样水平时，流体作水平渗流，Z1-Z2=0，则：,二、达西公式的推广,（一）达西公式的微分方程,式中，当Pr，L无限小时，可写成：,上式即为达西公式的微分形式，公式前面的负号代表压力 增加的方向与渗流距离增加的方向相反。即在渗流方向上， dPr/dL应该是负值。 由于Pr=PgZ 代入上式得：,这是达西定律的微分形式,（二）不可压缩液体渗流的达西公式表达式,前面介绍的公式是建立在一块岩心实验基础上的，并且认为这块岩心的孔隙介质由均质介质组成，流体在内部的渗流向一个方向。 实际上，地下流体的渗流是相当复杂的，下面主要讨论几种简单渗流方式的达西公式表达式。,1. 水平线性稳定渗流,水平线性稳定渗流的达西定律的基本表达式,从达西定律一般表达式推导，Z1=Z2（水平），代入一般表达式,从达西定律的微分形式推导， Z1=Z2（水平），代入达西定律微分形式的一般表达式,分离变量：,积分,可以看出： 由微分方程所导出的 水平线性稳定流表达 式与根据达西公式一 般表达式所得出的结 果是一致的。,2. 平面径向渗流,平面径向渗流的达西定律的基本表达式,式中： h地层厚度（m）； 外边界压力（Pa）； 内边界压力（m）； 外边界半径（m）； 内边界半径（m）。,参数的物理含义,（三）达西公式的修正 可压缩气体的达西公式,可压缩气体的最大特点是：当压力减小时，气体会发生 膨胀，温度一定时气体的膨胀服从波义尔定律：,只要将流量用平均流量代替即可,水平线性 稳定渗流,平面径向 稳定渗流,三、达西定律的适用范围,对大多数油田开发实践中，油气渗流一般服从达西定律，但对于高速流动的流体，尽管边界条件不变，但流型会变得瞬息万变，会产生涡旋，这种流速变大而导致的流型改变的转换可用“临界点”来加以描述。流速在该点以下时，流体以定常流的型式流动，称为层流，当流速超过“临界点”时，流线会变成非定向，不规则的流动型式，称为“紊流”（或湍流）。这二种不同的流动型式具有不同的渗流特性。,当渗流速度增大到一定值后，流速与压力梯度关系由线性转变为非线性，即流动型式从线性渗流转变为非线性渗流。达西定律就不适用了。,对于低渗透性致密岩石，在低速渗流时，由于流体与岩石之间存在吸附作用，或在粘土矿物表面形成水膜，当压力梯度很低时，流体不流动，因此存在一个启动压力梯度a，在低于该压力梯度范围内流速与压力梯度不呈线性关系,卡佳霍夫提出的判断指标雷诺数,式中：Re雷诺数，反映了惯性力与粘性力的比值，也反映了孔隙介质的特点； 流体密度（ ）；1750单位换算系数，与规定的各物理量的单位有关。,当 Re 0.2-0.3时渗流服从达西定律； 当 Re 0.2-0.3时，则渗流规律受到破坏，这时渗流速度和压差的关系如下式：,(Re),在已知岩石和流体物理参数如岩石孔隙度、渗透率和流体密度、粘度条件下，定义,一般，临界雷诺数Red=0.2-0.3,C取决于岩层和流体性质的系数； n渗流指数。,n=1时，为线性渗流，相当于管路水力学中的层流。这时渗流速度很小，即ReRed（临界雷诺数），液体在多孔介质内所产生的惯性力极小。 1/2n1时，为渗流过渡区。这时渗流速度已相当大（即ReRed），流体在多孔介质中惯性力已明显表示出来，故直线渗流定律已破坏。 n=1/2时，为渗流的平方区（类似于管路水力学中紊流平方区）。这时渗流速度已很大（ReRed），惯性力也很大。在这以后惯性力的增加与压力的下降又成不变的比例关系。,气体在致密岩石中低速渗流时会产生滑动效应 克林肯博格效应,2）同一岩石，同一平均压力，不同气体测得的渗透率不同,实验发现二点,1）同一岩石、同一种气体，在不同的平均压力下测得的气体渗透率不同，低平均压力下气体渗透率比较高，高平均压力下气体渗透率比较低（与实验相结合思考）,达西定律的修正,气体渗透率与平均压力的关系,同一岩石应该只有一个绝对渗透率，为什么测试条件不同（压力和气体类型）就会产生不同的Ka，应该选取哪一个Ka值作为岩石的绝对渗透率？,液体流动液测岩石渗透率的达西公式是建立在液体（确切讲是牛顿流体）渗流实验基础上的，认为液体的粘度不随流动状态改变，即所谓的粘性流动。液体流动时液体与管壁分子间会产生粘滞阻力，由于液-固间的分子力比液-液间的分子力大，在管壁附近表现出的粘滞阻力最大，可使得管壁处液体的流速为零，管道中心粘滞阻力最小，流速最大。管内流速分布是圆锥曲线。,克林博格从分析孔隙内气、液流速分布入手解释了这种现象,气体流动孔道壁表面的气体分子与孔道中心的分子流速几乎没有差别。气测渗透率时，除了气固间的分子作用力小以外，相邻层的气体分子还可以由于动能交换而使得管壁处的气体分子层与孔道中心的分子层的流速被不同程度均一化。管壁处的气体分子层流速不为零形成“气体滑脱效应”。克林博格发现了这种效应，因而也称为“克氏效应”。,气体的这种滑脱效应还与气体的性质有关，不同气体的分子量不同，分子直径不同，自由行程也就不同，使得滑脱系数（b）不同，分子量小，则b大，滑脱效应严重。,1、用气体测得的渗透率要比用液体测定的渗透率值高 气体法测定的渗透率更能反映岩石的真实渗透率 2、平均压力小，气体密度就小，气体分子间的碰撞就少，使得气体更易流动，气体滑脱现象越严重，所测的渗透率值越大。相反平均压力增大，气体滑脱效应就消失，渗透率就越小。如果压力增大到无穷大，气固间的作用力增大，管壁上的气膜逐渐趋于稳定，气体的流动性质接近液体的流动性质，这时的渗透率趋于一个常数，接近液测渗透率，故称该渗透率为等效液体渗透率或克氏渗透率。,气测渗透率与液测渗透率的差异：,同一岩石，不同气体测得的渗透率和平均压力的直线关系交纵座标于一点，该点的气体渗透率与同一岩石的液体渗透率是等价的。故称为等价液体渗透率，也称克林肯博格渗透率。,如果平均压力增至无穷大时，渗透率不再变化而趋于一个常数，这个数值一般接近于液测渗透率，这是因为压力无穷大时，气体的性质已经接近于液体的流动性质，气固间的作用力增大，气体的滑脱效应逐渐消失，管壁上的气膜逐渐稳定了。,1941年，克林博格给出了考虑气体滑脱效应的气测渗透率表达式,岩石本身渗透率的高低也将影响气体渗透率大小,我国多采用室内实验室测定过程中，改变几次平均压力，再按达西公式计算Kg，并编制Kg与平均压力的关系曲线，外推与Kg轴相交，便可得到克氏渗透率，以此作为岩石的绝对渗透率,四、储油（气）岩石渗透率的内在含义,绝对渗透率是岩石本身具有的固有性质，它只与岩石的孔隙结构有关，与通过岩石的流体性质无关。对于这一性质，可以从渗透率与裂缝宽度的关系、渗透率与毛细管孔隙半径的关系来理解。,孔隙度是岩石中孔隙体积与岩石外表体积的比值。这是岩石的内在性质。,一、渗透率与裂缝宽度的关系,对于单一裂隙来说，其长度和截面分别为L和lb，当压差为p时，可由Poiseulle方程求出流量q1:,对于n个裂缝来说，岩石截面积为A，流量为：,根据达西定律,二、渗透率与毛细管孔隙半径的关系,L,r,如果我们把孔隙看成是一根根等直径的毛细管束（n根），则可写成：,对于单根毛细管（半径为r 、长度为L），流体在其中流动时的流量，根据伯叶稷（Poisenille）定律可求得，即：,如果岩石长度为L，岩石截面积为A，溶孔孔隙度为孔，从孔隙度定义得：,L,r,按达西公式岩石的流量应为：,如果考虑孔隙结构的复杂性,是孔隙结构系数（曲折度）,五、渗透率的因次（量纲）与单位,1. 因次,目前在渗透率单位制中还得保留有渗流力学中过去常用的一套混合单位制，混合单位制中达西单位制用达西（D)表示，混合单位制中工程单位制不管公制或英制均用毫达西（mD）。,2. 单位,从因次中可知渗透率的单位是长度的平方，与面积的单位相同。,石油系统渗透率使用的绝对单位制：平方米太大，通常采用平方微米（ ）来表示。 也很大，最常用的单位是,所谓1达西（D）的含义是：以粘度为1厘泊（cp）的流体完全饱和于岩石孔隙中，在1个大气压的压差下，以层流方式通过截面积为1平方厘米，长度为1厘米的岩样时，若其流量为1立方厘米/s，则岩石的渗透率为1达西（英文字为Darcy,简写D）。,但这个数值太大，使用起来不方便，因而采用它的千分之一，称为毫达西（英文字母millidarcy简写mD）来作为渗透率的单位。,混合单位与绝对单位的换算,六、储油岩石渗透率的大小,渗透率是储油（气）岩物性的重要参数，不论对油气运移聚集，还是油（气）田开发都是基础数据。但其数值在不同的油（气）层中相差很大，由几个毫达西到几千个毫达西。,我国和世界一些油田的渗透率（据徐献中）,储层物性分级 （中国油气储层研究1994）,2 储层岩石渗透率的测定与计算,求取渗透率的方法包括二类： 1、以达西定律为基础的实验测定方法 2、间接法（测井方法或油藏工程方法）,1、实验法,分常规岩心分析和全直径岩心分析，在流动方向上分为水平、垂直和径向渗透率的测定,1）常规（小）岩心分析,指直径2.5cm（或3.8cm），长4-6cm的标准岩心柱。根据能源部石油行业标准SY5336-88标准，国内各油田主要采用气测法测定常规岩心渗透率。 分析仪器主要有：流量管气测渗透率仪、高低渗透率仪、全自动孔渗分析仪等,流量管空气渗透率法（或空气压水法测渗透率） （国内普遍采用的方法）,仪器由二部分组成：岩心夹持器和插于水池内的流量管。 岩心顶端通大气，用抽气球抽空时，迅速将水池水面吸升至某一高度，待抽空停止后，流量管内水面在静水压头作用下下落，水面下降造成负压从而建立起二端压差，外界空气进入通过岩样进入流量管。岩心出口压力是一个变量，取决于流量管内水柱静水压头h。水面下降快慢与岩心渗透率有关，只要测出流量管内水面下降某一段距离h0-hx的时间T便可求出K。,美国岩心公司的气测渗透率仪流程图及其使用该仪器测定渗透率的公式,2）全直径岩心分析,在钻井取出的岩心上截取一段柱状岩心，其直径一般为100-120mm，采用的仪器及其原理与常规岩心是一样的，唯一的区别是岩心夹持器，全直径岩心夹持器可用于测定岩心水平方向和垂直方向的渗透率,可在达西公式中引进一个附加的形状系数（筛网形状系数）E或G进行修正。该形状系数主要取决于弓形滤网所遮挡的角度和岩样的直径，可查曲线获得。,水平方向的渗透率计算,垂向渗透率计算公式与常规小岩心公式相同,岩心公司计算公式,全直径平面径向渗透率测定,在全直径岩心中心钻孔，上柱塞气孔与岩心中心孔相连。气体从环形空间流入岩心，并从岩心中心孔流出,2、间接法求取渗透率测井曲线法,利用测井曲线间接求取渗透率的原理是：束缚水饱和度Swi与岩心渗透率有关，岩石孔隙直径越小， Swi越高，渗透率K越低，同时对砂岩储层，孔隙度高的储层渗透性一般也好，因此可由测井得到的孔隙度和束缚水饱和度Swi资料估算储层渗透率。,2、间接法求取渗透率 基于毛管束渗流模型推导的计算模型,毛细管岩石模型的孔隙度为,毛细管岩石模型的比表面为,代入Kozeny公式,n为单位截面上毛细管数,3 非均质储油（气）层平均渗透率的计算,一、不均一地层平均渗透率的计算,已知地层总厚度为各层厚度之和， 即： h= h1+ h2+ h3 地层流体总流量Q为各层流量之和， 即： Q= Q1+ Q2+ Q3,1纵向不均一 （1）平面线性稳定渗流：,设有三个渗透率不等的地层，油气作平面线性渗