一种地质录井中随钻岩石风化势能指数确定方法及装置与流程

本发明涉及一种地质录井中随钻岩石风化势能指数确定方法及装置，属于岩石风化判别技术领域。

背景技术：

风化是指岩石在物理、化学、生物作用下发生破碎、不稳定矿物分解的过程，在石油地质研究中具有重要意义。风化过程中，随着不稳定矿物分解，岩石元素不断变化，首先是活跃的碱金属元素(na、k)和碱土金属元素(ca、mg)被溶解淋滤，继而是含p、s、si元素矿物分解，随着风化进程，由于其它组分的相对流失，造成较稳定元素fe、al、ti相对富集。因此，许多学者根据元素的流失和迁移特性，利用氧化物含量计算用以衡量风化程度的指数。

目前，石油钻井中风化指数主要是根据岩石矿物中氧化物摩尔数相对变化关系来表达，如威格特残积指数v＝(al2o3+k2o)/(mgo+cao+na2o)。为进一步精确刻画岩石风化程度，可根据母源矿物差异，优选反映风化程度的敏感矿物来风化指数，如母源矿物以硅酸岩为主的风化指数为帕克风化指数wip＝100*(2na2o/0.35+2k2o/0.25+mgo/0.9+cao*/0.7)；母源矿物以花岗岩为主的风化指数(weatheringindexofgranite)为wig＝100*(na2o+k2o+cao-10/3*p2o5)/(al2o3+fe2o3+tio2)。

例如，公布号为cn105717149a的中国发明专利申请文件，该文件公开了一种玄武岩风化程度的判别方法，该方法首先采集玄武岩样本，然后检测出玄武岩样本中八种氧化物的含量，并根据各氧化物的含量计算玄武岩样本的化学指数bwi的值，bwi＝(fe2o3+al2o3)/(sio2+k2o+na2o+cao+mgo+feo)，最后根据得到的化学指数来判断玄武岩石的风化程度。

上述方法均利用氧化物含量计算风化指数方法，需要在钻井岩心中采取岩石样品，并在实验室中化验分析获得氧化物含量。因钻井取心井段少，且不连续，所以获得的风化指数缺乏系统性；且因分析周期长，难以满足生产需要。

技术实现要素：

本发明的目的是提供一种地质录井中随钻岩石风化势能指数确定方法，以解决目前采用实验室中获取的岩石氧化物含量来表征岩石风化势能指数导致获取的风化势能指数不及时，不连续的问题；同时，本发明还提供了一种地质录井中随钻岩石风化势能指数确定装置。

本发明为解决上述技术问题而提供一种地质录井中随钻岩石风化势能指数确定方法，该确定方法包括以下步骤：

1)在研究区石油钻井现场对不同钻井深度中岩心、岩屑中的化学元素含量进行测定，得到不同深度化学元素含量；

2)将得到的不同深度化学元素含量带入元素风化势能指数计算模型，确定不同深度的元素风化势能指数，所述的元素风化势能指数包括帕克风化势能指数fwip和/或花岗岩风化势能指数fwig。

本发明根据石油钻井现场实时获取的岩石元素含量计算元素风化势能指数，得到元素风化势能指数时效性高，且得到的元素风化势能指数能够随井深连续变化，应用更广泛，克服了实验室得到的风化势能指数时效性低、不连续的问题。

进一步地，本发明给出了元素风化势能指数的计算模型，每种元素增加了相应的区域常数，计算的风化势能指数更能准确反映研究区岩石实际风化程度，母源矿物以硅酸盐岩矿物为主的帕克风化势能指数fwip计算模型为：

fwip＝100·(0.1242c1·na+0.1026c2·k+0.0463c3·mg+0.0357c4·ca-0.0768c5·p)

母源矿物以花岗岩为主的花岗岩风化势能指数fwig计算模型为：

