CN114813828A - 一种新型确定含水层热物性参数的微热试验方法 - Google Patents
一种新型确定含水层热物性参数的微热试验方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating 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处的热量相等,建立热量平衡方程;
径向流下,考虑热对流弥散项的热量运移方程:
根据热量守恒定律,井管内热量的变化与含水层r=rw处的热量相等:
初始条件
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为初始温度变化值;为地下水实际流速,B为含水层厚度,Q为注水流量,θ为含水层孔隙度;
当u=0时,即含水层不考虑水流流动时,式(1)变成如下形式,即热的传导方程:
式(1)至式(6)数学模型的方程形式比较复杂,且因公式(1)和(7)表达式不同,需将公式(1)方程分成考虑流速和不考虑流速两种条件的情况,分别计算,总体思路是采用拉普拉斯变换的方法对其进行求解,将上述方程、定解条件等通过无量纲因子和参变量转换为无量纲形式,依次对其进行拉普拉斯变换;
(1)不考虑流速情况下的理论模型求解:
通过无量纲变换将式(7)转换为如下形式:
再依次对式(2)-(6)做无量纲变化得:
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)热传导方程转化为:
井管内热量的变化与流过含水层r=rw处的热量相等,该式转化为
定解条件:
上述拉普拉斯转换式中:
为T*D的拉氏变换的象函数,为TD的拉氏变换的象函数,T*D为井孔中温度变化比值,为含水层中温度变化比值,为距井中心距离与井半径的比,井内导热管半径与井半径的比,δw=1,p为拉氏变换复变量,无量纲时间, 为与含水层固有性质相关的无量纲参数,为水的热容量与含水层的热容量的比。
根据n阶修正贝塞尔方程:
通解为:
y=C1Kn(βx)+C2In(βx) (20)
其中:In(x)和Kn(x)分别为第一类和第二类修正贝塞尔函数;
当n=0时,
式(14)的解为:
根据式(17)知:
又因为第一类修正贝塞尔函数I0(x)为指数增长型函数,具有性质:当x→∞时I0(x)→∞,因此C2=0,所以,
又根据0阶修正贝塞尔函数的性质知:
将式(22)代入式(15)和(18)中可得:
由(18)可得:
由此可以得到当含水层处于静止条件下,流速u=0时该理论模型在Laplace空间的解析解为:
公式(26)即为在不考虑流速情况下的拉普拉斯变换域中所得出的最终解,但是想要在实际实验中运用此解,还需要对公式(26)进行拉普拉斯逆变换,在拉普拉斯数值反演方法中,Stehfest(1970)方法以简单易行、计算时间短等优点被广泛应用:
将公式(26)代入公式(27)中进行拉普拉斯逆变换得:
综上可得:
T*D=f(α,γ,τ) (29)
由此可知,在不考虑流速情况下微热试验的理论模型中,控制测试井中无量纲温度变化T*D的无因次参数共3个,其中为水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,为与含水层热容量相关的参数,δs反映了井管内导热管半径大小和井孔半径大小的比值,在假设条件中,将忽略导热管和井孔间的空隙,所以δs为1,τ为无因次时间,由此可利用Matlab软件绘制关于α参量的无量纲温度变化随时间变化标准曲线;
(2)考虑流速情况下的理论模型求解:
通过无量纲变换可将式(1)转换成如下形式:
对无量纲因子τ做Laplace变换可得:
xy″+βy′-pxy=0 (32)
y′=x-mu′-mx-m-1u (33)
y″=x-mu″-2mx-m-1u′+m(m+1)x-m-2u (34)
将式(33)和(34)代入(32)得:
根据n阶修正贝塞尔方程的解,可得:
由此可得,当n=0时,
同样根据定解条件以及第一类修正贝塞尔函数I0(x)具有如下性质:当x→∞时I0(x)→∞,可得到C2=0。
再将式(39)和式(18)代入无量纲化后的热量守恒条件即式(15),得:
由此可以得到当含水层在考虑流速情况下,微热试验理论模型在Laplace空间的解析解为:
和不考虑流速情况下求解过程相同,采用精度更高,考虑了复变域的Talbot方法进行Laplace逆变换:
将公式(43)代入公式(44)中进行拉普拉斯逆变换,综上可得:
T*D=f(α,β,γ,τ) (45)
由此可以得到在考虑流速情况下的微热试验理论中,控制测试井中无量纲温度变化的无因次参数共4个;与不考虑流速情况下一致,α指的是水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,β的大小与注入含水层内部的水流量以及含水层的性质相关,且反映了热量在含水层的运移过程中随着水流量的增加,对流弥散作用明显增强,远远大于分子扩散;
再次对β进行分析:
其中,为热扩散率,反映了含水层的导热速率的快慢,为水的热容量和含水层的热容量之比;根据β的表达式可以看出,β的大小与含水层的导热系数息息相关,在考虑流速情况下含水层的导热系数与流速、介质的孔隙度以及渗透系数有关,且0<β<1;γ与参数α和δs密切相关,除了反映含水层本身的热物性外,还反映了微热试验中导热管半径大小对试验结果的影响,τ为无因次时间;由此可固定其他参数值,利用Matlab软件绘制关于α、β参量的无量纲温度变化随时间变化的标准曲线;
热对流弥散模型中含水层的有效导热系数λ是由滞止导热系数λ*和热机械弥散系数λv两部分组成;由此,含水层中有效导热系数λ可表示为:
λ=λ*+λv (46)
含水层滞止导热系数λ*指的是地下水不流动时的含水层导热系数,属含水层(包含含水层骨架和地下水)自身的热物理参数,与地下水流速、流量等因素无关,即在不考虑流速情况下的所需求解的参数,此时热的传递方式只有传导作用,同样含水层的滞止导热系数可由下式计算:
式中:θ为含水层孔隙度,λ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)中,微热试验是通过瞬间增加或减少井孔内的温度,记录井孔内温度随时间变化,研究其温度变化规律,进而求取含水层的热物性参数。
上述步骤3)中,求取含水层导热系数并通过公式推导其余热物性参数的步骤如下:①在Grapher软件中将不考虑流速和考虑流速情况下的实测曲线分别和相应条件下的标准曲线对应,置于同样模数的双对数坐标下,固定纵坐标在同一高度对横坐标进行平移,直至两曲线重合;②任取一匹配点,记录下匹配点的对应坐标值:[t],[τ],[α],其中,考虑流速情况下,还需要记录下对应标准曲线的值:[β]。已知水的热容量(ρc)f与井管半径rw,根据不考虑流速情况下标准曲线依据求出含水层的滞止导热系数根据考虑流速情况下标准曲线依据可以求出含水层的有效导热系数λ,依据λv=λ-λ*可以求出含水层的热机械弥散系数λv,依据求出热弥散度ax。同时,依据求出含水层的热容量(ρc)fs。此外,根据室内土工试验资料可已知孔隙度θ,同时,在求出含水层的热容量(ρc)fs基础上,根据(ρc)fs=θ(ρc)f+(1-θ)(ρc)s可求出岩土体固体骨架热容量(ρc)s,依据求出岩土体导热系数λ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所示,然后根据试验数据在相同模的透明双对数纸上作实测曲线,使实测曲线坐标原点与标准曲线坐标原点保持同一高度,通过平移实测曲线时间坐标轴t来匹配实测曲线与标准曲线;任取一匹配点,记下匹配点的对应坐标[τ]值和[t]值,记录拟合标准曲线的[α]值;根据公式计算出含水层滞止导热系数与含水层有效导热系数,并可以反推出含水层热机械弥散系数与含水层热弥散度。其中不考虑流速情况下参数计算的配线图如图8至图10所示,滞止导热系数、含水层热容量与根据式(47)反推的岩土体导热系数计算结果见表1和表2。考虑流速情况下参数计算的配线图如图11至图13所示,含水层有效导热系数、含水层热机械弥散系数与根据式(48)反推的含水层热弥散度计算结果见表3和表4。
表1不考虑流速情况下的微热试验模型计算含水层导热系数的结果
表2岩土体导热系数计算结果
从表1和表2可以看出,本次室内模拟试验的含水层孔隙度为θ=0.33,(ρc)f=4.18×106J/(m3·℃),计算求得岩土体导热系数的平均值为4.465W/(m·℃)。本次试验采用的是粒径在0.25mm-0.50mm的中细砂。根据查找相关手册和文献,中细砂的导热系数一般在1.7~5.0W/(m·℃)范围内,本次试验求取的岩土体导热系数刚好在此范围内,由此也可证明微热试验确定含水层导热系数以及岩土体导热系数的准确性和可行性。
表3考虑流速情况下的微热试验模型计算含水层导热系数的结果
表4含水层热弥散度计算结果
从表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处的热量相等,建立热量平衡方程;
径向流下,考虑热对流弥散项的热量运移方程:
根据热量守恒定律,井管内热量的变化与流过含水层r=rw处的热量相等:
初始条件
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为初始温度变化值;为地下水实际流速,B为含水层厚度,Q为注水流量,θ为含水层孔隙度;
当u=0时,即含水层不考虑流速时,式(1)变成如下形式,即热的传导方程:
式(1)至式(6)数学模型的方程形式比较复杂,且因公式(1)和(7)表达式不同,需将公式(1)方程分成考虑流速和不考虑流速两种条件的情况,分别计算,总体思路是采用拉普拉斯变换的方法对其进行求解,将上述方程、定解条件等通过无量纲因子和参变量转换为无量纲形式,依次对其进行拉普拉斯变换;
(1)不考虑流速情况下的理论模型求解:
通过无量纲变换将式(7)转换为如下形式::
再依次对式(2)-(6)做无量纲变化得:
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)热传导方程转化为:
井管内热量的变化与流过含水层r=rw处的热量相等,该式转化为
定解条件:
上述拉普拉斯转换式中:
为T*D的拉氏变换的象函数,为TD的拉氏变换的象函数,T*D为井孔中温度变化比值,为含水层中温度变化比值,为距井中心距离与井半径的比,井内导热管半径与井半径的比,δw=1,p为拉氏变换复变量,无量纲时间, 为与含水层固有性质相关的无量纲参数,为水的热容量与含水层的热容量的比;
根据n阶修正贝塞尔方程:
通解为:
y=C1Kn(βx)+C2In(βx) (20)
其中:In(x)和Kn(x)分别为第一类和第二类修正贝塞尔函数;
当n=0时,
式(14)的解为:
根据式(17)知:
又因为第一类修正贝塞尔函数I0(x)为指数增长型函数,具有性质:当x→∞时I0(x)→∞,因此C2=0,所以,
又根据0阶修正贝塞尔函数的性质知:
将式(22)代入式(15)和(18)中可得:
由(18)可得:
由此可以得到当含水层处于静止条件下,流速u=0时该理论模型在Laplace空间的解析解为:
公式(26)即为在不考虑流速条件下的拉普拉斯变换域中所得出的最终解,但是想要在实际实验中运用此解,还需要对公式(26)进行拉普拉斯逆变换:
将公式(26)代入公式(27)中进行拉普拉斯逆变换得:
综上可得:
T*D=f(α,γ,τ) (29)
由此可知,在不考虑流速情况下微热试验的理论模型中,控制测试井中无量纲温度变化T*D的无因次参数共3个,其中为水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,为与含水层热容量相关的参数,δs反映了井管内导热管半径大小和井孔半径大小的比值,在假设条件中,将忽略导热管和井孔间的空隙,所以δs为1,τ为无因次时间,由此可利用Matlab软件绘制关于α参量的无量纲温度变化随时间变化标准曲线;
(2)考虑流速情况下的理论模型求解:
通过无量纲变换可将式(1)转换成如下形式:
对无量纲因子τ做Laplace变换可得:
xy″+βy′-pxy=0 (32)
y′=x-mu′-mx-m-1u (33)
y″=x-mu″-2mx-m-1u′+m(m+1)x-m-2u (34)
将式(33)和(34)代入(32)得:
根据n阶修正贝塞尔方程的解,可得:
由此可得,当n=0时,
同样根据定解条件以及第一类修正贝塞尔函数I0(x)具有如下性质:当x→∞时I0(x)→∞,可得到C2=0;
再将式(39)和式(18)代入无量纲化后的热量守恒条件即式(15),得:
由此可以得到当含水层在考虑流速情况下,微热试验理论模型在Laplace空间的解析解为:
和不考虑流速情况下求解过程相同,采用精度更高,考虑了复变域的Talbot方法进行Laplace逆变换:
将公式(43)代入公式(44)中进行拉普拉斯逆变换,综上可得:
T*D=f(α,β,γ,τ) (45)
由此可以得到在考虑流速情况下的微热试验理论中,控制测试井中无量纲温度变化的无因次参数共4个;与不考虑流速情况下一致,α指的是水的热容量和含水层的热容量之比,反映了含水层相对于水的储热能量强弱,β的大小与注入含水层内部的水流量以及含水层的性质相关,且反映了热量在含水层的运移过程中随着水流量的增加,对流弥散作用明显增强,远远大于分子扩散;
再次对β进行分析:
其中,为热扩散率,反映了含水层的导热速率的快慢,为水的热容量和含水层的热容量之比;根据β的表达式可以看出,β的大小与含水层的导热系数息息相关,在考虑流速情况下含水层的导热系数与流速、介质的孔隙度以及渗透系数有关,且0<β<1;γ与参数α和δs密切相关,除了反映含水层本身的热物性外,还反映了微热试验中导热管半径大小对试验结果的影响,τ为无因次时间;由此可固定其他参数值,利用Matlab软件绘制关于α、β参量的无量纲温度变化随时间变化的标准曲线;
热对流弥散模型中含水层的有效导热系数λ是由滞止导热系数λ*和热机械弥散系数λv两部分组成;由此,含水层中有效导热系数系数λ可表示为:
λ=λ*+λv (46)
含水层滞止导热系数λ*指的是地下水不流动时的含水层导热系数,属含水层(包含含水层骨架和地下水)自身的热物理参数,与地下水流速、流量等因素无关,即在不考虑流速情况下的所需求解的参数,此时热的传递方式只有传导作用,同样含水层的滞止导热系数可由下式计算:
式中:θ为含水层孔隙度,λ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)中,微热试验是通过瞬间增加或减少井孔内的温度,记录井孔内温度随时间变化,研究其温度变化规律,进而求取含水层的热物性参数。
8.如权利要求1或2所述的新型确定含水层热物性参数的微热试验方法,其特征在于:步骤3)中,求取含水层导热系数并通过公式推导其余热物性参数的步骤如下:①在Grapher软件中将不考虑流速和考虑流速情况下的实测曲线分别和相应条件下的标准曲线对应,置于同样模数的双对数坐标下,固定纵坐标在同一高度对横坐标进行平移,直至两曲线重合;②任取一匹配点,记录下匹配点的对应坐标值:[t],[τ],[α],其中,考虑流速情况下,还需要记录对应标准曲线的值:[β];已知水的热容量(ρc)f与井管半径rw,根据不考虑流速情况下标准曲线依据求出含水层的滞止导热系数λ*;根据考虑流速情况下标准曲线依据可以求出含水层的有效导热系数λ,依据λv=λ-λ*可以求出含水层的热机械弥散系数λv,依据求出热弥散度ax;同时,依据求出含水层的热容量(ρc)fs;此外,根据室内土工试验资料可已知孔隙度θ,同时,在求出含水层的热容量(ρc)fs基础上,根据(ρc)fs=θ(ρc)f+(1-θ)(ρc)s可求出岩土体固体骨架热容量(ρc)s,依据求出岩土体导热系数λs。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115408846A (zh) * | 2022-08-25 | 2022-11-29 | 中国科学院武汉岩土力学研究所 | 珊瑚细砂热物理参数简化预测模型及其试验装置和方法 |
Citations (4)
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 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种多尺度热量运移同步监测试验系统和试验方法 |
-
2022
- 2022-04-25 CN CN202210440391.3A patent/CN114813828B/zh active Active
Patent Citations (4)
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)
Title |
---|
ZHOU YAN-ZHANG 等: "Simulation of Thermal Transport in Aquifer: A GWHP System in Chengdu, China", 《JOURNAL OF HYDRODYNAMICS, SER. B》 * |
刘国庆等: "小尺度含水层热量运移试验研究及热弥散效应评估", 《岩土力学》 * |
吕培等: "分层微水试验与数值模拟的正反演验证对比分析", 《勘察科学技术》 * |
Cited By (2)
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 |