一种深部地温场预测方法、设备及存储介质技术

本发明专利技术提供一种深部地温场预测方法、设备及存储介质，包括以下步骤：将研究区或相邻区域各测井温度‑电阻率数据对均等分成N段，且归一化处理各层段电阻率；推导不同井位不同层段的最优校正温度、本征归一化电阻率、温度校正系数，并将三者与对应深度构建数据集，进行回归分析，分别获取随深度的变化为T

A prediction method, equipment and storage medium of deep geothermal field

【技术实现步骤摘要】
一种深部地温场预测方法、设备及存储介质
本专利技术涉及地温场预测领域，尤其涉及一种深部地温场预测方法、设备及存储介质。
技术介绍
温度是地球内部的关键特征之一，对它的了解决定了我们研究基础地球科学问题和应用地热问题的能力。因此，最大限度地准确估计地下空间温度分布特征显得极为重要。目前，地球内部的温度估计通常是基于钻孔温度测井或热流梯度数据。前者根据地表不规则分布的几口井中测量的温度测井曲线进行空间插值，这常常导致相当大的估计错误。而后者对于温度模型的建立是基于对研究区域横向边界热流稳定状态的假设，以及对区域上下边界热流值(温度)基于先验信息的推测。由于这些数值通常只能得到非常粗略的估计(特别是在下边界)，因此在此基础上建立温度模型也会有相当大的误差。因此，现有的温度估计方法无法准确预测钻井未到达深度的温度，更无法有效提供井间空间的温度展布。而研究表明，使用随温度变化的替代参数可以克服这一困难，而在各种物性参数中，岩石电阻率随温度的变化相对最为敏感，因此，可使用电阻率去解决这一科学问题，且地下深部空间电阻率分布可以通过电磁法等地球物理手段获取。然而目前，利用电阻率进行温度预测是基于纯经验公式的使用，且其有效性被假设为不随空间位置的变化，即经验公式中的各项参数在任意地质环境及背景下均假设为定值，很明显这种方法是不合理的。
技术实现思路
针对目前地下温度预测尤其深部空间温度预测所建立的温度模型误差较大，且基于纯经验公式预测地温场的不合理性，本专利技术提出了一种基于最优温度系数的深部地温场预测方法，通过研究区或相邻区域钻孔测井电阻率-温度数据对，推导构建地下空间不同层段(不同深度)归一化电阻率与温度之间的精确关系表征，并结合电磁精细电阻率结构反演，实现地下空间温度场预测。为了实现上述目的，本专利技术提出了一种深部地温场预测方法，所述深部地温场预测方法包括以下步骤：获取研究区或相邻区域m(m>＝1的整数)口井的温度-电阻率数据对，依据深度将所述温度-电阻率数据对等分成N(N>＝2的整数)段，并将每口井各段电阻率进行归一化处理，得到多组归一化电阻率-温度数据集；根据所述多组归一化电阻率-温度数据集，通过两组归一化电阻率-温度关系约束方程，计算得到不同井位不同层段的最优校正温度、本征归一化电阻率和温度校正系数；将所述最优校正温度、所述本征归一化电阻率和所述温度校正系数在不同井位不同层段的取值分别与对应层段深度构建目标数据集，对所述目标数据集进行回归分析，得到所述最优校正温度、所述本征归一化电阻率和所述温度校正系数随深度的变化关系，分别表示为T0(z),ρNT0(z),αT0(z)；对研究区电磁数据体进行精细反演，获取地下空间各节点P(x,z)的反演电阻率ρinv(x,z)的分布特征，并将所述反演电阻率剖面随深度等分成M(M>＝2的整数)段，对每段的反演电阻率进行归一化处理，得到不同层段不同节点的归一化反演电阻率ρNinv(x,z)；其中，x为反演电阻率剖面的横向点位距离，z为反演电阻率剖面的纵向深度；基于所述T0(z),ρNT0(z),αT0(z)和所述ρNinv(x,z)，得到研究区地下不同层段不同节点的归一化反演电阻率与温度之间的关系表征，根据所述ρNinv(x,z)和所述关系表征预测研究区地下空间温度场的展布特征。优选地，所述将每口井各段电阻率进行归一化处理，其中，归一化处理公式为：其中，ρlog(i,j)为每口井第i层段内第j个测点电阻率测井数据，ρlog(i)为每口井第i层段内所有测点测井电阻率数据集，max[ρlog(i)]表示每口井第i层段内所有测井电阻率数据集的最大值，ρNorm(i,j)即为每口井第i层段内第j个测点归一化测井电阻率，i＝1,2,3…N。优选地，所述优选两组归一化电阻率-温度关系约束方程，根据所述多组归一化电阻率-温度数据集计算得到不同井位不同层段的最优校正温度、本征归一化电阻率和温度校正系数的步骤，包括：基于改进传统经验公式，得到与温度T(i)变化呈线性变化关系公式：其中，ρNorm(i)为每口井第i层段内所有测点归一化测井电阻率数据集，T(i)为每口井第i层段内所有测点测井温度数据集；对所述线性变化关系公式进行转换，得到转换公式：计算所述转换公式的斜率B和截距A，得到：根据岩石电阻率与温度的对数相关关系，在任意层段内，测量归一化电阻率ρNorm(i)与温度T(i)之间的关系，得到归一化电阻率与温度关系公式为：ρNorm(i)＝-Cln(T(i))+D(4)其中，C与D均为与岩层结构、地质背景、电性特征相关的系数；将T＝Tbest代入式(2)和式(4)进行折算，同时结合(3)式得到：将式(5)对Tbest求导，整理后得：求解式(6)，得到：上述A、B、C和D均由测井不同层段归一化电阻率ρNorm(i)与温度T(i)数据对获取，故在计算出每层段最优校正温度T0＝Ttest的基础上，由式(4)和式(3)分别计算得到每层段最优温度条件下的本征归一化电阻率ρNT0和温度校正系数αT0。优选地，所述对每段的反演电阻率进行归一化处理，得到不同层段不同节点的归一化反演电阻率ρNinv(x,z)，其中归一化处理公式为：其中，ρinv(k,x,z)为第k(k＝1,2,3,…M)层段中P(x,z)位置处反演电阻率数据；ρinv(k)为第k层段内所有反演电阻率数据集，max[ρinv(k)]为第k层段内所有反演电阻率数据集的最大值，ρNinv(k,x,z)为第k层段中P(x,z)位置处的归一化反演电阻率ρNinv(x,z)。优选地，所述基于所述T0(z),ρNT0(z),αT0(z)和所述ρNinv(x,z)，得到研究区地下不同层段不同节点的归一化反演电阻率与温度之间的关系表征，根据所述ρNinv(x,z)和所述关系表征预测研究区地下空间温度场的展布特征，其中，所述关系表征为：其中，T(x,z)为各节点温度预测值。此外，为了实现上述目的，本专利技术还提供了一种深部地温场预测设备，所述深部地温场预测设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的深部地温场预测程序，所述深部地温场预测程序被所述处理器执行时实现上述深部地温场预测方法的步骤。此外，为了实现上述目的，本专利技术还提供了一种存储介质，所述存储介质上存储有深部地温场预测程序，所述深部地温场预测程序被处理器执行时实现上述深部地温场预测方法的步骤。本专利技术所采取的技术方案带来的有益效果是：本专利技术基于最优温度系数建立研究区地下不同层段不同节点的归一化反演电阻率与温度之间的关系表征，根据关系表征将地下介质宏观电阻率特征精确转换为直观的温度场分布，实用性强，预测范围广且深，预测值的准确率达到了85.69％～97.75％。附图说明图1为本专利技术的一种基于最优温度系数的本文档来自技高网...

