CN114813828A - 一种新型确定含水层热物性参数的微热试验方法 - Google Patents

一种新型确定含水层热物性参数的微热试验方法 Download PDF

Info

Publication number
CN114813828A
CN114813828A CN202210440391.3A CN202210440391A CN114813828A CN 114813828 A CN114813828 A CN 114813828A CN 202210440391 A CN202210440391 A CN 202210440391A CN 114813828 A CN114813828 A CN 114813828A
Authority
CN
China
Prior art keywords
aquifer
heat
thermal
test
micro
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210440391.3A
Other languages
English (en)
Other versions
CN114813828B (zh
Inventor
赵燕容
魏裕丰
戎荣
杨义锴
王浩楠
董小松
张子民
张志豪
薛瑞丰
高秋爽
祁琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN202210440391.3A priority Critical patent/CN114813828B/zh
Publication of CN114813828A publication Critical patent/CN114813828A/zh
Application granted granted Critical
Publication of CN114813828B publication Critical patent/CN114813828B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/20Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

发明公开了一种新型确定含水层热物性参数的微热试验方法,包括如下步骤:1)以热径向对流弥散理论研究为基础,依据热量守恒原理建立微热试验的理论模型,通过数学变换,建立理论模型不同渗流条件、不同参数条件下的多组标准曲线;2)通过微热试验,进行不同渗流条件和不同激发强度下的微热试验,通过传感器获取实测数据,绘制实测数据曲线;3)采用配线法将标准曲线与微热试验实测数据曲线进行匹配,记录对应的坐标值,通过计算得到含水层导热系数,通过相关公式推导其岩土体导热系数、含水层热弥散度等其余热物性参数。上述方法试验周期短,操作方法简单,不仅可以得到含水层的导热系数,也可以反推得到岩土体的导热系数。

Description

