CN108798657B - 基于钻井液录井参数气测值的录井解释方法 - Google Patents

Abstract

本发明公开了一种基于钻井液录井参数气测值的录井解释方法，包括如下步骤：a、对钻井液录井采集到的气测值进行预处理；b、将预处理后的气测值、钻时、电导率和温度计算变化指数；c、通过变化指数计算含水指数W、含油指数O和含气指数G；d、建立以含水指数W、含油指数O和含气指数G相对比值为三边的三角图版；e、利用建立的三角图版进行储层含流体性质识别。本发明结合曲线法、图版法、自动识别技术和流体一体识别技术，集成度高、人工干预少且便于实现。

Description

基于钻井液录井参数气测值的录井解释方法

技术领域

本发明涉及一种基于钻井液录井参数利用集成三角图版和概率图道进行储层含流体性质判识的录井解释方法，属于机电领域石油天然气工业勘探开发范畴随钻录井资料处理技术应用领域。

背景技术

目前，针对钻井过程中钻遇地层储层及其含流体性质的录井解释方法国内外均有较为深入的研究，现分述如下：

国内方面，中石化倾向于新型装备的研发和实践，开展岩石热解、热蒸发烃、罐顶气和核磁共振等新型油气储层识别和发现装备的研究和应用，天然气方面，中石化较为推崇分析常规综合录井探测气组分NC₅之外的C₈以后的更复杂组分，值得注意的是，中石化在气测数据利用上重视定量脱气器采集、恒温恒压气样传输管路和地化色谱仪(GCMS)结合的高灵敏度探测样式，甚至开发了不同于中石油大规模列装的气相色谱和荧光色谱的离子色谱和拉曼色谱，油藏勘探方面，中石化大量采用新型地化技术并在定量荧光方面有一定的实践经验和成果。

中石油方面，长于天然气藏勘探服务的川庆地研院开发了一系列基于综合录井的录井解释方法，既有数值模拟类的“利用微钻时数据随钻计算录井孔隙度进行储层评价的方法(201210427291.3)、图版法类的“利用判别分析建立气测图版对储层流体进行判断的方法(201210427314.0)”，又有正演范畴的“利用随钻钻井液录井资料进行储层流体类型的判别方法(专利授权公告号：CN 103806911 B)”，此外，也有一批基于大数据处理的基于判别分析、神经网络和支持向量机的机器深度学习解释法，上述方法在利用录井参数进行流体性质识别方面取得了一定成效，但要实现自动化、智能化必须有所改进。

“利用微钻时数据随钻计算录井孔隙度进行储层评价的方法(20121 0427291.3)”基于大量测录井既有井的统计归纳和分析总结，系在大量工程和试验基础上的经验方法总结，应用过程中凸显出难以应用到新区块、仅能计算孔隙度对流体表征能力弱的问题；而“利用判别分析建立气测图版对储层流体进行判断的方法(201210427314.0)”基于图版法建立，关注气水关系识别，但对含油储层关注不够，此外，对于新区块、新层系需要大幅更新图版并调整判别区域，判别方式主观性较强，在探井区域资料井较少的情况下适应性较差。此外，值得注意的是，上述发明对于计算机编程实现没有做出针对性的设计，主要依靠人工统计和更新算法和图版，后期维护时间和人力成本较高。

综上，现将现有技术方案存在的问题归纳如下：

1、现有油井录井解释技术方案缺乏曲线法技术，北方油田服务公司长期以来图版法，如气测图版法、轻烃色谱图版法等进行分析，北方油田业主如长庆油田对此并不是很满意，由于录井资料随钻录取、干扰因素多、难以修正的特点，类似图版法的依赖数据具有较强统一性、系统误差较稳定的方法难以充分发挥录井原始资料作用；

2、现在推广大部分曲线法，通过趋势值分析进行定性识别，在川渝、长庆和塔里木取得一定进展，但其适应强、效率高的特性也带来了针对性偏弱、区块优化差的问题，因此，需要进行一定的修正和调整。

3、现有的图版法和曲线法共同的问题在于难以进行计算机应用，图版法需要频繁调整分离区间，曲线法需要人工定性判断，这都制约了录井解释自动化改进的进一步发展。

4、现有技术更多的关注气井或者油井的解释，鲜有油气水干层一体识别的技术。

发明内容

本发明的目的在于克服现有技术存在的上述问题，提供一种基于钻井液录井参数气测值的录井解释方法。本发明结合曲线法、图版法、自动识别技术和流体一体识别技术，集成度高、人工干预少且便于实现。

为实现上述目的，本发明采用的技术方案如下：

一种基于钻井液录井参数气测值的录井解释方法，其特征在于，包括如下步骤：

a、对钻井液录井采集到的气测值进行预处理；

b、将预处理后的气测值、钻时、电导率和温度计算变化指数；

c、通过变化指数计算含水指数W、含油指数O和含气指数G；

d、建立以含水指数W、含油指数O和含气指数G相对比值为三边的三角图版；

e、利用建立的三角图版进行储层含流体性质识别。

所述步骤a中，通过校正气测值，使单根气及频率高于单根气的权重较低或剔除干扰数据。

所述校正气测值为：设钻杆单根长度区间为[8,10]，则T单根∈[8,10]，f单根>0.125，以此为截止频率f_t，

a_i＝f_t/2πn_i＝0.1989,i＝1,2,3…j,i∈N

则用迭代公式计算

Y_i＝a_iX_i+(1-a_i)Y_i-1

其中，Y_i----------当前井深i点校正值；

X_i---------当前井深i点读值；

Y_i-1----------当前井深i点读值；

a_i-----------钻具系数,根据钻具长度取值，自然数组，无量纲；

n_i-----------钻具长度，m；

其中，a_i与单根长度有关。

所述气测值包括全烃组分值TG，单位％或ppm；甲烷组分值C₁，单位％或ppm；乙烷组分值C₂，单位％或ppm。

所述步骤b中，通过全烃组分值T、甲烷组分值C₁、乙烷组分值C₂、出口电导率值S、出口温度值T，分别计算甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil，再通过甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil计算含水指数W、含油指数O和含气指数G。

所述甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil分别通过下述公式计算得出：

其中，TG为色谱检测到的全烃组分值，单位％，C₁为色谱检测到的甲烷组分值，单位％，C₂为色谱检测到的乙烷组分值，单位％，WOB为钻压，单位kN·m，ROP为钻时，单位min/m，μ_C1为待求深度点所在井段对应甲烷组分平均值，μ_s为待求深度点所在井段对应出口电导率平均值，μ_T为待求深度点所在井段对应出口温度平均值，μ_TG-C1-C2为待求深度点所在井段对应出口C₁、C₂之外烃重组分平均值，△_C1为待求深度点所在井段对应出口电导率标准差，δ_s为待求深度点所在井段对应出口电导率标准差，△_T为待求深度点所在井段对应出口温度标准差，△_TG-C1-C2为待求深度点所在井段对应出口C₁、C₂之外烃重组分标准差。

所述步骤c中，含水指数W、含油指数O和含气指数G分别通过下述公式计算得出：

设

则含水指数W计算方式如下：

设

则含气指数G计算方式如下：

设

则含油指数O计算方式如下：

所述步骤d具体包括如下步骤：

d1、建立图版框

建立三角图版，边长L代表油Oil、气Gas、水Water组分占钻井液体积理论最大百分含量，即L最大值L_max＝1，将边长等分为三段，即有

L＝L₁+L₂+L₃

设定有一正圆O₁与此正三角形各边三等分点均相交，其半径：

有与该圆同心另一圆形O₂，其半径R₂

刻度为顺时针设置(0,1)，即按照顺时针方向，Oil边自下而上，为零点和1点，Water边自上而下，分别为零点和1点，Gas边自右及左，依次为零点和1点。

