CN107092036A - 一种基于储层真电阻率反演的火山岩地层流体识别方法及*** - Google Patents
一种基于储层真电阻率反演的火山岩地层流体识别方法及*** Download PDFInfo
- Publication number
- CN107092036A CN107092036A CN201710230144.XA CN201710230144A CN107092036A CN 107092036 A CN107092036 A CN 107092036A CN 201710230144 A CN201710230144 A CN 201710230144A CN 107092036 A CN107092036 A CN 107092036A
- Authority
- CN
- China
- Prior art keywords
- resistivity
- reservoir
- mrow
- msub
- stratum
- 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.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 135
- 239000012530 fluid Substances 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 87
- 230000005251 gamma ray Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 3
- 239000011148 porous material Substances 0.000 abstract description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 9
- 238000011161 development Methods 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000005755 formation reaction Methods 0.000 description 81
- 239000007789 gas Substances 0.000 description 28
- 238000011160 research Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 10
- 230000009286 beneficial effect Effects 0.000 description 10
- 230000008859 change Effects 0.000 description 10
- 230000004044 response Effects 0.000 description 10
- 230000009467 reduction Effects 0.000 description 8
- 229910052770 Uranium Inorganic materials 0.000 description 7
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 7
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 6
- 239000002253 acid Substances 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 239000001257 hydrogen Substances 0.000 description 6
- 229910052739 hydrogen Inorganic materials 0.000 description 6
- 230000009466 transformation Effects 0.000 description 6
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 230000005611 electricity Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000007935 neutral effect Effects 0.000 description 4
- 239000002734 clay mineral Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000000280 densification Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 229910052683 pyrite Inorganic materials 0.000 description 2
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 2
- 239000011028 pyrite Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于储层真电阻率反演的火山岩地层流体识别方法及***,涉及石油开发中的储层评价领域。建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;识别火山岩地层岩性;根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;根据获得系数因子数值后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。本发明建立消除岩性、裂缝、孔隙结构影响的反演模型,进而构造火山岩流体性质识别参数,确定该参数在油层、气层和水层中的界限值,以达到在实际火山岩油气藏的评价中提高储层流体性质的识别精度。
Description
技术领域
本发明涉及石油开发中的储层评价领域,尤其涉及一种基于储层真电阻率反演的火山岩地层流体识别领域。
背景技术
最近十几年随着整个油气产业的不断深入与细化,油气勘探的方向也发生了巨大的转变,从构造油气藏向复杂岩性油气藏的转变,从浅层油气藏向深层致密油气藏的转变,从正常电性、物性油气藏向“三低”油气藏的转变,从碎屑岩油气藏向火山岩油气藏的转变等,勘探目标的转变就必须有相应的新理论和新技术的支撑。目前,火山岩油气藏的勘探与开发已经进入技术攻坚阶段并且获得了可观的地质储量,近几年我国先后在渤海湾、二连、黄骅、准噶尔等盆地不断发现了火山岩油气储层,显示了火山岩地层油气勘探开发的巨大潜力。但是火山岩油气藏岩性复杂多变,岩相变化剧烈,非均质性强以及双重孔隙结构(孔隙+裂缝)等都给火山岩油气层的有效识别带来了很大的困难,常规的识别方法在火山岩含气层段基本上没有气体的响应特征,识别结果难以得到保证。
发明内容
本发明为了解决常规识别方法所存在的缺陷,提供了一种基于储层真电阻率反演的火山岩地层流体识别方法及***。
本发明解决上述技术问题的技术方案如下:一种基于储层真电阻率反演的火山岩地层流体识别方法,所述方法包括:
建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;
识别火山岩地层岩性;
根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;
根据获得系数因子数值后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。
本发明的有益效果是:本发明为了能够从测井信息中有效的提取火山岩储层流体的响应特征,建立消除岩性、裂缝、孔隙结构影响的反演模型,进而构造火山岩流体性质识别参数,确定该参数在油层、气层和水层中的界限值,以达到在实际火山岩油气藏的评价中提高储层流体性质的识别精度。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述储层真电阻率反演模型为:
其中,Rt为基岩无裂缝致密地层电阻率,单位为Ω·m;Rlld为储层深侧向电阻率,单位为Ω·m;φf为裂缝孔隙度;φ为地层有效孔隙度;φN1为基岩致密地层中子孔隙度;φN为储层中子孔隙度;c1和c2为系数因子。
采用上述进一步方案的有益效果是:上述储层真电阻率反演模型是基于双重孔隙介质模型建立的,电阻率的校正值由两部分构成,包括裂缝部分及孔隙流体部分,利用该模型能将裂缝对电阻率的影响消除,计算消除裂缝因素后的地层真电阻率,同时也能够反映出储层中含有不同流体时,孔隙流体性质对电阻率的影响。
进一步,通过自然伽马曲线值识别火山岩地层岩性。
采用上述进一步方案的有益效果是:为了提高储层真电阻率反演模型的识别精度,系数因子c1和c2的反演必须基于相同岩性的致密围岩与含流体储层,通过对研究区火山岩岩性资料与不同岩性测井响应特征的综合分析,发现自然伽马测井值对岩性比较敏感,随着伽马值的增大,火山岩的岩性从基性向酸性变化。因此,选择自然伽马曲线值来识别火山岩岩性。
进一步,所述系数因子c1的反演过程为:
根据FMI成像测井计算裂缝孔隙度;
选取有裂缝井段的电阻率作为储层深侧向电阻率Rlld,选取没有裂缝井段的电阻率作为基岩无裂缝致密地层电阻率Rt;
将储层深侧向电阻率Rlld按照裂缝孔隙度的大小从0.01%至0.09%依次取平均值;
根据平均值绘制Rt/Rlld与φf/φ的交会图,从而最终确定不同岩性的高角度裂缝地层、中角度裂缝地层和低角度裂缝地层的系数因子c1的数值。
采用上述进一步方案的有益效果是:由于不同角度的裂缝对电阻率的影响程度不尽相同,对系数因子c1的影响也同样存在差异,因此,在细分岩性的基础上,选取研究区有试油结论的典型井,每种岩性按照高角度裂缝、中角度裂缝和低角度裂缝三种情况分别确定系数因子c1的数值。
进一步,所述系数因子c2的反演过程为:
选取含流体层段电阻率作为储层深侧向电阻率Rlld,该电阻率所对应的中子孔隙度为储层中子孔隙度φN;
选取干层井段电阻率作为基岩无裂缝致密地层电阻率Rt,该电阻率所对应的中子孔隙度为基岩致密地层中子孔隙度φN1;
绘制Rt/Rlld与1/φN1-1/φN的交会图,从而最终确定不同岩性的系数因子c2的数值。
采用上述进一步方案的有益效果是:在火山岩岩性相同的情况下,含氢指数可以有效的识别地层流体性质,而含氢指数在测井曲线上的反映就是中子孔隙度,选取研究区有试油结论的典型井,分别确定不同岩性的系数因子c2的数值。
进一步,所述流体识别参数为:
采用上述进一步方案的有益效果是:储层真电阻率反演模型中的基岩无裂缝致密地层电阻率与储层深侧向电阻率的差值与储层深侧向电阻率的比值为火山岩地层的流体识别参数,其中分子反映的是裂缝和流体对电阻率降低幅度的影响,分母则是消除非均质的影响。
进一步,通过绘制储层深侧向电阻率Rlld与流体识别参数FI的交会图确定不同储层流体识别参数界限值。
采用上述进一步方案的有益效果是:利用火山岩地层储层真电阻率反演模型处理研究区内有试油结论的生产井,制作储层深侧向电阻率与流体识别参数的交会图,进而确定火山岩不同流体性质储层的流体识别参数界限值。
为了解决上述技术问题,本发明还提出了一种基于储层真电阻率反演的火山岩地层流体识别***,所述***包括:
模型建立模块,用于建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;
岩性识别模块,用于识别火山岩地层岩性;
反演模块,用于根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;
界限值计算模块,用于根据获得系数因子数值后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。
本发明的有益效果是:本发明为了能够从测井信息中有效的提取火山岩储层流体的响应特征,建立消除岩性、裂缝、孔隙结构影响的反演模型,进而构造火山岩流体性质识别参数,确定该参数在油层、气层和水层中的界限值,以达到在实际火山岩油气藏的评价中提高储层流体性质的识别精度。
进一步,所述储层真电阻率反演模型为:
其中,Rt为基岩无裂缝致密地层电阻率,单位为Ω·m;Rlld为储层深侧向电阻率,单位为Ω·m;φf为裂缝孔隙度;φ为地层有效孔隙度;φN1为基岩致密地层中子孔隙度;φN为储层中子孔隙度;c1和c2为系数因子。
采用上述进一步方案的有益效果是:上述储层真电阻率反演模型是基于双重孔隙介质模型建立的,电阻率的校正值由两部分构成,包括裂缝部分及孔隙流体部分,利用该模型能将裂缝对电阻率的影响消除,计算消除裂缝因素后的地层真电阻率,同时也能够反映出储层中含有不同流体时,孔隙流体性质对电阻率的影响。
进一步,所述岩性识别模块用于通过自然伽马曲线值识别火山岩地层岩性。
采用上述进一步方案的有益效果是:为了提高储层真电阻率反演模型的识别精度,系数因子c1和c2的反演必须基于相同岩性的致密围岩与含流体储层,通过对研究区火山岩岩性资料与不同岩性测井响应特征的综合分析,发现自然伽马测井值对岩性比较敏感,随着伽马值的增大,火山岩的岩性从基性向酸性变化。因此,选择自然伽马曲线值来识别火山岩岩性。
附图说明
图1为实施例1所述的火山岩地层流体识别方法流程图;
图2为英安岩高角度裂缝地层Rt/Rlld与φf/φ交会图;
图3为安山岩地层Rt/Rlld与1/φN1-1/φN交会图;
图4为实施例2所述的火山岩地层流体识别***原理图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
在上述背景技术中提及了目前常规的识别方法在火山岩含气层段基本上没有气体的响应特征,导致识别结果容易出现偏差的问题。基于上述问题,本发明提出如下具体实施例。
实施例1
如图1所示,本实施例提出了一种基于储层真电阻率反演的火山岩地层流体识别方法,所述方法包括:
建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;
识别火山岩地层岩性;
根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;
根据获得系数因子数值后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。
所述方法能够消除岩性及不同角度裂缝对火山岩油气藏流体性质识别的影响,提高储层流体性质评价的精度,具体的,所述储层真电阻率反演模型为:
其中,Rt为基岩无裂缝致密地层电阻率,单位为Ω·m;Rlld为储层深侧向电阻率,单位为Ω·m;φf为裂缝孔隙度;φ为地层有效孔隙度;φN1为基岩致密地层中子孔隙度;φN为储层中子孔隙度;c1和c2为系数因子。
本实施例所述方法为了从测井信息中有效的提取火山岩储层流体的响应特征,优选能够反映储层流体性质变化的测井参数,建立消除岩性、裂缝、孔隙结构影响的反演模型,进而准确判断储层流体类型。
电阻率是油气储层的一个重要参数,能够反映储层流体性质及不同流体饱和度的变化。但是,影响储层电阻率的因素又比较多,在基岩电阻率较高的硬地层中,岩性、裂缝、泥质含量、黏土矿物类型、流体性质、导电矿物等都可以导致储层电阻率的变化。在火山岩地层中出去泥质的影响,对于岩性相同的地层,储层电阻率与致密围岩电阻率的差异可以认为是裂缝及孔隙流体性质引起的。以准噶尔盆地陆东五彩湾地区石炭系火山岩油气藏为例,主要含流体地层的岩性为中性安山岩和中酸性英安岩,矿物成分分析表明地层中泥质含量非常少,并且没有黄铁矿等导电矿物,因此,导致火山岩地层电阻率降低的主要因素是裂缝及孔隙流体。
因此,火山岩地层流体性质的识别优选深侧向电阻率作为反演参数,在细分岩性的基础上构造一个消除裂缝和流体性质影响的地层真电阻率反演模型。
从上述储层真电阻率反演模型中可以看出,电阻率的校正值由两部分构成,包括裂缝部分以及孔隙流体部分,下面具体分析一下裂缝及孔隙流体性质对于电阻率的影响。
裂缝对电阻率影响分析:从上述储层真电阻率反演模型可以看出,地层中有裂缝时就会导致电阻率降低,电阻率降低幅度与裂缝孔隙度有关系,当裂缝孔隙度等于0时,即φf=0时,地层真电阻率等于储层深侧向电阻率Rlld;当裂缝孔隙度不等于0时,随着裂缝孔隙度增大,电阻率降低的幅度也越大,利用上述储层真电阻率反演模型则可以将裂缝对电阻率的影响消除,计算获得消除裂缝因素后的地层真电阻率。
孔隙流体性质对电阻率影响分析:孔隙流体校正的关键是流体性质的识别,在火山岩岩性相同也就是地层骨架对快中子减速能力基本相近的情况下,含氢指数可以有效地识别地层流体性质,当火山岩储层含有油、气、水时,由于各种流体含氢指数不同,对快中子的减速能力也不同,所以中子孔隙度比干层大,电阻率比致密围岩小。当储层是干层时,致密地层中子孔隙度等于储层中子孔隙度,即φN1=φN,基岩无裂缝致密地层电阻率等于储层深侧向电阻率;当储层中含有不同流体时,其对电阻率的降低幅度是不同的。
为了提高储层真电阻率反演模型的识别精度,系数因子c1和c2的反演必须基于相同岩性的致密围岩与含流体储层,通过对研究区火山岩岩性资料与不同岩性测井响应特征的综合分析,发现自然伽马测井值对岩性比较敏感,随着伽马值的增大,火山岩的岩性从基性向酸性变化,如表1所示。因此,选择自然伽马曲线值来识别火山岩岩性。
表1
在识别火山岩岩性之后,对系数因子c1和c2进行反演。
所述系数因子c1的反演过程为:
根据FMI成像测井计算裂缝孔隙度;
选取有裂缝井段的电阻率作为储层深侧向电阻率Rlld,选取没有裂缝井段的电阻率作为基岩无裂缝致密地层电阻率Rt;
将储层深侧向电阻率Rlld按照裂缝孔隙度的大小从0.01%至0.09%依次取平均值;
根据平均值绘制Rt/Rlld与φf/φ的交会图,从而最终确定不同岩性的高角度裂缝地层、中角度裂缝地层和低角度裂缝地层的系数因子c1的数值。
具体的,由于裂缝角度对系数因子c1影响较大,因此,系数因子c1的反演分高角度裂缝(≥75%)、中角度裂缝(45%至75%)和低角度裂缝(≤45%)三种情况进行。以高角度裂缝系数因子c1的反演为例,致密围岩和裂缝层段电阻率的选取都在英安岩地层中确定,裂缝孔隙度依据FMI成像测井计算,选取有裂缝井段的电阻率作为储层深侧向电阻率Rlld,选取没有裂缝井段的电阻率作为基岩无裂缝致密地层电阻率Rt;把储层深侧向电阻率Rlld按照裂缝孔隙度的大小从0.01%至0.09%依次取平均值,进而制作Rt/Rlld与φf/φ交会图,如图2所示,由交会图可以确定英安岩高角度裂缝地层系数因子c1的数值,c1=120.65,利用相同的方法可以确定英安岩中角度和低角度裂缝地层系数因子c1的数值,分别为c1=124.52和c1=136.32。由于研究区火山岩储层的岩性主要以中性安山岩和中酸性英安岩为主,按照英安岩地层系数因子c1的数值的反演原理,同样可以确定安山岩高角度、中角度和低角度裂缝地层系数因子c1的数值,分别为c1=118.27、c1=122.86和c1=135.16。
所述系数因子c2的反演过程为:
选取含流体层段电阻率作为储层深侧向电阻率Rlld,该电阻率所对应的中子孔隙度为储层中子孔隙度φN;
选取干层井段电阻率作为基岩无裂缝致密地层电阻率Rt,该电阻率所对应的中子孔隙度为基岩致密地层中子孔隙度φN1;
绘制Rt/Rlld与1/φN1-1/φN的交会图,从而最终确定不同岩性的系数因子c2的数值。
具体的,利用研究区典型井进行系数因子c2的反演,以安山岩地层为例,选取含流体层段电阻率作为储层深侧向电阻率Rlld,也就是由于流体的存在导致地层电阻率降低值,对应的中子孔隙度为储层中子孔隙度φN;选取干层井段电阻率作为基岩无裂缝致密地层电阻率Rt,对应的中子孔隙度为基岩致密地层中子孔隙度φN1;绘制Rt/Rlld与1/φN1-1/φN的交会图,如图3所示。由交会图可以确定安山岩地层系数因子c2的数值,c2=80.08。采用相同的原理,根据自然伽马测井值选取英安岩地层进行系数因子c2的反演,c2=73.68。
由于准噶尔盆地陆东五彩湾地区石炭系火山岩含流体地层的岩性主要以中性安山岩和中酸性英安岩为主,因此,地层真电阻率反演中系数因子c1和c2的选取如表2所示。
表2
结合表1中的岩性识别标准以及储层真电阻率反演模型,可以计算消除裂缝的流体性质影响后的火山岩地层真电阻率,模型反演的真电阻率与实际测量的深侧向电阻率的差值就是裂缝及孔隙流体所引起的,因此,获得流体识别参数即真电阻率与深侧向电阻率的差值与深侧向电阻率的比值,分子反映的是裂缝和流体对电阻率降低幅度的影响,分母则是消除非均质的影响。
采用火山岩地层储层真电阻率反演模型及方法处理研究区内有试油结论的生产井,制作储层深侧向电阻率Rlld与流体识别参数FI的交会图。火山岩储层由于非均质性严重及不同角度裂缝的分布导致深侧向电阻率变化比较大,流体性质对电阻率的影响远小于前两者对电阻率的贡献。因此,利用深侧向电阻率无法准确的识别火山岩地层不同流体性质的储层。流体识别参数消除了岩性、非均质性和裂缝对电阻率的影响,把不同流体性质对电阻率的贡献幅度提取出来,可以有效地识别火山岩裂缝地层不同流体性质的储层,其中油气同层的界限值大于0.3,气层为0.1至0.3,气水同层为0.035至0.1,水层小于0.035。
利用本实例建立的地层真电阻率反演模型及不同流体性质储层的识别标准评价准噶尔盆地陆东五彩湾地区火山岩裂缝地层的流体性质,取得了比较好的应用效果。从识别结果的综合分析中可以看出,在碎屑岩中识别气层比较好的三孔隙度曲线在火山岩裂缝地层基本上没有气层的响应特征;另外,相同岩性的电阻率曲线差异比较大,复杂岩性以及裂缝的发育对电阻率的贡献远大于流体性质,因此无法直接从深侧向电阻率曲线中提取不同流体的信息。利用真电阻率反演模型及流体性质识别界限值对该区火山岩地层进行评价,可以有效地识别含流体储层,并且能够准确判断储层的流体性质,为火山岩油气藏的后续开发提供了重要的地质依据。
实施例2
如图4所示,本实施例提出了一种基于储层真电阻率反演的火山岩地层流体识别***,所述***包括:
模型建立模块,用于建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;
岩性识别模块,用于识别火山岩地层岩性;
反演模块,用于根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;
界限值计算模块,用于根据获得系数因子数值后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。
优选的,所述储层真电阻率反演模型为:
其中,Rt为基岩无裂缝致密地层电阻率,单位为Ω·m;Rlld为储层深侧向电阻率,单位为Ω·m;φf为裂缝孔隙度;φ为地层有效孔隙度;φN1为基岩致密地层中子孔隙度;φN为储层中子孔隙度;c1和c2为系数因子。
电阻率是油气储层的一个重要参数,能够反映储层流体性质及不同流体饱和度的变化。但是,影响储层电阻率的因素又比较多,在基岩电阻率较高的硬地层中,岩性、裂缝、泥质含量、黏土矿物类型、流体性质、导电矿物等都可以导致储层电阻率的变化。在火山岩地层中出去泥质的影响,对于岩性相同的地层,储层电阻率与致密围岩电阻率的差异可以认为是裂缝及孔隙流体性质引起的。以准噶尔盆地陆东五彩湾地区石炭系火山岩油气藏为例,主要含流体地层的岩性为中性安山岩和中酸性英安岩,矿物成分分析表明地层中泥质含量非常少,并且没有黄铁矿等导电矿物,因此,导致火山岩地层电阻率降低的主要因素是裂缝及孔隙流体。
因此,火山岩地层流体性质的识别优选深侧向电阻率作为反演参数,在细分岩性的基础上构造一个消除裂缝和流体性质影响的地层真电阻率反演模型。
从上述储层真电阻率反演模型中可以看出,电阻率的校正值由两部分构成,包括裂缝部分以及孔隙流体部分,下面具体分析一下裂缝及孔隙流体性质对于电阻率的影响。
裂缝对电阻率影响分析:从上述储层真电阻率反演模型可以看出,地层中有裂缝时就会导致电阻率降低,电阻率降低幅度与裂缝孔隙度有关系,当裂缝孔隙度等于0时,即φf=0时,地层真电阻率等于储层深侧向电阻率Rlld;当裂缝孔隙度不等于0时,随着裂缝孔隙度增大,电阻率降低的幅度也越大,利用上述储层真电阻率反演模型则可以将裂缝对电阻率的影响消除,计算获得消除裂缝因素后的地层真电阻率。
孔隙流体性质对电阻率影响分析:孔隙流体校正的关键是流体性质的识别,在火山岩岩性相同也就是地层骨架对快中子减速能力基本相近的情况下,含氢指数可以有效地识别地层流体性质,当火山岩储层含有油、气、水时,由于各种流体含氢指数不同,对快中子的减速能力也不同,所以中子孔隙度比干层大,电阻率比致密围岩小。当储层是干层时,致密地层中子孔隙度等于储层中子孔隙度,即φN1=φN,基岩无裂缝致密地层电阻率等于储层深侧向电阻率;当储层中含有不同流体时,其对电阻率的降低幅度是不同的。
为了提高储层真电阻率反演模型的识别精度,系数因子c1和c2的反演必须基于相同岩性的致密围岩与含流体储层,通过对研究区火山岩岩性资料与不同岩性测井响应特征的综合分析,发现自然伽马测井值对岩性比较敏感,随着伽马值的增大,火山岩的岩性从基性向酸性变化。因此,选择自然伽马曲线值来识别火山岩岩性。
在识别火山岩岩性之后,需要针对具体的储层确定系数因子c1和c2的数值。
优选的,所述系数因子c1的反演过程为:
根据FMI成像测井计算裂缝孔隙度;
选取有裂缝井段的电阻率作为储层深侧向电阻率Rlld,选取没有裂缝井段的电阻率作为基岩无裂缝致密地层电阻率Rt;
将储层深侧向电阻率Rlld按照裂缝孔隙度的大小从0.01%至0.09%依次取平均值;
根据平均值绘制Rt/Rlld与φf/φ的交会图,从而最终确定不同岩性的高角度裂缝地层、中角度裂缝地层和低角度裂缝地层的系数因子c1的数值。
具体的,由于裂缝角度对系数因子c1影响较大,因此,系数因子c1的反演分高角度裂缝(≥75%)、中角度裂缝(45%至75%)和低角度裂缝(≤45%)三种情况进行。以高角度裂缝系数因子c1的反演为例,致密围岩和裂缝层段电阻率的选取都在英安岩地层中确定,裂缝孔隙度依据FMI成像测井计算,选取有裂缝井段的电阻率作为储层深侧向电阻率Rlld,选取没有裂缝井段的电阻率作为基岩无裂缝致密地层电阻率Rt;把储层深侧向电阻率Rlld按照裂缝孔隙度的大小从0.01%至0.09%依次取平均值,进而制作Rt/Rlld与φf/φ交会图,如图2所示,由交会图可以确定英安岩高角度裂缝地层系数因子c1的数值,c1=120.65,利用相同的方法可以确定英安岩中角度和低角度裂缝地层系数因子c1的数值,分别为c1=124.52和c1=136.32。由于研究区火山岩储层的岩性主要以中性安山岩和中酸性英安岩为主,按照英安岩地层系数因子c1的数值的反演原理,同样可以确定安山岩高角度、中角度和低角度裂缝地层系数因子c1的数值,分别为c1=118.27、c1=122.86和c1=135.16。
优选的,所述系数因子c2的反演过程为:
选取含流体层段电阻率作为储层深侧向电阻率Rlld,该电阻率所对应的中子孔隙度为储层中子孔隙度φN;
选取干层井段电阻率作为基岩无裂缝致密地层电阻率Rt,该电阻率所对应的中子孔隙度为基岩致密地层中子孔隙度φN1;
绘制Rt/Rlld与1/φN1-1/φN的交会图,从而最终确定不同岩性的系数因子c2的数值。
具体的,利用研究区典型井进行系数因子c2的反演,以安山岩地层为例,选取含流体层段电阻率作为储层深侧向电阻率Rlld,也就是由于流体的存在导致地层电阻率降低值,对应的中子孔隙度为储层中子孔隙度φN;选取干层井段电阻率作为基岩无裂缝致密地层电阻率Rt,对应的中子孔隙度为基岩致密地层中子孔隙度φN1;绘制Rt/Rlld与1/φN1-1/φN的交会图,如图3所示。由交会图可以确定安山岩地层系数因子c2的数值,c2=80.08。采用相同的原理,根据自然伽马测井值选取英安岩地层进行系数因子c2的反演,c2=73.68。
由于准噶尔盆地陆东五彩湾地区石炭系火山岩含流体地层的岩性主要以中性安山岩和中酸性英安岩为主,因此,地层真电阻率反演中系数因子c1和c2的选取如上述表2。
结合表1中的岩性识别标准以及储层真电阻率反演模型,可以计算消除裂缝的流体性质影响后的火山岩地层真电阻率,模型反演的真电阻率与实际测量的深侧向电阻率的差值就是裂缝及孔隙流体所引起的,因此,获得流体识别参数即真电阻率与深侧向电阻率的差值与深侧向电阻率的比值,分子反映的是裂缝和流体对电阻率降低幅度的影响,分母则是消除非均质的影响。
采用火山岩地层储层真电阻率反演模型及方法处理研究区内有试油结论的生产井,制作储层深侧向电阻率Rlld与流体识别参数FI的交会图。火山岩储层由于非均质性严重及不同角度裂缝的分布导致深侧向电阻率变化比较大,流体性质对电阻率的影响远小于前两者对电阻率的贡献。因此,利用深侧向电阻率无法准确的识别火山岩地层不同流体性质的储层。流体识别参数消除了岩性、非均质性和裂缝对电阻率的影响,把不同流体性质对电阻率的贡献幅度提取出来,可以有效地识别火山岩裂缝地层不同流体性质的储层,其中油气同层的界限值大于0.3,气层为0.1至0.3,气水同层为0.035至0.1,水层小于0.035。
本实施例所述的火山岩地层流体识别***在分析火山岩地层电阻率影响因素的基础上,建立了消除裂缝角度及孔隙流体性质的真电阻率反演模型,利用真电阻率与深侧向电阻率的差比值有效的反映出不同流体性质的储层在电阻率曲线上的差异,不仅仅提高了火山岩油藏储层流体性质识别的精度,为油、气储层的开发提供了地质依据。对于不同的火山岩油气藏,只要针对具体的储层确定系数因子c1和c2的数值,就可以应用到该类火山岩油气藏储层流体性质的识别中。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,所述方法包括:
建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;
识别火山岩地层岩性;
根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;
根据获得系数因子后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。
2.根据权利要求1所述的一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,所述储层真电阻率反演模型为:
<mrow>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>*</mo>
<mfrac>
<msub>
<mi>&phi;</mi>
<mi>f</mi>
</msub>
<mi>&phi;</mi>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>N</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&phi;</mi>
<mi>N</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>;</mo>
</mrow>
其中,Rt为基岩无裂缝致密地层电阻率,单位为Ω·m;Rlld为储层深侧向电阻率,单位为Ω·m;φf为裂缝孔隙度;φ为地层有效孔隙度;φN1为基岩致密地层中子孔隙度;φN为储层中子孔隙度;c1和c2为系数因子。
3.根据权利要求2所述的一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,通过自然伽马曲线值识别火山岩地层岩性。
4.根据权利要求2或3所述的一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,所述系数因子c1的反演过程为:
根据FMI成像测井计算裂缝孔隙度;
选取有裂缝井段的电阻率作为储层深侧向电阻率Rlld,选取没有裂缝井段的电阻率作为基岩无裂缝致密地层电阻率Rt;
将储层深侧向电阻率Rlld按照裂缝孔隙度的大小从0.01%至0.09%依次取平均值;
根据平均值绘制Rt/Rlld与φf/φ的交会图,从而最终确定不同岩性的高角度裂缝地层、中角度裂缝地层和低角度裂缝地层的系数因子c1的数值。
5.根据权利要求4所述的一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,所述系数因子c2的反演过程为:
选取含流体层段电阻率作为储层深侧向电阻率Rlld,该电阻率所对应的中子孔隙度为储层中子孔隙度φN;
选取干层井段电阻率作为基岩无裂缝致密地层电阻率Rt,该电阻率所对应的中子孔隙度为基岩致密地层中子孔隙度φN1;
绘制Rt/Rlld与1/φN1-1/φN的交会图,从而最终确定不同岩性的系数因子c2的数值。
6.根据权利要求5所述的一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,所述流体识别参数为:
<mrow>
<mi>F</mi>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
</mfrac>
<mo>.</mo>
</mrow>
7.根据权利要求6所述的一种基于储层真电阻率反演的火山岩地层流体识别方法,其特征在于,通过绘制储层深侧向电阻率Rlld与流体识别参数FI的交会图确定不同储层流体识别参数界限值。
8.一种基于储层真电阻率反演的火山岩地层流体识别***,其特征在于,所述***包括:
模型建立模块,用于建立储层真电阻率反演模型,所述储层真电阻率反演模型包含系数因子;
岩性识别模块,用于识别火山岩地层岩性;
反演模块,用于根据识别的火山岩地层岩性所对应的裂缝孔隙度和中子孔隙度反演获得系数因子数值;
界限值计算模块,用于根据获得系数因子数值后的储层真电阻率反演模型获得流体识别参数,并确定不同储层流体识别参数界限值。
9.根据权利要求8所述的一种基于储层真电阻率反演的火山岩地层流体识别***,其特征在于,所述储层真电阻率反演模型为:
<mrow>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>*</mo>
<mfrac>
<msub>
<mi>&phi;</mi>
<mi>f</mi>
</msub>
<mi>&phi;</mi>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>N</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&phi;</mi>
<mi>N</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>;</mo>
</mrow>
其中,Rt为基岩无裂缝致密地层电阻率,单位为Ω·m;Rlld为储层深侧向电阻率,单位为Ω·m;φf为裂缝孔隙度;φ为地层有效孔隙度;φN1为基岩致密地层中子孔隙度;φN为储层中子孔隙度;c1和c2为系数因子。
10.根据权利要求8或9所述的一种基于储层真电阻率反演的火山岩地层流体识别***,其特征在于,所述岩性识别模块用于通过自然伽马曲线值识别火山岩地层岩性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710230144.XA CN107092036A (zh) | 2017-04-10 | 2017-04-10 | 一种基于储层真电阻率反演的火山岩地层流体识别方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710230144.XA CN107092036A (zh) | 2017-04-10 | 2017-04-10 | 一种基于储层真电阻率反演的火山岩地层流体识别方法及*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107092036A true CN107092036A (zh) | 2017-08-25 |
Family
ID=59636918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710230144.XA Pending CN107092036A (zh) | 2017-04-10 | 2017-04-10 | 一种基于储层真电阻率反演的火山岩地层流体识别方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107092036A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861162A (zh) * | 2017-09-15 | 2018-03-30 | 北京旭日昌盛科技有限公司 | 基于微电极测井资料的天然裂缝识别方法及*** |
CN108830140A (zh) * | 2018-04-28 | 2018-11-16 | 中国石油大学(华东) | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
CN110873904A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
CN111827989A (zh) * | 2020-07-20 | 2020-10-27 | 中国石油天然气集团有限公司 | 一种沥青削减指数识别含沥青砾岩油层的方法 |
CN113031109A (zh) * | 2019-12-25 | 2021-06-25 | 中国石油天然气股份有限公司 | 岩石物理建模方法及装置 |
CN114427454A (zh) * | 2020-09-09 | 2022-05-03 | 中国石油化工股份有限公司 | 受泥浆侵入影响较大储层的真实地层电阻率矫正方法 |
-
2017
- 2017-04-10 CN CN201710230144.XA patent/CN107092036A/zh active Pending
Non-Patent Citations (3)
Title |
---|
康志勇: "《油气藏参数评价方法》", 30 April 2013, 石油工业出版社 * |
谭锋奇: "构造电阻率差比值法识别火山岩裂缝地层天然气层", 《吉林大学学报(地球科学版)》 * |
高成军: "《碳酸盐岩储层测井与录井评价技术》", 28 February 2007, 石油工业出版社 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861162A (zh) * | 2017-09-15 | 2018-03-30 | 北京旭日昌盛科技有限公司 | 基于微电极测井资料的天然裂缝识别方法及*** |
CN107861162B (zh) * | 2017-09-15 | 2019-12-24 | 北京旭日昌盛科技有限公司 | 基于微电极测井资料的天然裂缝识别方法及*** |
CN108830140A (zh) * | 2018-04-28 | 2018-11-16 | 中国石油大学(华东) | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
CN110873904A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
CN110873904B (zh) * | 2018-09-03 | 2022-02-01 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
CN113031109A (zh) * | 2019-12-25 | 2021-06-25 | 中国石油天然气股份有限公司 | 岩石物理建模方法及装置 |
CN113031109B (zh) * | 2019-12-25 | 2024-02-02 | 中国石油天然气股份有限公司 | 岩石物理建模方法及装置 |
CN111827989A (zh) * | 2020-07-20 | 2020-10-27 | 中国石油天然气集团有限公司 | 一种沥青削减指数识别含沥青砾岩油层的方法 |
CN111827989B (zh) * | 2020-07-20 | 2023-08-22 | 中国石油天然气集团有限公司 | 一种沥青削减指数识别含沥青砾岩油层的方法 |
CN114427454A (zh) * | 2020-09-09 | 2022-05-03 | 中国石油化工股份有限公司 | 受泥浆侵入影响较大储层的真实地层电阻率矫正方法 |
CN114427454B (zh) * | 2020-09-09 | 2024-05-31 | 中国石油化工股份有限公司 | 受泥浆侵入影响较大储层的真实地层电阻率矫正方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107092036A (zh) | 一种基于储层真电阻率反演的火山岩地层流体识别方法及*** | |
CN103867197B (zh) | 复杂岩性天然气层声波时差判别法 | |
CN109838230B (zh) | 油藏水淹层的定量评价方法 | |
CN103603659B (zh) | 利用常规测井资料识别复杂孔隙结构储层流体类型的方法 | |
CN110847901B (zh) | 一种变矿化度地层水下的致密砂岩储层流体识别方法 | |
CN104863574B (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN101980053A (zh) | 一种复杂礁滩储层预测方法 | |
CN105525917B (zh) | 识别储层的方法以及识别储层中流体类型的方法 | |
CN105275456A (zh) | 一种利用测井资料识别优质泥页岩的方法 | |
CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
CN110554441A (zh) | 一种基于钙质含量校正的致密砂岩含气饱和度计算方法 | |
CN103867198A (zh) | 碳酸盐岩天然气层地层密度判别法 | |
CN107515957A (zh) | 泥页岩层序地层划分方法 | |
CN111008451A (zh) | 一种快速识别油气层的测井解释方法 | |
CN112016753A (zh) | 基于三元耦合的变质岩潜山产能预测方法 | |
Zeng et al. | Tectonic stress field and its effect on hydrocarbon migration and accumulation in Mesozoic and Cenozoic in Kuqa depression, Tarim basin | |
Gao et al. | Fine characterization of large composite channel sandbody architecture and its control on remaining oil distribution: a case study of alkaline-surfactant-polymer (ASP) flooding test area in Xingshugang oilfield, China | |
CN114139242A (zh) | 基于岩相的水淹层测井评价方法 | |
Wang | Study on reservoir heterogeneity in block S | |
CN110688781A (zh) | 一种低渗非均质气藏储层测井解释方法 | |
Feng et al. | Lithology and oil-bearing properties of tight sandstone reservoirs: Chang 7 member of Upper Triassic Yanchang Formation, southwestern Ordos Basin, China | |
CN112784404A (zh) | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 | |
CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 | |
CN106990448A (zh) | 一种海相碎屑岩水平井隔夹层测井解释方法及装置 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170825 |