【技术保护点】
1.一种深部地温场预测方法，其特征在于，所述深部地温场预测方法包括以下步骤：/n获取研究区或相邻区域m口井的温度-电阻率数据对，依据深度将所述温度-电阻率数据对等分成N段，并将每口井各段电阻率进行归一化处理，得到多组归一化电阻率-温度数据集；/n根据所述多组归一化电阻率-温度数据集，通过两组归一化电阻率-温度关系约束方程，计算得到不同井位不同层段的最优校正温度、本征归一化电阻率和温度校正系数；/n将所述最优校正温度、所述本征归一化电阻率和所述温度校正系数在不同井位不同层段的取值分别与对应层段深度构建目标数据集，对所述目标数据集进行回归分析，得到所述最优校正温度、所述本征归一化电阻率和所述温度校正系数随深度的变化关系，分别表示为T

【技术特征摘要】
1.一种深部地温场预测方法，其特征在于，所述深部地温场预测方法包括以下步骤：
获取研究区或相邻区域m口井的温度-电阻率数据对，依据深度将所述温度-电阻率数据对等分成N段，并将每口井各段电阻率进行归一化处理，得到多组归一化电阻率-温度数据集；
根据所述多组归一化电阻率-温度数据集，通过两组归一化电阻率-温度关系约束方程，计算得到不同井位不同层段的最优校正温度、本征归一化电阻率和温度校正系数；
将所述最优校正温度、所述本征归一化电阻率和所述温度校正系数在不同井位不同层段的取值分别与对应层段深度构建目标数据集，对所述目标数据集进行回归分析，得到所述最优校正温度、所述本征归一化电阻率和所述温度校正系数随深度的变化关系，分别表示为T0(z),ρNT0(z),αT0(z)；
对研究区电磁数据体进行精细反演，获取地下空间各节点P(x,z)的反演电阻率ρinv(x,z)的分布特征，并将所述反演电阻率剖面随深度等分成M段，对每段的反演电阻率进行归一化处理，得到不同层段不同节点的归一化反演电阻率ρNinv(x,z)；其中，x为反演电阻率剖面的横向点位距离，z为反演电阻率剖面的纵向深度；
基于所述T0(z),ρNT0(z),αT0(z)和所述ρNinv(x,z)，得到研究区地下不同层段不同节点的归一化反演电阻率与温度之间的关系表征，根据所述ρNinv(x,z)和所述关系表征预测研究区地下空间温度场的展布特征。


2.如权利要求1所述的深部地温场预测方法，其特征在于，所述将每口井各段电阻率进行归一化处理，其中，归一化处理公式为：



其中，ρlog(i,j)为每口井第i层段内第j个测点电阻率测井数据，ρlog(i)为每口井第i层段内所有测点的测井电阻率数据集，max[ρlog(i)]表示每口井第i层段内所有测井电阻率数据集的最大值，ρNorm(i,j)即为每口井第i层段内第j个测点归一化测井电阻率，i＝1,2,3…N。


3.如权利要求1所述的深部地温场预测方法，其特征在于，所述优选两组归一化电阻率-温度关系约束方程，根据所述多组归一化电阻率-温度数据集计算得到不同井位不同层段的最优校正温度、本征归一化电阻率和温度校正系数的步骤，包括：
基于改进传统经验公式，得到与温度T(i)变化呈线性变化关系公式：



其中，ρNorm(i)为每口井第i层段内所有测点归一化测井电阻率数据集，T(i)为每口井第i层段内所有测点测井温度数据集；
对所述线性变化关系公式进行转换，得到转换公式：计算所述转换公式的斜率B和截...

【专利技术属性】
技术研发人员：黄国疏，胡祥云，蔡建超，马火林，陈斌，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：湖北;42