一种基于kozeny-carmen方程改进的渗透率预测新方法-张恒荣.pdf

《一种基于kozeny-carmen方程改进的渗透率预测新方法-张恒荣.pdf》由会员分享，可在线阅读，更多相关《一种基于kozeny-carmen方程改进的渗透率预测新方法-张恒荣.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第 卷第期 年 月 吉林大学学报(地球科学版)J o u r n a l o f J i l i n U n i v e r s i t y ( E a r t h S c i e n c e E d i t i o n ) V o l N o M a y 张恒荣,何胜林,吴进波,等一种基于K o z e n y C a r m e n方程改进的渗透率预测新方法吉林大学学报(地球科学版) , , ( ) : d o i : / j c n k i j j u e s e Z h a n g H e n g r o n g , H e S h e n g l i n , W u J i n b

2、o , e t a l A N e w M e t h o d f o r P r e d i c t i n g P e r m e a b i l i t y B a s e d o n M o d i f i e d K o z e n y C a r m e n E q u a t i o n J o u r n a l o f J i l i n U n i v e r s i t y ( E a r t h S c i e n c e E d i t i o n ) , , ( ) : d o i : / j c n k i j j u e s e 一种基于KozenyCarmen

3、方程改进的渗透率预测新方法张恒荣,何胜林,吴进波,吴一雄,梁玉楠中海石油(中国)有限公司湛江分公司,广东湛江 摘要:在计算复杂孔隙结构储层渗透率时,常规采用的孔渗指数方法或流动单元分类方法几乎很难准确评价渗透率.针对这一问题,本文提出一种引入修正迂曲度因子的改进的K o z e n y C a r m e n方程渗透率计算新方法.首先引入迂曲度因子修正K o z e n y C a r m e n方程,迂曲度因子可以表达为孔隙度与岩电参数的函数;然后对改进的K o z e n y C a r m e n方程进行推演变换,得到新的流动单元指数,能够更好地将储层进行分类;最后利用自适应神经模糊推理

4、系统建立取心段岩心渗透率与测井曲线的模型,并将此模型应用到非取心段的渗透率评价中.岩心渗透率与预测渗透率的对比验证了该方法的正确性与有效性,且渗透率计算精度较常规孔渗指数方法和流动单元分类方法有较大提高.该方法在南海西部海域莺歌海盆地东方气田储层评价中应用效果良好.关键词:渗透率; K o z e n y C a r m e n方程;流动单元指数( F Z I ) ;自适应神经模糊推理系统( A N F I S )doi: / j c n k i j j u e s e 中图分类号:P 文献标志码:A收稿日期: 作者简介:张恒荣( ) ,男,硕士,工程师,主要从事测井方法、解释的科研工作, E

5、 m a i l : z h a n g h r c n o o c c o m c n基金项目:国家科技重大专项项目( Z X )S u p p o r t e d b y N a t i o n a l S c i e n c e a n d T e c h n o l o g y M a j o r P r o j e c t ( Z X )ANewMethodforPredictingPermeabilityBasedonModifiedKozenyCarmenEquationZ h a n g H e n g r o n g , H e S h e n g l i n , W u J

6、 i n b o , W u Y i x i o n g , L i a n g Y u n a nZhanjiangBranchofCNOOCLtd,Zhanjiang ,Guangdong,ChinaAbstract: I n t h e c a l c u l a t i o n o f r e s e r v o i r p e r m e a b i l i t y o f c o m p l e x p o r e s t r u c t u r e , i t i s v e r y d i f f i c u l t t oa c c u r a t e l y e v a l

7、 u a t e t h e p e r m e a b i l i t y T o s o l v e t h i s p r o b l e m , t h i s p a p e r p r e s e n t s a n e w m e t h o d t oi m p r o v e t h e c a l c u l a t i o n o f p e r m e a b i l i t y b a s e d o n t h e K o z e n y C a r m e n e q u a t i o n b y i n t r o d u c i n g t h em o d

8、 i f i e d t o r t u o s i t y f a c t o r F i r s t l y , t h e K o z e n y C a r m e n e q u a t i o n i s m o d i f i e d b y a d d i n g t h e t o r t u o s i t yf a c t o r , w h i c h c a n b e e x p r e s s e d a s a f u n c t i o n o f p o r o s i t y a n d e l e c t r i c a l p a r a m e t