根据W系数、G系数和O系数结果，分别在三角图版各对应边分别标注，分别作过三点与其所在轴零刻度衔接的邻边之平行线，三线必然两两相交，围成一个正三角形，其重心落在图版上对应的区间则判定该异常井段的储层发育情况及流体性质；

d2、IGIC坐标变换计算

将二维笛卡尔坐标系转换到三角图版设定的三角坐标系，设含水指数为L_w，含油指数为L_o，含气指数L_g，

P_o(x_o,y_o):

P_w(x_w,x_w):

P_g(x_g,y_g):

x_g＝1-L_g

y_g＝0

由于平行于L_g边的辅助线必然斜率为0，故P_g和P_w所在辅助线的交点为三角形顶点，而后三角形中心P(x,y)的坐标分别为：

由此各深度点对应P坐标投影到三角解释集成图版中，或根据图版特征直接判断P点所在的流体识别区，将待解释井段的对应各深度点对应P点值投影到三角图版中根据散点的展布判断出该段可能钻遇某种流体，若散点分布在两个流体识别区则表明可能为两种流体并存。

采用本发明的优点在于：

本发明技术优势体现在三方面：

1、本发明基于计算机自动化设计和构建的计算方法和图版架构，集成三角图版和概率图道以及与其配套的处理方法、计算方法和判定方法都建立在计算机自动识别基础上，从预处理、计算准备、图版构建、判别模式设定到计算机自动化准备和实现的录井解释技术集成解决方案，集成度高、人工干预少以及便于计算机编程实现。

2、本发明集成三角图版系首次提出基于钻井液录井参数流体识别集成图版，从集成程度和信息归纳程度看，目前尚无能够集成电导率、温度、甲烷组分、重组分、钻时和钻压，综合对含油储层、含气储层、含水储层和干层进行综合判断的集成三角图版法。

3、本发明利用数据周期性特征分析对录井参数进行预处理和校正，根据目前文献调研和专利查询情况，目前尚无与本发明类似地利用周期性特征对钻井液录井参数进行预处理和校正的技术公开。

附图说明

图1为气测值分解示意图；

图2为三角图版建立框架示意图；

图3为三角图版效果图；

图4为三角图版效果(水层)图；

图5为三角图版效果(干层)图；

图6为三角图版效果(气水层)图；

图7为三角图版自动化准备计算坐标变换示意图；

图8为概率图道判别模式示意图。

具体实施方式

实施例1

一种基于钻井液录井参数气测值的录井解释方法，包括如下步骤：

a、对钻井液录井采集到的气测值进行预处理；

b、将预处理后的气测值、钻时、电导率和温度计算变化指数；

c、通过变化指数计算含水指数W、含油指数O和含气指数G；

d、建立以含水指数W、含油指数O和含气指数G相对比值为三边的三角图版；

e、利用建立的三角图版进行储层含流体性质识别。

所述步骤a中，通过校正气测值，使单根气及频率高于单根气的权重较低或剔除干扰数据。

所述校正气测值为：设钻杆单根长度区间为[8,10]，则T单根∈[8,10]，f单根>0.125，以此为截止频率f_t，

a_i＝f_t/2πn_i＝0.1989,i＝1,2,3…j,i∈N

则用迭代公式计算

Y_i＝a_iX_i+(1-a_i)Y_i-1

其中，Y_i----------当前井深i点校正值；

X_i---------当前井深i点读值；

Y_i-1----------当前井深i点读值；

a_i-----------钻具系数,根据钻具长度取值，自然数组，无量纲；

n_i-----------钻具长度，m；

其中，a_i与单根长度有关。

所述气测值包括全烃组分值TG，单位％或ppm；甲烷组分值C₁，单位％或ppm；乙烷组分值C₂，单位％或ppm。

