高温高压岩心电阻率标定的深部地温场预测方法及装置

1.本发明涉及地温场预测领域，尤其涉及一种高温高压岩心电阻率标定的深部地温场预测方法及装置。


背景技术：

2.温度是地球内部的关键特征之一，对深部温度属性的预测能力直接决定了我们研究基础地球科学和应用地热资源的能力。因此，准确快速估算地下空间温度分布特征显得极为重要。
3.目前，获取地球内部温度的方式主要有两大类：直接测量及间接预测。直接测量主要是通过钻孔测井获取沿深度方向温度特征，并基于不规则分布的钻孔测温进行空间插值。但直接测量存在两个问题：一方面，若钻孔太多则成本太高；另一方面，若钻孔太少则通过测井温度插值常常导致相当大的误差，尤其在地质构造复杂区域。而间接预测主要是基于地球化学或地球物理手段反推地温场。其中，地球化学手段主要是利用各种地球化学温标预测温度场，即通过收集地球化学同位素、流体或气体等数据反推温度特征，但地球化学地温计只能大致预估热储温度的分布范围，无法将预测温度与深度匹配，更无法获取区域性温度场展布。而基于地球物理探测的地温场预测方法，目前主要是通过搭建地球物理参数(如电阻率、波速等)与温度的耦合关系并基于地球物理探测反推地温场，而目前利用地球物理参数(如电阻率)预测温度场目前主要是使用纯经验公式，其有效性被假设为不随空间位置的变化，即经验公式中的各项参数在任意地质环境及深度下均假设为定值，但这种方法很明显是不合理的。因此，现有的地下深部温度场间接估算方法无法准确预测钻孔未到达深度的温度及井间空间的温度，更无法准确预测区域性深部地温场展布特征。


技术实现要素：

