本发明专利技术提供一种火成岩矿物成分确定方法、装置、计算机设备和存储介质，根据原始元素测井/XRF资料火成岩地层元素相对含量，利用纯矿物建立地层矿物质量与元素含量之间存在的系数矩阵，最终利用优化算法计算火成岩地层矿物重量含量，能评价出14种火成岩成分。反演过程实现了模块化操作，为适应于国内火成岩快速、高效开发，在形成自有知识产权的元素测井解释软件方面起到了积极作用。软件方面起到了积极作用。软件方面起到了积极作用。
【技术实现步骤摘要】
火成岩矿物成分确定方法、装置、计算机设备和存储介质
[0001]本专利技术涉及火成岩地层矿物类型及含量识别
，特别涉及一种火成岩矿物成分确定方法、装置、计算机设备和存储介质。
技术介绍
[0002]金云智，高楚桥，等(2018年)。在地层组分分析基础上建立火成岩地层元素测井资料与其各岩性间的关系模型；运用最优化算法，从地层元素测井Al、Ca、Fe、H、Si、S、Ti等地层元素含量出发，定量计算火成岩岩性组分的含量。
[0003]中国专利公开文本CN107194401 A公开了在火成岩QAPF等标准灰度图图版的基础上，通过火成岩矿物成分投点、投点位置灰度值获取、基于哈希表集合快速获取火成岩名称等环节，快速实现火成岩的自动化分类命名。中国专利公开文本CN107678072A公开了综合利用磁力数据、地震数据、岩心数据、测井数据预测火成岩储层的方法既能合理利用磁力识别火成岩，又能利用三维地震格架弥补其空间分辨率的不足，并且在钻井数据的约束下，使火成岩储层的预测效果达到最佳。中国专利文本CN101520518 A公开了一种利用重磁电异常的组合特征识别火成岩岩性的方法，采集记录岩芯物理特性；用k
‑
均值聚类对岩石的密度、磁化率、电阻率进行聚类分析，获得岩芯聚类编码；对重磁勘探结果进行重磁异常剥离获得异常，采用同样的聚类进行三维空间异常编码；将岩石物性聚类编码与重磁电异常编码相互结合反映探区每一处火成岩的岩性。
[0004]求取矿物质量含量解释方法归结为两种：一种是斯伦贝谢公司为代表的经验解释方法，另一种是哈里伯顿公司为代表的优化解释方法。斯伦贝谢公司解释方法是由经验关系式求取火成岩地层元素相对含量和矿物质量含量；哈里伯顿公司解释方法是直接由地层元素含量利用优化方法反演得到地层矿物质量含量。斯伦贝谢公司为代表的解释方法依赖大量的地区实验数据，对应一些新地层，可能需要重新建立模型，而且建立的模型没有公开，适用性需要进一步验证；哈里伯顿公司为代表的解释方法适应性较强，但存在反演矿物的不确定性。
技术实现思路
[0005]基于此，有必要针对上述技术问题，提供一种火成岩矿物成分确定方法、装置、计算机设备和存储介质。
[0006]一种火成岩矿物成分确定方法，包括：
[0007]测量火成岩所包含的纯矿物质量，根据所述纯矿物质量计算得到所述纯矿物的摩尔质量及相对分子量，计算所述纯矿物中的酸碱根的相对分子量；
[0008]根据所述纯矿物的摩尔质量及相对分子量确定所述纯矿物的分子式模型；
[0009]根据所述纯矿物的分子式模型，确定所述纯矿物中的至少一种主元素；
[0010]根据所述纯矿物的相对分子量和酸碱根的相对分子量，计算获得各所述主元素的相对分子量；
[0011]建立各所述主元素的相对分子量与所述纯矿物的分子式模型中的分配系数的关系式；
[0012]根据各所述关系式，采用优化算法计算获得所述分配系数的值；
[0013]检测所述分配系数的值是否符合预设条件；
[0014]当所述分配系数的值符合预设条件时，计算获得各所述主元素在所述纯矿物中的系数；
[0015]根据各所述主元素在所述纯矿物中的系数，确定所述纯矿物的矿物系数矩阵；
[0016]根据所述纯矿物的矿物系数矩阵，采用最优化方法计算所述纯矿物的含量；
[0017]根据所述纯矿物的含量，确定所述火成岩所含矿物成分。
[0018]在其中一个实施例中，所述根据所述纯矿物的矿物系数矩阵，计算所述纯矿物的矿物成分含量的步骤包括：
[0019]根据所述纯矿物的矿物系数矩阵，确定火成岩地层中所述主元素与纯矿物的关系式及误差目标函数；
[0020]根据所述主元素与纯矿物的关系式及误差目标函数，采用最优化方法计算所述纯矿物的含量。
[0021]在其中一个实施例中，所述检测所述分配系数的值是否符合预设条件的步骤之后还包括：
[0022]当所述分配系数的值不符合预设条件时，返回执行所述根据所述纯矿物的摩尔质量及相对分子量确定所述纯矿物的分子式模型的步骤。
[0023]在其中一个实施例中，所述根据各所述关系式，采用优化算法计算获得所述分配系数的值的步骤包括：
[0024]根据各所述关系式，采用最小二乘法计算获得所述分配系数的值。
[0025]在其中一个实施例中，所述当所述分配系数的值符合预设条件时，计算获得所述主元素在纯矿物中的系数的步骤包括：
[0026]当所述分配系数的值符合预设条件时，根据所述主元素的相对分子量与所述纯矿物的相对分子量之比，计算获得所述主元素在纯矿物中的系数。
[0027]在其中一个实施例中，所述建立各所述主元素的相对分子量与所述纯矿物的分子式模型中的分配系数的关系式的步骤中，
[0028]各所述主元素的相对分子量与所述纯矿物的分子式模型中的分配系数的关系式为：
[0029]Y＝ΣA
i
X
i
[0030]其中，Y为所述主元素的相对分子量，A分子式模型，为X为主元素在所述分子式模型中的分配系数。
[0031]在其中一个实施例中，所述根据所述纯矿物的摩尔质量及相对分子量确定所述纯矿物的分子式模型的步骤中，
[0032]所述纯矿物的分子式模型为A(X)，其中，A为所述纯矿物中包含的主元素，X为主元素的分配系数。
[0033]一种火成岩矿物成分确定装置，包括：
[0034]纯矿物分子量计算模块，用于测量火成岩所包含的纯矿物质量，根据所述纯矿物
质量计算得到所述纯矿物的摩尔质量及相对分子量，计算所述纯矿物中的酸碱根的相对分子量；
[0035]分子式模型确定模块，用于根据所述纯矿物的摩尔质量及相对分子量确定所述纯矿物的分子式模型；
[0036]主元素确定模块，用于根据所述纯矿物的分子式模型，确定所述纯矿物中的至少一种主元素；
[0037]相对分子量计算模块，用于根据所述纯矿物的相对分子量和酸碱根的相对分子量，计算获得各所述主元素的相对分子量；
[0038]分配关系式建立模块，用于建立各所述主元素的相对分子量与所述纯矿物的分子式模型中的分配系数的关系式；
[0039]分配系数计算模块，用于根据各所述关系式，采用优化算法计算获得所述分配系数的值；
[0040]预设条件符合检测模块，用于检测所述分配系数的值是否符合预设条件；
[0041]元素系数计算模块，用于当所述分配系数的值符合预设条件时，计算获得各所述主元素在所述纯矿物中的系数；
[0042]系数矩阵确定模块，用于根据各所述主元素在所述纯矿物中的系数，确定所述纯矿物的矿物系数矩阵；
[0043]矿物成分含量计算模块，用于根据所述纯矿物的矿物系数矩阵，采用最优化方法计算所述纯矿物的含量；
[0044]矿物成分确定模块，用于根据所述纯矿物的含量，确定所述火成岩所含矿物成分。
[0045]一本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种火成岩矿物成分确定方法，其特征在于，包括：测量火成岩所包含的纯矿物质量，根据所述纯矿物质量计算得到所述纯矿物的摩尔质量及相对分子量，计算所述纯矿物中的酸碱根的相对分子量；根据所述纯矿物的摩尔质量及相对分子量确定所述纯矿物的分子式模型；根据所述纯矿物的分子式模型，确定所述纯矿物中的至少一种主元素；根据所述纯矿物的相对分子量和酸碱根的相对分子量，计算获得各所述主元素的相对分子量；建立各所述主元素的相对分子量与所述纯矿物的分子式模型中的分配系数的关系式；根据各所述关系式，采用优化算法计算获得所述分配系数的值；检测所述分配系数的值是否符合预设条件；当所述分配系数的值符合预设条件时，计算获得各所述主元素在所述纯矿物中的系数；根据各所述主元素在所述纯矿物中的系数，确定所述纯矿物的矿物系数矩阵；根据所述纯矿物的矿物系数矩阵，采用最优化方法计算所述纯矿物的含量；根据所述纯矿物的含量，确定所述火成岩所含矿物成分。2.根据权利要求1所述的方法，其特征在于，所述根据所述纯矿物的矿物系数矩阵，计算所述纯矿物的矿物成分含量的步骤包括：根据所述纯矿物的矿物系数矩阵，确定火成岩地层中所述主元素与纯矿物的关系式及误差目标函数；根据所述主元素与纯矿物的关系式及误差目标函数，采用最优化方法计算所述纯矿物的含量。3.根据权利要求1所述的方法，其特征在于，所述检测所述分配系数的值是否符合预设条件的步骤之后还包括：当所述分配系数的值不符合预设条件时，返回执行所述根据所述纯矿物的摩尔质量及相对分子量确定所述纯矿物的分子式模型的步骤。4.根据权利要求1所述的方法，其特征在于，所述根据各所述关系式，采用优化算法计算获得所述分配系数的值的步骤包括：根据各所述关系式，采用最小二乘法计算获得所述分配系数的值。5.根据权利要求1所述的方法，其特征在于，所述当所述分配系数的值符合预设条件时，计算获得所述所有主元素在纯矿物中的系数的步骤包括：当所述分配系数的值符合预设条件时，根据所述所有主元素的相对分子量与所述纯矿物的相对分子量之比，计算获得所述所有主元素在纯矿物中的系数。6.根据权利要求1
‑
5任一项中所述的方法，其特征在于，所述建立各所述主元素的相对分子量与所述纯矿物的分子式模型中的...
【专利技术属性】
技术研发人员：廖东良，赵文杰，李永杰，吴海燕，王卫，
申请(专利权)人：中国石油化工股份有限公司石油工程技术研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。