所述步骤b中，通过全烃组分值T、甲烷组分值C₁、乙烷组分值C₂、出口电导率值S、出口温度值T，分别计算甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil，再通过甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil计算含水指数W、含油指数O和含气指数G。

所述甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil分别通过下述公式计算得出：

其中，TG为色谱检测到的全烃组分值，单位％，C₁为色谱检测到的甲烷组分值，单位％，C₂为色谱检测到的乙烷组分值，单位％，WOB为钻压，单位kN·m，ROP为钻时，单位min/m，μ_C1为待求深度点所在井段对应甲烷组分平均值，μ_s为待求深度点所在井段对应出口电导率平均值，μ_T为待求深度点所在井段对应出口温度平均值，μ_TG-C1-C2为待求深度点所在井段对应出口C₁、C₂之外烃重组分平均值，△_C1为待求深度点所在井段对应出口电导率标准差，δ_s为待求深度点所在井段对应出口电导率标准差，△_T为待求深度点所在井段对应出口温度标准差，△_TG-C1-C2为待求深度点所在井段对应出口C₁、C₂之外烃重组分标准差。

所述步骤c中，含水指数W、含油指数O和含气指数G分别通过下述公式计算得出：

设

则含水指数W计算方式如下：

设

则含气指数G计算方式如下：

设

则含油指数O计算方式如下：

所述步骤d具体包括如下步骤：

d1、建立图版框

建立三角图版，边长L代表油Oil、气Gas、水Water组分占钻井液体积理论最大百分含量，即L最大值L_max＝1，将边长等分为三段，即有

L＝L₁+L₂+L₃

设定有一正圆O₁与此正三角形各边三等分点均相交，其半径：

有与该圆同心另一圆形O₂，其半径R₂

刻度为顺时针设置(0,1)，即按照顺时针方向，Oil边自下而上，为零点和1点，Water边自上而下，分别为零点和1点，Gas边自右及左，依次为零点和1点。

根据W系数、G系数和O系数结果，分别在三角图版各对应边分别标注，分别作过三点与其所在轴零刻度衔接的邻边之平行线，三线必然两两相交，围成一个正三角形，其重心落在图版上对应的区间则判定该异常井段的储层发育情况及流体性质；

d2、IGIC坐标变换计算

将二维笛卡尔坐标系转换到三角图版设定的三角坐标系，设含水指数为L_w，含油指数为Lo，含气指数L_g，

P_o(x_o,y_o):

P_w(x_w,x_w):

P_g(x_g,y_g):

x_g＝1-L_g

y_g＝0

由于平行于L_g边的辅助线必然斜率为0，故P_g和P_w所在辅助线的交点为三角形顶点，而后三角形中心P(x,y)的坐标分别为：

由此各深度点对应P坐标投影到三角解释集成图版中，或根据图版特征直接判断P点所在的流体识别区，将待解释井段的对应各深度点对应P点值投影到三角图版中根据散点的展布判断出该段可能钻遇某种流体，若散点分布在两个流体识别区则表明可能为两种流体并存。

为了进一步提高识别的准确性，还可以进一步做如下处理：

建立基于三角图版的概率图道，利用概率图道进行储层含流体性质识别。

所述概率图道的建立方法如下：

A、单点概率图道的建立

对于单深度点的情形，设Wi为某深度点对应含水指数，Oi为某深度点含油指数，Gi为某深度点含气指数，对于任意序号为i的深度点，则有：

其中，i＝1,2,3,…，n，i∈Z

B、组合深度点的某井段概率图道的建立

设Wi为某深度点对应含水指数，Oi为某深度点含油指数，Gi为某深度点含气指数，设某井段内有m个深度点，即该井段长m米，则对应的概率图道建立方法如下：

DW＝∑DWi；DO＝∑DOi；DG＝∑DGi；

其中，

i＝1,2,3,…，m，i∈Z

DW：对应井段相对含水可能性，单位％

DO：对应井段相对含油可能性，单位％

DG：对应井段相对含气可能性，单位％。

