CN105651964B - 一种确定裂隙岩体表征单元体积的方法 - Google Patents

一种确定裂隙岩体表征单元体积的方法 Download PDF

Info

Publication number
CN105651964B
CN105651964B CN201511022607.0A CN201511022607A CN105651964B CN 105651964 B CN105651964 B CN 105651964B CN 201511022607 A CN201511022607 A CN 201511022607A CN 105651964 B CN105651964 B CN 105651964B
Authority
CN
China
Prior art keywords
crack rock
rock
volume
actual crack
actual
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.)
Active
Application number
CN201511022607.0A
Other languages
English (en)
Other versions
CN105651964A (zh
Inventor
巢志明
王环玲
徐卫亚
吉华
闫龙
向志鹏
王苏生
杨兰兰
赵恺
夏季
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201511022607.0A priority Critical patent/CN105651964B/zh
Publication of CN105651964A publication Critical patent/CN105651964A/zh
Application granted granted Critical
Publication of CN105651964B publication Critical patent/CN105651964B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Medicinal Chemistry (AREA)
  • Food Science & Technology (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • Dispersion Chemistry (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种确定裂隙岩体表征单元体积的方法,包括如下步骤:选取任意尺寸的实际裂隙岩体,绘制其三维数字模型;旋转该三维数字模型,使其与水平方向成不同倾角,通过3D打印制得该实际裂隙岩体不同倾角的相似材料;测定不同倾角相似材料的渗透率Ka,建立坐标轴,以不同倾角方向的值为点到原点的距离,做出一个封闭曲线;根据该封闭曲线与椭圆的相似性,定量分析裂隙岩体不同方向的渗透系数是否符合渗透张量的规律;根据定量分析结果,扩大或缩小选取的实际裂隙岩体的体积,重复上述步骤,确定该实际裂隙岩体的表征单元体积。本发明可充分模拟实际裂隙岩体,得出的REV值与实际裂隙岩体的REV差异很小,与数值方法相比,确定的岩体REV更精确。

Description

一种确定裂隙岩体表征单元体积的方法
技术领域
本发明涉及一种确定裂隙岩体表征单元体积的方法,特别涉及一种基于3D打印技术,通过室内模型试验确定裂隙岩体表征单元体积的方法,适用于不同类型含有复杂裂隙岩体REV的确定,属于岩土工程技术领域。
背景技术
自然岩体普遍存在着裂隙、节理、裙皱、断层等结构面,使得岩体的力学性质十分复杂,由于岩体的结构面对岩体的力学性质起着决定性的作用,因此岩体存在着尺寸效应,即随着岩体体积的变化,岩体的力学参数(强度、弹性模型、摩擦角、粘聚力和渗透系数等)会产生相应变化,但当岩体体积超过一定值时,岩体的力学参数基本保持稳定,不再随体积变化而变化,这个值称为岩体的REV(表征单元体积),即保持岩体的力学参数基本稳定时的岩体最小体积。
目前在岩石力学领域的研究中,岩体的力学模型主要采用等效连续介质模型和离散介质模型。在实际问题分析中,如果所研究岩体的REV存在,则可以将岩体视为等效连续介质,釆用等效介质模型进行计算;反之,如果岩体的REV不存在或者体积非常大,则说明岩体中裂隙离散程度较高,不适合采用等效连续介质模型,需要采用离散介质模型进行分析,因此,岩体REV是选取岩体力学介质类型、确定分析模型的基础,对确定岩体的力学参数至关重要,是岩石力学领域基本科学问题之一。由于裂隙岩体的REV一般较大,由于受现场条件、试验设备及试验费用的限制,在现场进行原位试验测量裂隙岩体的REV是十分困难甚至是不可能的,因此目前学者多采用解析法和数值模型的方法确定裂隙岩体的REV:
(1)文献“李亚军,姚军,黄朝琴,等.裂缝性油藏等效渗透率张量计算及表征单元体积研究[J].水动力学研究与进展A辑,2010,1:004.”利用等效渗透张量来确定岩体表征单元体积。
(2)文献“周创兵,於三大.论岩体表征单元体积REV—岩体力学参数取值的一个基本问题[J].工程地质学报,1999,7(4):332-336.”利用数学推导的方法体出了计算简单裂隙岩体单元体积的公式。
(3)文献“吴顺川,高艳华,高永涛,等.等效节理岩体表征单元体研究[J].中国矿业大学学报,2014,6:025.”采用有限元数值模拟方法确定岩体的REV。
(4)文献“宁宇,徐卫亚,郑文棠,等.柱状节理岩体随机模拟及其表征单元体尺度研究[J].岩石力学与工程学报,2008,27(6):1202-1208.”采用数值试验的方法确定岩体的REV。
(5)公开号为103197042A的中国发明专利“一种节理岩体表征单元体积的取值方法”采用数值试验的方法确定岩体的REV。
但是与通过试验确定裂隙岩体REV的方法相比,由解析法和数值模拟的方法确定裂隙岩体的REV难免会采用很多简化的假设,所得出的岩体模型与实际裂隙岩体出入较大,不能充分考虑实际现场裂隙岩体的特征,所计算出的裂隙岩体REV与实际裂隙岩体的REV差异较大,因此迫切需要一种比较有效的确定裂隙岩体REV的方法。
发明内容
发明目的:本发明的目的是提供一种确定裂隙岩体表征单元体积的方法,通过本方法能够精确确定不同类型含有复杂裂隙岩体的表征单元体积。
技术方案:
本发明涉及一种确定裂隙岩体表征单元体积的方法,包括如下步骤:
(1)选取任意体积的实际裂隙岩体,绘制其三维数字模型;
(2)旋转该三维数字模型,使其与水平方向成不同倾角,通过3D打印制得该实际裂隙岩体不同倾角的相似材料;
(3)测定不同倾角相似材料的渗透率Ka,建立坐标轴,以不同倾角方向的值为点到原点的距离,做出一个封闭曲线;
(4)根据该封闭曲线与椭圆的近似度,定量分析裂隙岩体不同方向的渗透系数是否符合渗透张量的规律;
(5)基于渗透张量性质,根据定量分析结果,扩大或缩小选取的实际裂隙岩体的体积,重复步骤(1)~(4),确定所测实际裂隙岩体的表征单元体积。
其中,步骤(4)中,当封闭曲线与椭圆的近似度大于90%时,说明裂隙岩体不同方向的渗透系数符合渗透张量的规律。
具体的,步骤(5)中,当定量分析结果为裂隙岩体不同方向的渗透系数符合渗透张量的规律时,缩小选取的实际裂隙岩体的体积,重复步骤(1)~(4),直至裂隙岩体不同方向的渗透系数不符合渗透张量的规律,最后两次选取的实际裂隙岩体体积的平均值即为所测实际裂隙岩体的表征单元体积。
而当定量分析结果为裂隙岩体不同方向的渗透系数不符合渗透张量的规律时,扩大选取的实际裂隙岩体的体积,重复步骤(1)~(4),直至裂隙岩体不同方向的渗透系数符合渗透张量的规律,最后两次选取的实际裂隙岩体体积的平均值即为所测实际裂隙岩体的表征单元体积。此时,当选取的实际裂隙岩体的体积扩大到不合理尺寸时,裂隙岩体不同方向的渗透系数仍然不符合渗透张量的规律,说明该实际裂隙岩体的表征单元体积不存在。
较优的,步骤(4)中,通过对封闭曲线进行函数拟合来确定其与椭圆的近似度。
上述确定封闭曲线与椭圆近似度的方法包括:通过迭代法算出所述封闭曲线对应的拟合度最高的椭圆方程,然后利用该椭圆方程对曲线进行拟合,得到的拟合度即为近似度。
具体的,步骤(3)中,确定第一象限内不同倾角方向上到原点距离为的点的位置,连接相应点形成开放曲线,根据对称性绘制其余象限内的开放曲线,相邻开放曲线在坐标轴上相交形成封闭曲线。确定柱状节理裂隙岩体的REV时,优选使用此方式绘制封闭曲线。
步骤(3)中,也可以通过确定第一象限和第二象限内不同倾角方向的值为点到原点的距离做出开放曲线,根据对称性绘制其余象限内的开放曲线,相邻开放曲线在坐标轴上相交形成封闭曲线。确定一般裂隙岩体的REV时,优选使用此方式绘制封闭曲线。
进一步的,步骤(2)中,采用所测实际裂隙岩体粉末为原材料进行3D打印。
更进一步的,模拟所测实际裂隙岩体的力学性能,配制3D打印的原材料。其相较于使用实际裂隙岩体粉末为原材料进行3D打印更为简单易行。
发明原理:根据渗透张量的性质,将以裂隙岩体各方向渗透系数倒数的开平方为点到原点的距离绘成曲线,如果曲线接近于椭圆,则该裂隙岩体可以采用等效连续介质模型分析其渗流性质,该裂隙岩体的体积可认为等于或大于表征单元体积。
有益效果:与现有技术相比,本发明的优点在于:
(1)本发明通过制备相似材料,并根据试验测得不同柱体倾角相似材料的渗透率,做出相似材料各个方向渗透率开平方倒数的曲线,并分析该曲线与椭圆的近似度,根据分析结果,相应放大或缩小相似材料所模拟的裂隙岩体的体积,重复判断,可以确定裂隙岩体的表征单元体积,本发明首次将该渗透张量性质用于确定裂隙岩体表征单元体积,开创了确定裂隙岩体表征单元体积REV的新方法;
(2)本方法充分考虑实际现场岩体的特征,首次将3D打印技术应用于裂隙岩体REV的确定,通过3D打印可以精确控制裂隙岩体的空间尺寸和岩体中裂隙的类型、大小、空间形态和表面粗糙度,打印出的裂隙岩体相似材料与实际裂隙岩体符合性较高,所制备的裂隙岩体相似材料与实际裂隙岩体极为接近,与通过解析法和数值模拟的方法确定岩体REV相比,计算出的岩体REV与实际裂隙岩体的REV差异很小;
(3)采用本发明的方法,不必进行原位试验,通过室内模型试验即可确定裂隙岩体的REV,在测得的REV精度高的同时,节省了大量物力和财力,实现简单、快捷、高效地确定裂隙岩体的REV。
附图说明
图1为本发明的确定裂隙岩体表征单元体积方法的流程图;
图2为本发明实施例中作出的封闭曲线的示意图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1,本发明的一种确定裂隙岩体表征单元体积的方法,其特征在于,包括如下步骤:
(1)选取任意体积的实际裂隙岩体,绘制其三维数字模型;
(2)旋转该三维数字模型,使其与水平方向成不同倾角,通过3D打印制得该实际裂隙岩体不同倾角的相似材料;
具体的,可以将该三维数字模型绕其中心点旋转成与水平方向成0°,15°,30°,45°,60°,75°,90°的倾斜角,然后分别通过3D打印得到对应倾角的相似材料;打印相似材料时,可以所测实际裂隙岩体为原材料进行3D打印,如将要确定REV的裂隙岩体制成粉末状作为打印岩块的材料,将其中裂隙的充填物制成粉末状作为打印岩体裂隙的材料;也可以模拟所测实际裂隙岩体的力学性能,配制水泥砂浆等原材料用于3D打印。
(3)测定不同倾角相似材料的渗透率Ka,建立坐标轴,以不同倾角方向的值为点到原点的距离,做出一个封闭曲线;
通过测定不同倾角相似材料的渗透率Ka,建立坐标轴,确定不同倾角方向上到原点距离为的点的位置,连接各象限内的相应点形成封闭曲线。
确定柱状节理裂隙岩体的REV时,可以先确定第一象限内不同倾角方向上到原点距离为的点的位置,连接相应点形成开放曲线,根据对称性绘制第二、第三、第四象限内的开放曲线,相邻开放曲线在坐标轴上相交形成封闭曲线;
确定一般裂隙岩体的REV时,可以先确定第一象限和第二象限内不同倾角方向的值为点到原点的距离做出开放曲线,根据对称性绘制其余象限内的开放曲线,相邻开放曲线在坐标轴上相交形成封闭曲线。
(4)根据该封闭曲线与椭圆的近似度,定量分析裂隙岩体不同方向的渗透系数是否符合渗透张量的规律;
当封闭曲线与椭圆的近似度大于90%时,说明裂隙岩体不同方向的渗透系数符合渗透张量的规律。
可以通过对封闭曲线进行函数拟合来确定其与椭圆的近似度,具体方法为:通过迭代法算出封闭曲线所对应的拟合度最高的椭圆方程,然后用得出的椭圆方程对曲线进行拟合,得到曲线对该椭圆方程的拟合度,即曲线与椭圆的近似程度。
(5)基于渗透张量性质,根据定量分析结果,扩大或缩小选取的实际裂隙岩体的体积,重复步骤(1)~(4),确定所测实际裂隙岩体的表征单元体积。
当定量分析结果为裂隙岩体不同方向的渗透系数符合渗透张量的规律时,缩小选取的实际裂隙岩体的体积,重复步骤(1)~(4),直至裂隙岩体不同方向的渗透系数不符合渗透张量的规律,最后两次选取的实际裂隙岩体体积的平均值即为所测实际裂隙岩体的表征单元体积。
当定量分析结果为裂隙岩体不同方向的渗透系数不符合渗透张量的规律时,扩大选取的实际裂隙岩体的体积,重复步骤(1)~(4),直至裂隙岩体不同方向的渗透系数符合渗透张量的规律,最后两次选取的实际裂隙岩体体积的平均值即为所测实际裂隙岩体的表征单元体积。此时,当选取的实际裂隙岩体的体积扩大到不合理尺寸时,裂隙岩体不同方向的渗透系数仍然不符合渗透张量的规律,说明该实际裂隙岩体的表征单元体积不存在。
实施例
以计算某水电站河床坝基处裂隙岩体的REV为例,本发明的确定裂隙岩体表征单元体积的方法,包括如下步骤:
(1)根据该水电站河床坝基处裂隙岩体的地质资料,确定选取体积为10m×10m×10m的裂隙岩体,设定比例尺缩小绘制其三维数字模型;
(2)将三维数字模型输入双喷头3D打印机,将该水电站河床坝基处的裂隙岩体制成粉末状作为打印岩块的材料,将裂隙岩体中裂隙的充填物制成粉末状作为打印岩体裂隙的材料,并装入双喷头3D打印机;
(3)将裂隙岩体的三维数字模型绕其中心点旋转成与水平方向成0°,15°,30°,45°,60°,75°,90°倾斜角,在不同的倾斜角下依次打印,分别打印出与水平方向成0°,15°,30°,45°,60°,75°,90°倾斜角的标准尺寸的模拟裂隙岩体的相似材料;
(4)利用岩石气体渗透率测试装置测出不同倾角裂隙岩体相似材料的渗透率Ka,以各方向的为点到原点的距离绘出曲线在坐标第一象限的部分,并根据对称性在其余各象限内补全曲线的剩余部分,如图2;
(5)通过迭代法算出裂隙岩体各方向绘成的曲线所对应的拟合度最高的椭圆方程是用该方程对曲线进行拟合,得到曲线对椭圆的拟合度是94.29%,即曲线对椭圆的近似程度是94.29%,各方向绘成的曲线对椭圆近似程度很高,说明裂隙岩体不同方向的渗透系数符合渗透张量的规律,则所选取裂隙岩体体积等于或大于REV;
(6)将选取的裂隙岩体实际尺寸缩小,重复上述步骤,直到裂隙岩体不同方向的渗透系数不符合渗透张量规律,取最后两次选取的裂隙岩体实际体积的平均值作为该裂隙岩体的REV,最终确定的岩体REV的值是3.2m×3.2m×3.2m。
采用本发明的方法得到某水电站河床坝基处裂隙岩体的REV为3.2m×3.2m×3.2m,与相关资料所得出的3m×3m×3m的结论大体一致。

Claims (5)

1.一种确定裂隙岩体表征单元体积的方法,其特征在于,包括如下步骤:
(1)选取任意体积的实际裂隙岩体,绘制其三维数字模型;
(2)旋转该三维数字模型,使其与水平方向成不同倾角,通过3D打印制得该实际裂隙岩体不同倾角的相似材料;
(3)测定不同倾角相似材料的渗透率Ka,建立坐标轴,以不同倾角方向的值为点到原点的距离,做出一个封闭曲线;
(4)根据该封闭曲线与椭圆的近似度,定量分析裂隙岩体不同方向的渗透系数是否符合渗透张量的规律;其中,通过对所述封闭曲线进行函数拟合来确定其与椭圆的近似度,确定所述封闭曲线与椭圆近似度的方法包括:通过迭代法算出所述封闭曲线对应的拟合度最高的椭圆方程,然后利用该椭圆方程对曲线进行拟合,得到的拟合度即为近似度;
(5)基于渗透张量性质,根据定量分析结果,扩大或缩小选取的实际裂隙岩体的体积,重复步骤(1)~(4),确定所测实际裂隙岩体的表征单元体积;
当定量分析结果为裂隙岩体不同方向的渗透系数符合渗透张量的规律时,缩小选取的实际裂隙岩体的体积,重复步骤(1)~(4),直至裂隙岩体不同方向的渗透系数不符合渗透张量的规律,最后两次选取的实际裂隙岩体体积的平均值即为所测实际裂隙岩体的表征单元体积;
当定量分析结果为裂隙岩体不同方向的渗透系数不符合渗透张量的规律时,扩大选取的实际裂隙岩体的体积,重复步骤(1)~(4),直至裂隙岩体不同方向的渗透系数符合渗透张量的规律,最后两次选取的实际裂隙岩体体积的平均值即为所测实际裂隙岩体的表征单元体积;
当选取的实际裂隙岩体的体积扩大到不合理尺寸时,裂隙岩体不同方向的渗透系数仍然不符合渗透张量的规律,说明该实际裂隙岩体的表征单元体积不存在。
2.根据权利要求1所述的确定裂隙岩体表征单元体积的方法,其特征在于,步骤(4)中,当所述封闭曲线与椭圆的近似度大于90%时,说明裂隙岩体不同方向的渗透系数符合渗透张量的规律。
3.根据权利要求1所述的确定裂隙岩体表征单元体积的方法,其特征在于,步骤(3)中,确定第一象限内不同倾角方向的值为点到原点的距离做出开放曲线,根据对称性绘制其余象限内的开放曲线,相邻开放曲线在坐标轴上相交形成封闭曲线。
4.根据权利要求1所述的确定裂隙岩体表征单元体积的方法,其特征在于,步骤(3)中,确定第一象限和第二象限内不同倾角方向的值为点到原点的距离做出开放曲线,根据对称性绘制其余象限内的开放曲线,相邻开放曲线在坐标轴上相交形成封闭曲线。
5.根据权利要求1所述的确定裂隙岩体表征单元体积的方法,其特征在于,步骤(2)中,采用所测实际裂隙岩体粉末为原材料进行3D打印,或者模拟所测实际裂隙岩体的力学性能,配制3D打印的原材料。
CN201511022607.0A 2015-12-29 2015-12-29 一种确定裂隙岩体表征单元体积的方法 Active CN105651964B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511022607.0A CN105651964B (zh) 2015-12-29 2015-12-29 一种确定裂隙岩体表征单元体积的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511022607.0A CN105651964B (zh) 2015-12-29 2015-12-29 一种确定裂隙岩体表征单元体积的方法

Publications (2)

Publication Number Publication Date
CN105651964A CN105651964A (zh) 2016-06-08
CN105651964B true CN105651964B (zh) 2017-11-03

Family

ID=56490873

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511022607.0A Active CN105651964B (zh) 2015-12-29 2015-12-29 一种确定裂隙岩体表征单元体积的方法

Country Status (1)

Country Link
CN (1) CN105651964B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106706393A (zh) * 2017-02-22 2017-05-24 河海大学 利用相似模型确定柱状节理岩体特征单元体rev的方法
CN110310363A (zh) * 2019-06-25 2019-10-08 上海电力学院 一种三维建模的表征单元体尺寸获取方法
CN115524275B (zh) * 2022-10-18 2023-08-22 中国水利水电科学研究院 考虑岩块渗流和裂隙渗流的裂隙岩体渗透张量确定方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130262028A1 (en) * 2012-03-30 2013-10-03 Ingrain, Inc. Efficient Method For Selecting Representative Elementary Volume In Digital Representations Of Porous Media
US9080946B2 (en) * 2012-06-15 2015-07-14 Ingrain, Inc. Digital rock analysis systems and methods with multiphase flow REV determination
CN202974813U (zh) * 2012-12-24 2013-06-05 长江水利委员会长江科学院 裂隙岩体渗流特性室内试验装置
CN103197042B (zh) * 2013-02-27 2015-07-15 北京科技大学 一种节理岩体表征单元体积的取值方法
CN105158039B (zh) * 2015-09-08 2018-05-08 河海大学 一种用于裂隙岩体渗流试验仿真裂隙制作方法

Also Published As

Publication number Publication date
CN105651964A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
Liu et al. Experimental investigation of the peak shear strength criterion based on three-dimensional surface description
Ghabraie et al. Application of 3D laser scanner, optical transducers and digital image processing techniques in physical modelling of mining-related strata movement
CN109211666B (zh) 基于ct扫描的预测应力加载条件下煤体渗透率的方法
Li et al. Evolution of the broken rock zone in the mixed ground tunnel based on the DSCM
Belayneh et al. Comparison of deterministic with stochastic fracture models in water-flooding numerical simulations
CN106837315B (zh) 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法
CN105651964B (zh) 一种确定裂隙岩体表征单元体积的方法
CN102288986B (zh) 一种地震尺度下碳酸盐岩储层弹性模量获取方法
Xiao et al. Combined effect of tortuosity and surface roughness on estimation of flow rate through a single rough joint
Dormieux et al. Homogenized models of stress-sensitive reservoir rocks
Marchand et al. A stochastic study of flow anisotropy and channelling in open rough fractures
Wang et al. A model of anisotropic property of seepage and stress for jointed rock mass
CN106707364B (zh) 一种雷达测井用标准井的施工方法
CN105298480A (zh) 缝洞性碳酸盐岩模型、其制备方法和用途
Wu et al. Fractal characteristics of low-permeability sandstone reservoirs
CN105717021B (zh) 一种定量确定裂隙岩体渗透各向异性的方法
CN105005667A (zh) 劈裂注浆复合土体等效弹性模量及泊松比计算方法
CN106248324B (zh) 长大软弱结构面强度参数的确定方法
CN110927041B (zh) 一种表征裂隙岩石渗透特性的层析成像方法
YAN et al. Research review of rock mechanics experiment and numerical simulation under THMC multi-field coupling
Hsu et al. Experimental study of the velocity field in solitary water waves
Iskander Transparent soils turn 25: Past, present, and future
Liu et al. Model test and 3 d numerical simulation study on excavation of double-arch tunnel.
CN105843992A (zh) 基于裂缝网络系统和渗流特征的页岩气数值模拟方法
Boldyrev et al. Particle image velocimetry and numeric analysis of sand deformations under a test plate

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant