[发明专利]一种深部地温场预测方法、设备及存储介质

说明书

本发明提供一种深部地温场预测方法、设备及存储介质，包括以下步骤：将研究区或相邻区域各测井温度‑电阻率数据对均等分成N段，且归一化处理各层段电阻率；推导不同井位不同层段的最优校正温度、本征归一化电阻率、温度校正系数，并将三者与对应深度构建数据集，进行回归分析，分别获取随深度的变化为T₀(z)，ρ_NT0(z)，α_T0(z)；对研究区电磁数据体进行电阻率反演并等分成M段，归一化处理不同层段的电阻率为ρ_Ninv(x，z)；基于T₀(z)，ρ_NT0(z)，α_T0(z)和ρ_Ninv(x，z)得到研究区地下不同层段不同节点的归一化反演电阻率与温度之间的关系表征，根据关系表征预测研究区地下空间温度场的展布特征。本发明能够将地下介质宏观电阻率特征精确转换为直观的温度场分布，实用性强，预测范围广且深。

技术领域

本发明涉及地温场预测领域，尤其涉及一种深部地温场预测方法、设备及存储介质。

背景技术

温度是地球内部的关键特征之一，对它的了解决定了我们研究基础地球科学问题和应用地热问题的能力。因此，最大限度地准确估计地下空间温度分布特征显得极为重要。

目前，地球内部的温度估计通常是基于钻孔温度测井或热流梯度数据。前者根据地表不规则分布的几口井中测量的温度测井曲线进行空间插值，这常常导致相当大的估计错误。而后者对于温度模型的建立是基于对研究区域横向边界热流稳定状态的假设，以及对区域上下边界热流值(温度)基于先验信息的推测。由于这些数值通常只能得到非常粗略的估计(特别是在下边界)，因此在此基础上建立温度模型也会有相当大的误差。

因此，现有的温度估计方法无法准确预测钻井未到达深度的温度，更无法有效提供井间空间的温度展布。而研究表明，使用随温度变化的替代参数可以克服这一困难，而在各种物性参数中，岩石电阻率随温度的变化相对最为敏感，因此，可使用电阻率去解决这一科学问题，且地下深部空间电阻率分布可以通过电磁法等地球物理手段获取。然而目前，利用电阻率进行温度预测是基于纯经验公式的使用，且其有效性被假设为不随空间位置的变化，即经验公式中的各项参数在任意地质环境及背景下均假设为定值，很明显这种方法是不合理的。

发明内容

针对目前地下温度预测尤其深部空间温度预测所建立的温度模型误差较大，且基于纯经验公式预测地温场的不合理性，本发明提出了一种基于最优温度系数的深部地温场预测方法，通过研究区或相邻区域钻孔测井电阻率-温度数据对，推导构建地下空间不同层段(不同深度)归一化电阻率与温度之间的精确关系表征，并结合电磁精细电阻率结构反演，实现地下空间温度场预测。

为了实现上述目的，本发明提出了一种深部地温场预测方法，所述深部地温场预测方法包括以下步骤：

获取研究区或相邻区域m(m＞＝1的整数)口井的温度-电阻率数据对，依据深度将所述温度-电阻率数据对等分成N(N＞＝2的整数)段，并将每口井各段电阻率进行归一化处理，得到多组归一化电阻率-温度数据集；

根据所述多组归一化电阻率-温度数据集，通过两组归一化电阻率-温度关系约束方程，计算得到不同井位不同层段的最优校正温度、本征归一化电阻率和温度校正系数；

将所述最优校正温度、所述本征归一化电阻率和所述温度校正系数在不同井位不同层段的取值分别与对应层段深度构建目标数据集，对所述目标数据集进行回归分析，得到所述最优校正温度、所述本征归一化电阻率和所述温度校正系数随深度的变化关系，分别表示为T₀(z)，ρ_NT0(z)，α_T0(z)；