利用概率图道进行储层含流体性质识别的判别原则如下：

设Wi为某深度点对应含水指数，Oi为某深度点含油指数，Gi为某深度点含气指数，则对应建立的概率图道，对于单深度点的情形，直接比较DOi、DWi和DGi的大小关系，且按照以下原则进行分类：

当DOi或者DWi或者DGi其一小于16％时，则不参与最终流体性质判别；

当DOi或者DWi或者DGi其一大于66％时，判断对应井深赋存对应流体类型的可能性最大，且以其为主进行定名。

所述步骤f中，利用概率图道进行储层含流体性质识别的判别模式如下：

1、含油层识别方面：

如DOi大于66％，DWi小于16％，DGi大于25％时，解释该层为油气同层；

如DOi大于66％，DWi大于25％，DGi小于16％时，解释该层为油水同层；

如DOi大于66％且小于75％，同时DWi、DGi均小于16％时，解释该层为差油层；

如DOi大于75％时，且DWi、DGi均小于16％时，解释该层为油层；

对于油层考虑到原油的物理性质特性，适当下调66％基准至40％～50％；

2、含气层识别方面：

如DGi大于66％，DWi小于16％，DOi大于25％时，解释该层为油气同层；

如DGi大于66％，DWi大于25％，DOi小于16％时，解释该层为气水同层；特殊地，当DGi与DWi均未大于66％，同时DOi小于16％的情形，解释为气水同层；

如DGi大于66％且小于75％，同时DWi、DOi均小于16％时，解释该层为差气层；

如DGi大于75％时，且DWi、DOi均小于16％时，解释该层为气层；

对于气层，上浮66％基准至75％～80％；

3、含水层识别方面：

如DWi大于66％，DGi小于16％，DOi大于25％时，解释该层为油水同层；

如DWi大于66％，DGi大于25％，DOi小于16％时，解释该层为含气水层；特殊地，当DGi与DWi均未大于66％，同时DOi小于16％的情形，解释为气水同层；考虑气水相渗中的滑脱效应现象，在DWi、DGi同时参与判定时，上调DGi基准10％～15％；

如DWi大于85％时，且DGi、DOi均小于16％时，解释该层为水层。

实施例2

本实施例结合附图对本发明做进一步说明。

一种基于钻井液录井参数气测值的录井解释方法，包括

1)数据预处理

如图1所示，录井参数采集、传输和输出周期跨度大、受控因素多，造成其难以分析利用。本发明利用随钻接钻具(俗称接单根)这一典型的近似周期地影响因素，引入周期性分析和校正计算方法。下面将以气测值为例进行预处理计算示例。

设钻杆单根长度区间为[8,10]，则T单根∈[8,10]，f单根>0.125；以此为截止频率f_t,

a_i＝f_t/2πn_i＝0.1989,i＝1,2,3…j,i∈N

则可用迭代公式计算

Y_i＝a_iX_i+(1-a_i)Y_i-1

其中，Y_i----------当前井深i点校正值；

X_i---------当前井深i点读值；

Y_i-1----------当前井深i点读值；

a_i-----------钻具系数,根据钻具长度取值，自然数组，无量纲；

n_i-----------钻具长度，m。

其中，a_i与单根长度有关。采用上述算法校正气测值，使单根气及频率高于单根气的权重较低或干扰数据大部分被剔除，实现一定程度上的去单根效应和去噪、降噪。

其他录井参数，如温度、电导率均出现类似地与随钻接单根有关的周期性读值重复、跃迁和突变现象，均可使用上述计算方法进行预处理。

2变化指数计算