9、e r s ; t h e n t h e m o d i f i e dK o z e n y C a r m e n e q u a t i o n i s d e d u c e d a n d t r a n s f o r m e d t o g e t t h e n e w f l o w u n i t i n d e x f o r b e t t e r r e s e r v o i rc l a s s i f i c a t i o n ; f i n a l l y , a m o d e l o f c o r e p r e m e a b i l i t y

10、a n d l o g c u r v e o f c o r i n g a n d l o g g i n g s e c t i o n u s i n g万方数据a d a p t i v e n e u r a l f u z z y i n f e r e n c e s y s t e m , a n d t h i s m o d e l i s a p p l i e d t o t h e n o n c o r i n g s e c t i o n i np e r m e a b i l i t y e v a l u a t i o n T h e c o r r

11、e c t n e s s a n d e f f e c t i v e n e s s o f t h e m e t h o d i s v e r i f i e d b y c o m p a r i n g t h ep e r m e a b i l i t y o f t h e c o r e a n d t h e p r e d i c t e d p e r m e a b i l i t y , a n d t h e a c c u r a c y o f t h e p e r m e a b i l i t y c a l c u l a t i o ni s

12、m u c h h i g h e r t h a n b o t h t h a t o f t h e c o n v e n t i o n a l p o r o s i t y e x p o n e n t i a l m e t h o d a n d t h e f l o w u n i tc l a s s i f i c a t i o n m e t h o d T h i s m e t h o d h a s g o o d e f f e c t s i n t h e g a s r e s e r v o i r e v a l u a t i o n o f

13、 Y i n g g e h a i b a s i ne a s t i n w e s t e r n S o u t h C h i n a S e a Keywords: p e r m e a b i l i t y ; K o z e n y C a r m e n e q u a t i o n ; f l o w z o n e i n d e x ( F Z I ) ; a d a p t i v e n e u r a l f u z z yi n f e r e n c e s y s t e m ( A N F I S )引言目前,地质学家和油藏工程师最大的挑战便是如何

14、开采地下的剩余石油;精确分析孔隙喉道属性和流体分布是提高采收率的关键;孔喉属性控制原始油、剩余油分布和流体流动.渗透率在油藏描述模型中是一个关键参数,储量估算和产量预测都需要它,而在非均质油藏中要想得到准确的渗透率是一项非常复杂的工作.渗透率数据通常从岩心数据中得来,然而大多数井往往不取心,这时最常采用的方法是在取心层段建立岩心渗透率与孔隙度模型,应用到那些没有取心的层段或者邻井.这种模型主要基于地质统计学,经典孔渗公式为l gKab(K为渗透率,为有效孔隙度,a和b为系数) .这种对数刻度的渗透率孔隙度交会图显然没有理论依据,不仅因为渗透率被刻意绘制为对数函数,而且渗透率与孔隙度之间没有必然

15、的因果关系 .孔隙度通常不依赖于粒度分布,而渗透率是非常依赖于粒度分布的.例如,在一个油藏中,孔隙度和渗透率一般成正比,但有时不同的渗透率可能对应相等的孔隙度值.因此,这种传统方法用来评价渗透率存在不可靠性,特别在复杂岩性储层中.对于常规均匀的孔隙岩石,国外学者K o z e n y 根据毛管理论提出一个公式, C a r m e n 对这个公式进行了改进,即K o z e n y C a r m e n公式. K o z e n y C a r m e n公式是一个理论公式,认为渗透率主要与孔隙迂曲度和孔隙几何形状有关. T i m u r 最早提出了渗透率和束缚水、孔隙度的公式. C o

16、a t e s等 将其改进变形,并认为束缚水可以通过核磁共振测井资料得到, C o a t e s渗透率模型已经成为核磁评价渗透率的常用计算公式. K r u m b e i n等 认为孔喉半径、孔径表面、粒度和分选、颗粒沉积的各向异性对渗透率存在至关重要的影响,提出了渗透率与以上几种岩石特性的关系. H e r r o n 提出了一个考虑矿物成分的渗透率模型,认为矿物成分往往随着沉积粒度和外观的改变而变化,种种要素均会影响岩石的孔径分布,从而导致岩石渗透率的变化,其中关键参数需要通过地球化学元素测井方法得到.渗透率主要考虑个方面的内容:黏土矿物因素、颗粒胶结物因素、孔隙结构因素.胶结矿物含量