一种新型确定含水层热物性参数的微热试验方法
技术领域
发明涉及一种新型确定含水层热物性参数的微热试验方法,属于微热试验技术领域。
背景技术
开发利用绿色可再生的浅层地热能不仅可以有效缓解能源紧张问题,更能大大减少环境污染,是我国目前势在必行的研究项目。在开发地热能的初期,需要对地热资源进行相关评价,评价内容包括浅层地热资源总容量、可开采资源量以及地层热量补给速率等等,而这些均与区域水文地质条件、含水层和岩土体的热物性密切相关。热物性中最为重要的参数是含水层、岩土体导热系数和比热容,其大小反映了含水层、岩土体的蓄热、导热能力和传热效率,在工程应用中,选择较大导热系数的含水层和岩土体进行开发利用往往可以有效提高热源的利用率,降低投资成本。因此对含水层和岩土体导热系数进行充分的研究和评估是十分必要的。
此外,含水层和岩土体的传热问题在地热资源的开发与利用中是一个十分重要的研究方向,其中有效导热系数作为表征含水层热运移的重要参数有着举足轻重的地位,其意义在于将包括对流弥散、传导作用在内的热量运移过程用有效导热系数这样一个参数来表示。采用一定方法精准确定含水层和岩土体的导热系数具有重要的理论意义和实际工程价值。
目前测试导热系数的方法通常是现场热响应试验,试验的原理是通过给与含水层(岩土体)一定热源,观测热量的运移分布从而推导出含水层和岩土体的导热系数,该类试验周期长,操作较为复杂。
发明内容
本发明基于热径向对流弥散理论研究之上,提出了一种新型确定含水层热物性参数的微热试验方法,具有试验周期短,操作简易,数据处理直观等优点,具有很好的推广应用价值。
为解决上述技术问题,发明所采用的技术方案如下:
一种新型确定含水层热物性参数的微热试验方法,包括如下步骤:
1)以热径向对流弥散理论研究为基础,依据热量守恒原理建立微热试验的理论模型,通过数学变换,建立理论模型不同渗流条件、不同参数条件下的标准曲线;
2)通过微热试验,进行不同渗流条件和不同激发强度下的微热试验,通过传感器获取实测数据,绘制实测数据曲线;
3)采用配线法将标准曲线与微热试验实测数据曲线进行匹配,记录对应的坐标值,通过计算得到含水层导热系数,通过相关公式推导得到岩土体导热系数、含水层热弥散度等其余热物性参数。
上述利用微热试验获取实测数据,绘制实测数据曲线,具有试验周期短、试验简单,结果精确等特点。
上述步骤1)中,根据微热试验理论得到的半解析解绘制无量纲温度变化T*D和无量纲时间τ的标准曲线;步骤2)中,根据采集到的井孔内温度变化值T与时间t的观测数据,绘制与标准曲线相同模数的温度变化T与时间t的实测曲线图,温度变化值T和时间t均取对数;步骤3)中,使实测曲线坐标原点与标准曲线坐标原点在坐标图层的同一高度,通过平移实测曲线时间坐标轴t来匹配实测曲线与标准曲线,记录拟合标准曲线的α值,无量纲时间τ值,时间t值。
上述步骤1)中,微热试验的理论模型,建立在热径向对流弥散理论研究之上,依据热量守恒原理建立的一维条件下,不考虑z方向上热量运移的理论公式,通过拉普拉斯变化及逆变换对理论公式进行求解,选用Stehfest算法和Talbot算法建立理论模型不同渗流条件、不同参数条件下的标准曲线。
上述步骤1)中,综合利用Matlab和Grapher绘图软件,绘制不考虑流速和考虑流速两种情况下无量纲温度变化T*D与无量纲时间τ关于α(水热容量与含水层热容量比值)参量的标准曲线,即前述选用Stehfest算法建立的不考虑流速情况标准曲线和Talbot算法建立的考虑流速情况标准曲线。
上述步骤1)中,微热试验的理论模型需要通过以下方程和定律建立:径向流下考虑热机械弥散项的热量运移方程;根据热量守恒定律,井管内热量的变化与流过含水层r=rw处的热量相等,建立热量平衡方程;
径向流下,考虑热对流弥散项的热量运移方程:
Figure BDA0003614869240000021
根据热量守恒定律,井管内热量的变化与含水层r=rw处的热量相等:
Figure BDA0003614869240000022
初始条件
T(r,t)|t=0=0,r>rw (3)
T*(0)=T0 (4)
边界条件
T(r,t)|r→∞=0,t>0 (5)
T*(t)=T(rw,t),t>0 (6)
以上式中:T为温度,T(rw,t)为t时刻井管壁处的温度变化值,T(r,t)表示含水层中距井中心为r时t任意时刻的温度的变化值,T*(t)为t时刻井孔内的温度变化值,r为含水层内某点到井中心的距离,rw为井管半径,rs为井管内导热管的半径;t为时间;λ(r)为含水层的有效导热系数;(ρc)f为地下水热容量,(ρc)fs为含水层的热容量;T0为初始温度变化值;
Figure BDA0003614869240000031
为地下水实际流速,B为含水层厚度,Q为注水流量,θ为含水层孔隙度;
当u=0时,即含水层不考虑水流流动时,式(1)变成如下形式,即热的传导方程:
Figure BDA0003614869240000032
式(1)至式(6)数学模型的方程形式比较复杂,且因公式(1)和(7)表达式不同,需将公式(1)方程分成考虑流速和不考虑流速两种条件的情况,分别计算,总体思路是采用拉普拉斯变换的方法对其进行求解,将上述方程、定解条件等通过无量纲因子和参变量转换为无量纲形式,依次对其进行拉普拉斯变换;
(1)不考虑流速情况下的理论模型求解:
通过无量纲变换将式(7)转换为如下形式:
Figure BDA0003614869240000033
其中
Figure BDA0003614869240000034
δw=1,
Figure BDA0003614869240000035
再依次对式(2)-(6)做无量纲变化得:
Figure BDA0003614869240000036
TD(δ,0)=0,0<δ<∞,τ=0 (10)
T*D(0)=TD,δ=1,τ=0 (11)
TD(δ,τ)=0,δ→∞,τ>0 (12)
T*D(τ)=TD(δ,τ),δ=1,τ>0 (13)
最后对式(8)τ做拉普拉斯变换,式(8)热传导方程转化为:
Figure BDA0003614869240000041
井管内热量的变化与流过含水层r=rw处的热量相等,该式转化为
Figure BDA0003614869240000042
定解条件:
Figure BDA0003614869240000043
Figure BDA0003614869240000044
Figure BDA0003614869240000045
上述拉普拉斯转换式中:
Figure BDA0003614869240000046
为T*D的拉氏变换的象函数,
Figure BDA0003614869240000047
为TD的拉氏变换的象函数,T*D为井孔中温度变化比值,
Figure BDA0003614869240000048
为含水层中温度变化比值,
Figure BDA0003614869240000049
为距井中心距离与井半径的比,
Figure BDA00036148692400000410
井内导热管半径与井半径的比,δw=1,p为拉氏变换复变量,
Figure BDA00036148692400000411
无量纲时间,
Figure BDA00036148692400000412
Figure BDA00036148692400000413
为与含水层固有性质相关的无量纲参数,
Figure BDA00036148692400000414
为水的热容量与含水层的热容量的比。
根据n阶修正贝塞尔方程:
Figure BDA00036148692400000415
通解为:
y=C1Kn(βx)+C2In(βx) (20)
其中:In(x)和Kn(x)分别为第一类和第二类修正贝塞尔函数;
当n=0时,
式(14)的解为:
Figure BDA00036148692400000416
根据式(17)知:
Figure BDA00036148692400000417
又因为第一类修正贝塞尔函数I0(x)为指数增长型函数,具有性质:当x→∞时I0(x)→∞,因此C2=0,所以,
Figure BDA00036148692400000418
又根据0阶修正贝塞尔函数的性质知:
Figure BDA0003614869240000051
Figure BDA0003614869240000052
将式(22)代入式(15)和(18)中可得:
Figure BDA0003614869240000053
由(18)可得:
Figure BDA0003614869240000054
由此可以得到当含水层处于静止条件下,流速u=0时该理论模型在Laplace空间的解析解为:
Figure BDA0003614869240000055
Figure BDA0003614869240000056
公式(26)即为在不考虑流速情况下的拉普拉斯变换域中所得出的最终解,但是想要在实际实验中运用此解,还需要对公式(26)进行拉普拉斯逆变换,在拉普拉斯数值反演方法中,Stehfest(1970)方法以简单易行、计算时间短等优点被广泛应用:
Figure BDA0003614869240000057
其中:f(t)为实空间的像原函数;F(s)为拉氏空间的像函数;N为求和项数;
Figure BDA0003614869240000058
Figure BDA0003614869240000059
将公式(26)代入公式(27)中进行拉普拉斯逆变换得:
Figure BDA00036148692400000510
综上可得:
T*D=f(α,γ,τ) (29)
由此可知,在不考虑流速情况下微热试验的理论模型中,控制测试井中无量纲温度变化T*D的无因次参数共3个,其中
Figure BDA00036148692400000511
为水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,
Figure BDA00036148692400000512
为与含水层热容量相关的参数,δs反映了井管内导热管半径大小和井孔半径大小的比值,在假设条件中,将忽略导热管和井孔间的空隙,所以δs为1,τ为无因次时间,由此可利用Matlab软件绘制关于α参量的无量纲温度变化随时间变化标准曲线;
(2)考虑流速情况下的理论模型求解:
通过无量纲变换可将式(1)转换成如下形式:
Figure BDA0003614869240000061
对无量纲因子τ做Laplace变换可得:
Figure BDA0003614869240000062
式中:
Figure BDA0003614869240000063
β为与含水层导热系数和试验过程中试验孔流量相关的参数,若在流量和孔隙度等变量固定的条件下,含水层的导热系数λ为定值,β为常数,而在研究含水层热量运移问题的过程中,水流对温度场的影响不可忽略;
又因为对比n阶修正贝塞尔方程和式(31),不可以直接套用,需对式(31)做如下变换,为方便书写和理解,令
Figure BDA0003614869240000064
x=δ,则式(31)变换为:
xy″+βy′-pxy=0 (32)
设:y=x-m·u,
Figure BDA0003614869240000065
则:
y′=x-mu′-mx-m-1u (33)
y″=x-mu″-2mx-m-1u′+m(m+1)x-m-2u (34)
将式(33)和(34)代入(32)得:
Figure BDA0003614869240000066
根据n阶修正贝塞尔方程的解,可得:
Figure BDA0003614869240000067
Figure BDA0003614869240000068
由此可得,当n=0时,
Figure BDA0003614869240000069
同样根据定解条件以及第一类修正贝塞尔函数I0(x)具有如下性质:当x→∞时I0(x)→∞,可得到C2=0。
Figure BDA00036148692400000610
再将式(39)和式(18)代入无量纲化后的热量守恒条件即式(15),得:
Figure BDA0003614869240000071
Figure BDA0003614869240000072
由此可以得到当含水层在考虑流速情况下,微热试验理论模型在Laplace空间的解析解为:
Figure BDA0003614869240000073
Figure BDA0003614869240000074
和不考虑流速情况下求解过程相同,采用精度更高,考虑了复变域的Talbot方法进行Laplace逆变换:
Figure BDA0003614869240000075
其中,f(t)为实空间的像原函数;F(p)为拉氏空间的像函数,p=2N/5t;s(θ)=pθ(cotθ+i);
Figure BDA0003614869240000076
σ(θ)=θ+(θcotθ-1)cotθ
将公式(43)代入公式(44)中进行拉普拉斯逆变换,综上可得:
T*D=f(α,β,γ,τ) (45)
其中,
Figure BDA0003614869240000077
为水的热容量和含水层的热容量之比,
Figure BDA0003614869240000078
为与含水层导热系数和试验过程中试验孔流量相关的参数,
Figure BDA0003614869240000079
为与α相关的参数,
Figure BDA00036148692400000710
为无量纲时间。
由此可以得到在考虑流速情况下的微热试验理论中,控制测试井中无量纲温度变化的无因次参数共4个;与不考虑流速情况下一致,α指的是水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,β的大小与注入含水层内部的水流量以及含水层的性质相关,且反映了热量在含水层的运移过程中随着水流量的增加,对流弥散作用明显增强,远远大于分子扩散;
再次对β进行分析:
Figure BDA00036148692400000711
其中,
Figure BDA0003614869240000081
为热扩散率,反映了含水层的导热速率的快慢,
Figure BDA0003614869240000082
为水的热容量和含水层的热容量之比;根据β的表达式可以看出,β的大小与含水层的导热系数息息相关,在考虑流速情况下含水层的导热系数与流速、介质的孔隙度以及渗透系数有关,且0<β<1;γ与参数α和δs密切相关,除了反映含水层本身的热物性外,还反映了微热试验中导热管半径大小对试验结果的影响,τ为无因次时间;由此可固定其他参数值,利用Matlab软件绘制关于α、β参量的无量纲温度变化随时间变化的标准曲线;
热对流弥散模型中含水层的有效导热系数λ是由滞止导热系数λ*和热机械弥散系数λv两部分组成;由此,含水层中有效导热系数λ可表示为:
λ=λ*v (46)
含水层滞止导热系数λ*指的是地下水不流动时的含水层导热系数,属含水层(包含含水层骨架和地下水)自身的热物理参数,与地下水流速、流量等因素无关,即在不考虑流速情况下的所需求解的参数,此时热的传递方式只有传导作用,同样含水层的滞止导热系数可由下式计算:
Figure BDA0003614869240000083
式中:θ为含水层孔隙度,λs为岩土体导热系数,(ρc)fs=θ(ρc)f+(1-θ)(ρc)s,(ρc)f和(ρc)s分别为水和固体骨架的热容量;可以通过不考虑流速情况下的微热试验求出含水层的滞止导热系数和含水层的热容量,从而可以利用式(47)反推出岩土体的导热系数;
热机械弥散系数λv则与含水层的流速密不可分,在孔隙含水层热量运移相关研究中某些学者认为纵向热机械弥散系数λv与流速u成正比,因此可表示为:
λv=ax(ρc)fu (48)
式中:ax为纵向弥散度;u为地下水实际流速,(ρc)f为水的热容量,热弥散系数从定义上看指的是单位温度梯度下通过单位多孔介质截面的热量,其大小展现了多孔介质通过热弥散作用传递热量的能力强弱。
上述步骤2)中,微热试验是通过瞬间增加或减少井孔内的温度,记录井孔内温度随时间变化,研究其温度变化规律,进而求取含水层的热物性参数。
进一步优选,步骤2)中,在室内承压含水层试验平台或现场试验井孔中进行微热试验,根据微热试验方法,先后获得两种不同流速情况下、温差比
Figure BDA0003614869240000084
和试验时间t试验数据,在Grapher绘图软件中绘制
Figure BDA0003614869240000091
关系实测曲线。
上述步骤3)中,求取含水层导热系数并通过公式推导其余热物性参数的步骤如下:①在Grapher软件中将不考虑流速和考虑流速情况下的实测曲线分别和相应条件下的标准曲线对应,置于同样模数的双对数坐标下,固定纵坐标在同一高度对横坐标进行平移,直至两曲线重合;②任取一匹配点,记录下匹配点的对应坐标值:[t],[τ],[α],其中,考虑流速情况下,还需要记录下对应标准曲线的值:[β]。已知水的热容量(ρc)f与井管半径rw,根据不考虑流速情况下标准曲线依据
Figure BDA0003614869240000092
求出含水层的滞止导热系数
Figure BDA0003614869240000093
根据考虑流速情况下标准曲线依据
Figure BDA0003614869240000094
可以求出含水层的有效导热系数λ,依据λv=λ-λ*可以求出含水层的热机械弥散系数λv,依据
Figure BDA0003614869240000095
求出热弥散度ax。同时,依据
Figure BDA0003614869240000096
求出含水层的热容量(ρc)fs。此外,根据室内土工试验资料可已知孔隙度θ,同时,在求出含水层的热容量(ρc)fs基础上,根据(ρc)fs=θ(ρc)f+(1-θ)(ρc)s可求出岩土体固体骨架热容量(ρc)s,依据
Figure BDA0003614869240000097
求出岩土体导热系数λs
上述方法,具有试验周期短、操作简单的优点;理论模型通过严格的数学推导得到,理论严谨,采用配线法确定含水层热物性参数,方法简单易掌握。
本发明未提及的技术均参照现有技术。
有益效果:与现有方法相比,本发明所提供的一种新型确定含水层热物性参数的微热试验方法,具有以下优点:
(1)本发明提出的微热试验理论公式基于热径向对流弥散理论研究之上,通过严格的数学推导,理论严谨,并结合室内试验进行了验证;
(2)热物性参数的确定采用配线法,方法简单,易于掌握;
(3)相较于传统的确定热物性参数的方法,微热试验试验周期短,操作方法简单,不仅可以得到含水层的有效导热系数、滞止导热系数、热机械弥散系数、热弥散度和热容量,也可以反推得到岩土体的导热系数,且准确性高,具有很好的推广价值。
附图说明
图1为本发明方法的操作流程图;
图2为根据微热试验理论,不考虑流速情况下,采用Stehfest算法得到的无量纲温度变化T*D和无量纲时间τ的标准曲线,即T*D-τ关系曲线;
图3为根据微热试验理论,在考虑流速情况下,采用Talbot算法得到的无量纲温度变化T*D和无量纲时间τ的标准曲线,即T*D-τ关系曲线;
图4为室内半圆柱形砂槽承压含水层试验平台主视图;
图5为室内半圆柱形砂槽承压含水层试验平台平面图;
图6为室内半圆柱形砂槽承压含水层试验平台AA’方向剖面图;
图7为室内半圆柱形砂槽承压含水层试验平台整体框架外观;
图8为微热试验不考虑流速情况下,试验编号为w65℃1时主孔试验数据采用微热试验方法进行线状热源激发获取主孔温度变化数据确定含水层热物性参数时的配线图;
图9为微热实验不考虑流速情况下,试验编号为w75℃1时主孔试验数据采用微热试验方法进行线状热源激发获取主孔温度变化数据确定含水层热物性参数时的配线图;
图10为微热实验不考虑流速情况下,试验编号为w85℃1时主孔试验数据采用微热试验方法进行线状热源激发获取主孔温度变化数据确定含水层热物性参数时的配线图;
图11为微热试验在考虑流速情况下,试验编号为y(△h=6)45℃1主孔试验数据采用微热试验方法进行线状热源激发获取主孔温度变化数据确定含水层热物性参数时的配线图;
图12为微热试验在考虑流速情况下,试验编号为y(△h=9)65℃1时主孔试验数据采用微热试验方法进行线状热源激发获取主孔温度变化数据确定含水层热物性参数时的配线图;
图13为微热试验在考虑流速情况下,试验编号为y(△h=12)85℃2时主孔试验数据采用微热试验方法进行线状热源激发获取主孔温度变化数据确定含水层热物性参数时的配线图;
图中,1为试验平台主井孔,2为径向布置的1-5号温度传感器,3为承压含水层,4为试验平台过滤层,5为试验平台外侧挡板,6为试验平台排水口,7为试验平台水位观测孔,8为试验平台9根测压管,9为计算机数据收集系统,10为隔水顶板与隔热层,11为沙袋,12为微热试验激发装置,13为隔水底板与隔热层,14为试验平台溢水口。
具体实施方式
为了更好地理解发明,下面结合实施例进一步阐明发明的内容,但发明的内容不仅仅局限于下面的实施例。
利用微热试验确定含水层热物性参数,首先在室内进行微热试验,如图4-图7所示,采用选定的多孔砂样介质装填的半圆柱形砂槽模拟承压含水层系统,系统由模型框架、承压含水层、温度观测装置构成。模型框架整体设计为高2.2m、半径2m的半圆柱状钢架,侧面和底面设置支撑架,预留排水孔和各观测孔外伸接口,在半圆柱圆弧侧边距底板1.5m处预留四个溢水孔以满足定水头要求。承压含水层模拟承压环境,厚0.75m,介质为0.25mm-0.50mm的中细砂,承压含水层的顶、底板铺设塑料膜和隔热膜进行防水隔热处理,顶板上铺设砂袋施加压力,承压含水层由侧边进水,底部排水。温度观测装置预置在主孔及主孔周围一定距离处,通过电缆与计算机相连进行数据传输。微热试验激发装置的结构及使用,参照本申请人同期申请的专利申请:一种确定含水层热物性参数的微热试验装置及方法。
如图1所示,一种新型确定含水层热物性参数的微热试验方法,包括如下步骤:
1)以热径向对流弥散理论研究为基础,依据热量守恒原理建立微热试验的理论模型,通过数学变换,建立理论模型不同渗流条件、不同参数条件下的标准曲线;
2)通过微热试验,在上述半圆柱形砂槽模拟承压含水层系统中进行不同渗流条件和不同激发强度下的微热试验,通过传感器获取实测数据,绘制实测数据曲线;
3)采用配线法将标准曲线与微热试验实测数据曲线进行匹配,记录对应的坐标值,通过计算得到含水层导热系数,通过相关公式推导得到岩土体导热系数、含水层热弥散度等其余热物性参数。
本次室内试验研究共分为:不考虑流速情况下、流速为9.3×10-3cm/s、6.82×10- 3cm/s、4.34×10-3cm/s四种流速条件,每种流速条件下又进行三种不同温度激发强度下的微热试验。注入热水量根据导热管的半径以及含水层厚度计算而来,约800ml。以主孔探头数据作为计算参数的数据,每组试验进行两组平行试验,共计24组。
根据微热试验计算含水层热物性参数的步骤和方法,首先利用Grapher绘图软件绘制T*D-τ关系的标准曲线,不考虑流速情况下的T*D-τ关系曲线如图2所示,在考虑流速情况下,即T*D-τ关系曲线如图3所示,然后根据试验数据在相同模的透明双对数纸上作
Figure BDA0003614869240000111
实测曲线,使实测曲线坐标原点与标准曲线坐标原点保持同一高度,通过平移实测曲线时间坐标轴t来匹配实测曲线与标准曲线;任取一匹配点,记下匹配点的对应坐标[τ]值和[t]值,记录拟合标准曲线的[α]值;根据公式计算出含水层滞止导热系数与含水层有效导热系数,并可以反推出含水层热机械弥散系数与含水层热弥散度。其中不考虑流速情况下参数计算的配线图如图8至图10所示,滞止导热系数、含水层热容量与根据式(47)反推的岩土体导热系数计算结果见表1和表2。考虑流速情况下参数计算的配线图如图11至图13所示,含水层有效导热系数、含水层热机械弥散系数与根据式(48)反推的含水层热弥散度计算结果见表3和表4。
表1不考虑流速情况下的微热试验模型计算含水层导热系数的结果
Figure BDA0003614869240000121
表2岩土体导热系数计算结果
Figure BDA0003614869240000122
表2中,水的导热系数和热容量根据经验值取得,含水层的热容量根据配线数据α求得,依据(ρc)fs=θ(ρc)f+(1-θ)(ρc)s求得岩土体的热容量(ρc)s,最后依据
Figure BDA0003614869240000123
Figure BDA0003614869240000124
求出岩土体导热系数λs
从表1和表2可以看出,本次室内模拟试验的含水层孔隙度为θ=0.33,(ρc)f=4.18×106J/(m3·℃),计算求得岩土体导热系数的平均值为4.465W/(m·℃)。本次试验采用的是粒径在0.25mm-0.50mm的中细砂。根据查找相关手册和文献,中细砂的导热系数一般在1.7~5.0W/(m·℃)范围内,本次试验求取的岩土体导热系数刚好在此范围内,由此也可证明微热试验确定含水层导热系数以及岩土体导热系数的准确性和可行性。
表3考虑流速情况下的微热试验模型计算含水层导热系数的结果
Figure BDA0003614869240000131
表4含水层热弥散度计算结果
Figure BDA0003614869240000132
从表3和表4可以看出,根据不考虑流速计算得到的滞止导热系数,可以反推出热机械弥散系数与热弥散度。通过微热试验理论进行室内微热试验,操作简便,计算结果精确,试验周期短,能够较好的求出含水层的热物性参数。