经过上述预处理计算，分别校正全烃组分值(TG，％或ppm)、甲烷组分值(C₁，单位％或者ppm)、乙烷组分值(C₂，％或ppm)出口电导率值(S,单位μs/cm或者μs/m)、出口温度值(T，单位℃或者℉)，并归一化、标准化；分别计算甲烷变化指数δ_C1、出口电导率变化指数δ_S、、出口温度变化指数δ_T、重组分权重系数δ_Oil。将甲烷数据数列、出口电导率数据数列和出口温度数据数列近似看作符合正态分布，显然，对于待求参数甲烷变化指数δ_C1、出口电导率变化指数δ_S、出口温度变化指数δ_T、出口钻井液含烃较重组分权重系数δ_Oil，有如下计算方式实现：

其中，TG为色谱检测到的全烃组分值，单位％，C₁为色谱检测到的甲烷组分值，单位％，C₂为色谱检测到的乙烷组分值，单位％，WOB为钻压，单位kN·m，ROP为钻时，单位min/m，μ_C1为待求深度点所在井段对应甲烷组分平均值，μ_s为待求深度点所在井段对应出口电导率平均值，μ_T为待求深度点所在井段对应出口温度平均值，μ_TG-C1-C2为待求深度点所在井段对应出口C₁、C₂之外烃重组分平均值，△_C1为待求深度点所在井段对应出口电导率标准差，δ_s为待求深度点所在井段对应出口电导率标准差，△_T为待求深度点所在井段对应出口温度标准差，△_TG-C1-C2为待求深度点所在井段对应出口C₁、C₂之外烃重组分标准差。

设若

则含水指数W计算方式如下：

设若

则含气指数G计算方式如下：

设若

则含油指数O计算方式如下：

3)集成三角图版构建和数字化

如图2所示，建立如图正三角形，边长L代表油(Oil)、气(Gas)、水(Water)组分占钻井液体积理论最大百分含量，即L最大值L_max＝1，将边长等分为三段，即有

L＝L₁-L₂+L₃

本发明设定有一正圆O₁与此正三角形各边三等分点均相交，易知其半径

更有与该圆同心另一圆形O₂，其半径R₂

如图3，刻度为顺时针(0,1)，即按照箭头方向，Oil边自下而上，为零点和1点，Water边自上而下，分别为零点和1点，Gas边自右及左，依次为零点和1点。

如图4-6所示，据W系数、G系数和O系数计算结果，分别在三角图版各对应边分别标注，分别作过三点与其所在轴零刻度衔接的邻边之平行线，三线必然两两相交成一正三角形，其重心(中心)落在图版上对应的区间则可判定该异常井段的储层发育情况及流体性质。

4)坐标变换

将二维笛卡尔坐标系转换到本发明设定的三角坐标系，设含水指数为L_w，含油指数为L_o，含气指数L_g，

P_o(x_o,y_o):

P_w(x_w，x_w)：

P_g(x_g，y_g)：

x_g＝1-L_g

y_g＝0

如图7所示，由于平行于L_g边的辅助线必然斜率为0，故P_g和P_w所在辅助线的交点为三角形顶点，而后三角形中心P(x，y)的坐标分别为：

由此可以方便将各深度点对应P坐标投影到三角图版中，也可以根据图版特征直接判断P点所在的流体识别区，将待解释井段的对应各深度点对应P点值投影到三角图版中即可根据散点的展布判断出该段可能钻遇某种流体，若散点分布在两个流体识别区则表明可能为两种流体并存。

5)基于含水指数W、含油指数O以及含气指数G的储层含流体概率图道建立

为了提高识别的准确性，还可以进一步做如下处理：

应用含水指数、含气指数和含油指数建立含流体性质相对概率图道，进一步提高识别的准确性。

该概率图道既可用多点井深组合的深度段进行总体识别，也可逐深度点识别。现分别声明如下：

A、单点概率图道的建立

对于单深度点的情形，设Wi为某深度点对应含水指数，Oi为某深度点含油指数，Gi为某深度点含气指数，对于任意序号为i的深度点，则有：

其中，i＝1,2,3,…，n，i∈Z

B、组合深度点的某井段概率图道的建立