fwig＝100·(0.0217d1·na+0.0128d2·k+0.0250d3·ca-0.0538d4·p)/(0.0185d5·al+0.0089d6·fe+0.0204d6·ti)其中al、k、mg、ca、na、p、fe和ti分别为化学元素铝、钾、镁、钙、钠、磷、铁和钛的质量百分含量，c1、c2、c3、c4、c5、d1、d2、d3、d4、d5和d6均为与母源矿物相关的区域常数，理想模型下均为1。

进一步地，本发明还给出了计算模型中区域常数的确定过程，使计算模型更能贴近实际，提高元素风化势能指数的准确率，所述区域常数的确定过程如下：

a.在研究区筛选风化程度不同层位的岩心为采样对象，获取各采样对象的化学元素含量和氧化矿物含量；

b.根据得到的化学元素含量确定理想模型下的元素风化势能指数，所述的元素风化势能指数包括帕克风化势能指数fwip和/或花岗岩风化势能指数fwig，根据氧化矿物含量确定氧化物风化指数；

c.将得到的各采样对象的氧化物风化指数与对应采样对象的元素风化势能指数进行拟合，对理想模型中的区域常数进行调整，使元素风化势能指数与氧化物风化指数两绝对相关系数达到设定值以上。

进一步地，本发明为实现对岩石化学元素含量的测定，给出了化学元素含量的具体测量方式，所述化学元素含量采用x射线元素录井仪测量得到。

进一步地，为了确定元素风化势能指数计算模型中的区域常数，本发明给出了拟合的方式，所述步骤b中的氧化物风化指数由样品岩心氧化矿物摩尔数百分含量计算得到，所述的氧化物风化指数包括帕克风化指数wip和/或花岗岩风化指数wig，所述步骤c在拟合时是将帕克风化指数wip与帕克风化势能指数fwip进行拟合，将花岗岩风化指数wig与花岗岩风化势能指数fwig进行拟合。

进一步地，所述步骤c中的设定值为0.9。

同时本发明还提供了一种地质录井中随钻岩石风化势能指数确定装置，该确定装置包括存储器和处理器，以及存储在所述存储器上并在所述处理器上运行的计算机程序，所述处理器与所述存储器相耦合，所述处理器执行所述计算机程序时实现以下步骤：

1)在研究区石油钻井现场对不同钻井深度中岩心、岩屑中的化学元素含量进行测定，得到不同深度化学元素含量；

2)将得到的不同深度化学元素含量带入元素风化势能指数计算模型，确定不同深度的元素风化势能指数，所述的元素风化势能指数包括帕克风化势能指数fwip和/或花岗岩风化势能指数fwig。

进一步地，母源矿物以硅酸盐岩矿物为主的帕克风化势能指数fwip计算模型为：

fwip＝100·(0.1242c1·na+0.1026c2·k+0.0463c3·mg+0.0357c4·ca-0.0768c5·p)

母源矿物以花岗岩为主的花岗岩风化势能指数fwig计算模型为：

fwig＝100·(0.0217d1·na+0.0128d2·k+0.0250d3·ca-0.0538d4·p)/(0.0185d5·al+0.0089d6·fe+0.0204d6·ti)其中al、k、mg、ca、na、p、fe和ti分别为化学元素铝、钾、镁、钙、钠、磷、铁和钛的质量百分含量，c1、c2、c3、c4、c5、d1、d2、d3、d4、d5和d6均为与母源矿物相关的区域常数，理想模型下均为1。

进一步地，所述区域常数的确定过程如下：

a.在研究区筛选风化程度不同层位的岩心为采样对象，获取各采样对象的化学元素含量和氧化矿物含量；

b.根据得到的化学元素含量确定理想模型下的元素风化势能指数，所述的元素风化势能指数包括帕克风化势能指数fwip和/或花岗岩风化势能指数fwig，根据氧化矿物含量确定氧化物风化指数；

