一种地层组分最优化确定方法及装置制造方法及图纸

本发明专利技术提供了一种地层组分最优化确定方法及装置，属于油气勘探测井技术领域。所述方法包括：根据待测地层的岩心分析资料及地质条件建立地层岩石组分模型，确定参与模型确定的测井曲线；确定参与模型确定的测井曲线对应的测井响应方程表达式；解析并记录存储测井响应方程表达式，建立最优化问题的目标函数，并通过迭代算法求解目标函数确定待测地层的最优组分含量。本发明专利技术通过建立地层岩石组分模型并确定对应的测井响应方程，再通过表达式解析法解析并记录存储测井响应方程表达式，然后建立最优化问题的目标函数，通过迭代算法求解获得待测地层的最优组分含量，不仅能够对用户自定义的测井响应方程进行最优化处理，而且具有较高的处理精度。

【技术实现步骤摘要】

本专利技术涉及一种地层组分最优化确定方法及装置，属于油气勘探测井

技术介绍
地层组分定量确定是储层测井解释的基础和关键。通过对地层矿物、流体含量的确定，并据此确定地层孔隙度、饱和度等参数，能够直接指导有效储层识别、油气层预测分析以及油田区域评价，对油气田勘探开发具有十分重要的意义。目前，利用测井资料确定地层组分含量通常采用的方法是测井最优化处理方法，它是将所有的测井信息、误差及地区地质经验综合成一个多维信息复合体，利用数学领域的最优化技术求取满足地层条件的最优解释结果。用于测井的最优化处理方法的基础是依据地层模型建立的各类测井响应方程，该测井响应方程表征了测井响应与地层特性参数之间的定量关系，通过测井响应方程能够确定最优化处理结果与实际地层的符合程度。但是在现有的最优化处理方法中，测井响应方程普遍只是采用预先定义的几类固定的形式。尽管在这些固定的形式在设计时都尽量为用户提供了多种模型方程，但在实际应用中仍存在以下问题：首先，已有的测井响应方程都是基于传统经典模型建立的，因而其适用范围有限，在很多复杂类型储层中应用效果不理想；其次，预定的模型公式无法在所有油田、地层都适用，因为不同研究区存在其特定的形式；最后，在测井解释中，通常还包括了很多根据岩石物理分析资料建立的经验公式，现有的最优化处理方法中欠缺对这些公式的处理能力，因而限制了该方法的应用范围和处理精度。
技术实现思路
本专利技术为解决现有的测井最优化处理方法存在的测井响应方程模型固定、适用范围有限、处理精度较差的问题，进而提出了一种地层组分最优化确定方法及装置，具体包括如下的技术方案：一种地层组分最优化确定方法，包括：根据待测地层的岩心分析资料及地质条件建立地层岩石组分模型，并确定参与模型确定的测井曲线；确定所述参与模型确定的测井曲线对应的测井响应方程表达式；解析并记录存储所述测井响应方程表达式；根据经过解析的所述响应方程表达式建立最优化问题的目标函数，并通过迭代算法求解所述目标函数，以确定所述待测地层的最优组分含量。在本专利技术所述的地层组分最优化确定方法中，所述测井曲线包括：自然伽马测井曲线、深侧向电阻率测井曲线、浅侧向电阻率测井曲线、密度测井曲线、中子测井曲线或元素俘获能谱测井曲线中的至少一种。在本专利技术所述的地层组分最优化确定方法中，所述测井响应方程为：其中，tci表示测井响应方程的值；表示待测地层的地层矿物及流体的组分含量；表示测井响应方程的表达式形式，并包括地层矿物和流体组分变量、数字、运算符以及参数符号。在本专利技术所述的地层组分最优化确定方法中所述解析并记录存储各响应方程表达式包括：对所述测井响应方程的表达式进行解析，并转义为用于计算机运算的后缀表达式；将所述后缀表达式中的每个元素通过数据结构体记录并保存为动态数组的存储结构；遍历所述后缀表达式的存储结构，并根据求导规则确定每个元素之间的组合关系，从而确定所述后缀表达式的偏导形式。在本专利技术所述的地层组分最优化确定方法中，所述优化问题的目标函数包括：其中，表示目标函数；v*表示使目标函数获得最小值时的取值；tci表示根据所述测井响应方程确定的各类测井方法的响应方程值；tmi表示实际测井测量响应值；wi表示所述测井曲线在最优化模型中的权重系数；n表示求解的测井曲线的数量。在本专利技术所述的地层组分最优化确定方法中，所述地层岩石组分模型还包括附加约束条件，所述附加约束条件为：其中，表示约束条件；Ck表示约束条件系数矩阵；表示待测地层的地层矿物及流体的组分含量；bk表示约束条件边界。在本专利技术所述的地层组分最优化确定方法中，通过迭代算法求解所述目标函数包括：通过惩罚函数法将所述目标函数转化为无约束问题的表达式：其中，M表示惩罚因子；并表示惩罚函数，当满足约束条件时，惩罚项当不满足约束条件时，惩罚项且随M的增大而增大；cn表示约束条件的数目。在本专利技术所述的地层组分最优化确定方法中，通过迭代算法求解所述目标函数还包括：通过Levenberg-Marquardt算法进行迭代增量可以表示为；(JTJ+μ·I)·h＝-JTR其中，J表示R的Jacobi矩阵；I表示单位矩阵；μ表示阻尼因子，在本专利技术所述的地层组分最优化确定方法中，通过迭代算法求解所述目标函数还包括：对每个处理深度点的一组初始设定的地层组分含量值进行第一次迭代，将各地层组分含量值代入解析存储的偏导表达式中，获得矩阵J；通过迭代增量求取当前迭代增量并转入到下一次迭代；当迭代地层组分值满足精度要求且满足附加约束条件时，即求解得到最优的地层组分含量。一种地层组分最优化确定装置，包括：测井曲线确定单元，用于根据待测地层的岩心分析资料及地质条件建立地层岩石组分模型，并确定参与模型确定的测井曲线；响应方程确定单元，用于确定所述参与模型确定的测井曲线对应的测井响应方程表达式；解析单元，用于解析并记录存储所述测井响应方程表达式；最优含量确定单元，根据经过解析的所述响应方程表达式建立最优化问题的目标函数，并通过迭代算法求解所述目标函数，以确定所述待测地层的最优组分含量。本专利技术的有益效果是：通过建立地层岩石组分模型并确定对应的测井响应方程，再通过表达式解析法解析并记录存储测井响应方程表达式，然后建立最优化问题的目标函数，通过迭代算法求解，从而获得待测地层的最优组分含量，不仅能够对用户自定义的测井响应方程进行最优化处理，适用范围较大，而且具有较高的处理精度。附图说明图1是以示例的方式示出了地层组分最优化确定方法的流程图。图2是实施例一所述的地层组分最优化确定方法的流程图。图3是实施例一提供的典型岩心含水饱和度-电阻增大率实验测量解结果示意图，其中的横坐标表示含税饱和度(Sw)，纵坐标表示电阻增大率(I)，圆点表示岩电实验结合，实曲线表示最优饱和度方程拟合结构，虚曲线表示Archie公式拟合结果。图4是实施例一提供的表达式解析的流程图。图5是实施例一提供的A井地层矿物及流体组分含量结果与实验分析结果对比图。图6是实施例二提供的地层组分最优化确定装置的结构图。具体实施方式结合图1所示，本实施例提供的地层组分最优化确定方法包括：步骤11，根据待测地层的岩心分析资料及地质条件建立地层岩石组分模型，并确定参与模型确定的测井曲线。对于预定的待测地层，可以获取的测井数据及实验分析数据包括：待测地层的常规、电成像等测井数据以及岩心、流体实验分析数据。根据待测地层的岩心、流体分析数据及地质条件，可以确定待测地层的岩石的主要矿物组成类型、微量矿物类型以及地层流体类型，从而建立用于测井解释的地层岩石组分模型。观察待处理层段测井曲线变化规律，优选出对地层矿物或流体变化具有明显响应特征的曲线参与模型最优化处理，并保证参与确定的曲线数量应尽量多于模型中矿物、流体组分的数量。步骤12，确定参与模型本文档来自技高网...

【技术保护点】
一种地层组分最优化确定方法，其特征在于，包括：根据待测地层的岩心分析资料及地质条件建立地层岩石组分模型，并确定参与模型确定的测井曲线；确定所述参与模型确定的测井曲线对应的测井响应方程表达式；解析并记录存储所述测井响应方程表达式；根据经过解析的所述响应方程表达式建立最优化问题的目标函数，并通过迭代算法求解所述目标函数，以确定所述待测地层的最优组分含量。

【技术特征摘要】
1.一种地层组分最优化确定方法，其特征在于，包括：
根据待测地层的岩心分析资料及地质条件建立地层岩石组分模型，并确定参与模
型确定的测井曲线；
确定所述参与模型确定的测井曲线对应的测井响应方程表达式；
解析并记录存储所述测井响应方程表达式；
根据经过解析的所述响应方程表达式建立最优化问题的目标函数，并通过迭代算
法求解所述目标函数，以确定所述待测地层的最优组分含量。
2.如权利要求1所述的方法，其特征在于，所述测井曲线包括：自然伽马测井
曲线、深侧向电阻率测井曲线、浅侧向电阻率测井曲线、密度测井曲线、中子测井曲
线或元素俘获能谱测井曲线中的至少一种。
3.如权利要求2所述的方法，其特征在于，所述测井响应方程为：
tci=fi(v&RightA...

【专利技术属性】
技术研发人员：冯周，李宁，武宏亮，王华峰，冯庆付，王克文，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11