CN104239703B - 页岩气储层多参数定量类比评价方法 - Google Patents
页岩气储层多参数定量类比评价方法 Download PDFInfo
- Publication number
- CN104239703B CN104239703B CN201410446356.8A CN201410446356A CN104239703B CN 104239703 B CN104239703 B CN 104239703B CN 201410446356 A CN201410446356 A CN 201410446356A CN 104239703 B CN104239703 B CN 104239703B
- Authority
- CN
- China
- Prior art keywords
- gas
- evaluation
- interval
- explained
- well
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及一种页岩气储层多参数定量类比评价方法，通过录井、测井资料获取待解释的目标层段孔隙度、有机碳含量、甲烷含量、总含气量、含气饱和度、脆性指数、厚度等页岩气储层评价关键参数与具有区域代表性的参考层的关键评价参数，绘制多参数类比图版，用目测法观察和对比，以直观定性评价待解释层的含气性；根据类比评价值＝目标层参数值/参考层参数值计算各参数的类比评价值，计算多参数定量类比评价指数Ia，绘制多参数定量类比相关性判别图版和计算多参数类比相关强度R2，绘制Ia—R2评价交会图，将待解释目标层的类比评价指数Ia—类比相关强度R2数据点置于类比评价指数Ia—类比相关强度R2交会图版中；输出结果。
Description
技术领域
本发明涉及一种基于录井与测井资料的页岩气储层多参数定量类比评价方法，用于页岩气勘探开发录井与测井资料解释评价。
背景技术
页岩气是一种新型的清洁能源，在国内勘探开发尚处于起步阶段，页岩气田发现较少，应用录井、测井资料评价页岩气储层的经验与方法相对缺乏，更缺少实用的页岩气测录井定量解释评价方法。
国内现有的页岩气储层录井、测井解释评价主要是以传统的常规油气层解释方法为主，引入了有机碳含量TOC、游离气含量Gf、吸附气含量Gs、总含气量Gt、脆性指数BRIT等页岩气储层评价参数，但在解释过程中，主要是用表格化的数据笼统定性对比评价，不能形成具有综合性的关键评价指标，解释方法相对单一，直观性不强，定量性差，具有多解性，有待通过实践研究，形成多参数定量类比评价方法，特别是图形化的评价方法，提高解释成果的可靠性，减少多解性。
页岩气勘探作为非常规油气资源勘探中的一个重要组成部分，已成为我国及全世界油气勘探的热门领域。准确、可靠的页岩气储层识别与评价结果能够为决策部门提供重要的决策依据。因此，在页岩气录井与测井资料解释评价方面，国内现场急需具有创新性的先进、实用技术。
发明内容
本发明的目的是针对上述技术现状，旨在提供一种可以直观展示录井、测井不同类别、不同计量单位的解释关键参数间的综合量化关系，适用于页岩气储层录井、测井解释评价的页岩气储层多参数定量类比评价方法。
本发明目的的实现方式为，页岩气储层多参数定量类比评价方法，具体步骤如下：
第一步：建立解释数据库
1)建立已知井对比解释层段关键参数数据库，
选择已试气并具有区域代表性的典型页岩气层作为对比层
选择区域经试气证实的页岩气层作为已知井对比层；选取对比层录井、测井解释关键参数孔隙度POR(％)，有机碳含量TOC(％)，全烃或甲烷含量QT(％)，总含气量Gt(％)，含气饱和度Sg(％)，脆性指数BRIT(％)，厚度H(m)、镜质体反射率Ro(％)、钻时ROP或钻时比值RROP、烃对比系数Kc、地层压力FP或地层压力梯度FPG中的3至10个，作为定量类比评价关键参数；
2)建立待解释井对比层解释关键参数数据库，
以待解释层段为对比对象，选择试气已证实为气层的井，在对比层段内按定步长读取录井、测井解释关键参数，以文本文件或WIS文件格式建立数据库；根据气层产量多选几个气层段对比，最好是选择一个达到接近工业气流产能的气层段；
第二步：读取解释关键参数
1)优选解释关键参数，
根据资料录取及地区地质特点优选解释关键参数；
2)划分待解释井页岩气异常显示层段，
依据录井地层岩性、甲烷、钻时、有机碳含量和测井自然伽马、密度曲线变化判断页岩气显示储层；以不扣除垂厚小于2m薄夹层为原则，划分页岩气储层评价井段，并作为待解释层段；
3)读取待解释层段录井与测井关键参数，
一个解释层段读取一组典型数据，参数分别为孔隙度POR2、有机碳含量TOC2、全烃或甲烷含量QT2、总含气量Gt2、含气饱和度Sg2、脆性指数BRIT2、储层厚度H2及钻时ROP2或钻时比值RROP2、地层压力梯度FPG2；
4)读取已知井对比层段录井与测井关键参数，
对应待解释层段取值方式，读取参数孔隙度POR1、有机碳含量TOC1、全烃或甲烷含量QT1、总含气量Gt1、含气饱和度Sg1、脆性指数BRIT1、储层厚度H1及钻时ROP1或钻时比值RROP1、地层压力梯度FPG1；
第三步：处理解释评价参数
1)计算类比评价值
类比评价值具体公式如下：
IPOR＝POR2/POR1，
ITOC＝TOC2/TOC1，
IQT＝QT2/QT1，
IGt＝Gt2/Gt1，
ISg＝Sg2/Sg1，
IBRIT＝BRIT2/BRIT1，
IH＝H2/H1,
IROP＝ROP1/ROP2,或IROP＝RROP2/RROP1
IFPG＝FPG2/FPG1。
式中，IPOR、ITOC、IQT、IGt、ISg、IBRIT、IH、IROP、IFPG分别为孔隙度、总有机碳含量、烃含量、总含气量、含气饱和度、脆性指数、储层厚度、钻时、地层压力梯度类比评价值，均用小数表示，无量纲；
2)计算类比评价指数Ia，采用几何平均数算法计算。
第四步：绘制解释评价图版
1)绘制多参数类比图版，
以孔隙度、有机碳含量、甲烷、总含气量、含气饱和度、脆性指数、厚度等不同类别的评价参数为横轴，评价参数数值为纵轴，绘制多参数类比图版；
分别将读取到的对比层、待解释层数值投射到多参数类比图上，用目测法观察和对比其相似性及差异性，直观定性评价待解释层段的含气性；
2)绘制多参数相关性判别图版及计算相关强度R2，
以对比层段参数值为横轴，待解释层段参数值为纵轴，绘制多参数相关性判别图版，并做线性回归分析，得出线性相关方程及相关强度R2；
3)绘制类比评价指数Ia—相关强度R2评价交会图，
以相关强度R2为横轴，类比评价指数Ia为纵轴，绘制Ia—R2交会图；
将区域内经过试气证实的已知井气层、含气层数据点绘到Ia—R2交会图上，根据数据统计原理确定解释评价标准，即确定不同储层含气类别界线AB与CD，界线AB与CD将图版分为页岩气层Ⅰ、含气层Ⅱ两个区域，用于表征页岩气储层评价结果；
将待解释层段的多参数类比评价指数Ia、相关强度R2数据点置于Ia—R2交会图版中；
第五步：输出评价结果
待解释层段数据交会点所处区域即为对应评价结果：Ⅰ区为气层，Ⅱ区为含气层。
本发明的实质是通过录井、测井资料获取待解释的目标层段孔隙度POR、有机碳含量TOC、全烃或甲烷含量QT、总含气量Gt、含气饱和度Sg、脆性指数BRIT、厚度H、镜质体反射率Ro、钻时比值ROPn/s、烃对比系数Kc等页岩气储层评价关键参数；一般选取3至10个，与经试气证实的并具有区域代表性的典型页岩气层进行多参数类比；通过计算多参数定量类比指数Ia和相关强度R2，建立类比评价指数Ia—类比相关强度R2交会图，用图形化直观显示、定量数据表征多参数类比评价结果。
本发明在四川盆地及周缘地区30多口井进行了生产应用，评价的页岩气层试气均获得了良好工业气流，符合率达到90％以上，解释精度得到大幅度提高，成为一种有效的页岩气储层录井、测井多参数定量评价方法，为页岩等非常规储层解释评价提供了崭新的手段和方法。
附图说明
图1为本发明工作流程框图，
图2为本发明工作原理图，
图3为多参数类比图版式样图，
图4为多参数定量类比相关性判别图式样图，
图5为定量类比指数Ia—相关强度R2评价交会图式样图，
图6为本发明JYHF-1井多参数类比图版，
图7为本发明JYHF-1井多参数定量类比相关性判别图版，
图8为本发明JYHF-1井定量类比指数Ia—相关强度R2评价交会图，
图9为本发明JY2井多参数类比图版，
图10为本发明JY2井多参数定量类比相关性判别图版，
图11为本发明JY2井定量类比指数Ia—相关强度R2评价交会图。
具体实施方式
参照图1、图2、图3、图4、图5，本发明的具体步骤(见图1、图2)为：
第一步：建立解释数据库
1)建立已知井对比解释层段关键参数数据库，
选择已试气并具有区域代表性的典型页岩气层作为对比层，
选择区域经试气证实的页岩气层作为已知井对比层；选取对比层录井、测井解释关键参数孔隙度POR(％)，有机碳含量TOC(％)，全烃或甲烷含量QT(％)，总含气量Gt(％)，含气饱和度Sg(％)，脆性指数BRIT(％)，厚度H(m)、镜质体反射率Ro(％)、钻时ROP或钻时比值RROP、烃对比系数Kc、地层压力FP或地层压力梯度FPG中的3至10个，作为定量类比评价关键参数。
2)建立待解释井对比层解释关键参数数据库，
以待解释层段为对比对象，选择试气已证实为气层的井，在对比层段内按一定步长读取录井、测井解释关键参数，以文本文件或WIS文件格式建立数据库；根据气层产量可多选几个气层段对比，最好是选择一个达到接近工业气流产能的气层段，一般不选用试气已证实为含气层的井。
第二步：读取解释关键参数
1)优选解释关键参数
根据资料录取及地区地质特点优选解释关键参数，在涪陵页岩气田或建南气田多选用孔隙度POR，有机碳含量TOC，全烃或甲烷含量QT，总含气量Gt，含气饱和度，脆性指数BRIT，厚度H等7个参数。
2)划分待解释井页岩气异常显示层段
依据录井地层岩性、甲烷、钻时、有机碳含量和测井自然伽马、密度曲线变化判断页岩气显示储层，作为待解释层段。储层岩性应为页岩或泥岩，具有纹理构造，甲烷高、有机碳含量高、钻时低、自然伽马高、密度低特征，粒径小于0.004mm，甲烷含量不低于1.0％，有机碳含量不低于0.5％。
以不扣除垂厚小于2m薄夹层为原则，划分页岩气储层评价井段，并作为待解释层段。
3)读取待解释层段录井与测井关键参数
页岩气储层待解释层段、已知井对比层段解释关键参数一般按层段取值，也可以在层段内按一定步长取值；
按层段取值：一般是一个解释层段读取一组平均数据，参数分别为孔隙度POR2、有机碳含量TOC2、全烃或甲烷含量QT2、总含气量Gt2、含气饱和度Sg2、脆性指数BRIT2、储层厚度H2及钻时ROP2或钻时比值RROP2、地层压力梯度FPG2等。
按步长取值：在同一层段内，一般按步长1m连续取值计算，该数值为步长内的平均值，参数分别为孔隙度POR2、有机碳含量TOC2、全烃或甲烷含量QT2、总含气量Gt2、含气饱和度Sg2、脆性指数BRIT2、储层厚度H2及钻时ROP2或钻时比值RROP2、地层压力梯度FPG2等。
4)读取已知井对比层段录井与测井关键参数
对应待解释层段取值方式，读取参数孔隙度POR1、有机碳含量TOC1、全烃或甲烷含量QT1、总含气量Gt1、含气饱和度Sg1、脆性指数BRIT1、储层厚度H1及钻时ROP1或钻时比值RROP1、地层压力梯度FPG1。
第三步：处理解释评价参数
1)计算类比评价值
类比评价值具体公式如下：
IPOR＝POR2/POR1，
ITOC＝TOC2/TOC1，
IQT＝QT2/QT1，
IGt＝Gt2/Gt1，
ISg＝Sg2/Sg1，
IBRIT＝BRIT2/BRIT1，
IH＝H2/H1,
IROP＝ROP1/ROP2,或IROP＝RROP2/RROP1
IFPG＝FPG2/FPG1。
式中，IPOR、ITOC、IQT、IGt、ISg、IBRIT、IH、IROP、IFPG分别为孔隙度、总有机碳含量、烃含量、总含气量、含气饱和度、脆性指数、储层厚度、钻时、地层压力梯度类比评价值，均用小数表示，无量纲。
一般多选用IPOR、ITOC、IQT、IGt、ISg、IBRIT、IH等7项参数进行解释，在同一区块钻时、地层压力梯度变化不明显，用不同区块资料类比时应增选地层压力梯度、钻时或钻时比值。
2)计算多参数定量类比评价指数Ia
类比评价指数选用几何平均数表示。所谓的几何平均数是N个变量值连乘积的n次方根，适用于计量单位不统一的各类值的平均值，选用7个参数变量时n＝7，以此类推。
7个参数类比评价指数Ia的计算公式为：
Ia＝(IPOR×ITOC×IQT×IGt×ISg×IBRIT×IH)1/7
第四步：绘制解释评价图版
1)绘制多参数类比图版
以孔隙度、有机碳含量、甲烷、总含气量、含气饱和度、脆性指数、厚度等7个不同类别的评价参数为横轴，评价参数数值为纵轴，可以选用线性、对数或双曲线等刻度形式，绘制多参数类比图版，一般选用线性刻度形式；为了方便对比，使各参数间的分布范围不宜过大，可以将各参数进行整数倍放大或缩小，一般选用10倍放大或缩小。
分别将读取到的已知井对比层段、待评价井解释层段数值投射到多参数类比图上，已知井的对比层段的数据用空心方块虚线标记，待评价井解释层段用实心方块实线标记，待解释层段与对比层段差值用涨跌柱线标记，待解释层段在已知井对比层段之上部分用上涨柱线实心填充标识(见图3)；待解释层段在已知井对比层段之下部分用下跌柱线空心标识，可以用目测法观察和对比其相似性及差异性，可以直观定性评价待解释层的含气性。
2)绘制多参数相关性判别图版及计算相关强度R2
以对比层段解释参数值为横轴，待解释层段参数值为纵轴，均选用线性刻度，绘制多参数相关性判别图版，并做线性回归分析，得出线性相关方程及相关强度R2(见图4)。
3)绘制类比评价指数Ia—相关强度R2评价交会图
以相关强度R2为横轴，类比评价指数Ia为纵轴，绘制Ia—R2交会图(见图5)；
将区域内经过试气证实的已知井气层、含气层数据点绘到Ia—R2交会图上，根据数据统计原理确定解释评价标准，即确定不同储层含气类别界线AB与CD，界线AB与CD将图版分为页岩气层(Ⅰ)、含气层(Ⅱ)两个区域，用于表征页岩气储层评价结果；
将待解释层段的Ia与R2数据点置于Ia—R2交会图版中。
统计四川盆地及周缘地区页岩气典型气层、含气层12口井30层段类比评价结果，Ⅰ、Ⅱ区的评价标准是：Ⅰ区气层，相关强度R2大于或等于0.5，同时类比指数Ia大于或等于0.8；Ⅱ区含气层，相关强度R2小于0.5，或类比指数Ia小于0.8。
第五步：输出评价结果
待解释层段数据交会点所处区域即为对应评价结果：交会点落在Ⅰ区为气层，试气能获得工业气流；落在Ⅱ区为含气层，试气不能获得工业气流。
下面用实例来说明本发明的实现方法与步骤：
实例一：川东鄂西地区JYHF-1井
JYHF-1井是川东鄂西地区某构造上的重点页岩气试验井，分别进行了测井和录井工作，完成导眼段钻井后，根据测录井解释结果进行了侧钻水平井。在导眼井侏罗系东岳庙段发现并解释页岩气层1层段厚52.0m。具体解释评价步骤是：
第一步：建立解释数据库
选择该区域经试气证实的侏罗系J111井的典型页岩气层作为已知对比层。J111井在侏罗系东岳庙段598.0～646.0m井段，厚48.0m，完井压裂试气，日产天然气2100m3～3900m3。获取已知井对比层解释关键参数孔隙度、有机碳含量、全烃含量、甲烷含量、总含气量、含气饱和度、脆性指数、地层压力梯度、厚度共9个；
POR1＝3.8％、TOC1＝1.55％、Ct1＝9.85％、QT1＝8.14％、Gt1＝1.28％、Sg1＝60.0％、BRIT1＝50.0％、FPG1＝1.07MPa/100m、H1＝48.0m。
第二步：读取解释关键参数
1)解释关键参数优选
优选POR1＝3.8％、TOC1＝1.55％、QT1＝8.14％、Gt1＝1.28％、Sg1＝60.0％、BRIT1＝50.0％、H1＝48.0m共7个。
2)划分待解释井页岩气异常显示层段
依据录井地层岩性、甲烷、钻时、有机碳含量和测井自然伽马、密度曲线变化判断JYHF-1井页岩气异常显示层段，判断侏罗系东岳庙段591.0～643.0m井段，厚52.0m，岩性主要为页岩，甲烷由基值0.03％上升到8.46％，有机碳含量平均1.6％；该段作为待解释层段。
3)读取待解释层录井与测井关键参数
待解释层段按层段读取一组平均数据，
对应读取JYHF-1井页岩异常显示待解释层段解释关键参数孔隙度、有机碳含量、甲烷、总含气量、含气饱和度、脆性指数、厚度共7个。
POR2＝3.6％、TOC2＝1.6％、QT2＝8.46％、Gt2＝1.27％、Sg2＝58.0％、BRIT2＝50.0％、H2＝52.0m。
4)读取已知井对比层段录井与测井关键参数
POR1＝3.8％、TOC1＝1.55％、QT1＝8.14％、Gt1＝1.28％、Sg1＝60.0％、BRIT1＝50.0％、H1＝48.0m共7个。
第三步：处理解释评价参数
1)计算类比评价值
根据公式
IPOR＝POR2/POR1＝3.6/3.8＝0.94，
ITOC＝TOC2/TOC1＝1.6/1.55＝1.03，
IQT＝QT2/QT1＝8.46/8.14＝1.04，
IGt＝Gt2/Gt1＝1.27/1.28＝0.99，
ISg＝Sg2/Sg1＝58.0/60.0＝0.96，
IBRIT＝BRIT2/BRIT1＝50.0/50.0＝1，
IH＝H2/H1＝52.0/48.0＝1.08,
2)计算多参数定量类比评价指数Ia
Ia＝(IPOR×ITOC×IQT×IGt×ISg×IBRIT×IH)1/7＝(0.94×1.03×1.04×0.99×0.96×1×1.08)1/7＝1.01
第四步：绘制解释评价图版
1)绘制多参数类比图版
JYHF-1井591.0～643.0m井段与J111井598.0～646.0m井段绘制多参数类比解释图版(见图6)，类比图版显示，JYHF-1井页岩气异常显示待解释层与已知井对比层形状大致接近，待解释层全烃QT值、厚度H值与已知井对比层相比略高，其余解释关键参数均与参考层接近。
2)绘制多参数相关性判别图版及计算相关强度R2
以对比层段解释参数值为横轴(x)，待解释层段参数值为纵轴(y)，绘制多参数相关性判别图版(见图7)，由线性回归分析，得到线性相关方程为：
y＝1.03x-0.79、相关强度R2＝0.99。
3)绘制类比评价指数Ia——相关强度R2评价交会图
以相关强度R2为横轴(x)，类比评价指数Ia为纵轴(y)，绘制Ia——R2交会图(见图8)；将区域内经过试气证实的已知井气层、含气层数据点绘到Ia——R2交会图上，将JYHF-1井591.0～643.0m井段的Ia与R2数据点置于Ia—R2交会图版中，待解释层交会点落在气层区(Ⅰ区)。
第五步：输出评价结果
JYHF-1井导眼段的591.0～643.0m井段解释为气层，厚52.0m，侧钻水平井能够获得工业气流。
完井后，对JYHF-1井导眼段气层侧钻水平井，水平段长1022.0m，分7段压裂试气，日产气1.16×104m3，测试结论为气层，说明本发明解释结果与试气验证符合。
实例二：中扬子地区JY2井
JY2井是中扬子地区川东南某构造上的重点页岩气预探井，分别进行了测井和录井工作，在导眼井志留系龙马溪组下部—奥陶系五峰组发现并解释页岩气层1层段厚98.0m。具体解释评价步骤是：
第一步：建立解释数据库
选择该区域经试气证实的JY1井志留系龙马溪组下部—奥陶系五峰组典型页岩气层作为已知对比层。JY1井2326.0～2415.0m井段，侧钻水平井，水平段长1008m，分15段压裂试气，日产天然气21.3×104m3。
获取已知井对比层解释关键参数孔隙度、有机碳含量、全烃含量、甲烷含量、总含气量、含气饱和度、脆性指数、地层压力梯度、厚度共9个；
POR1＝5.2％、TOC1＝2.52％、Ct1＝2.13％、QT1＝2.03％、Gt1＝2.99％、Sg1＝80.0％、BRIT1＝60.0％、FPG1＝1.45MPa/100m、H1＝89.0m。
第二步：读取解释关键参数
1)解释关键参数优选
优选POR1＝5.2％、TOC1＝2.52％、QT1＝2.03％、Gt1＝2.99％、Sg1＝80.0％、BRIT1＝60.0％、H1＝89.0m共7个。
2)划分待解释井页岩气异常显示层段
依据录井地层岩性、甲烷、钻时、有机碳含量和测井自然伽马、密度曲线变化判断JY2井页岩气异常显示层段，判断志留系龙马溪组下部—奥陶系五峰组2477.0～2575.0m井段，厚98.0m，岩性主要为碳质页岩，甲烷由基值0.10％上升到4.30％，该段有机碳含量平均2.23％；该段作为待解释层段。
3)读取待解释层录井与测井关键参数
待解释层段按层段读取一组平均数据，
对应读取JY2井页岩异常显示待解释层段解释关键参数孔隙度、有机碳含量、甲烷、总含气量、含气饱和度、脆性指数、厚度共7个。
POR2＝5.7％、TOC2＝2.23％、QT2＝4.30％、Gt2＝2.99％、Sg2＝82.0％、BRIT2＝60.0％、H2＝98.0m。
4)读取已知井对比层段录井与测井关键参数
POR1＝5.2％、TOC1＝2.52％、QT1＝2.03％、Gt1＝2.99％、Sg1＝80.0％、BRIT1＝60.0％、H1＝89.0m共7个。
第三步：处理解释评价参数
1)计算类比评价值
根据公式
IPOR＝POR2/POR1＝5.7/5.2＝1.09，
ITOC＝TOC2/TOC1＝2.23/2.52＝0.88，
IQT＝QT2/QT1＝4.30/2.03＝2.12，
IGt＝Gt2/Gt1＝2.99/2.99＝1，
ISg＝Sg2/Sg1＝82.0/80.0％＝1.02，
IBRIT＝BRIT2/BRIT1＝60.0/60.0＝1，
IH＝H2/H1＝98.0/89.0＝1.10,
2)计算多参数定量类比评价指数Ia
Ia＝(IPOR×ITOC×IQT×IGt×ISg×IBRIT×IH)1/7＝(1.09×0.88×2.12×1×1.02×1×1.10)1/7＝1.13
第四步：绘制解释评价图版
1)绘制多参数类比图版
JY2井2477.0～2575.0m井段与JY1井2326.0～2415.0m井段绘制多参数类比解释图版(见图9)，类比解释图版显示，JY2井页岩气异常显示待解释层与已知井对比层形状大致接近，待解释层孔隙度POR值、全烃QT值、含气饱和度Sg值、厚度H值与已知井对比层相比略高，有机碳含量TOC值与已知对比层相比略低，脆性指数BRIT值、总含气量Gt值与已知对比层数据一致。
2)绘制多参数相关性判别图版及计算相关强度R2
以对比层段解释参数值为横轴(x)，待解释层段参数值为纵轴(y)，绘制多参数相关性判别图版(见图10)，由线性回归分析，得到线性相关方程为：
y＝0.95x+7.58、相关强度R2＝0.90。
3)绘制类比评价指数Ia——相关强度R2评价交会图
以相关强度R2为横轴(x)，类比评价指数Ia为纵轴(y)，绘制Ia—R2交会图(见图11)；将区域内经过试气证实的已知井气层、含气层数据点绘到Ia—R2交会图上，将JY2井2477.0～2575.0m井段的Ia与R2数据点置于Ia—R2交会图版中，待解释层交会点落在气层区(Ⅰ区)。
第五步：输出评价结果
JY2井2477.0～2575.0m井段，厚98.0m，试气能获工业气流，气层产能高于JY1井。
随后，针对上述解释的气层中下部2450.0～2568.0m井段进行了水平井侧钻，水平段长1500.0m，分20段压裂试气，日产天然气35.0×104m3，说明本发明解释结果与试气验证符合。
Claims (5)
1.页岩气储层多参数定量类比评价方法，其特征在于具体步骤如下：
第一步：建立解释数据库
1)建立已知井对比解释层段关键参数数据库，
选择已试气并具有区域代表性的典型页岩气层作为对比层,
选择区域经试气证实的页岩气层作为已知井对比层；选取对比层录井、测井解释关键参数孔隙度POR，有机碳含量TOC，全烃或甲烷含量QT，总含气量Gt，含气饱和度Sg和脆性指数BRIT,其单位为％，厚度H、镜质体反射率Ro,单位为％,钻时ROP或钻时比值RROP、烃对比系数Kc、地层压力FP或地层压力梯度FPG中的3至10个，作为定量类比评价关键参数；
2)建立待解释井对比层解释关键参数数据库，
以待解释层段为对比对象，选择试气已证实为气层的井，在对比层段内按定步长读取录井、测井解释关键参数，以文本文件或WIS文件格式建立数据库；根据气层产量选择至少一个气层段对比；
第二步：读取解释关键参数
1)选择解释关键参数，
根据资料录取及地区地质特点选择解释关键参数；
2)划分待解释井页岩气异常显示层段，
依据录井地层岩性、甲烷、钻时、有机碳含量和测井自然伽马、密度曲线变化判断页岩气显示储层；以不扣除垂厚小于2m薄夹层为原则，划分页岩气储层评价井段，并作为待解释层段；
3)读取待解释层段录井与测井关键参数，
一个解释层段读取一组典型数据，参数分别为孔隙度POR2、有机碳含量TOC2、全烃或甲烷含量QT2、总含气量Gt2、含气饱和度Sg2、脆性指数BRIT2、储层厚度H2及钻时ROP2或钻时比值RROP2、地层压力梯度FPG2；
4)读取已知井对比层段录井与测井关键参数，
对应待解释层段取值方式，读取参数孔隙度POR1、有机碳含量TOC1、全烃或甲烷含量QT1、总含气量Gt1、含气饱和度Sg1、脆性指数BRIT1、储层厚度H1及钻时ROP1或钻时比值RROP1、地层压力梯度FPG1；
第三步：处理解释评价参数
1)计算类比评价值
类比评价值具体公式如下：
IPOR＝POR2/POR1，
ITOC＝TOC2/TOC1，
IQT＝QT2/QT1，
IGt＝Gt2/Gt1，
ISg＝Sg2/Sg1，
IBRIT＝BRIT2/BRIT1，
IH＝H2/H1,
IROP＝ROP1/ROP2,或IROP＝RROP2/RROP1
IFPG＝FPG2/FPG1；
式中，IPOR、ITOC、IQT、IGt、ISg、IBRIT、IH、IROP、IFPG分别为孔隙度、总有机碳含量、烃含量、总含气量、含气饱和度、脆性指数、储层厚度、钻时、地层压力梯度类比评价值，均用小数表示，无量纲；
2)计算类比评价指数Ia，采用几何平均数算法计算；
第四步：绘制解释评价图版
1)绘制多参数类比图版，
以孔隙度、有机碳含量、甲烷、总含气量、含气饱和度、脆性指数、厚度不同类别的评价参数为横轴，评价参数数值为纵轴，绘制多参数类比图版；
分别将读取到的对比层、待解释层数值投射到多参数类比图上，用目测法观察和对比其相似性及差异性，直观定性评价待解释层段的含气性；
2)绘制多参数相关性判别图版及计算相关强度R2，
以对比层段参数值为横轴，待解释层段参数值为纵轴，绘制多参数相关性判别图版，并做线性回归分析，得出线性相关方程及相关强度R2；
3)绘制类比评价指数Ia—相关强度R2评价交会图，
以相关强度R2为横轴，类比评价指数Ia为纵轴，绘制Ia—R2交会图；
将区域内经过试气证实的已知井气层、含气层数据点绘到Ia—R2交会图上，根据数据统计原理确定解释评价标准，即确定不同储层含气类别界线AB与CD，界线AB与CD将图版分为页岩气层Ⅰ、含气层Ⅱ两个区域，用于表征页岩气储层评价结果；
将待解释层段的多参数类比评价指数Ia、相关强度R2数据点置于Ia—R2交会图版中；
第五步：输出评价结果
待解释层段数据交会点所处区域即为对应评价结果：Ⅰ区为气层，Ⅱ区为含气层。
2.根据权利要求1所述的页岩气储层多参数定量类比评价方法，其特征在于根据气层产量选择一个达到接近工业气流产能的气层段。
3.根据权利要求1所述的页岩气储层多参数定量类比评价方法，其特征在于解释评价取值按层段取值，每一层段取一组平均数据；按步长取值，在同一层段内，按步长为1m连续取值计算。
4.根据权利要求1所述的页岩气储层多参数定量类比评价方法，其特征在于绘制多参数类比图版各参数10倍放大或缩小。
5.根据权利要求1所述的页岩气储层多参数定量类比评价方法，其特征在于储层岩性为页岩或泥岩，粒径小于0.004mm，甲烷含量不低于1.0％，有机碳含量不低于0.5％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410446356.8A CN104239703B (zh) | 2014-09-03 | 2014-09-03 | 页岩气储层多参数定量类比评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410446356.8A CN104239703B (zh) | 2014-09-03 | 2014-09-03 | 页岩气储层多参数定量类比评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104239703A CN104239703A (zh) | 2014-12-24 |
| CN104239703B true CN104239703B (zh) | 2017-04-12 |
Family
ID=52227753
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410446356.8A Active CN104239703B (zh) | 2014-09-03 | 2014-09-03 | 页岩气储层多参数定量类比评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104239703B (zh) |
Families Citing this family (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105986813B (zh) * | 2015-02-13 | 2019-05-10 | 中国石油化工股份有限公司 | 准致密储层快速评价方法及准致密储层多指标评价方法 |
| CN105989541B (zh) * | 2015-03-05 | 2019-07-05 | 中国石油化工股份有限公司 | 一种基于大数据的测井资料的处理方法 |
| CN105257287A (zh) * | 2015-11-05 | 2016-01-20 | 成都理工大学 | 油页岩石油储层参数评价方法 |
| CN106932836B (zh) * | 2015-12-30 | 2019-07-12 | 中国石油化工股份有限公司 | 一种用于评价页岩气含气丰度的方法及系统 |
| CN107367755A (zh) * | 2016-05-11 | 2017-11-21 | 中国石油化工股份有限公司 | 一种改进的多参数交会图绘制方法 |
| CN107451310B (zh) * | 2016-05-31 | 2020-09-04 | 中国石油化工股份有限公司 | 基于页岩源储相互关系的分类评价方法及装置 |
| CN107689013B (zh) * | 2016-08-03 | 2021-06-15 | 中国石油化工股份有限公司 | 圈闭资源量钻前评价参数校正方法及装置 |
| CN106646594A (zh) * | 2016-09-09 | 2017-05-10 | 中国海洋石油总公司 | 基于“断‑砂”耦合关系的浅层岩性油气藏定量预测方法 |
| CN106529199B (zh) * | 2016-12-15 | 2019-01-04 | 中国石油新疆油田分公司勘探开发研究院 | 一种砾岩油藏化学驱井距的确定方法 |
| CN108240952A (zh) * | 2016-12-24 | 2018-07-03 | 中石化石油工程技术服务有限公司 | 一种解析法计算页岩含气量的方法 |
| CN108661630B (zh) * | 2017-03-31 | 2021-05-11 | 中国石油化工股份有限公司 | 一种基于参数优选的地质甜点定量评价方法 |
| CN107701172B (zh) * | 2017-09-22 | 2020-07-24 | 中石化石油工程技术服务有限公司 | 基于线性模型的页岩气水平井初期最高产能的预测方法 |
| CN109580906B (zh) * | 2017-09-28 | 2021-08-24 | 中国石油化工股份有限公司 | 基于岩石物理的页岩脆性识别图版的制作方法及系统 |
| CN107977480B (zh) * | 2017-10-18 | 2021-04-30 | 中国石油化工集团有限公司 | 一种页岩气储层产气性快速评价方法 |
| CN108303510A (zh) * | 2017-12-25 | 2018-07-20 | 中国石油天然气股份有限公司 | 页岩气储层性能的评价方法、装置及计算机存储介质 |
| CN108182531A (zh) * | 2017-12-27 | 2018-06-19 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩气开发评价方法、装置及终端设备 |
| CN108363114B (zh) * | 2018-01-12 | 2019-11-08 | 中国石油天然气股份有限公司 | 致密油甜点区评价方法及装置 |
| CN108301823B (zh) * | 2018-01-19 | 2021-05-28 | 捷贝通石油技术集团股份有限公司 | 一种识别储层油气甜点的方法 |
| CN110322363A (zh) * | 2018-03-29 | 2019-10-11 | 中国石油化工股份有限公司 | 页岩气储层改造体积计算方法及系统 |
| CN108915678A (zh) * | 2018-08-15 | 2018-11-30 | 中海石油(中国)有限公司 | 一种大西洋两岸油气田开发指标区域类比方法 |
| CN109991123B (zh) * | 2019-03-28 | 2021-07-06 | 中国石油化工股份有限公司 | 页岩油资源可动性的地球化学评价方法 |
| CN110991933A (zh) * | 2019-12-19 | 2020-04-10 | 西南石油大学 | 一种山地页岩气资源评价方法及系统 |
| CN111475685B (zh) * | 2019-12-27 | 2020-10-27 | 北京国双科技有限公司 | 油气勘探方法、装置、存储介质及电子设备 |
| CN113358538B (zh) * | 2020-03-06 | 2022-11-01 | 中国石油天然气股份有限公司 | 确定岩石模态的方法及装置 |
| CN113627700B (zh) * | 2020-05-07 | 2024-04-26 | 中国石油化工股份有限公司 | 一种气井生产效果的评价方法及装置 |
| CN111997598B (zh) * | 2020-09-04 | 2024-03-26 | 长江大学 | 录井随钻储层评价方法、模型建立方法、装置及电子设备 |
| CN115965273A (zh) * | 2022-12-13 | 2023-04-14 | 中国石油大学(华东) | 一种页岩油水平井钻探过程中甜点评价方法 |
| CN116008512B (zh) * | 2023-03-02 | 2024-03-08 | 西南石油大学 | 一种针对未知页岩储层的判别含气情况的分析方法 |
| CN117216669B (zh) * | 2023-11-09 | 2024-03-22 | 贵州能源产业研究院有限公司 | 一种页岩储层分类评价图板的建立方法及应用 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102312671A (zh) * | 2011-09-16 | 2012-01-11 | 中国石油化工股份有限公司 | 能快速解释评价储层流体性质的方法 |
| CN102590889A (zh) * | 2012-02-17 | 2012-07-18 | 中国石油化工股份有限公司 | 基于雷达图与云模型的录井多参数油气解释方法 |
| CN103615242A (zh) * | 2013-12-17 | 2014-03-05 | 中国海洋石油总公司 | 实时地层流体录井多参数油气层综合解释评价方法 |
| MX2013001215A (es) * | 2013-01-30 | 2014-07-29 | Servicios Y Suministros En Informatica S A De C V | Método de exploración y manipulación tridimensional para simulaciones numéricas de modelos de flujo en yacimientos de aceite y gas lutita. |
-
2014
- 2014-09-03 CN CN201410446356.8A patent/CN104239703B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102312671A (zh) * | 2011-09-16 | 2012-01-11 | 中国石油化工股份有限公司 | 能快速解释评价储层流体性质的方法 |
| CN102590889A (zh) * | 2012-02-17 | 2012-07-18 | 中国石油化工股份有限公司 | 基于雷达图与云模型的录井多参数油气解释方法 |
| MX2013001215A (es) * | 2013-01-30 | 2014-07-29 | Servicios Y Suministros En Informatica S A De C V | Método de exploración y manipulación tridimensional para simulaciones numéricas de modelos de flujo en yacimientos de aceite y gas lutita. |
| CN103615242A (zh) * | 2013-12-17 | 2014-03-05 | 中国海洋石油总公司 | 实时地层流体录井多参数油气层综合解释评价方法 |
Non-Patent Citations (1)
| Title |
|---|
| 页岩气储层常规测井解释模型与应用实例;陈扬 等;《江汉石油职工大学学报》;20140131;第27卷(第1期);第20-24页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104239703A (zh) | 2014-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104239703B (zh) | 页岩气储层多参数定量类比评价方法 | |
| CN110644980B (zh) | 一种超低渗透油藏储层综合分类评价方法 | |
| CN102312671B (zh) | 能快速解释评价储层流体性质的方法 | |
| CN103912269B (zh) | 页岩气储层地层破裂压力梯度录井确定方法 | |
| CN104533400B (zh) | 一种重构测井曲线的方法 | |
| CN105930932B (zh) | 基于含气指数的页岩气层标准化无阻流量的获取方法 | |
| CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| CN103615243B (zh) | 一种利用地化衍生参数图版判断油质类型的方法 | |
| CN110130883A (zh) | 岩层参数的确定方法及装置 | |
| CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
| CN104695950A (zh) | 火山岩油藏产能预测方法 | |
| CN108374657B (zh) | 井断点自动识别方法 | |
| CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
| CN104569343A (zh) | 页岩储层脆性矿物测井定量表征方法 | |
| CN104428492A (zh) | 使用矿物学确定压裂刺激点的布置的系统、方法、和计算机程序产品 | |
| CN115030707A (zh) | 一种油页岩“甜点”快速评价方法 | |
| CN106560820A (zh) | 一种页岩气储层的测井评价方法 | |
| CN112145165A (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
| CN105350959A (zh) | 由测井岩性密度确定页岩气储层含气饱和度的方法 | |
| CN109667576B (zh) | 一种高矿化度成因低阻油层测井识别方法 | |
| CN104463686A (zh) | 利用判别分析法随钻识别页岩气储层的方法 | |
| CN105464655B (zh) | 一种流体测井识别方法 | |
| CN105064987A (zh) | 利用随钻录井q参数进行水层识别的解释评价方法 | |
| CN109738955B (zh) | 一种基于成分-结构分类下的变质岩岩性综合判别方法 | |
| CN112832738A (zh) | 一种碎屑岩累计生烃强度确定方法及甜点层的识别与评价方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180904 Address after: 433123 235 Xiangyang 57 Avenue, Qianjiang, Hubei Co-patentee after: SINOPEC JIANGHAN PETROLEUM ENGINEERING CO., LTD. Patentee after: Measurement Logging Co, Ltd. of SINOPEC Jianghan Petroleum Engineering Co., Ltd. Address before: 433123 Qianjiang City of administrative units directly under the control of the provinces of Hubei Province face south five or seven main road 235 Patentee before: Measurement Logging Co, Ltd. of SINOPEC Jianghan Petroleum Engineering Co., Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200717 Address after: 100728 Beijing, Chaoyangmen, North Street, No. 22, No. Co-patentee after: SINOPEC OILFIELD SERVICE Corp. Patentee after: SINOPEC Group Co-patentee after: SINOPEC JIANGHAN OILFIELD SERVICE Corp. Address before: 433123 No. 57, No. 235, Xiangyang Road, Hubei, Qianjiang Co-patentee before: SINOPEC JIANGHAN OILFIELD SERVICE Corp. Patentee before: Logging company of Sinopec Jianghan Petroleum Engineering Co.,Ltd. |