一种碳酸盐气藏储量动用状况的评价方法
【专利摘要】本发明涉及一种碳酸盐气藏储量动用状况的评价方法。本发明充分应用测井资料多元组合得到一种定量评价气藏储量动用状况的渗流特征参数，该参数充分考虑了储层的岩性、泥质含量、含气性、渗透性、孔隙结构类型及孔隙结构系数等，可以较为准确地反映碳酸盐气藏小层层间、层内渗流能力与储量动用程度的差异性。有效利用随机建模手段建立三维渗流特征参数连续地质模型，开展渗流特征参数井间预测，克服了传统插值方法（如克里金方法）对参数的平滑效应，真实地反映了碳酸盐气藏储层的非均质性，为气藏储量的有效动用提供技术支持，夯实了气藏持续稳产基础，达到提高气藏最终采收率的目的。
【专利说明】一种碳酸盐气藏储量动用状况的评价方法
【技术领域】
[0001]本发明涉及一种碳酸盐气藏储量动用状况的评价方法，属于碳酸盐气藏预测与评价【技术领域】。
【背景技术】
[0002]碳酸盐气藏高含硫、气藏埋藏深、储层厚度大、沉积和成岩相变快、储层纵横向上非均质性严重，气藏储量空间上动用不均匀，开发难度大；另外，长井段、非均质气井多层合采，且井筒积液现象普遍，液面以下还有气层生产，产剖测试精度受到影响，产剖资料显示，测试产气厚度小于实际产出厚度，液面以下气层产气量、出力层位置不准，纵向上小层储量的采出状况不清，给气藏动态分析、开发调整方案的制定带来难度，严重到影响气藏持续稳产。
[0003]常规评价碳酸盐气藏储量动用状况的主要方法:(I)采用井网控制程度、采收率、动态储量、采出程度等参数进行计算、评价，虽然这些参数能够反映气藏宏观的储量动用差异，但随着气藏开发的深入，对碳酸盐气藏储量动用评价精度要求越来越高，需要对气层层间、层内进行精细评价。由于气藏层间和层内的储量动用的影响因素较多，除了受储层物性的影响外，还受储集的有效性、微观孔隙结构的差异性、小层的开发条件等因素的影响，因此，采用常规方法显然不能满足目前碳酸盐气藏精细开发的要求；(2)碳酸盐气藏采用丛式井网，开发井距大，在100m以上，由于井间储层非均质强，采用常规方法评价气藏井间的储量动用状况准确率低；另外，采用常规方法只能通过二维图件来评价气藏储量动用状况，存在一定局限性，在一定程度上掩盖了储层纵横向上的非均质性，不能进行井间气层储量的三维定量化及可视化预测；(3)应用产剖测试技术能够直接定量评价气层纵向上层间的采出状况、储量动用程度，但受井筒积液、井筒测试条件不具备等因素的影响，该方法评价不够准确，产剖测试资料仅能作为气藏储量动用状况评价的参考。
【发明内容】
[0004]本发明的目的是针对现有方法的精度低、不能提供三维可视化直观图件、难以满足碳酸盐气藏精细开发的需要等缺陷，提供了一种利用综合性测井评价参数、三维可视化直观图件定量评价碳酸盐气藏的储量动用状况的方法。
[0005]本发明包括以下步骤:
[0006]1、确定产出层与非产层产剖特征
[0007]利用气田产剖测试资料，结合单井地质特性，研究分析测试井的产出层与非产层产剖特征，为气藏储量动用状况评价参数优选提供依据。
[0008]2、储层渗流特征参数的建立
[0009]2.1建立测试井的测井评价参数与产层、非产层特征的关系图版，确定出区分产层与非产层的评价单参数，包括补偿中子、深侧向电阻率参差比、含气饱和度、自然伽玛；
[0010]2.2利用优选出的单参数进行多元组合，确定能定量区分产层与非产层的渗流特征参数。
[0011]3、利用储层的渗流特征参数，确定产层与非产层
[0012]3.1建立渗流参数与产层、非产层特征的关系图版，确定出可动用层下限标准，区分出产层与非产层；
[0013]3.2、建立渗流特征参数与米采气指数的关系图版，确定出产层优劣的分类标准。
[0014]4、气藏储量动用状况评价
[0015]4.1根据渗流特征参数，确定各井点的渗流特征参数值，在孔隙度、渗透率等相关属性模型的三维约束下，采用地质建模的手段，建立储层的三维渗流特征参数连续地质模型；
[0016]4.2通过对储层的三维渗流特征参数连续地质模型的约束控制、雕刻，得到气藏储量动用状况评价剖面图以及不同层段的气藏储量动用状况评价平面图，实现对碳酸盐气藏的储量动用状况三维、可视化评价。
[0017]本发明充分应用测井资料多元组合得到一种定量评价气藏储量动用状况的渗流特征参数，该参数充分考虑了储层的岩性、泥质含量、含气性、渗透性、孔隙结构类型及孔隙结构系数等，可以准确、定量反映碳酸盐气藏小层层间、层内渗流能力与储量动用程度的差异性。本发明有效利用随机建模手段建立三维渗流特征参数连续地质模型，开展渗流特征参数及气藏储量动用程度的井间预测，真实地反映了碳酸盐气藏储层的非均质性，由此生成的三维、可视化储量动用状况评价直观图件，能够准确地评价碳酸盐气藏空间储量动用状况，为气藏稳产开发阶段动态分析和开发调整提供依据，为气藏储量的有效动用提供技术支持，夯实了气藏持续稳产基础，达到提高气藏最终采收率的目的。
【专利附图】
【附图说明】
[0018]图1是本发明技术方案流程框图；
[0019]图2是利用现有技术编制的两口中国川东北取芯井普光104-1、302_1井不同孔隙结构类型的产出状况图；
[0020]图3是利用现有技术编制的中国川东北普光104-1井产层电测成果综合柱状图；
[0021]图4是利用现有技术编制的13 口产剖井测井孔隙度、渗透率与产层、非产层特征关系图；
[0022]图5是利用现有技术编制的中国川东北普光气田13 口产剖井测井孔隙度、饱和度与产层、非产层特征关系图；
[0023]图6是利用现有技术编制的中国川东北普光气田13 口产剖井深侧向电阻率参差比与产层、非产层特征关系图；
[0024]图7是利用现有技术编制的中国川东北普光气田2 口全直径取芯井普光302-1、104-1井深电阻率与微观孔隙结构系数、孔隙度关系图；
[0025]图8是利用现有技术编制的中国川东北普光气田6 口取心井不同孔隙结构的深侧向电阻率与孔隙结构微观系数关系与孔渗关系对比图；
[0026]图9是利用本发明编制的中国川东北普光气田13 口产剖井渗流特征参数与产层、非产层特征关系图版；
[0027]图10是利用本发明编制的中国川东北普光气田13 口产剖井渗流特征参数与米采气指数交会图版；
[0028]图11是利用本发明建立中国川东北普光气田三维渗流特征参数连续地质模型；
[0029]图12是利用本发明生成的中国川东北普光气田普光107-1H?305-2井和普光3011-5?103-4井两个方向的储量动用状况评价剖面图；
[0030]图13是利用本发明编制的中国川东北普光气田飞一二段上、中、下三个层段的储量动用状况评价平面图。
【具体实施方式】
[0031]下面结合中国川东北普光气田储量动用状况的评价实例和附图，对本发明实施方式做进一步详细说明，由图1可知，本发明具体步骤如下:
[0032]1、确定产层与非产层产剖特征
[0033]利用普光气田产剖测试资料，结合单井地质特性，研究分析测试井的产层与非产层产剖特征，为气藏储量动用状况评价参数优选提供依据。
[0034]根据普光气田2011年以来的13 口产剖测试井资料评价及地质特性分析，确定测试井产剖特征如表I所示，包括:
[0035]表I普光气田2011-2013年13 口产剖测试井解释成果表 ^ZTZ ^^Wm^占射孔厚度比例^
(m)(1mVd) ;zr 亏__；zr1iL__Cm)__(ra)__^_V>_
30.58-26 TIfS-TIf 1-2 4823.0-4933.5266.7
50.56-26 TIfS-TIf 1-2 4786.0-4933.529.57.6
388.6 ------
68.24-26 TlfS-TIf 1-2 4775.0-4933.5 33.58,6
86.9 4-26 TIf3-TIfl-2 4775.0-4933.5 36.59.4
3017-31TIf 1-25634.5—5688.0 25.911.3
4516-31TIf 1-25629.5-5688.0 34.415
104-1 229.2 ------
6016-32TIf1-25634.5-5694.04218.3
___72__16-31__TIfl-2__5634.5-5694.044.5__19.4_
30__69-85 Tlfl-2~P2ch 5258.0-5385.739.8__21_
302-1189.56069—85 Tlfl-2_P3ch 5258.0-5385.739.821
9069-85 Tlfl-2-P4ch 5258.0-5385.739.S21
25.863-79TIf 1-25667.7-5768.934.418.3
3563-79TIf 1-25667； 7-5768.9 34.418.3
I 187 7 ------
45.5 63-79TIf 1—25667； 7-5768.9 34.418.3
59.263-79TIf 1-25667.7-5768.9 34.418.3
29.2 47-64__TIfl~2__5127.6-5260.5 92.82__63.8_
59.8 47-64TIf 1—25127.6-5260.5 92.8263.8
202-1 145 5 ______
'91.6 47-64TIf 1-25127.6-5260.5 92.8263.S
114.447-64Tlfl-25127.6-5260.5 V>2.8263.8
38.6 15-33__Tlfl-2__5337.5-5501.1 35.7__11.5_
54.8 15-33 TIf 1-2 5337.5-5501.1 35.7 11.5
ΙΟ」-，10------
76.411-33TI￡l-25316.3-5501.1 37.312.0
10011-33TIf 1-25316.3-5501.14113.2
45.216-30TIf 1-24^16.5-4961.<S 29.TS, 9
70.716-30TIf 1-24916.5-4961.8 29.78.9
__OQO-OOOO Q______
L0036」一95.7 16-30__Tlfl-2__4Q16.P.-4961.(S 29.T__S1B_
___104.9 16-30__Tlfl-2__4916.5-4961.8 29.7___
27.5 62-73TIf 1—250<S2.S-5133.8.1Π.98.5
39.962-73TIf 1—2 5082.8-5133.837.08.5
398 ------
50.862-73TIf 1-2β~η1?Λ.8WH8.5
68.3 62-73TIf 1-25082.8-5133.8 36.68.5
19.5 66-67TIf1-26054-605952.9
29.1 66-67 Tlfl-2 6054-6059 5 2.9
105-2173.2 ------
39.466-67Tlfl-26054-605952.9
49.566-67Tlfl-26054-605952.9
40.3 7-8 ΤΙΠ-2 5130-5141 7 1.4
H>2 J 519.459.27-8TIf 1—25130—5141T1,4
79.57-8TIf 1-25130-514171.4
93.67-8TIf1-25130-5141 71.4
3375-95TIf 1-25409.5-548127.511.8
201-2I48.575-79TIf 1-25409.5-548127.511.8
58.275-63TIf 1-25409.5-548127.511.8
82.375-47TIf1-25409.5-548127.511.8
307-8TI￡l-25607-569734.517.7
304-1 195.138.47-8TIf 1—25607—569734.517.7
46.47-8Tlfl-25607-569734.517.7
58.37-8TIf 1-25607-569734.517.7
26.339-40Tlfl-26073.5-609016.56.19
10^-433.5 39-40TIf 1—26073.5-609016.56.19
45.7 39-40TIf 1-26073,5-609016, 56.19
[0037]1.1测试井产出通道较薄，测试产气厚度小于实际产出厚度；
[0038]13 口测试井产出厚度介于5?92.8m,占测试井射孔厚度的1.4-63.8%。其中，普光202-1井产气厚度占射孔厚度的比例最大，达到63.8% ;普光302-2井产气厚度占射孔厚度最小，仅为1.4%。分析认为储层非均质性强，存在平面渗流及纵向窜流现象，另外，由于井筒积液造成，液面以下出气层气量、位置测量不准，造成测试产气厚度小于实际产出厚度。
[0039]1.2随着产气量增大，生产压差增加，产出通道有所增加；
[0040]如普光201-1井产气量由30.5增加到86.9X 104m3/d时，产气厚度由26.0m增加到36.5m，产气比例由6.7%上升至Ij 9.4% ;普光104-1井产气量由30.0增加到72.0X 104m3/d时，产气厚度由25.9m增加到44.5m，产气比例由11.3%上升到19.4%。说明随着气井生产压差的增大，部分相对较差的气层开始动用。
[0041]1.3产出层一般为物性和孔隙结构相对较好的层段；
[0042]如普光104-1、302_1两口产剖测试井的粒内孔隙产出比例较小，仅占粒内孔隙样本点的37.8% ;粒间孔隙、混合孔隙产出比例较大，产出层比例达到62.9-69.8%，表明孔隙结构相对较好的粒间、混合孔隙产气能力相对较强，见图2。
[0043]如普光104-1产剖测试井产出层主要为一、二类储层，即孔隙度>5%的储层。该井产气剖面测试解释产气层为飞一二段40.6m，排除积液面以下的射孔层，产气厚度占液面以上射孔层段厚度的29.2%。其中一、二类层产气厚度32.2m,占总产气层厚度的79.3% ;三类层8.4m,占总产气层厚度的20.7%，见图3。
[0044]2、产层与非产层渗流特征参数的建立
[0045]2.1建立测试井的测井评价参数与产层、非产层特征的关系图版，确定出区分产层与非产层的评价单参数；
[0046]常规的测井评价参数种类多，包括孔隙度、渗透率、含气饱和度、深、中、浅电阻率、深、中、浅侧向电阻率参差比、密度、声波、补偿中子、自然伽玛、自然电位、微梯度、感应电导率、微观孔隙结构系数等，通过建立测试井的各类测井参数与产层、非产层特征关系图，见图4、图5、图6，从中优选出反映孔隙度的补偿中子、反映渗透性和微观孔隙结构差异的深侧向电阻率参差比、反映含气性的含气饱和度、反映泥质含量和沉积环境的自然伽玛4个单参数。
[0047]其中深侧向电阻率参差比为深侧向电阻率/深侧向电阻率参差系数，参差系数定义为各样品测量值累积百分数偏离“完全均质线”的相对值的平均值，也就是储层电阻率厚度加权平均值。
[0048]结合普光302-1、104_1井全直径样品分析资料，两口井的孔隙结构系数、深侧向电阻率相关系数达到0.88，说明可以通过深侧向电阻率率研究储层孔隙结构特征，如图7所示；不同孔隙结构储层孔隙度-深侧向电阻率关系，与孔-渗关系对比具有较大相似性，说明深侧向电阻率参差比可以反映储层的渗透性，如图8所示。
[0049]微观孔隙结构系数是用来表征真实岩石孔隙特征与长度相同的平行柱状毛细管束模型之间的差别，它的数值是影响这种差别的各种综合因素的度量。孔隙结构参数的测量采用毛管压力分析一压汞法。其计算公式:
Φ*Γ2
[0050]φρ =-
8K
[0051]式中:r是平均孔喉半径(相对应于平均毛管压力)，单位是孔隙度，单位％;K是渗透率，单位um2 ;当Φρ=1时，即为毛管束模型。
[0052]2.2利用优选出的单参数进行多元组合，确定能定量区分产层与非产层的渗流特征参数；
[0053]对优选出来的补偿中子、深侧向电阻率参差比、含气饱和度、自然伽玛4个单参数进行多元组合，定量评价储量动用状况的渗流特征参数Fs,其方法如下:
P ( Rlld/Clld) *CNL*Sg
[0054]Fs =---
GR
[0055]式中:Rlld/Clld—深恻向电阻率参差比，无单位；
[0056]Rlld—深侧向电阻率，Ω.m ;
[0057]Clld—深侧向电阻率参差系数，Ω.m ;
[0058]CNL—补偿中子，％ ；
[0059]GR—自然伽玛，API;
[0060]Sg—含气饱和度，% ；
[0061]3、利用储层的渗流特征参数，确定产层与非产层
[0062]3.1建立渗流参数Fs与产层、非产层特征关系图版，确定出可动用层下限标准，区分出产层与非产层；
[0063]通过建立的渗流特征参数Fs与产层、非产层特征关系图版见图9，确定可动用层下限标准，即渗流特征参数值Fs=0.1。由图可知，渗流特征参数值Fs在0.1以上的储层为非产层、差层；渗流特征参数值Fs在0.1以下的储层为产层。
[0064]3.2、建立渗流特征参数Fs与米采气指数的关系图版，确定出产层优劣的分类标准；
[0065]根据规范，气藏产能划分标准如表2所示
[0066]表2气藏产能划分标准
储层分类I类(好)II类(中) mL'
[0067]^^
lotf/d_____5 50__u b
^o-9 αι ?。.9 。?。厂
[0068]通过绘制产剖测试井小层渗流特征参数Fs与米采气指数关系图看出，如图10所示，可以利用渗流特征参数Fs进一步比较小层的产气能力大小，渗流特征参数匕<0.01，米采气指数> 1.0，为产量高的层，小层储量动用相对要好；渗流特征参数Fs介于0.01?0.1之间，米采气指数介于0.1?1.0之间，为产量中等的层，小层储量动用程度中等；渗流特征参数Fs ^ 0.1，米采气指数介于O?0.1的为差层或非产层，小层储量动用差或未动用，见表3。
[0069]表3普光碳酸盐气田产层与非产层划分标准
~~S?差或不产渗征^0.010.01-0.1--θ.1
参数____
米采气指彳ηΓΜ I ηπ π ,
I,U0,1?1.0O?O,I
[0070]数____
绝对无阻流量>505?卬υ
1iUi/d
储量动用紅
好中差未动用
[0071]4、气藏储量动用状况评价
[0072]4.1根据渗流特征参数Fs，确定各井点的渗流特征参数值Fs，在现有的气藏孔隙度、渗透率等相关属性模型的三维约束下，采用随机地质建模的手段，建立储层的三维渗流特征参数连续地质模型，如图11所示；
[0073]4.2通过对储层的三维渗流特征参数连续地质模型进行约束控制、雕刻，得到任意层位或方向上的气藏储量动用状况评价剖面图以及不同层段的气藏储量动用状况评价平面图，实现对碳酸盐气藏的储量动用状况三维、可视化评价；
[0074]依据步骤3.1可动用层下限标准、步骤3.2中的产层划分标准，对步骤4.1中的三维渗流特征参数连续地质模型进行控制约束、雕刻，生成气藏储量动用状况评价剖面图、平面图，如图12、图13-1、13-2、13-3所示，开展对气田储量动用状况评价:
[0075]从图12看出，纵向上飞三段动用较差，可动用层仅零星分布；飞一二段动用最好，但仍有少部分未动用层夹杂在动用层中；长兴组动用比飞三段好，比飞一二段差。在从南到北方向，普光104平台井区、2011平台区域动用较充分，边部普光304平台?305平台区域动用较差。从东到西储层动用整体较充分，未动用层、差层交错分布在动用层中间。
[0076]从图13_1_1、13_1_2、13_1_3看出，平面上，飞一二段未动用层、差层渗流特征参数值在0.1?2和5?10两个区间。其中，飞一二段上段渗流特征参数值在0.1?2之间的未动用层、差层主要分布在有利沉积相带鮞粒滩沉积中心的普光102-1?201-1平台区域；渗流特征参数值在5?10之间的未动用层、差层主要分布在构造相变线附近303平台附近；中段未动用层、差层分布与上段相似；下段未动用层、差层在构造高部位较集中分布。
[0077]从图13-2看出，平面上，飞一二段动用好的层渗流特征参数值<0.01。其中，飞一二段上段动用好的层主要分布在北部鮞粒滩沉积中心普光104平台?107-1H平台区域；中段动用好的层分布与上段相似；飞一二下段动用好的层较集中在构造高部位的普光302平台、301平台?201平台区域。
[0078]从图13-3看出，平面上，飞一二段动用中等的层渗流特征参数值介于0.01?0.1，分布范围较动用好的层扩大。其中，飞一二段上段动用中等的层主要分布普光302?202?203平台一线以北的区域；中段动用中等的层主要分布与上段相似，不同的是中段在普光107-1H还有发育；飞一二下段动用中等的层分布范围较动用好的层大幅度缩小，零星分布在普光302平台、201平台区域。
【权利要求】
1.一种碳酸盐气藏储量动用状况的评价方法，其特征包括以下步骤: (1)确定产层与非产层产剖特征； (2)确定储层的渗流特征参数； (3)利用储层的渗流特征参数，确定产层与非产层； (4)气藏储量动用状况评价。
2.根据权利要求1所述的一种碳酸盐气藏储量动用状况的评价方法，其特征是确定储层的渗流特征参数，其步骤如下: (1)建立产剖测试井的测井评价参数与产层、非产层特征的关系图版，确定出区分产层与非产层的评价单参数，包括补偿中子、深侧向电阻率参差比、含气饱和度、自然伽玛参数； (2)利用优选出的单参数进行多元组合，确定能定量区分产层与非产层的渗流特征参数。
3.根据权利要求1所述的一种碳酸盐气藏储量动用状况的评价方法，其特征是对优选出来的单参数进行多元组合，确定能定量评价储量动用状况的渗流特征参数Fs，其计算公式如下: c ( Riid/Ciid) *CNL^S, rs == -
GR 式中:Rlld/Clld—深侧向电阻率参差比，无单位； Riid—深侧向电阻率，Ω.Π1 ； Clld一深侧向电阻率参差系数Ω.m ; CNL—补偿中子，% ； GR—自然伽玛，API ； Sg—含气饱和度，％。
4.根据权利要求1、2或3所述的一种碳酸盐气藏储量动用状况的评价方法，其特征是利用储层的渗流特征参数Fs，确定产层与非产层，其步骤如下: (1)建立渗流参数Fs与产层、非产层的关系图版，确定出可动用层下限标准，区分出产层与非产层； (2)建立渗流特征参数Fs与米采气指数的关系图版，确定出产层优劣的分类标准。
5.根据权利要求4所述的一种碳酸盐气藏储量动用状况的评价方法，其特征是气藏储量动用状况评价，其步骤如下: (1)根据渗流特征参数Fs，确定各井点的渗流特征参数值Fs，在气藏孔隙度、渗透率等相关属性模型的三维约束下，采用地质建模的手段，建立储层的三维渗流特征参数连续地质模型； (2)通过对储层的三维渗流特征连续地质模型的约束控制、雕刻，得到气藏储量动用状况评价剖面图以及不同层段的气藏储量动用状况评价平面图，实现对碳酸盐气藏的储量动用状况评价。
【文档编号】E21B49/00GK104196523SQ201310533681
【公开日】2014年12月10日 申请日期:2013年10月31日 优先权日:2013年10月29日
【发明者】曾大乾, 张雪松, 毕建霞, 宿亚仙, 刘红磊, 张纪喜, 陈建 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司中原油田分公司勘探开发科学研究院