资源描述
单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第三章 储油(气)岩石的渗透率,概念,在一定的压差下,岩石允许流体通过的性质称为岩石的,渗透性,。从数量上度量岩石渗透性的参数就叫岩石的,渗透率,。,渗透率,就是岩石允许流体通过的能力。,基本知识,孔隙度,-,度量岩石储存能力的参数,它是一个没有方向性的标量。,渗透率,-,度量岩石渗透能力的参数,是一个具有方向性的向量。,渗透率又可分为:绝对渗透率、相渗透率与相对渗透率。,本章着重讨论绝对渗透率,相渗透率与相对渗透率将在第三篇中介绍。,岩石中只有一种流体通过时,岩石允许该流体通过的能力称为,单相渗透率,。,绝对渗透率,是指当岩石中只有一种流体通过,且流体不与岩石发生任何物理和化学反应时,岩石允许该流体通过的能力。,实质上任何一种流体都会或多或少地与岩石发生物理和化学反应。绝对渗透率只是一个理论值。在实际应用中,只能选用一种与岩石反应非常少的流体的单相渗透率来近似代替绝对渗透率。,通常采用气体,氩气、氮气、空气,的渗透率作为绝对渗透率。,绝对渗透率是岩石本身具有的固有性质,它只与岩石的孔隙结构有关,与通过岩石的流体性质无关。,严格说,自然界中任何物质均具有一定渗透性。如致密钢板,在超高压条件下,也可以让气体通过。,但实际上自然界达不到那样高的压力,一些粘土岩、页岩和致密的碳酸盐岩等由于,孔喉直径太小(微毛细管孔隙),即使有部分流体进入,也由于颗粒表面张力作用而被牢固地吸附在颗粒表面,从而堵塞微毛细管,使流体无法通过。,渗透性与非渗透性是个相对的概念,1,达西定律及其表达式,一、达西定律及其表达式,达西定律是,1856,年法国亨利,达西在解决城市供水问题时,用直立均质未胶结砂柱做水流渗滤试验,得出的一个经验公式,后人为纪念他,把这一公式命名为达西公式或称达西定律。,在砂柱中,顶底分别用渗,透性铁丝网封住,紧靠砂柱顶,底分别与测压管相连接,当水,流通过砂柱时,水在测压管内,分别上升到相对于任一基准面,以上,h1,和,h2,的高度,,实验中发现,,无论砂柱中,砂层类型如何改变,流量总是,与测压管水柱高差、及砂柱横,截面积成正比,而与砂柱的长,度成反比。,式中;,Q,总流量;,A,截面积;,v,渗流速度,可以理解为单位时间内单位截面积的注入量(,cm/s,);,h,相对于某个基准面压力计的液面高差(,cm,);,k,比例常数,也叫介质的渗流系数(,cm,2,)。,1,)渗流的液体是均质的、不可压缩的水,水的粘度不变,,因此没有考虑粘度对渗流规律的影响;,2,)均质砂柱由极细小的细砂组成,具微小的连通孔隙通道,,(,达西改变砂子类型,实际上仅改变了,k,的大小,),;,3,)渗流速度较小,且变化不大;,4,)试验装置始终保持在垂直条件下;,在该项实验中,其边界条件如下:,之后,曾有他人在改变边界条件,4,(即将实验装置摆放成各种角度的倾斜位置)重复进行达西实验,结果发现不管装置倾斜程度如何,只要测验管水头差(,h1-h2,)相同,则流量相同,从而证明达西实验定律在地球重力场中与流动方向无关,其主要影响力是重力,是一种势能。,另外,人们通过改变边界条件,2,,用实际岩心代替砂柱进行实验,证明达西定律是成立的,但介质特性(,k,)对流量有影响;,当在改变边界条件,1,时,即用各种液体而不仅仅是水作实验时达定律仍成立,但发现流体粘度对流量有影响;,因此达西公式进一步表示为:,上述实验表明,不管如何改变边界条件,达西定律是成,立的。改变不同介质与流体所导致的对流量的影响主要是因,为渗流系数发生了改变。,因此原始达西公式中的,k,只代表了某种特定流体在特定介质条件下的渗流能力。,由此可看出,不同的流体、不同的介质条件,其渗流系,数是不同的。,达西公式中的,h1,和,h2,代表了渗流液体液面相对于某一基准面水柱的高度。,h1,h2,Z1,Z2,我们可将水头高,h1,、,h2,分别折算成液面,h,高度时的压力,Pr1,和,Pr2,(称为折算压力),即:,Pr1=,gh1,Pr2=,gh2,将上述折算压力代入达西公式,,即:,注:,Pr,的大小与选用的基准面有关,称为基准压力或折算压力,该公式实际上是以压力形式表示油层中各点液体所具有的总能,量),h1,=Pr1/,g,h2=Pr2/,g,达西公式,设,k=K,g K=k/,g,,则,此公式即为达西公式的折算压力表达式,由于总水压头(总能量),Pr,gh,P,(压力计压能),+,gZ,(势能),故,Pr1,g,h1,P1,gZ1,Pr2,g,h2,P2,gZ2,h1,h2,Z1,Z2,代入达西折算压力公式:,当,Z=0,时,即流体为水平流动时,达西定律的,一般表达式,h1,h2,Z1,Z2,变换上式,得:,K,即为岩石的渗透率(,cm,2,),当流体性质不变情况下,岩石渗透率仅仅是与多孔介质(岩石性质)有关的参数。,上述达西公式是均质孔隙介质中单相流体在作直线稳定渗流情况下推导出来的一个平均关系式。如不满足上述条件,岩石的渗透率就会降低。实际上,孔隙介质是不均匀的,流体在孔隙介质中的渗流也常常表现为非稳定的线性渗流。但经大量实验证明,很多渗流是符合达西定律的。,但对于高速流动的液体,以及速度极低或极高的气体,达西定律就不适用了。,对于实际中不均匀的孔隙介质,加上不均质的流体(即多相)流体同时渗流时,常作非平面、非稳定的线性渗流。大量实验证明,达西定律也是适用的。,达西公式的一般表达式为:,当岩样水平时,流体作水平渗流,,Z1-Z2=0,,则:,二、达西公式的推广,(一)达西公式的微分方程,式中,当,Pr,,,L,无限小时,可写成:,上式即为达西公式的微分形式,公式前面的,负号代表压力,增加的方向与渗流距离增加的方向相反,。即在渗流方向上,,dPr/dL,应该是负值。,由于,Pr=P,gZ,代入上式得:,这是达西微分方程的一般表达式,(二),不可压缩液体,渗流的达西公式表达式,前面介绍的公式是建立在一块岩心实验基础上的,并且认为这块岩心的孔隙介质由均质介质组成,流体在内部的渗流向一个方向。,实际上,地下流体的渗流是相当复杂的,下面主要讨论几种简单渗流方式的达西公式表达式。,1.,水平线性稳定渗流,水平线性稳定渗流的达西定律的基本表达式,从达西定律一般表达式推导,,,Z1=Z2,(水平),代入一般表达式,从达西定律的微分形式推导,,,Z1=Z2,(水平),代入达西定律微分形式的一般表达式,分离变量:,积分,可以看出:,由微分方程所导出的,水平线性稳定流表达,式与根据达西公式一,般表达式所得出的结,果是一致的。,2.,平面径向渗流,平面径向渗流,的达西定律的基本表达式,式中:,h,地层厚度(,m,);,外边界压力(,Pa,);,内边界压力(,m,);,外边界半径(,m,);,内边界半径(,m,)。,参数的物理含义,(三)达西公式的修正,可压缩气体的达西公式,可压缩气体的最大特点是:当压力减小时,气体会发生,膨胀,,温度一定,时气体的膨胀服从波义尔定律:,只要将流量用平均流量代替即可,水平线性稳定渗流,平面径向稳定渗流,三、达西定律的适用范围,对大多数油田开发实践中,油气渗流一般服从达西定律,但对于高速流动的流体,尽管边界条件不变,但流型会变得瞬息万变,会产生涡旋,这种流速变大而导致的流型改变的转换可用“临界点”来加以描述。流速在该点以下时,流体以定常流的型式流动,称为层流,当流速超过“临界点”时,流线会变成非定向,不规则的流动型式,称为“紊流”(或湍流)。这二种不同的流动型式具有不同的渗流特性。,当渗流速度增大到一定值后,流速与压力梯度关系由线性转变为非线性,即流动型式从线性渗流转变为非线性渗流。达西定律就不适用了。,对于低渗透性致密岩石,在低速渗流时,由于流体与岩石之间存在吸附作用,或在粘土矿物表面形成水膜,当压力梯度很低时,流体不流动,因此存在一个启动压力梯度,a,,在低于该压力梯度范围内流速与压力梯度不呈线性关系,b,a,卡佳霍夫提出的判断指标,雷诺数,式中:,R,e,雷诺数,反映了惯性力与粘性力的比值,也反映了孔隙介质的特点;,流体密度();,1750,单位换算系数,与规定的各物理量的单位有关。,当,R,e,0.2-0.3,时渗流服从达西定律;,当,R,e,0.2-0.3,时,则渗流规律受到破坏,这时渗流速度和压差的关系如下式:,(R,e,),在已知岩石和流体物理参数如岩石孔隙度、渗透率和流体密度、粘度条件下,定义,临界雷诺数一般,=0.2-0.3,C,取决于岩层和流体性质的系数;,n,渗流指数。,n=1,时,为线性渗流,相当于管路水力学中的层流。这时渗流速度很小,即,Re,Red,(临界雷诺数),液体在多孔介质内所产生的惯性力极小,可忽略不计。,1/2,n,1,时,为渗流过渡区。这时渗流速度已相当大(即,Re,Red,),流体在多孔介质中惯性力已明显表示出来,故直线渗流定律已破坏。,n=1/2,时,为渗流的平方区(类似于管路水力学中紊流平方区)。这时渗流速度已很大(,Re,Red,),惯性力也很大。在这以后惯性力的增加与压力的下降又成不变的比例关系。,四、渗透率的因次(量纲)与单位,1.,因次,石油系统渗透率使用的绝对单位制,2.,单位,从因次中可知渗透率的单位是长度的平方,与面积的单位相同。,平方米这个单位太大,通常采用平方微米()来表示。,也很大,最常用的单位是,石油系统渗透率使用的混合单位制,目前在渗透率单位制中还得保留有渗流力学中过去常用的一套混合单位制,混合单位制中达西单位制用达西,(,D),表示,混合单位制中工程单位制不管公制或英制均用毫达西(,m,d,)。,所谓,1,达西(,D,)的含义是:以粘度为,1,厘泊(,cp,)的流体完全饱和于岩石孔隙中,在,1,个大气压的压差下,以层流方式通过截面积为,1,平方厘米,长度为,1,厘米的岩样时,若其流量为,1,立方厘米,/s,,则岩石的渗透率为,1,达西(英文字为,Darcy,简写,D,)。,但这个数值太大,使用起来不方便,因而采用它的千分之一,称为毫达西(英文字母,millidarcy,简写,md,),来作为渗透率的单位。,五、储油岩石渗透率的大小,渗透率是储油(气)岩物性的重要参数,不论对油气运移聚集,还是油(气)田开发都是基础数据。但其数值在不同的油(气)层中相差很大,由,几个毫达西到几千个毫达西,。,我国和世界一些油田的渗透率(据徐献中),级别,孔隙度,(%),渗透率,(md),特高,30,2000,高,25-30,500-2000,中,15-25,100-500,低,10-15,10-100,特低,10,10,储层物性分级,(中国油气储层研究,1994,),1,、滑动效应,克林肯博格效应,2,)同一岩石,同一平均压力,不同气体测得的渗透率不同,气体渗透率与平均压力的关系,实验发现,1,)同一岩石、同一种气体,在不同的平均压力下测的的气体渗透率不同,低平均压力下气体渗透率比较高,高平均压力下气体渗透率比较低,孔道中的液体流动速度,孔道中的气体流动速度,油,气,液体流动,孔道中心的液体分子比孔道壁表面的分子流速高,气体流动时,孔道壁表面的气体分子与孔道中心的分子流速几乎没有差别由于液固间的分子力比液液间的分子力大,在管壁附近表现出的粘滞阻力最大,可使得管壁处液体得流速为零,管道中心粘滞阻力最小,流速最大。,克林博格从分析孔隙内气、液流速分布入手解释了这种现象,气测渗透率时,出了气固间得分子作用力小以外,相邻层的气体分子还可以由于动能交换而使得管壁处的气体分子层与孔道中心的分子层的流速被不同程度均一化。管壁处的气体分子层流速不为零,形成,“,气体滑脱效应,”,。克林博格发现了这种效应,因而也称为,“,克氏效应,”,。,气体的这种滑脱效应还与气体的性质有关,,不同气体的分子量不同,分子直径不同,自由行程也就不同,使得滑脱系数(,b,)不同,分子量大,则,b,大,滑脱效应严重。,同一岩石,不同气体测得的渗透率和平均压力的直线关系交纵座标于一点,该点的气体渗透率与同一岩石的液体渗透率是等价的。故称为等价液体渗透率,也称,克林肯博格渗透率,。,用气体测得的渗透率要比用液体测定的渗透率值高,气体法测定的渗透率更能反映岩石的真实渗透率,平均压力小,气体密度就小,气体间分子间的碰撞就少,使得气体更易流动,气体滑脱现象越严重,所测的渗透率值越大。相反平均压力增大,气体滑脱效应就消失,渗透率就越小。如果压力增大到无穷大,气体的流动性质接近液体的流动性质,气固间的作用力增大,管壁上的气膜逐渐趋于稳定,这时的渗透率趋于一个常数,接近液测渗透率,故称该渗透率为等效液体渗透率或克氏渗透率。,根据滑动效应,岩石本身渗透率的高低也将影响气体渗透率大小,1941,年,克林博格给出了考虑气体滑脱效应的气测渗透率表达式,第三章 储油(气)岩石的渗透率,2,储油(气)岩石渗透率的内在含义,绝对渗透率是岩石本身具有的固有性质,它只与岩石的孔隙结构有关,与通过岩石的流体性质无关。对于这一性质,可以从渗透率与裂缝宽度的关系、渗透率与毛细管孔隙半径的关系来理解。,孔隙度是岩石中孔隙体积与岩石外表体积的比值。这是岩石的内在性质。,一、渗透率与裂缝宽度的关系,对于单一裂隙来说,其长度和截面分别为,L,和,lb,,当压差为,p,时,可由,Poiseulle,方程求出流量,q,1,:,对于,n,个裂缝来说,岩石截面积为,A,,流量为:,根据达西定律,二、渗透率与毛细管孔隙半径的关系,L,r,如果我们把,孔隙,看成是一根根等直径的毛细管束(,n,根,),则可写成:,对于单根毛细管(半径为,r,、长度为,L,),流体在其中流动时的流量,根据伯叶稷(,Poisenille,)定律可求得,即:,如果岩石长度为,L,,岩石截面积为,A,,溶孔孔隙度为,孔,,从孔隙度定义得:,L,r,按达西公式岩石的流量应为:,如果考虑孔隙结构的复杂性,是孔隙结构系数(曲折度),3,储油(气)层,平均,渗透率的计算,一、不均一地层平均渗透率的计算,已知地层总厚度为各层厚度之和,,即:,h=h,1,+h,2,+h,3,地层流体总流量,Q,为各层流量之和,,即:,Q=Q,1,+Q,2,+Q,3,1,纵向不均一,(,1,)平面线性稳定渗流:,设有三个渗透率不等的地层,油气作平面线性渗流,求地层的平均渗透率。,根据达西公式:,h=h,1,+h,2,+h,3,Q=Q,1,+Q,2,+Q,3,加权平均,(,2,)平面径向稳定渗流:,设有三个渗透率不等的地层(,K,1,、,K,2,、,K,3,),油气作平面径向渗透,求地层平均的渗透率。,已知地层总厚度为各层厚度之和,,即:,h=h,1,+h,2,+h,3,地层流体总流量,Q,为各层流量之和,,即:,Q=Q,1,+Q,2,+Q,3,根据达西平面径向稳定渗流公式:,Q=Q,1,+Q,2,+Q,3,h=h,1,+h,2,+h,3,加权平均,2,横向不均一,(,1,)平面线性稳定渗流:,设地层横向不均一造成三个不同渗透率带,K,1,、,K,2,、,K,3,,地层的延伸长度为,L,1,、,L,2,、,L,3,,地层油气为平面线性渗流,求地层总的渗透率。,已知地层的总延伸长度为各段长度之和为,:,L,=,L,1,+,L,2,+,L,3,地层流体的总压力差,p,为各段压差之和为,:,p,=,p,1,+,p,2,+,p,3,Q,=,Q,1,=,Q,2,=,Q,3,,,h,=,h,1,=,h,2,=,h,3,W,=,W,1,=,W,2,=,W,3,根据达西定律总压差和各段压差如下:,Q,=,Q,1,=,Q,2,=,Q,3,,,h,=,h,1,=,h,2,=,h,3,L,=,L,1,+,L,2,+,L,3,调和平均,p,=,p,1,+,p,2,+,p,3,设地层在横向上有三个不同的渗透带,油气为平面径向渗流,求地层平均渗透率。,地层总压差为各段分压差之和,,p,=,p,1,+,p,2,+,p,3,如果井径为,R,w,,从井中心到各带的半径分别为,R,1,、,R,2,、,R,3,,井底压力为,p,w,,最外带压力为,p,e,。,Q,=,Q,1,=,Q,2,=,Q,3,,,h,=,h,1,=,h,2,=,h,3,(,2,)平面径向渗流:,根据达西定律的径向渗流公式可计算总压差和各带压差:,p,=,p,1,+,p,2,+,p,3,调和平均,二、,具有溶孔和裂缝地层总渗透率的计算,1.,具有溶孔的地层,K,总,=K,基,+K,孔,或,计算有溶蚀孔的地层渗透率时,我们可以把它看成由两部分组成。一部分是不含溶蚀孔的基质渗透率(,K,基,),另一部分是含溶蚀孔的岩石渗透率(,K,孔,),总渗透率为:,2.,具有裂缝的地层,对具有裂缝的地层计算其渗透率,也是把它分为两部分。一是不具裂缝的岩石基质的渗透率,一是裂缝的渗透率。,也可用帕尔逊斯方程加以描述,,书,p57,可以看出,该方程与上式基本上是一致的,差别在于:帕尔逊斯方程不仅考虑孔隙渗透率和裂隙渗透率,而且还考虑到了,裂隙的组数和裂隙与流体渗流方向之间的夹角,孔隙(或基质)渗透率,裂隙渗透率,2.,垂直线性稳定渗流,
展开阅读全文