设Wi为某深度点对应含水指数，Oi为某深度点含油指数，Gi为某深度点含气指数，设某井段内有m个深度点，即该井段长m米，则对应的概率图道建立方法如下：

DW＝∑DWi；DO＝∑DOi；DG＝∑DGi；

其中，

i＝1,2,3,…，m，i∈Z

DW：对应井段相对含水可能性，单位％

DO：对应井段相对含油可能性，单位％

DG：对应井段相对含气可能性，单位％

7)利用概率图道进行储层含流体性质识别的判别模式，如图8所示：

设Wi为某深度点对应含水指数，Oi为某深度点含油指数，Gi为某深度点含气指数，则对应的概率图道建立方法如下，对于单深度点的情形，直接比较DOi、DWi和DGi的大小关系，且按照以下原则进行分类：

当DOi或者DWi或者DGi其一小于16％时，则不参与最终流体性质判别；

当DOi或者DWi或者DGi其一大于66％时，判断对应井深赋存对应流体类型的可能性最大，且以其为主进行定名；

(1)含油层识别方面：

如DOi大于66％，DWi小于16％，DGi大于25％时，解释该层为油气同层；

如DOi大于66％，DWi大于25％，DGi小于16％时，解释该层为油水同层；

如DOi大于66％且小于75％，同时DWi、DGi均小于16％时，解释该层为差油层；

如DOi大于75％时，且DWi、DGi均小于16％时，解释该层为油层；

对于油层考虑到原油的物理性质特性，可适当下调66％基准至40％～50％。

(2)含气层识别方面：

如DGi大于66％，DWi小于16％，DOi大于25％时，解释该层为油气同层；

如DGi大于66％，DWi大于25％，DOi小于16％时，解释该层为气水同层；特殊地，当DGi与DWi均未大于66％且数值接近，同时DOi小于16％的情形，解释为气水同层亦可；

如DGi大于66％且小于75％，同时DWi、DOi均小于16％时，解释该层为差气层；

如DGi大于75％时，且DWi、DOi均小于16％时，解释该层为气层；

对于气层，建议上浮66％基准至75％～80％。

(3)含水层识别方面：

如DWi大于66％，DGi小于16％，DOi大于25％时，解释该层为油水同层；

如DWi大于66％，DGi大于25％，DOi小于16％时，解释该层为含气水层；特殊地，当DGi与DWi均未大于66％且数值接近，同时DOi小于16％的情形，解释为气水同层；考虑气水相渗中的滑脱效应等现象，在DWi、DGi同时参与判定时，适当上调DGi基准10％～15％。

如DWi大于85％时，且DGi、DOi均小于16％时，解释该层为水层。

Claims (6)

1.一种基于钻井液录井参数气测值的录井解释方法，其特征在于，包括如下步骤：

a、对钻井液录井采集到的气测值进行预处理；

b、将预处理后的气测值、钻时、电导率和温度计算变化指数；

c、通过变化指数计算含水指数W、含油指数O和含气指数G；

d、建立以含水指数W、含油指数O和含气指数G相对比值为三边的三角图版；

e、利用建立的三角图版进行储层含流体性质识别；

所述步骤a中，通过校正气测值，剔除干扰数据；校正气测值为：设钻杆单根长度区间为[8,10]，则T_单根∈[8,10]，单根频率f>0.125，以此为截止频率f_t，

a_i＝f_t/2πn_i＝0.1989,i＝1,2,3…j,i∈N

则用迭代公式计算

Y_i＝a_iX_i+(1-a_i)Y_i-1

其中，Y_i----------当前井深i点校正值；

X_i---------当前井深i点预处理前读值；

Y_i-1----------当前井深的i-1点预处理后读值；

a_i-----------钻具系数,根据钻具长度取值，自然数组，无量纲；

n_i-----------钻具长度，m；

其中，a_i与单根长度有关。