17、直接影响孔隙的大小和形态,一般会降低岩石的渗透率和孔隙度.一般碎屑岩储层用K o z e n y C a r m a n方程计算渗透率,综合考虑适用范围较广.随着油气勘探的深入,复杂孔隙结构低渗储层渗透率评价越来越困难.笔者引入迂曲度因子对K o z e n y C a r m e n方程进行修正,通过对方程进行推演变换,得到新的流动单元指数,以期更好地将储层进行分类;并利用自适应神经模糊推理系统建立取心段岩心渗透率与测井曲线的模型,将其推广到非取心段储层中.常规KozenyCarmen模型国内外学者已经从孔隙度和其他相关参数中提出了许多评价渗透率的模型,最早的是K o z e n y模型,其

18、主要关键参数是有效孔隙度、迂曲度和比表面积. K o z e n y模型认为多孔介质可等效为大量相同孔隙半径的迂曲毛管束,再结合达西定律求解泊松方程可得到渗透率:Kr . ( )式中:K为渗透率, m ;为迂曲度;为有效孔隙度;r为毛管束半径, m .方程( )后来被C a r m a n修改为如下比较通用的形式:KfsSg v( ) . ( )式中:fs为形状系数;Sg v为颗粒比表面, m .A m a e f u l e 等利用K o z e n y C a r m e n理论提出 吉林大学学报(地球科学版) 第 卷万方数据了一种基于平均水力半径概念的可以识别和表征地质相的流动单元方法.

19、流动单元被认为是一种有代表性的地层基本单元,它是岩石矿物学特征(如类型、丰度、形态学参数及孔喉的相对位置等)与结构(如颗粒大小、形状、分选性及接触方式等)的函数.而流动单元可以用“流动单元指数” (FZI)来描述和量化:FZIfsSg vRQIz. ( )式中:RQI为储层品质因子,RQI K/;z为归一孔隙度指数,z/ ( ) .对式( )两边取对数,得到l gRQIl gzl gFZI. ( )式( )可以绘制为一种双对数坐标下的交会图,RQI为纵坐标,z为横坐标.在RQI与z双对数绘图中,相等的FZI值对应的岩心数值落在斜率为的线上,其他的FZI值则分布在各自平行的斜线上,而FZI的数值

20、便是斜率方程在z 处的截距.FZI值相等的岩石具有孔径分布类似特征,这样的储层属于同一种流动单元.改进的KozenyCarmen模型K o z e n y C a r m e n方程是在孔喉半径平均值的概念上发展起来的,大量学者研究 认为渗透率与孔隙度、迂曲度的这种相关性对各向同性均匀多孔岩石非常实用且容易定量化,但该方程不能很好地适用于复杂各向异性孔隙岩石.本文提出改进的K o z e n y C a r m e n模型,主要是将模型中的迂曲度进行修改,适用于更宽泛的多孔岩石.岩石迂曲度可以较准确地从岩石的导电性质和有效孔隙度测量中得到.迂曲度定义为流体质点通过多孔介质的旅行长度La与实际多

21、孔介质长度L比值的平方 .它的数学表达形式为(LaL) . ( )与地层因素F有以下理论关系 :(F) . ( )其中,Fam. ( )式中:a是岩性指数;m是胶结指数.因此,迂曲度可以写成(am ) . ( )式( )表明孔隙度与迂曲度成非线性关系,当孔隙为直毛管束时,a和m值等于 ,也等于 .同样,假设为 ,也为 .随着m值增大,迂曲度和孔隙度的关系越来越非线性化(图 ) .图岩石迂曲度与孔隙度、胶结指数的关系FigRelationshipamongrocktortuosity,porosityandcementationindex迂曲度的理论公式已经被实验验证,将公式( )融合到K o

