基于全岩化学数据的花岗岩热导率极限提升梯度预测方法技术

本发明专利技术公开了基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，涉及地热资源开采技术领域。本发明专利技术包括以下步骤：S1：从实验室花岗岩热导率和全岩化学测试数据中提取用于训练模型的样本数据；S2：对数据进行预处理；S3：基于极限梯度提升的机器学习算法，训练预测模型；S4：利用模型评价指标，评估预测模型的性能，优选预测模型并保存；S5：利用保存的预测模型，实现对新数据热导率的预测。本发明专利技术结合机器学习极限梯度提升算法，建立了利用花岗岩岩屑全岩化学数据对热导率的回归模型，实现了对花岗岩岩屑热导率的快速和精准预测，降低了时间和经济成本，解决了现有的花岗岩岩屑热导率预测方法用量大，时间和经济成本高误差大的问题。问题。问题。

【技术实现步骤摘要】
基于全岩化学数据的花岗岩热导率极限提升梯度预测方法


[0001]本专利技术属于地热资源开采
，特别是涉及基于全岩化学数据的花岗岩热导率极限提升梯度预测方法。

技术介绍

[0002]花岗岩属于酸性岩浆岩中的侵入岩，主要矿物为石英、钾长石和酸性斜长石，次要矿物则为黑云母、角闪石，有时还有少量辉石。作为大陆的标志性岩石，花岗岩构成大陆上部地壳的基础，且花岗岩的形成过程通常与大陆的构造作用、变质作用和成矿作用密切相关。其热物理性质（热扩散率和热导率）对于探讨地壳热结构、定量计算各类地质事件中的温度分布及热状态等具有重要意义。
[0003]花岗岩热导率在干热岩地热资源数值模拟中意义重大。但是干热岩地热钻孔中仅在有限的深度层段进行取心工作，可用于获取热导率的只有随钻井液排出的花岗岩岩屑。如何利用岩屑获取较为准确的热导率值，一直是研究的热点和难点。岩屑热导率的获取可以分为实验室直接测量和间接估算两种方法。但是，由于实验室测量岩屑的热导率数值的方法存在一定的局限性：1、要求岩屑的量在100mL左右；2、在进行岩屑热导率测试之前需要对岩屑进一步处理，比如利用光学扫描法测量时，还需要对岩屑进行压实等额外操作；3、室内测试岩屑热导率效率低，线源法测试效率最快也需要15分钟。
[0004]对于岩屑量较少，难以进行实验室的热导率测量时，常用的方法是利用岩屑的矿物含量数据来预测起热导率数值，常用的计算模型有几何平均、调和平均模型，但是这些常用的数学模型的计算基础是基于理想矿物排列和接触假设条件下而推导的，与花岗岩的实际矿物空间排列和接触关系相差甚远，导致预测效果不佳。此外，利用花岗岩全岩化学进行热导率预测的实践较少，只有Jennings，et al；（2019）基于340个深成岩全岩化学和热导率数据，建立了经验公式，且近年来，机器学习算法逐渐广泛的应用于各种材料导热性能的探索和设计中，但是用于岩石热导率预测的方面还鲜有涉及。
[0005]为了克服上述现有技术的不足，本专利技术提供了基于全岩化学数据的花岗岩热导率极限提升梯度预测方法。

技术实现思路