2.根据权利要求1所述的基于钻井液录井参数气测值的录井解释方法，其特征在于：所述气测值包括全烃组分值TG，单位％或ppm；甲烷组分值C₁，单位％或ppm；乙烷组分值C₂，单位％或ppm。

3.根据权利要求2所述的基于钻井液录井参数气测值的录井解释方法，其特征在于：所述步骤b中，通过全烃组分值T、甲烷组分值C₁、乙烷组分值C₂、出口电导率值S、出口温度值T，分别计算甲烷变化指数δ_C1、出口电导率变化指数δ_s、出口温度变化指数δ_T、重组分权重系数δ_Oil，再通过甲烷变化指数δ_C1、出口电导率变化指数δ_S、出口温度变化指数δ_T、重组分权重系数δ_Oil计算含水指数W、含油指数O和含气指数G。

4.根据权利要求3所述的基于钻井液录井参数气测值的录井解释方法，其特征在于：所述甲烷变化指数δ_C1、出口电导率变化指数δ_S、出口温度变化指数δ_T、重组分权重系数δ_Oil分别通过下述公式计算得出：

其中，TG为色谱检测到的全烃组分值，单位％，C₁为色谱检测到的甲烷组分值，单位％，C₂为色谱检测到的乙烷组分值，单位％，WOB为钻压，单位kN·m，ROP为钻时，单位min/m，

为待求深度点所在井段对应甲烷组分平均值，μ_s为待求深度点所在井段对应出口电导率平均值，μ_T为待求深度点所在井段对应出口温度平均值，

为待求深度点所在井段对应出口C₁、C₂之外烃重组分平均值，

为待求深度点所在井段对应出口电导率标准差，Δ_S为待求深度点所在井段对应出口电导率标准差，Δ_T为待求深度点所在井段对应出口温度标准差，

为待求深度点所在井段对应出口C₁、C₂之外烃重组分标准差。

5.根据权利要求4所述的基于钻井液录井参数气测值的录井解释方法，其特征在于：所述步骤c中，含水指数W、含油指数O和含气指数G分别通过下述公式计算得出：

设

则含水指数W计算方式如下：

设

则含气指数G计算方式如下：

设

则含油指数O计算方式如下：

6.根据权利要求5所述的基于钻井液录井参数气测值的录井解释方法，其特征在于：所述步骤d具体包括如下步骤：

d1、建立图版框

建立三角图版，边长L代表油Oil、气Gas、水Water组分占钻井液体积理论最大百分含量，即L最大值L_max＝1，将边长等分为三段，即有

L＝L₁+L₂+L₃

设定有一正圆O₁与此正三角形各边三等分点均相交，其半径：

有与该圆同心另一圆形O₂，其半径R₂

刻度为顺时针设置(0,1)，即按照顺时针方向，Oil边自下而上，为零点和1点，Water边自上而下，分别为零点和1点，Gas边自右及左，依次为零点和1点；

根据W系数、G系数和O系数结果，分别在三角图版各对应边分别标注，分别作过三点与其所在轴零刻度衔接的邻边之平行线，三线必然两两相交，围成一个正三角形，其重心落在图版上对应的区间则判定井段的储层发育情况及流体性质；

d2、IGIC坐标变换计算

将二维笛卡尔坐标系转换到三角图版设定的三角坐标系，设含水指数为L_w，含油指数为L_o，含气指数L_g，

P_o(x_o,y_o):

P_w(x_w,y_w):

P_g(x_g,y_g):

x_g＝1-L_g

y_g＝0

由于平行于L_g边的辅助线必然斜率为0，故P_g和P_w所在辅助线的交点为三角形顶点，而后三角形中心P(x,y)的坐标分别为：

x：

y：

由此各深度点对应P坐标投影到三角解释集成图版中，或根据图版特征直接判断P点所在的流体识别区，将待解释井段的对应各深度点对应P点值投影到三角图版中根据散点的展布判断出该段可能钻遇某种流体，若散点分布在两个流体识别区则表明可能为两种流体并存。

Images