22、z e n y C a r m e n方程中,则改进的K o z e n y C a r m e n模型可以写成KfsaSg vm( ) . ( )改进的K o z e n y C a r m e n方程(式 )引入迂曲度改善了对渗透率的评价.对于直毛管束孔隙,a、m、均等于 ,说明常规K o z e n y C a r m e n模型是改进K o z e n y C a r m e n模型的一种特例.K o z e n y C a r m e n模型在非均质性弱的岩石中应用效果很好,但是实际情况中是非线性的,强烈依赖于m.在自然界岩石中,m值一般为 ,它反映了岩石的物理性质和几何性质,包括胶

23、结、颗粒大小、压实、分选和孔隙类型.因此,将m值引入改进的K o z e n y C a r m e n模型从本质来讲更适合渗透率评价.将式( )进行开根号和单位转换,得到 KfsaSg vm( ) . ( )豪达西( m D )为非法定计量单位, m D m ,下同.式中,K的单位变为m D .式( )左侧是储层品第期 张恒荣,等:一种基于K o z e n y C a r m e n方程改进的渗透率预测新方法万方数据质因子RQI, / (fsaSg v )是改进的流动单元指数FZIm ,同时,z/ ( ) ,则式( )可简写为RQIFZIm z m . ( )两边均取对数可得l gRQIl

24、 gzl gFZIm(m ) l g.( )式( )表明:如果m ,则式( )等效于常规K o z e n y C a r m e n模型;随着m值的增加,RQI与z m在双对数坐标下斜线倾斜越高.每组具有相同FZIm的岩石会构成一个流动单元.W i n s a u e r等 做了一份详细的超过 块砂岩样品的孔隙度、渗透率、地层因素、胶结指数等岩石物理实验.图描述了改进FZI(RQIz m )与常规FZI(RQIz )在双对数坐标下的关系,图中所有实线代表斜率为的斜线,该图表示不同流动单元在常规K o z e n y C a r m e n模型和改进K o z e n y C a r m e

25、n模型中的聚类.为了将岩心数据分为离散流动单元,按照以下公式进行分类 :HFUR o u n d l gFZIC . ( )式中:HFU为流动单元; R o u n d表示取整数;C一般取 .从图中可看到:常规K o z e n y C a r m e n模型分类样品点比较分散;改进K o z e n y C a r m e n模型分类更合理,每一类流动单元涵盖的样品点跨度更平均、代表性更强.这表明改进K o z e n y C a r m e n模型比常规K o z e n y C a r m e n模型更能表征流动单元.Anfis模糊神经网络方法岩心的改进K o z e n y C a r

26、 m e n模型建立后,需要建立岩心段FZI与对应深度段测井曲线的关系,将建立好的模型应用到非取心段,得到准确的渗透率逐点曲线.最常用的做法是多元数理统计方法 ,另外还有F i s h e r线性判别函数法 、 B a y e s概率判别函数法 等,最近几年神经网络系统也被很多学者应用到建模中.本文将采用A n f i s模糊神经推理系统 .它是模糊理论与神经网络相互作用、影响的结晶,包含了模糊理论与神经网络的特点,综合建模、预测、推广、应用等技术于一体. A n f i s是一种类似自适应网络的方法,其参数选择是以调整隶属函数的输入/输出数据来响应所有数据变化的,这种技术被称为神经自适应学习

27、.使用给定的输入/输出数据集,构造一个模糊隶属函数参数的推理系统,对一个反向传播算法或组合使用最小二乘法,然后利用该方法与推理系统联合计算的数据来建模.模糊推理系统输入和输出的规则库包含个模糊逻辑,主要是规则T a k a g i和S u g e n o型.本次研究的靶区为南海西部莺歌海盆地东方区气田,将所有取心样品点的化验分析孔隙度、渗透率、岩石电阻率参数利用新FZI方法计算得到FZI值,然后选取取心段的常规测井曲线如伽马、密度、中子和声波等,建立岩心FZI与对应常规曲线的A n f i s模型,并用此模型进行回判.图为A n f i s模糊神经网络方法与普通B P神经网络方法计算的FZI值

