Dynamic model for water-rock interface of softening of soft rock and its evolution law
-
摘要: 软岩的软化机制本质上是水岩界面改变引起的微细观结构的劣化,对水岩界面的研究具有十分重要的意义。考虑到软岩的一般性特征,从胶结成岩过程的逆过程出发,对不同浸水时间的粉砂质泥岩浸泡液进行了化学元素测试。利用偏光显微薄片研究了软岩微细观结构的演化过程,提出了界面胶结连接结构的概念,建立了软岩软化过程的水岩界面演化模型。基于扩散理论,推导了描述水岩界面演化的理论方程,并与试验结果进行比对。推导和提出了软岩细观软化损伤因子,并将其引入块体离散元数值计算中。研究表明:软岩软化过程中伴随着颗粒的脱落和悬浮以及可溶性物质的溶解,软岩的软化具有明显的非线性动力学特征,理论计算值与试验结果的拟合程度较高,验证了理论方程的可靠性和合理性。数值模拟结果与试验结果吻合度较高。Abstract: The softening mechanism of soft rock is essentially the deterioration of microstructure caused by change of water-rock interface. The research on the water-rock interface is of great importance. Considering the general characteristics of the soft rock, from the idea of the inverse process of cemented diagenesis, the chemical element analysis of silty mudstone soak solution with different immersion time is carried out. The evolution process of microstructure of the soft rock is studied by using the polarized light micro-section. The concept of interfacial cemented bonding structure is proposed to investigate the established model for evolution of water-rock interface. Based on the diffusion theory, the theoretical equation describing the evolution of water-rock interface is derived and compared with the experimental results. The meso-softening damage factor of soft rock is proposed and deduced, and introduced into block discrete element simulation. The results show that the softening process of the soft rock is accompanied by the shedding and suspension of particles and the dissolution of soluble substances. Besides, the softening of the soft rock has obvious nonlinear dynamic characteristics. The fitting degree between the calculated values and the experimental results is relatively high, which verifies the reliability and rationality of the theoretical equation. The numerical results are in good agreement with the experimental ones.
-
0. 引言
页岩油是一种在中国存储量较丰富且仅次于石油的油气资源[1]。为了满足中国日益增长的能源需求,大力发展页岩产业和页岩开采技术是十分必要的。射流破岩技术能够大面积且有导向地碎裂储层岩体,方便页岩气和煤层气等油气资源的解吸与渗流,已被广泛应用于页岩的开采领域[2-3]。
目前,国内外学者已对磨料水射流破岩问题展开了大量的研究[4-6],但由于射流和岩体具有不同的物理属性和力学性能,尤其是岩体中还存在着大量的孔隙[7],再加上磨料颗粒对岩体的冲蚀作用,使得岩体的变形破损机制变得极其复杂[8]。因此,开展磨料水射流冲击破岩的研究,并对含有孔隙结构岩体的破损过程进行分析,能够为射流破岩开采技术的优化和完善提供一定的理论依据。
随着计算机技术的发展,数值模拟方法已经成为研究射流破岩问题的有效手段之一[9-10]。黄中伟等[11]基于有限元法模拟了液氮磨料射流冲击岩体的过程,发现冲击是岩体破损的主要因素。庄欠伟等[12]应用有限元方法模拟研究了磨料水射流冲击钢筋和混凝土的过程。穆朝民等[13]使用网格法构建了岩体,并分析了在受到磨料水射流冲击后岩体的损伤演化过程。由于网格类的数值方法在模拟大变形等问题时存在着网格畸变的情况,因此,近年来无网格方法得到了大量的关注并广泛应用于射流破岩问题的研究。
光滑粒子流体动力学(smoothed particle hydro- dynamics,SPH)方法作为一种典型的无网格粒子方法,在处理大变形、应力局部化等问题中具有一定的优势,被广泛采纳并应用于冲击大变形和流固耦合等领域。刘勇等[14]运用SPH方法模拟分析了磨料空气射流破岩的冲蚀过程。李世杰等[15]采用SPH方法模研究了不同直径的水射流对土壤冲蚀坑的影响。赵健等[16]将SPH方法和有限元方法进行耦合,模拟分析了磨料水射流冲击下岩体的损伤演化机制。虽然目前有许多学者已采用SPH方法模拟研究了不同的射流破岩问题,但针对射流冲击下含孔隙岩体的破损过程以及孔隙结构对破岩效果影响的研究相对较少。鉴于实际岩体中含有孔隙结构,本文在岩体模型的内部引入了孔隙,模拟研究了磨料水射流冲击含孔隙岩体的破损过程,并分析了多种孔隙特性对射流破岩效果的影响。
本文依据SPH方法的基本理论,分别引入描述射流和岩体力学特性的NULL型和H-J-C型本构模型,建立磨料水冲击含孔隙岩体的数值模型。通过模拟磨料水射流破岩问题,验证了方法和模型的有效性。在此基础上,加入随机混合的孔隙,构建出磨料水射流冲击含孔隙结构岩体的动力学模型,对射流破岩的动态过程进行模拟研究,分析了岩体的孔隙特性(孔隙尺寸、孔隙度和孔隙填充物)以及射流参数(射流速度和磨料浓度)对岩体破损效果的影响,为深入理解射流开采中岩体的破损机理和提高开采效率提供了一定的理论依据。
1. 模型的建立
1.1 SPH方法的基础理论
SPH方法是一种通过利用离散点来构造近似函数的纯拉格朗日方法。对于磨料水射流冲击含孔隙岩体的问题而言,可将含磨料颗粒的水射流和具有孔隙结构的岩体离散为一系列的粒子,通过差值求解非稳态N-S控制方程即可得到粒子的物理属性。因此,宏观变量可通过计算域内离散点的积分插值得到,质点的近似函数为[17]
f(x)=∫Ωf(x′)W(x−x′,h)dx′。 (1) 式中f(x)为三维坐标向量x的函数;Ω为x的支持域;x−x′为粒子间距离;h为SPH粒子光滑长度,光滑长度随时间和空间变化,如图 1所示;W(x−x′,h)为核函数,用辅助函数θ(x−x′)定义:
W(x−x′,h)=1h(x−x′)dθ(x−x′), (2) 式中,d为空间维数。用粒子近似法将连续形式积分方程换算为离散形式的方程,即
f(x)=n∑i=1miρi1h(x−x′)dθ(x−x′), (3) 式中,ρi为粒子i的密度,mi为粒子i的质量。
SPH方法下的N-S方程中质量守恒方程为
ρi=n∑j=1mjWijN∑j=1mjρjWij。 (4) 动量守恒方程为
DvαiDt=N∑j=1mj(σαβiρ2i+σαβjρ2j)∂Wij∂xβi。 (5) 能量守恒方程为
DeiDt=12N∑j=1mj(Pi+Pjρiρj)vβij∂Wij∂xβi+μi2ρiεαβiεαβj, (6) 式中,v为速度,e为热能,σ为应力张量,P为压力,m为质量,ρ为密度,μ为动力黏度系数,ε为应变率,N为质点总数,x为坐标轴符号,α,β代表不同坐标轴,i,j代表不同质点。
为了保证计算中的精度和稳定性,提高计算效率,研究采用了应用较广的三次B样条核函数[17]。
1.2 磨料和水的本构模型
研究采用了LS-DYNA中的NULL本构模型描述水和磨料,塑性材料是一种在外力的作用下,能够产生显著变形而不被破坏的材料,由于在射流破岩问题中,磨料和水较符合塑性材料的特性[2],因此本文在仿真模拟中将两者都视为完全塑性材料并赋予Gruneisen状态方程[18],即
p=ρ0C2μ[1+(1−γ02)μ−a2μ2][1−(S1−1)μ−S2μ2μ+1−S3μ3(μ+1)2]+(γ0+aμ)Ea。 (7) 式中C为冲击波速度与质点速度变化曲线截距;S1,S2,S3为冲击波与质点速度变化曲线的斜率系数;γ0为Gruneisen常数;a为γ0和μ=ρ/ρ0−1的一阶体积修正量。水和磨料颗粒的本构模型参数[16]如表 1所示。为了实现磨料粒子在流场中的真实分布情况,从射流粒子中进行随机抽样,并赋予磨料的属性,实现射流中磨料颗粒的随机混合。
表 1 水和磨料本构模型的参数Table 1. Parameters of constitutive model for water and abrasive名称 ρ0 /(g·cm-3) C/(m·s-1) S1 S2 S3 水 1.0 1480 2.56 -1.986 0.228 磨料 3.5 4569 1.49 0 0.460 1.3 岩体的本构模型
在射流的冲击作用下,不同性质的岩体会呈现出不同程度的破损。研究选用页岩材料,并使用H-J-C本构模型描述,页岩本构模型的参数[19]如表 2所示。H-J-C本构模型综合考虑了岩石材料损伤、应变率和静水压力对屈服力的影响,是一种能够在岩体大变形、高应变率和高压等情况下应用的模型,其等效强度通过应力、应变速率和损伤状态来描述,屈服面的方程为[20]
σ∗=[A(1−D)+Bp∗N](1−Clnˉε∗)。 (8) 表 2 页岩本构模型的参数Table 2. Parameters of constitutive model for shale岩石 ρ/(g·cm-3) G/GPa fc/GPa T/GPa μcrush 页岩 2.44 13.6 0.0048 0.004 0.001 式中σ∗为岩体冲击载荷作用下等效应力与静态屈服强度之比;P∗为岩体实际压力与静态屈服强度之比;¯ε∗为冲击载荷作用下岩体应变率与静态应变率之比;A,B,N和c为材料的强度参数。
损伤因子D通过等效塑性应变和塑性体积应变累计得到
D=∑Δεp+ΔμpD1(p∗+T∗)D2。 (9) 式中Δεp和Δμp分别为一个计算循环内的等效塑性应变和塑性体积应变;D1与D2为损伤常数;T∗=T/fe,T为材料的最大拉伸强度,fe为材料静态抗压强度。
1.4 孔隙的构建
为了构建与天然岩体更接近并且能够进行数值计算的孔隙结构,研究中参考最大球算法[21]建立孔隙的过程和随机堆积球模型[22-24],将孔隙离散为球体,允许孔隙相互重叠,重叠程度与孔隙自身相比较小,如图 2所示。
在使用SPH方法将岩体与孔隙的离散后,空间中的任意一个代表孔隙中心位置的粒子,与最近岩体粒子的欧氏距离,为球心与岩体表面相切的最大球的半径,计算公式如下:
R=√(x1−x2)2+(y1−y2)2+(z1−z2)2, (10) 式中,(x1,y1,z1)为孔隙中心点的坐标,(x2,y2,z2)为岩体点的坐标。
根据岩体与孔隙的质量比α能够确定出页岩的孔隙度β,之后通过β即可确定出孔隙的个数:
α=mamw+ma=VaρaVwρ0+Vaρa, (11) β=VwVa=αρ01−αρ0, (12) na=3dw3β4π R3。 (13) 式中ma,mw分别为岩体与孔隙的质量;Va,Vw分别为岩体与孔隙的体积;ρa,ρ0分别为岩体与孔隙的密度;na为页岩模型中孔隙的个数;dw为岩体的边长,R为孔隙最大内切球的半径。
页岩等沉积岩的形成过程十分复杂,因此在建立孔隙结构需要有以下特征[25]:①随机性。在岩石的沉积过程中,孔隙在岩体中随机生成,且在任何点都有均等形成的机会。②独立性。孔隙生成在岩体中达到平衡且静止的状态,不会受到其他孔隙的影响。③长周期性。在岩体沉积、形成的过程中,孔隙生成的方向为重力势能梯度最大的方向,选取的孔隙的位置需要符合此规律。研究中采用了蒙特卡洛法随机生成孔隙,对于空间内任意一个SPH粒子,具有相同的概率被采样,能够使需要生成的na个孔隙均匀的分布在岩体内部,被选取的单个孔隙位置表示为Pi(x,y,z),其描述如下:
Pi(x,y,z)={x,y,z=Random(L)},i⩽na, (14) 式中,Random(L)是能够满足孔隙分布的随机函数,x,y,z为坐标信息。
在随机得到的Pi(x,y,z)位置上生成一个SPH粒子之后,沿着X,Y,Z三个方向继续增加粒子,最终生成单个孔隙,演化过程如图 3所示。
本文离散球体的规则与SILIN等[21]和文献[22,23]一致。为了比较不同尺寸的孔隙之间的差异,建立球体大小相等的孔隙模型,如图 4所示。孔隙的本构模型及状态方程与水相同
经过对比,发现在相同参数下的构建的岩体仿真模型与真实岩体的孔隙分布[26]基本相同,并且通过三轴压缩试验,得到了在载荷作用下孔隙岩体模型的应力应变曲线,如图 5所示。经计算得知,该模型的弹性模量与文献[27]中真实孔隙岩体相差1%,其他力学参数基本一致,说明了构建的孔隙岩体模型能够非常接近真实岩体孔隙结构。
1.5 磨料水射流冲击孔隙岩体模型
针对所研究的磨料水射流冲击含孔隙岩体的问题,进行如下的简化假设:①磨料水射流破岩过程中忽略气相环境的影响,只包含水射流、磨料颗粒、岩体和孔隙4种物质;②磨料颗粒与水粒子为相同直径的球体。
为了保证计算效率,采用了1/4的对称模型,对称面设置为SPH_SYMMETRY_PLANE对称约束,并在对称面的基础上增加SPH_NON_REFLECTING约束来模拟岩石无限大的情况。对称后岩体为50 mm×50 mm×25 mm的立方体,孔隙的直径介于0.3~2 mm,射流喷距为1 mm,模型共由424395个光滑粒子构成,其中2394个为射流粒子,包含了水与磨料。磨料水射流冲击孔隙岩体的几何模型如图 6所示。
研究中选取了影响岩体破损效果的以下因素进行分析:孔隙特性(孔隙尺寸ϕ、孔隙度β和孔隙填充物)以及射流参数(射流速度v和磨料浓度δ)。
1.6 模型的验证
为了验证所建模型的有效性,对李井慧[28]所研究的磨料水射流冲击大理石过程进行了模拟,各项参数均与文献[28]保持一致,图 7为不同射流压力下破损坑截面演化图,由图 7可以看出,随着射流压力的增加,大理石的侵彻深度在逐渐的增大。
图 8为试验与模拟结果对比图。由图 8可知,由于真实试验的冲击时间大于仿真模拟,且试验中的磨料水射流能够持续冲击岩体,导致试验切割深度整体大于模拟结果,但岩体的切割深度随射流压力的变化趋势基本相同,均随着射流压力的增加而增大,验证了本文所构建的磨料水射流模型的有效性。
本文还对林晓东等[29]所研究的磨料水射流冲击灰岩问题进行了模拟,对岩体破损形状进行分析和验证。射流直径为4 mm,射流高度为75 mm,磨料浓度为30%,均与文献[29]保持一致。射流冲击灰岩破损坑截面演化过程如图 9所示,当射流撞击岩体的瞬间,受到冲击的区域出现了破损坑,随着冲击时间的增加,破损坑逐渐扩大,最终接近子弹型。由于林晓东等[29]采用网格法描述岩体,受到冲蚀作用的岩石网格单元破坏失效后删除,而本文采用SPH方法对岩体进行了离散,使得破损坑的形状略有差异,但冲孔截面形态基本相同,与试验的冲孔扫描结果[29]吻合,说明使用SPH方法替代网格法构建岩体是可行的。
图 10给出了射流速度模拟结果对比,由图 10可知,当射流速度小于200 m/s时,冲击速度对破损坑半径的影响较小;当射流速度超过200 m/s时,破损坑半径随着射流速度的增加急剧增大,不同模拟方法所得的岩体冲孔半径随冲击速度的变化趋势较为相似,进一步验证了本研究基于SPH方法建立的数值模型的有效性。后续关于磨料水射流冲击含孔隙岩体的研究以此模型为基础,为了更好的观察破损坑形态的变化,射流高度调整为25 mm。
2. 破岩指标参数定义
为了全面分析岩体的破损情况,对不同工况下的破岩过程进行了模拟,并提取了岩体破损图像的特征,研究表明破损坑深和破损面积的变化规律基本一致,因此本文选取了表面坑径(r/R)和破损面积(s/S)为岩体破损指标,并进行了无量纲处理,如图 11所示,其中r,s分别为岩体破损坑的表面坑径和纵截面积,R,S分别是岩体的宽度和纵截面积。
3. Von-Mises应力分析
选取孔隙度β=10%,孔隙尺寸为1 mm3,孔隙填充物为水,射流速度v=500 m/s,磨料浓度δ=5%的数值模型进行Von-Mises应力分析,图 12为射流冲击过程的Von-Mises应力演化图。
由图 12可知,射流撞击到岩体的瞬间,受到冲击的岩体区域出现了较大的应力集中现象,当Von-Mises应力超过岩体的强度极限时,会导致岩体与射流的接触区域破损。由于岩体部分采用了脆性较高的页岩,当t=0.03 ms时,应力波开始呈现出裂纹状;距离射流较近的孔隙边缘出现了应力集中现象,并且孔隙的形状发生了改变,因孔隙填充物为水,导致孔隙内部应力较小。随着冲击时间的增加,冲击区域两侧的孔隙持续受到岩体挤压。当t≥0.06 ms时,孔隙被完全压碎,导致无法看出孔隙具体的形状和位置,岩体的破损程度不断加剧,应力波在岩体内部向四周传播并且以裂纹状继续扩散。
4. 孔隙特性及射流参数的影响
4.1 孔隙尺寸
孔隙的尺寸会影响射流破岩的效率和岩体的破损程度。研究中取孔隙尺寸ϕ=1~6(以ϕ=3为例说明,孔隙直径离散为3个粒子),如图 13所示。孔隙度β=15%,孔隙填充物为水,磨料体积浓度δ=5%,射流速度v=500 m/s。
图 14给出了不同孔隙尺寸下t=0.1 ms时岩体破损区域与破损指标随冲击时间的变化曲线。由图 14(a)可看出,在相同的孔隙度下,当孔隙尺寸较小时(ϕ≤4),岩体内孔隙数量较多,能够均匀的分布在岩体中,随着孔隙体积的增大,岩体破损加剧;当孔隙尺寸ϕ≥5时,岩体的破损坑的形状变化较大。由图 14(b)可以看出,随着射流冲击时间的增加,岩体破损坑的表面坑径和破损面积都逐渐增大。当孔隙较小时(ϕ≤4),随着孔隙尺寸的增加,破损面积不断增大,但表面坑径之间的差异较小;当孔隙尺寸ϕ≥5时,表面坑径之间出现了较大的差异。
为了探究岩体的破损是否受到孔隙分布位置的影响,研究中对孔隙尺寸ϕ≤4的4种情况进行了多次模拟并进行了对比,发现在相同孔隙尺寸下,岩体的破损情况受到孔隙分布位置的影响较小。在孔隙尺寸较大的两种情况下(ϕ≥5),分别增加3种不同的孔隙随机分布。图 15为不同的孔隙分布位置下t=0.1 ms时岩体破损区域与破损指标随冲击时间的变化曲线。结合图 14,15可以得知,当孔隙较大并且在相同的孔隙度下,不同孔隙分布位置的破损坑形状有显著的差别。当孔隙尺寸ϕ≥5时,不同的孔隙分布位置下破损坑的表面坑径以及破损面积差异较大,说明了在磨料水射流冲击的过程中含大孔隙岩体的破损受到孔隙随机分布位置的影响,因此后续的研究中孔隙尺寸选取为1 mm3(ϕ=3)。
4.2 孔隙度
孔隙度是影响岩体破损的重要原因。研究中取孔隙度β=0~25%,孔隙尺寸为1 mm3(ϕ=3),孔隙填充物为水,射流速度v=500 m/s,磨料体积浓度δ=5%。图 16给出了不同孔隙度下t=0.1 ms时岩体破损区域与破损指标随冲击时间的变化曲线。从图 16可以看出,随着岩体内部孔隙的增多,在磨料水射流的冲击作用下岩体破损坑的表面坑径和破损面积均增大,尤其是破损面积急剧增加,说明孔隙率越高的岩体越容易破碎。当岩体内不含孔隙(β=0%)时,岩体的破损面积略大于含孔隙较少的岩体(0%<β≤10%),这也说明少量的孔隙及其填充物可能会对磨料水射流产生干扰并起到一定的缓冲作用,降低了射流的冲击破岩能力。
4.3 孔隙填充物
不同区域的页岩孔隙中蕴藏着不同的填充物质。研究选取了水、轻油、重油和空气4种孔隙填充物,由于空气的密度与液体相比非常小,可近似忽略,因此将孔隙粒子删除视为空气,填充物的参数[30]如表 3所示,孔隙尺寸为1 mm3(ϕ=3),孔隙度β=15%,孔隙的分布位置完全相同;射流速度v=500 m/s,磨料体积浓度δ=5%。
表 3 孔隙填充物的参数Table 3. Parameters of pore fillers材料 ρ /(g·cm-3) 动力黏度系数/(pa·s) 水 1.00 1.009×10-3 重油 1.10 2.5×10-3 轻油 0.85 2.5×10-3 图 17给出了不同填充物下t=0.1 ms时岩体破损区域与破损指标随冲击时间的变化曲线。由图 17可以看出,当填充物为液体时,岩体破损坑的形状较为相似,表面坑径以及破损面积差异较小;当孔隙填充物为空气时,岩体破损坑的面积最大,但表面坑径最小,这也说明与其他填充物相比,含气体的岩体更容易沿着磨料水射流的冲击方向发生破损。
4.4 射流速度
射流速度是决定开采破损效率的重要因素。射流速度取v=100~500 m/s,磨料浓度δ=5%,孔隙尺寸为1 mm3(ϕ=3),孔隙度β=15%,孔隙填充物为水。图 18给出了不同射流速度下t=0.1 ms时岩体的破损区域与破损指标随冲击时间的变化曲线。由图 18可以看出,岩体破损坑的表面坑径和破损面积随着射流速度的增大而增加。当射流速度v>200 m/s时,表面坑径和破损面积急剧增加,说明磨料水射流速度越高,破损孔隙岩体能力越强。
4.5 磨料浓度
磨料浓度也是影响射流破岩效率和破损精度的重要因素之一。根据实际磨料水射流破岩过程中的磨料颗粒体积占比,磨料浓度δ=0~25%,射流速度v=500 m/s,孔隙尺寸为1 mm3(ϕ=3),孔隙度β=15%,孔隙填充物为水。图 19给出了不同磨料浓度下t=0.1 ms时岩体的破损区域与破损指标随冲击时间的变化曲线。
由图 19可知,随着磨料浓度的增加,岩体破损坑的纵截面形状发生了明显的改变。当水射流中添加磨料颗粒后,岩体破损坑的表面坑径急剧减小,磨料浓度越大,岩体破损坑的表面坑径和破损面积越小;当磨料浓度增大到10%以上时,各浓度之间的表面坑径以及破损面积差异较小。
为了深入分析磨料浓度对岩体破损效果的影响,给出了底部坑径(l/R)指标,如图 20所示。
图 21为不同磨料浓度下底部坑径随冲击时间的变化曲线,由图 21可知,当水射流中混有磨料颗粒时,岩体破损坑的底部坑径急剧增加,并且随着磨料浓度的增加而增大。结合图 19(a)可知,磨料颗粒主要集中在破损坑四周的表面上并产生冲蚀作用,使破损坑的表面坑径和底部坑径相差较小,说明添加磨料后能够有效的提高射流破岩的精准性。
5. 结论
(1)模拟分析了射流冲击岩体过程的破损形态和指标,通过与已有研究结果进行对比,验证了本研究所建模型的有效性。
(2)文中采用PYTHON语言编写离散后的球体孔隙,并实现随机分布,采用该方法可以对数值模拟中的孔隙特性实现精准的控制。
(3)模拟研究了孔隙特性和射流参数对岩体破损的影响。孔隙尺寸和孔隙分布情况会显著影响岩体的破损形状;当孔隙度介于0~25%时,相同的孔隙尺寸下,岩体孔隙度越大,岩体破损坑的表面坑径和破损面积越大;当孔隙填充物之间的密度以及动力黏度系数相对接近时,填充物对岩体破损的影响相对较小,其中当填充物为气体时,岩体则更容易沿着射流冲击方向发生破损变形。射流的初始速度和磨料浓度也是影响岩体破损情况的主要因素,射流初速度越大,射流的冲击能力越强,岩体破损坑的表面坑径和破损面积也随之增大;对于多孔岩体而言,冲击过程中的磨料颗粒主要集中在破损坑四周的壁面上,随着磨料浓度的增大,磨料颗粒对破损坑壁面的冲蚀效果越强,坑底部的直径明显增大。
-
表 1 泰勒级数展开式未知数的取值
Table 1 Values of unknowns in Taylor series expansion
元素 C0/(μg·mL-1) Cl/(μg·mL-1) D/(cm2·s-1) r/cm Al 0.162 0.383 0.115 1.547 Ca 3.5×10-15 0.350 0.235 1.615 Fe 0.042 0.192 0.144 1.491 K 0.076 0.245 0.141 1.248 Mg 0.017 0.233 0.117 1.263 Na 1.0×10-14 1.384 0.199 1.640 Si 1.0×10-14 2.360 0.167 1.605 表 2 Voronoi多面体块细观参数
Table 2 Mesoscopic parameters of Voronoi blocks
时间/d 密度/(kg·m-3) 弹性模量/GPa 泊松比 0 2400 1.31 0.27 1 2400 1.31 0.27 4 2400 1.31 0.27 8 2400 1.31 0.27 15 2400 1.31 0.27 表 3 接触单元细观参数
Table 3 Mesoscopic parameters of contact elements
时间/d 法向接触刚度/(1012 Pa·m-1) 切向接触刚度/(1012 Pa·m-1) 黏聚力/MPa 内摩擦角
/(°)抗拉强度/MPa 0 1.35 1.31 1.75 30 0.76 1 0.80 0.80 0.77 30 0.34 4 0.36 0.35 0.11 30 0.05 8 0.32 0.31 0.08 30 0.04 15 0.29 0.28 0.06 30 0.03 -
[1] 朱俊杰. 滇中红层软岩水–岩作用机理及时效性变形特性研究[D]. 成都: 成都理工大学, 2019. ZHU Jun-jie. Study on Water-Rock Mechanism and Time-Dependent Deformation Characteristics of Red Beds Soft Rock in Central Yunnan[D]. Chengdu: Chengdu University of Technology, 2019. (in Chinese)
[2] 何满潮. 软岩工程力学[M]. 北京: 科学出版社, 2002. HE Man-chao. Engineering Mechanics of Soft Rock[M]. Beijing: Science Press, 2002. (in Chinese)
[3] 冯启言, 韩宝平, 隋旺华. 鲁西南地区红层软岩水岩作用特征与工程应用[J]. 工程地质学报, 1999, 7(3): 266–271. doi: 10.3969/j.issn.1004-9665.1999.03.012 FENG Qi-yan, HAN Bao-ping, SUI Wang-hua. Characteristics of water rock interaction of red beds and its application to engineering in southwestern Shandong[J]. Journal of Engineering Geology, 1999, 7(3): 266–271. (in Chinese) doi: 10.3969/j.issn.1004-9665.1999.03.012
[4] 周翠英, 谭祥韶, 邓毅梅, 等. 特殊软岩软化的微观机制研究[J]. 岩石力学与工程学报, 2005, 24(3): 394–400. doi: 10.3321/j.issn:1000-6915.2005.03.006 ZHOU Cui-ying, TAN Xiang-shao, DENG Yi-mei, et al. Research on softening micro-mechanism of special soft rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(3): 394–400. (in Chinese) doi: 10.3321/j.issn:1000-6915.2005.03.006
[5] 程强. 红层软岩开挖边坡致灾机理及防治技术研究[D]. 成都: 西南交通大学, 2008. CHENG Qiang. Study on Hazard Mechanism and Prevention Technology of Cutting Slope in Red Beds Soft Rock[D]. Chengdu: Southwest Jiaotong University, 2008. (in Chinese)
[6] 刘光廷, 胡昱, 李鹏辉. 软岩遇水软化膨胀特性及其对拱坝的影响[J]. 岩石力学与工程学报, 2006, 25(9): 1729–1734. doi: 10.3321/j.issn:1000-6915.2006.09.001 LIU Guang-ting, HU Yu, LI Peng-hui. Behavior of soaking rock and its effects on design of arch dam[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(9): 1729–1734. (in Chinese) doi: 10.3321/j.issn:1000-6915.2006.09.001
[7] 李洪志, 何满潮. 膨胀型软岩力学化学性质研究[J]. 煤, 1995(6): 9–12. https://www.cnki.com.cn/Article/CJFDTOTAL-MEIA506.002.htm LI Hong-zhi, HE Man-chao. Study on mechanical and chemical properties of swelled soft rock[J]. Coal, 1995(6): 9–12. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MEIA506.002.htm
[8] MORGENSTERN N R, EIGENBROD K D. Classification of argillaceous soils and rocks[J]. Journal of the Geotechnical Engineering Division, 1974, 100(10): 1137–1156. doi: 10.1061/AJGEB6.0000106
[9] 吴道祥, 刘宏杰, 王国强. 红层软岩崩解性室内试验研究[J]. 岩石力学与工程学报, 2010, 29(增刊2): 4173–4179. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S2104.htm WU Dao-xiang, LIU Hong-jie, WANG Guo-qiang. Laboratory experimental study of slaking characteristics of red-bed soft rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 4173–4179. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S2104.htm
[10] 周翠英, 彭泽英, 尚伟, 等. 论岩土工程中水-岩相互作用研究的焦点问题——特殊软岩的力学变异性[J]. 岩土力学, 2002, 23(1): 124–128. doi: 10.3969/j.issn.1000-7598.2002.01.028 ZHOU Cui-ying, PENG Ze-ying, SHANG Wei, et al. On the key problem of the water-rock interaction in geoengineering: mechanical variability of special weak rocks and some development trends[J]. Rock and Soil Mechanics, 2002, 23(1): 124–128. (in Chinese) doi: 10.3969/j.issn.1000-7598.2002.01.028
[11] 张丹, 陈安强, 刘刚才. 紫色泥岩水热条件下崩解过程的分维特性[J]. 岩土力学, 2012, 33(5): 1341–1346. doi: 10.3969/j.issn.1000-7598.2012.05.010 ZHANG Dan, CHEN An-qiang, LIU Gang-cai. Disintegration characteristics of purple mudstone based on fractal dimension under hydrothermal condition[J]. Rock and Soil Mechanics, 2012, 33(5): 1341–1346. (in Chinese) doi: 10.3969/j.issn.1000-7598.2012.05.010
[12] 黄明, 詹金武. 酸碱溶液环境中软岩的崩解试验及能量耗散特征研究[J]. 岩土力学, 2015, 36(9): 2607–2612, 2623. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201509023.htm HUANG Ming, ZHAN Jin-wu. Disintegration tests and energy dissipation characteristics of soft rock in acid and alkali solution[J]. Rock and Soil Mechanics, 2015, 36(9): 2607–2612, 2623. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201509023.htm
[13] 杨建林, 王来贵, 李喜林, 等. 遇水–风干循环作用下泥岩断裂的微观机制研究[J]. 岩石力学与工程学报, 2014, 33(增刊2): 3606–3612. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2014S2028.htm YANG Jian-lin, WANG Lai-gui, LI Xi-lin, et al. Research on micro-fracture mechanism of mudstone after wet-dry cycles[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S2): 3606–3612. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2014S2028.htm
[14] 沈照理, 王焰新. 水–岩相互作用研究的回顾与展望[J]. 地球科学, 2002, 27(2): 127–133. doi: 10.3321/j.issn:1000-2383.2002.02.001 SHEN Zhao-li, WANG Yan-xin. Review and outlook of water-rock interaction studies[J]. Earth Science, 2002, 27(2): 127–133. (in Chinese) doi: 10.3321/j.issn:1000-2383.2002.02.001
[15] 周翠英, 黄思宇, 刘镇, 等. 红层软岩软化的界面过程及其动力学模型[J]. 岩土力学, 2019, 40(8): 3189–3196, 3206. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201908036.htm ZHOU Cui-ying, HUANG Si-yu, LIU Zhen, et al. The interface process and its dynamic model of red-bed soft rock softening[J]. Rock and Soil Mechanics, 2019, 40(8): 3189–3196, 3206. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201908036.htm
[16] 潘艺, 刘镇, 周翠英. 红层软岩遇水崩解特性试验及其界面模型[J]. 岩土力学, 2017, 38(11): 3231–3239. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201711020.htm PAN Yi, LIU Zhen, ZHOU Cui-ying. Experimental study of disintegration characteristics of red-bed soft rock within water and its interface model[J]. Rock and Soil Mechanics, 2017, 38(11): 3231–3239. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201711020.htm
[17] 陈瑜, 曹平, 蒲成志, 等. 水–岩作用对岩石表面微观形貌影响的试验研究[J]. 岩土力学, 2010, 31(11): 3452–3458. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201011017.htm CHEN Yu, CAO Ping, PU Cheng-zhi, et al. Experimental study of effect of water-rock interaction on micto-topography of rock surface[J]. Rock and Soil Mechanics, 2010, 31(11): 3452–3458. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201011017.htm
[18] CIANTIA M O, CASTELLANZA R, CROSTA G B, et al. Effects of mineral suspension and dissolution on strength and compressibility of soft carbonate rocks[J]. Engineering Geology, 2015, 184: 1–18.
[19] 卢良兆, 许文良. 岩石学[M]. 北京: 地质出版社, 2011. LU Liang-zhao, XU Wen-liang. Lithology[M]. Beijing: Geological Publishing House, 2011. (in Chinese)
[20] GU D M, HUANG D A, ZHANG W G, et al. A 2D DEM-based approach for modeling water-induced degradation of carbonate rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2020, 126: 104188–104188.
[21] 付腾飞, 徐涛, 朱万成, 等. 基于多晶离散元法的砂岩三轴压缩损伤特性[J]. 东北大学学报(自然科学版), 2020, 41(7): 968–974. https://www.cnki.com.cn/Article/CJFDTOTAL-DBDX202007010.htm FU Teng-fei, XU Tao, ZHU Wan-cheng, et al. Damage compression based on polycrystalline discrete element method[J]. Journal of Northeastern University (Natural Science), 2020, 41(7): 968–974. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DBDX202007010.htm
[22] GHAZVINIAN E, DIEDERICHS M S, QUEY R. 3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2014, 6(6): 506–521.
-
期刊类型引用(8)
1. 王泽志,王杰,马小刚,李帆,范新亚,陈燕. 用磁性微磨料射流技术光整加工交叉深孔内壁. 金刚石与磨料磨具工程. 2025(02): 245-255 . 百度学术
2. 王高辉,孔维伟,卢文波,潘鑫豪,舒奕展. 高拱坝坝后侵彻爆炸毁伤效应分析. 武汉大学学报(工学版). 2024(07): 853-862 . 百度学术
3. 王义江,郁东旭,孙立鹏,朱启银,王建州. 激光照射岩石热裂特性与裂隙分布研究. 岩土工程学报. 2024(09): 1809-1819 . 本站查看
4. 韩明兴,余锴,段宏兵,熊利荣,徐琨,李淼,刘琦. 基于SPH-FEM的油菜茎秆磨粒气体射流切割仿真与试验. 农业工程学报. 2024(17): 82-92 . 百度学术
5. 李新旺,王汉青,程立朝,张学栋,温学君,申鹏飞. 高压水射流冲击下煤体初始失效形式分析. 中国科技论文. 2023(01): 64-70+90 . 百度学术
6. 胡阳,曹安业,秦续峰,薛成春,郭文豪,刘耀琪,彭雨杰. 深部冲击地压煤层水射流割缝卸压参数优化研究. 岩土工程学报. 2023(07): 1509-1516 . 本站查看
7. 林琳,赵志磊,蒋东岑,张云朋,尤晖. 尖边喷嘴磨料水射流特性分析及对去除函数轮廓的影响. 液压与气动. 2023(09): 116-129 . 百度学术
8. 王璐,徐绯,杨扬. 完全拉格朗日SPH在冲击问题中的改进和应用. 力学学报. 2022(12): 3297-3309 . 百度学术
其他类型引用(9)
-
其他相关附件
-
本文附件音频
论文音频介绍 点击下载(901KB) -
DOCX格式
论文补充图文材料 点击下载(2152KB) -
PDF格式
审稿意见与作者答复 点击下载(387KB)
-