c.将得到的各采样对象的氧化物风化指数与对应采样对象的元素风化势能指数进行拟合，对理想模型中的区域常数进行调整，使元素风化势能指数与氧化物风化指数两绝对相关系数达到设定值以上。

进一步地，所述化学元素含量采用x射线元素录井仪测量得到。

进一步地，所述步骤b中的氧化物风化指数由样品岩心氧化矿物摩尔数百分含量计算得到，所述的氧化物风化指数包括帕克风化指数wip和/或花岗岩风化指数wig，所述步骤c在拟合时是将帕克风化指数wip与帕克风化势能指数fwip进行拟合，将花岗岩风化指数wig与花岗岩风化势能指数fwig进行拟合。

进一步地，所述步骤c中的设定值为0.9。

附图说明

图1是本发明华北地区某油田不同风化程度层位样品元素风化势能指数fwip与帕克风化指数wip拟合图；

图2是本发明华北地区某油田不同风化程度层位样品元素风化势能指数fwig与花岗岩风化指数wig拟合图；

图3是本发明华北地区某油田某井x射元素录井风化势能指数综合分析图。

具体实施方式

下面结合附图对本发明的具体实施方式进行详细说明。

本发明地质录井中随钻岩石风化势能指数确定方法的实施例

本发明首先在研究区石油钻井现场对不同钻井深度中岩心、岩屑中的化学元素含量进行测定，得到不同深度化学元素含量；然后将得到的不同深度化学元素含量带入元素风化势能指数计算模型，确定不同深度的元素风化势能指数，其中元素风化势能指数包括帕克风化势能指数fwip和花岗岩风化势能指数fwig。下面对各步骤的具体实现方式进行详细说明。

1.在研究区筛选风化程度不同层位的岩心为采样对象，获取各采样对象的化学元素含量和氧化矿物含量。

在区域地质研究的基础上，筛选风化程度不同层位具有代表性的岩心为采样对象。采用x射线元素录井仪(或具有相同功能的x射线荧光分析仪)对采样对象的化学元素进行分析，得到各化学元素含量，其中测定的化学元素包括但不限于以下元素：硅、铝、铁、钛、钠、钾、钙、镁、磷；利用实验室分析各采样对象的氧化矿物摩尔数百分含量。

2.确定采样对象的氧化物风化指数和理想模型下元素风化势能指数。

本发明的元素风化势能指数包括帕克风化势能指数fwip和花岗岩风化势能指数fwig。本实施例选取包括帕克风化势能指数fwip和花岗岩风化势能指数fwig进行风化势能指数计算。

根据研究区母源类型优选出反映风化程度的敏感元素种类，利用敏感元素组合优化风化势能指数计算模型，母源矿物以硅酸盐岩矿物为主的帕克风化势能指数fwip计算模型为：

fwip＝100·(0.1242c1·na+0.1026c2·k+0.0463c3·mg+0.0357c4·ca-0.0768c5·p)

其中k、mg、ca、na、p分别为化学元素钾、镁、钙、钠、磷的质量百分含量，c1、c2、c3、c4、c5均为与母源矿物相关的区域常数，理想模型下均为1。fwip数值越大代表越容易被风化，即岩石风化程度越低，0代表岩石完全风化。

根据研究区母源类型优选出反映风化程度的敏感元素种类，利用敏感元素组合优化风化指数计算模型，母源矿物以花岗岩为主的花岗岩风化势能指数fwig计算模型为：

fwig＝100·(0.0217d1·na+0.0128d2·k+0.0250d3·ca-0.0538d4·p)/(0.0185d5·al+0.0089d6·fe+0.0204d6·ti)

其中al、k、ca、na、p、fe和ti分别为化学元素铝、钾、钙、钠、磷、铁和钛的质量百分含量；d1、d2、d3、d4、d5和d6均为与母源矿物相关的区域常数，理想模型下均为1。fwig数值越大代表越容易被风化，即岩石风化程度越低，0代表岩石完全风化。

通过上述公式即可计算出理想模型下各采样对象的风化势能指数。

利用实验室分析的氧化矿物摩尔数百分含量，根据帕克风化指数wip和花岗岩风化指数wig公式，得到相应的氧化物风化指数，具体计算公式已在背景技术中介绍，这里不再详述。

3.进行数据拟合，调整元素风化势能指数中的区域常数。

将实验室获取的各采样对象的帕克风化指数wip与帕克风化势能指数fwip进行拟合，根据统计分析规律对帕克风化势能指数计算模型中的区域常数微调，使帕克风化势能指数与实验室氧化矿物计算的帕克风化指数wip两绝对相关系数达到设定值(本实施例为0.9)以上，此时，帕克风化势能指数fwip计算模型中的区域常数即为针对该研究区的区域常数，将得到的区域常数带入帕克风化势能指数fwip计算模型，即可得到优化的帕克风化势能指数fwip计算模型。

同理，将实验室获取的各采样对象的花岗岩风化指数wig与花岗岩风化势能指数fwig进行拟合，根据统计分析规律对风化势能指数计算模型中的区域常数微调，使花岗岩风化势能指数fwig与实验室氧化矿物计算的花岗岩风化指数wig两绝对相关系数达到设定值(本实施例为0.9)以上，此时，花岗岩风化势能指数fwig计算模型中的区域常数即为针对该研究区的区域常数，将得到的区域常数带入花岗岩风化势能指数fwig计算模型，即可得到优化的花岗岩风化势能指数fwig计算模型。

4.根据优化后的计算模型计算现场元素风化指数。

采用x射线元素录井仪(或具有相同功能的x射线荧光分析仪)在石油钻井现场随钻井深度对岩心、岩屑中的化学元素进行连续分析，得出随深度变化的元素含量数据，其中化学元素包括但不限于以下元素：硅、铝、铁、钛、钠、钾、钙、镁、磷。将现场获取的岩石元素含量带入到优化后的帕克风化势能指数fwip计算模型和花岗岩风化势能指数fwig计算模型计算相应的元素风化势能指数，并绘制出风化势能指数曲线，利用得到的风化势能指数值对岩石风化程度定量评价，根据风化势能指数值和曲线变化对地层进行划分对比，确定风化壳风化程度对储层物性和油气富集规律研究。

下面以华北地区某油田为例对本发明的效果进行验证。

1.首先在该区域地质研究的基础上，筛选出风化程度不同层位具有代表性的岩心为采样对象，取1吋小岩心柱从中间平分，编写相同编号，一半样品做x射线元素分析，获得岩生石元素质量百分比数据，如表1所示；另一半送实验室做氧化矿物含量分析，获得岩石氧化矿物摩尔数百分含量数据，如表2所示。

2.利用x射线元素录井仪分析各采样对象的元素质量百分含量，根据理想模型下元帕克风化势能指数fwip和花岗岩风化势能指数fwig计算模型(即计算模型中区域常数为1)，初步计算岩石风化指数；利用实验室分析的氧化矿物摩尔数百分含量，根据帕克风化指数wip和花岗岩风化指数wig公式，计算岩石风化程度。

3.将实验室获取的各采样对象的帕克风化指数wip与帕克风化势能指数fwip进行拟合，将实验室获取的各采样对象的风化指数wig与花岗岩风化势能指数fwig进行拟合，根据统计分析规律对风化势能指数计算模型中的区域常数微调，使元素风化势能指数与实验室氧化矿物计算的风化指数两绝对相关系达0.9以上(见图1、图2所示)，此时，元素风化势能指数计算模型中的区域常数即为针对研究区的区域常数，根据得到区域常数确定优化后的元素风化势能指数计算模型。