28、对比图,可见A n f i s模糊神经网络方法与岩心数据吻合更好,相比普通B P神经网络方法,此方法更为精确.图改进KozenyCarmen模型(a)与常规KozenyCarmen模型(b)效果对比图FigComparisonofeffectofimprovedKozenyCarmenmodel(a)andconventionalKozenyCarmenmodel(b) 吉林大学学报(地球科学版) 第 卷万方数据应用效果分析基于以上研究,将改进的K o z e n y C a r m e n模型与A n f i s预测方法应用在莺歌海盆地东方区气田进行渗透率评价.该气田岩心资料与测井资料丰富,

29、测井资料主要是斯伦贝谢公司的随钻测井V i s i o n系列.东方区气田主要目的层为黄流组,岩性为海底扇粉砂岩,泥质较重,岩心分析孔隙度为 ,渗透率为 m D ,属于典型的低孔隙度低渗透率储层 .根据该方法,利用在东方区气田所有取心段孔隙度和渗透率数据建立的模型如图所示,此气田利用FZI指数可分为类流动单元(图 a ) ,对应种不同的岩心孔渗幂函数关系(图 b ) .将东方区岩心FZI和A n f i s模糊神经网络模型应用到本气田D F A井所有井段(包括取心段) ,将得到的每类FZI平均值带入如下公式:K FZI mm( ) , ( )于是可得到逐点渗透率.图第道为计算FZI与岩心FZI

30、的对比信息,可见与岩心吻合较好,同时也证明了A n f i s模糊神经网络方法的回判率非常高;第道为岩心渗透率、常规渗透率、改进渗透率的对比信息.从中子、密度曲线发现, x x m段储层孔隙度几乎一致,而岩心显示储层渗透率上下不一致;说明储层上下两段虽然孔隙度基本相同,但喉道还是有较大差异 .流动单元定义是孔喉分布相似的具有同一FZI值的储层,常规K o z e n y C a r m e n模型虽然也能反映孔喉分布,但只适用于较均质条件下的储层类别,而改进的K o z e n y C a r m e n模型对迂曲度进行修改适用于更宽泛的多孔岩石,可见改进的K o z e n y C a r

31、m e n模型更能反映储层的微观结构,计算的渗透率也更准确.图Anfis模糊神经网络(a)与BP神经网络(b)结果对比图FigComparisonofAnfisfuzzyneuralnetworkandBPneuralnetwork图中每种颜色代表一类流动单元.图东方区岩心FZI分类图FigClassificationmapofFZIinDFcore第期 张恒荣,等:一种基于K o z e n y C a r m e n方程改进的渗透率预测新方法万方数据英尺( f t ) 、英寸( i n )为非法定计量单位, f t i n m ,下同.图DFA井FZI计算渗透率与岩心渗透率对比成果图Fig

32、ComparisonofFZIpermeabilityandcorepermeabilityinDFAwell图东方气田改进渗透率与常规渗透率计算结果误差分析FigErroranalysisofimprovedpermeabilityandconventionalpermeabilitycalculationresultsinDFgasfield对莺歌海盆地东方区气田所有的预测渗透率与岩心渗透率进行应用效果对比分析,结果图所示,中间的蓝色实线表示与岩心渗透率完全相同,上下两条红色虚线内的区域表示个数量级的误差范围.由图可以看出,常规K o z e n y C a r m e n模型计算的数据有

33、些已经落在一个数量级误差范围外了,而改进K o z e n y C a r m e n模型计算的渗透率大部分落在半个数量级内,较常规公式精度有较大提高.结论 )基于K o z e n y C a r m e n方程改进的流动单元指数,相比于传统的流动单元指数方法,引入胶结指数能更好地对流动单元进行分类,更适合于非均质储层. 吉林大学学报(地球科学版) 第 卷万方数据 )应用A n f i s模糊神经网络系统建立取心段测井曲线与岩心流动单元指数的模型,该方法比常用的B P神经网络方法更精确. )本方法很好地解决了复杂储层渗透率求取的问题,使测井渗透率的预测精度有所提高,可推广到其他碎屑岩或碳酸盐

34、岩储层.参考文献(References): T i m u r A A n I n v e s t i g a t i o n o f P e r m e a b i l i t y , P o r o s i t ya n d R e s i d u a l W a t e r S a t u r a t i o n R e l a t i o n s h i p s f o rS a n d s t o n e R e s e r v o i r s J T h e L o g A n a l y s t , , ( ) : D u b r u l e O , H a l d o r s