Claims (8)

1.一种新型确定含水层热物性参数的微热试验方法,其特征在于:包括如下步骤:
1)以热径向对流弥散理论为基础,依据热量守恒原理建立微热试验的理论模型,通过数学变换,建立理论模型不同渗流条件、不同参数条件下的标准曲线;
2)通过微热试验,进行不同渗流条件和不同激发强度下的微热试验,通过传感器获取实测数据,绘制实测数据曲线;
3)采用配线法将标准曲线与微热试验实测数据曲线进行匹配,记录对应的坐标值,通过计算得到含水层导热系数,通过相关公式推导得到其余热物性参数。
2.如权利要求1所述的新型确定含水层热物性参数的试验方法,其特征在于:步骤1)中,根据微热试验理论得到的半解析解绘制无量纲温度变化T*D和无量纲时间τ的标准曲线;步骤2)中,根据采集到的井孔内温度变化值T与时间t的观测数据,绘制与标准曲线相同模数的温度变化T与时间t的实测曲线图,温度变化值T和时间t均取对数;步骤3)中,使实测曲线坐标原点与标准曲线坐标原点在坐标图层的同一高度,通过平移实测曲线时间坐标轴t来匹配实测曲线与标准曲线,记录拟合标准曲线的α值、无量纲时间τ值和时间t值。
3.如权利要求1或2所述的新型确定含水层热物性参数的试验方法,其特征在于:步骤1)中,微热试验的理论模型,建立在热径向对流弥散理论研究之上,依据热量守恒原理建立的一维条件下,不考虑z方向上热量运移的理论公式,通过拉普拉斯变化及逆变换对理论公式进行求解,选用Stehfest算法和Talbot算法建立理论模型不同渗流条件、不同参数条件下的标准曲线。
4.如权利要求3所述的新型确定含水层热物性参数的试验方法,其特征在于:步骤1)中,综合利用Matlab和Grapher绘图软件,绘制不考虑流速和考虑流速两种条件下无量纲温度变化T*D与无量纲时间τ关于α参量的标准曲线,其中,α为水热容量与含水层热容量比值。
5.如权利要求4所述的新型确定含水层热物性参数的试验方法,其特征在于:步骤1)中,微热试验的理论模型需要通过以下方程和定律建立:径向流下考虑热机械弥散项的热量运移方程;根据热量守恒定律,井管内热量的变化与流过含水层r=rw处的热量相等,建立热量平衡方程;
径向流下,考虑热对流弥散项的热量运移方程:
Figure FDA0003614869230000011
根据热量守恒定律,井管内热量的变化与流过含水层r=rw处的热量相等:
Figure FDA0003614869230000021
初始条件
T(r,t)|t=0=0,r>rw (3)
T*(0)=T0 (4)
边界条件
T(r,t)|r→∞=0,t>0 (5)
T*(t)=T(rw,t),t>0 (6)
以上式中:T为温度,T(rw,t)为t时刻井管壁处的温度变化值,T(r,t)表示含水层中距井中心为r时t任意时刻的温度的变化值,T*(t)为t时刻井孔内的温度变化值,r为含水层内某点到井中心的距离,rw为井管半径,rs为井管内导热管的半径;t为时间;λ(r)为含水层的有效导热系数;(ρc)f为地下水热容量,(ρc)fs为含水层的热容量;T0为初始温度变化值;
Figure FDA0003614869230000022
为地下水实际流速,B为含水层厚度,Q为注水流量,θ为含水层孔隙度;
当u=0时,即含水层不考虑流速时,式(1)变成如下形式,即热的传导方程:
Figure FDA0003614869230000023
式(1)至式(6)数学模型的方程形式比较复杂,且因公式(1)和(7)表达式不同,需将公式(1)方程分成考虑流速和不考虑流速两种条件的情况,分别计算,总体思路是采用拉普拉斯变换的方法对其进行求解,将上述方程、定解条件等通过无量纲因子和参变量转换为无量纲形式,依次对其进行拉普拉斯变换;
(1)不考虑流速情况下的理论模型求解:
通过无量纲变换将式(7)转换为如下形式::
Figure FDA0003614869230000024
其中
Figure FDA0003614869230000025
δw=1,
Figure FDA0003614869230000026
再依次对式(2)-(6)做无量纲变化得:
Figure FDA0003614869230000027
TD(δ,0)=0,0<δ<∞,τ=0 (10)
T*D(0)=TD,δ=1,τ=0 (11)
TD(δ,τ)=0,δ→∞,τ>0 (12)
T*D(τ)=TD(δ,τ),δ=1,τ>0 (13)
最后对式(8)τ做拉普拉斯变换,式(8)热传导方程转化为:
Figure FDA0003614869230000031
井管内热量的变化与流过含水层r=rw处的热量相等,该式转化为
Figure FDA0003614869230000032
定解条件:
Figure FDA0003614869230000033
Figure FDA0003614869230000034
Figure FDA0003614869230000035
上述拉普拉斯转换式中:
Figure FDA0003614869230000036
为T*D的拉氏变换的象函数,
Figure FDA0003614869230000037
为TD的拉氏变换的象函数,T*D为井孔中温度变化比值,
Figure FDA0003614869230000038
为含水层中温度变化比值,
Figure FDA0003614869230000039
为距井中心距离与井半径的比,
Figure FDA00036148692300000310
井内导热管半径与井半径的比,δw=1,p为拉氏变换复变量,
Figure FDA00036148692300000311
无量纲时间,
Figure FDA00036148692300000312
Figure FDA00036148692300000313
为与含水层固有性质相关的无量纲参数,
Figure FDA00036148692300000314
为水的热容量与含水层的热容量的比;
根据n阶修正贝塞尔方程:
Figure FDA00036148692300000315
通解为:
y=C1Kn(βx)+C2In(βx) (20)
其中:In(x)和Kn(x)分别为第一类和第二类修正贝塞尔函数;
当n=0时,
式(14)的解为:
Figure FDA00036148692300000316
根据式(17)知:
Figure FDA00036148692300000317
又因为第一类修正贝塞尔函数I0(x)为指数增长型函数,具有性质:当x→∞时I0(x)→∞,因此C2=0,所以,
Figure FDA0003614869230000041
又根据0阶修正贝塞尔函数的性质知:
Figure FDA0003614869230000042
Figure FDA0003614869230000043
将式(22)代入式(15)和(18)中可得:
Figure FDA0003614869230000044
由(18)可得:
Figure FDA0003614869230000045
由此可以得到当含水层处于静止条件下,流速u=0时该理论模型在Laplace空间的解析解为:
Figure FDA0003614869230000046
Figure FDA0003614869230000047
公式(26)即为在不考虑流速条件下的拉普拉斯变换域中所得出的最终解,但是想要在实际实验中运用此解,还需要对公式(26)进行拉普拉斯逆变换:
Figure FDA0003614869230000048
其中:f(t)为实空间的像原函数;F(s)为拉氏空间的像函数;N为求和项数;
Figure FDA0003614869230000049
Figure FDA00036148692300000410
将公式(26)代入公式(27)中进行拉普拉斯逆变换得:
Figure FDA00036148692300000411
综上可得:
T*D=f(α,γ,τ) (29)
由此可知,在不考虑流速情况下微热试验的理论模型中,控制测试井中无量纲温度变化T*D的无因次参数共3个,其中
Figure FDA0003614869230000051
为水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,
Figure FDA0003614869230000052
为与含水层热容量相关的参数,δs反映了井管内导热管半径大小和井孔半径大小的比值,在假设条件中,将忽略导热管和井孔间的空隙,所以δs为1,τ为无因次时间,由此可利用Matlab软件绘制关于α参量的无量纲温度变化随时间变化标准曲线;
(2)考虑流速情况下的理论模型求解:
通过无量纲变换可将式(1)转换成如下形式:
Figure FDA0003614869230000053
对无量纲因子τ做Laplace变换可得:
Figure FDA0003614869230000054
式中:
Figure FDA0003614869230000055
β为与含水层导热系数和试验过程中试验孔流量相关的参数,若在流量和孔隙度等变量固定的条件下,含水层的导热系数λ为定值,β为常数,而在研究含水层热量运移问题的过程中,水流对温度场的影响不可忽略;
又因为对比n阶修正贝塞尔方程和式(31),不可以直接套用,需对式(31)做如下变换,为方便书写和理解,令
Figure FDA0003614869230000056
x=δ,则式(31)变换为:
xy″+βy′-pxy=0 (32)
设:y=x-m·u,
Figure FDA0003614869230000057
则:
y′=x-mu′-mx-m-1u (33)
y″=x-mu″-2mx-m-1u′+m(m+1)x-m-2u (34)
将式(33)和(34)代入(32)得:
Figure FDA0003614869230000058
根据n阶修正贝塞尔方程的解,可得:
Figure FDA0003614869230000059
Figure FDA00036148692300000510
由此可得,当n=0时,
Figure FDA00036148692300000511
同样根据定解条件以及第一类修正贝塞尔函数I0(x)具有如下性质:当x→∞时I0(x)→∞,可得到C2=0;
Figure FDA0003614869230000061
再将式(39)和式(18)代入无量纲化后的热量守恒条件即式(15),得:
Figure FDA0003614869230000062
Figure FDA0003614869230000063
由此可以得到当含水层在考虑流速情况下,微热试验理论模型在Laplace空间的解析解为:
Figure FDA0003614869230000064
Figure FDA0003614869230000065
和不考虑流速情况下求解过程相同,采用精度更高,考虑了复变域的Talbot方法进行Laplace逆变换:
Figure FDA0003614869230000066
其中,f(t)为实空间的像原函数;F(p)为拉氏空间的像函数,p=2N/5t;s(θ)=pθ(cotθ+i);
Figure FDA0003614869230000067
σ(θ)=θ+(θcotθ-1)cotθ
将公式(43)代入公式(44)中进行拉普拉斯逆变换,综上可得:
T*D=f(α,β,γ,τ) (45)
其中,
Figure FDA0003614869230000068
为水的热容量和含水层的热容量之比,
Figure FDA0003614869230000069
为与含水层导热系数和试验过程中试验孔流量相关的参数,
Figure FDA00036148692300000610
为与α相关的参数,
Figure FDA00036148692300000611
为无量纲时间;
由此可以得到在考虑流速情况下的微热试验理论中,控制测试井中无量纲温度变化的无因次参数共4个;与不考虑流速情况下一致,α指的是水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,β的大小与注入含水层内部的水流量以及含水层的性质相关,且反映了热量在含水层的运移过程中随着水流量的增加,对流弥散作用明显增强,远远大于分子扩散;
再次对β进行分析:
Figure FDA0003614869230000071
其中,
Figure FDA0003614869230000072
为热扩散率,反映了含水层的导热速率的快慢,
Figure FDA0003614869230000073
为水的热容量和含水层的热容量之比;根据β的表达式可以看出,β的大小与含水层的导热系数息息相关,在考虑流速情况下含水层的导热系数与流速、介质的孔隙度以及渗透系数有关,且0<β<1;γ与参数α和δs密切相关,除了反映含水层本身的热物性外,还反映了微热试验中导热管半径大小对试验结果的影响,τ为无因次时间;由此可固定其他参数值,利用Matlab软件绘制关于α、β参量的无量纲温度变化随时间变化的标准曲线;
热对流弥散模型中含水层的有效导热系数λ是由滞止导热系数λ*和热机械弥散系数λv两部分组成;由此,含水层中有效导热系数系数λ可表示为:
λ=λ*v (46)
含水层滞止导热系数λ*指的是地下水不流动时的含水层导热系数,属含水层(包含含水层骨架和地下水)自身的热物理参数,与地下水流速、流量等因素无关,即在不考虑流速情况下的所需求解的参数,此时热的传递方式只有传导作用,同样含水层的滞止导热系数可由下式计算:
Figure FDA0003614869230000074
式中:θ为含水层孔隙度,λs为岩土体导热系数,(ρc)fs=θ(ρc)f+(1-θ)(ρc)s,(ρc)f和(ρc)s分别为水和固体骨架的热容量;可以通过不考虑流速情况下的微热试验求出含水层的滞止导热系数和含水层的热容量,从而可以利用式(47)反推出岩土体的导热系数;
热机械弥散系数λv则与含水层的流速密不可分,在孔隙含水层热量运移相关研究中某些学者认为纵向热机械弥散系数λv与流速u成正比,因此可表示为:
λv=ax(ρc)fu (48)
式中:ax为纵向弥散度;u为地下水实际流速,(ρc)f为水的热容量,热弥散系数从定义上看指的是单位温度梯度下通过单位多孔介质截面的热量,其大小展现了多孔介质通过热弥散作用传递热量的能力强弱。
6.如权利要求1或2所述的新型确定含水层热物性参数的微热试验方法,其特征在于:步骤2)中,微热试验是通过瞬间增加或减少井孔内的温度,记录井孔内温度随时间变化,研究其温度变化规律,进而求取含水层的热物性参数。
7.如权利要求6所述的新型确定含水层热物性参数的试验方法,其特征在于:步骤2)中,在室内承压含水层试验平台或现场试验井孔中进行微热试验,根据微热试验方法,先后获得两种不同流速情况下、温差比
Figure FDA0003614869230000081
和试验时间t试验数据,在Grapher绘图软件中绘制
Figure FDA0003614869230000082
关系实测曲线。
8.如权利要求1或2所述的新型确定含水层热物性参数的微热试验方法,其特征在于:步骤3)中,求取含水层导热系数并通过公式推导其余热物性参数的步骤如下:①在Grapher软件中将不考虑流速和考虑流速情况下的实测曲线分别和相应条件下的标准曲线对应,置于同样模数的双对数坐标下,固定纵坐标在同一高度对横坐标进行平移,直至两曲线重合;②任取一匹配点,记录下匹配点的对应坐标值:[t],[τ],[α],其中,考虑流速情况下,还需要记录对应标准曲线的值:[β];已知水的热容量(ρc)f与井管半径rw,根据不考虑流速情况下标准曲线依据
Figure FDA0003614869230000083
求出含水层的滞止导热系数λ*;根据考虑流速情况下标准曲线依据
Figure FDA0003614869230000084
可以求出含水层的有效导热系数λ,依据λv=λ-λ*可以求出含水层的热机械弥散系数λv,依据
Figure FDA0003614869230000085
求出热弥散度ax;同时,依据
Figure FDA0003614869230000086
求出含水层的热容量(ρc)fs;此外,根据室内土工试验资料可已知孔隙度θ,同时,在求出含水层的热容量(ρc)fs基础上,根据(ρc)fs=θ(ρc)f+(1-θ)(ρc)s可求出岩土体固体骨架热容量(ρc)s,依据
Figure FDA0003614869230000087
求出岩土体导热系数λs
CN202210440391.3A 2022-04-25 2022-04-25 一种确定含水层热物性参数的微热试验方法 Active CN114813828B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210440391.3A CN114813828B (zh) 2022-04-25 2022-04-25 一种确定含水层热物性参数的微热试验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210440391.3A CN114813828B (zh) 2022-04-25 2022-04-25 一种确定含水层热物性参数的微热试验方法