[0006]本专利技术的目的在于提供基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，本专利技术结合机器学习极限梯度提升算法，建立了利用花岗岩岩屑全岩化学数据对热导率的回归模型，实现了对花岗岩岩屑热导率的快速和精准预测，降低了时间和经济成本，解决了现有的常规花岗岩岩屑热导率预测方法用量大，时间和经济成本高以及误差大的问题。
[0007]为解决上述技术问题，本专利技术是通过以下技术方案实现的：
本专利技术为基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，包括以下步骤：S1：从实验室花岗岩热导率和全岩化学测试数据中提取用于训练模型的样本数据；S2：对数据进行预处理；S3：基于极限梯度提升的机器学习算法，训练多个花岗岩热导率预测模型；S4：利用模型评价指标，评估预测模型的性能，优选预测模型并保存；S5：利用保存的预测模型，实现对新数据热导率的预测。
[0008]进一步地，S1中提取样本数据的具体方法为，从实验室测试数据中获取花岗岩热导率和全岩化学数据组成数据集，以热导率作为有监督学习的目标集，并以全岩化学数据作为特征集。
[0009]进一步地，S2具体包括以下步骤：S21：对数据进行清洗和统计分析，得到特征数据；S22：对特征数据进行归一化处理；S23：对归一化后的特征数据集进行降维处理。
[0010]进一步地，S22中归一化处理公式为：其中，为归一化后的数据，为这一特征数据中的最小值，为这一特征数据中的最大值。
[0011]进一步地，S23中降维处理的具体方法为，通过对各个特征数据之间的相关性进行分析，以及对各个特征数据的主成分进行分析，从而对特征数据集进行降维判断和处理。
[0012]进一步地，S3的具体方法为：S31：将归一化和降维后的数据随机抽取80%数据用作模型训练部分，剩余的20%数据作为模型测试集；S32：利用极限梯度提升的机器学习算法，通过多次训练和测试，得到花岗岩热导率预测模型。
[0013]进一步地，S4具体为，利用模型评价指标，对多个花岗岩热导率预测模型的性能进行评估，对性能最好的花岗岩热导率预测模型进行保存。
[0014]进一步地，S4中模型评价指标包括均方误差（MSE）、平均绝对误差（MAPE）、均方根误差（RMSE）和相关系数（R2）。
[0015]进一步地，均方误差MSE的计算公式为：平均绝对误差MAPE的计算公式为：
均方根误差RMSE的计算公式为：相关系数R2的计算公式为：其中，为测试集中样品室内试验测试得到的花岗岩热导率，为模型预测的花岗岩热导率，为模型预测的花岗岩平均值。
[0016]进一步地，S5具体为，利用保存好的最佳花岗岩热导率预测模型对新加入的花岗岩全岩化学数据进行热导率预测。
[0017]本专利技术具有以下有益效果：1、本专利技术将极限提升梯度学习算法应用于岩石热导率预测方面，机器学习能够通过相关可靠数据来快速、自动地创建精确预测模型，通过机器学习算法来精确预测花岗岩岩屑的热导率，将有利于提高热导率的预测效率和精度，减少时间和经济成本，且本专利技术为进一步研究建立岩石热导率预测模型提供参考。
[0018]2、本专利技术利用多项模型评价指标，评估预测模型的性能，通过均方误差MSE、平均绝对误差MAPE、均方根误差RMSE和相关系数R2的综合考量，可得到最优的花岗岩热导率预测模型，进而大大提高预测的准确性。
附图说明
[0019]为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例描述所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0020]图1为流程示意图；图2为本专利技术极限梯度提升算法和现有经验公式热导率预测效果对比图。
具体实施方式
[0021]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述。
[0022]请参阅图1所示，本专利技术为基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，包括以下步骤：
S1：从实验室花岗岩热导率和全岩化学测试数据中提取用于训练模型的样本数据，具体为，从实验室测试数据中获取花岗岩热导率和全岩化学数据组成数据集，以热导率作为有监督学习的目标集，并以全岩化学数据作为特征集；S2：对数据进行预处理，具体包括以下步骤：S21：对数据进行清洗和统计分析，得到特征数据；S22：对特征数据进行归一化处理；S23：对归一化后的特征数据集进行降维处理；S22中归一化处理公式为：其中，为归一化后的数据，为这一特征数据中的最小值，为这一特征数据中的最大值；S23中降维处理的具体方法为，通过对各个特征数据之间的相关性进行分析，以及对各个特征数据的主成分进行分析，从而对特征数据集进行降维判断和处理。
[0023]S3：基于极限梯度提升的机器学习算本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，其特征在于：包括以下步骤：S1：从实验室花岗岩热导率和全岩化学测试数据中提取用于训练模型的样本数据；S2：对数据进行预处理；S3：基于极限梯度提升的机器学习算法，训练多个花岗岩热导率预测模型；S4：利用模型评价指标，评估预测模型的性能，优选预测模型并保存；S5：利用保存的预测模型，实现对新数据热导率的预测。2.根据权利要求1所述的基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，其特征在于，所述S1中提取样本数据的具体方法为，从实验室测试数据中获取花岗岩热导率和全岩化学数据组成数据集，以热导率作为有监督学习的目标集，并以全岩化学数据作为特征集。3.根据权利要求2所述的基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，其特征在于，所述S2具体包括以下步骤：S21：对数据进行清洗和统计分析，得到特征数据；S22：对特征数据进行归一化处理；S23：对归一化后的特征数据集进行降维处理。4.根据权利要求3所述的基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，其特征在于，所述S22中归一化处理公式为：其中，为归一化后的数据，为这一特征数据中的最小值，为这一特征数据中的最大值。5.根据权利要求4所述的基于全岩化学数据的花岗岩热导率极限提升梯度预测方法，其特征在于，所述S23中降维处理的具体方法为，通过对各个特征数据之间的相关性进行分析，以及对各个特征数据的主成分进行分析，从而对特征数据集进行...

【专利技术属性】
技术研发人员：汪虎，李朦，徐玳笠，陈康，李勇，余如洋，
申请(专利权)人：数皮科技湖北有限公司，
类型：发明
国别省市：