35、e n R H G e o s t a t i s t i c s f o r P e r m e a b i l i t y E s t i m a t i o n J R e s e r v o i r C h a r a c t e r i z a t i o n , , ( ) : A n s e l m e t t i F S , L u t h i S Q u a n t i t a t i v e C h a r a c t e r i z a t i o no f C a r b o n a t e P o r e S y s t e m s b y D i g i t a l

36、 I m a g e A n a l y s i s J A A P G B u l l e t i o n , , ( ) : K o z e n y J U e b e r K a p i l l a r e L e i t u n g d e s W a s s e r s i mB o d e n , S i t z u n g s b e r i c h t e J R o y a l A c a d e m y o fS c i e n c e , V i e n n a , P r o c C l a s s I , , ( ) : C a r m a n P C F l u i

37、 d F l o w T h r o u g h G r a n u l a r B e d s J T r a n s a c t i o n s I n s t i t u t i o n o f C h e m i c a l E n g i n e e r s , , ( ) : T i m u r A P r o d u c i b l e P o r o s i t y a n d P e r m e a b i l i t y o fS a n d s t o n e I n v e s t i g a t e d T h r o u g h N u c l e a r M a

38、g n e t i cR e s o n a n c e P r i n c i p l e s J T h e L o g A n a l y s t , , ( ) : C o a t e s G R , D u m a n o i r J L A N e w A p p r o a c h t oI m p r o v e d L o g d e r i v e d P e r m e a b i l i t y J T h e L o gA n a l y s t , , ( ) : K r u m b e i n W C , M o n k G D P e r m e a b i l

39、 i t y a s a F u n c t i o no f t h e S i z e P a r a m e t e r s o f U n c o n s o l i d a t e d S a n d J T r a n s a c t i o n s o f t h e A i m e , , ( ) : H e r r o n M M E s t i m a t i o n g t h e I n t r i n s i c P e r m e a b i l i t y o fC l a s s i c S e d i m e n t s f r o m G e o c h e

40、 m i c a l D a t a C / /S P W L A t h A n n u a l S y m p o s i u m L o n d o n : s n , : A m a e f u l e J O , A l t u n b a y M E n h a n c e d R e s e r v o i r D e s c r i p t i o n : U s i n g C o r e a n d L o g D a t a t o I d e n t i f yH y d r a u l i c ( F l o w ) U n i t s a n d P r e d i

41、 c t P e r m e a b i l i t y i nU n c o r e d I n t e r v a l s / W e l l s C / / A n n u a l T e c h n i c a lC o n f e r e n c e a n d E x h i b i t i o n o f t h e S o c i e t y o f P e t r o l e u mE n g i n e e r s H o u s t o n : S P E , : A h m e d U , C r a r y S F , C o a t e s G R P e r m

42、e a b i l i t yE s t i m a t i o n : T h e V a r i o u s S o u r c e s a n d T h e i rI n t e r r e l a t i o n s h i p s J I n t e r n a t i o n a l J o u r n a l o f R o c kM e c h a n i c s & M i n i n g S c i e n c e & G e o m e c h a n i c sA b s t r a c t s , , ( ) : B a b a d a g l i T , A l

43、S a l m i S A R e v i e w o f P e r m e a b i l i t yP p r e d i c t i o n M e t h o d s f o r C a r b o n a t e R e s e r v o i r s U s i n gW e l l L o g D a t a J S p e R e s e r v o i r E v a l u a t i o n &E n g i n e e r i n g , , ( ) : F r a n c i s c o J V P , A l b e r t o J O T V a l i d i

44、 t y o f t h eP e r m e a b i l i t y C a r m a n K o z e n y E q u a t i o n : A V o l u m eA v e r a g i n g A p p r o a c h J P h y s , , A : R o s e W , B r u c e W A E v a l u a t i o n o f C a p i l l a r y C h a r a c t e r i n P e t r o l e u m R e s e r v o i r R o c k J J o u r n a l o fP e t r o l e u m T e c h n o l o g y , , ( ) :