Publications (2)

Publication Number Publication Date
CN114813828A true CN114813828A (zh) 2022-07-29
CN114813828B CN114813828B (zh) 2022-09-20

Family

ID=82507056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210440391.3A Active CN114813828B (zh) 2022-04-25 2022-04-25 一种确定含水层热物性参数的微热试验方法

Country Status (1)

Country Link
CN (1) CN114813828B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115408846A (zh) * 2022-08-25 2022-11-29 中国科学院武汉岩土力学研究所 珊瑚细砂热物理参数简化预测模型及其试验装置和方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007263957A (ja) * 2006-03-01 2007-10-11 Nippon Steel Engineering Co Ltd 地下水流れ状態推定装置及び方法、土壌の有効熱伝導率推定方法、並びに地盤調査方法
CN102323294A (zh) * 2011-08-05 2012-01-18 江亚斌 一种岩土热响应测试方法
US20130211771A1 (en) * 2012-02-13 2013-08-15 Kambix Innovations, LLC. Manipulating heat flux bifurcation & dispersion inside porous media for heat transfer control
CN106248546A (zh) * 2016-06-16 2016-12-21 水利部交通运输部国家能源局南京水利科学研究院 一种多尺度热量运移同步监测试验系统和试验方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007263957A (ja) * 2006-03-01 2007-10-11 Nippon Steel Engineering Co Ltd 地下水流れ状態推定装置及び方法、土壌の有効熱伝導率推定方法、並びに地盤調査方法
CN102323294A (zh) * 2011-08-05 2012-01-18 江亚斌 一种岩土热响应测试方法
US20130211771A1 (en) * 2012-02-13 2013-08-15 Kambix Innovations, LLC. Manipulating heat flux bifurcation & dispersion inside porous media for heat transfer control
CN106248546A (zh) * 2016-06-16 2016-12-21 水利部交通运输部国家能源局南京水利科学研究院 一种多尺度热量运移同步监测试验系统和试验方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHOU YAN-ZHANG 等: "Simulation of Thermal Transport in Aquifer: A GWHP System in Chengdu, China", 《JOURNAL OF HYDRODYNAMICS, SER. B》 *
刘国庆等: "小尺度含水层热量运移试验研究及热弥散效应评估", 《岩土力学》 *
吕培等: "分层微水试验与数值模拟的正反演验证对比分析", 《勘察科学技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115408846A (zh) * 2022-08-25 2022-11-29 中国科学院武汉岩土力学研究所 珊瑚细砂热物理参数简化预测模型及其试验装置和方法
CN115408846B (zh) * 2022-08-25 2023-08-15 中国科学院武汉岩土力学研究所 珊瑚细砂热物理参数简化预测模型及其试验装置和方法

Also Published As

Publication number Publication date
CN114813828B (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
Bauer et al. Transient 3D analysis of borehole heat exchanger modeling
Luo et al. Analysis on performance of borehole heat exchanger in a layered subsurface
Signorelli et al. Numerical evaluation of thermal response tests
Stähli et al. Soil moisture redistribution and infiltration in frozen sandy soils
Raymond et al. Numerical analysis of thermal response tests with a groundwater flow and heat transfer model
Zhang et al. Comparison of test methods for shallow layered rock thermal conductivity between in situ distributed thermal response tests and laboratory test based on drilling in northeast China
Durner et al. Determining soil hydraulic properties
Ren et al. A review on using heat as a tool for studying groundwater–surface water interactions
CN114813828B (zh) 一种确定含水层热物性参数的微热试验方法
Tan et al. An innovative method to evaluate hydraulic conductivity of a single rock fracture based on geometric characteristics
Yang et al. Optimal control of salinity boundary condition in a tidal model using a variational inverse method
Zhou et al. An analytical model for heat transfer outside a single borehole heat exchanger considering convection at ground surface and advection of vertical water flow
Song et al. Seepage monitoring of an embankment dam based on hydro-thermal coupled analysis
Liu et al. Development and validation of a method to predict the soil thermal conductivity using thermal piezocone penetration testing (T-CPTU)
Kang et al. Feasibility analysis of extreme learning machine for predicting thermal conductivity of rocks
Ren et al. A comparison of numerical and Lu modeling of water flow and heat transport with laboratory experiments
Courtois et al. Estimating thermal properties of phase change material from heat flux measurements
Wang et al. Comparative study on effects of macroscopic and microscopic fracture structures on the performance of enhanced geothermal systems
Krisdani et al. Use of instantaneous profile and statistical methods to determine permeability functions of unsaturated soils
CN105572161B (zh) 一种非恒态导热系数测试的方法及所用测试装置
CN110728064B (zh) 基于comsol数学模块的河岸潜流带水热耦合建模方法
Zhang et al. Experimental study on the effect of initial water content and temperature gradient on soil column segregation frost heave
Zhang et al. Comparison of 1-D analytical solutions and a numerical model for quantifying hyporheic exchange flux using the temperature tracer method
Liu et al. Analysis and experimental study on thermal dispersion effect of small scale saturated porous aquifer
CN106153293A (zh) 一种基于温度示踪的潜流通量的测量装置及测量方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20220729

Assignee: Yangzhou Runbai Engineering Construction Co.,Ltd.

Assignor: HOHAI University

Contract record no.: X2022320000310

Denomination of invention: A Microcalorimetric Test Method for Determining the Thermophysical Parameters of Aquifers

Granted publication date: 20220920

License type: Common License

Record date: 20221210

EE01 Entry into force of recordation of patent licensing contract