CN116401966A - 一种多孔裂隙岩体时变浆液扩散过程的模拟方法 - Google Patents
一种多孔裂隙岩体时变浆液扩散过程的模拟方法 Download PDFInfo
- Publication number
- CN116401966A CN116401966A CN202310208388.3A CN202310208388A CN116401966A CN 116401966 A CN116401966 A CN 116401966A CN 202310208388 A CN202310208388 A CN 202310208388A CN 116401966 A CN116401966 A CN 116401966A
- Authority
- CN
- China
- Prior art keywords
- slurry
- water
- fracture
- matrix
- rock mass
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000002002 slurry Substances 0.000 title claims abstract description 95
- 239000011435 rock Substances 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000009792 diffusion process Methods 0.000 title claims abstract description 23
- 238000004088 simulation Methods 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 44
- 239000011159 matrix material Substances 0.000 claims description 39
- 239000012530 fluid Substances 0.000 claims description 17
- 230000035699 permeability Effects 0.000 claims description 12
- 239000003795 chemical substances by application Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 230000036571 hydration Effects 0.000 claims description 5
- 238000006703 hydration reaction Methods 0.000 claims description 5
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000005429 filling process Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 239000000758 substrate Substances 0.000 claims description 3
- 238000003756 stirring Methods 0.000 claims description 2
- 206010057175 Mass conditions Diseases 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 15
- 239000006185 dispersion Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000004593 Epoxy Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001802 infusion Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
本发明公开了一种多孔裂隙岩体时变浆液扩散过程的模拟方法，该方法包括步骤一：使用数值仿真软件建立多孔介质模型；步骤二：在数值仿真软件中设置多孔介质模型相关方程组和关键参数，方程组包括质量守恒、力学本构、水力本构方程，关键参数包括岩体参数和浆液参数；步骤三：在数值仿真软件中构建网格；步骤四：使用数值仿真软件自动进行扩散计算并生成结果文件。本申请在测定岩体参数和浆液参数后，利用数值仿真软件对方程组进行联立求解，即可获得模拟结果，模拟结果可以反映不同注浆压力下时变性浆液在多孔介质中扩散过程的主要特征，进而可将此方法用于不同岩体条件下的注浆方案比选和优化，并最终为实现安全高效的岩体注浆提供技术指导。
Description
技术领域
本发明属于地下工程技术领域，具体涉及一种多孔裂隙岩体时变浆液扩散过程的模拟方法。
背景技术
注浆作为处理微裂隙渗漏水问题的有效技术手段，在地下工程中得到了广泛应用与极大发展。然而由于涉及到化学、流体力学、工程地质学、岩土力学以及机械技术等多门学科的注浆技术极其复杂，因此注浆理论至今仍远落后于工程实践。
目前，岩体灌浆的数值模型大多未考虑浆液时变性与岩体多孔特性的影响效应。值得注意的是，由于注浆所用水泥浆液的流变参数在灌注过程中并非保持不变的，同时裂隙的存在也将岩体介质分为了岩隙(裂隙体)和岩基(孔隙体)两个组成部分，因此现有的数值模型不能准确地描述时变性流体在多孔介质中的运移过程。针对上述问题，建立一个时变浆液灌注过程的多孔介质模型以描述多孔裂隙岩体中时变性浆液的扩散过程是一种可行的方法。
发明内容
为解决上述背景技术中提出的问题，本发明提供一种多孔裂隙岩体时变浆液扩散过程的模拟方法，以解决现有岩体灌浆的数值模型大多未考虑浆液时变性与岩体多孔特性的影响效应的问题。
为实现上述目的，本发明提供如下技术方案：
一种多孔裂隙岩体时变浆液扩散过程的模拟方法，该方法包括：
步骤一：使用数值仿真软件建立多孔介质模型；
步骤二：在数值仿真软件中设置多孔介质模型相关方程组和关键参数，方程组包括质量守恒、力学本构、水力本构方程，关键参数包括岩体参数和浆液参数；
步骤三：在数值仿真软件中构建网格；
步骤四：使用数值仿真软件自动进行扩散计算并生成结果文件。
优选地，使用的数值仿真软件为COMSOL Multiphysics。
优选地，步骤二中的多孔介质模型相关方程组具体表达为：
裂隙中水与浆液的质量守恒方程表达为：
其中ρfw和ρff是裂隙中水和裂隙的密度，φf是裂隙的孔隙度，
和/>
是裂隙中水和浆液的达西流速，Sfw和Sff分别是裂隙中水和浆液的饱和度，εfv是裂隙的体积应变，pfw和pff分别是裂隙中水和浆液的压力，χw和χf分别是水和浆液的压缩系数，/>
是存储系数，Se是有效饱和度，pc＝pff-pfw是毛细压力，гfw和гff分别是水和浆液的交换项，上标小点表示时间导数算子；
多孔基质中水和浆液的质量守恒方程表达为：
其中ρpw和ρpf分别是基质中水和浆液的密度，αp是基质的等效毕奥有效应力系数，φp是基质的孔隙度，εpv是基质的体积应变，ppw和ppf分别是基质中水和浆液的压力，Ks是固体颗粒的体积模量，
和/>
是基质中水和浆液的达西流速；
交换项гfw和гff表达为：
其中
是裂隙的平均固有渗透率，krw和krf分别是水和浆液的相对渗透率，ρw和ρf分别是水和浆液的密度，μw和μf分别是浆液的动力黏度，/>
是交换系数；
交换系数定义为：
其中β是骨料形状的几何系数，γw是一个与骨料几何形状和初始流体压力无关的经验系数，a是裂隙的平均间距；
多孔基质的有效应力σ’p表达为：
σ’p＝σ+αpppw
其中σ是总应力；
基质的应变εp表达为：
εp＝Sp:σ'p
其中Sp是基质的柔度张量；
裂隙的有效应力σ’f表达为：
其中
是裂隙中的平均水压力；
裂隙的应变εf表达为：
εf＝Sf:σ'f
其中Sf是裂隙的柔度张量；
多孔介质的总应变ε表达为：
ε＝εp+εf；
裂隙的法向刚度表达为：
其中Δun是i方向上力学孔径的变化，Kni是初始裂隙刚度，Vm是最大裂隙开度；
流体在流场中的平流运动受一般达西定律的支配,公式表达为：
其中kf和kp分别是裂隙和基质的固有渗透性张量，krπ和μπ分别是流体π的相对渗透率和动力黏度，pfπ和ρfπ分别是裂隙中流体π的压力与密度，ppπ和ρpπ分别是基质中流体π的压力与密度；
浆液的水化时间t主要为：
对主剂和固化剂以及促进剂成一定比例的浆材进行室内测试，进而得到在拌制后一个时间段内该浆材的切力黏度曲线并拟合为浆液水化时间t的公式；
以上方程组构成了时变浆液灌注过程的多孔介质模型，在测定岩体参数和浆液参数后，利用COMSOL Multiphysics对多孔介质模型相关方程组进行联立同时求解，即可获得不同注浆压力下浆液运移过程的演化规律。
优选地，步骤三中构建网格时设置序列类型为物理场控制网格，以及设置超细化的单元大小。
与现有技术相比，本发明的有益效果是：
本申请提出一种多孔裂隙岩体时变浆液扩散过程的模拟方法，在测定岩体参数和浆液参数后，利用数值仿真软件对设置的方程组进行联立求解，即可获得模拟结果，模拟结果可以反映不同注浆压力下时变性浆液在多孔介质中扩散过程的主要特征，进而可将此方法用于不同岩体条件下的注浆方案比选和优化，并最终为实现安全高效的岩体注浆提供技术指导。
附图说明
图1为一种多孔裂隙岩体时变浆液扩散过程的模拟方法的步骤示意图。
图2为数值仿真软件的参数取值界面示意图；
图3为数值仿真软件的方程设置界面示意图；
图4为数值仿真软件的网格设置界面示意图；
图5为数值仿真软件的结果处理界面示意图；
图6为多孔介质模型结构示意图；
图7为t＝30s、注浆压力为1Mpa时的裂隙示意图；
图8为t＝30s、注浆压力为1Mpa时的基质示意图；
图9为t＝60s、注浆压力为1Mpa时的裂隙示意图；
图10为t＝60s、注浆压力为1Mpa时的基质示意图；
图11为t＝90s、注浆压力为1Mpa时的裂隙示意图；
图12为t＝90s、注浆压力为1Mpa时的基质示意图；
图13为t＝90s、注浆压力为2Mpa时的裂隙示意图；
图14为t＝90s、注浆压力为2Mpa时的基质示意图；
图15为t＝90s、注浆压力为3Mpa时的裂隙示意图；
图16为t＝90s、注浆压力为3Mpa时的基质示意图；
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
如图1所示，一种多孔裂隙岩体时变浆液扩散过程的模拟方法，该方法包括：
步骤一：使用数值仿真软件建立多孔介质模型；
步骤二：在数值仿真软件中设置多孔介质模型相关方程组和关键参数，方程组包括质量守恒、力学本构、水力本构方程，关键参数包括岩体参数和浆液参数；
步骤三：在数值仿真软件中构建网格；
步骤四：使用数值仿真软件自动进行扩散计算并生成结果文件。
在本实施例中，本申请提出一种多孔裂隙岩体时变浆液扩散过程的模拟方法，在测定岩体参数和浆液参数后，利用数值仿真软件对设置的方程组进行联立求解，即可获得模拟结果，模拟结果可以反映不同注浆压力下时变性浆液在多孔介质中扩散过程的主要特征，进而可将此方法用于不同岩体条件下的注浆方案比选和优化，并最终为实现安全高效的岩体注浆提供技术指导。
实施例2
本实施例与实施例1的区别在于：使用的数值仿真软件为COMSOL Multiphysics的5.5版本。
实施例3
本实施例与实施例2的区别在于：如图3所示，步骤二中的多孔介质模型相关方程组具体表达为：
裂隙中水与浆液的质量守恒方程表达为：
其中ρfw和ρff是裂隙中水和裂隙的密度，φf是裂隙的孔隙度，
和/>
是裂隙中水和浆液的达西流速，Sfw和Sff分别是裂隙中水和浆液的饱和度，εfv是裂隙的体积应变，pfw和pff分别是裂隙中水和浆液的压力，χw和χf分别是水和浆液的压缩系数，/>
是存储系数，Se是有效饱和度，pc＝pff-pfw是毛细压力，гfw和гff分别是水和浆液的交换项，上标小点表示时间导数算子；
多孔基质中水和浆液的质量守恒方程表达为：
其中ρpw和ρpf分别是基质中水和浆液的密度，αp是基质的等效毕奥有效应力系数，φp是基质的孔隙度，εpv是基质的体积应变，ppw和ppf分别是基质中水和浆液的压力，Ks是固体颗粒的体积模量，
和/>
是基质中水和浆液的达西流速；
交换项гfw和гff表达为：
其中
是裂隙的平均固有渗透率，krw和krf分别是水和浆液的相对渗透率，ρw和ρf分别是水和浆液的密度，μw和μf分别是浆液的动力黏度，/>
是交换系数；
交换系数定义为：
其中β是骨料形状的几何系数，γw是一个与骨料几何形状和初始流体压力无关的经验系数，a是裂隙的平均间距；
多孔基质的有效应力σ’p表达为：
σ’p＝σ+αpppw
其中σ是总应力；
基质的应变εp表达为：
εp＝Sp:σ'p
其中Sp是基质的柔度张量；
裂隙的有效应力σ’f表达为：
其中
是裂隙中的平均水压力；
裂隙的应变εf表达为：
εf＝Sf:σ'f
其中Sf是裂隙的柔度张量；
多孔介质的总应变ε表达为：
ε＝εp+εf；
裂隙的法向刚度表达为：
其中Δun是i方向上力学孔径的变化，Kni是初始裂隙刚度，Vm是最大裂隙开度；
流体在流场中的平流运动受一般达西定律的支配,公式表达为：
其中kf和kp分别是裂隙和基质的固有渗透性张量，krπ和μπ分别是流体π的相对渗透率和动力黏度，pfπ和ρfπ分别是裂隙中流体π的压力与密度，ppπ和ρpπ分别是基质中流体π的压力与密度；
对主剂∶固化剂∶促进剂＝100∶6.0∶1.5的JD-1改性环氧浆材进行室内测试，进而得到浆液在拌制后6h内的切力黏度曲线并拟合成如下公式：
η(t)＝1×10-20t5-4×10-16t4+4×10-12t3-2×10-8t2+3×10-5t+0.0025
其中t是浆液的水化时间；
以上方程组通过质量守恒、力学本构、水力本构方程建立时变浆液灌注过程的多孔介质模型，在测定岩体参数和浆液参数后，利用COMSOL Multiphysics对多孔介质模型相关方程组进行联立同时求解，即可获得不同注浆压力下浆液运移过程的演化规律。
实施例4
本实施例与实施例1的区别在于：如图4所示，步骤三中构建网格时设置序列类型为物理场控制网格，以及设置超细化的单元大小。
实施例5
如图6所示，多孔介质模型包括6864个域单元和228个边界单元。
实施例6
如图2、图5所示，岩石力学参数及浆液流变参数的确定是保证数值计算准确性及合理性的关键。通常来说，室内试验获得的岩石力学参数不能直接运用于岩体工程数值计算，因为试样为含有少量微裂隙的岩块，而实际工程岩体节理裂隙较多且受力环境更加复杂，其力学性质差异较大。室内试验结果表明岩石力学参数受到加载应力及水压的影响，再加上岩石试样的尺寸效应问题，故室内试验得到的岩石参数需根据现场工程地质条件及其他类似工程数值计算数据进行适当修正。本次岩体灌浆数值模拟计算参数如表1所示，
表1数值模型计算参数
注:实际参数与工程案例有关，参数的标定可通过敏感性分析及物理实验等方法进行测量，不同工程应用中参数的取值具有较大差异。
如图7、图8所示，图为当t＝30s，注浆压力为1Mpa时浆液的扩散范围；
如图9、图10所示，图为当t＝60s，注浆压力为1Mpa时浆液的扩散范围；
如图11、图12所示，图为当t＝90s，注浆压力为1Mpa时浆液的扩散范围；
如图13、图14所示，图为当t＝90s，注浆压力为2Mpa时浆液的扩散范围；
如图15、图16所示，图为当t＝90s，注浆压力为3Mpa时浆液的扩散范围。
在本实施例中，模拟结果可以反映不同注浆压力下时变性浆液在多孔介质中扩散过程的主要特征，进而可将此方法用于不同岩体条件下的注浆方案比选和优化，并最终为实现安全高效的岩体注浆提供技术指导。
Claims (4)
1.一种多孔裂隙岩体时变浆液扩散过程的模拟方法，其特征在于，该方法包括：
步骤一：使用数值仿真软件建立多孔介质模型；
步骤二：在数值仿真软件中设置多孔介质模型相关方程组和关键参数，方程组包括质量守恒、力学本构、水力本构方程，关键参数包括岩体参数和浆液参数；
步骤三：在数值仿真软件中构建网格；
步骤四：使用数值仿真软件自动进行扩散计算并生成结果文件。
2.根据权利要求1所述的一种多孔裂隙岩体时变浆液扩散过程的模拟方法，其特征在于，使用的数值仿真软件为COMSOL Multiphysics。
3.根据权利要求2所述的一种多孔裂隙岩体时变浆液扩散过程的模拟方法，其特征在于，步骤二中的多孔介质模型相关方程组具体表达为：
裂隙中水与浆液的质量守恒方程表达为：
其中ρfw和ρff是裂隙中水和裂隙的密度，φf是裂隙的孔隙度，
和/>
是裂隙中水和浆液的达西流速，Sfw和Sff分别是裂隙中水和浆液的饱和度，εfv是裂隙的体积应变，pfw和pff分别是裂隙中水和浆液的压力，χw和χf分别是水和浆液的压缩系数，/>
是存储系数，Se是有效饱和度，pc＝pff-pfw是毛细压力，гfw和гff分别是水和浆液的交换项，上标小点表示时间导数算子；
多孔基质中水和浆液的质量守恒方程表达为：
其中ρpw和ρpf分别是基质中水和浆液的密度，αp是基质的等效毕奥有效应力系数，φp是基质的孔隙度，εpv是基质的体积应变，ppw和ppf分别是基质中水和浆液的压力，Ks是固体颗粒的体积模量，
和/>
是基质中水和浆液的达西流速；
交换项гfw和гff表达为：
其中
是裂隙的平均固有渗透率，krw和krf分别是水和浆液的相对渗透率，ρw和ρf分别是水和浆液的密度，μw和μf分别是浆液的动力黏度，/>
是交换系数；
交换系数定义为：
其中β是骨料形状的几何系数，γw是一个与骨料几何形状和初始流体压力无关的经验系数，a是裂隙的平均间距；
多孔基质的有效应力σ’p表达为：
σ′p＝σ+αpppw
其中σ是总应力；
基质的应变εp表达为：
εp＝Sp:σ'p
其中Sp是基质的柔度张量；
裂隙的有效应力σ’f表达为：
其中
是裂隙中的平均水压力；
裂隙的应变εf表达为：
εf＝Sf:σ'f
其中Sf是裂隙的柔度张量；
多孔介质的总应变ε表达为：
ε＝εp+εf；
裂隙的法向刚度表达为：
其中Δun是i方向上力学孔径的变化，Kni是初始裂隙刚度，Vm是最大裂隙开度；
流体在流场中的平流运动受一般达西定律的支配,公式表达为：
其中kf和kp分别是裂隙和基质的固有渗透性张量，krπ和μπ分别是流体π的相对渗透率和动力黏度，pfπ和ρfπ分别是裂隙中流体π的压力与密度，ppπ和ρpπ分别是基质中流体π的压力与密度；
浆液的水化时间t主要为：
对主剂和固化剂以及促进剂成一定比例的浆材进行室内测试，进而得到在拌制后一个时间段内该浆材的切力黏度曲线并拟合为浆液水化时间t的公式；
以上方程组构成了时变浆液灌注过程的多孔介质模型，在测定岩体参数和浆液参数后，利用COMSOL Multiphysics对多孔介质模型相关方程组进行联立同时求解，即可获得不同注浆压力下浆液运移过程的演化规律。
4.根据权利要求1所述的一种多孔裂隙岩体时变浆液扩散过程的模拟方法，其特征在于，步骤三中构建网格时设置序列类型为物理场控制网格，以及设置超细化的单元大小。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310208388.3A CN116401966A (zh) | 2023-03-07 | 2023-03-07 | 一种多孔裂隙岩体时变浆液扩散过程的模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310208388.3A CN116401966A (zh) | 2023-03-07 | 2023-03-07 | 一种多孔裂隙岩体时变浆液扩散过程的模拟方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116401966A true CN116401966A (zh) | 2023-07-07 |
Family
ID=87013258
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310208388.3A Pending CN116401966A (zh) | 2023-03-07 | 2023-03-07 | 一种多孔裂隙岩体时变浆液扩散过程的模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116401966A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117252033A (zh) * | 2023-11-13 | 2023-12-19 | 山东科技大学 | 一种考虑热流固耦合的粗糙裂隙浆液扩散过程描述方法 |
| CN117408191A (zh) * | 2023-12-15 | 2024-01-16 | 山东大学 | 一种裂隙岩体渗透-流动注浆模拟方法及系统 |
| CN117421934A (zh) * | 2023-12-15 | 2024-01-19 | 山东大学 | 一种考虑温度效应的多孔介质注浆模拟方法及系统 |
-
2023
- 2023-03-07 CN CN202310208388.3A patent/CN116401966A/zh active Pending
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117252033A (zh) * | 2023-11-13 | 2023-12-19 | 山东科技大学 | 一种考虑热流固耦合的粗糙裂隙浆液扩散过程描述方法 |
| CN117252033B (zh) * | 2023-11-13 | 2024-03-15 | 山东科技大学 | 一种考虑热流固耦合的粗糙裂隙浆液扩散过程描述方法 |
| CN117408191A (zh) * | 2023-12-15 | 2024-01-16 | 山东大学 | 一种裂隙岩体渗透-流动注浆模拟方法及系统 |
| CN117421934A (zh) * | 2023-12-15 | 2024-01-19 | 山东大学 | 一种考虑温度效应的多孔介质注浆模拟方法及系统 |
| CN117421934B (zh) * | 2023-12-15 | 2024-03-15 | 山东大学 | 一种考虑温度效应的多孔介质注浆模拟方法及系统 |
| CN117408191B (zh) * | 2023-12-15 | 2024-04-02 | 山东大学 | 一种裂隙岩体渗透-流动注浆模拟方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116401966A (zh) | 一种多孔裂隙岩体时变浆液扩散过程的模拟方法 | |
| CN106599449B (zh) | 一种溶洞体积计算的试井解释方法 | |
| CN110296892B (zh) | 基于能量分析的岩石损伤演化过程中特征应力确定方法 | |
| Wu et al. | Long‐term strength determination and nonlinear creep damage constitutive model of salt rock based on multistage creep test: implications for underground natural gas storage in salt cavern | |
| Nelson | An experimental study of fracture permeability in porous rock | |
| CN111257202A (zh) | 一种含吸附气条件下页岩压裂液强制渗吸及返排实验方法 | |
| Bohnsack et al. | Stress sensitivity of porosity and permeability under varying hydrostatic stress conditions for different carbonate rock types of the geothermal Malm reservoir in Southern Germany | |
| CN104101564A (zh) | 一种非稳态高温高压测试低渗透岩心启动压力梯度的装置及方法 | |
| Al-Ibadi et al. | Experimental investigation and correlation of treatment in weak and high-permeability formations by use of gel particles | |
| Peng et al. | Quantifying the influence of fractures for more-accurate laboratory measurement of shale matrix permeability using a modified gas-expansion method | |
| CN105424466A (zh) | 一种评价波浪作用下砂土对埋置海底管线抗力的方法 | |
| CN104374684A (zh) | 用于测试采动过程卸荷煤岩体渗透率的系统及其应用 | |
| Chen et al. | Simulation of damage–permeability coupling for mortar under dynamic loads | |
| Abdulhadi et al. | Experimental study of wellbore instability in clays | |
| Kong et al. | The variation of grain size distribution in rock granular material in seepage process considering the mechanical–hydrological–chemical coupling effect: An experimental research | |
| Yan et al. | Experimental study of hydraulic fracturing for unconsolidated reservoirs | |
| Wang et al. | Investigation of the effect of cemented fractures on fracturing network propagation in model block with discrete orthogonal fractures | |
| CN204165873U (zh) | 用于测试采动过程卸荷煤岩体渗透率的系统 | |
| Niu et al. | Hydromechanical behaviour of columnar jointed rock masses under true triaxial conditions: An experimental and theoretical investigation | |
| CN104675396B (zh) | 一种层状硬脆性泥页岩水化特性的评价装置 | |
| Brown | Recent progress in HDR reservoir engineering | |
| CN112001095B (zh) | 一种建立固井质量评价指标的方法及固井质量评价方法 | |
| Pan et al. | Re-evaluation of the porosity measurements under different confining pressures: A better appraisal of reservoir porosity | |
| CN209148314U (zh) | 一种应用3d技术全真模拟地层评价装置 | |
| Gao et al. | Application of dilation-recompaction model in fracturing optimisation in tight oil reservoir |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |