Numerical simulation of abrasive water jet impacting porous rock based on SPH method
-
摘要: 针对磨料水射流冲击含孔隙岩体的动态过程，依据光滑粒子流体动力学（smoothed particle hydrodynamics, SPH）方法，将含磨料颗粒的水射流和具有孔隙的岩体离散为一系列的粒子，分别引入NULL型和H-J-C型本构关系描述磨料水射流和岩体的力学特性，建立了磨料水射流冲击孔隙岩体模型，通过模拟射流冲击岩体问题，并与已有研究结果进行对比验证了所建模型的有效性，在此基础上，着重对磨料水射流冲击含孔隙岩体的动态破损过程进行模拟研究，分析了岩体的孔隙特性和射流参数对岩体破损效果的影响。计算结果表明，孔隙的尺寸和分布情况会影响岩体破损坑的形状，孔隙度对岩体的破损程度影响较大，其他因素对岩体的破损效果具有不同程度的影响。研究结果可为工程中提高射流破岩效率提供理论依据。Abstract: Shale contains pores. The abrasive water jet technology is used to extract oil and gas resources inside the shale, and the impact process is affected by the pores. For the damage process of abrasive water jet impacting the pore-containing rock, according to the basic principle of smoothed particle hydrodynamics (SPH) method, the abrasive water jet and the pore-containing rock are discretized into a series of particles, and a dynamic model for the abrasive water jet impacting the pore-containing rock is established by the constitutive relationship of NULL and H-J-C to describe the mechanical properties of abrasive water jets and rock respectively. The effectiveness of the established model is verified and analyzed by simulating the impact of abrasive water jet on the rock and comparing with the existing research results. On this basis, the dynamic damage process of the abrasive water jet impacting the pore-bearing rock is simulated, and the effects of pore characteristics and jet parameters on the damage degree of rock are also analyzed. The simulated results show that the velocity of abrasive water jet and porosity of rock have significant effects on the breaking process of the rock, while the size and distribution of the pores can affect the breaking shape of the rock, and other factors have different influences on the damage degree of rock. The research conclusion can provide theoretical support for improving the rock-breaking efficiency of water jet in engineering.
-
0. 引言
长距离引调水工程是一项解决中国北方部分地区水资源匮乏所采用的规模巨大且意义深远的水利基础建设项目[1]。近年来,输水渠道和堤坝工程规模越来越大,从几公里到几十、几百甚至是上千公里,工程条件复杂、环境恶劣[2]。引额济乌一期项目位于新疆北部,是一项跨流域、长距离的调水工程。输水渠道在冻融循环、季节性输水、高纬度强光照等影响因素下可能会存在渗漏安全隐患,降低渠堤的渗流稳定性,导致溃堤等安全事故,将会对沿线地区的经济和人民财产安全造成不可估量的损失。
渠堤渗漏监测通常采用渗压计、测压管、量水堰等监测仪器对构筑物的代表性点进行监测,通过单点的测值来反映整体渗流情况。但出于成本考虑,渠道监测断面布置稀疏,受仪器和人员随机误差影响大,极易出现较长距离的漏测、误测的情况,给长距离渠道的安全运行管理带来了挑战。基于分布式光纤温度传感技术,通过监测渠堤不同位置的温度变化值来推算其渗流情况是一种新型的渗流监测技术[3-4]。与传统渗漏监测方法相比较,分布式光纤测渗技术具有抗高压、抗电磁干扰、耐腐蚀、耐高温、防雷击、成本低廉、施工过程简单等特点,非常适合在岩土体实现长距离分布式监测。随着分布式光纤解析技术的日益完善,大量的学者开展了各类基于分布式光纤的测温系统研究[5-7],促进了该技术日趋成熟。
本文提出了一种基于拉曼光时域反射原理的新型渗流监测光缆型式,并研发了可用于长距离输水渠道的分布式光纤拉曼渗流监测技术,通过理论研究、室内试验以及现场应用,验证了该技术的适用性。
1. 高寒区渠道渗漏传感光缆设计
现有的DTS技术（分布式光纤温度传感技术）是将发热材料产生的热量直接传递给测温光纤,当外界环境存在渗流时,部分热量随渗流向外散失,导致测温光纤接收到的热量减少,测点温度低于无渗流的情形,由此来判断渗流的存在。但在实际工程中,由于测温光纤在渗流存在的情况下依旧能从发热材料接收到较多的热量,温度变化程度较小,敏感度低[8]。本文设计并制作了采用特殊封装技术和护套材料的测量光缆。该光缆结构为带加热功能双层四纤非铠装拉曼光缆（光缆结构如图1）,该结构可以抵御北疆输水干渠现场零下20℃的低温,保证光缆安全过冬。
图1中1为隔热橡胶圆柱、21为1号光纤套管、22为2号光纤套管、31为1号发热导线、32为2号发热导线、4为氟塑料层、5为钢带、6为导热外护套。
2. 光纤渗漏监测室内试验
2.1 室内试验设计
室内试验设计将光缆置于不同的预设环境（包括：空气、干沙、含水率为5%,10%,15%,20%,30%和饱和的砂土、纯水）中以验证光缆的适用性。试验装置如图2所示。DTS解调仪是解析光缆温度信息,并将数据传输给电脑的重要仪器设备；环境模拟箱是大型水箱,用来放置试验所需的模拟环境。
试验时将模拟环境加入环境模拟箱,再将光缆埋设于模拟环境中,开展各种环境下的光缆升温试验。
2.2 室内试验结果分析
综合对比与分析1号光纤与2号光纤的温升状态,如图3与图4所示。可以看出光纤在不同环境下的温度变化值均不相同。以2号光纤的温升为例,在空气环境中光纤温升幅度最大,超过了13℃；在砂环境中光纤温升幅度较在空气环境中小一点,但也超过了11℃；在纯水环境中光纤温升幅度最小,温度升高了6.53℃；而在各种不同含水率的砂环境中光纤温升幅度相差不大,但存在差异。这些数据表明,在不同的环境中该新型光缆均可反映出环境特征,验证了光缆用于堤坝渗漏监测的可行性。
![]()
图 3 1号光纤在不同环境下温升对比图Figure 3. Comparison of temperature rises of fiber No. 1 under different environments![]()
图 4 2号光纤在不同环境下温升对比图Figure 4. Comparison of temperature rises of fiber No. 2 under different environments3. 光纤渗漏监测现场原型试验
为研究新型分布式光纤渗流监测技术的适用性与仪器现场安装埋设技术,依托引额济乌一期工程开展了渠堤渗漏现场试验研究与验证。
3.1 试验段布置方案
图5为新疆某输水渠道照片。选择在引额济乌总干渠+70断面附近进行现场试验,该处由于铺设横向排水管的需要,对渠坡原有防渗结构进行了拆除,方便现场安装分布式渗漏监测仪器。考虑到渠道的对称型式,只在渠道一侧的渠堤上布设渗漏光缆（布置剖面见图6）,用于测量渠堤防渗膜后的实际渗漏水位高度,为后期研究渠坡变形破坏提供实测数据。
![]()
图 5 非运行期时新疆某输水渠道图片Figure 5. Picture of a water transfer in Xinjiang during non-operational period为验证分布式监测技术的可行性和精度,在试验段的局部位置安装埋设了传统的孔隙水压力计,并加装自动化测读装置,通过物联网将测量数据上传至位于南京水利科学研究院的服务器内。如图7所示,在渠底防渗膜下部布设两只孔隙水压力计,在坡顶以下11 m位置,布设3只孔隙水压力计,孔隙水压力计的安装要求参照《土石坝安全监测技术规范（SL551— 2012）》。
![]()
图 7 渗漏光缆和孔隙水压力计布设位置示意图Figure 7. Schematic diagram of placement of pore water pressure gauge and cable3.2 试验段测温光纤安装埋设技术
现场光缆布设方案见图8。
埋设步骤：①将要铺设的光缆一头熔接跳线,另一头相互熔接成一体；②依图6和7在渠坡土体中开挖宽约15 cm、深约50 cm的光缆沟,光缆沟间隔2 m,依次开挖10条,相互连通；③在光缆沟底部摊铺、压实厚约10 cm的剔除过大块砾石的原岸坡土料；④在光缆沟内铺设好光缆,光缆检查合格后回填光缆沟,小型振动碾压实；⑤光缆的光损检查合格后,铺设防渗膜和衬砌板。
经现场实际踏勘、放样,依据光缆布设方案安装光缆后将光缆引出至渠顶。同沟埋设孔隙水压力计与光缆,以保证测值的对比性,安装时记录仪器实际安装高程,并将电缆线呈S型顺光缆沟引出至堤顶。
采用光功率计测试光缆的工作状态,检查合格后,采用人工结合小型机械,进行渠坡的削坡、整平和压实工作。待摊铺一层厚约10 cm的素砂浆找平层后,铺设两布一膜,采用防水膜焊接机,将新铺设的两布一膜与老膜连成一体,检测合格后,在两布一膜上部铺设六棱块衬砌板。现场施工完毕后,将孔隙水压力计的电缆线汇集到光缆在渠顶的露头位置,并安装由南京水利科学研究院研制的自动化采集设备,开始监测渠道膜后水位随渠道水位的变化情况。
3.3 光纤测点定位分析
根据渠道运行的规律,分别在渠道通水后、渠道通水期间、渠道停水前逐渐降低运行水位时以及渠道停水后等有代表性的时间点内进行测量。为了更好的发挥渗漏监测光缆的优点,在正式测量前开展了点加热法的光纤定位工作,得到了光缆1号纤和2号纤的起始测量位置和终了测量位置。
如图9所示,为了测量方便,测量前将1号纤和2号纤串联成一根光纤,试验过程中通过拉曼光纤解调仪同时测量,保证了测量数据的时效性和可对比性,能有效减少其他试验因素的干扰。从图9中可以明显看出,在光纤距离向22.5,189.5,205.5和370.5这4个点上,光纤测值有明显的变化,说明点加热定位技术在现场条件下能有效定位光纤测点的位置,将光纤测量数据与光纤实际安装位置一一对应。
依据笔者已有的监测经验,为有效消除系统测量误差对数据的干扰,本次试验选择加热系统未工作前6次测值的平均值作为各个测点的基准值,待光纤温升数值逐渐趋于稳定,取稳定后6次测值的平均值作为各个测点加热后的温度值,前后两值相减,即可得到光纤各点的温升值。
3.4 监测结果分析
为更好地通过数值来判断光缆所处位置的渗流情况,对温升曲线进行了拟合,得到相应的量化指标,温升曲线拟合公式如下所示：
ΔT=Kt2, (1) 式中,K为环境特征值,反映温度的增长速率。
根据光纤渗漏监测室内试验的成果,采用式（1）分析光纤每个点的温升曲线,可得各个点的K值,绘制曲线如图10所示。
由于光缆呈S型布置,从曲线形式上可以看出,渠道膜后渗水对光缆各点K值有较明显的影响,即渗水区域的光纤K值均在一定区间波动。依据前期室内试验的研究成果,确定温差3℃的点为膜后水位的临界点判断标准,K值大于0.7为膜后水位的临界点判断标准,现场代表监测数据分析结果如表1所示。
表 1 光纤测渗漏现场测值分析表Table 1. Analysis of measured values of optical fiber leakage日期 光纤在水位以上的长度/m 水面距离孔压的高度/m 孔压相应的测值/m 2019-07-22 温差法 56.5~66 3.25 3.175 85.5~93.5 3.61 116.5~122.5 4.03 144.5~152.5 3.62 2019-07-22 K值法 56~65.5 3.25 84.5~93 3.50 116.5~122.5 4.03 144.5~152 3.73 2019-09-19 温差法 50.5~70 1.02 0.685 82~101.5 1.04 109~130.5 0.57 140~157.5 1.50 2019-09-19 K值法 51~70 1.13 82.5~101 1.26 111~130 1.13 141.5~160.5 1.16 上述的两种分析判断方法,其结果基本一致,与孔隙水压力的测值较为接近。由于本次光纤测量以0.5 m为距离单位进行测量的,可以形象的认为光纤上每隔0.5 m有一个传感器。根据光纤埋设安装图纸可知,当膜后水位位于两个光纤测点之间时,采用光纤数据判断水位的高度,会与孔隙水压力计测值有一定差距,该差值最大值为0.5 m。由于光纤安装呈S型分布,上述对水位的分析判断,往往采用两个边界点,因此该种测量方法与孔隙水压力计测值的最大差值为1 m。光纤测量的平均值与孔压测值相比,误差在0.5 m以内,完全达到了光纤测量精度,说明采用光纤渗漏监测技术监测现场膜后水位,测值准确,能够适应现场测量条件,证明了光纤渗漏监测技术在高寒区长距离输水渠道应用的可行性和可靠性。
4. 结论
本文基于拉曼光时域反射原理研发了分布式光纤拉曼渗流监测技术,提高了光缆测温敏感性和抗冻性,可用于高寒区长距离输水渠道渗漏监测。结论如下：
（1）室内模拟试验结果表明新型光缆在不同环境中均可反映出环境特征,验证了光缆用于堤坝渗漏监测的可行性。
（2）现场采用点加热定位技术能有效定位光纤测点的位置,实现了光纤测量数据与光纤实际安装位置的一一对应。
（3）本次测量光纤采样点间距为0.5 m,光纤测量的平均值与孔隙水压力计测值差小于0.5 m,差值在理论范围内,说明光纤渗漏监测技术监测现场膜后水位测值准确,完全能够适应现场测量条件。
分布式光纤拉曼渗流监测技术克服了传统方法间隔长、可靠性差、施工成本高等缺点；可实现长距离分布式测量、施工简便、适应性强,适合于推广应用。
-
![]()
图 21 不同磨料浓度下底部坑径随冲击时间的变化曲线
Figure 21. Variation curves of diameter of bottom pit with impact time under different abrasive volume concentrations
表 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 
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
-
[1] 侯冰, 武安安, 常智, 等. 页岩油储层多甜点压裂裂缝垂向扩展试验研究[J]. 岩土工程学报, 2021, 43(7): 1322–1330. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202107024.htm HOU Bing, WU An-an, CHANG Zhi, et al. Experimental study on vertical propagation of fractures of multi-sweet of spots shale oil reservoir[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(7): 1322–1330. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202107024.htm
[2] 司鹄, 谢延明, 杨春和. 磨料水射流作用下岩石损伤场的数值模拟[J]. 岩土力学, 2011, 32(3): 935–940. doi: 10.3969/j.issn.1000-7598.2011.03.048 SI Hu, XIE Yan-ming, YANG Chun-he. Numerical simulation of rock damage field under abrasive water jet[J]. Rock and Soil Mechanics, 2011, 32(3): 935–940. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.03.048
[3] ZHANG T T, LIU Z Y, XIE Z X, et al. Numerical simulation of oil shale in situ mining using fluid-thermo-solid coupling[J]. Global Geology, 2020, 23(4): 247–254.
[4] 夏彬伟, 刘浪, 彭子烨, 等. 致密砂岩水平井多裂缝扩展及转向规律研究[J]. 岩土工程学报, 2020, 42(8): 1549–1555. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202008027.htm XIA Bin-wei, LIU Lang, PENG Zi-ye, et al. Multi-fracture propagation and deflection laws of horizontal wells in tight sandstone[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(8): 1549–1555. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202008027.htm
[5] LIU W, BELL A, WANG Z M, et al. Evaluation of the slurry erosion resistance of the body materials of oil & gas drill bits with a modified abrasive waterjet[J]. Wear, 2020, 456/457: 203364. doi: 10.1016/j.wear.2020.203364
[6] QIANG C H, WANG F C, GUO C W. Study on cutting speed and energy utilization rate in processing stainless steel with abrasive water jet[J]. The International Journal of Advanced Manufacturing Technology, 2020, 108(5/6): 1875–1886.
[7] 陈国庆, 简大华, 陈宇航, 等. 不同含水率冻融后红砂岩剪切蠕变特性[J]. 岩土工程学报, 2021, 43(4): 661–669. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202104010.htm CHEN Guo-qing, JIAN Da-hua, CHEN Yu-hang, et al. Shear creep characteristics of red sandstone after freeze-thaw with different water contents[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(4): 661–669. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202104010.htm
[8] 程晓泽, 任福深, 方天成, 等. 粒子射流耦合冲击破岩实验[J]. 石油学报, 2018, 39(2): 232–239. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB201802012.htm CHENG Xiao-ze, REN Fu-shen, FANG Tian-cheng, et al. Experiment of rock breaking by particle-jet coupling impact[J]. Acta Petrolei Sinica, 2018, 39(2): 232–239. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB201802012.htm
[9] MA G W, WANG X J, REN F. Numerical simulation of compressive failure of heterogeneous rock-like materials using SPH method[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(3): 353–363. doi: 10.1016/j.ijrmms.2011.02.001
[10] POZZETTI G, PETERS B. A numerical approach for the evaluation of particle-induced erosion in an abrasive waterjet focusing tube[J]. Powder Technology, 2018, 333: 229–242. doi: 10.1016/j.powtec.2018.04.006
[11] 黄中伟, 张世昆, 李根生, 等. 液氮磨料射流破碎高温花岗岩机理[J]. 石油学报, 2020, 41(5): 604–614. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB202005010.htm HUANG Zhong-wei, ZHANG Shi-kun, LI Gen-sheng, et al. Breakage mechanism of high-temperature granite by abrasive liquid nitrogen jet[J]. Acta Petrolei Sinica, 2020, 41(5): 604–614. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB202005010.htm
[12] 庄欠伟, 袁一翔, 徐天明, 等. 射流联合盾构切削钢筋混凝土仿真与试验[J]. 岩土工程学报, 2020, 42(10): 1817–1824. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202010009.htm ZHUANG Qian-wei, YUAN Yi-xiang, XU Tian-ming, et al. Simulation and experiment on cutting reinforced concrete with jet combined shield method[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(10): 1817–1824. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202010009.htm
[13] 穆朝民, 戎立帆. 磨料射流冲击岩石损伤机制的数值分析[J]. 岩土力学, 2014, 35(5): 1475–1481. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201405038.htm MU Chao-min, RONG Li-fan. Numerical simulation of damage mechanism of abrasive water jet impaction on rock[J]. Rock and Soil Mechanics, 2014, 35(5): 1475–1481. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201405038.htm
[14] 刘勇, 李志飞, 魏建平, 等. 磨料空气射流破煤冲蚀模型研究[J]. 煤炭学报, 2020, 45(5): 1733–1742. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202005019.htm LIU Yong, LI Zhi-fei, WEI Jian-ping, et al. Erosion model of abrasive air jet used in coal breaking[J]. Journal of China Coal Society, 2020, 45(5): 1733–1742. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202005019.htm
[15] 李世杰, 王艾伦, 刘向军, 等. 基于SPH算法土壤水射流冲击演化数值仿真研究[J]. 计算机仿真, 2019, 36(3): 243–247, 384. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJZ201903049.htm LI Shi-jie, WANG Ai-lun, LIU Xiang-jun, et al. Numerical simulation of soil water jet impact evolution based on SPH algorithm[J]. Computer Simulation, 2019, 36(3): 243–247, 384. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSJZ201903049.htm
[16] 赵健, 张贵才, 徐依吉, 等. 基于SPH方法粒子射流破岩数值模拟与实验研究[J]. 爆炸与冲击, 2017, 37(3): 479–486. https://www.cnki.com.cn/Article/CJFDTOTAL-BZCJ201703014.htm ZHAO Jian, ZHANG Gui-cai, XU Yi-ji, et al. SPH-based numerical simulation and experimental study on rock breaking by particle impact[J]. Explosion and Shock Waves, 2017, 37(3): 479–486. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BZCJ201703014.htm
[17] LIU M B, LIU G R. Smoothed particle hydrodynamics (SPH): an overview andRecent developments[J]. Archives of Computational Methods in Engineering, 2010, 17(1): 25–76. doi: 10.1007/s11831-010-9040-7
[18] HEUZÉ O. General form of the Mie-Grüneisen equation of state[J]. Comptes Rendus Mécanique, 2012, 340(10): 679-687. doi: 10.1016/j.crme.2012.10.044
[19] 许传波. 页岩气藏地层特性分析及爆炸压裂适应性研究[D]. 东营: 中国石油大学(华东), 2014. XU Chuan-bo. The Study of Shale Gas Reservoir Characteristics and Adaptability of Explosive Fracturing[D]. Dongying: China University of Petroleum (Huadong), 2014. (in Chinese)
[20] IQBAL M A, KUMAR V, MITTAL A K. Experimental and numerical studies on the drop impact resistance of prestressed concrete plates[J]. International Journal of Impact Engineering, 2019, 123: 98–117. doi: 10.1016/j.ijimpeng.2018.09.013
[21] SILIN D, PATZEK T. Pore space morphology analysis using maximal inscribed spheres[J]. Physica A: Statistical Mechanics and Its Applications, 2006, 371(2): 336–360. doi: 10.1016/j.physa.2006.04.048
[22] 饶登宇, 白冰. 孔隙尺度下三维多孔介质扩散迂曲度的SPH计算[J]. 岩土工程学报, 2020, 42(5): 961–967. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202005024.htm RAO Deng-yu, BAI Bing. Pore-scale SPH simulations of diffusive tortuosity in 3-D porous media[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 961–967. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202005024.htm
[23] RAO D Y, BAI B. Study of the factors influencing diffusive tortuosity based on pore-scale SPH simulation of granular soil[J]. Transport in Porous Media, 2020, 132(2): 333–353. doi: 10.1007/s11242-020-01394-0
[24] BAI B, RAO D Y, XU T, et al. SPH-FDM boundary for the analysis of thermal process in homogeneous media with a discontinuous interface[J]. International Journal of Heat and Mass Transfer, 2018, 117: 517–526. doi: 10.1016/j.ijheatmasstransfer.2017.10.004
[25] 王冬欣. 基于Micro-CT图像的数字岩心孔隙级网络建模研究[D]. 长春: 吉林大学, 2015. WANG Dong-xin. The Research of Digital Core Network Extraction Based on Micro-CT Images[D]. Changchun: Jilin University, 2015. (in Chinese)
[26] 张鹏伟, 胡黎明, Jay N Meegoda, 等. 基于岩土介质三维孔隙结构的两相流模型[J]. 岩土工程学报, 2020, 42(1): 37–45. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202001008.htm ZHANG Peng-wei, HU Li-ming, MEEGODA J N, et al. Two-phase flow model based on 3D pore structure of geomaterials[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 37–45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202001008.htm
[27] 王金波. 岩石孔隙结构三维重构及微细观渗流的数值模拟研究[D]. 北京: 中国矿业大学(北京), 2014. WANG Jin-bo. 3d Reconstruction of Porous Rock and Numerical Simulations of Fluid Flow at Mesoscale Levels[D]. Beijing: China University of Mining & Technology, Beijing, 2014. (in Chinese)
[28] 李井慧. 磨料水射流切割大理石的试验及仿真研究[D]. 哈尔滨: 哈尔滨理工大学, 2015. LI Jing-hui. Research on Simulation and Experiment of AWJ Cutting Marble[D]. Harbin: Harbin University of Science and Technology, 2015. (in Chinese)
[29] 林晓东, 卢义玉, 汤积仁, 等. 基于SPH-FEM耦合算法的磨料水射流破岩数值模拟[J]. 振动与冲击, 2014, 33(18): 170–176. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ202116017.htm LIN Xiao-dong, LU Yi-yu, TANG Ji-ren, et al. Numerical simulation of abrasive water jet breaking rock with SPH-FEM coupling algorithm[J]. Journal of Vibration and Shock, 2014, 33(18): 170–176. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ202116017.htm
[30] 原油和石油产品密度在638 kg/m3到1074 kg/m3范围内的烃压缩系数: GB/T 21450—2008[S]. 2008. Crude Petroleum and Petroleum Products-Compressibility Factors for Hydrocarbons in the Range 638 kg/m3 to 1074 kg/m3: GB/T 21450—2008[S]. 2008. (in Chinese)
-
期刊类型引用(2)
1. 张子洋，汪波，刘锦超，刘金炜，杨凯. 后注浆预应力树脂锚杆应力损失后注浆段效用分析. 铁道科学与工程学报. 2025(02): 783-794 . 
百度学术
2. 谈博海，姚囝，杜键，孙明伟，刘海，关文超. 膨胀型浆体注浆锚杆拉拔力学特性及损伤失效机制研究. 岩石力学与工程学报. 2024(12): 3058-3069 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 335
- HTML全文浏览量: 33
- PDF下载量: 113
- 被引次数: 2
