CN110415167A - 一种基于数字图像技术的粗糙面裂隙生成方法及试验*** - Google Patents
一种基于数字图像技术的粗糙面裂隙生成方法及试验*** Download PDFInfo
- Publication number
- CN110415167A CN110415167A CN201910710425.4A CN201910710425A CN110415167A CN 110415167 A CN110415167 A CN 110415167A CN 201910710425 A CN201910710425 A CN 201910710425A CN 110415167 A CN110415167 A CN 110415167A
- Authority
- CN
- China
- Prior art keywords
- rough surface
- oil cylinder
- crack
- telescopic oil
- surface crack
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005516 engineering process Methods 0.000 title claims abstract description 33
- 239000011435 rock Substances 0.000 claims abstract description 72
- 238000004088 simulation Methods 0.000 claims abstract description 31
- 238000012360 testing method Methods 0.000 claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000002474 experimental method Methods 0.000 claims abstract description 6
- 230000009467 reduction Effects 0.000 claims abstract description 6
- BVPWJMCABCPUQY-UHFFFAOYSA-N 4-amino-5-chloro-2-methoxy-N-[1-(phenylmethyl)-4-piperidinyl]benzamide Chemical compound COC1=CC(N)=C(Cl)C=C1C(=O)NC1CCN(CC=2C=CC=CC=2)CC1 BVPWJMCABCPUQY-UHFFFAOYSA-N 0.000 claims abstract 3
- 238000010606 normalization Methods 0.000 claims abstract 3
- 230000011218 segmentation Effects 0.000 claims abstract 2
- 208000013668 Facial cleft Diseases 0.000 claims 1
- 238000007788 roughening Methods 0.000 abstract description 4
- 239000012530 fluid Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 239000000463 material Substances 0.000 description 7
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 4
- 239000010931 gold Substances 0.000 description 4
- 229910052737 gold Inorganic materials 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 2
- 206010000117 Abnormal behaviour Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/04—Context-preserving transformations, e.g. by using an importance map
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0256—Triaxial, i.e. the forces being applied along three normal axes of the specimen
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10008—Still image; Photographic image from scanner, fax or copier
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种基于数字图像技术的粗糙面裂隙生成方法及试验***，涉及岩石力学试验技术领域，步骤包括：A.加载岩石试件，试件内形成贯穿裂隙，清除试件表面碎屑；B.扫描仪对粗糙裂隙表面进行扫描，获得粗糙裂隙表面的数字图像；C.处理所述数字图像，包括去除噪声处理和阀值分割处理；D.对数字图像RGB值进行归一化处理，获得归一化矩阵，并建立空间坐标系，获取图像数字矩阵还原比例；E.计算数字还原矩阵；F.完成粗糙面裂隙的模拟，还原完整的空间裂隙；其试验***包括壳体、立柱、伸缩油缸阵列、控制***和液压回路。本发明提供了一种研究岩石粗糙裂隙渗流物理力学规律的方法，准确描述了岩石裂隙的粗糙性结构。
Description
技术领域
本发明涉及岩石力学试验技术领域，尤其是一种基于数字图像技术的粗糙面裂隙生成方法，以及模拟粗糙面裂隙的试验***。
背景技术
岩体的天然裂隙结构和渗流行为异常复杂,岩体裂隙渗流机制与定量描述一直是岩土、矿业、地质、石油及天然气工程高度关注的难点问题，并且岩体中的节理、裂隙等是地下水渗流的主要通道，岩体裂隙表面粗糙性对地下水渗流的影响较大。天然岩体中的裂隙长期受复杂的地质作用影响，其表面的粗糙特性千差万别，在数值模拟中很难直观的表征岩体表面的粗糙性，因此将裂隙表面简化为两个平行板，假设裂隙表面光滑，或者应用统计数中的概率密度函数假定一种概率分布来描述来处理岩体裂隙表面的非均匀性，在一定程度上研究了表面粗糙度对渗流的影响，但对于特定的岩体表面裂隙，目前仍然无法反映其实际的表面起伏结构。岩石各组分的强度不同使得岩石破裂面的粗糙度发生较大差别，只有准确的描述岩石裂隙表面的起伏粗糙性结构才能更客观的反映岩石应力-渗流耦合特征。
数字图像处理技术作为一种材料细观尺度上的空间结构精确测量和数字表述手段，已经广泛应用于岩土、混凝土的细观结构定量分析中。数字图像技术通过将数码设备获得的图形信息进行分类提取，通过归一化处理转化成0-1之间的数值，再根据相应的物理参数实测信息，对物理参数进行还原。该技术现阶段的主要应用包括岩石裂隙宽度非接触测量和数字表述岩体非均匀性等。
现有的天然岩体粗糙裂隙表面模拟，大多采用树脂材料复制岩石表面特征，利用相似材料重构该裂隙面，这种方法无法避免试验中的不确定因素对裂隙重构的干扰，并且在流固耦合试验中无法对相同裂隙情况的两个试件分别进行不同的试验，并且无法准确的获得节理面水力开度，从细观尺度研究岩石的裂隙岩体流固耦合问题缺乏必要的研究方法及试验手段；为此需要对现有的试验方法及试验***做进一步的改进。
发明内容
为了提供一种研究岩石粗糙裂隙渗流物理力学规律的方法，准确描述岩石裂隙的粗糙性结构，从细观尺度研究岩石的裂隙岩体流固耦合问题提供方便，本发明提供了一种基于数字图像技术的粗糙面裂隙生成方法及试验***，具体技术方案如下。
一种基于数字图像技术的粗糙面裂隙生成方法，步骤包括：
步骤A.制作岩石试件，并在岩石力学试验机上加载岩石试件，试件内形成贯穿裂隙，然后清除岩石试件表面碎屑；
步骤B.通过扫描仪对粗糙裂隙表面进行扫描，获得粗糙裂隙表面的数字图像；
步骤C.处理所述数字图像，包括去除噪声处理和阀值分割处理；
步骤D.对处理后的数字图像RGB值进行归一化处理，获得数字图像表面的归一化矩阵，建立空间坐标系，测量坐标系原点的起伏高度，获取图像数字矩阵还原比例；
步骤E.计算得到数字还原矩阵，包括用所述还原比例乘以所述归一化矩阵；
步骤F.利用数字还原矩阵和基于数字图像技术的粗糙面裂隙生成试验***，完成粗糙面裂隙的模拟，还原完整的空间裂隙。
优选的是，岩石试件包括长方体试件、正方体试件和圆柱形试件，岩石力学试验机包括单轴加载试验机和三轴加载试验机；岩石试件表面的碎屑使用干燥的毛刷清理。扫描仪为激光扫描仪或CT扫描仪，数字图像为位图或向量图。
优选的是，空间坐标系的建立选取裂隙面上的确定点作为原点。
进一步优选的是，归一化处理包括对亮度值I进行归一化处理，区分岩体表面裂隙是下凹还是上凸；当以岩石表面下凹最低点为水平面，则数字图像的I值均大于0。
进一步优选的是，原点的I值为hm(x)，选取岩体表面裂隙下凹最低点所在平面为基准面，实测原点与基准面之间的高程值hs(x)，获取还原比例因子。
还优选的是，基于数字图像技术的粗糙面裂隙生成试验***，包括壳体、立柱、伸缩油缸阵列、控制***和液压回路，壳体包括上壳体和下壳体，立柱连接上壳体和下壳体，伸缩油缸阵列包括上部伸缩油缸组和下部伸缩油缸组，上部伸缩油缸组和下部伸缩油缸组均设置有多个紧密排列的伸缩油缸，控制***控制液压回路工作，液压回路控制伸缩油缸的伸缩。
还进一步优选的是，上壳体和下壳体的边缘设置有多个立柱，立柱为伸缩结构，立柱同时伸缩调整上壳体和下壳体之间的距离。
还进一步优选的是，上部伸缩油缸组和下部伸缩油缸组的边缘设置有挡板，并且上部伸缩油缸组和下部伸缩油缸组的伸缩油缸相对布置并相互配合；伸缩油缸的截面为矩形。
更进一步优选的是，粗糙面裂隙的模拟包括上部粗糙面裂隙的模拟和下部粗糙面裂隙的模拟；
具体是，下部粗糙面裂隙的模拟中，先获取下部粗糙面裂隙的数字还原矩阵记为[A]，对[A]进行归一化处理得到通过Matlab确定出中一个最小值点和一个最大值点的坐标位置，确定下部裂隙的基准面与岩石试件底面的距离h1，同时确定最小值点和最大值点相对于基准面的高度分别记为ξ1和ξ2，得到下部粗糙面裂隙表面任意一点的在坐标系中的位态函数：
上部粗糙面裂隙的模拟中，先获取下部粗糙面裂隙的数字还原矩阵记为[B]，对[B]进行归一化处理得到通过Matlab确定出中一个最小值点和一个最大值点的坐标位置，确定上部裂隙的基准面与岩石试件底面的距离h2，同时确定最小值点和最大值点相对于基准面的高度分别记为η1和η2，得到上部粗糙面裂隙表面任意一点的在坐标系中的位态函数：
将下部粗糙面裂隙和上部粗糙面裂隙合并，在岩石试件底面选定原点建立新的空间坐标系，则裂隙开度空间分布函数ω(x)为：
将数字还原矩阵通过控制***的电脑输入，基于数字图像技术的粗糙面裂隙生成试验***还原完整的空间裂隙，得到裂隙空间分布模型；
裂隙空间分布模型应用于粗糙面裂隙岩体渗流模拟试验。
一种用于粗糙面裂隙生成方法的试验***，包括壳体、立柱、伸缩油缸阵列、控制***和液压回路，所述壳体包括上壳体和下壳体，所述立柱连接上壳体和下壳体，所述伸缩油缸阵列包括上部伸缩油缸组和下部伸缩油缸组，上部伸缩油缸组和下部伸缩油缸组均设置有多个紧密排列的伸缩油缸，所述控制***控制液压回路工作，所述液压回路控制伸缩油缸的伸缩。
本发明的有益效果包括：
(1)提供了一种基于数字图像技术的粗糙面裂隙生成方法，该方法通过对数字图像的处理，进而精确的重现岩石试件在拉伸、扭转或剪切破坏的粗糙面裂隙形貌特征，并真实的模拟岩体节理粗糙裂隙面起伏变化特征。
(2)本方法通过获取数字图像，并对数字图像进行数据处理，得到数字还原矩阵；利用该矩阵在基于数字图像技术的粗糙面裂隙生成试验***中，模拟还原粗糙裂隙面；通过该***还原还具有准确性高、成本低廉、可重复性好、操作简便等优点；并且该***可以重复多次还原多个不同的粗糙裂隙面，为室内岩石力学试验提供了更为便利的试验方法。
(3)该方法对裂隙表面图像进行处理，很好的还原了岩体裂隙的分布，利用基于数字图像技术的粗糙面裂隙生成试验***，还原完整的空间裂隙得到的裂隙空间分布模型，应用于渗流试验中，可以准确的获得节理面水力开度，有效的模拟裂隙流流动过程，避免了岩石粗糙裂隙离散性大，可重复性差的缺陷，为从细观尺度研究岩石的裂隙岩体流固耦合问题提供参考。另外本方法及试验***能够完成不同条件下结构面物理力学性能测试的试验，还在渗流试验中具有广阔的应用前景。
附图说明
图1为基于数字图像技术的粗糙面裂隙生成方法流程图；
图2为粗糙裂隙表面数字图像和归一化图像示意图；
图3为采集的粗糙裂隙表面数字图像的灰度分布规律图；
图4为粗糙裂隙表面结构剖面示意图；
图5为构建空间裂隙分布特征而建立的空间坐标系示意图；
图6为上部粗糙面裂隙的模拟和下部粗糙面裂隙的模拟原理示意图；
图7为裂隙模拟三维模型图；
图8为基于数字图像技术的粗糙面裂隙生成试验***示意图；
图9为试验***外部结构示意图；
图10为上部伸缩油缸组和下部伸缩油缸组结构示意图；
图11为伸缩油缸阵列的局部放大图；
图12为模拟粗糙裂隙表面截面拆分图；
图中：1-壳体；11-上壳体；12-下壳体；2-立柱；3-伸缩油缸阵列；4-控制***；5-液压回路。
具体实施方式
结合图1至图12所示，本发明提供的一种基于数字图像技术的粗糙面裂隙生成方法及试验***具体实施方式如下。
岩体粗糙裂隙表面的模拟，大多采用树脂材料复制岩石表面特征，再利用相似材料重构该裂隙面，这种方法无法避免试验中的多次复制及制作过程中不确定因素对裂隙重构的干扰，在流固耦合试验中无法对相同裂隙情况的两个试件分别进行不同的试验，并且难以准确的获得节理面水力开度，从细观尺度研究岩石的裂隙岩体流固耦合问题缺乏必要的研究方法及试验手段。本发明通过获取裂隙面的数字图像，建立数字模型，利用岩石裂隙表面某一点的实际测量结果，得到裂隙面对应点的数据值，获得图像和实际裂隙高差的比例因子对表面起伏信息进行还原，通过数字还原矩阵得到裂隙空间分布模型，裂隙空间分布模型进行流固耦合模拟试验，结果表明裂隙流流动速度明显小于相似材料的模拟，并且岩石裂隙数值模拟中能准确的获得节理面水力开度，有效的模拟裂隙流流动过程。
一种基于数字图像技术的粗糙面裂隙生成方法，如图1所示，其步骤包括：
步骤A.
制作岩石试件，并在岩石力学试验机上加载岩石试件，试件内形成贯穿裂隙，然后清除岩石试件表面碎屑。
岩石试件包括长方体试件、正方体试件和圆柱形试件等多种形状及尺寸的试件，试件由现场岩石取样制作而成，岩石力学试验机包括单轴加载试验机和三轴加载试验机等试验机，试验机加载试件得到贯穿裂隙。岩石试件表面的碎屑使用干燥的毛刷清理，从而消除碎屑对裂隙表面成像的干扰。
步骤B.
通过扫描仪对粗糙裂隙的表面进行扫描，获得粗糙裂隙表面的数字图像，图像为.jpg格式或.bmp格式及其他可处理的图像格式。其中扫描仪为激光扫描仪、岩石扫描仪或CT扫描仪，数字图像为位图或向量图。
步骤C.
处理所述数字图像，包括去除噪声处理和阀值分割处理，从而保证裂隙面数字图像质量。其中，得到的数字图像具有明显的峰值，可以作为确定的分割阀值。图像在生成和传输过程中常因受到各种噪声的干扰和影响而使图像降质，噪声种类很多包括电噪声、机械噪声、信道噪声和其他噪声，这对后续图像的处理(如分割、压缩和图像理解等)将产生不利影响；常用去除图像噪声的方法包括折叠均值滤波器去噪、折叠自适应维纳滤波器去噪、折叠中值滤波器去噪、形态学噪声滤除器去噪和小波去噪等。
步骤D.
对处理后的数字图像RGB值进行归一化处理，获得数字图像表面的归一化矩阵，建立空间坐标系，测量坐标系原点的起伏高度，获取图像数字矩阵还原比例。
图像RGB值即亮度，并使用整数来表示。通常情况下，RGB各有256级亮度，用数字表示为从0～255。灰度值是指指黑白图像中点的颜色深度，范围一般从0～255，白色为255，黑色为0。因此获取的数字图像可以为黑白图像或彩色图像，满足该处理条件的图像质量要求即可。
其中，空间坐标系的建立选取裂隙面上的确定点作为原点。归一化处理包括对亮度值I进行归一化处理，区分岩体表面裂隙是下凹还是上凸；当以岩石表面下凹最低点为水平面，则数字图像的I值均大于0。原点的I值为hm(x)，选取岩体表面裂隙下凹最低点所在平面为基准面，实测原点与基准面之间的高程值hs(x)，获取还原比例因子。
以具体的试验为例，粗糙裂隙表面数字图像灰度图像分布频率分布如图3所示，对图像的亮度值I进行归一化处理，用来区分岩体表面裂隙是下凹还是上凸。图2给出了三组裂隙表面图像I值归一化分布图，以岩石表面下凹最低点为水平面，该水平面作为参考平面的选取如图4所示，所以图中I值均大于0，I值分布介于0.1-0.9之间，对比图2中给定的裂隙表面和归一化数字图像，得到给定的岩石裂隙表面图像与生成的I值图像完全对等，说明基于该阀值和图像处理技术的I值分割方法在岩石裂隙表面起伏性的表征中表现良好。
其中，对获取的I值进行归一化处理后，得到图像数据矩阵Fig(x)，本试验中通过建立300mm×150mm×150mm的计算模型。然后获取图像中某点的I值hm(x)，以模型中的下凹最深度点所在水平面为基准，通过测量方法获得图像中选取点对应的高程实际值hs(x)。
步骤E.
计算得到数字还原矩阵，包括用所述还原比例乘以所述归一化矩阵。
以上述试验为例在计算出还原比例因子后，便可以根据获得的I值函数Fig(x)对裂隙表面图像进行还原，得到裂隙表面的实际函数H(x)，粗糙裂隙表面如图4所示。
步骤F.
利用数字还原矩阵和基于数字图像技术的粗糙面裂隙生成试验***，完成粗糙面裂隙的模拟，还原完整的空间裂隙。
其中，粗糙面裂隙的模拟包括上部粗糙面裂隙的模拟和下部粗糙面裂隙的模拟。
下部粗糙面裂隙的模拟中，先获取下部粗糙面裂隙的数字还原矩阵记为[A]，对[A]进行归一化处理得到通过Matlab确定出中一个最小值点和一个最大值点的坐标位置，确定下部裂隙的基准面与岩石试件底面的距离h1，同时确定最小值点和最大值点相对于基准面的高度分别记为ξ1和ξ2，得到下部粗糙面裂隙表面任意一点的在坐标系中的位态函数：
上部粗糙面裂隙的模拟中，先获取下部粗糙面裂隙的数字还原矩阵记为[B]，对[B]进行归一化处理得到通过Matlab确定出中一个最小值点和一个最大值点的坐标位置，确定上部裂隙的基准面与岩石试件底面的距离h2，同时确定最小值点和最大值点相对于基准面的高度分别记为η1和η2，得到上部粗糙面裂隙表面任意一点的在坐标系中的位态函数：
将下部粗糙面裂隙和上部粗糙面裂隙合并，在岩石试件底面选定原点建立新的空间坐标系，如图5至图7所示，则裂隙开度空间分布函数ω(x)为：
将数字还原矩阵通过控制***的电脑输入，基于数字图像技术的粗糙面裂隙生成试验***还原完整的空间裂隙，得到裂隙空间分布模型；该裂隙空间分布模型可以应用于粗糙面裂隙岩体渗流模拟试验。
上述基于数字图像技术的粗糙面裂隙生成方法，通过对数字图像的处理，进而精确的重现岩石试件在拉伸、扭转或剪切破坏的粗糙面裂隙形貌特征，并真实的模拟岩体节理粗糙裂隙面起伏变化特征；还通过获取数字图像，并对数字图像进行数据处理，得到数字还原矩阵；利用该矩阵在基于数字图像技术的粗糙面裂隙生成试验***中，模拟还原粗糙裂隙面；通过该***还原还具有准确性高、成本低廉、可重复性好、操作简便等优点；并且该***可以重复多次还原多个不同的粗糙裂隙面，为室内岩石力学试验提供了更为便利的试验方法。
如图8至图12所示，基于数字图像技术的粗糙面裂隙生成试验***，具体结构包括壳体1、立柱2、伸缩油缸阵列3、控制***4和液压回路5，其壳体1包括上壳体11和下壳体12，立柱2连接上壳体11和下壳体12，伸缩油缸阵列3包括上部伸缩油缸组和下部伸缩油缸组，上部伸缩油缸组和下部伸缩油缸组均设置有多个紧密排列的伸缩油缸，控制***4控制液压回路5工作，液压回路5控制伸缩油缸的伸缩。壳体1具有稳定的结构，上下壳体组合的结构包覆整个试验装置从而避免该***受到外界环境的干扰；立柱保证上下壳体的连接和组合，伸缩缸阵列在其面积和体积范围内可以模拟任意组合的粗糙面裂隙，控制***和控制回路用于伸缩油缸的精准控制。
其中上壳体11和下壳体12的边缘设置有多个立柱2，多个立柱2可以沿上壳体11和下壳体12的边缘对称布置，其立柱2可以为伸缩结构，立柱2同时伸缩调整上壳体和下壳体之间的距离，从而方便裂隙模拟的操作和观察。上部伸缩油缸组和下部伸缩油缸组的边缘设置有挡板，对上部伸缩油缸组和下部伸缩油缸组进行限位，并且上部伸缩油缸组和下部伸缩油缸组的伸缩油缸相对布置并相互配合，如图10和图12所示，保证上部伸缩油缸组和下部伸缩油缸组的伸缩油缸能够逐个整齐对应。伸缩油缸的截面为矩形，相连的伸缩油缸之间还需要做放水处理，方便进行粗糙面裂隙岩体渗流模拟试验。
该方法对裂隙表面图像进行处理，很好的还原了岩体裂隙的分布，利用基于数字图像技术的粗糙面裂隙生成试验***，还原完整的空间裂隙得到的裂隙空间分布模型，若应用于渗流试验中，可以准确的获得节理面水力开度，有效的模拟裂隙流流动过程，避免了岩石粗糙裂隙离散性大，可重复性差的缺陷，为从细观尺度研究岩石的裂隙岩体流固耦合问题提供参考。另外本方法及试验***能够完成不同条件下结构面物理力学性能测试的试验，还在渗流试验中具有广阔的应用前景。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (10)
1.一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，步骤包括：
步骤A.制作岩石试件，并在岩石力学试验机上加载岩石试件，试件内形成贯穿裂隙，然后清除岩石试件表面碎屑；
步骤B.通过扫描仪对粗糙裂隙表面进行扫描，获得粗糙裂隙表面的数字图像；
步骤C.处理所述数字图像，包括去除噪声处理和阀值分割处理；
步骤D.对处理后的数字图像RGB值进行归一化处理，获得数字图像表面的归一化矩阵，建立空间坐标系，测量坐标系原点的起伏高度，获取图像数字矩阵还原比例；
步骤E.计算得到数字还原矩阵，包括用所述还原比例乘以所述归一化矩阵；
步骤F.利用数字还原矩阵和基于数字图像技术的粗糙面裂隙生成试验***，完成粗糙面裂隙的模拟，还原完整的空间裂隙。
2.根据权利要求1所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述岩石试件包括长方体试件、正方体试件和圆柱形试件，所述岩石力学试验机包括单轴加载试验机和三轴加载试验机；岩石试件表面的碎屑使用干燥的毛刷清理；所述扫描仪为激光扫描仪或CT扫描仪，所述数字图像为位图或向量图。
3.根据权利要求1所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述空间坐标系的建立选取裂隙面上的确定点作为原点。
4.根据权利要求3所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述归一化处理包括对亮度值I进行归一化处理，区分岩体表面裂隙是下凹还是上凸；当以岩石表面下凹最低点为水平面，则数字图像的I值均大于0。
5.根据权利要求4所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述原点的I值为hm(x)，选取岩体表面裂隙下凹最低点所在平面为基准面，实测原点与基准面之间的高程值hs(x)，获取还原比例因子。
6.根据权利要求1-5任一项所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述基于数字图像技术的粗糙面裂隙生成试验***，包括壳体、立柱、伸缩油缸阵列、控制***和液压回路，所述壳体包括上壳体和下壳体，所述立柱连接上壳体和下壳体，所述伸缩油缸阵列包括上部伸缩油缸组和下部伸缩油缸组，上部伸缩油缸组和下部伸缩油缸组均设置有多个紧密排列的伸缩油缸，所述控制***控制液压回路工作，所述液压回路控制伸缩油缸的伸缩。
7.根据权利要求7所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述上壳体和下壳体的边缘设置有多个立柱，所述立柱为伸缩结构，立柱同时伸缩调整上壳体和下壳体之间的距离。
8.根据权利要求7所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述上部伸缩油缸组和下部伸缩油缸组的边缘设置有挡板，并且上部伸缩油缸组和下部伸缩油缸组的伸缩油缸相对布置并相互配合；所述伸缩油缸的截面为矩形。
9.根据权利要求7或8所述的一种基于数字图像技术的粗糙面裂隙生成方法，其特征在于，所述粗糙面裂隙的模拟包括上部粗糙面裂隙的模拟和下部粗糙面裂隙的模拟；
具体是，下部粗糙面裂隙的模拟中，先获取下部粗糙面裂隙的数字还原矩阵记为[A]，对[A]进行归一化处理得到通过Matlab确定出中一个最小值点和一个最大值点的坐标位置，确定下部裂隙的基准面与岩石试件底面的距离h1，同时确定最小值点和最大值点相对于基准面的高度分别记为ξ1和ξ2，得到下部粗糙面裂隙表面任意一点的在坐标系中的位态函数：
上部粗糙面裂隙的模拟中，先获取下部粗糙面裂隙的数字还原矩阵记为[B]，对[B]进行归一化处理得到通过Matlab确定出中一个最小值点和一个最大值点的坐标位置，确定上部裂隙的基准面与岩石试件底面的距离h2，同时确定最小值点和最大值点相对于基准面的高度分别记为η1和η2，得到上部粗糙面裂隙表面任意一点的在坐标系中的位态函数：
将下部粗糙面裂隙和上部粗糙面裂隙合并，在岩石试件底面选定原点建立新的空间坐标系，则裂隙开度空间分布函数ω(x)为：
将数字还原矩阵通过控制***的电脑输入，基于数字图像技术的粗糙面裂隙生成试验***还原完整的空间裂隙，得到裂隙空间分布模型；
所述裂隙空间分布模型应用于粗糙面裂隙岩体渗流模拟试验。
10.一种用于粗糙面裂隙生成方法的试验***，其特征在于，包括壳体、立柱、伸缩油缸阵列、控制***和液压回路，所述壳体包括上壳体和下壳体，所述立柱连接上壳体和下壳体，所述伸缩油缸阵列包括上部伸缩油缸组和下部伸缩油缸组，上部伸缩油缸组和下部伸缩油缸组均设置有多个紧密排列的伸缩油缸，所述控制***控制液压回路工作，所述液压回路控制伸缩油缸的伸缩。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910710425.4A CN110415167B (zh) | 2019-08-02 | 2019-08-02 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验*** |
| PCT/CN2019/099350 WO2021022465A1 (zh) | 2019-08-02 | 2019-08-06 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910710425.4A CN110415167B (zh) | 2019-08-02 | 2019-08-02 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110415167A true CN110415167A (zh) | 2019-11-05 |
| CN110415167B CN110415167B (zh) | 2020-03-06 |
Family
ID=68365375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910710425.4A Active CN110415167B (zh) | 2019-08-02 | 2019-08-02 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验*** |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110415167B (zh) |
| WO (1) | WO2021022465A1 (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111507988A (zh) * | 2020-04-14 | 2020-08-07 | 山东科技大学 | 一种基于数字图像处理的破碎岩体建模及渗流试验方法 |
| CN111739149A (zh) * | 2020-06-15 | 2020-10-02 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
| CN112131642A (zh) * | 2020-09-17 | 2020-12-25 | 四川大学 | 一种基于波数法的岩体粗糙离散裂隙网络生成方法 |
| CN112435217A (zh) * | 2020-11-02 | 2021-03-02 | 上海祐云信息技术有限公司 | 图像识别预制钢筋混泥土粗糙面的算法 |
| CN116167119A (zh) * | 2023-04-24 | 2023-05-26 | 河海大学 | 二级配混凝土破裂面分形维数计算方法 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113899767B (zh) * | 2021-10-08 | 2024-03-12 | 山东科技大学 | 基于3d打印技术的巷道围岩实验模型的成型实验方法 |
| CN113920141B (zh) * | 2021-12-15 | 2022-02-11 | 中南大学 | 一种岩体完整性系数的计算方法、装置及存储介质 |
| CN114675004B (zh) * | 2022-03-22 | 2023-03-24 | 中国矿业大学 | 一种采动三维裂缝模拟板制备装置及突水溃砂模拟方法 |
| CN115015251B (zh) * | 2022-06-01 | 2024-03-12 | 青岛理工大学 | 一种多作用力下的可视化三维裂隙注浆实验***及方法 |
| CN115370419A (zh) * | 2022-08-01 | 2022-11-22 | 浙江大学 | 一种考虑岩石裂隙生物加固-渗流-温度耦合作用的可视化试验***与方法 |
| CN117313589B (zh) * | 2023-11-30 | 2024-02-06 | 山东科技大学 | 三维粗糙离散裂隙网络构建、评价与渗流模拟的方法 |
| CN118111808B (zh) * | 2024-04-30 | 2024-07-09 | 煤炭科学研究总院有限公司 | 一种岩体稳定性分析的试验方法及*** |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419303A (zh) * | 2011-08-15 | 2012-04-18 | 山东科技大学 | 复杂条件下的裂隙注浆可视化试验装置 |
| CN102607950A (zh) * | 2012-03-05 | 2012-07-25 | 山东科技大学 | 一种岩石剪切渗流耦合真三轴试验*** |
| CN102855806A (zh) * | 2012-09-18 | 2013-01-02 | 中国水利水电科学研究院 | 一种物理模型快速搭建***和方法 |
| CN103698171A (zh) * | 2013-11-12 | 2014-04-02 | 山东科技大学 | 底板岩层固流耦合材料的裂隙演化试验方法 |
| CN103954511A (zh) * | 2014-04-09 | 2014-07-30 | 北京工业大学 | 一种裂隙网络岩石剪切-渗流耦合实验方法 |
| CN105486581A (zh) * | 2015-12-14 | 2016-04-13 | 山东科技大学 | 一种反演岩石预制裂隙发育过程的方法 |
| CN106228559A (zh) * | 2016-07-26 | 2016-12-14 | 绍兴文理学院 | 基于分解‑重构建立裂隙岩体三维数值模型的方法 |
| CN106872241A (zh) * | 2017-04-28 | 2017-06-20 | 山东科技大学 | 一种自动调节的岩石或混凝土圆柱试样制作装置 |
| CN108613872A (zh) * | 2018-06-13 | 2018-10-02 | 山东科技大学 | 透明类岩石试件裂隙扩展监测***及方法 |
| CN108763650A (zh) * | 2018-04-28 | 2018-11-06 | 湘潭大学 | 一种覆岩采动裂隙网络模型构建方法 |
| CN109186935A (zh) * | 2018-08-08 | 2019-01-11 | 武汉大学 | 一种模拟交叉裂隙渗流的可视化试验装置及方法 |
| US20190212460A1 (en) * | 2018-01-08 | 2019-07-11 | Dagang Oil Field Of Cnpc | Method for secondary exploration of old oil area in fault subsidence basin |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5297420A (en) * | 1993-05-19 | 1994-03-29 | Mobil Oil Corporation | Apparatus and method for measuring relative permeability and capillary pressure of porous rock |
| WO2012118868A2 (en) * | 2011-02-28 | 2012-09-07 | Schlumberger Technology Corporation | Petrographic image analysis for determining capillary pressure in porous media |
| CN102621010B (zh) * | 2012-04-13 | 2014-01-29 | 山东大学 | 深部围岩锚固性能多功能试验机 |
| CN103871064B (zh) * | 2014-03-25 | 2017-02-08 | 中国石油大学(华东) | 一种火山岩ct图像的预处理和确定分割阈值的方法 |
| CN104677726B (zh) * | 2015-02-09 | 2017-05-31 | 四川大学 | 三爪卡盘限位拉压循环试验岩石试样固定装置 |
| CN108519264B (zh) * | 2018-04-09 | 2020-10-02 | 中国科学院武汉岩土力学研究所 | 一种裂隙试样制备方法 |
-
2019
- 2019-08-02 CN CN201910710425.4A patent/CN110415167B/zh active Active
- 2019-08-06 WO PCT/CN2019/099350 patent/WO2021022465A1/zh active Application Filing
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419303A (zh) * | 2011-08-15 | 2012-04-18 | 山东科技大学 | 复杂条件下的裂隙注浆可视化试验装置 |
| CN102607950A (zh) * | 2012-03-05 | 2012-07-25 | 山东科技大学 | 一种岩石剪切渗流耦合真三轴试验*** |
| CN102855806A (zh) * | 2012-09-18 | 2013-01-02 | 中国水利水电科学研究院 | 一种物理模型快速搭建***和方法 |
| CN103698171A (zh) * | 2013-11-12 | 2014-04-02 | 山东科技大学 | 底板岩层固流耦合材料的裂隙演化试验方法 |
| CN103954511A (zh) * | 2014-04-09 | 2014-07-30 | 北京工业大学 | 一种裂隙网络岩石剪切-渗流耦合实验方法 |
| CN105486581A (zh) * | 2015-12-14 | 2016-04-13 | 山东科技大学 | 一种反演岩石预制裂隙发育过程的方法 |
| CN106228559A (zh) * | 2016-07-26 | 2016-12-14 | 绍兴文理学院 | 基于分解‑重构建立裂隙岩体三维数值模型的方法 |
| CN106872241A (zh) * | 2017-04-28 | 2017-06-20 | 山东科技大学 | 一种自动调节的岩石或混凝土圆柱试样制作装置 |
| US20190212460A1 (en) * | 2018-01-08 | 2019-07-11 | Dagang Oil Field Of Cnpc | Method for secondary exploration of old oil area in fault subsidence basin |
| CN108763650A (zh) * | 2018-04-28 | 2018-11-06 | 湘潭大学 | 一种覆岩采动裂隙网络模型构建方法 |
| CN108613872A (zh) * | 2018-06-13 | 2018-10-02 | 山东科技大学 | 透明类岩石试件裂隙扩展监测***及方法 |
| CN109186935A (zh) * | 2018-08-08 | 2019-01-11 | 武汉大学 | 一种模拟交叉裂隙渗流的可视化试验装置及方法 |
Non-Patent Citations (2)
| Title |
|---|
| 尹立明 等: ""岩石应力–渗流耦合真三轴试验***的研制与应用"", 《岩石力学与工程学报》 * |
| 陈军涛 等: ""深部开采底板裂隙扩展演化规律试验研究"", 《岩石力学与工程学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111507988A (zh) * | 2020-04-14 | 2020-08-07 | 山东科技大学 | 一种基于数字图像处理的破碎岩体建模及渗流试验方法 |
| CN111739149A (zh) * | 2020-06-15 | 2020-10-02 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
| CN111739149B (zh) * | 2020-06-15 | 2023-09-01 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
| CN112131642A (zh) * | 2020-09-17 | 2020-12-25 | 四川大学 | 一种基于波数法的岩体粗糙离散裂隙网络生成方法 |
| CN112131642B (zh) * | 2020-09-17 | 2021-06-01 | 四川大学 | 一种基于波数法的岩体粗糙离散裂隙网络生成方法 |
| CN112435217A (zh) * | 2020-11-02 | 2021-03-02 | 上海祐云信息技术有限公司 | 图像识别预制钢筋混泥土粗糙面的算法 |
| CN112435217B (zh) * | 2020-11-02 | 2022-01-11 | 上海祐云信息技术有限公司 | 图像识别预制钢筋混凝土粗糙面的方法 |
| CN116167119A (zh) * | 2023-04-24 | 2023-05-26 | 河海大学 | 二级配混凝土破裂面分形维数计算方法 |
| CN116167119B (zh) * | 2023-04-24 | 2023-11-28 | 河海大学 | 二级配混凝土破裂面分形维数计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110415167B (zh) | 2020-03-06 |
| WO2021022465A1 (zh) | 2021-02-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110415167A (zh) | 一种基于数字图像技术的粗糙面裂隙生成方法及试验*** | |
| Ramandi et al. | Digital rock analysis for accurate prediction of fractured media permeability | |
| Shan et al. | Influence of CT scanning parameters on rock and soil images | |
| Cheng | A new model for quantifying anisotropic scale invariance and for decomposition of mixing patterns | |
| Shirole et al. | Experimental investigation of multi-scale strain-field heterogeneity in rocks | |
| Rupil et al. | Fatigue microcrack detection with digital image correlation | |
| Promma et al. | Application of the virtual fields method to mechanical characterization of elastomeric materials | |
| US10223782B2 (en) | Digital rock physics-based trend determination and usage for upscaling | |
| CN108693572B (zh) | 一种基于三维建模的地应力场反演方法 | |
| US20140019054A1 (en) | Digital rock analysis systems and methods with reliable multiphase permeability determination | |
| NO20121394A1 (no) | Fremgangsmåte for å innhente konsistente og integrerte fysiske egenskaper for porøse medier | |
| Zhao et al. | Digital measurement of 2D and 3D cracks in sandstones through improved pseudo color image enhancement and 3D reconstruction method | |
| CN111624146A (zh) | 一种定量表征致密储层全尺度孔径分布特征的方法 | |
| CN109063606B (zh) | 矿化蚀变遥感信息提取方法及装置 | |
| Wei et al. | Experimental techniques for the study of the cracking mechanisms in drying clays | |
| Wu et al. | Microfabric evolution of coral sand foundations during seismic liquefaction using 3D images | |
| CN105975990B (zh) | 一种复垦土壤剖面的精细表达方法 | |
| Zhang et al. | The development of a 3D mesoscopic model of metallic foam based on an improved watershed algorithm | |
| CN110660051B (zh) | 一种基于导航金字塔的张量投票处理方法 | |
| Kahl et al. | Microfabric and anisotropy of elastic waves in sandstone–An observation using high-resolution X-ray microtomography | |
| Eichler | C-Axis analysis of the NEEM ice core–An approach based on digital image processing | |
| Lyu et al. | An image processing approach to measuring the sphericity and roundness of fracturing proppants | |
| Dipova | Automated strain measurements in uniaxial testing via computer vision | |
| Moshfeghi et al. | Validating the use of material point method and SANISAND model for relating the state parameter with cone tip resistance | |
| Qi et al. | Thermal pyrolysis-induced shale micro-cracks: 3D characterization and implication for reservoir stimulation |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |