考虑土塞效应的开口管桩沉桩与承载全过程离散元分析

李立辰; 刘卓; 刘浩; 吴文兵; 罗仑博; 蒋国盛; 梅国雄

doi:10.11779/CJGE20230340

考虑土塞效应的开口管桩沉桩与承载全过程离散元分析 English Version

李立辰^1,,
刘卓¹,
刘浩¹,
吴文兵^1, ,,
罗仑博²,
蒋国盛¹,
梅国雄^{1, 3}

1.
中国地质大学(武汉)工程学院岩土钻掘与防护教育部工程研究中心, 湖北武汉 430074
2.
中国长江三峡集团有限公司科学技术研究院，北京 101100
3.
浙江大学海洋学院，浙江舟山 316021

基金项目:

国家自然科学基金项目 52178371

国家自然科学基金项目 52308383

湖北省自然科学基金青年项目 2023AFB008

中国博士后科学基金面上项目 2022M712964

详细信息

作者简介:
李立辰(1994—)，男，博士，助理研究员，主要从事桩基工程理论与技术方面的教学与科研工作。E-mail: 20121003645@cug.edu.cn

通讯作者:
吴文兵, E-mail: zjuwwb1126@163.com

中图分类号: TU470
计量
- 文章访问数: 474
- HTML全文浏览量: 65
- PDF下载量: 144
出版历程
- 收稿日期: 2023-04-18
- 网络出版日期: 2023-11-29
- 刊出日期: 2024-06-30

DEM analysis of installation and bearing process of open-ended piles considering plugging effects

1.
Faculty of Engineering, Engineering Research Centre of Rock-Soil Drilling & Excavation and Protection, Ministry of Education, China University of Geosciences, Wuhan 430074, China
2.
Institute of Science and Technology, China Three Gorges Corporation, Beijing 101100, China
3.
Ocean College, Zhejiang University, Zhoushan 316021, China

摘要

摘要: 开口管桩贯入过程中会产生复杂的土塞效应，显著影响管桩的沉桩过程与承载特性。考虑到试验方法成本高、周期长、且难以直观反映土塞与管桩的相互作用特性，基于离散元方法开展砂土地基中静压开口管桩沉桩与承载全过程模拟分析。基于室内试验数据对离散元模型参数进行标定，并采用周期性单元复制法生成均质且应力沿深度线性分布的地基模型。基于土体变形特征，揭示管桩中土塞的主动拱与被动拱分布特征。进一步地，结合土塞参数的演化规律，得出土塞增长率（IFR）对土塞状态变化更为敏感，且土塞阻力的发挥趋势与IFR值变化趋势有较高的相似性。基于管桩静载试验模拟，得出静载下开口管桩呈完全闭塞状态，土塞的荷载-位移响应呈弹塑性，并最终建立考虑CPT锥尖阻力、土塞效应、管桩尺寸的静压开口管桩桩端阻力经验计算公式。本文研究成果揭示了土塞的形成演化过程与承载机制，对于考虑土塞效应的开口管桩沉桩阻力预测与承载力计算具有一定参考价值。
- 开口管桩 /
- 土塞效应 /
- 离散元法 /
- 桩端阻力
Abstract: The penetration process of open-ended piles will cause complex plugging effects, which further influence the installation and bearing characteristics of piles. Considering the high cost, long period and limitations in reflecting the interaction between piles and soil plugs of the traditional experimental approach, the jacking installation and subsequent loading process of open-ended piles in sand are simulated based on the discrete element method (DEM). The periodic cell replication method is used to generate a homogeneous foundation sample with linear stress distribution. Based on the analysis of the soil deformation characteristics, the characteristics of the active arch and passive arch distribution of the soil plugs are revealed. Based on the analysis of the evolution of soil plugging parameters and plug bearing characteristics, it is found that the incremental filling ratio (IFR) is more sensitive to the plug behavior, and the development of plug resistance shows good consistency with the trend of the IFR value. Also, the open-ended piles are in the fully plugged mode under the static loads, and the load-displacement response of soil plugs is elastic-plastic. Finally, the method for calculating the end bearing capacity of open-ended jacked piles is established by considering the CPT cone resistances, plugging effects and pile size effects. The above results reveal the evolution process and bearing mechanism of soil plugs and can provide reference value for the prediction of installation resistance and calculation of bearing capacity of open-ended piles considering soil plugging effects.
- open-ended pile /
- plugging effect /
- discrete element method /
- end bearing capacity

HTML全文

0. 引言

随着全球经济的迅速发展,人类对能源的需求日益增加。煤层气作为一种清洁能源,正逐步取代传统常规能源成为人们生活中必不可少的一部分^[1]。煤是典型的双重孔隙结构,由裂隙系统和基质系统组成。裂隙系统为气体流动的主要通道,基质系统为主要的储气空间。在煤层气开采过程中,气体压力降低将会引起有效应力的增加,从而压缩裂隙来减少储层的采收率。与此同时,煤基质发生解吸,使得基质收缩导致裂隙开度的增加,进而提高采收率^[2]。长期的煤层气开采取决于这种基质-裂隙相互作用^[3-4]。因此,正确认识基质-裂隙相互作用过程对于预测煤层气开采至关重要。

由于煤储层复杂的双重孔隙结构及“低渗透”、“强吸附”的特点,导致煤层渗透率难以预测。目前大量的渗透率模型依据吸附膨胀和有效应力来表征渗透率的演化,这些模型通常认为气体吸附引起的应变完全由裂隙开度闭合所适应。Gray^[5]在渗透率模型建立中首次考虑了煤基质的膨胀和收缩,从而对煤层储层渗透率的动态波动进行了解释。随后,Seidle等^[6]提出了火柴棍几何模型,概念性地解释了煤渗透率随应力增加而减少这一现象,并在之后的研究中用实验测量了甲烷和二氧化碳的收缩系数,进而推导了与吸附/收缩有关的渗透率模型^[7]。Palmer等^[8]（P-M模型）通过直接类比温度对弹性多孔介质的影响,将基质收缩作为储层压力的函数,建立了一种结合有效应力和基质收缩效应计算孔隙体积压缩性和渗透率的理论模型。Shi等^[9]（S-D模型）同样采用火柴棍几何模型的概念,在Gray模型的基础上认为基质收缩与解吸气体体积而不是等效吸附压力成正比,这一点与Seidle等^[7]的观点相同,并在之后提出了一种将解吸气体含量与体积基质应变联系起来的公式,从而更新了渗透率模型^[10]。不同于S-D模型的理论推导,Cui等^[11]（C-B模型）以Piceance盆地煤的等温吸附和体积应变为约束条件,通过实验量化储层压力和气体吸附引起的体积应变对煤层渗透率的影响,建立起一个应力依赖性的渗透率模型。Connell^[12]在现有概念模型的基础上,增加了孔压、围压和吸附应变的独立效应,对渗透率的指数模型进行了扩展。

尽管上述模型对储层渗透率都有较好的预测,但是在研究吸附引起的应变影响时,通常认为基质吸附引起的变形完全挤压裂隙开度,这可能会高估吸附膨胀的影响。Robertson^[13]首先提出运用应变修正系数对模型进行修正,提高了P-M模型和S-D模型与数据的拟合效果。Liu等^[14]引入了内部膨胀系数的概念,在单轴应变和恒定围压条件下建立了一个新的渗透率模型,认为只有部分吸附应变会引起裂隙开度的变形。在此之后,Chen等^[3]对于膨胀应变效应定义了分配比来考虑基质-裂隙之间的相互作用,从而合理解释了恒定有效应力条件下的渗透率实验。Guo等^[15]在模型建立中同样认为只有一部分基质变形对裂隙变形有贡献,并定量测试了采用N₂和CO₂测试Sulcis煤样的基质变形贡献。Zhou等^[16]考虑吸附膨胀对裂隙开度的影响时,认为内膨胀系数是个常数,结合变Klinkenberg效应和裂隙开度建立起一个新的渗透率模型。Liu等^[17]重新定义了内部膨胀系数f,考虑基质和裂隙之间气体压力的差异来减少评估基质裂隙相互作用时的错误,并对系数f的控制机理和影响因素进行了系统的研究。Jiang等^[18]则系统性的将渗透率模型分为4类：①吸附应变只影响裂隙应变；②吸附应变只影响煤体应变；③吸附引起体应变等于裂隙应变；④吸附引起基质应变、裂隙应变和体应变之比为常数,并着重讨论了第④种分类,将经验公式扩展到理论模型,从而明确的考虑了不同应变的影响。

除了基质裂隙相互作用外,应力对储层变形也有很大的影响^[19]。一般认为,煤储层受自身重力、裂隙孔隙压力、基质孔隙压力以及吸附引起的体应力的影响。基质吸附应变一部分造成裂隙开度的变形,其余部分对裂隙煤层的体应变产生影响^{[14, 17, 20]},这导致只有部分吸附变形对体应力产生影响。然而,上述模型忽略了对吸附应力的修正,因此在对储层渗透率进行预测时会产生一定的误差。此外,基质和裂隙在开采过程中是动态变化的^[21],但前人的模型中常常仅考虑有效应力及吸附膨胀对裂隙开度的影响,把基质当成刚体,忽略了基质本身在有效应力作用下的变形。因此仅考虑裂隙开度变化的渗透率模型难以准确地对储层的渗透率进行预测。

不可否认,之前的研究极大的提高了对于非常规储层气体运移过程的理解。然而,在建立模型预测储层渗透率时仍然存在着一些问题没有解决。在本文中推导了一个新的渗透率模型从3个方面去解决这些问题：①同时考虑基质和裂隙的变形,建立一个可以反映煤层气开采过程中动态变化的模型；②采用基质膨胀部分（f）造成裂隙开度变化,其余部分（1-f）产生煤岩体应变的假设,在传统的应力控制方程吸附项前引入系数（1-f）对外部应力进行修正；③考虑基质和裂隙之间的压力差异来减少对基质-裂隙相互作用评估的误差。此外,通过对新模型进行简化,得到仅考虑裂隙开度变化的渗透率模型以及未考虑应力修正的渗透率模型,分别用现场数据和实验室数据在不同边界条件下对这3个模型进行相互比较,分析考虑基质本身变形及应力修正对渗透率演化的影响。并且进一步地将不同边界条件下的新模型与4种传统的模型和考虑内膨胀系数的模型进行比较,验证了新模型对渗透率数据预测的优越性。

1. 模型的建立

在立方定律中,裂隙的渗透率k_f可以用基质宽度s和裂隙开度b来定义^[22]：

$k_{f} = \frac{b^{3}}{12 s}$ 。

(1)

式（1）表示的渗透率计算方式可直接由已知的裂隙开度和基质宽度进行计算,然而在开采过程中,由于复杂的应力作用和吸附膨胀,裂隙开度和基质宽度是动态变化的,因此为了反映变应力状态下储层真实的渗透率,将裂隙开度和基质宽度当作变量,对式（1）进行求导可得如下所示的储层在变应力状态下渗透率的导数表达式：

$d k_{f} = k_{f} (\frac{3 Δ b}{b} - \frac{Δ s}{s})$ 。

(2)

式（2）表明裂隙的渗透率不仅取决于裂隙开度的变化,基质变形同样会影响渗透率的值。因此只单独考虑裂隙开度的变化,进而直接利用立方定律得到的渗透率模型在预测储层开采量时会存在一定的误差。下面分别考虑基质和裂隙在有效应力和吸附膨胀作用下产生的变形,从而建立一个更加精确的储层渗透率模型。

1.1 力学变形

由于煤岩基质的不均质性及致密性,通常将其理想化为以基质为代表的中、小孔隙和以裂隙为代表的大孔隙的双重孔隙度介质^[23-24]。由于裂隙和基质的导流能力差异很大,导致气体在裂隙中的流动远远高于基质内的扩散速率,进而造成基质和裂隙在相同的外部应力作用下,承受不同的有效应力^[25]。然而,在之前的模型中,尽管认识到基质和裂隙中不同的孔隙压力,却往往忽略基质本身的力学变形^[17]。本文将裂隙开度和基质宽度视为变量,在模型建立中同时考虑基质和裂隙在力学作用下产生的变形,假设基质和裂隙压缩为负,膨胀为正,则基质产生的力学变形可由有效应力原理和胡克定律得到

$Δ s_{e} = \frac{- s Δ {σ^{'}}_{m}}{K_{m}} = - \frac{s (Δ σ - α Δ p_{m})}{K_{m}}$ ,

(3)

式中, $Δ s_{e}$ 为基质在有效应力作用下产生的变形, ${σ^{'}}_{m}$ 为基质承受的有效应力, $K_{m}$ 为基质的体积模量, $σ$ 为外部应力, $p_{m}$ 为基质内的孔隙压力, $α$ 定义为基质的Biot系数,

$α = 1 - \frac{K}{K_{m}}$ 。

(4)

式中K为煤块的体积模量,K=E/3(1-2ν)；E为弹性模量； $ν$ 为泊松比。同理裂隙在有效应力作用下产生的变形为

$Δ b_{e} = \frac{- b Δ {σ^{'}}_{f}}{K_{f}} = - \frac{b (Δ σ - β Δ p_{f})}{K_{f}}$ ,

(5)

其中, ${σ^{'}}_{f}$ 为裂隙承受的有效应力,K_f为裂隙的体积模量, $p_{f}$ 为裂隙内的孔隙压力, $β$ 为裂隙的Biot系数,

$β = 1 - \frac{K_{f}}{K_{m}}$ 。

(6)

1.2 吸附产生的变形

煤岩作为一种双重孔隙介质,吸附只发生在基质中^[26]。由于裂隙中的破碎岩屑和压裂液中的支撑剂影响,典型的裂隙包括两个表面接触的粗糙区,如图1（a）所示,并不能将基质块完全分割开来^{[3, 17, 27]}。由朗缪尔等温吸附理论可计算出基质的吸附变形为

图 1 基质和裂隙的相关关系图

Figure 1. Interrelationship of matrix and fracture

下载: 全尺寸图片幻灯片

$Δ s_{s} = s Δ ε_{s} = s ε_{L} \frac{p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}$ ,

(7)

式中, $ε_{s}$ 为吸附引起的体积应变, $ε_{L}$ 为朗缪尔体积应变常数,定义为孔隙压力为无穷大时的吸附应变, $p_{L}$ 为朗缪尔压力,定义为测量的体积应变为0.5 $ε_{L}$ 时的孔隙压力, $p_{m0}$ 为初始的基质孔隙压力。

假设基质膨胀是各向同性的,且只有部分吸附应变对裂隙开度造成影响,则由于基质吸附变形造成的裂隙开度变化可由下式计算：

$Δ b_{s} = - f Δ s_{s} = - f s ε_{L} \frac{p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}$ ,

(8)

式中,f为内部膨胀系数,认为是一个变量,随着不同的储层取值不同,取值为0～1。

1.3 动态的渗透率模型

前两部分考虑了基质和裂隙在力学作用下和吸附膨胀作用下产生的变形,其变形过程可总结为图1（b）所示。结合式（5）,（8）可知,裂隙在有效应力和吸附共同作用下产生的变形为

$Δ b = - \frac{b (Δ σ - β Δ p_{f})}{K_{f}} - f s ε_{L} \frac{p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}$ 。

(9)

将式（3）,（9）代入式（2）可得裂隙渗透率的导数形式：

$\begin{matrix} d k_{f} = k_{f} [\frac{(Δ σ - α Δ p_{m})}{K_{m}} - \frac{3 (Δ σ - β Δ p_{f})}{K_{f}} - \\ 3 f \frac{s}{b} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}] \end{matrix}$ 。

(10)

将式（10）两端进行积分并同时除以初始孔隙率 $k_{f0}$ 可得在变应力作用下渗透率的表达式：

$\begin{matrix} \frac{k_{f}}{k_{f0}} = exp [\frac{(σ - σ_{0}) - α (p_{m} - p_{m0})}{K_{m}} - \frac{3 (σ - σ_{0}) - 3 β (p_{f} - p_{f0})}{K_{f}} - \\ 3 f \frac{s_{0}}{b_{0}} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}] \end{matrix}$ 。

(11)

式（11）所构建的渗透率模型与传统的渗透率模型相比,其最主要的区别在于考虑了基质本身在力学作用下的变形,并认为吸附应变对裂隙开度只有部分影响。为了进一步凸显基质变形本身对储层渗透率演化的影响,建立了一个仅考虑裂隙开度变化,进而直接运用立方定律得到的一种简化的渗透率模型,具体表达式为

$\begin{array}{l} \frac{k_{f}}{k_{f0}} = {(1 + \frac{Δ b}{b})}^{3} \\ = {[1 - \frac{(Δ σ - β Δ p_{f})}{K_{f}} - f ε_{L} \frac{s}{b} \frac{p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}]}^{3} \end{array}$ 。

(12)

为方便表述,将式（11）所表示的渗透率模型定义为模型1,式（12）所表示的渗透率模型定义为模型2,至此渗透率模型已经建立完成,下一部分将对模型中的外部应力做进一步的研究。

1.4 外部应力的修正

当不考虑惯性力的作用时,页岩储层变形的平衡方程为^[13]

$σ_{i j, j} + F_{i} = 0 (i = 1, 2, 3; j = 1, 2, 3)$ ,

(13)

式中, $σ_{i j}$ 为储层的总应力分量,F_i为体积力分量。

基于页岩的均质性和各向同性假设,储层变形的几何方程为

$ε_{i j} = \frac{1}{2} (u_{i, j} + u_{j, i})$ ,

(14)

式中, $ε_{i j}$ 为储层总的应变分量,u_i,u_j为位移分量。

Liu等^{[14, 17]}认为基质吸附应变只有部分对裂隙开度有影响,剩下的部分对岩体应变造成影响,因此在研究储层整体应变的本构方程时,在吸附项前应该引入系数（1-f）对其进行修正,根据Shi等^[9]和Liu等^[28]提出的应变本构方程,修正后的表达式为

$\begin{array}{l} ε_{i j} = \frac{1}{E} [(1 + ν) σ_{i j} - ν σ_{k k} δ_{i j}] + \frac{α}{3 K} p_{m} δ_{i j} + \\ \frac{β}{3 K} p_{f} δ_{i j} + \frac{(1 - f) ε_{s}}{3} δ_{i j} \end{array}$ ,

(15)

式中, $σ_{k k}$ 为体积应力, $σ_{k k}$ = $σ_{x x}$ + $σ_{y y}$ + $σ_{z z}$ , $δ_{i j}$ 为Kronecker符号（i=j,δ_ij=1；i≠j,δ_ij=0）。

求解式（15）可得修正后的储层应力表达式

$σ_{i j} = λ ε_{k k} δ_{i j} + 2 G ε_{i j} - α p_{m} δ_{i j} - β p_{f} δ_{i j} - (1 - f) K ε_{s} δ_{i j}$ 。

(16)

式中 $λ$ 为拉梅常数, $λ$ =Eν/(1+ν)(1-2ν)； $ε_{k k}$ 为体积应变, $ε_{k k}$ = $ε_{x x}$ + $ε_{y y}$ + $ε_{z z}$ ；G为剪切模量,G=E/2 (1+ $ν$ )。

将式（14）代入到式（16）中可得各个方向上的有效应力分量,

$\begin{matrix} {σ^{'}}_{x x} = λ ε_{k k} + 2 G ε_{x x} - (1 - f) K ε_{s} ， \\ {σ^{'}}_{y y} = λ ε_{k k} + 2 G ε_{y y} - (1 - f) K ε_{s} ， \\ {σ^{'}}_{z z} = λ ε_{k k} + 2 G ε_{z z} - (1 - f) K ε_{s} ， \end{matrix}}$

(17)

式中, ${σ^{'}}_{x x}$ , ${σ^{'}}_{y y}$ , ${σ^{'}}_{z z}$ 为x,y,z 3个方向的有效应力分量方程。

值得一提的是,若不考虑此修正,计算所得的各个方向的有效应力分量为

$\begin{matrix} {σ^{'}}_{x x} = λ ε_{k k} + 2 G ε_{x x} - K ε_{s} ， \\ {σ^{'}}_{y y} = λ ε_{k k} + 2 G ε_{y y} - K ε_{s} ， \\ {σ^{'}}_{z z} = λ ε_{k k} + 2 G ε_{z z} - K ε_{s} 。 \end{matrix}}$

(18)

与式（17）相比,式（18）没有对应力进行修正,所得吸附应变项前的系数为1,因此可以用于下面部分讨论考虑应力修正与不考虑应力修正的区别。

1.5 不同边界条件下的渗透率模型

通常认为煤岩气藏在原位应力条件下开采时可近似认为单轴应变边界条件^[8]或常体积边界条件^[29],而实验室测量煤岩渗透率时普遍采用恒定围压边界条件^[13]。在本部分中,将对所建立的渗透率模型在不同边界条件下的演变作进一步的研究。

（1）单轴应变条件下渗透率模型

单轴应变条件指水平方向应变为0,且假设在竖直方向承受恒定的应力,即 $ε_{x x}$ = $ε_{y y}$ =0, $ε_{z z}$ ≠0,如图2（a）所示,则式（17）可以重写为

$\begin{array}{l} {σ^{'}}_{x x} = \frac{E ν}{(1 + ν) (1 - 2 ν)} ε_{x x} - \frac{(1 - f) E}{3 (1 - 2 ν)} ε_{s} ， \\ {σ^{'}}_{y y} = \frac{E ν}{(1 + ν) (1 - 2 ν)} ε_{y y} - \frac{(1 - f) E}{3 (1 - 2 ν)} ε_{s} ， \\ {σ^{'}}_{z z} = \frac{E ν}{(1 + ν) (1 - 2 ν)} ε_{z z} + \frac{E}{(1 + ν)} ε_{z z} - \frac{(1 - f) E}{3 (1 - 2 ν)} ε_{s} 。 \end{array}}$

(19)

图 2 不同边界条件原理图

Figure 2. Diagram of different boundary conditions

下载: 全尺寸图片幻灯片

求解式（19）可得水平方向的有效应力,

${σ^{'}}_{x x} = {σ^{'}}_{y y} = \frac{ν}{1 - ν} {σ^{'}}_{z z} - \frac{(1 - f) E}{3 (1 - ν)} ε_{s}$ 。

(20)

式（20）所示的表达式和已有文献[10,30]求得的公式类似,其不同之处在于吸附项前多出一项（1-f）,体现出部分膨胀应变的作用。

前面提到,基质和裂隙对于流体的渗透性能不同,因此流经基质和裂隙中的孔隙压力会存在一定的区别,对于双重孔隙度竖直方向的有效应力可用下式计算^[31-33]：

${σ^{'}}_{z z} = σ_{z z} - (β p_{f} + α p_{m})$ 。

(21)

根据式（20）,（21）可得煤岩储层的平均有效应力为

$\begin{matrix} Δ σ = σ - σ_{0} \\ = \frac{1}{3} (Δ {σ^{'}}_{x x} + Δ {σ^{'}}_{y y} + Δ {σ^{'}}_{z z}^{}) \\ = - \frac{(1 + ν)}{3 (1 - ν)} [β (p_{f} - p_{f0}) + α (p_{m} - p_{m0})] - \\ \frac{2 (1 - f) E}{9 (1 - ν)} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})} \end{matrix}$ 。

(22)

当不考虑应力修正时,只需将式（17）替换为式（18）,重复上述计算步骤,所得到的不考虑应力修正的页岩储层平均有效应力表达式为

$\begin{matrix} Δ σ = σ - σ_{0} \\ = \frac{1}{3} (Δ {σ^{'}}_{x x} + Δ {σ^{'}}_{y y} + Δ {σ^{'}}_{z z}) \\ = - \frac{(1 + ν)}{3 (1 - ν)} [β (p_{f} - p_{f0}) + α (p_{m} - p_{m0})] - \\ \frac{2 E}{9 (1 - ν)} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})} \end{matrix}$ 。

(23)

将式（22）,（23）分别代入式（11）可得单轴应变条件下储层平衡方程考虑应力修正和不考虑应力修正的裂隙渗透率模型

$\begin{array}{l} \frac{k_{f}}{k_{f0}} = \exp {[\frac{2 (2 - ν)}{K_{f} (1 - ν)} - \frac{(1 + ν)}{3 K_{m} (1 - ν)}] β (p_{f} - p_{f 0}) + \\ [\frac{(1 + ν)}{K_{f} (1 - ν)} - \frac{2 (2 - ν)}{3 K_{m} (1 - ν)}] α (p_{m} - p_{m0}) + \\ [\frac{2 E (1 - f) (3 K_{m} - K_{f})}{9 K_{m} K_{f} (1 - ν)} - 3 f \frac{s_{0}}{b_{0}}] \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}} \end{array}$ ,

(24)

$\begin{matrix} \frac{k_{f}}{k_{f0}} = \exp {[\frac{2 (2 - v)}{K_{f} (1 - v)} - \frac{(1 + v)}{3 K_{m} (1 - v)}] β (p_{f} - p_{f0}) + \\ [\frac{(1 + v)}{K_{f} (1 - v)} - \frac{2 (2 - v)}{3 K_{m} (1 - v)}] α (p_{m} - p_{m0}) + \\ [\frac{2 E (3 K_{m} - K_{f})}{9 K_{m} K_{f} (1 - ν)} - 3 f \frac{s_{0}}{b_{0}}] \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}} \end{matrix}$ 。

(25)

进一步地,为了分析基质本身变形对渗透率演化的影响,将式（23）代入式（12）可得单轴应变条件下简化的渗透率模型

$\begin{matrix} \frac{k_{f}}{k_{f0}} = {1 + \frac{2 （ 2 - ν)}{3 K_{f} (1 - ν)} β (p_{f} - p_{f0}) - \frac{(1 + ν)}{3 K_{f} (1 - ν)} α (p_{m} - p_{m0}) + \\ {[\frac{2 (1 - f) E}{9 K_{f} (1 - ν)} - f \frac{s}{b}] \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}}}^{3} \end{matrix}$ 。

(26)

式（24）在渗透率模型建立过程中同时考虑基质本身变形和裂隙开度变形,并对应力项进行了修正,为方便表述,定义为模型1-11；式（25）同样考虑基质变形和裂隙开度变形,但是没有对应力项进行修正,定义为模型1-12；式（26）为在模型建立过程中考虑了应力项的修正,但仅考虑裂隙开度变化,之后直接运用立方定律得到的一种简化的渗透率模型,定义为模型2-1；这3个模型均是在单轴应变条件下推导而来,将用于下一部分的验证。

（2）常体积条件下的渗透率模型

煤岩气藏在原位应力开采的另一种近似条件为常体积条件,如图2（b）所示。这种条件指的是四周边界固定,各个方向的应变为0,即 $ε_{x x}$ = $ε_{y y}$ = $ε_{z z}$ =0。在常体积条件下式（17）所示的各个方向的有效应力分量可重写为

$\begin{matrix} {σ^{'}}_{x x} = - \frac{(1 - f) E}{3 (1 - 2 ν)} ε_{s} ， \\ {σ^{'}}_{y y} = - \frac{(1 - f) E}{3 (1 - 2 ν)} ε_{s} ， \\ {σ^{'}}_{z z} = - \frac{(1 - f) E}{3 (1 - 2 ν)} ε_{s} 。 \end{matrix}}$

(27)

由式（27）可以计算常体积条件下的平均有效应力为

$\begin{matrix} Δ σ = σ - σ_{0} \\ = \frac{1}{3} (Δ {σ^{'}}_{x x} + Δ {σ^{'}}_{y y} + Δ {σ^{'}}_{z z}) \\ = - \frac{(1 - f) E}{3 (1 - 2 ν)} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})} \end{matrix}$ 。

(28)

同理,当不考虑应力修正项时,重写式（18）并按上述计算步骤可得常体积条件下的平均有效应力为

$\begin{matrix} Δ σ = σ - σ_{0} \\ = \frac{1}{3} (Δ {σ^{'}}_{x x} + Δ {σ^{'}}_{y y} + Δ {σ^{'}}_{z z}) \\ = - \frac{E}{3 (1 - 2 ν)} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})} \end{matrix}$ 。

(29)

将式（28）,（29）分别代入到式（11）可得常体积条件下考虑应力修正和不考虑应力修正的裂隙渗透率模型：

$\begin{matrix} \frac{k_{f}}{k_{f0}} = \exp {- \frac{1}{K_{m}} - α (p_{m} - p_{m0}) + \frac{3}{K_{f}} β (p_{f} - p_{f0}) + \\ [\frac{(3 K_{m} - K_{f}) (1 - f) E}{3 K_{m} K_{f} (1 - 2 ν)} - 3 f \frac{s_{0}}{b_{0}}] \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}} \end{matrix}$ ,

(30)

$\begin{matrix} \frac{k_{f}}{k_{f0}} = \exp {- \frac{1}{K_{m}} α (p_{m} - p_{m0}) + \frac{3}{K_{f}} β (p_{f} - p_{f0}) + \\ [\frac{(3 K_{m} - K_{f}) E}{3 K_{m} K_{f} (1 - 2 ν)} - 3 f \frac{s_{0}}{b_{0}}] \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}} \end{matrix}$ 。

(31)

将式（28）代入式（12）可得常体积条件下仅考虑裂隙开度变化的简化的渗透率模型：

$\begin{array}{l} \frac{k_{f}}{k_{f0}} = {1 + \frac{1}{K_{f}} β (p_{f} - p_{f0}) + \\ {[\frac{(1 - f) E}{3 K_{f} (1 - 2 ν)} - f \frac{s}{b}] \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}}}^{3} \end{array}$ 。

(32)

与单轴应变条件下渗透率模型相似,式（30）～（32）均是在常体积条件下建立的渗透率模型,其中式（30）建立的渗透率模型同时考虑了基质本身变形和裂隙开度变形,并对应力项进行了修正,定义为模型1-21,与式（30）相比,式（31）建立的渗透率模型没有考虑应力项的修正,定义为模型1-22；式（32）是在对应力项进行修正的基础上,只考虑裂隙开度变化,并直接运用立方定律得到的一种简化的渗透率模型,定义为模型2-2。

（3）恒定围压条件下渗透率模型

在实验室测量岩芯渗透率时,常采用如图2（c）所示的恒定围压边界条件,因此建立恒定围压条件下的渗透率模型对于实验室岩芯测量至关重要。在这种条件下,平均应力的增量为0,即

$σ - σ_{0} = 0$ 。

(33)

此时不存在应力项修正的问题,将式（33）代入式（11）可得页岩储层渗透率模型为

$\begin{matrix} \frac{k_{f}}{k_{f0}} = exp [\frac{- α (p_{m} - p_{m0})}{K_{m}} - \frac{- 3 β (p_{f} - p_{f0})}{K_{f}} - \\ 3 f \frac{s_{0}}{b_{0}} \frac{ε_{L} p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}] \end{matrix}$ 。

(34)

同理,将式（33）代入式（12）可得简化的储层渗透率模型为

$\frac{k_{f}}{k_{f0}} = {[1+ \frac{β (p_{f} - p_{f0})}{K_{f}} - f ε_{L} \frac{s}{b} \frac{p_{L} (p_{m} - p_{m0})}{(p_{m} + p_{L}) (p_{m0} + p_{L})}]}^{3}$ 。

(35)

式（34）,（35）为在恒定围压条件下建立的渗透率模型,其中式（34）建立的模型同时考虑基质本身变形和裂隙开度对储层渗透率演化的影响,定义为模型1-3；式（35）为仅考虑裂隙开度变化所建立的一种简化的渗透率模型,定义为模型2-3。不同边界条件下的渗透率模型已经建立完成,图3为具体建立模型的流程,下一部分将对已经建立的渗透率模型进行验证。

图 3 渗透率建模的工作流程

Figure 3. Workflow used in permeability modeling

下载: 全尺寸图片幻灯片

2. 模型的验证

在本部分,将在不同边界条件下对前面部分建立的渗透率模型进行验证,主要为用现场测井数据验证单轴应变条件和常体积条件下的渗透率模型,用实验室数据验证恒定围压条件下的渗透率模型。现有的室内测试渗透率实验一般是在吸附平衡的条件下进行的^[34],而现场条件下某一测井的数据也是在储层局部平衡的基础上测试获取。在这一平衡状态,假设基质和裂隙的孔隙压力相等,由于此时基质与裂缝之间的质量传递是处于平衡状态,故不需要考虑气体在基质中的扩散作用。

2.1 单轴应变条件下渗透率模型的验证

单轴应变条件下渗透率的模型为模型1-11,1-12和2-1,分别对应于式（24）～（26）。用于验证的现场数据为圣胡安内Fruitland煤层的6个测井^[8,21]。表1给出了这6个测井的初始渗透率以及不同测点处的储层压力和对应的渗透率。采用matlab中的cftool进行拟合,其原理为非线性最小二乘法中的信赖域算法。模型拟合所用的各个参数如表2所示,拟合结果如图4所示。综合6个测井的拟合结果可知,与简化的模型2-1相比,模型1-11和模型1-12的变化趋势基本相同,除测井2的数据外,模型1-12和模型1-11在低压处的渗透率差异高于在高压处的差异,表明修正应力项对渗透率演化的影响在孔隙压力较低时更为明显。在图4中,当孔隙压力较高时,测井2,3,5中模型2-1的拟合曲线均呈现出明显的脱离现场数据的趋势。此外,在6个测井的拟合曲线中,孔隙压力较低时,模型2-1的渗透率在3个模型中最低。正如前文中所述,由于模型2-1仅考虑了裂隙开度的变化,没有考虑基质本身对于渗透率演变的影响,因此模型2-1在预测煤层气开采时渗透率的演变效果较模型1-11和1-12更差。

表 1 圣胡安盆地Fruitland 煤层测量的现场数据^[8]

Table 1. Field data measured in Fruitland coal seam in San Juan Basin

测井	k ₀/mD	储层压力(MPa)/渗透率(mD)
测井	k ₀/mD	测点1	测点2	测点3
井1	1.8	4.85/7.2	3.12/16.5	2.43/20.6
井2	6.2	4.61/1.7	3.52/2.7	3.04/2.9
井3	5.3	4.52/4	3.03/18.5	2.07/28.3
井4	2.1	3.38/9.3	2.88/9.6	2.21/16.2
井5	9.1	3.2/10.6	2.48/14.5	2.05/23.5
井6	2.1	2.98/11.5	2.11/24.3	1.96/28.8

下载: 导出CSV

| 显示表格

表 2 拟合Fruitland煤层渗透率数据所用的参数

Table 2. Parameters used for matching permeability data of Fruitland coal seam

符号	值	含义
E	2902^[9]	煤的杨氏模量(MPa)
E_m	8143^[21]	基质的杨氏模量(MPa)
ν	0.35^[9]	煤的泊松比
p_L	4.3 ^[9]	朗缪尔压力(MPa)
ε_L	0.01266^[9]	朗缪尔体积应变
b₀	1×10^-6	裂隙开度(m)
s₀	0.01^[21]	基质宽度(m)

下载: 导出CSV

| 显示表格

图 4 单轴应变条件下现场数据与模型的匹配图

Figure 4. Matching curves of field data with model under uniaxial strain

下载: 全尺寸图片幻灯片

2.2 常体积条件下渗透率模型的验证

通常认为常体积和单轴应变均能很好地近似模拟储层开采时的边界条件。因此,同单轴应变条件下模型验证一样,在本部分仍采用圣胡安盆地内Fruitland煤层的6个测井数据对常体积条件下的渗透率模型进行验证。图5为3个模型与现场数据的拟合结果,与图4相似,模型1-21与模型1-22的变化趋势基本相同,且在低压处渗透率的差异高于高压处渗透率的差异,进一步表明修正应力项在孔隙压力较低时更加明显。然而与单轴应变条件下拟合结果相比,模型1-21和模型1-22在低压处与高压处的渗透率差异更小。此外,模型2-2在孔隙压力较高时渗透率的下降速度更快,在孔隙压力较低时渗透率明显低于其他两个模型,对现场数据的拟合效果也明显更差。

图 5 常体积条件下现场数据与模型的匹配图

Figure 5. Matching curves of field data with model under constant volume

下载: 全尺寸图片幻灯片

图4,5拟合所得的内部膨胀系数f如表3所示。在单轴应变条件下模型1-12拟合所得的f值平均比模型1-11拟合的值大1.46%,然而这两个模型所得的f值之差明显小于模型2-1与模型1-11两者之间的平均差值7.87%。常体积条件下的3个模型对6个测井数据拟合所得的f值与单轴应变条件的总体规律相同,即模型2-2所得f值平均比模型1-21大16.3%,远远高于模型1-22与模型1-21所得的f的平均差值0.16%。

表 3 单轴应变和常体积条件下不同测井的f值

Table 3. Values of f of different wells under uniaxial strain and constant volume

测井	单轴应变条件			常体积条件
测井	模型 1-11	模型 1-12	模型 2-1	模型 1-21	模型 1-22	模型 2-2
井1	0.1020	0.1040	0.1995	0.0800	0.0800	0.1545
井2	0.3525	0.3795	0.4310	0.0720	0.0745	0.0930
井3	0.1875	0.1940	0.2475	0.1025	0.1040	0.1520
井4	0.1425	0.1460	0.2250	0.1035	0.1045	0.1630
井5	0.2575	0.2710	0.2900	0.1300	0.1320	0.1450
井6	0.1995	0.2065	0.3070	0.1455	0.1470	0.2725

下载: 导出CSV

| 显示表格

2.3 恒定围压条件下渗透率模型的验证

基于恒定围压建立的模型为1-3和模型2-3,采用Robertson^[13]在实验室内用甲烷气体测得的Anderson 01和Gilson 02岩样的相关数据进行验证。同样的,假设试验状态下的基质和裂隙内孔隙压力相等。Anderson 01和Gilson 02岩样的拟合结果分别如图6,7所示,拟合所用的参数值如表4所列。从图6,7可知,模型1-3对实验数据的拟合很好,然而模型2-3却无法准确地拟合实验室内的测量数据。拟合结果表明模型1-3能够很好的预测在恒定围压下的渗透率演变。

图 6 用甲烷所测Anderson 01 渗透率数据与模型的拟合关系

Figure 6. Fitting relationship between model and experiment permeability data for core Anderson 01 using pure CH4 gas

下载: 全尺寸图片幻灯片

图 7 用甲烷所测Gilson 02 渗透率数据和模型的拟合关系

Figure 7. Fitting relationship between model and experiment permeability data for core Gilson 02 using pure CH4 gas

下载: 全尺寸图片幻灯片

表 4 模型中用到的拟合参数^{[3, 13, 29]}

Table 4. Fitting parameters used in models

煤样	煤杨氏模量E/MPa	朗缪尔压力 $p_{L}$ /MPa	朗缪尔体积应变 $ε_{L}$	煤的泊松比 $ν$	裂隙开度/(10^-5m)	基质宽度/m
Anderson 01	1379	6.11	0.00931	0.35	1	0.01
Gilson 02	1379	6.11	0.00765	0.35	1	0.01
Sulcis 煤	1120	7.25	0.04900	0.26	1	0.01

下载: 导出CSV

| 显示表格

3. 结果和讨论

本节对前面模型验证的结果做了比较分析,并进一步用现场数据和实验室数据对模型1和经典的模型进行了比较研究。通过不同朗缪尔体积应变常数,不同基质宽度和不同裂隙开度对模型1做了敏感性分析,从而研究不同气体类型,不同裂隙发育程度以及不同裂隙张开度对渗透率演变的影响。

3.1 模型1和模型2的比较

在第三部分用现场测井数据和实验室数据分别对单轴应变、常体积以及恒定围压条件下的模型进行了验证。验证结果表明,当引入内部膨胀系数f,并对应力项进行修正时,修正效果在单轴应变条件下更为明显,如图4,5所示。对于常体积近似储层边界条件,由于模型1-21和模型1-22较小的拟合曲线差异（图5）,以及较小的f值差异（表3）,认为系数f对应力项的影响可以忽略不计。且通过分析图4,5可知,模型2-1与模型2-2相对于模型1-11与模型1-21,虽也能拟合测井数据,但是在孔隙压力较低时拟合所得的渗透率偏低,且在孔隙压力较高时,有偏离实际情况的趋势。结合图6,7可知,模型2-3无法准确的拟合实验室测得的渗透率数据,进一步表明仅考虑裂隙开度变形建立的渗透率模型其适用范围和准确程度均小于考虑基质本身变形所建立的模型。因此分析模型1对于深入研究多孔介质渗透率的演变规律更具有价值。

3.2 模型1和经典模型的比较

一般在拟合现场尺度煤层条件时通常采用单轴应变和常体积条件,拟合实验室测量条件采用恒定围压条件。采用Shi等^[10]使用的典型的圣胡安盆地的现场数据和Pini等^[35]测量的实验室数据对不同边界条件下的新模型与传统不考虑内膨胀系数模型的C-B模型,S-D模型,P-M模型做了比较。为了进一步说明新模型的优越性,还选取了引用量最高的考虑内膨胀系数的模型进行对比,定义为L-R模型^[14]。所采用圣胡安盆地的基本参数和拟合参数见表2,其余模型的物理参数不变。Pini等^[35]测量的参数见表4的Sulcis煤,初始压力值为3.76 MPa。通过调整S-D模型,C-B模型的裂隙压缩性、P-M模型的孔隙度以及L-R模型的裂隙应变变化量和裂隙压缩性来获取最佳的拟合效果。拟合所得的参数见表5,比较结果如图8,9所示。

表 5 4种模型在不同边界条件下的拟合参数

Table 5. Fitting parameters of four models under different boundary conditions

测试环境	S-D模型	C-B模型	P-M模型	L-R模型
现场	C_f=0.0967 MPa^-1	C_f=0.2339 MPa^-1	ϕ_f=9.16×10^-4	C_f=0.02031 MPa^-1
现场	C_f=0.0967 MPa^-1	C_f=0.2339 MPa^-1	ϕ_f=9.16×10^-4	Δε_f=0.0039
室内	C_f=0.01 MPa^-1	C_f=0.001 MPa^-1	ϕ_f=0.01	C_f=0.2926 MPa^-1
室内	C_f=0.01 MPa^-1	C_f=0.001 MPa^-1	ϕ_f=0.01	Δε_f=1.69×10^-5

下载: 导出CSV

| 显示表格

图 8 现场条件下4种经典模型与新模型之间的比较

Figure 8. Comparison of four common models with new models under field conditions

下载: 全尺寸图片幻灯片

图 9 实验室条件下4种经典模型与新模型之间的比较

Figure 9. Comparison of four common models with new models under field conditions

下载: 全尺寸图片幻灯片

图8为几种模型对现场数据的拟合结果,由图可知,考虑内膨胀系数的模型1-11,1-12和L-R模型对数据的拟合程度明显优于其他模型。此外,3个模型对现场数据拟合方程的R²分别达到了0.9616,0.9594,0.9567。造成拟合程度差异的原因如下：3个模型均考虑内膨胀系数的影响,但是模型1-11在此基础上同时考虑基质在有效应力下的变形和吸附膨胀对总应力的修正,而模型1-12却忽略了吸附对应力项造成的调整,至于L-R模型则同时忽略基质在应力下的变形和吸附带来的应力项修正。图9为恒定围压条件下4种模型对实验数据的拟合比较,很明显模型1-3对数据的拟合效果最好,R²达到了0.9972。L-R模型对数据的拟合程度次之,R²达到了0.9885。C-B模型和S-D模型均只能拟合低压范围内的数据,P-M在整个压力范围内的拟合结果均较差。当结合图8,9可知,考虑内膨胀系数的模型同时适用于现场条件和实验室条件下的数据预测,而传统的C-B模型,S-D模型和P-M模型虽也能初粗略的适用现场条件,但是无法适用于实验室条件下渗透率数据的预测。

4. 结论

建立了一个新的描述储层渗透率演化的模型,该模型考虑了基质和裂隙在力学作用和吸附膨胀作用下的变形,并对储层应力的表达式进行了修正。随后对渗透率模型进行简化,并采用现场测井数据和实验室数据在不同边界条件下对这几个模型进行比较验证,之后与传统的考虑内膨胀系数和不考虑内膨胀系数的渗透率模型进行了比较。基于模型验证和分析结果,可以得到3点结论。

（1）新的模型可以同时适用现场条件和实验室条件测得的数据预测,对数据的拟合程度均很好。当简化为不考虑基质变形的模型2时,在现场条件下的拟合效果变差,且无法拟合恒定围压条件下测得的数据。此外,与新模型相比,模型2在单轴应变条件和常体积条件下所得的内部膨胀系数较大。

（2）对内部膨胀系数量化分析后发现,当用单轴应变近似储层条件时,应力修正对渗透率演化的影响比常体积条件更为明显,且拟合所得的f值通常也更大。此外,应力项的修正对渗透率演化的影响在孔隙压力较低时更为明显。

（3）考虑内部膨胀系数的模型对现场和实验室数据的拟合程度均优于不考虑内膨胀系数的模型,且考虑基质在有效应力下的变形和吸附膨胀对应力的修正均能不同程度的提高对数据的匹配效果。

图 1 应力应变曲线标定结果

Figure 1. Calibrated results of stress-strain response

下载: 全尺寸图片幻灯片

图 2 地基模型示意图

Figure 2. Schematic view of soil foundation

下载: 全尺寸图片幻灯片

图 3 初始平衡状态下地基参数沿深度分布特征

Figure 3. Distribution of foundation parameters at balance

下载: 全尺寸图片幻灯片

图 4 不同贯入速度下P2沉桩总阻力发展特征

Figure 4. Development of installation resistance for P2 under different penetration speeds

下载: 全尺寸图片幻灯片

图 5 沉桩结束时桩周土体变形特征

Figure 5. Soil deformations around piles at the end of pile installation

下载: 全尺寸图片幻灯片

图 6 沉桩过程中桩端土体运移轨迹特征

Figure 6. Particle movements relative to pile tip during pile Installation

下载: 全尺寸图片幻灯片

图 7 土塞参数随贯入深度变化特征

Figure 7. Development of soil plugging parameter with penetration.depth

下载: 全尺寸图片幻灯片

图 8 土塞增长率与土塞率经验拟合关系

Figure 8. Relationship between IFR and PLR

下载: 全尺寸图片幻灯片

图 9 沉桩阻力随贯入深度变化特征

Figure 9. Development of pile installation resistance with penetration depth

下载: 全尺寸图片幻灯片

图 10 单位内摩阻力随贯入深度变化特征

Figure 10. Development of unit inner shaft resistance with penetration depth

下载: 全尺寸图片幻灯片

图 11 内摩阻力沿桩身分布特征

Figure 11. Distribution of inner shaft resistance

下载: 全尺寸图片幻灯片

图 12 承载结束时桩周土体变形特征

Figure 12. Soil deformations around pile at the end of static loading

下载: 全尺寸图片幻灯片

图 13 桩顶位移对贯入阻力的影响

Figure 13. Development of penetration resistance with displacement of pile head

下载: 全尺寸图片幻灯片

图 14 q_plug/q_c和q_ann/q_c与IFR经验拟合关系

Figure 14. Relationship between q_plug/q_c, q_ann/q_c and IFR

下载: 全尺寸图片幻灯片

图 15 归一化桩端阻力(q_b/q_c)经验拟合关系

Figure 15. Normalized relationship between q_b/q_c and A_{rb, eff}

下载: 全尺寸图片幻灯片

表 1 UWA硅质砂及其等效离散砂样参数

Table 1 Properties of UWA superfine sand and its DEM analogue

参数	G_s	颗粒尺寸/mm			C_u	C_c	${\varphi '_{{\text{cs}}}}$ /(°)
参数	G_s	d₁₀	d₅₀	d₆₀	C_u	C_c	${\varphi '_{{\text{cs}}}}$ /(°)
试验值	2.67	0.12	0.18	0.20	1.67	1.020	33
模拟值	2.67	1.80	2.70	3.00	1.67	0.975	32
注： ${\varphi '_{{\text{cs}}}}$ 为极限状态摩擦角，C_u，C_c分别为不均匀系数和曲率系数。

下载: 导出CSV

表 2 离散元模型参数

Table 2 Parameters for discrete model

地基土样								墙		桩
k_n/(N·m^-1)	k_s/(N·m^-1)	μ	$\gamma$	μ_r	β_n	β_s	M_d	k_n/(N·m^-1)	k_s/(N·m^-1)	μ_w	k_n/(N·m^-1)	k_s/(N·m^-1)	μ_p
1×10⁸	2×10⁷	0.5	0.1	0.45	0.1	0.1	3	1×10¹⁰	1×10¹⁰	0	1.5×10⁹	1.5×10⁹	0.5
注：k_n，k_s分别为法向刚度和剪切刚度；μ为摩擦系数； $\gamma$ 为阻尼比；μ_r为滚动摩擦系数；β_n为法向极限阻尼比；β_s为剪切极限阻尼比；M_d为黏壶模式；μ_w为墙摩擦系数；μ_p为桩摩擦系数。

下载: 导出CSV

参考文献(32)

[1]	陆昭球, 高倚山, 宋铭栋. 关于开口钢管桩工作性状的几点认识[J]. 岩土工程学报, 1999, 21(1): 111-114. doi: 10.3321/j.issn:1000-4548.1999.01.024 LU Zhaoqiu, GAO Yishan, SONG Mingdong. Some knowledge about working performance open-end steel pipe pile[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 111-114. (in Chinese) doi: 10.3321/j.issn:1000-4548.1999.01.024
[2]	李立辰. 考虑土塞效应的管桩与复合单桩基础承载特性试验及数值研究[D]. 武汉: 中国地质大学(武汉), 2021. LI Lichen. Experimental and Numerical Research on the Bearing Characteristics of Open-ended Pile and Hybrid Pile Foundation Considering Plugging Effect[D]. Wuhan: China University of Geosciences (Wuhan), 2021. (in Chinese)
[3]	张忠苗, 刘俊伟, 俞峰, 等. 静压预应力混凝土管桩土塞效应试验研究[J]. 岩土力学, 2011, 32(8): 2274-2280. doi: 10.3969/j.issn.1000-7598.2011.08.005 ZHANG Zhongmiao, LIU Junwei, YU Feng, et al. Research on plugging effect of jacked prestressed concrete pipe pile[J]. Rock and Soil Mechanics, 2011, 32(8): 2274-2280. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.08.005
[4]	吴文兵, 蒋国盛, 王奎华, 等. 土塞效应对管桩纵向动力特性的影响研究[J]. 岩土工程学报, 2014, 36(6): 1129-1141. doi: 10.11779/CJGE201406019 WU Wenbing, JIANG Guosheng, WANG Kuihua, et al. Influence of soil plug effect on vertical dynamic response of pipe piles[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(6): 1129-1141. (in Chinese) doi: 10.11779/CJGE201406019
[5]	RANDOLPH M F, LEONG E C, HOULSBY G T. One-dimensional analysis of soil plugs in pipe piles[J]. Géotechnique, 1991, 41(4): 587-598. doi: 10.1680/geot.1991.41.4.587
[6]	郑俊杰, 聂重军, 鲁燕儿. 基于土塞效应的柱形孔扩张问题解析解[J]. 岩石力学与工程学报, 2006, 25(增刊2): 4004-4008. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2006S2104.htm ZHENG Junjie, NIE Chongjun, LU Yaner. Analytical solutions of cylindrical cavity expansion problems considering plugging effects[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(S2): 4004-4008. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2006S2104.htm
[7]	刘润, 郭绍曾, 周龙, 等. 大直径钢管桩土塞效应的拟静力判断方法[J]. 地震工程学报, 2017, 39(1): 20-27. doi: 10.3969/j.issn.1000-0844.2017.01.0020 LIU Run, GUO Shaozeng, ZHOU Long, et al. Quasi-static method to evaluate the soil plug effect on a large-diameter steel pipe pile[J]. China Earthquake Engineering Journal, 2017, 39(1): 20-27. (in Chinese) doi: 10.3969/j.issn.1000-0844.2017.01.0020
[8]	LI L C, LIU X, LIU H, et al. Experimental and numerical study on the static lateral performance of monopile and hybrid pile foundation[J]. Ocean Engineering, 2022, 255: 111461. doi: 10.1016/j.oceaneng.2022.111461
[9]	刘俊伟. 静压开口混凝土管桩施工效应试验及理论研究[D]. 杭州: 浙江大学, 2012. LIU Junwei. Experimental and Theoretical Study on Construction Effect of Static Pressure Open Concrete Pipe Pile[D]. Hangzhou: Zhejiang University, 2012. (in Chinese)
[10]	俞峰, 张忠苗. 混凝土开口管桩竖向承载力的经验参数法设计模型[J]. 土木工程学报, 2011, 44(7): 100-110. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201107016.htm YU Feng, ZHANG Zhongmiao. A design framework for evaluating the vertical bearing capacity of open-ended concrete pipe pile from empirical correlations[J]. China Civil Engineering Journal, 2011, 44(7): 100-110. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201107016.htm
[11]	HAN F, GANJU E S, PREZZI M, et al. Axial resistance of open-ended pipe pile driven in gravelly sand[J]. Géotechnique, 2020, 70(2): 138-152. doi: 10.1680/jgeot.18.P.117
[12]	曹兆虎, 孔纲强, 刘汉龙, 等. 基于透明土的管桩贯入特性模型试验研究[J]. 岩土工程学报, 2014, 36(8): 1564-1568. doi: 10.11779/CJGE201408025 CAO Zhaohu, KONG Gangqiang, LIU Hanlong, et al. Model tests on pipe pile penetration by using transparent soils[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(8): 1564-1568. (in Chinese) doi: 10.11779/CJGE201408025
[13]	KIKUCHI Y, SATO T, MIZUTANI T, et al. Plugging mechanism of open‐ended piles[J]. Advances in Computed Tomography for Geomaterials, 2010: 406-413.
[14]	周健, 陈小亮, 周凯敏, 等. 静压开口管桩沉桩过程模型试验及数值模拟[J]. 岩石力学与工程学报, 2010, 29(增刊2): 3839-3846. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S2057.htm ZHOU Jian, CHEN Xiaoliang, ZHOU Kaimin, et al. Model test and numerical simulation of driving process of open-ended jacked pipe piles[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 3839-3846. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S2057.htm
[15]	LI L C, WU W B, LIU H, et al. DEM analysis of the plugging effect of open-ended pile during the installation process[J]. Ocean Engineering, 2021, 220: 108375. doi: 10.1016/j.oceaneng.2020.108375
[16]	LI L C, ZHENG M Y, LIU X, et al. Numerical analysis of the cyclic loading behavior of monopile and hybrid pile foundation[J]. Computers and Geotechnics, 2022, 144: 104635. doi: 10.1016/j.compgeo.2022.104635
[17]	Itasca Consulting Group, Inc. PFC2D User's Manual, Version 5.0. 64[R]. Miuneapoils: Itasca Consulting, Group, 2014.
[18]	刘清秉, 项伟, BMLEHANE, 等. 颗粒形状对砂土抗剪强度及桩端阻力影响机制试验研究[J]. 岩石力学与工程学报, 2011, 30(2): 400-410. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201102026.htm LIU Qingbing, XIANG Wei, LEHANE B M, et al. Experimental study of effect of particle shapes on shear strength of sand and tip resistance of driven piles[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(2): 400-410. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201102026.htm
[19]	蒋明镜, 肖俞, 陈双林, 等. 砂土中单桩竖向抗压承载机制的离散元分析[J]. 岩土力学, 2010, 31(增刊2): 366-372. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2010S2066.htm JIANG Mingjing, XIAO Yu, CHEN Shuanglin, et al. Discrete element analysis of bearing mechanism of single pile in sand under vertical load[J]. Rock and Soil Mechanics, 2010, 31(S2): 366-372. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2010S2066.htm
[20]	ZHANG Z T, WANG Y H. Three-dimensional DEM simulations of monotonic jacking in sand[J]. Granul Matter. 2015, 17: 359-376. doi: 10.1007/s10035-015-0562-4
[21]	CIANTIA M, BOSCHI K, SHIRE T, et al. Numerical techniques for fast generation of large discrete-element models[J]. Engineering and Computational Mechanics, 2018, 171(4): 147-161.
[22]	郦建俊, 黄茂松, 王卫东, 等. 开挖条件下抗拔桩承载力的离心模型试验[J]. 岩土工程学报, 2010, 32(3): 388-396. http://www.cgejournal.com/cn/article/id/12408 LI Jianjun, HUANG Maosong, WANG Weidong, et al. Centrifugal model tests on bearing capacity of uplift piles under deep excavation[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(3): 388-396. (in Chinese) http://www.cgejournal.com/cn/article/id/12408
[23]	杨济铭, 张小勇, 张福友, 等. 砂土中桩-土-承台协同作用下桩基承载特性细观研究[J]. 岩土力学, 2020, 41(7): 2271-2282. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202007014.htm YANG Jiming, ZHANG Xiaoyong, ZHANG Fuyou, et al. Mesoscopic study on bearing characteristics of pile foundation under pile-soil-cap combined interaction in sand[J]. Rock and Soil Mechanics, 2020, 41(7): 2271-2282. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202007014.htm
[24]	LEHANE B M, SCHNEIDER J A, XU X. The UWA-05 method for prediction of axial capacity of driven piles in sand[C]// Frontiers in Offshore Geotechnics: Proceedings of the International Symposium on Frontiers in Offshore Geotechnics. London, 2015.
[25]	XU X T, SCHNEIDER J A, LEHANE B M. Cone penetration test (CPT) methods for end-bearing assessment of open- and closed-ended driven piles in siliceous sand[J]. Canadian Geotechnical Journal, 2008, 45(8): 1130-1141. doi: 10.1139/T08-035
[26]	LIU J W, DUAN N, CUI L, et al. DEM investigation of installation responses of jacked open-ended piles[J]. Acta Geotechnica, 2019, 14(6): 1805-1819. doi: 10.1007/s11440-019-00817-7
[27]	PAIK K, SALGADO R. Determination of bearing capacity of open-ended piles in sand[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(1): 46-57. doi: 10.1061/(ASCE)1090-0241(2003)129:1(46)
[28]	FELLENIUS B H. The response of a "plug" in an open-toe pile[J]. Geotechnical Engineering Journal of the SEAGS & AGSSEA, 2015, 46(2): 82-86.
[29]	童立元, 李洪江, 刘松玉, 等. 基于静力触探试验的基坑开挖卸荷单桩水平承载力损失预测研究[J]. 岩土工程学报, 2019, 41(3): 501-508. doi: 10.11779/CJGE201903012 TONG Liyuan, LI Hongjiang, LIU Songyu, et al. Prediction of lateral capacity losses of a single pile adjacent to excavation of foundation pits based on CPT tests[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(3): 501-508. (in Chinese) doi: 10.11779/CJGE201903012
[30]	俞峰, 杨峻. 砂土中钢管桩承载力的静力触探设计方法[J]. 岩土工程学报, 2011, 33(增刊2): 349-354. http://www.cgejournal.com/cn/article/id/14379 YU Feng, YANG Jun. Design methods for bearing capacity of steel pipe piles driven in sand by means of cone penetration tests[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(S2): 349-354. (in Chinese) http://www.cgejournal.com/cn/article/id/14379
[31]	YU F, YANG J. Base capacity of open-ended steel pipe piles in sand[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2011, 138(9): 1116-1128.
[32]	LEHANE B M, GAVIN K G. Base resistance of jacked pipe piles in sand[J]. Journal Geotechnical Geoenvironmental Engineering, 2001, 127(6): 473-480.

施引文献(14)

期刊类型引用(4)

1.	冉志红，董国华，杨子富，林帆. 基于经验频率分布的偏态概率模型可靠度计算方法研究. 四川建筑科学研究. 2024(01): 52-59 . 百度学术
2.	易顺，潘家军，王艳丽，徐晗，白强强，杨志勇. 基坑地表沉降的偏态分布函数应用. 长江科学院院报. 2024(08): 135-141+179 . 百度学术
3.	周勇，王延凯，王正振，王宁. 考虑局部超载下的桩锚支护结构变形分析. 科学技术与工程. 2021(14): 5943-5950 . 百度学术
4.	赵平，韩治勇. 深基坑开挖引起的地表竖向变形分析. 安徽科技学院学报. 2021(06): 112-118 . 百度学术