4.针对目前地下温度场预测尤其深部空间温度场预测模型误差较大，且基于纯经验公式预测地温场的不合理性。本发明提出了一种高温高压岩心电阻率标定的深部地温场预测方法，通过高温高压岩心电阻率测试推导构建不同岩性不同深度的岩心归一化电阻率与温度、压力关系表征；反演研究区三维密度分布特征并换算成上覆地层压力分布，同时反演电磁数据获取三维电阻率并基于最优化理论优选出岩心尺度电阻率标定场地尺度电阻率的最佳转换公式，将反演电阻率归一化处理；最终结合研究区不同深度归一化电阻率与温度、压力之间精确关系表征及上覆地层压力与归一化电阻率实现地下深部温度场预测。
5.为了实现上述目的，本发明提出了一种高温高压岩心电阻率标定的深部地温场预测方法及装置。
6.根据本发明的一方面，一种高温高压岩心电阻率标定的深部地温场预测方法包括以下步骤：
7.选取研究区或邻区已有的n
l
口钻孔内不同深度、不同地层、不同岩性的岩心m颗，对m颗岩心电阻率开展高温高压试验，获取不同温度、不同压力下岩心电阻率并进行归一化
处理，形成不同岩心归一化电阻率-温度-压力数据集；
8.根据所述不同岩心归一化电阻率-温度-压力数据集，构造归一化电阻率-温度-压力耦合关系约束方程f
rtp
，通过多元回归分析计算得到不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的压力控制系数a(i)、温度控制系数b(i)及常系数c(i)，其中，i表示第i颗岩心；
9.提取岩心样品对应钻孔处测井温度tb、上覆地层压力p
plog
随深度h的变化特征，分别视作地层背景温度tb(h)及地层背景压力pb(h)，计算深度zi处的背景温度tb(h＝zi)及背景压力pb(h＝zi)，将所述背景温度tb(h＝zi)及背景压力pb(h＝zi)代入相应岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)，获取该深度zi处标准归一化电阻率r
nstand
(zi)；
10.对研究区电磁数据体进行三维精细反演，获取地下空间各节点p(x,y,z)的反演电阻率r
inv
(x,y,z)的分布特征，根据r
inv
(x,y,zi)所处位置深度zi及地层，选取与地层具有相似岩性的岩心，并提取相应岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率r
inv
(xb,yb,zi)，基于最优化理论选取岩心尺度电阻率标定场地尺度电阻率的最佳转换公式，结合r
inv
(x,y,zi)、r
inv
(xb,yb,zi)、r
nstand
(zi)，逐点求取深度zi处归一化电阻率r
ninv
(x,y,zi)，进而逐层计算不同深度处归一化电阻率r
ninv
(x,y,z)，同时对研究区重力数据体开展三维精细反演，获取三维密度分布特征ρg(x,y,z)并计算上覆地层压力分布特征p
over
(x,y,z)；
11.根据所述不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中包含的压力控制系数a(i)、温度控制系数b(i)与常系数c(i)，及所述p
over
(x,y,z)，得到研究区不同节点处p(x,y,z)归一化反演电阻率与温度之间的精确关系表征，根据所述精确关系表征及r
ninv
(x,y,z)，逐层逐点计算不同深度处地温场t
cal
(x,y,z)展布特征。
12.优选地，所述获取不同温度、不同压力下岩心电阻率并进行归一化处理，形成不同岩心归一化电阻率-温度-压力数据集，其中，岩心电阻率归一化处理公式为：
[0013][0014]
i＝1,2,3,
…
m，r
ntest
(i,tj,pj)为第i颗岩心某一测试温度tj、压力pj下岩心归一化电阻率；r
test
(i,tj,pj)为第i颗岩心某一测试温度tj、压力pj下岩心电阻率值；r
test
(i,t,p)为第i颗岩心所有测试温度、压力下岩心电阻率数据集；max[r
test
(i,t,p)]为第i颗岩心所有测试温度、压力下岩心电阻率数据集的最大值。
[0015]
优选地，所述通过多元回归分析计算得到不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的压力控制系数a(i)、温度控制系数b(i)及常系数c(i)，其中，不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)的表达式为：
[0016]frtp
(i)＝f(r
ntest
(i,t,p),t(i),p(i))
[0017]
其中，r
ntest
(i,t,p)、t(i)及p(i)三者之间耦合关系表示为：
[0018][0019]
i＝1,2,3,
…
m，f
rtp
(i)表示第i颗岩心归一化电阻率-温度-压力耦合关系约束方程，为r
ntest
(i,t,p)、t(i)及p(i)的函数，表示为f(p
ntest
(i,t,p),t(i),p(i))；r
ntest
(i,t,p)表示第i颗岩心在不同测试温度及压力下的归一化电阻率数据集；p(i)表示第i颗岩心高温高压测试中压力数据集；t(i)表示第i颗岩心高温高压测试中温度数据集；p0表示常压。
[0020]
优选地，所述提取岩心样品对应钻孔处测井温度tb、上覆地层压力p
plog
随深度h的变化特征，分别视作地层背景温度tb(h)及地层背景压力pb(h)，其中，p
plog
、tb(h)及pb(h)的表达式分别表示为：
[0021][0022]
tb(h)＝f
t
(h)
[0023]
pb(h)＝f
p
(h)
[0024]
p
plog
为上覆地层压力；ρ
log
为地层测井密度；g为重力加速度；h为地层深度；δh为测井数据间距；n
p
为深度h处至地表处测井数据节点数；f
t
(h)表示地层背景温度随深度的映射；f
p
(h)表示地层背景温度随深度的映射。
[0025]
优选地，所述计算深度zi处的背景温度tb(h＝zi)及背景压力pb(h＝zi)，将所述背景温度tb(h＝zi)及背景压力pb(h＝zi)代入相应岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)，获取该深度zi处标准归一化电阻率r
nstand
(zi)，其中，标准归一化电阻率r
nstand
(zi)的表达式表示为：
[0026][0027][0028]
i＝1,2,3
…
m，r
nstand
(zi)表示研究区zi深度处标准归一化电阻率；a
ibest
表示第i颗岩心f
rtp
(i)中的压力系数；b
ibest
表示第i颗岩心f
rtp
(i)中的温度系数；c
ibest
表示第i颗岩心f
rtp
(i)中的常系数；h
control
(i)表示利用第i颗岩心计算标准归一化电阻率r
nstand
(zi)时的适用深度范围；z
cmin
(i)≤h
control
(i)≤z
cmin
(i)表示所述适用深度范围在z
cmin
(i)至z
cmax
(i)之间；z
cmin
(i)表示所述适用深度范围的最小深度；z
cmax
(i)表示所述适用深度范围的最大深度。
[0029]
优选地，所述逐点求取深度zi处归一化电阻率r
ninv
(x,y,zi)，具体计算公式为：
[0030][0031]
k1为归一化电阻率扩大转换系数，k2为归一化电阻率缩小转换系数，两者均定义为岩心尺度归一化电阻率标定场地尺度归一化电阻率的转换系数；r
inv
(x,y,zi)表示电磁反演电阻率r
inv
(x,y,z)在深度zi处的电阻率数据集；r
inv
(xb,yb,zi)表示计算研究区归一化电阻率时，相应标定岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率；x为地下三维空
间节点横向距离，y为地下三维空间节点纵向距离，z为地下三维空间节点沿垂向深度。
[0032]
优选地，所述对研究区重力数据体开展三维精细反演，获取三维密度分布特征ρg(x,y,z)并计算上覆地层压力分布特征p
over
(x,y,z)，其中，上覆地层压力分布特征p
over
(x,y,z)的计算公式表示为：
[0033][0034]
p
over
(x,y,z)为地下空间各节点p(x,y,z)位置处基于反演密度计算的上覆地层压力；ρg(x,y,z)为地下空间各节点p(x,y,z)位置处重力数据反演地层密度；g为重力加速度；δz
inv
为反演密度数据体沿深度z方向间距；n
p
为深度z处至地表处反演密度数据节点数。
[0035]
优选地，所述根据所述精确关系表征及r
ninv
(x,y,z)，逐层逐点计算不同深度处地温场t
cal
(x,y,z)展布特征，计算公式为：
[0036][0037]
i＝1,2,3
…
m，a
ical
表示计算深度z处的温度场时，适用的第i颗岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的压力系数；b
ical
表示计算深度z处的温度场时，适用的第i颗岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的温度系数；c
ical
表示计算深度z处的温度场时，适用的第i颗岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)的常系数；h
cal
(i)表示利用第i颗岩心标定计算温度场时，其适用深度范围为z
cmin
(i)至z
cmax
(i)之间；z
cmin
(i)表示所述适用深度范围的最小深度；z
cmax
(i)表示所述适用深度范围的最大深度。
[0038]
优选地，基于最优化理论，求取最优化电阻率归一化转化系数k1＝k
best1
及k2＝k
best2
，使得利用所述最优化电阻率归一化转化系数k
best1
及k
best2
换算归一化电阻率r
ninv
(x,y,z)并反推温度场t
cal
(x,y,z)＝t
best
(x,y,z)时，提取的m颗岩心对应n
l
口钻孔位置处的预测温度与实测温度拟合度最佳，t
best
(x,y,z)即为研究区现今深部地温场分布特征。
[0039]
根据本发明的另一方面，一种基于岩心高温高压电阻率标定的深部地温场预测装置，包括以下模块：
[0040]
数据集获取模块，用于选取研究区或邻区已有的n
l
口钻孔内不同深度、不同地层、不同岩性的岩心m颗，对m颗岩心电阻率开展高温高压试验，获取不同温度、不同压力下岩心电阻率并进行归一化处理，形成不同岩心归一化电阻率-温度-压力数据集；
[0041]
多元回归分析模块，用于根据所述不同岩心归一化电阻率-温度-压力数据集，构造归一化电阻率-温度-压力耦合关系约束方程f
rtp
，通过多元回归分析计算得到不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的压力控制系数a(i)、温度控制系数b(i)及常系数c(i)，其中，i表示第i颗岩心；
[0042]
标准归一化电阻率计算模块，用于提取岩心样品对应钻孔处测井温度tb、上覆地层压力p
plog
随深度h的变化特征，分别视作地层背景温度tb(h)及地层背景压力pb(h)，计算深度zi处的背景温度tb(h＝zi)及背景压力pb(h＝zi)，将所述背景温度tb(h＝zi)及背景压
力pb(h＝zi)代入相应岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)，获取深度zi处标准归一化电阻率r
nstand
(zi)；
[0043]
深部电阻率归一化处理模块，用于对研究区电磁数据体进行三维精细反演，获取地下空间各节点p(x,y,z)的反演电阻率r
inv
(x,y,z)的分布特征，根据r
inv
(x,y,zi)所处位置深度zi及地层，选取与地层具有相似岩性的岩心，并提取相应岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率r
inv
(xb,yb,zi)，基于最优化理论选取岩心尺度电阻率标定场地尺度电阻率的最佳转换公式，结合r
inv
(x,y,zi)、r
inv
(xb,yb,zi)、r
nstand
(zi)，逐点求取深度zi处归一化电阻率r
ninv
(x,y,zi)，进而逐层计算不同深度处归一化电阻率r
ninv
(x,y,z)；
[0044]
上覆地层压力分布特征计算模块，用于对研究区重力数据体开展三维精细反演，获取三维密度分布特征ρg(x,y,z)并计算上覆地层压力分布特征p
over
(x,y,z)；
[0045]
地温场分布预测模块，用于根据所述不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中包含的压力控制系数a(i)、温度控制系数b(i)与常系数c(i)，及所述p
over
(x,y,z)，得到研究区不同节点处p(x,y,z)归一化反演电阻率与温度之间的精确关系表征，根据所述精确关系表征及r
ninv
(x,y,z)，逐层逐点计算深部地温场t
cal
(x,y,z)展布特征。
[0046]
本发明所采取的技术方案带来的有益效果是：本发明基于不同深度不同岩性的岩心高温高压电阻率试验，结合岩心尺度电阻率校正场地尺度电阻率转换公式，构建研究区不同深度归一化反演电阻率-压力-温度之间的精确关系表征，融合上覆地层压力属性，根据所述精确关系表征将岩心标定的地下介质宏观电阻率特征转换为直观的温度场分布，预测值的准确率达到了86.33％～97.28％。实用性强，预测范围广且深，对于地热资源勘探评价、可持续开发利用及地下温度监测具有重要意义。
附图说明
[0047]
下面将结合附图及实施例对本发明作进一步说明，附图中：
[0048]
图1为本发明的一种高温高压岩心电阻率标定的深部地温场预测流程图；
[0049]
图2为本发明的不同深度不同岩性岩心归一化电阻率-温度-压力关系分布图；
[0050]
图3为本发明的研究区地下空间反演电阻率剖图；
[0051]
图4为本发明的研究区地下空间上覆地层压力剖图；
[0052]
图5为本发明的研究区地下空间温度场预测剖面图；
[0053]
图6为本发明的研究区预测温度与实测温度对比验证图；
[0054]
图7为本发明的一种高温高压岩心电阻率标定的深部地温场预测装置结构图。
具体实施方式
[0055]
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
[0056]
为了使本发明的目的、技术方案和效果更加清楚的理解，现将结合附图对本发明实施方式作进一步的描述。
[0057]
实施例一：参考图1，图1为本发明的一种高温高压岩心电阻率标定的深部地温场预测方法流程图，实施例提供的一种高温高压岩心电阻率标定的深部地温场预测方法，包
括以下步骤：
[0058]
s1、选取研究区或邻区已有的n
l
(n
l
》＝1)口钻孔内不同深度、不同地层、不同岩性的岩心m颗。对m颗岩心电阻率开展高温高压试验，获取不同温度、不同压力下岩心电阻率并进行归一化处理，构建不同岩心归一化电阻率-温度-压力数据集。
[0059]
在本实施例中，步骤s1具体为：优选研究区或邻区已有的n
l
＝2口钻孔内不同深度、不同地层、不同岩性的岩心m＝6颗。对m＝6颗岩心电阻率开展高温高压试验，获取不同温度、不同压力下岩心电阻率并进行归一化处理，构建不同岩心归一化电阻率-温度-压力数据集。其中，岩心电阻率归一化处理公式为：
[0060][0061]rntest
(i,tj,pj)为第i(i＝1,2,3,
…
m)颗岩心在测试温度tj、压力pj下岩心归一化电阻率；r
test
(i,tj,pj)为第i(i＝1,2,3,
…
m)颗岩心在测试温度tj、压力pj下岩心电阻率值；r
test
(i,t,p)为第i(i＝1,2,3,
…
m)颗岩心在所有测试温度、压力下的岩心电阻率数据集；max[r
test
(i,t,p)]为第i(i＝1,2,3,
…
m)颗岩心在所有测试温度、压力下的岩心电阻率数据集的最大值。
[0062]
s2、基于所述不同岩心归一化电阻率-温度-压力数据集，构造归一化电阻率-温度-压力耦合关系约束方程f
rtp
，通过多元回归分析计算得到不同岩心(适用地下不同深度范围)f
rtp
(i)中的压力控制系数a(i)、温度控制系数b(i)及常系数c(i)，其中i表示第i颗岩心。
[0063]
在本实施例中，步骤s2具体为：基于步骤s1中不同岩心归一化电阻率-温度-压力数据集，构造归一化电阻率-温度-压力耦合关系约束方程f
rtp
，通过多元回归分析计算得到步骤s1中m＝6颗岩心(适用地下不同深度范围)归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)(参考图2)中的压力控制系数a(i)、温度控制系数b(i)及常系数c(i)。其中，不同岩心f
rtp
(i)表示为：
[0064]frtp
(i)＝f(r
ntest
(i,t,p),t(i),p(i))
ꢀꢀ
(2)
[0065]
其中，r
ntest
(i,t,p)、t(i)及p(i)三者之间耦合关系表示为：
[0066][0067]frtp
(i)表示第i(i＝1,2,3
…
m)颗岩心归一化电阻率-温度-压力耦合约束关系方程，为r
ntest
(i,t,p)、t(i)及p(i)的函数，表示为f(r
ntest
(i,t,p),t(i),p(i))；r
ntest
(i,t,p)表示第i(i＝1,2,3
…
m)颗岩心在不同测试温度及压力下的归一化电阻率数据集；p(i)表示第i(i＝1,2,3
…
m)颗岩心所有测试压力数据集；t(i)表示第i(i＝1,2,3
…
m)颗岩心所有测试温度数据集；p0表示常压，取0.1mpa。
[0068]
s3、提取所述岩心样品对应钻孔处测井温度tb、上覆地层压力p
plog
(基于测井密度计算获取)随深度h的变化特征，分别视作地层背景温度tb(h)及地层背景压力pb(h)。计算深度zi处的背景温度tb(h＝zi)及背景压力pb(h＝zi)。然后将所述背景温度tb(h＝zi)及背景压力pb(h＝zi)代入所述相应岩心(适用深度范围包含zi)的归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)，获取深度zi处标准归一化电阻率r
nstand
(zi)。
[0069]
在本实施例中，步骤s3具体为：提取步骤s1中m＝6颗岩心样品对应n
l
＝2口钻孔位
置处测井温度tb、上覆地层压力p
plog
(利用公式(4)基于测井密度计算获取)随深度h的变化特征，分别视作地层背景温度tb(h)及地层背景压力pb(h)。计算深度zi处的背景温度tb(h＝zi)及背景压力pb(h＝zi)。然后将所述背景温度tb(h＝zi)及背景压力pb(h＝zi)代入步骤s2中相应岩心(适用深度范围包含zi)的归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)，获取深度zi处标准归一化电阻率r
nstand
(zi)。其中，地层背景温度tb(h)及地层背景压力pb(h)，分别表示为(5)式和(6)式：
[0070][0071]
tb(h)＝f
t
(h)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0072]
pb(h)＝f
p
(h)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0073]
p
plog
为上覆地层压力；ρ
log
为地层测井密度；g为重力加速度；h为地层深度；δh为测井数据间距；n
p
为深度h处至地表处测井数据节点数。f
t
(h)表示地层背景温度随深度的映射(函数)。f
p
(h)表示地层背景温度随深度的映射(函数)。
[0074]
而标准归一化电阻率r
nstand
(zi)。表示为：
[0075][0076][0077]rnstand
(zi)表示研究区zi深度处标准归一化电阻率；a
ibest
表示第i(i＝1,2,3
…
m)颗岩心f
rtp
(i)中的压力系数；b
ibest
表示第i(i＝1,2,3
…
m)颗岩心f
rtp
(i)中的温度系数；c
ibest
表示第i(i＝1,2,3
…
m)颗岩心f
rtp
(i)中的常系数；h
control
(i)表示利用第i(i＝1,2,3,
…
m)颗岩心计算标准归一化电阻率r
nstand
(zi)时的适用深度范围；z
cmin
(i)≤h
control
(i)≤z
cmax
(i)表示所述适用深度范围在z
cmin
(i)至z
cmax
(i)之间；z
cmin
(i)表示所述适用深度范围的最小深度；z
cmax
(i)表示所述适用深度范围的最大深度。
[0078]
s4、对研究区电磁数据体进行三维精细反演，获取地下空间各节点p(x,y,z)的反演电阻率r
inv
(x,y,z)的分布特征，根据r
inv
(x,y,zi)所处位置深度zi及地层，优选与所述地层具有相似岩性的岩心，并提取相应岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率p
inv
(xb,yb,zi)。基于岩心尺度电阻率标定场地尺度电阻率的转换公式，结合r
inv
(x,y,zi)、r
inv
(xb,yb,zi)、r
nstand
(zi)，逐点求取深度zi处归一化电阻率r
ninv
(x,y,zi)，进而逐层计算不同深度处归一化电阻率r
ninv
(x,y,z)。
[0079]
在本实施例中，步骤s4具体为：对研究区电磁数据体进行三维精细反演，获取地下空间各节点p(x,y,z)的反演电阻率r
inv
(x,y,z)的分布特征(参考图3，是对三维反演电阻率ρ
inv
(x,y,z)进行切片，获取某条电磁测线对应的电阻率剖面)。根据r
inv
(x,y,zi)所处位置深度zi及地层，优选m＝6颗岩心中与所述地层具有相似岩性的岩心，并提取相应岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率r
inv
(xb,yb,zi)。基于岩心尺度电阻率标定场地尺度电阻率的转换公式，结合r
inv
(x,y,zi)、r
inv
(xb,yb,zi)、r
nstand
(zi)，逐点求取深度zi处归一化电阻率r
ninv
(x,y,zi)，进而逐层计算不同深度处归一化电阻率r
ninv
(x,y,z)。其中，深
度zi处基于岩心尺度电阻率标定场地尺度电阻率的转换公式表示为：
[0080][0081]
k1为归一化电阻率扩大转换系数，k2为归一化电阻率缩小转换系数，两者均定义为岩心尺度归一化电阻率标定场地尺度归一化电阻率的转换系数；r
inv
(x,y,zi)表示电磁反演电阻率r
inv
(x,y,z)在深度zi处的电阻率数据集；r
inv
(xb,yb,zi)表示计算研究区归一化电阻率时，相应标定岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率；x为地下三维空间节点横向距离，y为地下三维空间节点纵向距离，z为地下三维空间节点沿垂向深度。
[0082]
s5、对研究区重力数据体开展三维精细反演，获取三维密度分布特征ρg(x,y,z)并计算上覆地层压力分布特征p
over
(x,y,z)。
[0083]
在本实施例中，步骤s5具体为：对研究区重力数据体开展三维精细反演，获取三维密度分布特征ρg(x,y,z)并计算上覆地层压力分布特征p
over
(x,y,z)(参考图4，是对三维上覆地层压力p
pre
(x,y,z)进行切片，获取步骤s4中电磁测线对应剖面的上覆地层压力分布)。其中，上覆地层压力计算公式表示为：
[0084][0085]
p
over
(x,y,z)为地下空间各节点p(x,y,z)位置处基于反演密度计算的上覆地层压力；ρg(x,y,z)为地下空间各节点p(x,y,z)位置处重力数据反演地层密度；g为重力加速度；δz
inv
为反演密度数据体沿深度z方向间距；n
p
为深度z处至地表处反演密度数据节点数。
[0086]
s6、基于所述不同岩心(适用地下不同深度范围)归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中包含的压力控制系数a(i)、温度控制系数b(i)与常系数c(i)及所述p
over
(x,y,z)，可得到研究区不同节点处p(x,y,z)归一化反演电阻率与温度之间的精确关系表征，根据所述精确关系表征及r
ninv
(x,y,z)，可逐层逐点计算不同深度处地温场t
cal
(x,y,z)展布特征。同时，基于最优化理论，求取最优归一化电阻率转化系数k1＝k
best1
及k2＝k
best2
，使得利用所述电阻率转换系数k
best1
及k
best2
换算归一化电阻率r
ninv
(x,y,z)并反推温度场t
cal
(x,y,z)＝t
best
(x,y,z)时，提取的m颗岩心对应n
l
口钻孔位置处的预测温度与实测温度拟合度最佳，t
best
(x,y,z)即为研究区现今深部地温场分布特征。
[0087]
在本实施例中，步骤s6具体为：基于步骤s2中所述m＝6颗不同岩心(适用地下不同深度范围)归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中包含的压力控制系数a(i)、温度控制系数b(i)与常系数c(i)及步骤s5中所述p
over
(x,y,z)，可得到研究区不同节点处p(x,y,z)归一化反演电阻率与温度之间的精确关系表征；同时，基于最优化理论，求取最优归一化电阻率转化系数k1＝k
best1
及k2＝k
best2
，使得利用所述电阻率转换系数k
best1
＝0.0035及k
best2
＝0.0030换算归一化电阻率r
ninv
(x,y,z)并结合所述精确关系表征反推温度场t
cal
(x,y,z)＝t
best
(x,y,z)时，提取的m＝6颗岩心对应n
l
＝2口钻孔位置处的预测温度与实测温度拟合度最佳，t
best
(x,y,z)即为研究区现今深部地温场分布特征(参考图5，是将地
下深部预测三维温度场t
best
(x,y,z)进行切片，获取步骤s4中电磁测线对应剖面的地下温度场分布特征)。其中，计算t
cal
(x,y,z)可表示为：
[0088][0089][0090]
i＝1,2,3
…
m，a
ical
表示计算深度z处的温度场时，适用的第i颗岩心归一化电阻率-温度-压力耦合关系约束方程f
ptp
(i)中的压力系数；b
ical
表示计算深度z处的温度场时，适用的第i颗岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的温度系数；c
ical
表示计算深度z处的温度场时，适用的第i颗岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)的常系数；h
cal
(i)表示利用第i颗岩心标定计算温度场时，其适用深度范围为z
cmin
(i)至z
cmax
(i)之间；z
cmin
(i)表示所述适用深度范围的最小深度；z
cmax
(i)表示所述适用深度范围的最大深度。
[0091]
参考图6，图6为本发明的研究区预测温度与实测温度对比验证图；在本实施例中，提取研究区预测温度剖面中钻孔测温点位处温度随深度的变化曲线，并与实际测井温度进行对比。实验结果表明，d16井与d17井预测温度与实测温度的拟合优度分别为r2＝0.8633及r2＝0.9728，由此表明，预测温度准确率在86.33％～97.28％之间。
[0092]
实施例二：参考图7，图7为本发明的一种高温高压岩心电阻率标定的深部地温场预测装置结构图。实施例提供的一种岩心高温高压电阻率标定的深部地温场预测装置，包括以下模块：
[0093]
数据集获取模块1，用于选取研究区或邻区已有的n
l
口钻孔内不同深度、不同地层、不同岩性的岩心m颗，对m颗岩心电阻率开展高温高压试验，获取不同温度、不同压力下岩心电阻率并进行归一化处理，形成不同岩心归一化电阻率-温度-压力数据集；
[0094]
多元回归分析模块2，用于根据所述不同岩心归一化电阻率-温度-压力数据集，构造归一化电阻率-温度-压力耦合关系约束方程f
rtp
，通过多元回归分析计算得到不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中的压力控制系数a(i)、温度控制系数b(i)及常系数c(i)，其中，i表示第i颗岩心；
[0095]
标准归一化电阻率计算模块3，用于提取岩心样品对应钻孔处测井温度tb、上覆地层压力p
plog
随深度h的变化特征，分别视作地层背景温度tb(h)及地层背景压力pb(h)，计算深度zi处的背景温度tb(h＝zi)及背景压力pb(h＝zi)，将所述背景温度tb(h＝zi)及背景压力pb(h＝zi)代入相应岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)，获取该深度zi处标准归一化电阻率r
nstand
(zi)；
[0096]
深部电阻率归一化处理模块4，用于对研究区电磁数据体进行三维精细反演，获取地下空间各节点p(x,y,z)的反演电阻率r
inv
(x,y,z)的分布特征，根据r
inv
(x,y,zi)所处位置深度zi及地层，选取与地层具有相似岩性的岩心，并提取相应岩心对应钻孔位置p(xb,yb)深度zi处的电磁反演电阻率r
inv
(xb,yb,zi)，基于最优化理论选取岩心尺度电阻率标定场地尺度电阻率的最佳转换公式，结合r
inv
(x,y,zi)、r
inv
(xb,yb,zi)、r
nstand
(zi)，逐点求取深度zi处归一化电阻率r
ninv
(x,y,zi)，进而逐层计算不同深度处归一化电阻率r
ninv
(x,y,z)；
[0097]
上覆地层压力分布特征计算模块5，用于对研究区重力数据体开展三维精细反演，获取三维密度分布特征ρg(x,y,z)并计算上覆地层压力分布特征p
over
(x,y,z)；
[0098]
地温场分布预测模块6，用于根据所述不同岩心归一化电阻率-温度-压力耦合关系约束方程f
rtp
(i)中包含的压力控制系数a(i)、温度控制系数b(i)与常系数c(i)，及所述p
over
(x,y,z)，得到研究区不同节点处p(x,y,z)归一化反演电阻率与温度之间的精确关系表征，根据所述精确关系表征及r
ninv
(x,y,z)，逐层逐点计算深部地温场t
cal
(x,y,z)展布特征。
[0099]
本发明实施例提供一种高温高压岩心电阻率标定的深部地温场预测方法及装置，包括：开展m颗高温高压岩心电阻率测试并归一化处理电阻率；构造电阻率-温度-压力耦合关系约束方程，推导不同深度不同岩心归一化电阻率与温度、压力精确关系表征；反演研究区三维密度并换算成上覆地层压力分布；反演电磁数据获取三维电阻率，基于最优化理论选取岩心尺度电阻率标定场地尺度电阻率的最佳转换公式，将三维电阻率归一化；结合精确关系表征、三维归一化电阻率及上覆地层压力逐层逐点反推地下深部温度场展布特征。本发明能够基于高温高压岩心电阻率试验及岩心标定将地下介质宏观电阻率特征精确转换为直观温度场展布，预测范围广且深。
[0100]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
[0101]
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。词语第一、第二、以及第三等的使用不表示任何顺序，可将这些词语解释为标识。
[0102]
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。