4.采用x射元素分析仪在石油钻井现场随钻井深度对岩心、岩屑中的化学元素含量进行连续分析，得出随深度变化的x射线元素含量；利用x射线元素含量数据，根据优化后的元素风化指数计算模型计算岩石帕克风化势能指数fwip和花岗岩风化势能指数fwig，计算结果见表3。

5.根据岩石帕克风化势能指数fwip和花岗岩风化势能指数fwig数据、岩心段实验室风化指数数据以及常规地质录井数据绘制风化指数综合图。

图3为华北地区某油田某钻井根据x射元素录井分析的元素含量计算的帕克风化势能指数fwip和花岗岩风化势能指数fwig与岩心段实验室分析计算的帕克风化指数wip和花岗岩风化指数wig对比分析综合图。

从图3中看出，实验室分析数据计算的风化指数与石油钻井现场x射线元素分析数据计算的风化指数吻合性好。四个取心段(50205028m，5080-5084，5115-5124m，5180-5184m)实验室计算的岩石风化程度(帕克风化指数wip和花岗岩风化指数wig)分布在元素风化势能指数曲线(元素风化势能指数fwip和元素风化势能指数fwig)两侧，证明地质录井中随钻岩石风化势能指数计算方法获得的岩石风化指数准确可靠。

同时，从图3中还可以看出，风化壳上覆地层徐庄组属于强风化地层，元素风化势能指数fwip和元素风化势能指数fwig大幅负异常，说明风化程度高；其下伏海相地层毛庄组属还原环境沉积，元素风化势能指数fwip和元素风化势能指数fwig大幅正异常异常，说明风化程度低；上部陆相地层沙河街组元素风化势能指数fwip和元素风化势能指数fwig居中，说明元素风化势能指数fwip和元素风化势能指数fwig真实反映了地层岩石风化程度。

因此，本发明能够在石油钻井现场实时计算风化指数，时效性高，且能获得连续景深的风化指数，应用更广泛，克服了实验室得到的风化指数不连续的问题；此外，每种元素增加了相应的区域常数，计算的风化指数更能准确反映研究区岩石实际风化程度。

本发明地质录井中随钻岩石风化势能指数确定装置的实施例

上述方法可以作为一种计算机程序，存储在地质录井中随钻岩石风化势能指数确定装置中的存储器中并可在地录井中随钻岩石风化势能指数确定装置中的处理器上运行。该装置中的处理器可以采用单片机、dsp、plc或mcu等实现，存储器可以采用ram存储器、闪存、rom存储器、eprom存储器、eeprom存储器、寄存器、硬盘、移动磁盘、cd-rom或者本领域已知的任何其他形式的存储介质，可以将该存储介质耦接至处理器，使处理器能够从该存储介质读取信息，或者该存储介质可以是处理器的组成部分。该装置中所执行的步骤方法已在方法的实施例中进行了详细说明，这里不再赘述。

上述各实施例仅用于说明本发明，其中方法的各实施步骤等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之内。

表1：华北地区某油田x射线元素录井分析与风化势能计算统计表

续表1：华北地区某油田x射线元素录井分析与风化指数计算统计表

续表1：华北地区某油田x射线元素录井分析与风化指数计算统计表

表2：华北地区某油田岩心氧化矿物分析与风化指数计算统计表

续表2：华北地区某油田岩心氧化矿物分析与风化指数计算统计表

续表2：华北地区某油田岩心氧化矿物分析与风化指数计算统计表

表3：华北地区某油田某井x射线元素录井分析与风化势能计算表

续表3：华北地区某油田某井x射线元素录井分析与风化势能计算表

续表3：华北地区某油田某井x射线元素录井分析与风化势能计算表

续表3：华北地区某油田某井x射线元素录井分析与风化势能计算表

续表3：华北地区某油田某井x射线元素录井分析与风化势能计算表

续表3：华北地区某油田某井x射线元素录井分析与风化势能计算表

续表3：华北地区某油田某井x射线元素录井分析与风化势能计算表