CN102174888B - 一种基于储集岩岩性参数的地层数据处理方法 - Google Patents
一种基于储集岩岩性参数的地层数据处理方法 Download PDFInfo
- Publication number
- CN102174888B CN102174888B CN2011100504160A CN201110050416A CN102174888B CN 102174888 B CN102174888 B CN 102174888B CN 2011100504160 A CN2011100504160 A CN 2011100504160A CN 201110050416 A CN201110050416 A CN 201110050416A CN 102174888 B CN102174888 B CN 102174888B
- Authority
- CN
- China
- Prior art keywords
- oil
- volume
- reservoir
- rock
- saturation
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 130
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title abstract description 35
- 238000012545 processing Methods 0.000 title abstract description 3
- 239000012530 fluid Substances 0.000 claims abstract description 35
- 239000011148 porous material Substances 0.000 claims abstract description 24
- 238000011161 development Methods 0.000 claims abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 67
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 claims description 25
- 238000003672 processing method Methods 0.000 claims description 11
- 102000003712 Complement factor B Human genes 0.000 claims description 10
- 108090000056 Complement factor B Proteins 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 abstract description 19
- 238000011156 evaluation Methods 0.000 abstract description 7
- 239000003208 petroleum Substances 0.000 abstract description 5
- 239000003921 oil Substances 0.000 description 231
- 238000004458 analytical method Methods 0.000 description 22
- 229920006395 saturated elastomer Polymers 0.000 description 10
- 239000000460 chlorine Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 239000008398 formation water Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013213 extrapolation Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 3
- 229910052753 mercury Inorganic materials 0.000 description 3
- 239000003129 oil well Substances 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000009736 wetting Methods 0.000 description 2
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000007872 degassing Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000011888 foil Substances 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000003350 kerosene Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供一种基于储集岩岩性参数的地层数据处理方法,其属于储集岩岩性参数和地层数据的地球物理勘探与测井评价开发领域,包括利用地质勘探及开发设备采集地层岩性和流体信息;通过分析设备对采集的地层岩性及流体信息进行分析获得计算储集岩孔隙流体体积所需的基础参数,通过选取相应孔隙度和含油饱和度等储层参数,建立了储集岩孔隙中石油体积和石油地质储量计算方程。本发明适用于探明储量、控制储量和预测储量的储集岩石油资源评价,属于容积法储量计算范畴。本发明建立的容积法储量计算方程,消除了原容积法储量计算公式本身造成的***误差,使石油地质储量计算精度提高了6%~21%。
Description
技术领域
本发明涉及地球物理勘探和测井评价开发方法,具体为一种基于储集岩岩性参数的地层数据处理方法,更进一步地通过基于储集岩岩性参数的地层数据处理来计算获得储集岩孔隙流体体积。
背景技术
石油是储集岩孔隙流体中的一种,计算石油地质储量有多种方法,其中计算的地质储量包括探明地质储量、控制地质储量和预测地质储量。探明地质储量是指在油藏评价阶段,经评价钻探证实油藏可提供开采并能获得经济效益后,估算求得的、确定性很大的地质储量。此外,客观准确地评价油藏地质储量规模一直是油田勘探开发科技人员要解决的技术难题。
目前计算油藏石油地质储量主要有三种方法:一是容积法,根据《石油天然气储量计算规范》提供的公式计算油藏石油地质储量;二是动态法,通过对生产一定阶段的油藏动态信息测试分析,获取地层孔隙流体的各类信息及地层温度和地层压力信息,依据对获取信息的规律性分析,采用相应技术方法计算油藏石油地质储量;三是概率法,通过获取的地质信息分析获取油藏石油地质储量。
众所周知,概率法确定石油地质储量主要应用于预测地质储量,精度要求低,应用范围有限;动态法计算油藏石油地质储量要求油藏的生产动态资料能反映整个油藏的生产能力,油藏生产到中后期时计算的石油地质储量才具有代表性。容积法计算石油地质储量适用于勘探开发的各个阶段,特别是油藏评价阶段和期其以后的各阶段。,例如在中国专利申请CN101937108A中使用的容积法储量计算公式N=100AohφSoi/Boi,其要求计算的是理论静态石油地质储量,但是其计算的石油地质储量N是油藏实际石油地质储量的(100φ/φt)%或(100Soi/Soe)%,与实际石油地质储量存在***误差。《石油天然气储量计算规范》是经过专家学者多年努力总结的研究成果,也是勘探开发各个阶段都适用的油藏地质储量评价技术,但在应用中发现规范中的石油储量计算公式仍存在一些问题,需要进一步完善和改进,如地层有效孔隙度φ与原始含油饱和度Soi的乘积并不等于储集岩孔隙中的石油体积百分比,再乘以含油面积和油层有效厚度,也不等于地层条件下储集岩孔隙中的石油体积,而是小于储集岩孔隙中的石油体积;由于应用习惯和管理上的周期性等诸多原因,目前一直在延用这一尚不完善的容积法石油地质储量计算公式。
综上所述,为了准确计算油藏石油地质储量,需要对现有技术中提供的容积法储量计算公式进行修正,进而克服因计算方法而造成的计算石油地质储量的***误差。
发明内容
为了解决上面背景技术中存在的问题,本发明提出一种基于储集岩岩性参数的地层数据处理方法,更进一步地通过基于储集岩岩性参数的地层数据处理来计算获得储集岩孔隙流体体积。
依据本发明提出的技术方案,提供一种基于储集岩岩性参数的地层数据处理方法,其包括以下步骤:
1)利用地质勘探开发设备,可以采集地层岩性、孔隙流体及测井信息;
2)基于采集到的地层岩性数据,通过岩性数据分析设备获得储集岩的地层总孔隙度φt、地层有效孔隙度φ、束缚水饱和度Swi、有效含水饱和度Swe、有效束缚水饱和度Swie和含水饱和度Sw;
3)基于采集到的孔隙流体和测井信息,结合岩性数据分析设备获得油层有效厚度h和含油面积Ao;
4)基于采集到的地层孔隙流体信息,通过流体分析设备获得原始原油体积系数Boi;
5)根据在步骤2)中获得的储集岩地层总孔隙度φt、地层有效孔隙度φ、束缚水饱和度Swi、有效含水饱和度Swe、有效束缚水饱和度Swie和含水饱和度Sw,来获得储集岩孔隙中的石油体积百分比OR(单位f);
其中,OR为石油体积百分比,指岩样孔隙中的石油所占据的孔隙体积与岩样总体积的比值,单位f;
6)根据得到的石油体积百分比OR和在步骤3)中获得的油层有效厚度h(单位m)和含油面积Ao(单位km2),来获得含油储集岩孔隙中的石油体积Vf(单位104m3);
7)根据得到的石油体积Vf和在步骤4)中获得的原始原油体积系数Boi(无因次),来获得油藏石油地质储量N′(单位104m3),石油地质储量N′分别为下面几种情形;
对于纯油层:N′=100AohφtSoi/Boi;或N′=100AohφSoie/Boi;
对于含可动水油层:N′=100AohφtSo/Boi;或N′=100AohφSoe/Boi。
优选地,石油体积百分比OR分别为下列情况:
纯油层为:OR=φt(1-Swi);或OR=φ(1-Swie);或OR=φtSoi;或OR=φSoie;
含可动水油层为:OR=φt(1-Sw);或OR=φ(1-Swe);或OR=φtSo;或OR=φSoe;
其中,φ为地层有效孔隙度(指岩样的有效孔隙体积与岩样总体积的比值),单位f;φt为地层总孔隙度(指岩样的总孔隙体积与岩样总体积的比值),单位f;Soi为在纯油层中对应地层总孔隙度φt的原始含油饱和度(在纯油层中,指岩样中石油所占据的孔隙体积与岩样总孔隙体积的比值),并有Soi=1-Swi,单位f;Soie为在纯油层中对应地层有效孔隙度φ的有效原始含油饱和度(在纯油层中,指岩样中石油所占据的孔隙体积与岩样有效孔隙体积的比值),并有Soie=1-Swie,单位f;So为在含可动水油层中对应地层总孔隙度φt的含油饱和度(在含可动水油层中,指岩样中石油所占据的孔隙体积与岩样总孔隙体积的比值),并有So=1-Sw,单位f;Soe为在含可动水油层中对应地层有效孔隙度φ的有效含油饱和度(在含可动水油层中,指岩样中石油所占据的孔隙体积与岩样有效孔隙体积的比值),并有Soe=1-Swe,单位f。
更优选地,储集岩孔隙中的石油体积Vf分别为下列情况:纯油层时为Vf=100AohφtSoi;或Vf=100AohφSoie;含可动水油层时为Vf=100AohφtSo;或Vf=100AohφSoe。
使用本发明所述的基于储集岩岩性参数的地层数据处理方法,利用计算储集岩孔隙流体体积来计算油藏石油地质储量,实现了储量计算的客观性和合理性。本发明消除了原有容积法储量公式本身产生的***误差,提高了探明石油地质储量计算精度。
附图简要说明
图1为饱和砂岩油藏地层脱气原油密度与地层油密度关系示意图;
图2为饱和砂岩油藏原始溶解气油比与与地层油密度关系示意图;
图3为单位岩石体积模型示意图;
图4为饱和砂岩油藏体积系数与原始溶解气油比关系示意图。
具体实施方式
依据本发明所述的一种基于储集岩岩性参数的地层数据处理方法,利用计算储集岩孔隙流体体积来计算油藏石油地质储量,具有客观、合理、无***误差的优点。该地层岩性、孔隙流体信息采集处理方法和地层测井信息主要基于以下事实进行研究来获得的。
在相同油藏地质条件下,利用地层有效孔隙度和有效含油饱和度(或有效原始含油饱和度),或利用地层总孔隙度和原始含油饱和度(或含油饱和度)计算储集岩孔隙中的石油体积是恒定的。也就是说,针对某一油藏而言,无论采用地层有效孔隙度还是采用地层总孔隙度计算石油地质储量,只要选用的含油饱和度与地层孔隙度相匹配,计算的石油地质储量是恒定不变的。
本发明采用如下技术方案来计算油藏石油地质储量,即一种基于储集岩岩性参数的地层数据处理方法,利用计算储集岩孔隙流体体积来计算油藏石油地质储量,其包括以下步骤:
1)利用地质勘探开发设备,可以采集地层岩性、孔隙流体及测井信息;
2)基于采集到的地层岩性数据,通过岩性数据分析设备获得储集岩的地层总孔隙度、地层有效孔隙度、束缚水饱和度、有效含水饱和度、有效束缚水饱和度和含水饱和度;
3)基于采集到的孔隙流体和测井信息,结合岩性数据分析设备获得油层有效厚度和含油面积;
4)基于采集到的地层孔隙流体信息,通过流体分析设备获得原始原油体积系数;
5)根据在步骤2)中获得的储集岩地层总孔隙度、地层有效孔隙度、束缚水饱和度、有效含水饱和度、有效束缚水饱和度和含水饱和度,来获得储集岩孔隙中的石油体积百分比;
在纯油层中:
OR=φt(1-Swi) (1)
或OR=φ(1-Swie) (2)
或OR=φtSoi (3)
或OR=φSoie (4)
在含可动水油层中:
OR=φt(1-Sw) (5)
或OR=φ(1-Swe) (6)
或OR=φtSo (7)
或OR=φSoe (8)
式中OR-石油体积百分比(指岩样孔隙中的石油所占据的孔隙体积与岩样总体积的比值),f;φ-地层有效孔隙度(指岩样的有效孔隙体积与岩样总体积的比值),f;φt-地层总孔隙度(指岩样的总孔隙体积与岩样总体积的比值),f;Soi-在纯油层中对应地层总孔隙度φt的原始含油饱和度,Soi=1-Sw1,f;Soie-在纯油层中对应地层有效孔隙度φ的有效原始含油饱和度,Soie=1-Swie,f;So-在含可动水油层中对应地层总孔隙度φt的含油饱和度,So=1-Sw,f;Soe-在含可动水油层中对应地层有效孔隙度φ的有效含油饱和度,Soe=1-Swe,f。
6)根据上面得到的石油体积百分比OR(单位,f)和在步骤3)中获得的油层有效厚度h(单位m)和含油面积Ao(单位km2),来获得含油储集岩孔隙中的石油体积Vf(单位104m3);
在纯油层中:
Vf=100AohφtSoi (9)
或Vf=100AohφSoie (10)
在可动水油层中:
Vf=100AohφtSo (11)
或Vf=100AohφSoe (12)
7)根据上面得到的石油体积Vf(单位104m3))和在步骤4)中获得的原始原油体积系数Boi(无因次),来获得油藏石油地质储量N′(单位104m3);
在纯油层中:
N′=100AohφtSoi/Boi (13)
或N′=100AohφSoie/Boi (14)
在可动水油层中:
N′=100AohφtSo/Boi (15)
或N′=100AohφSoe/Boi (16)
为了更进一步详细说明本发明的技术方案,即一种基于储集岩岩性参数的地层数据处理方法,利用计算储集岩孔隙流体体积来计算油藏石油地质储量,下面详细描述如何得到本发明中采用的地层岩性、孔隙流体和测井信息相关参数。但是本领域技术人员应当清楚,下面对各种相应参数的解释和说明仅仅是示例性的;不应当将本发明限制到具体实施例,本领域技术人员也可以经过创造性劳动,通过其他途径来获得储集岩地层总孔隙度、地层有效孔隙度、束缚水饱和度、有效束缚水饱和度、有效含水饱和度、含水饱和度、含油面积、油层有效厚度和原始原油体积系数等。
1地层参数
1.1岩石孔隙度(或称为,地面孔隙度或岩层孔隙度或地层孔隙度)
岩石孔隙度分有效孔隙度和总孔隙度,有效孔隙度是指岩样的有效孔隙体积与岩样总体积的比值;总孔隙度是指岩样的总孔隙体积与岩样总体积的比值。有效孔隙度又包括地面有效孔隙度和地层有效孔隙度;总孔隙度又包括地面总孔隙度和地层总孔隙度。
(1)地层总孔隙度φt
对固结岩石而言,可直接采用岩心孔隙压缩性样品分析结果确定地层总孔隙度;还可根据1986年J.P.Martin等人提出的声波地层因素方程求取地层总孔隙度:
式中:φt-地层总孔隙度,f;Δt-测井声波时差,μs/m;Δtma-骨架时差,μs/m;x-岩性系数。
根据固结岩石有限的岩心分析地层总孔隙度φt和对应的实测声波时差Δt,由(17)式变形后回归求得骨架时差Δtma和岩性系数x。
logΔt=logΔtma-xlog(1-φt) (18)
再根据声波时差Δt,由(17)式计算固结岩石的地层总孔隙度φt。
对疏松岩石而言,实验室分析的地面总孔隙度φa需要校正为地层总孔隙度,分岩性后的孔隙度校正方程分别为:
细砾状砂岩:
φt=φa[0.381+0.619(Pe+1)-0 1628] (19)
不等粒砂岩:
φt=φa[0.381+0.619(Pe+1)-0 1169] (20)
中-粗砂岩:
φt=φa[0.381+0.619(Pe+1)-0 0907] (21)
中-细砂岩:
φt=φa[0.381+0.619(Pe+1)-0 0675] (22)
细砂岩:
φt=φa[0.381+0.619(Pe+1)-0 0481] (23)
粉-细砂岩:
φt=φa[0.381+0.619(Pe+1)-0 0298] (24)
净上覆岩层压力Pe(单位MPa)可根据样品埋深D,由下式求得:
Pe=0.01Dρb-Pf (25)
式中ρb-上覆岩石平均密度,取2.30g/cm3;Pf-实测地层孔隙流体压力,MPa。
疏松岩石采用(19)~(25)式即可求得地层总孔隙度φt。
(2)地层有效孔隙度φ
在实验室采用煤油法分析的孔隙度为地面有效孔隙度。对固结岩石而言,地层有效孔隙度可采用孔隙压缩性岩心分析结果直接计算石油地质储量;对疏松岩石而言,地面有效孔隙度需经孔隙压缩性研究后由下列方程形式回归求取地层有效孔隙度φ:
φ=A+Bφe (26)
式中:A、B-回归系数;φe-实验室分析地面有效孔隙度,f。
1.2含油饱和度
含油饱和度分有效含油饱和度和原始含油饱和度,原始含油饱和度是指岩样中石油所占据的孔隙体积与岩样总孔隙体积的比值;有效含油饱和度是指岩样中石油所占据的孔隙体积与岩样有效孔隙体积的比值。
(1)原始含油饱和度
原始含油饱和度与束缚水饱和度的关系:
Soi=1-Swi (27)
原始含油饱和度利用压汞资料求取。首先利用J(Sw)函数对毛管压力曲线进行分类和平均:
式中:σ,θ和Pc-分别为实验室内的界面张力、接触角和毛管压力;K-砂岩渗透率,mD;φ-地层有效孔隙度,f。
J(Sw)函数确定原始含油饱和度的具体步骤:
第一步,计算每块样品的C值。
对某一样品的毛管压力曲线而言,σ,cosθ,K和φ均为常数。为计算方便,令(28)式中的常数项为C,即:
第二步,计算岩样每个测压点的J(Sw)函数值。
J(Sw)=CPc (30)
第三步,统计不同饱和度区间的平均汞饱和度、平均J(Sw)函数和平均压力数据。
第四步,以纵坐标表示J(Sw)函数,横坐标表示汞饱和度SHg,点上J(Sw)点子。
如果点子集中,说明这些样品同属于一种孔隙结构类型,可以综合为一条代表该类储层的J(Sw)函数曲线。
第五步,求取该类储层的平均毛管压力曲线。
已知该类储层的平均渗透率和平均孔隙度,利用(29)式算出平均而后将J(Sw)函数曲线上任一含汞饱和度对应的J(Sw)值乘以平均值的倒数,即可得到该点对应的平均毛管压力值进而得出该类储层的平均毛管压力曲线。
第六步,将实验室条件下的平均毛管压力曲线换算为油藏条件下的毛管压力曲线。
在实验室内测量的毛管压力曲线,虽然使用的样品是实际岩心,但非润湿相和润湿相流体不能直接采用油藏条件下的流体,而不同组合的流体界面张力和润湿接触角不同,测得的毛管压力也各不相同。因此无论采用油-水或水银-空气那个***测得的毛管压力曲线,在计算原始含油饱和度时,必须校正为油藏条件下的平均毛管压力曲线。
实验室毛管压力表达式和油藏毛管压力表达式分别为:
式中:σ,θ和Pc-分别为实验室内的界面张力、接触角和毛管压力;σR,θR和PcR-分别为油藏条件下的界面张力、接触角和毛管压力。
由(32)、(33)式联立得:
将实验室分析的参数(表1为界面张力、接触角实验室分析成果)代入(34)式得:
表1界面张力、接触角实验室分析成果
第七步,将油藏条件下的毛管压力换算为油柱高度。
油藏的毛管压力由油水的重力差来平衡,而毛管压力可表示为:
PcR=10-3H(ρwf-ρof)g (36)
式中:H-油藏的自由水面以上高度,m;g-重力加速度,9.80m/s2。
将(36)式简化后得:
式中:PcR-油藏毛管压力,MPa;ρwf、ρof-分别为油藏条件下地层水密度和地层油密度,g/cm3。
第八步,求取原始含油饱和度。
将(35)式代入(37)式得:
由(38)式和油藏平均毛管压力曲线即可得到原始含油饱和度与油藏高度关系,进而求得不同油藏高度下的原始含油饱和度和束缚水饱和度。上式中的地层油密度ρof和地层水密度ρwf采用下列方法求取。
1)地层油密度计算方法
地层油密度是指在原始地层条件下,单位体积石油的质量(g/cm3)。
根据辽河油区已录取的数据,将地面脱气原油密度转换为地层温度下的地层脱气原油密度,从而使地层油密度与地层脱气原油密度具有统一的对比基础,提高了统计分析的相关性。分别建立了不同油藏地质条件下地层油密度统计方程。
饱和砂岩油藏地层油密度统计方程(图1为饱和砂岩油藏地层脱气原油密度与地层油密度关系示意图、图2为饱和砂岩油藏原始溶解气油比与与地层油密度关系示意图)
ρof=1.7527ρot-0.7298,n=280,r=0.9380 (39)
ρof=1.3486-0.3116logRsi,n=280,r=0.9962 (40)
前第三系及高凝油不饱和油藏地层油密度ρof统计方程:
ρof=0.7997ρot+0.1855,n=44,r=0.90896 (41)
其中(39)式和(41)式中的地层脱气原油密度ρot(指地层温度t)是通过地面分析20℃时脱气原油密度ρo求取:
ρot=ρo+(13.560-0.191ξ)×10-3-(63.900-0.870ξ)×10-5t102 (42)
ξ=INT[100(ρo-0.8001)] (43)
2)地层水密度计算方法
采用S-K系列方程计算地层水密度ρwf:
当S·C≤12964.28mg/L时:
ρwf=1.00249ρ′wf+1.3094×10-6S·C-4.7710×10-7ρ′wfS·C-0.00236 (45)
当S·C>12964.28mg/L时:
ρwf=1.02827ρ′wf+5.3106×10-8S·C+6.3135×10-7ρ′wfS·C-0.02376 (46)
当Cl-1≤8000mg/L时:
ρwf=1.00249ρ′wf+2.1219×10-6Cl-1-7.7316×10-7ρ′wfCl--0.00236 (47)
当Cl-1>8000mg/L时:
ρwf=1.64215ρ′wf-8.6060×10-8Cl-1-1.0231×10-6ρ′wfCl-+0.02788 (48)
式中的地层水矿化度S·C(单位mg/L)为实测值;ρ′wf为纯水密度(单位g/cm3);Cl-1为氯离子浓度(单位mg/L);地层温度t(单位℃)是根据研究地区的实测地层温度与垂直深度D(单位m)建立的统计方程计算;地层压力Pi(单位MPa)是根据实测地层压力与垂直深度D建立的统计方程计算。
Pi=a+bD (50)
式中a、b是回归系数。
(2)有效含油饱和度
有效含油饱和度与有效含水饱和度关系式为:
Soe=1-Swe (51)
根据原始含油饱和度与有效含油饱和度的关系式求取有效含油饱和度。
储集岩含油体积是指具有工业性开采价值的含油储集岩总体积,为含油面积与油层有效厚度的乘积;而石油体积是指在地层条件下含油储集岩中石油所占据的孔隙体积。石油体积与原始原油体积系数的比值即为地面脱气状态下的石油地质储量。
由单位岩石体积模型(图3为单位岩石体积模型示意图)可推导出地层孔隙中石油体积与地层总孔隙度、地层有效孔隙度及相应饱和度间的关系方程,推导过程如下:
Vma=Vsma+Vshma (52)
Vp=Vshwi+Vswi+Vsb+Vsf (53)
V=Vma+Vp=1 (54)
式中V-单位岩石体积,V=1m3;Vma-单位岩石骨架体积,m3;Vp-单位岩石孔隙体积,m3;Vshma-单位岩石泥质骨架体积,m3;Vsma-单位岩石砂岩骨架体积,m3;Vshwi-单位岩石泥质或粘土中的流体体积,m3;Vswi-单位岩石泥质以外无效孔隙中的束缚水体积,m3;Vsb-单位岩石有效孔喉中薄膜水体积,m3;Vsf-单位岩石有效孔喉中可动流体体积,m3;Vf-地层条件下的石油体积,m3。
根据单位岩石体积模型,结合孔隙度定义,建立地层总孔隙度φt和地层有效孔隙度φ理论表达式:
对纯油层而言,与地层总孔隙度φt对应的束缚水饱和度Swi理论表达式,以及与原始含油饱和度Soi的理论关系式为:
Soi=1-Swi (58)
根据地层总孔隙度和原始含油饱和度可导出对应纯油层石油体积百分比OR(单位f):
对纯油层而言,与地层有效孔隙度对应的有效束缚水饱和度Swie理论表达式,以及与有效原始含油饱和度Soie的理论关系式为:
Soie=1-Swie (61)
根据地层有效孔隙度φ和有效原始含油饱和度Soie可导出对应纯油层石油体积百分比:
由(58)、(59)、(61)、(62)式联立得纯油层地层总孔隙度φt与地层有效孔隙度φ理论关系式:
φt(1-Swi)=φ(1-Swie) (63)
或φtSoi=φSoie (64)
对含可动水油层而言,与地层总孔隙度φt对应的含水饱和度Sw表达式,以及与含油饱和度So和可动水饱和度Swf等饱和度参数的理论关系式为:
So=1-Sw (66)
Sw=Swi+Swf (67)
So=Soi-Swf (68)
式中的含油饱和度So可通过实验室岩心分析直接获取。
根据地层总孔隙度φt和含油饱和度So可导出对应含可动水油层石油体积百分比:
对含可动水油层而言,与地层有效孔隙度φ对应的有效含水饱和度Swe定义,以及与有效含油饱和度Soe和有效可动水饱和度Swfe等饱和度参数的理论关系式为:
Swe+Soe=1 (71)
Swe=Swie+Swfe (72)
Soe=Soie-Swfe (73)
根据地层有效孔隙度和有效含油饱和度可导出对应含可动水油层石油体积百分比:
由(66)、(69)、(71)、(74)式联立得含可动水油层地层总孔隙度φt与地层有效孔隙度φ理论关系式:
φt(1-Sw)=φ(1-Swe) (75)
或φtSo=φSoe (76)
1.3油层有效厚度
(1)有效厚度标准确定
分别制定油层划分和夹层扣除标准;以岩心分析资料和测井解释资料为基础,测试资料为依据,在研究岩性、物性、电性与含油性关系后,确定其有效厚度划分的岩性、物性、电性下限标准;储层特征和流体性质相近的多个小型油藏,可分别制定统一的标准;借用邻近油藏下限标准应论证类比依据和标明参考文献;有效厚度标准图版符合率大于80%。
(2)有效厚度划分
以测井解释资料划分有效厚度时,应对有关测井曲线进行必要的井筒环境(如井径变化、泥浆侵入等)校正和不同测井系列的归一化处理;以岩心分析资料划分有效厚度时,油气层段应取全岩心,收获率不低于80%;有效厚度的起算厚度为0.2~0.4m,夹层起扣厚度为0.2m。
(3)计算单元有效厚度取值
计算储量需要取得代表储量计算单元的油层有效厚度。计算单元油层有效厚度取值方法与油田地质条件和井点分布情况有关。归纳起来有两种方法:算术平均和面积权衡。
1)算术平均值
井点分布较均匀的勘探开发区,储量计算单元油层有效厚度采用算术平均值:
式中计算单元有效厚度平均值,m;hk-第k口井的油层有效厚度,m;n-计算单元井数,口。
2)面积权衡值
井点分布不均匀的评价钻探区,油层有效厚度采用面积权衡值,分以下两种形式:
①等值线面积权衡
将井与井之间的油层厚度视为线性变化,采用直线内插法绘制油层有效厚度等值图。
式中Ai-第i条与第i+1条等值线间的含油面积,km2;hi、hi+1-第i条与第i+1条等值线的油层有效厚度,m。
②井点面积权衡
这一方法的基础是:单井控制面积为该井至邻井距离二分之一范围内的面积,各井所能控制的面积大小随井距变化而异,每口井所钻遇的油层有效厚度代表该井控制面积内的油层有效厚度。
具体做法:首先将邻近的井点依次连成三角网,然后取锐角三角形的中垂线(钝角三角形取中点连线)划分的单井控制含油面积。油水过渡带厚度取邻井厚度之半,岩性边界至油井距离二分之一处取相邻油井厚度之半,纯油区和油层厚度稳定区域平均油层有效厚度取井点面积权衡值,计算公式为:
式中Ak-第k口井控制的含油面积,km2。
1.4含油面积
已开发探明储量的含油面积,根据生产井静态和动态资料综合圈定。
未开发探明储量的含油面积,各种边界的确定需达到以下条件:
1)用以圈定含油面积的流体界面,应经测井或测试资料,或钻井取心资料证实,或可靠的压力测试资料确定。
2)未查明流体界面的油藏,应以测试证实的最低的出油层(或井段)底界,或有效厚度累计值或集中段高度外推圈定含油面积。
3)油藏断层(或地层)遮挡边界,宜以油层顶(底)面与断层(或地层不整合)面相交的外含油边界圈定含油面积。
4)油藏储层岩性(或物性)遮挡边界,用有效厚度零线或渗透储层厚度线圈定含油面积;未查明边界时以开发井距的1~1.5倍外推划计算线。
5)在储层厚度和埋藏深度等条件下,高分辩率地震解释预测的流体界面和岩性边界,经钻井资料约束解释时,可作为圈定含油面积的依据。
6)在确定的含油边界内,边部油井到含油边界的距离大于3倍开发井距时,可按照油藏开发井距的1~1.5倍外推划计算线。
1.5原始原油体积系数
原始原油体积系数Boi是指在原始地层条件下的石油体积Vf(单位m3)与地面标准条件下脱气原油体积Vs(单位m3)之比:
确定计算单元原油体积系数有三种方法:
一是取样分析法,这是最准确、最有效的方法,但对多数计算单元而言都难以实现。
二是类比法,通过与已知油藏的地质条件和原油性质等相关参数类比分析,直接选择类似油藏的原始原油体积系数,这种方法更适合于欠饱和油藏及特殊类型油藏,当然饱和油藏也适用。
三是经验方程计算法,这种方法不适合特殊类型油藏。
辽河油区依据已录取的数百口井的原油高压物性分析数据,分岩性、分区带建立了原油体积系数统计方程。
①辽河油区饱和砂岩油藏原油体积系数Boi统计方程:
②辽河油区前第三系不饱和油藏原油体积系数Boi统计方程:
③西部凹陷潜山(Pt、Ar)油藏原油体积系数Boi统计方程:
④大民屯凹陷潜山(Pt、Ar)油藏原油体积系数Boi统计方程:
⑤辽河油区饱和砂岩油藏体积系数Boi与原始溶解气油比Rsi统计方程(图4为饱和砂岩油藏体积系数与原始溶解气油比关系示意图):
Rsi=321(Boi-1)1 02,n=307,r=0.9868 (85)
⑥辽河油区前第三系不饱和油藏体积系数Boi与原始溶解气油比Rsi统计方程:
Rsi=280(Boi-1)1 02,n=43,r=0.9524 (86)
2石油体积及石油地质储量计算
2.1储集岩孔隙中的石油体积Vf(单位104m3)
纯油层储集岩孔隙中的石油体积:
Vf=100AohφtSoi
或Vf=100AohφSoie
含可动水油层储集岩孔隙中的石油体积:
Vf=100AohφtSo
或Vf=100AohφSoe
2.2地面脱气石油地质储量N′
纯油层石油地质储量计算方程:
N′=100AohφtSoi/Boi
或N′=100AohφSoie/Boi
含可动水油层石油地质储量计算方程:
N′=100AohφtSo/Boi
或N′=100AohφSoe/Boi
本发明技术方案中所述方法,经过实际资料验证,取得了明显的效果,具体为:上述系列容积法储量计算方程在辽河油田十个储量计算单元应用后,效果显著,使探明地质储量计算精度在原来基础上提高了6%~21%。下面列举实例加以说明。
实例
在辽河油区优选了10个不同孔隙度、饱和度区间的探明储量计算单元,并根据实际分析的孔隙度、饱和度数据,采用本发明和原来的容积法储量公式计算探明储量,其相对误差为-6.0%~-21.3%,平均为-11.0%(表2),属***误差。
表2不同孔隙度和饱和度区间计算的地质储量误差分析
如上所述,已经清楚详细地描述了本发明提出的一种基于储集岩岩性参数的地层数据处理方法来计算油藏石油地质储量。尽管本发明的优选实施例详细描述并解释了本发明,但是本领域普通的技术人员可以理解,在不背离所附权利要求定义的本发明的精神和范围的情况下,可以在形式和细节中做出多种修改。
Claims (1)
1.一种基于储集岩岩性参数的地层数据处理方法,包括以下步骤:
1)利用地质勘探开发设备,可以采集地层岩性、孔隙流体及测井信息;
2)基于采集到的地层岩性数据,通过岩性数据分析设备获得储集岩的地层总孔隙度φt、地层有效孔隙度φ、束缚水饱和度Swi、有效含水饱和度Swe、有效束缚水饱和度Swie和含水饱和度Sw;
3)基于采集到的孔隙流体和测井信息,结合岩性数据分析设备获得油层有效厚度h和含油面积Ao,油层有效厚度h的单位为m,含油面积Ao的单位为km2;
4)基于采集到的地层孔隙流体信息,通过流体数据分析设备获得原始原油体积系数Boi;
5)根据在步骤2)中获得的储集岩地层总孔隙度φt、地层有效孔隙度φ、束缚水饱和度Swi、有效含水饱和度Swe、有效束缚水饱和度Swie和含水饱和度Sw,来获得储集岩孔隙中的石油体积百分比OR;
其中,OR为石油体积百分比,指岩样孔隙中的石油所占据的孔隙体积与岩样总体积的比值,单位f;
6)根据得到的石油体积百分比OR和在步骤3)中获得的油层有效厚度h和含油面积Ao,来获得含油储集岩孔隙中的石油体积Vf,单位104m3;
7)根据得到的石油体积Vf和在步骤4)中获得的原始原油体积系数Boi,来获得油藏石油地质储量N′,单位104m3,石油地质储量N′分别为下面几种情形;
对于纯油层:N′=100AohφtSoi/Boi;或N′=100AohφSoie/Boi;
对于含可动水油层:N′=100AohφtSo/Boi;或N′=100AohφSoe/Boi;
其中,石油体积百分比OR分别为下列情况:
纯油层为:OR=φt(1-Swi);或OR=φ(1-Swie);或OR=φtSoi;或OR=φSoie;
含可动水油层为:OR=φt(1-Sw);或OR=φ(1-Swe);或OR=φtSo;或OR=φSoe;
其中,φ为地层有效孔隙度,是指岩样的有效孔隙体积与岩样体积的比值,单位f;φt为地层总孔隙度,是指岩样的总孔隙体积与岩样体积的比值,单位f;Soi为在纯油层中对应地层总孔隙度φt的原始含油饱和度,在纯油层中原始含油饱和度是指岩样中石油所占据的孔隙体积与岩样总孔隙体积的比值,并有Soi=1-Swi,单位f;Soie为在纯油层中对应地层有效孔隙度φ的有效原始含油饱和度,在纯油层中有效原始含油饱和度是指岩样中石油所占据的孔隙体积与岩样有效孔隙体积的比值,并有Soie=1-Swie,单位f;So为在含可动水油层中对应地层总孔隙度φt的含油饱和度,在含可动水油层中含油饱和度是指岩样中石油所占据的孔隙体积与岩样总孔隙体积的比值,并有So=1-Sw,单位f;Soe为在含可动水油层中对应地层有效孔隙度φ的有效含油饱和度,在含可动水油层中有效含油饱和度是指岩样中石油所占据的孔隙体积与岩样有效孔隙体积的比值,并有Soe=1-Swe,单位f;
储集岩孔隙中的石油体积Vf分别为下列情况:纯油层时为Vf=100AohφtSoi;或Vf=100AohφSoie;含可动水油层时为Vf=100AohφtSo;或Vf=100AohφSoe。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100504160A CN102174888B (zh) | 2011-03-03 | 2011-03-03 | 一种基于储集岩岩性参数的地层数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100504160A CN102174888B (zh) | 2011-03-03 | 2011-03-03 | 一种基于储集岩岩性参数的地层数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102174888A CN102174888A (zh) | 2011-09-07 |
CN102174888B true CN102174888B (zh) | 2012-05-23 |
Family
ID=44518108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100504160A Active CN102174888B (zh) | 2011-03-03 | 2011-03-03 | 一种基于储集岩岩性参数的地层数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102174888B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989233A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种地层数据处理方法 |
US11847391B2 (en) | 2020-06-29 | 2023-12-19 | Dassault Systemes Simulia Corp. | Computer system for simulating physical processes using surface algorithm |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631907A (zh) * | 2014-10-29 | 2016-06-01 | 中国石油天然气股份有限公司 | 油气藏龟背图绘制方法及装置 |
CN104933246A (zh) * | 2015-06-15 | 2015-09-23 | 中国海洋石油总公司 | 一种定量表征气田储量参数对储量变化贡献的方法 |
CN105221144B (zh) * | 2015-09-21 | 2018-04-06 | 中国石油天然气股份有限公司 | 确定油藏储量的方法及装置 |
CN106772664A (zh) * | 2015-11-24 | 2017-05-31 | 中国石油化工股份有限公司 | 海上新区滚动勘探开发的方法 |
CN106845667A (zh) * | 2016-11-23 | 2017-06-13 | 红有软件股份有限公司 | 石油、天然气储量全生命周期管理*** |
CN110573842B (zh) | 2017-01-26 | 2021-05-11 | 达索***西姆利亚公司 | 基于流体占据时间的多相流可视化图 |
CN107120106B (zh) * | 2017-05-25 | 2020-06-23 | 中石化石油工程技术服务有限公司 | 基于有机孔隙度和总有机碳含量的页岩品质评价方法 |
US11714040B2 (en) * | 2018-01-10 | 2023-08-01 | Dassault Systemes Simulia Corp. | Determining fluid flow characteristics of porous mediums |
CN110032744B (zh) * | 2018-01-11 | 2021-11-09 | 中国石油化工股份有限公司 | 一种非常规油气储量计算方法及*** |
US11530598B2 (en) | 2018-08-21 | 2022-12-20 | Dassault Systemes Simulia Corp. | Determination of oil removed by gas via miscible displacement in reservoir rock |
CN109343145B (zh) * | 2018-11-28 | 2020-04-14 | 程立华 | 基于修正储量丰度的低渗致密砂岩气藏富集区确定方法 |
CN111521532B (zh) * | 2019-02-01 | 2023-02-10 | 中国石油天然气股份有限公司 | 一种致密储层含气饱和度的测定方法 |
CN110130884B (zh) * | 2019-04-28 | 2022-05-10 | 中国石油天然气股份有限公司 | 一种气井剩余潜力确定方法 |
CN112558153B (zh) * | 2019-09-25 | 2022-03-29 | 中国石油天然气股份有限公司 | 一种双相介质的油气储层预测方法及装置 |
CN110929383A (zh) * | 2019-10-28 | 2020-03-27 | 中国石油大港油田勘探开发研究院 | 一种源下油气柱高度的计算方法 |
CN111127644B (zh) * | 2019-12-12 | 2021-03-02 | 北京龙软科技股份有限公司 | 一种露天煤矿绘制煤层厚度等值线的方法和装置 |
CN114428166B (zh) * | 2020-09-29 | 2024-06-04 | 中国石油化工股份有限公司 | 一种地层含油气性的评价方法及*** |
US11907625B2 (en) | 2020-12-29 | 2024-02-20 | Dassault Systemes Americas Corp. | Computer simulation of multi-phase and multi-component fluid flows including physics of under-resolved porous structures |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4841490A (en) * | 1987-06-19 | 1989-06-20 | Schlumberger Technology Corporation | Inversion of a vertical seismic profile by minimization of an entropy like function |
CN1164651A (zh) * | 1996-10-29 | 1997-11-12 | 中国石油天然气总公司西北地质研究所 | 地震-测井多信息储层参数反演方法 |
CN101078776A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 核测井数据多尺度处理方法 |
CN101787884A (zh) * | 2010-01-28 | 2010-07-28 | 中国石油集团川庆钻探工程有限公司 | 声波孔隙度和中子孔隙度差值储层流体类型判别方法 |
-
2011
- 2011-03-03 CN CN2011100504160A patent/CN102174888B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4841490A (en) * | 1987-06-19 | 1989-06-20 | Schlumberger Technology Corporation | Inversion of a vertical seismic profile by minimization of an entropy like function |
CN1164651A (zh) * | 1996-10-29 | 1997-11-12 | 中国石油天然气总公司西北地质研究所 | 地震-测井多信息储层参数反演方法 |
CN101078776A (zh) * | 2007-06-29 | 2007-11-28 | 西安奥华电子仪器有限责任公司 | 核测井数据多尺度处理方法 |
CN101787884A (zh) * | 2010-01-28 | 2010-07-28 | 中国石油集团川庆钻探工程有限公司 | 声波孔隙度和中子孔隙度差值储层流体类型判别方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989233A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种地层数据处理方法 |
US11847391B2 (en) | 2020-06-29 | 2023-12-19 | Dassault Systemes Simulia Corp. | Computer system for simulating physical processes using surface algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN102174888A (zh) | 2011-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102174888B (zh) | 一种基于储集岩岩性参数的地层数据处理方法 | |
CN101487390B (zh) | 一种确定油层原始含油饱和度的阿尔奇模式方法 | |
CN105134195A (zh) | 一种基于测井资料的页岩气储层品质评价方法 | |
CN101775983B (zh) | 基于砂岩储层水层电阻率的地层数据处理方法 | |
CN102707333A (zh) | 页岩气资源量/储量的测量方法 | |
Harris | The role of geology in reservoir simulation studies | |
Alpak et al. | Simplified modeling of turbidite channel reservoirs | |
CN104834003B (zh) | 一种非常规致密储层的相控压缩系数地震预测方法 | |
CN102621586B (zh) | 一种识别地层属性的地层数据处理方法 | |
CN110424956A (zh) | 页岩油资源量计算中评价单元保存系数权重量化赋值方法 | |
Anderson | Coring and core analysis handbook | |
Deng et al. | An integrated workflow for reservoir modeling and flow simulation of the nikanassin tight gas reservoir in the western Canada sedimentary basin | |
Yang et al. | Addressing microseismic uncertainty from geological aspects to improve accuracy of estimating stimulated reservoir volumes | |
Osborne | Neural networks provide more accurate reservoir permeability | |
Dubois et al. | Multiscale geologic and petrophysical modeling of the giant Hugoton gas field (Permian), Kansas and Oklahoma, USA | |
CN102134994B (zh) | 一种基于基岩油藏油水层电阻率的地层数据处理方法 | |
Eikeland et al. | Dry gas reinjection in a strong waterdrive gas/condensate field increases condensate recovery—case study: The Sleipner Øst Ty field, South Viking Graben, Norwegian North Sea | |
Cherian et al. | Understanding development drivers in horizontal wellbores in the Midland Basin | |
Osborne | Permeability estimation using a neural network: A case study from the Roberts unit, Wasson field, Yoakum County, Texas | |
Worthington et al. | Optimizing the value of reservoir simulation through quality-assured initialization | |
Zhang et al. | Multi-parameters logging identifying method for sand body architectures of tight sandstones: A case from the Triassic Chang 9 Member, Longdong area, Ordos Basin, NW China | |
CN114059999A (zh) | 一种重力流沉积成因测井识别方法 | |
Almasoodi et al. | Stochastic-based Coupling of Static and Dynamic Models: An Example From the Meremac Formation in the STACK Play | |
Wang et al. | 4D Geomechanical Research in Hydraulic Fracturing of Unconventional Resources-A Case Study From Junggar Basin, China | |
Famubode et al. | Integrating Geochemistry with Other Reservoir Diagnostic Tools–Implications for Infill Planning and Development in the Midland Basin |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C56 | Change in the name or address of the patentee |
Owner name: CHINA PETROLEUM +. NATURAL GAS CO., LTD. Free format text: FORMER NAME: KANG ZHIYONG |
|
CP03 | Change of name, title or address |
Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: PetroChina Company Limited Address before: 124010 development and synthesis Institute, exploration and Development Institute, No. 95 petroleum street, Xinglongtai District, Liaoning, Panjin Patentee before: Kang Zhiyong |