Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
  • Scopus数据库收录期刊
  • 中国科技核心期刊
  • 全国中文核心期刊
  • 美国工程索引(EI)收录期刊

基于多元概率分布模型的珠海黏土多参数预测

汪海林, 刘航宇, 顾晓强, 宋许根

汪海林, 刘航宇, 顾晓强, 宋许根. 基于多元概率分布模型的珠海黏土多参数预测[J]. 岩土工程学报, 2021, 43(S2): 193-196. DOI: 10.11779/CJGE2021S2046
引用本文: 汪海林, 刘航宇, 顾晓强, 宋许根. 基于多元概率分布模型的珠海黏土多参数预测[J]. 岩土工程学报, 2021, 43(S2): 193-196. DOI: 10.11779/CJGE2021S2046
WANG Hai-lin, LIU Hang-yu, GU Xiao-qiang, SONG Xu-gen. Multi-parameter prediction of Zhuhai clay based on multivariate probability distribution model[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(S2): 193-196. DOI: 10.11779/CJGE2021S2046
Citation: WANG Hai-lin, LIU Hang-yu, GU Xiao-qiang, SONG Xu-gen. Multi-parameter prediction of Zhuhai clay based on multivariate probability distribution model[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(S2): 193-196. DOI: 10.11779/CJGE2021S2046

基于多元概率分布模型的珠海黏土多参数预测  English Version

基金项目: 

国家自然科学基金项目 51822809

详细信息
    作者简介:

    汪海林(1998— ),男,博士研究生,主要从事岩土基本特性大数据统计分析方面的研究工作。E-mail: wanghailin@tongji.edu.cn

    通讯作者:

    顾晓强, E-mail: guxiaoqiang@tongji.edu.cn

  • 中图分类号: TU43

Multi-parameter prediction of Zhuhai clay based on multivariate probability distribution model

  • 摘要: 在岩土工程实践中,一个重要任务是根据场地勘察报告获取实际设计中所需要的设计参数。土体参数之间的相互转换对于岩土工程实践具有重要意义。本文在前人研究的基础上,对珠海市某城区的软土的室内试验及现场数据进行了分析,并建立了9个土体参数的多元概率分布模型以及相应的贝叶斯更新模型。使用该贝叶斯模型可以方便地使用多元土体参数信息对未知参数进行更新预测,研究表明,多参数更新预测模型能够综合所有已知土体参数的信息,从而有效提升对于未知参数预测的准确性。
    Abstract: In geotechnical engineering practice, an important task is to obtain the design parameters required based on the site survey report. The mutual conversion of soil parameters is of great significance to the practice of geotechnical engineering. Based on the previous studies, the laboratory test and field test data of soft soil in a certain urban area of Zhuhai are analyzed. A multivariate probability distribution model of 9 soil parameters and a Bayesian updating model are established. Using this Bayesian updating model can easily make use of the information of multiple soil parameters to predict unknown parameters. The results show that the multi-parameter prediction model can synthesize the information of all known parameters, thereby effectively improving the accuracy of predicting unknown parameters.
  • 煤层气作为一种高效且清洁的能源受到越来越多学者的关注[1]。在煤层气开采过程中会引起煤岩的力学变形和吸附解吸变形[2]。煤岩渗透率变化是这两种作用的竞争结果,裂隙与基质孔隙之间的气体压差将加剧这两种过程的复杂性,从而增加渗透率预测的难度。因此研究基质与裂隙内气体压差对正确认识储层内渗透率的动态演化具有重要意义。

    常用的煤岩渗透率理论模型可分为单孔隙模型与双重孔隙模型两大类。由于煤裂隙系统的渗透率远大于基质系统的渗透率[2],早期许多研究中将煤岩视为单孔隙模型[3-5],不存在气体压差的概念。双重孔隙介质模型对煤岩孔隙结构有很好的描述,区分了裂隙压力和基质孔隙压力[6-9],然而对孔隙压力的理解仅停留在不同介质系统具有不同的孔隙压力上,忽略了气体压差对煤岩变形的影响。这些传统模型未能解释应力控制边界条件下的实验室测量结果,在匹配现场数据方面也仅有有限的成功[10]

    上述研究中使用的煤岩理想化模型有火柴杆模型与立方体模型,模型中假设裂隙将基质完全隔开。而在实际情况中裂隙往往不能贯穿各基质块,基质之间通过基质桥连接,因此需要考虑基质与裂隙的相互作用[11]。这种相互作用使得基质吸附变形只有一部分影响裂隙开度变化,另一部分则影响煤整体变形。为了解决这一问题,一些研究通过引入一个值为0~1之间的内部膨胀系数f来表征基质吸附变形对裂隙开度变化的贡献[11-14],并通过数据拟合的方式将系数f确定为一个常量。然而,由于煤岩基质孔隙的致密性,抽注时裂隙压力与基质孔隙压力无法瞬间达到平衡。吸附性气体在裂隙与基质之间的长期扩散过程中,煤岩经历初始平衡、非平衡和最终平衡3个状态[15-16],在低压条件下甚至永远无法达到最终平衡[17]。非平衡态演化过程中基质吸附变形对裂隙开度变化的贡献是随时间变化的,演化初期裂隙附近的基质非均匀变形(局部变形)几乎完全影响裂隙开度变化(f≈1),而在最终平衡阶段,基质均匀变形(整体变形)几乎不影响裂隙开度变化(f≈0)。因此f被看作一个常量是不合理的。

    过去渗透率测量的实验室观察均在短期内完成,得到的试验数据无法反映煤岩的长期非平衡演化过程。Wei等[18]对煤岩试样进行了超过80 d的CO2注气试验,使用改进的脉冲衰减技术进行了渗透率连续测量,研究了煤岩从初始平衡到最终平衡的裂隙渗透率演化。试验结果显示,煤岩渗透率随时间先迅速增大后降低,最后再缓慢恢复。为了解释这一结果,近期许多研究基于非平衡理论进行了建模。Wei等[19-20]通过区分裂隙应变的方式描述了局部变形与整体变形之间的关系,又通过划分基质区域的方式区分了不同位置处的基质变形对裂隙开度变化和煤体变形的影响。Zeng等[21-22]引入气体侵入/排出体积与基质整体体积的比值来量化了局部膨胀/收缩到全局膨胀/收缩的过程。Huang等[23]假设最靠近裂隙的基质点压力与裂隙压力相同,使用该点膨胀应变与整体膨胀应变的比值来表示非平衡因子,以描述基质的非均匀变形。Liu等[24]根据通过薄膜的非稳态扩散解析解计算了基质中气体浓度的分布,定义了新的非平衡因子。这些研究都推出了更符合长期试验数据的理论成果。

    上述关于非平衡的研究在确定空间或时间变化因子方面存在很高的难度,或者对基质裂隙相互作用的力学解释不够充分。裂隙与基质的气体压差除了造成煤岩长期处于非平衡的吸附状态之外,还会造成基质和裂隙系统内有效应力的变化[25]。因此,本文考虑了气体压差引起的有效应力变化,并定义了一个新的与气体压差相关的非平衡因子,以描述非平衡演化过程中基质吸附变形对裂隙开度变化与煤体变形的影响。同时,本文考虑了基质桥的吸附变形,对基质块与基质桥吸附变形对裂隙开度变化的共同作用做出了解释。在此基础上,建立了一个气体压差影响下双重孔隙煤岩非平衡演化过程的渗透率理论模型,并进行了该理论模型控制下的煤岩试样注气与原位煤储层采气的数值模拟。研究结果对煤层气开采工程具有一定理论指导作用。

    一般情况下,气体在基质中的储存形式包含游离态和吸附态,而裂隙中只有游离态。即气体吸附仅发生在基质系统中,不发生在裂隙系统中。假设吸附引起的应变在各个方向的分量相等,双重孔隙煤岩变形弹性力学本构关系可以定义为[6]

    εij=12Gbσij(16Gb19Kb)σkkδij+α3Kbpmδij+β3Kbpfδij+εbs3δij (1)

    式中:下标m,f,b分别代表基质、裂隙与煤体;G为剪切模量,G=E/2(1+ν),E为杨氏弹性模量(MPa);ν为泊松比;K为体积模量,K=E/3(1-2ν);αβ分别为基质与裂隙的Biot系数,表达式分别为α=1-Km/ Ksβ=1-Kb/KmKs为煤基质骨架(煤颗粒)的体积模量;p为孔隙压力(MPa);总应力σkk=σ11+σ22+σ33,MPa;δij为Kronecker符号;εbs为基质吸附变形产生的煤体应变。根据Langmuir等温吸附方程,基质的吸附应变εms可以定义为

    εms=εLmpmpm+PLm (2)

    式中:εLmPLm分别为煤基质的Langmuir应变常数与压力常数。

    将式(1)与弹性力学平衡方程、几何方程联立,可得基质与裂隙系统的Navier型变形控制方程:

    Gbui,kk+Gb12vuk,kiαPm,iβPf,iKbεbs,i+Fi=0 (3)

    式中:u为位移矢量,F为体积力。

    向煤岩注气后,气体侵入裂隙壁附近的基质区域,引起基质局部膨胀,进而压缩裂隙,渗透率减小。随着气体向基质内部扩散,膨胀变形逐渐均匀,局部变形对裂隙开度的影响减弱,渗透率逐渐恢复。采气过程则与之相反。因此,引入一个数值介于0~1的非平衡因子Ue来量化基质吸附变形对裂隙开度变化的贡献程度。局部应变是由气体压差引起的基体与裂隙之间的应变差产生的[19],因此将该系数与气体压差建立起联系。

    Ue=|pfpm|max (4)

    Ue从0增大至1后再减小最终回到0,当Ue=0时,基质的吸附变形不会对裂隙开度变化产生影响,当Ue=1时,基质的吸附变形完全影响裂隙开度变化。 {p_{\text{f}}} {p_{\text{m}}} 适用于注气情况; {p_{\text{m}}} {p_{\text{f}}} 适用于采气情况。

    煤裂隙是气体的主要流动通道,而煤基质则为气体提供主要的储存空间,裂隙渗透率远大于基质渗透率[6-9],因此,此处只研究裂隙渗透率模型。煤裂隙孔隙率的表达式为

    \phi_{\mathrm{f}}=\frac{V_{\mathrm{f}}}{V_{\mathrm{b}}}。 (5)

    式中:Vf为裂隙体积;Vb为煤体体积。对式(5)两边取微分,可得裂隙孔隙率的改变量为

    {\text{d}}{\phi _{\text{f}}} = {\text{d}}\left( {\frac{{{V_{\text{f}}}}}{{{V_{\text{b}}}}}} \right) = {\phi _{\text{f}}}\left( {\frac{{{\text{d}}{V_{\text{f}}}}}{{{V_{\text{f}}}}} - \frac{{{\text{d}}{V_{\text{b}}}}}{{{V_{\text{b}}}}}} \right)。 (6)

    根据式(1)可得裂隙体应变增量与煤岩体应变增量分别为[26]

    \frac{\text{d}{V}_{\text{b}}}{{V}_{\text{b}}}=-\frac{1}{{K}_{\text{b}}}\left[\text{d}\overline{\sigma }-\alpha \text{d}{p}_{\text{m}}-\beta \text{d}{p}_{\text{f}}\right]+\text{d}{\varepsilon }_{\text{bs}}\text{ }。 (7)

    式中:\overline \sigma 为平均应力,\overline \sigma = - {\sigma _{kk}}/3。煤体、裂隙的变形均由两部分组成:一部分由有效应力变化引起,另一部分由气体吸附引起[27],类似的,裂隙应变增量则可表示为[25]

    \frac{\text{d}{V}_{\text{f}}}{{V}_{\text{f}}}=-\frac{1}{{K}_{\text{f}}}\left[\text{d}\overline{\sigma }-\beta \text{d}{p}_{\text{f}}-\text{d}p(t)\right]+\text{d}{\varepsilon }_{\text{fs}}\text{ }。 (8)

    式中:p(t)为裂隙与基质之间的压差表达式,其值随注气或采气过程的时间而变化,p(t)= {p_{\text{f}}} - {p_{\text{m}}} {\varepsilon _{{\text{fs}}}} 为由基质吸附变形与基质桥吸附变形共同产生的裂隙应变。方程(7),(8)等式右侧第一项为有效应力变化产生的应变增量,第二项为吸附应变增量。

    将式(7),(8)代入式(6)可得

    \begin{aligned} \frac{\mathrm{d} \phi_{\mathrm{f}}}{\phi_{\mathrm{f}}}=\left(\frac{1}{K_{\mathrm{b}}}-\right. & \left.\frac{1}{K_{\mathrm{f}}}\right)\left(\mathrm{d} \bar{\sigma}-\beta \mathrm{d} p_{\mathrm{f}}\right)- \\ & \frac{\alpha}{K_{\mathrm{b}}} \mathrm{~d} p_{\mathrm{m}}+\frac{1}{K_{\mathrm{f}}} \mathrm{~d} p(t)+\mathrm{d} \varepsilon_{\mathrm{fs}}-\mathrm{d} \varepsilon_{\mathrm{bs}} \end{aligned}。 (9)

    根据Betti-Maxwell互易定理可得Kf[28]

    {K}_{\text{f}}=\frac{{\varphi }_{\text{f}}}{\beta }{K}_{\text{b}}\text{ }。 (10)

    考虑基质-裂隙相互作用时,常用的煤岩理想化模型为基质桥模型[11]图 1为基质桥模型代表性单元,向煤岩注气时,基质吸附变形与基质桥吸附变形共同影响裂隙开度变化。图 1中深色区域表示吸附膨胀变形前的体积,浅色区域为吸附膨胀变形产生的体积差,红色箭头表示基质或基质桥的膨胀应变方向。基质桥与裂隙之间的接触比表面积较大,且基质桥自身体积较小,因此可以假设在注气后裂隙与基质桥之间的传质速度很快,即基质桥内的气体压力能在很短的时间内与裂隙压力达到平衡。基质由于体积较大,在与裂隙的传质过程中存在长期的非平衡态演化。因此,推导中假设基质桥中的气体吸附压力为与裂隙压力相等,基质中的气体压力则与裂隙压力存在压差。

    图  1  基质桥模型及吸附变形示意图
    Figure  1.  Schematic diagram of matrix bridge model and adsorption deformation

    基质桥是由于裂隙没有完全贯穿基质而产生的连接各基质块的部分,因此可以认为基质桥属于基质的一部分,其吸附特性与基质相同,可采用相同的Langmuir应变常数与压力常数。因此,吸附导致的基质桥膨胀应变为

    {\varepsilon }_{\text{ms}}^{\text{r}}={\varepsilon }_{\text{Lm}}\frac{{p}_{\text{f}}}{{p}_{\text{f}}+{P}_{\text{Lm}}}\text{ }。 (11)

    在注气过程中,基质桥吸附膨胀会增大裂隙开度,与基质吸附膨胀而压缩裂隙起到相反的作用。基于各向同性假设,由吸附膨胀导致的各个方向线应变为{\varepsilon _{{\text{s}}x}} = {\varepsilon _{{\text{s}}y}} = {\varepsilon _{{\text{s}}z}} = {\varepsilon _{\text{s}}}/3,以裂隙压缩方向为负,在水平方向,由基质吸附变形与基质桥吸附变形相互作用导致的裂隙开度变化为

    \Delta {b}_{\text{s}}=b\cdot \frac{\Delta {\varepsilon }_{\text{ms}}^{\text{r}}}{3}-{U}_{\text{e}}\cdot a\cdot \frac{\Delta {\varepsilon }_{\text{ms}}}{3}\text{ }。 (12)

    式中:b为裂隙开度(m);a为基质宽度(m)。水平方向上吸附导致的裂隙线应变增量为

    \Delta {\varepsilon }_{\text{fs}}^{x}\text{=}\frac{\Delta {b}_{\text{s}}}{b}=\frac{\Delta {\varepsilon }_{\text{ms}}^{\text{r}}}{3}-{U}_{\text{e}}\cdot \frac{a}{b}\cdot \frac{\Delta {\varepsilon }_{\text{ms}}}{3}\text{ }。 (13)

    因此,基质和基质桥吸附变形的共同影响的裂隙应变增量为

    \Delta {\varepsilon }_{\text{fs}}\text{=}3{\varepsilon }_{\text{fs}}^{x}=\Delta {\varepsilon }_{\text{ms}}^{\text{r}}-{U}_{\text{e}}\cdot \frac{a}{b}\cdot \Delta {\varepsilon }_{\text{ms}}\text{ }。 (14)

    由于基质吸附变形一部分导致裂隙开度变化,另一部分导致煤体变形,因此同理式(14),可得基质吸附变形导致的煤体应变增量为

    \Delta {\varepsilon _{{\text{bs}}}} = {\text{(}}1 - {U_{\text{e}}}{\text{)}} \cdot \frac{a}{{a + b}} \cdot {\varepsilon _{{\text{ms}}}}。 (15)

    由于裂隙开度b远小于基质宽度a,即a±ba,式(15)可简化为

    \Delta {\varepsilon }_{\text{bs}}=(1-{U}_{\text{e}})\cdot {\varepsilon }_{\text{ms}}\text{ }。 (16)

    联立式(2),(11),(14),(16)代入式(9),并对等式两侧进行积分、简化,得裂隙孔隙度模型为

    \begin{gathered} \frac{\phi_{\mathrm{f}}}{\phi_{\mathrm{f} 0}}=\exp \left\{\left(\frac{1}{K_{\mathrm{b}}}-\frac{1}{K_{\mathrm{f}}}\right)\left[\left(\bar{\sigma}-\bar{\sigma}_0\right)-\beta\left(p_{\mathrm{f}}-p_{\mathrm{f} 0}\right)\right]-\frac{\alpha}{K_{\mathrm{b}}}\left(p_{\mathrm{m}}-p_{\mathrm{m} 0}\right)+\right. \\ \frac{1}{K_{\mathrm{f}}}\left[p(t)-p\left(t_0\right)\right]+\varepsilon_{\mathrm{Lm}} P_{\mathrm{Lm}} \cdot\left[\frac{p_{\mathrm{f}}-p_{\mathrm{f} 0}}{\left(p_{\mathrm{f}}+P_{\mathrm{Lm}}\right)\left(p_{\mathrm{f} 0}+P_{\mathrm{Lm}}\right)}-\right. \\ \left.\left.\quad\left(U_{\mathrm{e}} \cdot \frac{a}{b}+1\right) \cdot \frac{p_{\mathrm{m}}-p_{\mathrm{m} 0}}{\left(p_{\mathrm{m}}+P_{\mathrm{Lm}}\right)\left(p_{\mathrm{m} 0}+P_{\mathrm{Lm}}\right)}\right]\right\}。 \end{gathered} (17)

    根据渗透率与孔隙率的立方关系

    \frac{{k}_{\text{f}}}{{k}_{\text{f}0}}={\left(\frac{{\phi }_{\text{f}}}{{\phi }_{\text{f}0}}\right)}^{3}\text{ }。 (18)

    将式(17)代入式(18),得裂隙渗透率比值为

    \begin{aligned} \frac{k_{\mathrm{f}}}{k_{\mathrm{f} 0}}= & \exp \left\{\left(\frac{3}{K_{\mathrm{b}}}-\frac{3}{K_{\mathrm{f}}}\right)\left[\left(\bar{\sigma}-\bar{\sigma}_0\right)-\beta\left(p_{\mathrm{f}}-p_{\mathrm{f} 0}\right)\right]-\frac{3 \alpha}{K_{\mathrm{b}}}\left(p_{\mathrm{m}}-p_{\mathrm{m} 0}\right)+\right. \\ & \frac{3}{K_{\mathrm{f}}}\left[p(t)-p\left(t_0\right)\right]+3 \varepsilon_{\mathrm{Lm}} P_{\mathrm{Lm}} \cdot\left[\frac{p_{\mathrm{f}}-p_{\mathrm{f} 0}}{\left(p_{\mathrm{f}}+P_{\mathrm{Lm}}\right)\left(p_{\mathrm{f} 0}+P_{\mathrm{Lm}}\right)}-\right. \\ & \left.\left.\left(U_{\mathrm{e}} \cdot \frac{a}{b}+1\right) \cdot \frac{p_{\mathrm{m}}-p_{\mathrm{m} 0}}{\left(p_{\mathrm{m}}+P_{\mathrm{Lm}}\right)\left(p_{\mathrm{m} 0}+P_{\mathrm{Lm}}\right)}\right]\right\}。 \end{aligned} (19)

    式(17),(19)即为本文裂隙孔隙率与渗透率模型的一般表达式。

    在实验室中,通常对煤岩试样应用恒围压边界条件[10, 18],即d\overline \sigma =0。根据式(19),并由Kb \gg Kf,可得恒围压条件下的渗透率模型为

    \begin{gathered} \frac{k_{\mathrm{f}}}{k_{\mathrm{f} 0}}=\exp \left\{\frac{3 \beta}{K_{\mathrm{f}}}\left(p_{\mathrm{f}}-p_{\mathrm{f} 0}\right)-\frac{3 \alpha}{K_{\mathrm{b}}}\left(p_{\mathrm{m}}-p_{\mathrm{m} 0}\right)+\frac{3}{K_{\mathrm{f}}}[p(t)-\right. \\ \\ \left.p\left(t_0\right)\right]+3 \varepsilon_{\mathrm{Lm}} P_{\mathrm{Lm}} \cdot\left[\frac{p_{\mathrm{f}}-p_{\mathrm{f} 0}}{\left(p_{\mathrm{f}}+P_{\mathrm{Lm}}\right)\left(p_{\mathrm{f} 0}+P_{\mathrm{Lm}}\right)}-\right. \\ \left.\left.\quad\left(U_{\mathrm{e}} \cdot \frac{a}{b}+1\right) \cdot \frac{p_{\mathrm{m}}-p_{\mathrm{m} 0}}{\left(p_{\mathrm{m}}+P_{\mathrm{Lm}}\right)\left(p_{\mathrm{m} 0}+P_{\mathrm{Lm}}\right)}\right]\right\} 。 \end{gathered} 。 (20)

    单轴应变条件通常被作为原位煤层的边界条件,在煤层气开采研究中被广泛应用[12, 26]。结合边界条件( {\varepsilon _{11}} = {\varepsilon _{22}} = 0 \Delta {\sigma _{33}} = 0 ),根据式(1)可得

    \begin{aligned} \Delta \sigma_{11}=\Delta \sigma_{22} & =\frac{v}{1-v} \Delta \sigma_{33}+ \\ & \frac{2 v-1}{1-v}\left(\alpha \Delta p_{\mathrm{m}}+\beta \Delta p_{\mathrm{f}}\right)-\frac{E_{\mathrm{b}}}{3(1-v)} \Delta \varepsilon_{\mathrm{bs}} \end{aligned} 。 (21)

    根据平均应力的表达式,可得平均应力增量为

    \overline{\sigma }-{\overline{\sigma }}_{0}\text{=}-\frac{1}{3}({\sigma }_{11}-{\sigma }_{110}+{\sigma }_{22}-{\sigma }_{220}+{\sigma }_{33}-{\sigma }_{330})\text{ }。 (22)

    将式(16),(21)代入式(22)可得

    \begin{array}{r} \bar{\sigma}-\bar{\sigma}_0=\frac{2-4 v}{3(1-v)}\left[\alpha\left(p_{\mathrm{m}}-p_{\mathrm{m} 0}\right)+\beta\left(p_{\mathrm{f}}-p_{\mathrm{f} 0}\right)\right]+ \\ \frac{2 E_{\mathrm{b}}}{9(1-v)}\left(1-U_{\mathrm{e}}\right)\left(\varepsilon_{\mathrm{ms}}-\varepsilon_{\mathrm{ms} 0}\right) \end{array}。 (23)

    将式(23)代入式(19),并由Kb \gg Kf,可得单轴应变条件下的裂隙渗透率模型为

    \begin{gathered} \frac{{{k_{\text{f}}}}}{{{k_{{\text{f}}0}}}} = \exp \left\{ {\frac{1}{{{K_{\text{f}}}}}\left[ {\frac{{1 + \nu }}{{1 - \nu }}\beta ({p_{\text{f}}} - {p_{{\text{f0}}}}) + 3(p(t) - p({t_0}))} \right] - \frac{\alpha }{{{K_{\text{b}}}}}\frac{{1 + \nu }}{{1 - \nu }} \cdot } \right. \hfill \\ {\text{ }}({p_{\text{m}}} - {p_{{\text{m}}0}}) + 3{\varepsilon _{{\text{Lm}}}}{P_{{\text{Lm}}}} \cdot \left[ {\frac{{{p_{\text{f}}} - {p_{{\text{f}}0}}}}{{({p_{\text{f}}} + {P_{{\text{Lm}}}})({p_{{\text{f}}0}} + {P_{{\text{Lm}}}})}}} \right. - \hfill \\ \left. {\left( {\frac{{2{E_{\text{b}}}(1 - {U_{\text{e}}})}}{{9{K_{\text{f}}}(1 - \nu )}} + \frac{a}{b}{U_{\text{e}}} + 1} \right) \cdot \frac{{{p_{\text{m}}} - {p_{{\text{m}}0}}}}{{({p_{\text{m}}} + {P_{{\text{Lm}}}})({p_{{\text{m0}}}} + {P_{{\text{Lm}}}})}}} \right\} \end{gathered}。 (24)

    气体在基质中的运移遵循Fick定律,在裂隙中的运移遵循Darcy定律。结合气体运移质量守恒定律,气体在基质与裂隙中的运移方程为[6]

    \frac{{\partial {c_{\text{m}}}}}{{\partial t}} + \nabla \cdot ( - D\nabla {c_{\text{m}}}) = - {Q_{\text{s}}}\text{,} (25)
    \frac{\partial {m}_{\text{f}}}{\partial t}+\nabla \cdot \left({\rho }_{\text{g}}\overrightarrow{{v}_{\text{g}}}\right)={Q}_{\text{s}}\text{ }。 (26)

    式中:cm为基质中的气体体积浓度(kg/m3);mf为裂隙中的气体含量(kg/m3);D为气体扩散系数(m2/s);ρg为气体密度(kg/m3);\overrightarrow {{v_{\text{g}}}} 为达西流速矢量(m/s);Qs为基质与裂隙之间的气体交换速率(kg/(m3·s))。

    基质中包含游离态和吸附态气体,气体吸附量可由Langmuir方程计算;而裂隙中只有游离态气体,所以在基质与裂隙中气体含量表达式分别为

    {c}_{\text{m}}={\rho }_{\text{gm}}{\phi }_{\text{m}}+{\rho }_{\text{c}}{\rho }_{\text{ga}}\frac{{p}_{\text{m}}{V}_{\text{Lm}}}{{p}_{\text{m}}+{P}_{\text{Lm}}}\text{ }\text{,} (27)
    {m}_{\text{f}}={\rho }_{\text{gf}}{\phi }_{\text{f}}\text{ }。 (28)

    式中:ρc为煤岩的密度(kg/m3);ρga为标准状况下的气体密度(kg/m3);VLm为Langmuir体积常数(m3/kg)。

    根据理想气体状态方程可得气体密度表达式为

    {\rho }_{\text{g}}=\frac{M}{RT}p\text{ }。 (29)

    式中:M为气体的摩尔质量(kg/mol);R为摩尔气体常数(J/(mol·K));T为温度(K)。

    基质与裂隙之间的气体交换速率表达式为

    {Q_{\text{s}}} = D \cdot w({\rho _{{\text{gm}}}} - {\rho _{{\text{gf}}}}) = - \frac{{DwM}}{{RT}}p{\text{(}}t{\text{)}}。 (30)

    式中:w为形状因子(m-2)。

    裂隙中气体达西流速矢量表达式为

    \overrightarrow{{v}_{\text{g}}}=-\frac{{k}_{\text{f}}}{\mu }\nabla {p}_{\text{f}}\text{ }。 (31)

    式中:μ为气体的动力学黏度(Pa·s)。

    联立式(25)~(31),可得煤基质与裂隙中的气体运移耦合控制方程

    \begin{aligned} & {\left[\phi_{\mathrm{m}}+\frac{\rho_{\mathrm{c}} p_{\mathrm{a}} V_{\mathrm{Lm}} P_{\mathrm{Lm}}}{\left(p_{\mathrm{m}}+P_{\mathrm{Lm}}\right)^2}\right] \frac{\partial p_{\mathrm{m}}}{\partial t}+p_{\mathrm{m}} \frac{\partial \phi_{\mathrm{m}}}{\partial t}+} \\ & \nabla \cdot\left[-D\left[\phi_{\mathrm{m}}+\frac{\rho_{\mathrm{c}} p_{\mathrm{a}} V_{\mathrm{Lm}} P_{\mathrm{Lm}}}{\left(p_{\mathrm{m}}+P_{\mathrm{Lm}}\right)^2}\right] \nabla p_{\mathrm{m}}\right]=D w \Delta p(t), \end{aligned} (32)
    {\phi }_{\text{f}}\frac{\partial {p}_{\text{f}}}{\partial t}+{p}_{\text{f}}\frac{\partial {\phi }_{\text{f}}}{\partial t}+\nabla \cdot \left(-{p}_{\text{f}}\frac{{k}_{\text{f}}}{\mu }\nabla {p}_{\text{f}}\right)=-Dw\Delta p(t)\text{ }。 (33)

    煤层气开采过程中气体流动与扩散会导致气体压力变化,进而使煤岩发生变形,而煤岩变形会引起其孔隙率与渗透率的改变,又反过来影响气体流动与扩散。这些物理场的交叉耦合关系如图 2所示。涉及多物理场耦合的煤岩渗透率演化过程可于COMSOL Multiphysics中进行数值模拟求解。

    图  2  煤岩流固耦合关系示意图
    Figure  2.  Diagram of fluid-solid coupling relationship in coal

    在过去关于基质与裂隙相互作用的研究中,建立的渗透率模型往往能够顺利匹配“渗透率-压力”试验数据[12-14],但无法成功匹配“渗透率-时间”试验数据。这是因为这些模型通过拟合出一个定值f来描述基质与裂隙相互作用,而实际上基质吸附变形对裂隙开度变化的贡献程度与时间有关,定值f无法描述非平衡演化过程。本部分使用COMSOL模拟煤试样注气过程,以得到式(20)控制下的渗透率演化预测结果。将预测曲线与Wei等[18]试验获取的长期渗透率数据进行对比,以验证本文渗透率模型的合理性。该试验在恒围压条件下进行,根据岩芯形状与其受力的几何对称性,在数值模拟中可将圆柱形煤试样简化为二维轴对称模型,如图 3所示。几何模型底部边界为辊支撑,其余边界设置6 MPa恒定压力,上边界设置3 MPa的注入压力,煤试样内初始气体压力为0.1 MPa。

    图  3  实验室注气模拟模型
    Figure  3.  Simulation model for gas injection in laboratory

    煤体弹性模量E取3 GPa,煤基质弹性模型Em取10 GPa[20],泊松比 \nu 取0.34,则可分别计算得到煤体体积模量Kb与煤基质体积模量Km;煤裂隙体积模量Kf由式(10)计算;煤基质骨架(颗粒)体积模量Ks可使用下式计算[31]

    {K_{\text{s}}}{\text{ = }}\frac{{{K_{\text{m}}}}}{{1 - 3{\phi _{{\text{m0}}}}(1 - \nu )/\left[ {2(1 - 2\nu )} \right]}} 。 (34)

    根据各类体积模量可计算得到Biot系数αβ的值。其他相关参数取自Wei等[20]的研究,将模拟所需的所有参数汇总于表 1

    表  1  用于理论模型验证的参数
    Table  1.  Parameters used for verification of theoretical model
    变量 物理意义
    E 3000 MPa 煤岩弹性模量
    Em 10000 MPa[20] 煤基质弹性模量
    \nu 0.34 泊松比
    ρc 1400 kg/m³ 煤岩密度
    ρga 1.98 kg/m³ 标况下CO2密度
    M 0.044 kg/mol CO2的摩尔质量
    μ 1.84×10-5 Pa·s CO2动力黏度
    R 8.414 J/(mol·K) 摩尔气体常数
    T 293 K 温度
    pa 0.1 MPa 大气压
    pinj 3 MPa[20] 注入压力
    σ0 6 MPa[20] 围压
    D 2×10-8 m²/s[20] CO2扩散系数
    w 12 m-2[20] 形状因子
    εLm 0.0518[20] Langmuir应变常数
    VLm 0.017 m3/kg[20] Langmuir体积常数
    PLm 4.48 MPa[20] Langmuir压力常数
    φf0 1.2%[20] 初始裂隙孔隙率
    φm0 2%[20] 初始基质孔隙率
    kf0 1×10-17 m2[20] 初始裂隙渗透率
    下载: 导出CSV 
    | 显示表格

    模型验证结果如下图 4所示,并与Lu等[12]建立的模型进行比较。Lu等[12]模型考虑了基质裂隙相互作用,f取值为1和0时可以表示渗透率演化的上下限。f=1时,基质吸附变形完全影响裂隙开度的变化,此时的预测曲线能够成功匹配前期的试验数据,这是因为在前期基质发生局部变形,主要影响裂隙开度变化。f=0时,基质吸附变形对裂隙开度变化没有影响,此时的预测曲线则在后期接近试验数据,这是因为后期逐渐接近最终平衡状态,煤岩发生整体变形,对裂隙渗透率影响不大。但Lu等[12]模型无法描述从初始平衡至最终平衡演化的过程。相比之下,本文模型可以很好地匹配整个试验过程的渗透率数据,有效解释了非平衡演化过程的渗透率趋势。

    图  4  理论模型验证与对比
    Figure  4.  Validation of theoretical model and its comparison

    煤层产气过程与气体注入煤样过程相反,但渗透率的演化机制相同[19]。因此,本文理论模型也可用于解释煤层气开采过程中渗透率的演化。本部分采用COMSOL研究式(24)控制下的煤层气开采过程的渗透率演化规律。如图 5所示,建立一个200 m×200 m的煤层几何模型,中心位置为坐标原点。模型下边界设置为固定约束,左右边界设置为辊支撑,上边界设置竖直向下的恒定边界荷载Fy,以此来模拟单轴应变条件。

    图  5  煤层气开采的模拟模型
    Figure  5.  Simulation model for extraction of coal bed methane

    煤层初始状态下的基质孔隙压力与裂隙压力相等,均设为4 MPa[24]。在中心位置处钻一半径为0.5 m的井口,井口边界压力为0.3 MPa[24]。取A点(50 m,50 m)以研究采气过程中气体压力与渗透率的演化规律。煤层气开采数值模拟中所用的参数列于表 2,其他参数与表 1所列相同。

    表  2  用于现场开采煤层气数值模拟的参数
    Table  2.  Parameters used for numerical simulation of CBM
    符号 物理意义
    ρga 0.717 kg/m³ 标况下CH4密度
    M 0.016 kg/mol CH4的摩尔质量
    μ 1.15×10-5 Pa·s CH4动力黏度
    D 2×10-9 m²/s[20] CH4扩散系数
    pbh 0.3 MPa[24] 井底压力
    p0 4 MPa[24] 储层初始压力
    εLm 0.01[24] Langmuir应变常数
    VLm 0.017 m3/kg[24] Langmuir体积常数
    PLm 2.5 MPa[24] Langmuir压力常数
    kf0 2 mD[24] 初始裂隙渗透率
    下载: 导出CSV 
    | 显示表格

    图 6绘制了气体压力随时间的分布规律。开采前,A点的基质孔隙压力等于裂隙压力,均为初始值4 MPa。钻井后进行开采时,由于裂隙渗透率远大于基质渗透率,裂隙压力首先迅速降低,初始平衡状态被打破,煤储层开始非平衡态演化。而基质孔隙压力在103 h时才出现较为明显的下降趋势,且在后续的很长时间的采气过程中均大于裂隙压力。相应地,非平衡因子Ue随着气体压差的变化从0迅速增大后缓慢降低。在106 h之后才恢复至0,煤储层到达最终平衡状态。非平衡状态的演化持续了很长的时间。

    图  6  气体压力与非平衡因子演化规律
    Figure  6.  Evolution of gas pressure and non-equilibrium factor

    图 7描述了与井点同一水平位置处(x轴处)的气体压力空间分布规律。对比图 7(a),(b)可知,裂隙压力与基质孔隙压力均从远离井点的位置向井点处不断降低,横坐标x=0处的气体压力为井底压力0.3 MPa。但裂隙压力随空间的降低幅度大于基质孔隙压力的降低幅度。这是因为裂隙中的气体不断从远离井点处向近井点处运输,近井点处的气体更快排出井外,因此近井点处的裂隙压力较低。而基质孔隙中的气体则主要是向裂隙中传质,传质过程发生在储层的各个位置,因此基质孔隙压力曲线趋于水平。图 6显示在非平衡演化过程中,储层A点处的裂隙压力在任意时间均小于基质孔隙压力。图 7则显示在非平衡演化过程中,裂隙压力在同一时间的任意位置均小于基质孔隙压力。

    图  7  气体压力的空间分布
    Figure  7.  Spatial distribution of gas pressure

    为了分析采气过程中的裂隙渗透率演化规律,绘制了其随时间的演化曲线,由图 8可知,渗透率演化同样可分为5个阶段,经历了先减小后增大再减小的趋势。其演化原理与注气过程相同,但过程相反。阶段Ⅰ代表初始平衡状态,裂隙压力与基质孔隙压力均为初始值。阶段Ⅱ中,采气使裂隙压力迅速减小,裂隙渗透率因有效应力的增大而减小。同时基质桥产生解吸收缩会对裂隙两侧的基质产生一种类似“牵拉”的力,同样会使裂隙开度减小。在此阶段,初始平衡状态被打破,非平衡演化过程开始。阶段Ⅲ中,由于裂隙与基质之间的气体压差较大,气体解吸发生在靠近裂隙的基质区域,产生的局部收缩使裂隙开度增大,渗透率增大。阶段Ⅳ中,随着基质内部气体的排出,裂隙与基质的气体压差减小,基质局部收缩向整体收缩演化,基质吸附变形从主要增大裂隙开度变为主要减小煤体体积,渗透率逐渐降低。阶段Ⅴ中,基质孔隙压力与裂隙压力再次相等,煤岩进入最终平衡状态。因此,在开采前期,有效应力与.基质桥的收缩对渗透率的影响占主导作用;随着开采的进行,主导渗透率演化的因素变为基质局部收缩与整体收缩效应。

    图  8  煤层气开采时渗透率长期演化规律
    Figure  8.  Long-term evolution of permeability in CBM

    图 9显示了裂隙渗透率随孔隙压力的变化规律。采气过程为降压过程,孔隙压力从4 MPa降至0.3 MPa的整个过程中,渗透率比值先降低后升高再降低,其随裂隙压力和基质孔隙压力变化的趋势相同,但变化位置不同。渗透率比值从1降至最小值时,基质孔隙压力几乎不变,而裂隙压力降低至1.5 MPa左右。渗透率比值从最小值增大到最大值时,基质孔隙压力降低至0.75 MPa左右,裂隙压力降低至0.4 MPa左右。渗透率比值从最大值降至最终平衡状态下的比值时,基质孔隙压力与裂隙压力从不同的压力值分别降低至0.3 MPa。由图 9也可看出,基质孔隙压力的降低始终滞后于裂隙压力,且其二者共同决定渗透率的演化过程。值得注意的是,根据图 6的注气曲线可得,最终平衡状态下的渗透率因裂隙压力的升高而大于初始平衡状态下的渗透率,因此图 89采气曲线的最终平衡状态下的渗透率应该由于裂隙压力的降低而小于初始平衡状态下的渗透率,但曲线却呈现出了相反的结果,这可能是由适用于现场的单轴应变条件导致。图 10绘制了不同条件下的渗透率演化规律,针对这一现象做出了解释。

    图  9  煤层气开采时渗透率随孔隙压力的变化
    Figure  9.  Variation of permeability with pore pressure in CBM
    图  10  不同情况下的渗透率演化规律
    Figure  10.  Evolution of permeability under different conditions

    图 10可得,恒围压条件下的和单轴应变条件下的渗透率演化趋势相同,但变化幅度不同,且最终平衡状态下的渗透率大小也不同。在恒围压条件下,最终平衡状态的渗透率小于初始平衡状态的渗透率。而单轴应变条件下,最终平衡状态的渗透率大于初始平衡状态的渗透率,这是因为单轴应变条件限制了xy方向的应变,使得基质收缩更多的影响了裂隙开度变化而不是煤整体变形,因此渗透率变化幅度更大。此外,对比单轴应变条件下不同基质宽度与基质桥宽度之比(a/b)的3条曲线可得,a/b值不会影响最终平衡状态下的渗透率,但会影响渗透率的变化幅度。a/b越大,基质占比越大,基质桥占比越小,基质吸附变形带来的影响越大,基质桥对基质变形的阻碍作用越小,渗透率变化的幅度越大。

    建立了一个考虑裂隙与基质气体压差和基质桥吸附变形的煤岩非平衡演化全过程渗透率模型;模型中定义了一个与气体压差相关的非平衡因子,以表征煤岩抽注过程中基质吸附变形在不同时刻对裂隙开度变化的影响程度;使用渗透率长期演化的试验数据对该模型进行了验证;并采用数值模拟进一步分析了原位开采煤层气时的气体压力时空分布与渗透率演化规律,具体得到以下3点结论。

    (1)煤岩抽注过程中基质孔隙压力的变化总是滞后于裂隙压力的变化,裂隙与基质之间则会产生气体压差,进而造成煤结构的变形。这种变形既包括因有效应力改变而产生的力学变形,也包括非平衡状态下的基质非均匀吸附变形。

    (2)煤试样注气和煤层采气过程的渗透率长期演化的原理相同,过程相反。其演化可分为5个阶段,首尾两阶段代表初始平衡与最终平衡状态,中间3个阶段代表非平衡状态,非平衡演化过程中煤岩渗透率受基质非均匀吸附变形的影响而发生非单调变化。

    (3)用于煤储层的不同边界条件会影响渗透率演化曲线的变化幅度和最终数值,单轴应变条件下的渗透率变化幅度和最终数值均大于恒定围压条件下的渗透率变化幅度和最终数值。不同基质块与基质桥宽度之比会影响渗透率曲线的变化幅度,但不影响其最终数值。基质块与基质桥宽度之比越大,渗透率变化幅度增大。

  • 图  1   土体参数的直方图、散点图以及相关系数

    Figure  1.   Histogram, scatter plot and correlation coefficients of soil parameters

    图  2   压缩模量的更新结果

    Figure  2.   Updating results of compressive modulus

    图  3   压缩模量的预测均值与实测值对比

    Figure  3.   Comparison of predicted and measured values of compressive modulus

    表  1   9个土体参数的统计信息

    Table  1   Statistics of 9 soil parameters

    参数n均值最小值最大值STDCOV
    w/%88565.1139.4098.6010.250.16
    e8851.751.002.620.270.15
    wL88556.5432.9087.309.290.16
    IP88525.4912.8044.705.320.21
    IL8851.350.942.580.260.19
    α 7421.560.533.740.530.34
    Es 7421.750.773.520.490.28
    σv 88598.305.00368.9469.970.71
    Su 21610.532.1054.007.130.68
    下载: 导出CSV

    表  2   分布类型与参数

    Table  2   Distribution types and parameters

    参数类型aXbXaYbYp
    wSU0.850.015.9165.190.24
    eSU1.06-0.190.221.700.48
    wLSB1.190.6553.2936.300.67
    IPSB1.430.7134.4512.090.21
    ILSB0.730.891.201.010.06
    α SU2.18-2.350.630.650.89
    Es SU1.65-1.020.611.280.97
    σv SB0.981.20394.85-7.320.77
    Su SB0.991.0226.441.860.99
    下载: 导出CSV
  • [1]

    CHING J, PHOON K K. Transformations and correlations among some clay parameters—the global database[J]. Canadian Geotechnical Journal, 2014, 51(6): 663-685. doi: 10.1139/cgj-2013-0262

    [2]

    PHOON K K, KULHAWY F H. Evaluation of geotechnical property variability[J]. Canadian Geotechnical Journal, 1999, 36(4): 625-639. doi: 10.1139/t99-039

    [3]

    PHOON K K, KULHAWY F H. Characterization of geotechnical variability[J]. Canadian Geotechnical Journal, 1999, 36(4): 612-624. doi: 10.1139/t99-038

    [4] 高大钊. 地基土力学性质指标的可靠性分析与取值[J]. 同济大学学报, 1985, 13(4): 59-68. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ198504007.htm

    GAO Da-zhao. Reliability analysis and evaluation of soil-property parameters[J]. Journal of Tongji University, 1985, 13(4): 59-68. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ198504007.htm

    [5] 姚耀武, 陈东伟. 土坡稳定可靠度分析[J]. 岩土工程学报, 1994, 16(2): 80-87. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC402.009.htm

    YAO Yao-wu, CHEN Dong-wei. Reliability analysis of slope stability[J]. Chinese Journal of Geotechnical Engineering, 1994, 16(2): 80-87. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC402.009.htm

    [6] 倪红, 刘新宇, 秦玉. 土性参数概率特性对地基承载力可靠度的影响[J]. 解放军理工大学学报(自然科学版), 2004, 5(3): 67-69. https://www.cnki.com.cn/Article/CJFDTOTAL-JFJL200403016.htm

    NI Hong, LIU Xin-yu, QIN Yu. Reliability analysis of bearing capacity of soil[J]. Journal of PLA University of Science and Technology (Natural Science), 2004, 5(3): 67-69. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JFJL200403016.htm

    [7] 顾晓强, 吴瑞拓, 梁发云, 等. 上海土体小应变硬化模型整套参数取值方法及工程验证[J]. 岩土力学, 2021, 42(3): 833-845. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202103026.htm

    GU Xiao-qiang, WU Rui-tuo, LIANG Fa-yun, et al. On HSS model parameters for Shanghai soils with engineering verification[J]. Rock and Soil Mechanics, 2021, 42(3): 833-845. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202103026.htm

    [8]

    CHING J, PHOON K K. Modeling parameters of structured clays as a multivariate normal distribution[J]. Canadian Geotechnical Journal, 2012, 49(5): 522-545. doi: 10.1139/t2012-015

    [9]

    CHING J, PHOON K K. Correlations among some clay parameters—the multivariate distribution[J]. Canadian Geotechnical Journal, 2014, 51(6): 686-704. doi: 10.1139/cgj-2013-0353

    [10]

    ZOU H F, LIU S Y, CAI G J, et al. Multivariate correlation analysis of seismic piezocone penetration (SCPTU) parameters and design properties of Jiangsu quaternary cohesive soils[J]. Engineering Geology, 2017, 228: 11-38. doi: 10.1016/j.enggeo.2017.07.005

    [11]

    PHOON K K. Modeling and simulation of stochastic data[C]//GeoCongress 2006. Atlanta, 2006: 1-17.

    [12]

    PHOON K K, CHING J. Multivariate model for soil parameters based on Johnson distributions[C]//Geo- Congress 2013. 2013, San Diego.

    [13]

    SLIFKER J F, SHAPIRO S S. The Johnson system: selection and parameter estimation[J]. Technometrics, 1980, 22(2): 239-246.

    [14]

    HÄRDLE W K, SIMAR L. Multivariate distributions[M]//Applied Multivariate Statistical Analysis. Berlin: Springer Berlin Heidelberg, 2012: 107-165.

图(3)  /  表(2)
计量
  • 文章访问数:  286
  • HTML全文浏览量:  29
  • PDF下载量:  179
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-08-15
  • 网络出版日期:  2022-12-05
  • 刊出日期:  2021-10-31

目录

/

返回文章
返回