应用多参数划分层内流动单元
国内外流动单元划分方法
１９ ９４正 １９ ９５正
冯 晓宏 姚 光庆
特征角度看 ，岩石物理相就是 “ 水力单元” １９ 正 ９５ １９ 正 ９６ 焦养泉、李思 田 裘亦楠 沉积体系内部按水动力条件进一步划分的建筑块体 由于储层的非均质性、隔挡和窜流旁通条件 ，注入水沿
着地 质 结构 引起 的一 定途径 驱 油 、 自然 形成 的流体 流动
法
半径 、 流动 带 尼 一 曼孔 渗 关 系公 Ａｂ ａａｅ 断块 卡 ｂｚｄｎ 指标 式 划 分 取 心 关 键 井 的 魏 斌 １ 阿拉 斯加 北斜 坡恩 迪
岩 心 、 传导系数、 储 流动单元， 非取心井 科特油 田储层 测 井 资 存系数 、 净砂 流动单元的划分， 采 Ｇ ａｇ ｉ ２渤海湾高尚堡油 田 ｕｎｍ ｎ ｇ
特征 角度研究流动单元， 出岩石物理相的概念及研究方法；冯晓宏等根据 提 取心井资料， 油藏角度研 究流动单元的特征； 从 焦养泉等研 究流动单元的边
界和规模， 将其归纳为建筑结构 的一部分； 裘亦楠和穆 龙新等从成 因角度探 讨流动单元 的特征； 阎长辉等在岩石物理性质基础上， 考虑流体性质 的动态
段 为
孔 隙 组 合 类
１美 国得 克萨斯 海 相碳 酸 岩储 层 ２辽 河 油区西 部 凹 陷 ２ ６
主 的 孔 隙 几 型 、 孔 喉 半 理 论基础 ： 平均 水 力 Ｄａｉｓ ｖｅ 储 层 何 学 划 径 、 均流 动 半 径 的概 念 和 柯 兹 Ｖ ｓｅｌ 平 ｅｓｌ 参 数 分 分 析
岩 相 及 岩石 类 型 、 沉 研 究 沉积 学 ， 垂 向 在 宏 观 岩 积构 造 、 粒度 上 划 分 为 若 干 成 因
石 物 理 中值 、 泥质含 单元， 根据孔渗等物 Ｈ ａ ｅｒ ｎ 参 数 划 量 、 砂 层 厚 参 进 一 步 细 分 成 因 Ｈａ ｎ ｍｉ
《油气田地下地质学》复习思考题(15级资工)
《油气田地下地质学》提纲第一章1、名词解释：地质井、参数井、预探井、评价井、开发井、调整井、定向井、丛式井、水平井、CT值、井斜角、井斜方位角、井号编排、钻时、钻时录井、岩心、岩心收获率、岩心编号、岩屑、岩屑录井、岩屑迟到时间、捞砂时间、钻达时间、套管程序、方入、进尺、补心高、补心海拔。
2、录井方法一般包括哪几种？3、影响钻时的因素有哪些？如何根据钻时来判断岩性？4、现场上常用的荧光录井方法有哪几种？5、如何划分碎屑岩的含油级别？6、为什么要进行岩心归位？简述岩心归位的原则和步骤。
7、如何获取有代表性的岩屑？常用的测定迟到时间的方法有哪几种？8、在钻井中泥浆的功能是什么？泥浆性能包括有哪些？9、什么是泥浆的失水量和泥饼？钻井过程中对其作何要求？为什么？10、如何根据井号编排判断井别：渔浅1井、荆参2井、浩4 -3井、陵1井、沙36井。
11、泥浆显示分为哪几类？12、完井方法因地质条件不同可分为哪几类？13、如何根据泥浆性能的变化来判断油、气、水层和其它特殊岩层？14、通过岩心录井及岩心分析可获得哪些资料及信息？第二章1、概念：油气水的综合判断、束缚水、可动水、含油饱和度、相渗透率、增阻侵入、减阻侵入、地层测试、中途测试、跨隔测试、测试垫。
2、在进行油气水层的判断时，为什么对低渗透性砂岩油气层的含油性解释偏低？3、在进行油气水层的判断时，为什么对高渗透性砂岩油气层的含油性解释偏高？4、简述在碳酸盐岩双重孔隙结构中，基质孔隙系统和裂缝系统的主要区别。
5、说明钻柱测试压力卡片上不同压力段测试阀、旁通阀、封隔器所处的工作状态，标注压力卡片上各点所表示的压力。
6、满足什么样条件的压力卡片才能供我们解释分析用？7、对几张典型的压力卡片进行初步分析。
8、简述低阻油层的成因。
9、简述水淹层的地质特征。
第三章1、名词解释：有效厚度、沉积旋回、细分沉积相、标准层、标准剖面、含油层系、油层组、砂岩组、单油层、测井相、地震相。
储层流动单元划分方法与应用
储层流动单元划分及其对产能动态的影响———以鄂尔多斯盆地元 54区长 1油层组为例宋广寿 1 , 2,孙 卫 1,高 辉 1,莫建武 1( 1. 西北大学 地质学系 /大陆动力学国家重点实验室 ,陕西 西安710069; 2. 中国石油长庆油田公司 ,陕西 西安 710021) 摘要 :目的 对有效提高产能和注水开发效果的措施进行分析 ,提出合理化建议 。
方法 通过沉积微相 、储层综合定量分析 、产层生产动态资料进行综合评价 。
结果 本次划分 ,选取了有效厚度 、孔 隙度 、渗透率 、突进系数 、流动带指数等参数 ,以元城油田元 54 区长 1 油层组为例 ,将储层划分为E 、G 、M 、P 4类流动单元 。
结论 各类流动单元与储层沉积微相 、岩性 、物性和油层初始产能具有很好的对应关系 ,能够真实客观地反映特低渗砂岩储层物性差 、非均质性强的地质特点 。
关 键 词 :鄂尔多斯盆地 ;特低渗透砂岩储层 ;流动单元划分 ;产能动态评价 中图分类号 : TE33 + 1 文献标识码 : A 文章编号 : 1000 2274 Ⅹ ( 2010 ) 02 20299 205The cl a s si f i c a t i o n of fl ow u n i ts and i ts i n fluence o n deli v er abili ty performance : A ca s e of the Ch ang 1 forma t i o n of Y u an 54 Area s i n O r dos Ba si nS ON G Guang 2shou 1 , 2 , S UN W e i 1 , G AO H u i 1 , M O J ian 2wu 1( 1. D e p a r tm e n t of G eo l o g y / The Sta t e Key L a b o r a t o r y of C on t inen t a l D y nam ic s , No r thwe st U n i ve r sity, X i ′an 710069 , Ch i na;2. Chan g q ing O ilfie ld Comp a ny, Pe t r och i na, X i ′an 710021, Ch i na )T o ana l yze and p r opo s e the ra t i o na l m e a s u r e s of enhanc i ng the p r oduc t i o n and fl o od i ng deve l o p 2A b s tra c t : A i m m e n t effec t effec t i ve l y i n a l l k i nd s of fl o w un its . M e thod sB a s ed on comp rehen s i ve quan t ita t i ve ana l ysis and eva l u 2a ti o n of sed i m en ta ry m i c r ofac i e s and re se rvo ir . Re su lts Tak i ng Yuan 54 A rea s Chang 1 f o r m a ti o n i n Y uancheng o i lfi e l d a s an examp le, fi ve p a ram e te rs, i nc l ud i ng effi c i ency th ickne ss, po r o sity, p e r m eab ility, da sh fac t o r and fl o w 2z one i ndex a re se l ec ted f o r . The f o r m a ti o n is fi ne l y d i vi ded i n t o 4 typ e s of fl o w un its, i nc l ud i ng E, G, M , P . C o n 2c lu s i on The re su lts can be ve rifi ed w ith p r oduc ti o n behavi o r da ta each o the r and show tha t the se l ec ted va l u e s ofp a ram e te rs, comp rehen si ve eva l ua ti o n func ti o n s and i ndex e s f o r the fl o w un it d i vis i o n can tru l y and ob j ec t i v e l y r e 2flec t the geo l o gi c fea tu re s of u ltra 2l o w p e r m eab ility re se rvo ir, such a s poo r p r op e rty cha rac te risti c s and str o ng he t e r 2ogene i ty .Key word s :O rdo s B a s i n; u ltra 2l o w p e r m e ab i lity sand re s e r vo i r; fl o w un its subd i visi o n; de l i ve r ab i lity p e r f o r m a nce eva l ua t i o n1319 % ,平均空气渗透率 7111 ×10 - 3μm 2,孔隙类型以粒间孔和溶蚀孔为主 ,属成岩性低孔 、特低渗透油 藏 。
浅水三角洲储层流动单元划分
[收稿日期]20220223[基金项目]国家重大科技专项 渤海油田加密调整及提高采收率油藏工程技术示范 (2016Z X 05058001)㊂ [第一作者]吴穹螈(1990),男,硕士,工程师,现主要从事油气田开发地质研究,w u _q y90@126.c o m ㊂ 吴穹螈,王少鹏,张岚,等.浅水三角洲储层流动单元划分[J ].长江大学学报(自然科学版),2023,20(2):84-95.WU Q Y ,WA N GSP ,Z HA N G L ,e t a l .C l a s s i f i c a t i o no f f l o wu n i t s i ns h a l l o w w a t e r d e l t a r e s e r v o i r [J ].J o u r n a l o fY a n g t z eU n i v e r s i t y(N a t u r a l S c i e n c eE d i t i o n ),2023,20(2):84-95.浅水三角洲储层流动单元划分吴穹螈,王少鹏,张岚,翟上奇,常会江中海石油(中国)有限公司天津分公司,天津300452[摘要]为准确反映储层内部非均质性,以指导油田后期开发调整,有必要对砂体开展流动单元研究㊂综合应用测井㊁地震㊁分析化验和动态资料,采用熵权T o p s i s 法对渤海B 油田开展浅水三角洲储层流动单元划分㊂研究表明,研究区浅水三角洲储层非均质性较强,孔隙度㊁渗透率㊁泥质含量和饱和度中值压力是影响储层储集渗流能力的主要因素,根据以上4类参数可将储层划分为Ⅰ㊁Ⅱ㊁Ⅲ㊁Ⅳ四类流动单元㊂流动单元的分布受到沉积相的控制,Ⅰ㊁Ⅱ类流动单元主要分布在分流河道和分流砂坝中,连通性较好;Ⅲ㊁Ⅳ类流动单元主要分布在分流河道边部和溢岸砂中,连通性较差㊂流动单元的组合样式对连通性具有明显控制作用,不同流动单元间形成渗流差异界面,使得砂体内部弱连通或不连通,导致剩余油大量富集㊂[关键词]浅水三角洲;流动单元;熵权;T o p s i s 法;剩余油[中图分类号]T E 321,T E 122.2[文献标志码]A [文章编号]16731409(2023)02008412C l a s s i f i c a t i o no f f l o wu n i t s i n s h a l l o ww a t e r d e l t a r e s e r v o i rWU Q i o n g y u a n ,WA N GS h a o p e n g ,Z HA N GL a n ,Z HA I S h a n g q i ,C HA N G H u i j i a n gT i a n j i nB r a n c h ,C N O O CL i m i t e d (C h i n a ),T i a n ji n300452A b s t r a c t :I no r d e r t oa c c u r a t e l y r e v e a l t h e i n t e r n a lh e t e r o g e n e i t y o f r e s e r v o i r sa n d g u i d et h e l a t e rd e v e l o pm e n ta n d a d j u s t m e n t o f o i l f i e l d s ,i t i s n e c e s s a r y t o c a r r y o u t t h e r e s e a r c h o n t h e f l o wu n i t s o f s a n db o d y .B a s e d o n t h e c o m p r e h e n s i v e a p p l i c a t i o n o f l o g g i n g ,s e i s m i c ,a n a l y t i c a l t e s t a n dd y n a m i cd a t a ,t h e e n t r o p y w e i g h tT o p s i s m e t h o dw a su s e d t oc l a s s i f yt h e f l o wu n i t s o f s h a l l o ww a t e r d e l t a r e s e r v o i r i nB o h a iO i l f i e l dB .T h e s t u d y s h o w s t h a t t h e s h a l l o ww a t e r d e l t a r e s e r v o i r i n t h e s t u d y a r e a h a s s t r o n g h e t e r o g e n e i t y ,a n d p o r o s i t y ,p e r m e a b i l i t y ,s h a l e c o n t e n t a n d s a t u r a t i o nm e d i a n p r e s s u r e a r e t h em a i n f a c t o r s a f f e c t i n g t h e r e s e r v o i r s e e p a g e c a p a c i t y .A c c o r d i n g t o t h e a b o v e f o u r p a r a m e t e r s ,t h e r e s e r v o i r c a nb e d i v i d e d i n t o f o u r t y p e s o f f l o wu n i t s :Ⅰ,Ⅱ,Ⅲa n d Ⅳ.T h ed i s t r i b u t i o no f f l o wu n i t s i s c o n t r o l l e db y s e d i m e n t a r yf a c i e s ;t y p e I a n d Ⅱf l o wu n i t s a r em a i n l y d i s t r i b u t e d i nd i s t r i b u t a r y c h a n n e l s a n dd i s t r i b u t a r y s a n db a r s ,w i t hg o o d c o n n e c t i v i t y ,whi l e t y p eⅢa n dⅣf l o wu n i t s a r em a i n l y d i s t r i b u t e d i n t h e e d g e o f d i s t r i b u t a r y c h a n n e l s a n d o v e r b a n k s a n d ,w i t h p o o rc o n n e c t i v i t y .T h ec o m b i n a t i o n p a t t e r no f f l o w u n i t sh a sas i g n i f i c a n tc o n t r o lo nc o n n e c t i v i t y ,a n d d i f f e r e n t f l o wu n i t s f o r ms e e p a g e d i f f e r e n t i a l i n t e r f a c e s ,w h i c h m a k e s t h e s a n db o d y w e a k l y c o n n e c t e do r d i s c o n n e c t e d ,r e s u l t i n g i na l a r g e a m o u n t o f r e m a i n i n g oi l e n r i c h m e n t .K e yw o r d s :s h a l l o w w a t e r d e l t a ;f l o wu n i t ;e n t r o p h y w e i g h t ;T o p s i sm e t h o d ;r e m a i n i n g o i l ㊃48㊃长江大学学报(自然科学版) 2023年第20卷第2期J o u r n a l o fY a n g t z eU n i v e r s i t y (N a t u r a l S c i e n c eE d i t i o n ) 2023,V o l .20N o .2随着油田开发阶段的深入,油气藏生产情况也呈现出明显的变化,多表现为油气产量递减大㊁含水上升快以及注采不平衡等特征㊂因此,为准确反映储层内部非均质性,提高油藏精细描述精度,以指导油田后期开发调整,有必要将油藏描述的尺度从小层㊁单层细化到单砂体内部流动单元㊂流动单元为纵横向连续的孔隙度㊁渗透率以及层理特征相似的储集体,多用于表征油藏非均质性,流动单元的划分实际上就是对储层质量进行分类[1-4]㊂近年来,随着油藏精细描述在油田开发生产中的广泛应用,流动单元对储层连通性及剩余油的控制作用逐渐被国内外学者所关注[5-7],国内外学者对流动单元开展了大量研究,认为流动单元的研究是非均质性研究的延伸,流动单元内部渗流特征相似并随着开发呈现出动态变化[6];不同流动单元内部渗流特征差异对注入水波极具有重要影响作用[7];流动单元分类是流动单元研究的核心内容,并提出了流动带指数法㊁聚类分析法㊁综合参数法㊁孔隙结构参数法㊁生产动态参数法等多种流动单元的定性及定量划分方法[8-10]㊂然而,目前被广泛使用的流动单元定量划分方法都需要对评判指标进行权重的确定,而权重的确定难以避免主观性或繁琐的计算量,缺乏一种简洁且能客观对评判指标权重赋值的识别方法[11]㊂随着近年来智能化理念的深入,信息熵逐渐被引入到石油领域研究中,信息熵是系统衡量信息不确定性的指标,其对数据的非均一性具有良好的反映,能够直观体现出评判参数的差异性㊂然而采用信息熵进行判别仍然面临评判参数权重难以客观赋予的问题㊂探寻一种客观简洁评判流动单元的方法,对于明确储层内部非均质性差异及连通性具有重大指导意义㊂而T o p s i s 法是一种逼近于理想解的排序方法,能够高效简洁地对不同评判参数权重进行客观赋值㊂将信息熵和T o p s i s 法相结合进行流动单元研究能够客观高效地筛选出主控因素,并对储层内部非均质性进行分析㊂因此,以渤海B 油田为例,针对储层内部非均质性差异认识难度大的问题,将熵权与T o ps i s 法相结合,对研究区储层进行流动单元划分,进而指导油田后续调整方案的制定㊂1 研究区概况B 油田位于渤海南部海域黄河口凹陷中南部,北靠渤南低凸起,南依莱北低凸起,西傍渤中33-1次断隆,东抵黄河口凹陷东部次洼(见图1)㊂该油田为一继承性发育并被断层复杂化的断块构造,受南北两组掉向相反的北东向正断层控制,断裂系统发育㊂主力含油层系为明化镇组下段(以下简称明下段),油藏类型为岩性构造油气藏,平均孔隙度31.1%,平均渗透率1455m D ,为典型高孔特高渗储层㊂图1 研究区构造位置及地层特征F i g .1 S t r u c t u r a l l o c a t i o na n d s t r a t i g r a p h i c c h a r a c t e r i s t i c s o f t h e s t u d y ar e a ㊃58㊃第20卷第2期吴穹螈等:浅水三角洲储层流动单元划分明下段发育浅水三角洲前缘沉积,以分流砂坝和分流河道为主(见图2),砂体连续性极好,平面均方根属性连片,高值区呈分区坨状分布㊂分流砂坝岩性以细砂岩为主,可见中砂岩,分选磨圆较好,泥质含量较少,砂体厚度多在2.5~7.0m ,垂向韵律以均质韵律和反韵律为主(见图2(a)),分流砂坝测井曲线形态主要为中高幅漏斗形㊁箱形㊂粒度概率曲线为以跳跃和悬浮组分为主的多段式(见图2(b )),岩心可见波状层理㊁块状层理㊁槽状交错层理等沉积构造(见图2(c )㊁图2(d )㊁图2(e))㊂这些特征反映了研究区储层受到波浪改造的牵引流特点,与河口环境处河流㊁波浪交替的水动力条件相匹配㊂该类型浅水三角洲沉积骨架为分流砂坝,分流河道体系主要作为砂质供给通道,这是因为浅水湖盆频繁的长距离岸线迁移以及低坡度河口区摩擦力占主导,使得沉积主体部位在空间上快速迁移,早期分流河道发生迅速废弃与改道[12-14]㊂B 油田于2009年初投产,钻井100余口,目前井距150~350m ,综合含水率已达82%,但采出程度仅为24%,剩余油仍有较大潜力㊂为了给油田提供剩余油挖潜依据,迫切需要开展储层流动单元识别与划分㊂图2 浅水三角洲沉积特征F i g .2 S e d i m e n t a r y ch a r a c t e r i s t i c s o f s h a l l o ww a t e r d e l t a 2 熵权T o ps i s 法原理笔者采用基于熵权的T o p s i s 法对研究区开展储层流动单元的识别与划分㊂T o ps i s 法是根据有限个评价对象与理想化目标的接近程度进行排序的方法,是对现有的对象进行相对优劣的评价[15]㊂其基本原理是通过检测评价对象与正理想解㊁负理想解的距离来进行排序,若评价对象最靠近正理想解同时又最远离负理想解,则为最好㊂T o p s i s 法包括对评判参数构建归一化初始矩阵㊁确定权重㊁确定决策矩阵㊁确定正理想解和负理想解㊁确定贴近程度5个步骤㊂2.1 构建归一化初始矩阵假设共有m 个待评价样本,每个对象都有n 个指标,则原始数据矩阵X 为: X =x 11 x 1j x 1n ︙︙x i 1x i jx i n ︙︙x m 1 x mjx m n éëêêêêêêêùûúúúúúúú(1)㊃68㊃长江大学学报(自然科学版)2023年3月式中:m 为待评价样本的数量;n 为指标的数量;x i j为第i 个样本的第j 个指标的数值㊂由于不同参与评判的指标,其单位量纲不尽相同,导致类似孔隙度和渗透率等指标绝对数量值差异较大,因此需要先对评判指标进行归一化:Z i j =x i jmi =1xi j(2)归一化后得到初始矩阵:Z =z 11 z 1j z 1n ︙︙z i 1z i j z i n ︙︙z m 1 z mjz m néëêêêêêêêùûúúúúúúú(3)2.2 确定评判指标权重评判指标权重的确定是流动单元划分的关键,其直接决定了评判模型以及划分结果㊂研究通过引入熵 的概念来对权重进行客观定量赋值㊂熵反映的是信息不确定性的情况,即衡量一个系统里不确定性的程度㊂熵越大,表明该系统或者指标分布越分散,其不确定性越强;熵越小,表明其不确定性越弱㊂在权重分析中,若一个评判指标其熵值越大,则说明该指标权重度越大㊂根据熵值e j的定义: e j =-ðmi =1z i j l n z i j (4)计算指标差异度h j : h j=1-e jl n m (5)在参与评判的n 个指标中,第j 项指标的权重为:w j =h jðni =1hj(6)2.3 确定决策矩阵决策矩阵V 为归一化后的指标值乘以其对应的权重值:V =w 1z 11 w n z 1n ︙︙w 1z m 1w n z m n éëêêêêùûúúúú=v 11 v 1j v 1n ︙︙v i 1v i jv i n ︙︙v m 1v mjv m néëêêêêêêêùûúúúúúúú(7)式中:v i j为第i 个样本的第j 个指标的权重与归一化值的乘积㊂2.4 确定正理想解和负理想解正理想解代表每个参与评判的指标最理想的值,反之负理想解为其最不理想的值㊂在求解的过程中,需要注意评判参数数值的物理意义与正负理想解的对应关系;比如孔隙度和渗透率等效益型参数其值越大,则代表储层质量越好,而泥质含量等成本型参数值越大则代表储层质量越差㊂V +=(m a x iv i j |j ɪJ 1),(m i n v i j |j ɪJ 2),i =1,2, ,m {}(8) V -=(m i n iv i j |j ɪJ 1),(m a x v i j |j ɪJ 2),i =1,2, ,m {}(9)式中:J 1为效益型参数;J 2为成本型参数㊂㊃78㊃第20卷第2期吴穹螈等:浅水三角洲储层流动单元划分2.5 确定贴近程度计算各评价对象与正㊁负理想解之间的距离㊂距离正理想解的距离越近,相对贴近程度越近,则流动单元的储层质量越好㊂样本与正理想解之间的距离: d +i=nj =1(v i j-v +j )2i =1,2, ,m (10)样本与负理想解之间的距离: d -i=nj =1(v i j -v -j )2i =1,2, ,m (11)相对贴近程度D i 为:D i =d +id +i +d -i(12) 图3 吸水强度与孔喉半径㊁退汞效率㊁饱和度中值压力拟合关系F i g .3 T h e f i t t i n g r e l a t i o n s h i p b e t w e e nw a t e r a b s o r pt i o n i n t e n s i t y a n d p o r e t h r o a t r a d i u s ,m e r c u r yr e m o v a l e f f i c i e n c y an d s a t u r a t i o nm e d i a n p r e s s u r e 其中,0ɤD i ɤ1,D i 越无限逼近1,表明评价对象越优㊂3 流动单元划分泛连通体内部的渗流差异是由其储层质量的差异而导致,渗流差异的判别实际上就是对相对连续的渗流能力的截断分类[16-17]㊂研究应用熵权T o ps i s 法,利用优选的参数划分研究区流动单元,明确其渗流差异界面㊂3.1 参数优选流动单元的分布受到沉积作用㊁成岩作用㊁岩石宏观物性及微观孔隙结构等参数的控制[18-21]㊂因此在选取储层流动单元划分参数时,应当遵循全面且适当的原则,充分考虑沉积成岩㊁宏观微观㊁动态静态等因素;同时还要结合研究区的地质油藏特征以及资料的丰富程度,确保选取的参数能够客观真实反映储层渗流能力,并且在研究区易于获取,具有良好的可复制性㊂流动单元的分布与沉积环境具有密切的关系,泥质含量和粒度中值是最常用的反映沉积环境的参数㊂研究区为典型浅水三角洲前缘沉积,目前仅几口探井和少量开发井进行了取心分析化验,由于粒度中值数据较少,难以获取,考虑资料丰度,选取泥质含量作为反映沉积环境的参数进行流动单元的划分㊂储层的储集和渗流能力对流动单元的划分起着决定性的作用,同一沉积相带内部,其垂向和侧向上的物性差异是明显的,而孔隙度㊁渗透率作为最具有代表性的参数,其具有良好的物性指示作用且易于获取,因此选择孔隙度和渗透率作为反映物性条件的参数来进行评判㊂除了宏观沉积环境㊁成㊃88㊃长江大学学报(自然科学版)2023年3月岩条件以及物性外,岩石微观孔隙结构也是影响流体运移的重要因素㊂孔喉半径㊁退汞效率㊁饱和度中值压力是表征岩石微观孔隙结构最常用的参数,其能够直接反映储层的连通性㊂研究区孔喉半径㊁退汞效率和饱和度中值压力分析化验资料较少,孔喉半径和退汞效率与渗透率和注水井吸水强度相关性较差,而饱和度中值压力与渗透率和注水井吸水强度相关性较好(见图3),表明饱和度中值压力反映了岩石微观孔隙结构的差异,进而影响了岩石的渗流能力,因此选用饱和度中值压力作为反映岩石微观孔隙机构的参数来进行流动单元的划分㊂3.2划分结果图4 饱和度中值压力与渗透率拟合关系F i g .4 T h e f i t t i n g r e l a t i o n s h i p be t w e e n s a t u r a t i o n m e d i a n p r e s s u r e a n d p e r m e a b i l i t y通过参数的优选,研究筛选出孔隙度㊁渗透率㊁泥质含量和饱和度中值压力4个参数进行流动单元的划分㊂在划分之前需要对参数的类型进行划分以便于后续正负理想解的评判,4个参数中孔隙度和渗透率为效益型参数,参数的数值越大,代表储层的质量和渗流能力越高;泥质含量和饱和度中值压力为成本型参数,参数的数值越大,代表储层的质量和渗流能力越差㊂因此,在求取正负理想解的时候,需要将参数的最大最小值与其物理意义进行匹配㊂以取心井为例,运用熵权T o ps i s 法对储层流动单元进行划分㊂取心井的孔隙度㊁渗透率和泥质含量参数可以通过测井曲线获取㊂研究表明,取心井分析化验所得到的饱和度中值压力数据和渗透率具有良好的相关性(见图4),可通过式(13)进行获取: P C 50=6.8339K -0.589(13)式中:P C 50为饱和度中值压力;K 为渗透率㊂评判中,选取样本为工区120口井的数据,共计768000个样本点,开展流动单元划分㊂由于样本数据量较大,下文以部分数据为例,选取7口井的4个评判参数进行划分,故取m =7,n =4,由式(1)得到初始矩阵X 为:0.22746.3080.2310.7410.22028.6500.2860.9470.24266.3790.2650.5770.24990.8160.2350.4800.25099.9750.2170.4510.24995.6990.2030.4390.23353.7410.2430.654éëêêêêêêêêêêùûúúúúúúúúúú对上述矩阵进行归一化处理,消除量纲的差异,得到标准化后的矩阵,并按照式(4)~(7)进行计算得到决策矩阵V :9.3ˑ10-79.3ˑ10-72.6ˑ10-53.6ˑ10-59.0ˑ10-76.2ˑ10-73.2ˑ10-54.9ˑ10-510.0ˑ10-714.4ˑ10-73.0ˑ10-53.0ˑ10-510.2ˑ10-719.8ˑ10-72.6ˑ10-52.5ˑ10-510.2ˑ10-722.0ˑ10-72.4ˑ10-52.3ˑ10-510.2ˑ10-723.0ˑ10-72.3ˑ10-52.3ˑ10-59.6ˑ10-711.7ˑ10-72.7ˑ10-53.4ˑ10-5éëêêêêêêêêêêêùûúúúúúúúúúúú在决策矩阵获取的基础上,根据式(8)~(9)进行计算得到理想解集合:㊃98㊃第20卷第2期吴穹螈等:浅水三角洲储层流动单元划分d +i =7.9ˑ10-64.1ˑ10-68.7ˑ10-32.6ˑ10-3[] d -i =7.9ˑ10-64.1ˑ10-68.7ˑ10-32.6ˑ10-3[] 图5 砂泥岩波阻抗与自然伽马关系F i g .5 R e l a t i o n s h i p b e t w e e nw a v e i m pe d a n c e a n d n a t u r a l ga m m a i n s a n da n dm u d s t o n e 在理想解集合求取的基础上计算指标相对贴近程度D i 的分布㊂研究中,非取心井的参数获取与取心井一致,主要利用测井资料进行计算获取,而井间的参数则主要依靠地震资料进行预测㊂研究区明下段为一套浅水三角洲沉积,整体上具有较低的砂地比㊂三维地震资料同相轴连续性好,信噪比与分辨率较高(主频约45H z),砂岩速度为2600m /s ,密度为2100k g /m 3,泥岩速度为3100m /s ,密度为2350k g/m 3,砂泥岩波阻抗差异明显(见图5),资料品质较好,整体上看,研究区的地震资料能够有效定量表征储层的非均质性㊂开发井实钻后储层参数(孔隙度㊁渗透率㊁泥质含量)和地震振幅参数具有良好的相关性,因此井间的 图6 D i 正态分布曲线F i g.6 N o r m a l d i s t r i b u t i o n c u r v e o f D i 评判参数则主要利用地震振幅属性构建孔隙度场㊁渗透率场和泥质含量场,进而再根据分析化验得到的渗透率和饱和度中值压力的相关性参数构建饱和度中值压力场,来实现参数的获取㊂在D i 分布获取的基础上,根据数据构形理论对储层流动单元完成划分,其原理为储层渗流能力近似的样本其正态分布曲线图上斜率近似相同(见图6),据此将研究区储层流动单元划分为Ⅰ㊁Ⅱ㊁Ⅲ㊁Ⅳ四类,单井流动单元划分结果如图7所示,Ⅰ类流动单元D i 小于0.3,Ⅱ类流动单元D i 介于0.3~0.5,Ⅲ类流动单元D i 介于0.5~0.7,Ⅳ类流动单元D i 大于0.7㊂4 流动单元特征流动单元的划分以平面上泛连通的砂体为约束,根据邻井和无井区域的划分结果将同类流动单元划分为一类,最终得到流动单元的平面分布,进而分析泛连通体内部的侧向连通性㊂4.1 平面分布特征Ⅰ类流动单元储层质量最好,是研究区储集和渗流能力最好的储层㊂研究区Ⅰ类流动单元占比15.4%㊂该类流动单元指标逼近程度小于0.3,砂体厚度大于10m ,孔隙度大于30%,渗透率大于1800m D ,泥质含量小于8%,饱和度中值压力小于0.03M P a (见图8),孔隙结构以大孔粗喉为主(见图9(a ))㊂Ⅱ类流动单元储层质量较好,研究区Ⅱ类流动单元占比32.4%㊂该类流动单元砂体厚度介于6~10m ,指标逼近程度介于0.3~0.5,孔隙度介于26%~30%,渗透率介于1200~1800m D ,泥质含量介于8%~14%,饱和度中值压力介于0.03~0.08M P a ,孔隙结构以大中孔粗中喉为主(见图9(b ))㊂Ⅲ类流动单元储层质量中等,研究区Ⅲ类流动单元最为发育,占比38.6%,该类流动单元砂体厚度介于3~6m ,指标逼近程度介于0.5~07,孔隙度介于20%~26%,渗透率介于400~㊃09㊃长江大学学报(自然科学版)2023年3月1200m D ,泥质含量介于14%~20%,饱和度中值压力介于0.08~0.15M P a,孔隙结构以中小孔中细喉为主(见图9(c ))㊂Ⅳ类流动单元储层质量最差,研究区Ⅳ类流动单元占比13.6%㊂该类流动单元砂体厚度小于3m ,指标逼近程度大于0.7,孔隙度小于20%,渗透率小于400m D ,泥质含量大于20%,饱和度中值压力大于0.15M P a ,孔隙结构以小孔细微喉为主(见图9(d))㊂图7 单井流动单元识别划分F i g .7 F l o wu n i t i d e n t i f i c a t i o na n d c l a s s i f i c a t i o no f s i n gl ew e l ls 图8 取心井不同流动单元统计特征F i g.8 S t a t i s t i c a l c h a r a c t e r i s t i c s o f d i f f e r e n t f l o wu n i t s i n c o r ew e l l s 流动单元的平面分布和沉积微相具有良好的匹配关系㊂研究区为浅水三角洲前缘沉积,以分流砂坝和分流河道为主要沉积骨架,其中分流砂坝为优势沉积单元㊂Ⅰ㊁Ⅱ类流动单元只分布在储层质量较好㊃19㊃第20卷第2期吴穹螈等:浅水三角洲储层流动单元划分的分流砂坝以及分流河道的中心部位;Ⅲ类流动单元在分流河道和溢岸砂中普遍发育,Ⅳ类流动单元则主要分布在厚度薄且物性差的溢岸砂以及分流间湾中(见图10)㊂图9 不同流动单元微观孔隙结构特征F i g .9 M i c r o s c o pi c p o r e s t r u c t u r e c h a r a c t e r i s t i c s o f d i f f e r e n t f l o wu n i ts 图10 流动单元及沉积微相平面分布F i g .10 P l a n e d i s t r i b u t i o no f f l o wu n i t s a n d s e d i m e n t a r y mi c r o f a c i e s 4.2 连通性分析由于沉积环境以及宏观物性㊁微观孔隙结构的差异,同一泛连通体内部不同流动单元之间往往存在渗流差异界面,进而对泛连通体内部侧向连通性产生重要影响[22-25]㊂研究区具有丰富的生产动态资料,产液剖面和吸水剖面资料表明,流动单元组合样式对侧向连通性具有明显的控制作用㊂分流砂坝和分流㊃29㊃长江大学学报(自然科学版)2023年3月河道内部发育多种流动单元,不同流动单元之间呈切叠式连续分布,这也反映出复合砂体内部垂向上是由多期单一成因单元砂体叠置而成(见图11)㊂图11 流动单元剖面分布F i g.11 P r o f i l e d i s t r i b u t i o no f f l o wu n i t s Ⅰ类流动单元之间连通性极好,生产动态表现为注采迅速受效,注入水短期内迅速突进,采油井无水采油期多小于3个月,单井产油量大于100m 3/d ;Ⅰ㊁Ⅱ类之间或者Ⅱ类之间接触,连通性较好,注采受效,注入水缓慢推进,含水率上升较为缓慢,采油井无水采油期多介于3~8个月,单井产油量多介于60~100m 3/d ;Ⅰ㊁Ⅱ类与Ⅲ类之间接触,连通性中等,注采受效缓慢,采油井无水采油期多介于8~16个月,单井产油量多介于30~60m 3/d ;Ⅰ㊁Ⅱ㊁Ⅲ类与Ⅳ类之间接触,连通性差,注采受效困难,采油井无水采油期多大于16个月,单井产油量多小于30m 3/d㊂4.3 挖潜实例笔者根据流动单元划分结果及分布特征,有针对性地提出了挖潜措施,即对不同流动单元接触区部署调整井来挖掘内部剩余油㊂2021年底完钻了采油井X 11井,该井位于注水井X 6和注水井X 7之间,流动单元呈现出 Ⅱ类-Ⅲ类-Ⅱ类 接触样式,分析认为X 11井区域为Ⅲ类流动单元,注采受效缓慢,剩余油较为富集㊂实钻结果表明,X 11井在目的层钻遇大套油层,水淹比例较低(见图12),初期日产油40m 3,无水采油期为9个月,目前含水率较低,仅为7.8%㊂由此可见,通过熵权T o p s i s 法对储层进行流动单元划分,能够较好地反映储层内部非均质性情况,为挖潜高含水期储层内部剩余油提供了科学的指导㊂5 结论1)与传统流动带指数法㊁聚类分析法㊁综合参数法㊁孔隙结构参数法等方法相比,熵权T o ps i s 法能够定量客观赋予评判参数权重,受主观因素影响较小,且计算简便,具有较好的实用性和可操作性,是划分储层流动单元的一种有效方法㊂2)孔隙度㊁渗透率㊁泥质含量和饱和度中值压力4种参数能够较好地反映沉积环境以及储层宏观物性和微观孔隙结构,可用于研究区储层流动单元划分,研究区发育4类流动单元,以Ⅱ㊁Ⅲ类为主㊂3)流动单元的平面分布受到沉积相的控制,流动单元的组合样式控制了连通性,不同流动单元间形成渗流差异界面,使得复合砂体内部弱连通或不连通,导致剩余油富集㊂㊃39㊃第20卷第2期吴穹螈等:浅水三角洲储层流动单元划分㊃49㊃长江大学学报(自然科学版)2023年3月图12X11井实钻情况F i g.12A c t u a l d r i l l i n g o fw e l l X11参考文献:[1]H E A R NCL,E B A N K S WJ J,T Y ERS.G e o l o g i c a l f a c t o r s i n f l u e n c i n g r e s e r v o i r p e r f o r m a n c eo f t h eH a r t z o g D r a wf i e l d,W y o m i n g[J]. P e t r o lT e c h,1984,36:1335-1344.[2]康立明,任战利.多参数定量研究流动单元的方法 以鄂尔多斯盆地W93井区为例[J].吉林大学学报(地球科学版),2008, 38(5):749-756.K A N GL M,R E NZL.S t u d y o n M u l t i-P a r a m e t e r sD i s c r i m i n a t i o n M e t h o d f o rF l o w U n i t s:T a k i n g W93W e l l b l o c k i nO r d o sB a s i na sA nE x a m p l e[J].J o u r n a l o f J i l i nU n i v e r s i t y(E a r t hS c i e n c eE d i t i o n),2007,28(4):99-103.[3]袁丙龙,张辉,叶青,等.基于三角洲复合砂体构型的流动单元划分及剩余油分布模式[J].沉积学报,2021,39(5):1253-1263.Y U A NBL,Z H A N G H,Y E Q,e ta l.F l o w-u n i t c l a s s i f i c a t i o nb a s e do nc o m p o u n ds a n d-b o d y a r c h i t e c t u r eo fd e l t aa n dd i s t r i b u t i o n p a t t e r no f r e m a i n i n g o i l[J].A c t aS e d i m e n t o l o g i c aS i n i c a,2021,39(5):1253-1263.[4]王石,万琼华,陈玉琨,等.基于辫状河储层构型的流动单元划分及其分布规律[J].油气地质与采收率,2015,22(5):47-51. WA N GS,WA N Q H,C H E N Y K,e t a l.F l o wu n i t s d i v i s i o na n d t h e i r d i s t r i b u t i o n l a wb a s e do nb r a i d e d r i v e r r e s e r v o i r a r c h i t e c t u r[J]. P e t r o l e u m G e o l o g y a n dR e c o v e r y E f f i c i e n c y,2015,22(5):47-51.[5]商建霞,张乔良,李华,等.基于储层构型的渗流单元刻画方法及应用 以珠江口盆地文昌C油田为例[J].沉积学报,2021, 39(4):1020-1030.S HA N GJX,Z HA N G QL,L IH,e t a l.C h a r a c t e r i z a t i o nm e t h o d a n d a p p l i c a t i o no f s e e p a g e u n i t s b a s e do n r e s e r v o i r c o n f i g u r a t i o n:ac a s e s t ud y o f t h eWe n C h a n g O i lf i e l d i n t h eP e a r lR i v e rM o u t hB a s i n[J].A c t aS e d i m e n t o l og i c aS i n i c a,,2021,39(4):1020-1030.[6]吕明针,林承焰,张宪国,等.储层流动单元划分方法评价及优选[J].岩性油气藏,2015,27(1):74-80,88.L Y U M Z,L I NCY,Z H A N GXG,e t a l.E v a l u a t i o n a n d o p t i m i z a t i o n o f f l o wu n i t d i v i s i o nm e t h o d s[J].L i t h o l o g i cR e s e r v o i r s,2015, 27(1):74-80,88.[7]贾胜彬.古城油田泌浅10区储层流动单元划分及储集质量评价[J].长江大学学报(自然科学版),2017,17(3):25-28,33. J I ASB.F l o wu n i t c l a s s i f i c a t i o na n d r e s e r v o i r q u a l i t y e v a l u a t i o no f f a u l t b l o c kB i q i a n10i nG u c h e n g O i l f i e l d[J].J o u r n a l o fY a n g t z e U n i v e r s i t y(N a t u r a l S c i e n c eE d i t i o n),2017,17(3):25-28,33.[8]高伟,王允诚,刘宏.流动单元划分方法的讨论与研究[J].西南石油大学学报(自然科学版),2009,31(6):37-40.G A O W,WA N G YC,L I U H.T h e d i s c u s s i o n a n d s t u d y o n t h e d i v i s i o no f f l o wu n i t[J].J o u r n a l o f S o u t h w e s t P e t r o l e u m U n i v e r s i t y (S c i e n c e a n dT e c h n o l o g y E d i t i o n),2009,31(6):37-40.[9]C A N A S JA,M A L I KZA,W uCH,C h a r a c t e r i z a t i o n o f f l o wu n i t s i n s a n d s t o n e r e s e r v o i r s:L a C i r a F i e l dC o l o m b i a S o u t hA m e r i c a[R].S P E 27732,1994.[10]T IG M,O G B E D O,MU N L Y W,e t a l.U s eo f f l o wu n i t s a sa t o o l f o r r e s e r v o i rd e s c r i p t i o n:Ac a s es t u d y[J].S P EF o r m a t i o nCopyright©博看网. All Rights Reserved.E v a l u a t i o n ,1995,10(2):122-128.[11]A H R I MA N K O S H M ,K A S I R IN ,MO U S A V IS M .I m p r o v e d p e r m e a b i l i t yp r e d i c t i o no f ah e t e r o g e n e o u sc a r b o n a t er e s e r v o i ru s i n ga r t i f i c i a l n e u r a l n e t w o r k sb a s e do nt h ef l o wz o n e i n d e xa p p r o ac h [J ].P e t r o l e u m S c i e n c ea nd Te c h n o l o g y,2011,29(23):2494-2506.[12]朱永进,尹太举,刘玲利.浅水型三角洲沉积研究进展及问题讨论[J ].石油天然气学报,2011,33(3):22-26.Z HU YJ ,Y I N TJ ,L I U L L .P r o g r e s sa n dd i s c u s s i o no ns h a l l o w -w a t e rd e l t as e d i m e n t [J ].J o u r n a lo fO i l a n dG a s t e c h n o l o g y,2011,33(3):22-26.[13]尹太举,张昌民,朱永进,等.叠覆式三角洲 一种特殊的浅水三角洲[J ].地质学报,2014,88(2):263-272.Y I N TJ ,Z H A N GC M ,Z HU YJ ,e ta l .O v e r l a p p i n g d e l t a :A n e ws p e c i a l t y p eo fd e l t af o r m e db y o v e r l a p pe dl o b e s [J ].A c t a G e o l o g i c a l S i n i c a ,2014,88(2):263-272.[14]段冬平,侯加根,刘钰铭,等.河控三角洲前缘沉积体系定量研究 以鄱阳湖三角洲为例[J ].沉积学报,2014,32(2):270-277.D U A N DP ,H O UJG ,L I U Y M ,e t a l .Q u a n t i t a t i v e r e s e a r c hof f l u v i a l -d o m i n a t e dd e l t a f r o n t s e d i m e n t a r y s y s t e m :Ac a s e s t u d y of P o y a ng L a k e d e l t a [J ].A c t aS e d m e n t o l o g i c aS i n i c a ,2014,32(2):270-277.[15]董凤娟,卢学飞,琚惠娇,等.基于熵权T o ps i s 法的低渗透砂岩流动单元划分[J ].地质科技情报,2012,31(6):124-128.D O N GFJ ,L U XF ,J U HJ ,e t a l .D i v i s i o no f f l o wu n i t so f t h e l o w p e r m e a b i l i t y s a n d s t o n e r e s e r v o i rb a s e do ne n t r o p h y -w e i g h t i n gT o p s i sm e t h o d [J ].G e o l o g i c a l S c i e n c e a n dT e c h n o l o g y I n f o r m a t i o n ,2012,31(6):124-128.[16]郑占,吴胜和,许长福,等.克拉玛依油田六区克下组冲积扇岩石相及储层质量差异[J ].石油与天然气地质,2010,31(4):463-471.Z H E N GZ ,WUSH ,X U CF ,e t a l .L i t h o f a c i e sa n dr e s e r v o i r so f a l l l u v i a l f a n i nt h eL o w e rK e r a m a y F o r m a t i o n i nt h eb l o c k -6o f K a r a m a y O i l f i e l d ,t h e J u n g g a rB a s i n [J ].O i l&G a sG e o l o g y,2010,31(4):463-471.[17]秦国省,胡文瑞,邹存友,等.基准面旋回对冲积扇沉积演化及储层质量的控制作用 以百口泉油田百21井区百口泉组为例[J ].油气地质与采收率,2018,25(3):1-10.Q I N GS ,HU W R ,Z O U CY ,e t a l .E f f e c t s o f b a s e -l e v e l c y c l e so na l l u v i a l f a ne v o l u t i o na n d r e s e r v o i r q u a l i t y :ac a s e s t u d y of t h e B a i k o u q u a nF o r m a t i o n i nB a i 21B l o c k ,B a i k o u q u a nO i l f i e l d [J ].P e t r o l e u m G e o l og y a n dR e c o v e r y E f f i c i e n c y ,2018,25(3):1-10.[18]吴胜和,王仲林.陆相储层流动单元研究的新思路[J ].沉积学报,1999,17(2):252-257.WUSH ,WA N GZL .An e w m e t h o do fn o n -m a r i n er e s e r v o i r f l o wu n i t s t u d y [J ].A c t aS e d i m e n t o l o gi c aS i n i c a ,1999,17(2):252-257.[19]施玉娇,高振东,王起琮,等.碎屑岩储层流动单元划分及特征 以陕北富昌地区延长组长2段储层为例[J ].岩性油气藏,2009,21(4):99-104.S H IYJ ,G A OZD ,WA N G QZ ,e t a l .T h e d i v i s i o n a n d c h a r a c t e r o f c l a s t i c r e s e r v o i r f l o w i n g u n i t :A n e x a m p l e f r o mC h a n g 2r e s e r v o i r o f Y a n c h a n g F o r m a t i o n i nF u c h a n g a r e a ,N o r t h e r nS h a n x i [J ].L i t h o l o g i cR e s e r v o i r s ,2009,21(4):99-104.[20]万琼华,罗伟,梁杰,等.基于储层构型的流动单元渗流屏障级次研究[J ].西南石油大学学报(自然科学版),2019,41(1):77-84.WA N Q H ,L U O W ,L I A N GJ ,e ta l .R e s e r v o i ra r c h i t e c t u r e -b a s e dc l a s s i f i c a t i o no f s e e p a geb a r r i e r so f f l o w u n i t [J ].J o u r n a lo f S o u t h w e s t P e t r o l e u m U n i v e r s i t y (S c i e n c e a n dT e c h n o l o g y E d i t i o n ),2019,41(1):77-84.[21]刘联群,李勇,张建国,等.多参数流动单元在储层评价中的应用[J ].西北大学学报(自然科学版),2009,39(1):114-120.L I U LQ ,L IY ,Z H A N GJG ,e t a l .A p p l i c a t i o n o fm u l t i -pa r a m e t e r s f l o wu n i t s i n e v a l u a t i o n o f r e s e r v o i r [J ].J o u r n a l o fN o r t h w e s t U n i v e r s i t y (N a t u r a l S c i e n c eE d i t i o n ),2009,39(1):114-120.[22]王志章,何刚.储层流动单元划分方法与应用[J ].天然气地球科学,2010,21(3):362-366.WA N GZZ ,H EG .D i v i s i o no f r e s e r v o i r f l o wu n i t a n d i t s a p pl i c a t i o n [J ].N a t u r a lG a sG e o s c i e n c e ,2010,21(3):362-366.[23]周路路,刘俊刚,舒广强,等.松辽盆地北部三肇凹陷升554断块特低渗油气储层流动单元研究[J ].中国地质,2021,48(3):807-819.Z H O U LL ,L I UJG ,S HU G Q ,e t a l .R e s e a r c ho n f l o wu n i t s o f u l t r a -l o w p e r m e a b i l i t y o i l -g a s r e s e r v o i r s i nS h e n g 554f a u l t b l o c k o f S a n z h a od e p r e s s i o n ,n o r t h e r nS o n g l i a oB a s i n [J ].G e o l o g y i nC h i n a ,2021,48(3):807-819.[24]李照永.三肇凹陷低渗透扶杨油层流动单元划分[J ].断块油气田,2011,18(1):66-69.L I ZY .F l o wu n i t d i v i s i o n o f F u y a n g r e s e r v o i rw i t h l o w p e r m e a b i l i t y i n S a n z h a o d e pr e s s i o n [J ].F a u l t -b l o c kO i l a n dG a s F i e l d ,2011,18(1):66-69.[25]陈雷,胡望水,李韵龙.流动单元研究及其在油气开发中的应用[J ].重庆科技学院学报(自然科学版),2013,15(6):53-58.C H E NL ,HU W S ,L IYL .T h e s t u d y o f f l o wu n i t a n d i t s a p p l i c a t i o n i no i l a n d g a s d e v e l o p m e n t [J ].J o u r n a l o f C h o n g q i n g U n i v e r s i t y of S c i e n c e a n dT e c h n o l og y (N a t u r a l S c i e n c e sE d i t i o n ),2013,15(6):53-58.[编辑] 帅群㊃59㊃第20卷第2期吴穹螈等:浅水三角洲储层流动单元划分 Copyright ©博看网. All Rights Reserved.。
流动单元定量论文流动单元研究论文
流动单元定量论文流动单元研究论文储层流动单元定量研究方法与注意问题摘要：自C．L．Hearn等人于1984年提出流动单元（Flow Unit）的概念[1]，目前人们一致性的认识是：“横向和垂向上连续的，具有相似渗透率、孔隙度及流动特征的储集带”。
流动单元作为一种油气田地质研究理论和表征方法，已经得到快速推广与发展。
流动单元的研究应以石油地质学、沉积岩石学、层序地层学、油层物理以及渗流力学等学科的理论为基础，综合应用油田动、静态生产资料，以储集体的非均质性特征和渗流特征为主线，对地下油气藏进行分单元的综合研究和评价。
关键词：流动单元；沉积相带；高分辨率层序1流动单元研究方法与关键技术1.1流动单元定量化研究进展流动单元研究实质上是对油气藏进行深入细致的综合研究和表征。
特别是对于复杂和长期注水开发的油气藏，要正确揭示地下油气藏内流动单元的分布规律，必须采用多学科的理论为指导，应用多种手段，多种信息，才能做好流动单元的识别、划分和评价。
本文通过文献调研，把国内外储层流动单元研究的基本思路和步骤大体分为两种类型：1.1.1 以数学手段为主的储层参数分析法广泛应用储层中的各种地质参数，通过单井中密集取样的聚类分析寻找划分流动单元的有效参数和定量界限，然后直接在整套储层中定量划分流动单元。
在这类方法中，仅仅需要做少量的地层对比和沉积学研究，隔、夹层的分布可作为流动单元划分的分界依据，最终建立以流动单元为基础的三维定量地质模型。
上述方法存在以下不足：强调了成因单元（或沉积相带）内影响流体渗流地质参数的差异，虽应用多种参数进行流动单元划分，但对成因单元本身的分布、单元间渗流屏障（沉积屏障、成岩胶结带和断层遮挡）及各种地质界面的研究很难用统一的定量标准进行刻画。
而陆相储层砂体时、空分布的复杂性、渗流屏障及微地质界面的分布状况，在某种程度上对地下流体运动的影响更为重要。
此外，某些研究者过分强调流动单元在垂向上的分层性，甚至将流动单元看作是更细微的“地层单元”，忽视了平面上渗流差异的研究，而平面上严重的非均质性正是陆相储层的又一特色。
一种流动单元自动分层新方法及其应用
模糊数学方法在划分流动单元中的应用
１ 模糊 综 合 评 价 的原 理
设 为评 价 因素集 ， 素就是 研 究 对象 的各种 因 属性 或数 量 。设 为评 语 集 ， 表示 研 究对 象 的评 价 结果 。当考 虑 的因 素很 多 、 因素之 间又 有层 次 之 各
级的综 合模糊 评价 ¨ 。
为 了使各项 因素 都参 与评 价 ， 须建 立从 到 必 （ 的模 糊 映射 ， ） 即诱 导 出模糊 变换矩 阵
２ 流动 单 元 的 模 糊综 合 评 价
２１ 评 价方 法 ．
ｒ … ｒ
＝
ｒ ］
（） ＝ ； ｒ ｆ
Ｌ … ｒ
Ｉ
ｒｎ ｍｊ
（） １
流动 元的 价参 多， 括有 厚度、 单 评 数较 包 效 孔隙
度、 渗透 率 和 泥 质 含 量 等 ３，－ ￣ 如 何 综 合 利 用 １ ，１ １ 。 －７ １２
收稿 日期 ：０ ０ ０ — ０； ２ １－ ６ ３ 改 Ｉ Ｈ期 ：０Ｏ ０ — ３ ２１－８ １。 作 者 简 介 ： 长春 ， ，ｌ ９ ， 郭 男 ： 从事 油 …开 发 地 质 产 及科 研 工作 。 联 系 电 活 ：０ ４ ） ７５ ２ ， －／ｌ 程 ｒ ｌ ｉ ｉ 一 （５ ６ ８ １４ ７ Ｅ Ｉ ｉ ｌ ； 皋 金 项 日 ： 因石 油 化 Ｔ集 公 司科 技 攻 关项 目“ 油 层 挖 潜 配 套技 术 ” Ｐ 】３ ） 巾 ： 厚 （ ｏ ０ ５
式 中 ： 为各 因素在研 究 评价 中 的权重 向量 ； Ａ ａ 为第 ｋ个 因素 的权 重 系数 。 由权 重 向量 和模 糊变换 矩 阵经过模 糊运 算得到 综合评 价 指数 Ｂ ＝［ ｌ６ ， ， ｂ ，２ … ｂ ］ 式 中： 为综 合评 价指数 ； 第 ｋ个评 价对 象 ｂ为
利用渗流系数划分特低渗储层流动单元
（ ） 沉 积特 征 ， 且 具有 相 同 （ ） 流体流 动 （ 似 的 而 似 的 渗
收稿 日期 ：Ｏ Ｏ ３ ５ ２ １ 一０ —２
的储 层非 均 质单 元 ， 对储 层 结 构 模型 （ 积 模型 ） 是 沉 的进 一 步 细 划 和 定 量 表 征 ［。 ６ ］
作者简介： 贺顺义（９４ ）男， １７－ ， 甘肃静宁人， 在读博士， 石油地质专业， 现就读于中国地质大学（ 北京）Ｃ ｉ ｎ ｅ ｉ （ｈ ａ ｉ ｒ． ｎ Ｕ ｖ ｓ
位 的 基 础 上 ， 行 了Ｃ４ 进 ＋ ’ 层 “ 性 ” 系研 究 ， 心 分 析 结 果 标 明 孔 隙 度 一 渗 透 率 相 关 性 好 。在 细 分 油 四 关 岩
单 砂 层 和 沉 积 微 相 的 基 础 上 ， 用 单 砂 体 厚 度 加 权 平 均 方 法 ， 用 声 波 时 差 曲 线 求 取 孔 隙 度 ， 而 求 取 采 利 进
及 剩 余 油 分 布 研 究 等提 供 了地 质 依 据 。 关 键 词 ： 动 单 元 ； 流 系数 ； 低 渗 储 层 ； 质 统 计 ｝ 尔 多 斯 盆 地 流 渗 特 地 鄂
中图分 类号 ： ＴＥ８ １
文 献标 识码 ： Ａ
文 章 编 号 ：０ ６ ７ ８ （ Ｏ Ｏ Ｏ 一 Ｏ １ 一 Ｏ １０－ ９ １２１ ） ９ Ｏ ４ ３ 流 ） 征 ， 样 就 可 为 油 藏 数 值 模 拟 时 合 理 的层 系 划 特 这
１ ４
内 蒙古 石 油 化 工
２１ 年第 ９ ００ 期
利 用 渗 流 系数 划 分特 低 渗储 层 流 动单 元
贺 顺 义 师永 民 谢 楠。方 国庆 ， ， ，
（． 国地 质 大 学 ， 京 １中 北 １０ ８ ；２ 京 大 学 石 油 天 然 气研 究 中 心 ， 京 ０ ０ ３ ．北 北 １０ ７ ０８１
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
性 、渗 流等属 性参 数 的判别 分 析 ，建立 了研 究 区流
动单元 划 分模 型 。
１ 定 量 化参 数 的确 定
流 动单 元 划分 参 数 的选取 要 体 现 宏 观 与微 观 、 构造 与沉 积 、岩 石 骨架 与流 体性 质等 各个 方 面 的特 征 。本 次研 究划 分 流动单 元的 参数 主要包 括沉 积微 相 、孔 隙度 、渗 透 率 、油 藏 品质 指 数 ＲＱＪ及 流 动 带 指标 Ｆ ＺＩ等 。 （ ）反 映沉 积环 境 的特 征 参 数— — 沉 积 微 相 。 １ 储层 物性 的好 坏 和流体 运 动特 征在 很 大程度 上受 控 于沉 积特 征 ，不 同沉积 微 相控 制各 种成 因砂 体 的外 部形 态 和 内部 物性 的分 布特 征 ，在应 用 沉积 微相 作 为划 分流 动单 元 的参数 时 ，必须 将各 种微 相类 型转
４ ２ 流 动 带 指 数 分 析 ．
综 合 以上分析 ，Ｅ类 流动 单元 的平 均渗 透率 为
１ ５ ７ ×１ 一 肚 ，平均 孔 隙度 为 ２ ． ６ ，平 均 ２．４ ０ ｍ ４ １ Ｒ １为 ｏ ６ ，平均 Ｆ Ｉ为 ２ １ ，分布 最广 ，属 于 Ｑ ．９ Ｚ ．７ 高孔 高渗型储 层 ，主要 为河道 主体 带与河 口坝 砂体
２ 聚 类 分 析 法
聚类 分 析是 最常 用 的数理 统计 方法 之一 ，它是
按照 客体 在性 质 上或 成 因上 的亲疏 关 系 ，对 客体 进 行定 量分 类 的多元 统计 分 析方 法 。这种 分类 方法 不 仅综合 考 虑 全 部 因素 ，且 不 受 已有 分 类 结 构 的影 响 ，只是 以某 种分 类统 计量 为依 据 ，对 客体 进行分
关键 词 ：流 动单元 ；聚类分 析 ；剩余 油 分布 ；参 数
引 言
流动单 元 （Ｆｏ Ｕｎｔ） 又 称 水 力 单 元 （ ｌｗ ｉ ｓ Ｈｙ
—
（ ）反 映 岩 石 孔 隙 结 构 的 特 征 参 数 —— Ｆ Ｉ ４ Ｚ
（ 动带 指标 ） 流 。
ＦＺＩ一 一 （） ２
Ｄ ａ ｌ ｉ ） ｒｕｉ Ｕｎｔ ，是 指 影 响 流 体 流 动 的 岩 相 、岩 ｃ ｓ
编 号 为 ２ 河 道 砂 （ 括 分 流 河 道 、水 下 分 流 河 ； 包 道） 对应编 号为 ３ ， ；河 口 坝 ，对 应 编 号 为 ４ 。在 此
元 样品数为 ９ ，说 明 Ｐ ｘ小层 的流 动性 能普遍 较 个 ｘ 好 。各个 流动单元 的储层参数 统计结果见表 １ 。
表 １ Ｐ ｘ小 层 聚 类 分析 结 果 参数 统 计 ｘ
藏 品质 指 数 ＲＱ Ｉ及 流 动 带 指 标 Ｆ ＺＩ等 ５种 参 数 ，
首先通 过灰色 理论 确定各 个参数 的权 值 ，然后 以各
参数 自身评分 与各 参数 的权系数 相乘 得到 各参数 单
项 得分 ，最终 把各参 数单 项得分 累积后 按百 分制计 算 ，得 到流动 单元 的最终 综合得 分 ，依 此反 映地 下 地 质的 真实情 况 。
１中 国石油 勘探 开发 科学 研究 院 ２Ｓ ｈｕ ｅｇ ｒ ｈｎ ． ． 京 ３中国地 质大 学能 源学 院 ｃ ｌｍｂ ｒｅ ｉａＳ Ａ 北 Ｃ
摘 要 ：大 庆油 田某 井 区储 集 层流 动单 元采 用聚 类 分析 法可 划 分为 Ｅ、Ｇ 和 Ｐ共 ３种 类型 。Ｅ
类 主要 对应 河道 主体 带 与河 口坝砂 体 ，属 高孔 高渗 型储 层 ，综合 评 分 ６ ５左 右 ；Ｇ 类主 要 对 应 决 口扇 与分流 间席 状砂 体 ，属 高孔 中渗 型储 层 ，综合 评 分 ３ ８左 右 ；Ｐ类 绝 大部 分 对应 分 流 间席状 砂体 ，属 中孔 低 渗型储 层 ，综合评 分 １ ０左 右 。３类 流动 单 元 的 挖 潜 方 向 ：Ｅ类 箱 型 或 钟 型砂 体 的 顶部 易形成 剩 余 油富 集 区；Ｇ 类和 Ｐ类作 为二 类储 层 ，是 下一 步挖 潜 的有利 区域 。
仅 以 Ｐ ｘ小层 为例进 行 阐述 。 ｘ
４ １ 储 层 参 数 分 析 ．
通过 聚类分 析结果 ，对 不 同流 动单元 进行 孑 隙 Ｌ
度 与渗 透率 、Ｆ Ｉ与 Ｒ Ｚ ＱＪ散 点 图分 析 表 明 ，Ｐ ｘ ｘ 小层 的 Ｅ类 流 动 单 元 孔 、渗 性 能 最 好 ，Ｇ 类 流 动
直线。
本 次 研 究 同 时 对 各 个 样 品 加 入 了 Ｆ 属 性 参 ＺＩ
通 过 聚类 分析 结果 ，对不 同流 动单元 进行平 面 投点作 图 ，得到不 同单元 在平 面上 的分 布 ；再 同该 小层 的沉积微 相 图进行叠 合 ，就得 出流动 单元 与沉 积微相 的关 系分布 图 ，见图 １ 、图 ２ 。
由式 （ ） １ 、式 （ ） 可 知 ，Ｒ ２ ＱＩ一 ： Ｚ 或 Ｆ
ｌｇ Ｒ ） ｏ （￣ ＋ ｌｇ Ｆ Ｉ ｏ （ ＱＩ ＝ｌｇ ｑ ） ｏ （ Ｚ ），表 明 在 Ｒ 和 ＱＪ
流 特征 ［ ］ ４ 。层 内流动 单元 的划 分 则 是在 确 定地 层
的垂 向层序 、隔 和夹 层分 布等 地质 特 征之后 ，依 据 各 小层 层 内的岩 石物 理特 征 ，尤其 是孔 隙度 和 渗透 率 的分 布特征 ，应 用 多参 数量 化分 析 的方法 进行 平 面分 区研 究 来 进 行 。通 过 对 储 集 层 内 的沉 积 、物
４ ３ 流 动 单 元 聚 类 分 析 ．
３ 多参 数 的流 动单 元 评 分
本 次研究 选取 沉积 微相 、孔 隙度 、渗 透率 、油
通过 以上 分 析 ，可 以确 定 将 Ｐ ｘ层 大 致 划 分 ｘ 出 ３类流动 单元 。因此 ，在 对各个 样 品点聚类 分析 时 ，指定分 类个 数为 ３大类 。在对 多种 聚类方 法分 类结 果 比较 后 ，选 用 衡 量样 品亲疏 程 度 算 法 （ ｕ Ｅ— ｅｄ ｌ ）和样 品 聚类 算法 （ ｅａ ｅ ｉ Ａｖ ｒｇ ）进 行分类 。 Ｐ ｘ 层流 动单元 聚类 分析 的有效样 品点 总数 ｘ 小 为 ２ ２个 ，聚类结 果表 明 ，Ｅ类 流 动单 元样 品 数为 ８ ２ ０个 ，Ｇ类流动 单元 样 品数 为 ５ ２ ３个 ，Ｐ类 流 动单
４ 研 究 实 例
研究 工 区 目的层 段为 白垩系 姚家组 一段 葡萄 花 油层 组 ，岩性 为灰 白色 、灰 绿 色砂 岩 夹 有灰 绿 色 、 棕 红色泥 岩 ， 目的层段 厚 度 ０ ０ ｍ。本 区在 葡 １ ～６ 段发 育三 角洲平 原与三 角洲前 缘沉 积 ，北 部发 育三 角洲 平原 ，中、南部发 育三 角洲前 缘沉 积 ；主要微 相类 型 包 括 泛 滥 平 原 、分 流河 道 、决 口扇 、分 流 间 、水下分 流河 道 、河 口坝 、席状 砂等 ；分流 河道 北西 一南东走 向 ；河 口坝分 ３个带 发育 。砂岩单 层 厚度 ４ ０ｍ，层 内和平面非 均质性 较 强 。将砂 体 ～１ 类型综 合划 分为 ４大类 ：河 间砂 （ 括席 状砂 、泛 包 滥平 原 、分 流 问 ） ，对 应 编 号 为 １ ；决 口扇 ，对 应
数 分 析 ，样 品颜 色 值 越 红 ，Ｆ Ｉ值 越 高 。从 实 际 Ｚ 数据 投点 可 以发现 ，Ｐ ｘ小层 Ｆ Ｉ高值 样 品点 多 ， ｘ Ｚ 而低 值 较 少 ，说 明 其 总 体 孔 、渗 性 能 和流 动 性 能 好 ，这 与该 层 分流河 道 、水 下分 流河道 及河 口坝砂 体 发育 密切有 关 。
层 流 动 单 元 储 层 参 数 分 析 可 见 ，Ｐ ｘ小 层 在 ＲＱＩ ｘ
性差 异较大 ，而 Ｇ 类 流 动单 元 与 Ｅ类 流 动单 元 物
性差异 小 。 ４ ４ 流 动 单 元 平 面 分 布 规 律 ．
和 的 双 对 数 坐 标 上 ，有 ３组 斜 率 近 似 平 行 的
亲密 的 聚集合 并 为一 类 ，把另 一些 彼此 之 间亲 近的
换为 代码来 参 与运算 。 （ ）反 映储层物性特征 的参数—— Ｐ Ｒ 与 Ｋ 。 ２ Ｏ
（ ）反 映油 藏质 量 特征 的参 数 —— Ｒ ３ ＱＩ（ 油藏 品质指 数 ） 。
ｒ
客体 集合 为另 一类 。在客 体 聚合为 类 的基础 上 ，根
类 。根 据 分类对 象 不 同 ，聚类分 析分 为样 品 聚类分 析 和变 量 聚类分 析 ２种类 型 。
流 动单 元 的分类 是根 据 被观察 对 象 的变量 值来 进行 ，即采 用样 品 聚类分 析 法 。其 原 理是在 开 始时 每个 客 体 自成 一 类 ，然后 以某 种表 示 客体 亲疏关 系 的分 类 统计量 为 分类 依据 ，把 一些 彼此 之 间关 系最
一
的双对 数 坐 标 上 ，具 有 相 近 Ｆ Ｉ值 的样 品将 落 在 Ｚ 条斜 率 为 ｌ的直 线 上 ，具 有不 同 Ｆ ＺＩ值 的样 品 将 落在 一组 平行 线 直线 上 。而 同一直 线上 的样 品具 有 相似 的孔 喉 （ Ｌ 孑 隙结构 ） 特征 ，从 而构 成一 个水 力 流动 单元 。
Ｒ Ｉ ｏｏ１√ Ｑ — ．３４
一 芒
（ １ ）
据类 之 间 的亲疏 程度继 续 合并 ，直 到全 部客 体 聚为
一
类 为止 。
基金 论 文 ：高 等 学 校 博 士 学 科 点 专 项 科 研 基 金 （０ ５ ４ １ ０ ） ２ ０ ０ ９ ０ １
３ ６
《 气田地 面工程》第 ３ 卷第 １ ２１．１ 试验研究） 油 ０ 期（０ １０ ）＜
《 气田地 面工程》第 ３ 卷第 １ ２１．１ 试验研究＞ 油 ｏ 期（０ １０ ）（
ｄ ｉ １ ． ９ ９ ｊ ｉｓ ． ０ ６６ ９ ． ０ １ １ ０ ４ ｏ ： ０ ３ ６ ／ ． ｓｎ １ ０ － ８ ６ ２ １ ． ． １