Processing math: 0%

不同油田废弃井型转地热开发规律与参数

石宇, 彭俊岚, 孙文超, 张宏源, 刘斌, 李冰

石宇, 彭俊岚, 孙文超, 张宏源, 刘斌, 李冰. 不同油田废弃井型转地热开发规律与参数[J]. 华南师范大学学报(自然科学版), 2024, 56(4): 27-38. DOI: 10.6054/j.jscnun.2024047
引用本文: 石宇, 彭俊岚, 孙文超, 张宏源, 刘斌, 李冰. 不同油田废弃井型转地热开发规律与参数[J]. 华南师范大学学报(自然科学版), 2024, 56(4): 27-38. DOI: 10.6054/j.jscnun.2024047
SHI Yu, PENG Junlan, SUN Wenchao, ZHANG Hongyuan, LIU Bin, LI Bing. The Development Laws and Parameters of Geothermal Extraction from Different Abandoned Well Types in Oil Fields[J]. Journal of South China Normal University (Natural Science Edition), 2024, 56(4): 27-38. DOI: 10.6054/j.jscnun.2024047
Citation: SHI Yu, PENG Junlan, SUN Wenchao, ZHANG Hongyuan, LIU Bin, LI Bing. The Development Laws and Parameters of Geothermal Extraction from Different Abandoned Well Types in Oil Fields[J]. Journal of South China Normal University (Natural Science Edition), 2024, 56(4): 27-38. DOI: 10.6054/j.jscnun.2024047

不同油田废弃井型转地热开发规律与参数

基金项目: 

国家自然科学基金面上项目 52374010

中国石油天然气集团有限公司科学研究与技术开发项目 2020A-4811(JT)

中深层地热井换热技术研究开发项目 2021DJ5502(JT)

地热资源开发利用关键技术研究 

西南交通大学新型交叉学科培育基金项目 2682023ZTPY030

西南交通大学新型交叉学科培育基金项目 2682022KJ034

详细信息
    通讯作者:

    石宇,Email: shiyu@swjtu.edu.cn

  • 中图分类号: TK529

The Development Laws and Parameters of Geothermal Extraction from Different Abandoned Well Types in Oil Fields

  • 摘要:

    为解决传统地下热水直接利用方式存在的潜在地质灾害和环境污染等问题,提出了一种“取热不取水”的、以同轴套管换热系统为关键系统的新型地热开发方式,建立了直井与多分支井同轴套管闭式换热系统数值模型,对比分析取热效果并开展了参数敏感性分析。结果表明:相同条件下,多分支井系统较直井系统具有更高的出口温度和取热功率,但导致注入压力显著上升。两种井型对于各参数的响应基本一致,仅程度不同,其中注入温度、流量、保温管长度、地层和水泥导热系数会显著影响系统的取热效果。该结果对深层地热开采的井型选择具有指导意义。

    Abstract:

    To address the potential geological hazards and environmental pollution associated with traditional direct utilization methods of geothermal water, a new geothermal development approach, emphasizing "taking heat without water" and utilizing a coaxial casing heat exchange system as the core component, is proposed. Numerical models of the coaxial casing closed-loop heat transfer system are established for both straight wells and multi-branch wells, comparing heat extraction efficiency and conducting a sensitivity analysis of the parameters involved. The results demonstrate that the multi-branch well system achieves higher outlet temperatures and heat extraction power than the straight well system under identical conditions; however, this also leads to a significant increase in injection pressure. The different well types show similar responses to varying parameters, with differences primarily in degree. Factors that significantly influence the system's heat extraction efficiency include the injection temperature, flow rate, length of the insulation pipe, and thermal conductivity of the formation and cement. The study provides valuable guidance on selecting well types for deep geothermal extraction.

  • 在拓扑动力系统中,拓扑熵[1]是一个拓扑共轭不变量,可用来刻画系统的复杂程度。非紧集的拓扑熵[2-3]与测度熵、Lyapunov指数等密切相关[4-11],可用来研究系统中的重分形谱[12]、饱和集[6]和非正则集等[5, 7]

    学者们从多个角度对拓扑熵进行了深入研究,发现系统中满足一定条件的零测集具有与整个系统相同的拓扑熵或Hausdorff维数。如:BARREIRA和SCHMELING[4]证明了针对拓扑混合的有限型子转移,如果Hölder连续函数g对0是非上同调的,那么g的Birkhoff平均不规则集上的拓扑熵等于这个有限型子转移的拓扑熵; PESIN[5]证明了针对符号系统,非generic点集上的拓扑熵等于整个系统的拓扑熵; CHEN等[6]证明了满足specification性质的非正则集上的拓扑熵等于整个系统的拓扑熵; TIAN[7]证明了系统中若干周期类点集的差集上的拓扑熵等于整个系统的拓扑熵; ZHU和MA[8]将文献[6-7]的结论推广到了自由半群作用的系统中。

    在迭代函数系统(IFS)中,ELTON[13]得到了一个遍历定理(Elton定理):对任意连续函数,符号空间中几乎所有的点,IFS的时间平均等于空间平均。受文献[4-8]的启发,本文探讨符号空间中Elton定理不成立的点构成的集合的性质及其拓扑熵、Hausdorff维数与整个系统的拓扑熵、Hausdorff维数的关系。

    设(X, d)是紧致度量空间。映射T: XX称为一个压缩映射,如果d(T(x),T(y)), \forall x, y \in X,其中θ(0 < θ < 1)是一个常数。如果{T0, T1, …, Tm-1}是X上的压缩映射族,那么称(X, T0, T1, …, Tm-1)为迭代函数系统(IFS)。由文献[14]可知,对于迭代函数系统(X, T0, T1, …, Tm-1),如果存在唯一的非空子集E \subset X,使得E=\bigcup\limits_{j=0}^{m-1} T_j(E),则E被称为\left(X, T_0, T_1, \cdots, T_{m-1}\right)的吸引子。

    设(X, T0, T1, …, Tm-1)是一个迭代函数系统,(p0, p1, …, pm-1)是一个概率向量,即p_i>0(0 \leqslant i \leqslant m-1), \sum_{i=0}^{m-1} p_i=1。HUTCHINSON[14]证明了存在唯一的Borel概率测度ν,使得对所有的Borel集Z \subset X,有

    \nu(Z)=\sum\limits_{j=0}^{m-1} p_j \nu\left(T_j^{-1} Z\right)

    且supp ν是(X, T0, T1, …, Tm-1)的吸引子E, 这个概率测度ν称为是(X, T0, T1, …, Tm-1)的关于概率向量(p0, p1, …, pm-1)的不变测度。

    \varSigma=\prod\limits_{n=1}^{\infty}\{0, 1, \cdots, m-1\}为符号空间,(p0, p1, …, pm-1)是一个概率向量,P为由每一个因子的测度p({i})=pi诱导的乘积测度。

    1987年,ELTON[13]得到了迭代函数系统(X, T0, T1, …, Tm-1)中的一个遍历定理:

    定理A[13]   (Elton定理)对于X上的任意连续函数f和任意x \in X, 有

    \lim _{n \rightarrow \infty} \frac{1}{n} \sum\limits_{i=0}^{n-1} f\left(T_{j_i} \circ T_{j_{i-1}} \circ \cdots \circ T_{j_0}(x)\right)=\int_X f \mathrm{~d} \nu,

    相对于乘积测度P而言,对几乎所有的(j0, j1, …)\in Σ成立。

    对任意f\in C(X)和任意x \in X,记f的Elton定理不成立的点构成的集合为

    \begin{aligned} &B_1(f, x)=\left\{\left(j_0, j_1, \cdots\right) \in \varSigma \mid \lim _{n \rightarrow \infty} \frac{1}{n} \sum\limits_{i=0}^{n-1} f\left(T_{j_i} \circ T_{j_{i-1}} \circ \cdots \circ\right.\right. \\ &\left.T_{j_0}(x)\right) \text { 不存在\}。 } \end{aligned}

    根据定理A,对Σ上的每一个乘积测度P,均有P(B1(f, x))=0,B1(f, x)是σ-不变的,其中σΣ上的推移映射。

    本文研究集合B1(f, x)的拓扑熵、Hausdorff维数与整个系统的拓扑熵、Hausdorff维数的关系。

    Am×m阶矩阵,A的每个元素aij取值为0或1。设\varSigma_A \subset \varSigma是紧致的σ-不变子集,使得任意(i0, i1, …)\in ΣA,满足ain in+1 =1,n≥0,映射σ|ΣA称为由矩阵A决定的有限型子转移。称σ|ΣA拓扑混合当且仅当存在一个正整数k,使得Ak的所有元素的值都是正的。σ|ΣA的拓扑熵为h(σ|ΣA)=log ρ(A)[15],其中ρ(A)表示A的谱半径。

    如果存在连续函数\psi: \varSigma_A \rightarrow \mathbb{R},存在c \in \mathbb{R},使得g_1-g_2=\psi-\psi \circ \sigma+c,则称2个函数g_1, g_2: \varSigma_A \rightarrow \mathbb{R} 是上同调的。对于每一个连续函数g: \varSigma_A \rightarrow \mathbb{R},定义g的Birkhoff平均的不规则集为

    B(g)=\left\{\left(j_0, j_1, \cdots\right) \in \varSigma_A \mid \lim _{n \rightarrow \infty} \frac{1}{n} \sum\limits_{i=1}^n g\left(\sigma^k\left(j_0, j_1, \cdots\right)\right)\right. \text { 不 } \\ 存在\}。

    注意,集合B(g)是σ-不变的。根据Birkhoff遍历定理[15],对于ΣA上每一个σ-不变测度μ,有μ(B(g))=0。1997年,BARREIRA和SCHMELING[4]证明了满足一定条件的Birkhoff平均不规则集上的拓扑熵等于这个有限型子转移的拓扑熵。

    定理B[4]   对于拓扑混合的有限型子转移,如果Hölder连续函数g对0是非上同调的,那么

    h(\sigma \mid B(g))=h\left(\sigma \mid \varSigma_A\right),

    其中,h(σ|B(g))、h(σ|ΣA)分别表示B(g)、ΣA上的σ拓扑熵,B(g)表示g的Birkhoff平均不规则集,ΣA表示由矩阵A确定的符号空间Σσ-不变子集。

    下面给出本文定理证明需要用到的一个结论:

    定理C[12]   对符号系统(Σ, σ),对任意Z \subset \Sigma,有

    h(\sigma \mid Z)=\mathrm{HD}(Z) \cdot \log m,

    其中,Σm个符号生成,h(σ|Z)表示Z上的拓扑熵,HD(Z)表示Z上的Hausdorff维数。

    本文主要结论如下:

    定理1   设(X, T0, T1, …, Tm-1)是一个迭代函数系统,对任意g\in C1(X)和任意x \in X,这里C_1(X)=\{g \in \left.C(X) \mid \int_X g \mathrm{~d} \nu_1 \neq \int_X g \mathrm{~d} \nu_2, \forall \nu_1, \nu_2 \in M(X), \nu_1 \neq \nu_2\right\}, M(X)表示(X, T0, T1, …, Tm-1)中关于所有概率向量的不变测度的集合,有

    h\left(\sigma \mid B_1(g, x)\right)=h(\sigma)=\log m,

    B1(g, x)是满熵的,且具有满Hausdorff维数。

    证明   仿照定理B的证明方法,设(p0, p1, …, pm-1)和(q0, q1, …, qm-1)是任何概率向量,且

    \left(p_0, p_1, \cdots, p_{m-1}\right) \neq\left(q_0, q_1, \cdots, q_{m-1}\right)。

    P1P2分别表示由(p0, p1, …, pm-1)、(q0, q1, …, qm-1)诱导的乘积测度,ν1ν2分别表示(X, T0, T1, …, Tm-1)关于(p0, p1, …, pm-1)、(q0, q1, …, qm-1)的不变测度,则有P1P2, ν1ν2。给定ε>0,存在概率向量(p0, p1, …, pm-1)和(q0, q1, …, qm-1),使得(p0, p1, …, pm-1)≠(q0, q1, …, qm-1),hPi(σ)>h(σ)-ε(i=1, 2)。其中hPi(σ)表示关于Piσ测度论熵。因为ν1ν2,所以\int_X g \mathrm{~d} \nu_1 \neq \int_X g \mathrm{~d} \nu_2, g \in C_1(X)。选择δ\in (0, ε), 使得\left|\int_X g \mathrm{~d} \nu_1-\int_X g \mathrm{~d} \nu_2\right|>4 \delta。固定x \in X, 设Γil(i=1, 2;l≥1)表示Σ中满足下面不等式的点(j0, j1, …)构成的集合:对任意nl,满足

    \left|\frac{1}{n} \sum\limits_{i=0}^{n-1} g\left(T_{j_i} \circ T_{j_{i-1}} \circ \cdots \circ T_{j_0}(x)\right)-\int_X g \mathrm{~d} \nu_i\right| <\delta。

    p_s=s(\bmod 2), l_s是正整数的递增序列,使得对每一个整数s \geqslant 1, 有P_{p_s}\left(\Gamma_{P_s}^{l_s}\right)>1-\varepsilon / 2^s。令m_1= n_1=l_1, m_s=\left(n_{s-1}+1+l_{s+1}\right)!, n_s=n_{s-1}+1+m_s。构造如下柱集族:\mathscr{C}_s=\left\{C_{m_s}\left(j_0, j_1, \cdots\right) \mid\left(j_0, j_1, \cdots\right) \in \Gamma_{p_s}^{l_s}\right\}, \mathscr{D}_1= \mathscr{C}_1\mathscr{D}_s=\left\{\underline{C} C \bar{C} \mid \underline{C} \in \mathscr{D}_{s-1}, \bar{C} \in \mathscr{C}_s\right., 柱集的长度|C|=1}。

    \varLambda=\bigcap\limits_{s \geqslant 1} \cup C ,定义Λ上的测度μ如下:

    \begin{gathered} \mu(C)=P_1(C), C \in \mathscr{D}_1=\mathscr{C}_1, \\ \mu(\underline{C} C \bar{C})=\mu(\underline{C}) P_{p_s}(\bar{C}), \underline{C} C \bar{C} \in \mathscr{D}_s, s>1 。 \end{gathered}

    对于每个可测集A \subset \varSigma,令\mu(A)=\mu(A \cap \varLambda),则μ扩展到了Σ上。

    如果s>1, \underline{C} \in \mathscr{D}_{s-1},那么

    \mu\left(\bigcup\limits_{\bar{C} \in \mathscr{D}_s} \underline{C} \cap \bar{C}\right) \geqslant \mu(\underline{C})\left(1-\frac{\varepsilon}{2^s}\right)。

    如果ε < 2,那么

    \mu(\varLambda) \geqslant \prod\limits_{s=1}^{\infty}\left(1-\frac{\varepsilon}{2^s}\right)>0。

    \left(j_0, j_1, \cdots\right) \in C \in \mathscr{D}_s ,由于\left(j_{|C|-m_s}, j_{|C|-m_s+1}, \cdots\right) \in \Gamma_{p_s}^{l_s},当s \rightarrow \infty 时,有\frac{|C|}{m_s} \rightarrow 1足够大,那么

    \begin{aligned} &\left|\frac{1}{|C|} \sum\limits_{i=0}^{|C|} g\left(T_{j_i} \circ \cdots \circ T_{j_0}(x)\right)-\int_X g \mathrm{~d} \nu_{p_s}\right| \leqslant\\ &\left|\frac{1}{m_s} \sum\limits_{i=0}^{m_s} g\left(T_{j_{\mid C 1-m_s+i}} \circ \cdots \circ T_{j_1 \mid 1-m_s}(x)\right)-\int_X g \mathrm{~d} \nu_{p_s}\right| \times\\ &m_s \sum\limits_{i=0}^{|C|} g\left(T_{j_i} \circ \cdots \circ T_{j_0}(x)\right)\\ &|C| \sum\limits_{i=0}^{m_s} g\left(T_{j_{\mid C l}-m_s+i} \circ \cdots \circ T_{j_{|C|-m_s}}(x)\right) \end{aligned}
    \begin{aligned} &\left|1-\frac{m_s \sum\limits_{i=0}^{m_s} g\left(T_{j_i} \circ \cdots \circ T_{j_0}(x)\right)}{|C| \sum\limits_{i=0}^{m_s} g\left(T_{j_{|C|}-m_s+i} \circ \cdots \circ T_{\mathrm{j}_{|C|-m_s}}(x)\right)}\right| \times\\ &\left|\int_X g \mathrm{~d} \mu_{v_{p_s}}\right| <2 \delta。 \end{aligned}

    由于\delta \in(0, \varepsilon),可知\varLambda \subset B_1(g, x)。现在设\left(j_0\right., \left.j_1, \cdots\right) \in \varLambda,如果q足够大,那么

    -\frac{\log \left(\mu\left(C_q\left(j_0, j_1, \cdots\right)\right)\right)}{q} \geqslant h(\sigma)-\eta

    对某\eta \in(0, 2 \varepsilon)成立。这意味着

    h\left(\sigma \mid B_1(g, x)\right) \geqslant h(\sigma \mid \varLambda) \geqslant h_{\mu \mid \varLambda}(\sigma) \geqslant h(\sigma)-2 \varepsilon

    由于ε是任意的,则h\left(\sigma \mid B_1(g, x)\right)=h(\sigma)

    此外,由定理C可得

    \mathrm{HD}\left(B_1(g, x)\right)=\frac{h\left(\sigma \mid B_1(g, x)\right)}{\log m}=\frac{h(\sigma)}{\log m}=\mathrm{HD}(\varSigma)。

    证毕。

    例1   设X=[0, 1], T_0: x \mapsto \frac{1}{3} x, T_1: x \mapsto \frac{1}{3} x+ \frac{2}{3} \circ\left(X, T_0, T_1\right)的吸引子是Cantor集E

    \begin{aligned} E=& T_0(E) \cup T_1(E)=\left\{x_{i_0 i_1 \cdots}:=\bigcap\limits_{n=0}^{\infty} T_{i_0} \circ T_{i_1} \circ \cdots \circ T_{i_n}(E) \mid\right.\\ &\left.\left(i_0, i_1, \cdots\right) \in \prod\limits_{n=1}^{\infty}\{0, 1\}\right\}。 \end{aligned}

    g\in C(X), 使得g(x)=a, x\in [0, 1/3], g(x)=b, x\in [2/3, 1], 这里ab。对于任意x \in X,令E(g, x)={xi0 i1… |(i0, i1, …)\in B1(g, x)} E。那么以下结论成立:

    (1) 对于任意的x \in X,有

    h\left(\sigma \mid B_1(g, x)\right)=h(\sigma)=\log 2, \mathrm{HD}\left(B_1(g, x)\right)=\mathrm{HD}(\varSigma),

    这里\varSigma=\prod\limits_{n=1}^{\infty}\{0, 1\}

    (2) 对于任意x \in X,有

    \mathrm{HD}(E(g, x))=\mathrm{HD}(E)=\log 2 / \log 3。

    例1中的结论(1)可由定理1直接得出。结论(2)中集合E(g, x)是Cantor集E的一个子集,由于系统\left(\prod\limits_{n=1}^{\infty}\{0, 1\}, \sigma\right)与(E, ψ)是拓扑共轭的,E中任一集合与其在\prod\limits_{n=1}^{\infty}\{0, 1\}中的象集有相同的拓扑熵。针对系统(E, ψ),集合E(g, x)上的Hausdorff维数与其拓扑熵有关系,因此,下面给出结论(2)证明中需用到的一个定理:

    定理D[11]   对紧致系统(X, f),若存在α>1,对任意充分小的正数ε,当d(x, y) < \varepsilon时,有d(f(x), f(y))=\alpha d(x, y),则对任意Z \subset X,有

    h(f \mid Z)=\mathrm{HD}(Z) \cdot \log \alpha。

    例1的证明   (1)设(p0, p1)和(q0, q1)为任意概率向量,且(p0, p1)≠(q0, q1)。那么p0q0, p1q1。设ν1ν2分别表示(X, T0, T1)关于(p0, p1)、(q0, q1)的不变测度, 则ν1ν2,且

    \int_X g \mathrm{~d} \nu_1=p_0 a+p_1 b, \int_X g \mathrm{~d} \nu_2=q_0 a+q_1 b。

    因此

    \begin{aligned} &\left(p_0 a+p_1 b\right)-\left(q_0 a+q_1 b\right)=a\left(p_0-q_0\right)+b\left(p_1-q_1\right)= \\ &\quad\left(p_0-q_0\right)(a-b) \neq 0 \end{aligned}

    从而有\int_X g \mathrm{~d} \nu_1 \neq \int_X g \mathrm{~d} \nu_2。根据定理1,对任意x \in X,有

    h\left(\sigma \mid B_1(g, x)\right)=h(\sigma)=\log 2, \operatorname{HD}\left(B_1(g, x)\right)=\mathrm{HD}(\varSigma) 。

    (2) 对于迭代函数系统(X, T0, T1),可以定义一个连续的自映射

    \psi: E \rightarrow E, \\ \psi(x)=T_j^{-1}(x), x \in T_j(E)(j=0, 1),

    那么(E, ψ)是一个离散动力系统。定义 \varphi: \prod\limits_{n=0}^{\infty}\{0, 1\} \rightarrow E, \varphi\left(i_0, i_1, \cdots\right)=\bigcap\limits_{n=0}^{\infty} T_{i_0} \circ T_{i_1} \circ \cdots \circ T_{i_n}(E): =x_{i_0 i_1 \cdots}, \forall\left(i_0, \right. \left.i_1, \cdots\right) \in \prod\limits_{n=0}^{\infty}\{0, 1\},则φE \prod\limits_{n=0}^{\infty}\{0, 1\}的同胚,有\varphi \circ \sigma=\psi \circ \varphi, 也就是说,φ是拓扑共轭,σ \prod\limits_{n=0}^{\infty}\{0, 1\}上的推移映射。易见E(g, x)=φ(B1(g, x))。由定理C、定理D和(1),可得

    \mathrm{HD}(E(g, x))=\frac{h(\psi \mid E(g, x))}{\log 3}=\frac{h\left(\sigma \mid B_1(g, x)\right)}{\log 3}=\frac{\log 2}{\log 3}。
  • 图  1   单井同轴套管闭式取热示意图

    Figure  1.   Schematic of closed-loop heat extraction for a single well with coaxial casing

    图  2   模型耦合过程

    Figure  2.   Model coupling process

    图  3   几何模型及网格划分方案

    Figure  3.   Geometric model and meshing scheme

    图  4   120天的出口温度随单元数量的变化曲线

    Figure  4.   Outlet temperature variation over 120 days with unit count

    图  5   120天井底与垂直剖面的温度云图

    Figure  5.   120-day bottom and vertical cross-section temperature distribution

    图  6   直井与多分支井系统取热效果对比

    Figure  6.   Heat extraction comparison between vertical and multi-branch wells

    图  7   直井与多分支井系统取热效果随注入温度的变化曲线

    Figure  7.   Heat extraction effect variation with injection temperature

    图  8   直井与多分支井系统取热效果随注入流量的变化曲线

    Figure  8.   Heat extraction effect variation with injection flow rate

    图  9   直井与多分支井系统取热效果随保温管长度的变化曲线

    Figure  9.   Heat extraction effect variation with insulation pipe length

    图  10   直井与多分支井系统取热效果随套管尺寸-内管内径的变化曲线

    Figure  10.   Heat extraction effect variation with casing size-inner diameter of inner pipe

    图  11   直井与多分支井系统取热效果随水泥导热系数的变化曲线

    Figure  11.   Heat extraction effect variation with cement thermal conductivity

    图  12   直井与多分支井系统取热效果随地层导热系数的变化曲线

    Figure  12.   Heat extraction effect variation with formation thermal conductivity

    表  1   基础模型参数

    Table  1   Fundamental model parameters

    参数 参数值
    地层导热系数/(W·m-1·K-1) 2
    水泥导热系数/(W·m-1·K-1) 0.7
    套管导热系数/(W·m-1·K-1) 43.75
    内管导热系数/(W·m-1·K-1) 0.025
    内管内径/m 0.062
    外管内径/m 0.114 3
    外管外径/m 0.201 2
    套管外径/m 0.219 1
    井筒直径/m 0.269 9
    下载: 导出CSV

    表  2   地层物性参数

    Table  2   Formation physical properties

    地层参数 地层1 地层2 地层3 地层4
    孔隙度 0.55 0.257 0.186 0.105
    密度/(kg·m-3) 2 700 2 700 2 650 2 650
    热容/ (J/(kg-1·K-3) 950 950 1 000 1 000
    下载: 导出CSV

    表  3   分支井眼尺寸

    Table  3   Lateral wellbore dimensions

    参数 参数值
    内管内径/m 0.050 6
    环空内径/m 0.088 9
    环空外径/m 0.159 4
    套管外径/m 0.178 8
    井筒直径/m 0.215 9
    下载: 导出CSV

    表  4   敏感性分析参数设置表

    Table  4   Sensitivity analysis parameters

    参数 参数值
    注入温度/K 288,293,298,303,308
    注入流量/(m3·h-1) 15,20,25,30,35
    保温管长度/m 0,400,850,1300,1 750,2 200,2 650,3 100,3 500
    地层导热系数/(W·m-1·K-1) 1.5,2.0,2.5,3.0,3.5
    水泥导热系数/(W·m-1·K-1) 0.4,0.7,1.0,1.3,1.6
    下载: 导出CSV

    表  5   套管尺寸

    Table  5   Casing dimensions m

    参数 参数值
    内管内径 0.057 4,0.062 0,0.076 0,0.100 5,0.118 6
    环空内径 0.114 3,0.114 3,0.139 7,0.139 7,0.177 8
    环空外径 0.201 2,0.201 2,0.216 8,0.216 8,0.258 9
    套管外径 0.219 1,0.219 1,0.244 5,0.244 5,0.273 0
    井筒直径 0.269 9,0.269 9,0.269 9,0.269 9,0.298 5
    下载: 导出CSV
  • [1] 贾艳雨, 常青, 王俞文, 等. 我国地热能开发利用现状及双碳背景下的发展趋势[J]. 石油石化绿色低碳, 2021, 6: 5-9.

    JIA Y Y, CHANG Q, WANG Y W, et al. Geothermal energy development amid carbon peak and neutrality in China[J]. Green Petroleum & Petrochemicals, 2021, 6: 5-9.

    [2] 张聚斌, 邢晓森, 李振雄. 关于开发利用中深层地热取热不取水技术的思考[J]. 工程地球物理学报, 2021, 18: 400-408.

    ZHANG J B, XING X S, LI Z X. Reflections on the development and utilization of middle and deep geothermal energy without water extraction[J]. Chinese Journal of Engineering Geophysics, 2021, 18: 400-408.

    [3] 任虎俊. 水热型地热能同轴管换热技术研究——以河北省邯郸地区为例[J]. 中国煤炭地质, 2018, 30: 105-108.

    REN H J. Study on hydrothermal type geothermal energy coaxial tube heat transfer technology-a case study of Handan area[J]. Coal Geology of China, 2018, 30: 105-108.

    [4]

    BEIER R A. Transient heat transfer in a U-tube borehole heat exchanger[J]. Applied Thermal Engineering, 2013: 256-266.

    [5]

    LV Z, SONG X, LI G, et al. Numerical analysis of characteristics of a single U-tube downhole heat exchanger in the borehole for geothermal wells[J]. Energy, 2017: 186-196.

    [6] 宋先知, 张逸群, 李根生, 等. 雄安新区地热井同轴套管闭式循环取热技术研究[J]. 天津大学学报(自然科学与工程技术版), 2021, 54: 971-981.

    SONG X Z, ZHANG Y Q, LI G S, et al. Performance study of the downhole coaxial closed-loop heat exchange technology in Xiong'an New Area[J]. Journal of Tianjin University(Science and Technology), 2021, 54: 971-981.

    [7]

    HORNE R N. Design considerations of a down-hole coaxial geothermal heat exchanger[Z]. Geothermal Resource Council Annual Meeting. Salt Lake City: [S. n. ]. 1980: 569-572.

    [8]

    BEIER R A, ACUñA J, MOGENSEN P, et al. Transient heat transfer in a coaxial borehole heat exchanger[J]. Applied Thermal Engineering: Design, Processes, Equipment, Economics, 2014, 62: 256-266.

    [9]

    ZANCHINI E, LAZZARI S, PRIARONE A. Effects of flow direction and thermal short-circuiting on the performance of small coaxial ground heat exchangers[J]. Renewable Energy, 2010, 35(6): 1255/1-11.

    [10]

    ZHAO J, WANG H, LI X, et al. Experimental investigation and theoretical model of heat transfer of saturated soil around coaxial ground coupled heat exchanger[J]. Applied Thermal Engineering, 2008, 28(2): 116-125.

    [11]

    HOLMBERG H, ACUÑA J, NÆSS E, et al. Thermal evaluation of coaxial deep borehole heat exchangers[J]. Renewable Energy, 2016, 97: 65-76. doi: 10.1016/j.renene.2016.05.048

    [12]

    CHENG W L, LI T T, NIAN Y L, et al. An analysis of insulation of abandoned oil wells reused for geothermal power generation[J]. Energy Procedia, 2014, 61: 607-610. doi: 10.1016/j.egypro.2014.11.1181

    [13]

    GUO T, ZHANG Y, HE J, et al. Research on geothermal development model of abandoned high temperature oil reservoir in North China oilfield[J]. Renewable Energy, 2021, 177: 1-12. doi: 10.1016/j.renene.2021.05.128

    [14]

    ALIMONTI C, SOLDO E. Study of geothermal power generation from a very deep oil well with a wellbore heat exchanger[J]. Renewable Energy, 2016, 86: 292-301. doi: 10.1016/j.renene.2015.08.031

    [15] 张承良, 杨钊, 马庆利. 连续油管侧钻井的最新进展及新技术应用[J]. 国外油田工程, 2009, 25: 34-36.

    ZHANG C L, YANG Z, MA Q L. Latest developments and new technology applications of continuous lateral drilling[J]. Energy Conservation in Petroleum Petrochemical Industry, 2009, 25: 34-36.

    [16] 迟焕鹏, 李根生, 黄中伟, 等. 水力喷射径向水平井技术研究现状及分析[J]. 钻采工艺, 2013, 36: 119-124.

    CHI H P, LI G S, HUANG Z W, et al. Research status and analysis of water jet radial horizontal well technology[J]. Drilling & Production Technology, 2013, 36: 119-124.

    [17]

    WANG G, SONG X, SHI Y, et al. Heat extraction analysis of a novel multilateral-well coaxial closed-loop geothermal system[J]. Renewable Energy, 2021, 163: 974-986.

    [18] 刘豪, 王英民. 用地质统计学方法确定储集层物性参数的空间分布——以辽河油田沈35区块井间孔隙度插值为例[J]. 新疆石油地质, 1999(1): 36-39.
    [19] 王伟明, 卢双舫, 田伟超, 等. 吸附水膜厚度确定致密油储层物性下限新方法——以辽河油田大民屯凹陷为例[J]. 石油与天然气地质, 2016, 37: 135-140.

    WANG W M, LU S F, TIAN W C, et al. A new method to determine porosity and permeability cutoffs of tight oil reservoirs by using thickness of adsorption water film: a case study from the Damintun Sag, Liaohe oifield[J]. Oil & Gas Geology, 2016, 37: 135-140.

    [20] 汪少勇, 李建忠, 王社教, 等. 辽河坳陷雷家地区沙四段致密储层孔隙结构及物性下限[J]. 东北石油大学学报, 2016, 40(1): 51-61.

    WANG S Y, LI J Z, WANG S J, et al. Pore-structure and physical property cut-off of tight reservoirs in the fourth Shahejie formation, Leijia district, Liaohe depression[J]. Journal of Northeast Petroleum University, 2016, 40(1): 51-61.

    [21] 蒋宏伟, 王克雄, 翟应虎, 等. 辽河油田小龙湾地区粗面岩油藏岩石力学特性研究[J]. 石油钻探技术, 2007(1): 26-28.

    JIANG H W, WANG K X, ZHAI Y H, et al. Study of mechanics characteristics of xiaolongwan trachyte reservoir in Liaohe oilfield[J]. Petroleum & Drilling Techniques, 2007(1): 26-28.

    [22]

    WEI W, ZHU X, MENG Y, et al. Porosity model and its application in tight gas sandstone reservoir in the southern part of West Depression, Liaohe Basin, China[J]. Journal of Petroleum Science and Engineering, 2016, 141: 24-37.

    [23]

    HAALAND S E. Simple and explicit formulas for the friction factor in turbulent pipe flow[J]. Journal of Fluids Engineering, 1983, 105(1): 89-90.

    [24]

    INCROPERA F P, WITT D P D. Fundamentals of heat and mass transfer[M]. 6th edition. New York: Wiley, 2007.

    [25]

    WANG G, SONG X, SHI Y, et al. Numerical investigation on heat extraction performance of an open loop geothermal system in a single well[J]. Geothermics, 2019, 80: 170-184.

    [26] 杨世铭, 陶文铨. 传热学[M]. 第4版. 北京: 高等教育出版社, 2006.
    [27]

    GNIELINSKI V. New equations for heat and mass transfer in the turbulent flow in pipes and channels[M]. Washington DC: NASA STI, 1975.

    [28]

    CHEN J, JIANG F. A numerical study of EGS heat extraction process based on a thermal non-equilibrium model for heat transfer in subsurface porous heat reservoir[J]. Heat and Mass Transfer, 2016, 52(2): 255-267.

图(12)  /  表(5)
计量
  • 文章访问数:  36
  • HTML全文浏览量:  9
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-01-30
  • 刊出日期:  2024-08-24

目录

/

返回文章
返回