CN110617062B - 一种储层剩余气分布特征的分析方法 - Google Patents
一种储层剩余气分布特征的分析方法 Download PDFInfo
- Publication number
- CN110617062B CN110617062B CN201910972008.7A CN201910972008A CN110617062B CN 110617062 B CN110617062 B CN 110617062B CN 201910972008 A CN201910972008 A CN 201910972008A CN 110617062 B CN110617062 B CN 110617062B
- Authority
- CN
- China
- Prior art keywords
- water
- low
- degree
- namely
- gas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 33
- 238000004458 analytical method Methods 0.000 title abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 196
- 239000007789 gas Substances 0.000 claims abstract description 143
- 238000004519 manufacturing process Methods 0.000 claims abstract description 106
- 238000000034 method Methods 0.000 claims abstract description 85
- 230000035699 permeability Effects 0.000 claims abstract description 58
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 16
- 238000011160 research Methods 0.000 claims abstract description 13
- 238000005065 mining Methods 0.000 claims abstract description 10
- 230000003068 static effect Effects 0.000 claims abstract description 9
- 230000007246 mechanism Effects 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000010586 diagram Methods 0.000 claims abstract description 3
- 230000009545 invasion Effects 0.000 claims description 117
- 238000000605 extraction Methods 0.000 claims description 97
- 239000010410 layer Substances 0.000 claims description 32
- 239000011435 rock Substances 0.000 claims description 17
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 12
- 230000035945 sensitivity Effects 0.000 claims description 12
- 239000011229 interlayer Substances 0.000 claims description 11
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 230000003247 decreasing effect Effects 0.000 claims description 8
- 230000000704 physical effect Effects 0.000 claims description 7
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 6
- 239000003345 natural gas Substances 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 5
- 239000008398 formation water Substances 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 230000033764 rhythmic process Effects 0.000 claims description 3
- 239000002253 acid Substances 0.000 claims description 2
- 239000003513 alkali Substances 0.000 claims description 2
- 239000004568 cement Substances 0.000 claims description 2
- 239000002131 composite material Substances 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 239000012466 permeate Substances 0.000 claims description 2
- 239000011148 porous material Substances 0.000 claims description 2
- 150000003839 salts Chemical class 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 239000000758 substrate Substances 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
Abstract
本发明涉及一种储层剩余气分布特征的分析方法,其特征在于,包括以下步骤:步骤1,收集整理研究区块的静态资料和动态资料;步骤2,计算出单井的储量、水体能量、采气速度、渗透率变异系数;步骤3,作出步骤2计算结果的等值图;步骤4,在步骤3中的等值图上圈出单一储层符合某一标准的区域,并进行面积叠加,将研究区块划分为不同的小区域,分类研究不同剩余气分布特征及其形成机理。本发明优点在于,能够结合现场资料,将大的研究区块划分为剩余气分布特征不同的小区域,为制定不同剩余气分布特征的开采对策提供基础,实用性强。
Description
技术领域
本发明涉及一种储层剩余气分布特征的分析方法,属于天然气开发领域。
背景技术
当气田开采一定阶段之后,研究剩余气分布特征成为了一个十分重要的课题。
关于剩余气分布特征的分析方法,目前大多数学者仅对剩余油分布特征进行了相关的研究。并且现有的方法大多没有结合生产资料,仅从油藏岩石、流体固有的物理性质对剩余油进行分析,准确度不够。申请号为201010234800.1的专利申请《一种分析缝洞型油藏剩余油分布的方法》通过数值模拟技术找准缝洞型油藏剩余油的分布位置。申请号为201920143726.9的专利申请《一种特高含水期剩余油饱和度分布电法测定实验装置》、申请号为201310151735.X的专利申请《一种三维模拟油藏剩余油饱和度测量方法及装置》采用实验的方法确定不同位置的剩余油饱和度含量。因此迫切需要一个***的、全面的分析方法来分类研究不同剩余气分布特征及其形成机理,为气田综合调整,制定不同剩余气分布特征的开采对策,提高采收率提供科学依据和理论支持。
发明内容
本发明的目的是:分类研究不同剩余气分布特征及其形成机理,为制定不同剩余气分布特征的开采对策提供基础。
为达到上述目的,本发明采用以下技术方案:
一种储层剩余气分布特征的分析方法,所述方法包括以下步骤:
步骤1,收集整理研究区块的静态资料和动态资料;
步骤2,根据步骤1收集到的资料计算出单井的F储量、G水体能量、H单井年采气速度、I渗透率变异系数;
步骤3,结合Am含气边界、Al井位坐标和步骤2计算得到的单井各项数据,作出F储量、G水体能量、H单井年采气速度、I渗透率变异系数的等值图;
步骤4,在步骤3中的等值图上圈出单一储层符合某一标准的区域,并进行面积叠加,将研究区块划分为不同的小区域,分类研究不同剩余气分布特征及其形成机理。
进一步地,步骤1中静态资料包括:A基础数据、B储层特征、C岩石物性、D流体物性。动态资料包括:E生产参数。
进一步地,静态资料中的A基础数据包括:Aa孔隙度、Ab原始含气饱和度、Ac储层有效厚度、Ad含气面积、Ae渗透率、Af原始地层压力、Ag体积系数、Ah储层深度、Ai储层温度、Aj射孔层位、Ak射孔深度、Al井位坐标、Am含气外边界、An井距;B储层特征包括:Ba圈闭、Bb储层非均质性、Bc储层潜在伤害;C岩石物性包括:Ca孔隙度、Cb敏感性、Cc胶结类型、Cd渗透性、Ce饱和度、Cf相对渗透率;D流体物性包括:Da天然气、Db水。动态资料中的E生产参数包括:Ea生产日期、Eb压力、Ec产气量、Ed产水量、Ee气水比。
进一步地,Ak射孔深度可分为:Ak1射开小层上部、Ak2射开小层中部、Ak3射开小层下部;An井距可分为:An1大于经济极限井距、An2小于经济极限井距;Ba圈闭可分为:Ba1构造圈闭、Ba2地层圈闭、Ba3岩性圈闭、Ba4复合圈闭;Bb储层非均质性可分为:Bb1层间非均质性、Bb2平面非均质性、Bb3层内非均质性、Bb4微观非均质性、Bb5各向异性;Bc储层潜在伤害可分为:Bc1吼道堵塞、Bc2水封、Bc3应力敏感;Ca孔隙度可分为:Ca1致密(砂岩1~10%,碳酸盐岩1~5%)、Ca2中等(砂岩10~20%,碳酸盐岩5~10%)、Ca3好的(砂岩20~35%,碳酸盐岩10~20%);Cb敏感性可分为:Cb1水敏、Cb2盐敏、Cb3酸敏、Cb4碱敏、Cb5速敏;Cc胶结类型可分为:Cc1基底胶结、Cc2孔隙胶结、Cc3接触胶结;Cd渗透性可分为:Cd1渗透性极好(>1000mD)、Cd2渗透性好(100~1000mD)、Cd3渗透性中等(10~100mD)、Cd4渗透性微弱(1~10mD)、Cd5非渗透性的(<1mD);Ce饱和度可分为:Ce1含气饱和度、Ce2束缚水饱和度;Da天然气可分为:Da1干气、Da2湿气、Da3凝析气;Db水可分为:Db1地层水、Db2外来水。Eb压力可分为:Eb1油压、Eb2套压、Eb3油套压差;Ec产气量可分为:Ec1日产气量、Ec2年产气量、Ec3累计产气量;Ed产水量可分为:Ed1日产水量、Ed2年产水量、Ed3累计产水量;Ee气水比可分为:Ee1不变、Ee2增大、Ee3减小。
进一步地,Ba1构造圈闭可分为:Ba11背斜、Ba12断层、Ba13裂缝、Ba14岩性穿刺;Ba2地层圈闭可分为:Ba21地层不整合、Ba22地层超覆;Bb3层内非均质性可分为:Bb31正韵律、Bb32反韵律、Bb33复合韵律;Bc1吼道堵塞可分为:Bc11出砂、Bc12水锁;Cb1水敏可分为:Cb11膨胀性、Cb12崩解性、Cb13软化性;Db1地层水可分为:Db11边水、Db12底水、Db13层间水、Db14层内水、Db15束缚水。Eb3油套压差可分为:Eb31不变、Eb32增大、Eb33减小;Ec1日产气量可分为:Ec11不变、Ec12增大、Ec13减小;Ed1日产水量可分为:Ed11不变、Ed12增大、Ed13减小。
进一步地,步骤2中的F储量包括:Fa地质储量、Fb剩余地质储量丰度、Fc动态储量、Fd采出程度;G水体能量包括:Ga水侵常数、Gb存水体积系数、Gc水侵替换系数、Gd累计水侵量、Ge水驱指数、Gf水侵速度、Gg水体倍数、Gh水侵时间、Gi水侵方向;H单井年采气速度包括:Ha低(<0.1%)、Hb高(>0.1%);I渗透率变异系数包括:Ia层间渗透率变异系数、Ib平面渗透率变异系数、Ic层内渗透率变异系数。
进一步地,Fb剩余地质储量丰度可分为:Fb1低(<0.2×108m3/km2)、Fb2高(>0.2×108m3/km2);Fd采出程度可分为:Fd1低(<40%)、Fd2高(>40%);Ga水侵常数可分为:Ga1低(<6)、Ga2高(>6);Ia层间渗透率变异系数可分为:Ia1小(<0.3)、Ia2大(>0.3);Ib平面渗透率变异系数可分为:Ib1小(<0.3)、Ib2大(>0.3);Ic层内渗透率变异系数可分为:Ic1小(<0.3)、Ic2大(>0.3)。
进一步地,步骤2中的计算过程可分为以下7步:
(1)结合Ak井位坐标、Al含气外边界将小层含气面积劈分得到单井控制面积;
(2)计算单井控制的Fa地质储量、Fb剩余地质储量丰度;
(3)结合Eb1油压、Eb2套压、Ec1日产气量、Ed1日产水量、Ea生产日期计算得到Fc动态储量;
(4)计算单井的Fd采出程度、Ga水侵常数、Gb存水体积系数、Gc水侵替换系数、Gd累计水侵量、Ge水驱指数、Gf水侵速度、Gg水体倍数;
(5)作出Eb1油压、Eb2套压、Ec1日产气量、Ed1日产水量与Ea生产日期的关系曲线,分析曲线,得知油套压差增大,Ed1日产水量上升,Ec1日产气量下降的日期即为Gh水侵时间;
(6)计算H单井年采气速度;
(7)计算单井的Ia层间渗透率变异系数、Ib平面渗透率变异系数、Ic层内渗透率变异系数。
进一步地,步骤5中不同剩余气分布特征及其形成机理表现出的特征有以下40种:
(1)特征1为Fb1,即剩余地质储量丰度低,;
(2)特征2为Fb2+Fd2,即剩余地质储量丰度高,采出程度高;
(3)特征3为Fb2+Fd1+Ga1+Hb,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高;
(4)特征4为Fb2+Fd1+Ga1+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面非均质性严重;
(5)特征5为Fb2+Fd1+Ga1+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,垂向非均质性严重;
(6)特征6为Fb2+Fd1+Ga1+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,井距大;
(7)特征7为Fb2+Fd1+Ga1+Hb+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面非均质性严重;
(8)特征8为Fb2+Fd1+Ga1+Hb+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,垂向非均质性严重;
(9)特征9为Fb2+Fd1+Ga1+Hb+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,井距大;
(10)特征10为Fb2+Fd1+Ga1+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面与垂向非均质性严重;
(11)特征11为Fb2+Fd1+Ga1+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面非均质性严重,井距大;
(12)特征12为Fb2+Fd1+Ga1+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,垂向非均质性严重,井距大;
(13)特征13为Fb2+Fd1+Ga1+Hb+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面与垂向非均质性严重;
(14)特征14为Fb2+Fd1+Ga1+Hb+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面非均质性严重,井距大;
(15)特征15为Fb2+Fd1+Ga1+Hb+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,垂向非均质性严重,井距大;
(16)特征16为Fb2+Fd1+Ga1+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面与垂向非均质性严重,井距大;
(17)特征17为Fb2+Fd1+Ga1+Hb+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面与垂向非均质性严重,井距大;
(18)特征18为Fb2+Fd1+Ga1+Cb1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,储层岩石水敏;
(19)特征19为Fb2+Fd1+Ga2+Ha,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低;
(20)特征20为Fb2+Fd1+Ga2+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面非均质性严重;
(21)特征21为Fb2+Fd1+Ga2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,垂向非均质性严重;
(22)特征22为Fb2+Fd1+Ga2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,井距大;
(23)特征23为Fb2+Fd1+Ga2+Ha+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面非均质性严重;
(24)特征24为Fb2+Fd1+Ga2+Ha+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,垂向非均质性严重;
(25)特征25为Fb2+Fd1+Ga2+Ha+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,井距大;
(26)特征26为Fb2+Fd1+Ga2+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面与垂向非均质性严重;
(27)特征27为Fb2+Fd1+Ga2+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面非均质性严重,井距大;
(28)特征28为Fb2+Fd1+Ga2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,垂向非均质性严重,井距大;
(29)特征29为Fb2+Fd1+Ga2+Ha+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面与垂向非均质性严重;
(30)特征30为Fb2+Fd1+Ga2+Ha+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面非均质性严重,井距大;
(31)特征31为Fb2+Fd1+Ga2+Ha+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,垂向非均质性严重,井距大;
(32)特征32为Fb2+Fd1+Ga2+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面与垂向非均质性严重,井距大;
(33)特征33为Fb2+Fd1+Ga2+Ha+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面与垂向非均质性严重,井距大;
(34)特征34为Fb2+Fd1+Ga2+Bc3+Cc1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,储层岩石应力敏感强且胶结物含量多;
(35)特征35为Fb2+Fd1+Ga2+Cd4/Cd5,即剩余地质储量丰度高,采出程度低,水侵影响程度小,渗透率低;
(36)特征36为Fb2+Fd1+Ga2+Da3,即剩余地质储量丰度高,采出程度低,水侵影响程度小,且气藏为凝析气藏;
(37)特征37为Fb2+Fd1+Cb+Db2,即剩余地质储量丰度高,采出程度低,储层岩石敏感性强且有外来水渗入储层;
(38)特征38为Fb2+Fd1+Ak1+Bb31(/+Ak3+Bb32),即剩余地质储量丰度高,采出程度低,射孔时射开小层上部且为正韵律(或射孔时射开小层下部且为反韵律);
(39)特征39为Fb2+Fd1+Ba12,即剩余地质储量丰度高,采出程度低,气藏为断层气藏;
(40)特征40为Fb2+Fd1+Ba13,即剩余地质储量丰度高,采出程度低,气藏为裂缝气藏。
与现有技术相比,本发明具有以下有益效果:能结合现场实际,方便、迅速地将大的研究区块划分为剩余气分布特征不同的小区块,分类研究不同剩余气分布特征及其形成机理,为制定不同剩余气分布特征的开采对策提供基础,实用性强。
附图说明
图1剩余气分布特征的分析方法;
图2静态资料详细内容;
图3动态资料详细内容;
图4计算结果与分析详细内容;
图5本发明实施例剩余气分布特征划分结果。
1为剩余地质储量丰度低;2为剩余地质储量丰度高,采出程度高;6为剩余地质储量丰度高,采出程度低,水侵影响程度大,井距大;7为剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面非均质性严重;11为剩余地质储量丰度高,采出程度低,水侵影响程度大,平面非均质性严重,井距大;12为剩余地质储量丰度高,采出程度低,水侵影响程度大,垂向非均质性严重,井距大;16为剩余地质储量丰度高,采出程度低,水侵影响程度大,平面与垂向非均质性严重,井距大;17为剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面与垂向非均质性严重,井距大;25为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,井距大;27为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,井距大;30为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面非均质性严重,井距大;31为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,垂向非均质性严重,井距大;32为剩余地质储量丰度高,采出程度低,水侵影响程度小,平面与垂向非均质性严重,井距大;33为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面与垂向非均质性严重,井距大。
具体实施方式
结合附图对本发明作进一步说明。
实施例如下所示,附图1所示的剩余气分布特征的分析方法,针对国内某气田某小层,实施以下技术步骤:
第一步,收集整理研究区块的静态资料和动态资料。由测井资料得到各井在该研究小层的孔隙度、原始含气饱和度、储层有效厚度、渗透率、储层深度。根据现场实际得到各井井位坐标、含气外边界、井距。根据射孔层位统计表可以得到射孔层位、射孔深度。
根据生产日报可以得到每天的油压、套压、产气量、产水量。
第二步,计算相关结果。经过一系列的计算,得到该小层各井的相关结果如表1所示。
表1
第三步,作图并进行面积叠加。结合含气边界数据、井位数据和表1中的单井计算结果,作出相应的等值图,并分别在每张等值图上圈出该小层剩余地质储量丰度<0.2×108m3/km2、采出程度>40%、水侵常数>6、单井年采气速度>0.1%、井距<400m、平面渗透率变异系数>0.3、垂向渗透率变异系数>0.3的区域。
第四步,面积划分与剩余气分布特征分析。对圈出的区域进行面积叠加,将小层划分为不同的小区域,如附图5所示。该小层共划分了56个小区域,包含14种剩余气分布特征。每个区域中标注的数字即为对应的特征,如下:
1为剩余地质储量丰度低;
2为剩余地质储量丰度高,采出程度高;
6为剩余地质储量丰度高,采出程度低,水侵影响程度大,井距大;
7为剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面非均质性严重;
11为剩余地质储量丰度高,采出程度低,水侵影响程度大,平面非均质性严重,井距大;
12为剩余地质储量丰度高,采出程度低,水侵影响程度大,垂向非均质性严重,井距大;
16为剩余地质储量丰度高,采出程度低,水侵影响程度大,平面与垂向非均质性严重,井距大;
17为剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面与垂向非均质性严重,井距大;
25为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,井距大;
27为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,井距大;
30为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面非均质性严重,井距大;
31为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,垂向非均质性严重,井距大;
32为剩余地质储量丰度高,采出程度低,水侵影响程度小,平面与垂向非均质性严重,井距大;
33为剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面与垂向非均质性严重,井距大。
相应的,不同的剩余气分布特征应当制定不同的开采对策。如特征1区域,剩余地质储量丰度低,在现有技术条件下可以暂不进行开采。特征2区域,剩余地质储量丰度高,采出程度大,天然气有效动用,可以加强生产监控,合理论证井距。特征6、7、11、12、16、17区域水侵影响程度大,应当在东南部位,即s3-35、s3-13处布置强排井。在水侵严重时,采速过大会加剧边水的窜入以及底水的锥进,即特征7、17区域应当适当降低采速;而在水侵影响较小时,采速过低会导致经济效益低,如特征25、27、30、31、33区域应当提高采速。特征6、11、12、16、17、25、27、30、31、32、33区域,井距过大,远井区域采出程度低,可以合理加密井距,提高井控范围。特征7、11、16、17、30、32、33区域平面非均质性严重,应当调整工作制度,改善储层渗流。特征12、16、17、31、32、33区域垂向非均质性严重,应当分层开采,降低垂向非均质性的影响。
对于本领域的技术人员而言,基于本发明的思想可以做出各种添加、替换、组合,这些添加、替换、组合也均属于本发明的保护范围。
Claims (8)
1.一种储层剩余气分布特征的分析方法,其特征在于,所述方法包括以下步骤:
步骤1,收集整理研究区块的静态资料和动态资料;
步骤2,根据步骤1收集到的资料计算出单井的F储量、G水体能量、H单井年采气速度、I渗透率变异系数;
步骤2中的计算过程可分为以下7步:
(1)结合Ak井位坐标、A1含气外边界将小层含气面积劈分得到单井控制面积;
(2)计算单井控制的Fa地质储量、Fb剩余地质储量丰度;
(3)结合Ebl油压、Eb2套压、Ec1日产气量、Ed1日产水量、Ea 生产日期计算得到Fc动态储量;
(4)计算单井的Fd采出程度、Ga水侵常数、Gb存水体积系数、Gc水侵替换系数、Gd累计水侵量、Ge水驱指数、Gf水侵速度、Gg 水体倍数;
(5)作出 Eb1油压、Eb2套压、Ec1日产气量、Ed1日产水量与Ea 生产日期的关系曲线,分析曲线,得知油套压差增大,Ed1日产水量上升,Ec1日产气量下降的日期即为Gh水侵时间;
(6)计算H单井年采气速度;
(7)计算单井的Ia层间渗透率变异系数、Ib平面渗透率变异系数、Ic层内渗透率变异系数;
步骤3,结合Am含气边界、Al井位数据和步骤2计算得到的单井各项数据,作出F储量、G水体能量、H单井年采气速度、I渗透率变异系数的等值图;
步骤4,对步骤3中的等值图上圈出单一储层符合某一标准的区域,并进行面积叠加,将研究区块划分为不同的小区域,分类研究不同剩余气分布特征及其形成机理。
2.根据权利要求1所述的方法,其特征在于,步骤1中静态资料包括:A基础数据、B储层特征、C岩石物性、D流体物性;动态资料包括:E生产参数。
3.根据权利要求2所述的方法,其特征在于,静态资料中的A基础数据包括:Aa孔隙度、Ab原始含气饱和度、Ac储层有效厚度、Ad含气面积、Ae渗透率、Af原始地层压力、Ag体积系数、Ah储层深度、Ai储层温度、Aj射孔层位、Ak射孔深度、Al井位坐标、Am含气外边界、An井距;B储层特征包括:Ba圈闭、Bb储层非均质性、Bc储层潜在伤害;C岩石物性包括:Ca孔隙度、Cb敏感性、Cc胶结类型、Cd渗透性、Ce饱和度、Cf相对渗透率;D流体物性包括:Da天然气、Db水;动态资料中的E生产参数包括:Ea生产日期、Eb压力、Ec产气量、Ed产水量、Ee气水比。
4.根据权利要求3所述的方法,其特征在于,Ak射孔深度可分为:Ak1射开小层上部、Ak2射开小层中部、Ak3射开小层下部;An井距可分为:An1大于经济极限井距、An2小于经济极限井距;Ba圈闭可分为:Ba1构造圈闭、Ba2地层圈闭、Ba3岩性圈闭、Ba4复合圈闭;Bb储层非均质性可分为:Bb1层间非均质性、Bb2平面非均质性、Bb3层内非均质性、Bb4微观非均质性、Bb5各向异性;Bc储层潜在伤害可分为:Bc1吼道堵塞、Bc2水封、Bc3应力敏感;Ca孔隙度可分为:砂岩孔隙度范围为1~10%或碳酸盐岩孔隙度范围为1~5%即为Ca1致密、砂岩孔隙度范围为10~20%或碳酸盐岩孔隙度范围为5~10%即为Ca2中等、砂岩孔隙度范围为20~35%或碳酸盐岩孔隙度范围为10~20%即为Ca3好的;Cb敏感性可分为:Cb1水敏、Cb2盐敏、Cb3酸敏、Cb4碱敏、Cb5速敏;Cc胶结类型可分为:Cc1基底胶结、Cc2孔隙胶结、Cc3接触胶结;Cd渗透性可分为:渗透率>1000mD即为Cd1渗透性极好、渗透率范围为100~1000mD 即为Cd2渗透性好、渗透率范围为10~100mD 即为Cd3渗透性中等、渗透率范围为1~10mD 即为Cd4渗透性微弱、渗透率<1mD即为Cd5非渗透性的;Ce饱和度可分为:Ce1含气饱和度、Ce2束缚水饱和度;Da天然气可分为:Da1干气、Da2湿气、Da3凝析气;Db水可分为:Db1地层水、Db2外来水;Eb压力可分为:Eb1油压、Eb2套压、Eb3油套压差;Ec产气量可分为:Ec1日产气量、Ec2年产气量、Ec3累计产气量;Ed产水量可分为:Ed1日产水量、Ed2年产水量、Ed3累计产水量;Ee气水比可分为:Ee1不变、Ee2增大、Ee3减小。
5.根据权利要求4所述的方法,其特征在于,Ba1构造圈闭可分为:Ba11背斜、Ba12断层、Ba13裂缝、Ba14岩性穿刺;Ba2地层圈闭可分为:Ba21地层不整合、Ba22地层超覆;Bb3层内非均质性可分为:Bb31正韵律、Bb32反韵律、Bb33复合韵律;Bc1吼道堵塞可分为:Bc11出砂、Bc12水锁;Cb1水敏可分为:Cb11膨胀性、Cb12崩解性、Cb13软化性;Db1地层水可分为:Db11边水、Db12底水、Db13层间水、Db14层内水、Db15束缚水;Eb3油套压差可分为:Eb31不变、Eb32增大、Eb33减小;Ec1日产气量可分为:Ec11不变、Ec12增大、Ec13减小;Ed1日产水量可分为:Ed11不变、Ed12增大、Ed13减小。
6.根据权利要求1所述的方法,其特征在于,步骤2中的F储量包括:Fa地质储量、Fb剩余地质储量丰度、Fc动态储量、Fd采出程度;G水体能量包括:Ga水侵常数、Gb存水体积系数、Gc水侵替换系数、Gd累计水侵量、Ge水驱指数、Gf水侵速度、Gg水体倍数、Gh水侵时间、Gi水侵方向;H单井年采气速度包括:采气速度<0.1%即为Ha低、采气速度>0.1%即为Hb高;I渗透率变异系数包括:Ia层间渗透率变异系数、Ib平面渗透率变异系数、Ic层内渗透率变异系数。
7.根据权利要求6所述的方法,其特征在于,Fb剩余地质储量丰度可分为:剩余地质储量丰度<0.2×108m3/km2即为Fb1低、剩余地质储量丰度>0.2×108m3/km2即为Fb2高;Fd采出程度可分为:采出程度<40%即为Fd1低、采出程度>40%即为Fd2高;Ga水侵常数可分为:水侵常数<6即为Ga1低、水侵常数>6即为Ga2高;Ia层间渗透率变异系数可分为:层间渗透率变异系数<0.3即为Ia1小、层间渗透率变异系数>0.3即为Ia2大;Ib平面渗透率变异系数可分为:平面渗透率变异系数<0.3即为Ib1小、平面渗透率变异系数>0.3即为Ib2大;Ic层内渗透率变异系数可分为:层内渗透率变异系数<0.3即为Ic1小、层内渗透率变异系数>0.3即为Ic2大。
8.根据权利要求1所述的方法,其特征在于,步骤4中不同剩余气分布特征及其形成机理表现出的特征有以下40种:
(1)特征1为Fb1,即剩余地质储量丰度低;
(2)特征2为Fb2+Fd2,即剩余地质储量丰度高,采出程度高;
(3)特征3为Fb2+Fd1+Ga1+Hb,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高;
(4)特征4为Fb2+Fd1+Ga1+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面非均质性严重;
(5)特征5为Fb2+Fd1+Ga1+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,垂向非均质性严重;
(6)特征6为Fb2+Fd1+Ga1+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,井距大;
(7)特征7为Fb2+Fd1+Ga1+Hb+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面非均质性严重;
(8)特征8为Fb2+Fd1+Ga1+Hb+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,垂向非均质性严重;
(9)特征9为Fb2+Fd1+Ga1+Hb+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,井距大;
(10)特征10为Fb2+Fd1+Ga1+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面与垂向非均质性严重;
(11)特征11为Fb2+Fd1+Ga1+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面非均质性严重,井距大;
(12)特征12为Fb2+Fd1+Ga1+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,垂向非均质性严重,井距大;
(13)特征13为Fb2+Fd1+Ga1+Hb+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面与垂向非均质性严重;
(14)特征14为Fb2+Fd1+Ga1+Hb+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面非均质性严重,井距大;
(15)特征15为Fb2+Fd1+Ga1+Hb+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,垂向非均质性严重,井距大;
(16)特征16为Fb2+Fd1+Ga1+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,平面与垂向非均质性严重,井距大;
(17)特征17为Fb2+Fd1+Ga1+Hb+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,采气速度高,平面与垂向非均质性严重,井距大;
(18)特征18为Fb2+Fd1+Ga1+Cb1,即剩余地质储量丰度高,采出程度低,水侵影响程度大,储层岩石水敏;
(19)特征19为Fb2+Fd1+Ga2+Ha,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低;
(20)特征20为Fb2+Fd1+Ga2+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面非均质性严重;
(21)特征21为Fb2+Fd1+Ga2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,垂向非均质性严重;
(22)特征22为Fb2+Fd1+Ga2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,井距大;
(23)特征23为Fb2+Fd1+Ga2+Ha+Ib2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面非均质性严重;
(24)特征24为Fb2+Fd1+Ga2+Ha+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,垂向非均质性严重;
(25)特征25为Fb2+Fd1+Ga2+Ha+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,井距大;
(26)特征26为Fb2+Fd1+Ga2+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面与垂向非均质性严重;
(27)特征27为Fb2+Fd1+Ga2+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面非均质性严重,井距大;
(28)特征28为Fb2+Fd1+Ga2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,垂向非均质性严重,井距大;
(29)特征29为Fb2+Fd1+Ga2+Ha+Ib2+Ia2/Ic2,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面与垂向非均质性严重;
(30)特征30为Fb2+Fd1+Ga2+Ha+Ib2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面非均质性严重,井距大;
(31)特征31为Fb2+Fd1+Ga2+Ha+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,垂向非均质性严重,井距大;
(32)特征32为Fb2+Fd1+Ga2+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,平面与垂向非均质性严重,井距大;
(33)特征33为Fb2+Fd1+Ga2+Ha+Ib2+Ia2/Ic2+An1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,采气速度低,平面与垂向非均质性严重,井距大;
(34)特征34为Fb2+Fd1+Ga2+Bc3+Cc1,即剩余地质储量丰度高,采出程度低,水侵影响程度小,储层岩石应力敏感强且胶结物含量多;
(35)特征35为Fb2+Fd1+Ga2+Cd4/Cd5,即剩余地质储量丰度高,采出程度低,水侵影响程度小,渗透率低;
(36)特征36为Fb2+Fd1+Ga2+Da3,即剩余地质储量丰度高,采出程度低,水侵影响程度小,且气藏为凝析气藏;
(37)特征37为Fb2+Fd1+Cb+Db2,即剩余地质储量丰度高,采出程度低,储层岩石敏感性强且有外来水渗入储层;
(38)特征38为Fb2+Fd1+Ak1+Bb31,即剩余地质储量丰度高,采出程度低,射孔时射开小层上部且为正韵律;
(39)特征39为Fb2+Fd1+Ba12,即剩余地质储量丰度高,采出程度低,气藏为断层气藏;
(40)特征40为Fb2+Fd1+Ba13,即剩余地质储量丰度高,采出程度低,气藏为裂缝气藏。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910972008.7A CN110617062B (zh) | 2019-10-14 | 2019-10-14 | 一种储层剩余气分布特征的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910972008.7A CN110617062B (zh) | 2019-10-14 | 2019-10-14 | 一种储层剩余气分布特征的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110617062A CN110617062A (zh) | 2019-12-27 |
CN110617062B true CN110617062B (zh) | 2021-02-23 |
Family
ID=68925377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910972008.7A Active CN110617062B (zh) | 2019-10-14 | 2019-10-14 | 一种储层剩余气分布特征的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110617062B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111472722B (zh) * | 2020-03-12 | 2022-02-01 | 中国石油天然气股份有限公司 | 煤层气合采井分层产气能力预测方法以及装置 |
CN111222261B (zh) * | 2020-04-22 | 2020-07-14 | 西南石油大学 | 动静参数结合的产量劈分方法 |
CN111734357B (zh) * | 2020-08-10 | 2020-11-13 | 西南石油大学 | 一种基于人工智能的气藏单井控制面积劈分方法 |
CN114427432B (zh) * | 2020-09-11 | 2024-04-26 | 中国石油化工股份有限公司 | 一种气藏剩余气开发潜力确定方法 |
CN112595653B (zh) * | 2021-03-05 | 2021-05-07 | 西南石油大学 | 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法 |
CN115788385B (zh) * | 2022-12-29 | 2024-05-24 | 西南石油大学 | 一种高温凝析气藏凝析水气比的确定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104453834A (zh) * | 2014-10-31 | 2015-03-25 | 中国石油化工股份有限公司 | 一种井组注采关系优化调整方法 |
CN105822269A (zh) * | 2015-04-09 | 2016-08-03 | 中国石油化工股份有限公司 | 多层油藏二元驱后剩余油分布模式及开发调整方法 |
JP2016223064A (ja) * | 2015-05-27 | 2016-12-28 | 国立大学法人東北大学 | 二酸化炭素低排出発電方法および二酸化炭素低排出発電システム |
CN108104806A (zh) * | 2017-12-14 | 2018-06-01 | 中国石油化工股份有限公司 | 剩余油分布规律定量分析方法 |
CN109886559A (zh) * | 2019-01-25 | 2019-06-14 | 中国石油天然气集团有限公司 | 一种油田水淹层精细解释及剩余油综合评估方法 |
-
2019
- 2019-10-14 CN CN201910972008.7A patent/CN110617062B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104453834A (zh) * | 2014-10-31 | 2015-03-25 | 中国石油化工股份有限公司 | 一种井组注采关系优化调整方法 |
CN105822269A (zh) * | 2015-04-09 | 2016-08-03 | 中国石油化工股份有限公司 | 多层油藏二元驱后剩余油分布模式及开发调整方法 |
JP2016223064A (ja) * | 2015-05-27 | 2016-12-28 | 国立大学法人東北大学 | 二酸化炭素低排出発電方法および二酸化炭素低排出発電システム |
CN108104806A (zh) * | 2017-12-14 | 2018-06-01 | 中国石油化工股份有限公司 | 剩余油分布规律定量分析方法 |
CN109886559A (zh) * | 2019-01-25 | 2019-06-14 | 中国石油天然气集团有限公司 | 一种油田水淹层精细解释及剩余油综合评估方法 |
Non-Patent Citations (2)
Title |
---|
LW3-1气田总体开发方案-气藏工程研究;魏芳;《中国硕士学位论文全文数据库 工程科技Ⅰ辑》;20180516(第06期);B19-18页 * |
松南盆地孤家子气田气藏精细描述及剩余气分布规律研究;郭奇;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20171016(第11期);B19-2页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110617062A (zh) | 2019-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110617062B (zh) | 一种储层剩余气分布特征的分析方法 | |
Qin | Research progress of symbiotic accumulation of coal measure gas in China | |
CN103993862B (zh) | 复杂断块稀油油藏分层开发方法 | |
WO2019091048A1 (zh) | "五图-三带-两分区"保水采煤方法 | |
CN107664029B (zh) | 砂岩热储地热资源可循环利用最优井网布局方法 | |
CN111191849B (zh) | 一种西部矿区深埋工作面涌水量预测方法 | |
CN102747991A (zh) | 一种确定合采井单层产量的方法 | |
Taotao et al. | Characteristics of shale pores, fractal dimension and their controlling factors determined by nitrogen adsorption and mercury injection methods | |
CN105735960A (zh) | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 | |
CN110359895B (zh) | 一种非均质巨厚砂岩水平井分段压裂的探放水方法 | |
CN103352693A (zh) | 一种基于含油气含量的产能预测***及其方法 | |
CN108894777A (zh) | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 | |
CN104295292A (zh) | 多层叠置煤层气***开采井设计方法 | |
CN107237627A (zh) | 一种岩溶型碳酸盐岩储层布井方法 | |
CN109798149A (zh) | 一种坚硬顶板厚煤层采场来压强度分级预测方法 | |
Bi et al. | Geological characteristics and co-exploration and co-production methods of Upper Permian Longtan coal measure gas in Yangmeishu Syncline, Western Guizhou Province, China | |
CN105678082B (zh) | 一种识别油气井酸压沟通储层类型的双压降法 | |
CN108197421B (zh) | 一种致密气和煤层气联合开发有利区量化评价方法 | |
Fu et al. | Gas accumulation conditions and key technologies for exploration & development of Sulige gasfield | |
Chen et al. | In situ coal permeability and favorable development methods for coalbed methane (CBM) extraction in China: From real data | |
Yang et al. | Reservoir accumulation conditions and key exploration & development technologies for Keshen gas field in Tarim Basin | |
Dou et al. | The coalbed methane production potential method for optimization of wells location selection | |
Liu et al. | Effect of sandstone and mudstone thickness on artificial fracturing for hydrocarbon extraction from low-permeability reservoirs | |
Tian-xuan et al. | Optimization of goaf gas drainage parameters based on numerical simulation studying fracture in overlying strata | |
Zhang et al. | Mining thickness determination of upper and lower protective coal seam in pressure relief mining |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |