CN108061697B - 土体三维孔隙率计算方法 - Google Patents
土体三维孔隙率计算方法 Download PDFInfo
- Publication number
- CN108061697B CN108061697B CN201711278075.6A CN201711278075A CN108061697B CN 108061697 B CN108061697 B CN 108061697B CN 201711278075 A CN201711278075 A CN 201711278075A CN 108061697 B CN108061697 B CN 108061697B
- Authority
- CN
- China
- Prior art keywords
- soil
- threshold value
- porosity
- calculating
- image
- 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
Links
- 239000002689 soil Substances 0.000 title claims abstract description 190
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000011148 porous material Substances 0.000 claims abstract description 65
- 239000002245 particle Substances 0.000 claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000008859 change Effects 0.000 claims description 15
- 238000005520 cutting process Methods 0.000 claims description 14
- 238000005056 compaction Methods 0.000 claims description 4
- 239000002390 adhesive tape Substances 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 238000005498 polishing Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 abstract description 11
- 230000008569 process Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 3
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000011780 sodium chloride Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000007596 consolidation process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- NUHSROFQTUXZQQ-UHFFFAOYSA-N isopentenyl diphosphate Chemical compound CC(=C)CCO[P@](O)(=O)OP(O)(O)=O NUHSROFQTUXZQQ-UHFFFAOYSA-N 0.000 description 1
- 238000001000 micrograph Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N2015/0846—Investigating permeability, pore-volume, or surface area of porous materials by use of radiation, e.g. transmitted or reflected light
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供了一种土体三维孔隙率计算方法，属于室内土体孔隙率的测试领域，包括制备土样；利用扫描电镜获取土样断面的图像；设置二值化阈值序列，通过二值化阈值序列中的二值化阈值分别对图像进行二值化处理，并获取二值化处理后图像中的灰度面积；根据二值化阈值以及二值化阈值对应的灰度面积拟合灰度面积随二值化阈值变化的曲线；对灰度面积随二值化阈值变化的曲线方程进行积分，获得图像中土体表面微观结构的孔隙体积；获取图像中土体表面颗粒与孔隙整体体积；计算获得孔隙率。本发明提供的土体三维孔隙率计算方法的误差能够限定在合理的误差范围之内，测算结果精确度高，且测算过程简单，有利于提高测算效率。
Description
技术领域
本发明属于室内土体孔隙率的测试技术领域，更具体地说，是涉及一种土体三维孔隙率计算方法。
背景技术
土体是一种多相介质，具有孔隙、土骨架、土颗粒以及水分等复杂结构特征。在解释土体宏观特性方面，从土体微观结构研究土体的宏观性质是一种科学有效的手段。土体的孔隙大小及其分布是导致土体变形的主要原因之一，如地基沉降、边坡失稳等都与土体孔隙结构发生变化或孔隙与水等其他介质相互作用有关，因此观察土体的孔隙分布特征就显得尤为重要。
目前，对于研究土体孔隙的变化多是通过土的常规固结试验进行的，从微观结构观察孔隙的变化的方法研究尚不成熟。对于土体孔隙的微观观察采用的方式大多为某一确定阈值的方法。这种方法多是通过某一阈值下，通过扫描图像得到的计算孔隙率与常规室内试验标定的孔隙率一致时，确定观察土体断面图像的阈值，从而进一步研究在这一阈值下该图像呈现出的土体微观孔隙的分布。此方法首先在试验步骤上过于繁多冗杂，效率不高；其次对于在不方便确定室内试验标定时，不能精确地定量反映土体微观孔隙的数量多少和分布情况，存在阈值变化引起的观测系统误差。
发明内容
本发明的目的在于提供一种土体三维孔隙率计算方法，以解决现有技术中存在的土体孔隙的微观观察方法操作复杂，且精确度低的技术问题。
为实现上述目的，本发明采用的技术方案是：提供一种土体三维孔隙率计算方法，包括如下步骤：
制备土样；
利用扫描电镜获取土样断面的图像；
设置二值化阈值序列，通过所述二值化阈值序列中的二值化阈值分别对所述图像进行二值化处理，并获取二值化处理后图像中的灰度面积；
根据所述二值化阈值以及所述二值化阈值对应的灰度面积拟合灰度面积随二值化阈值变化的曲线；
对所述灰度面积随二值化阈值变化的曲线方程进行积分，获得所述图像中土体表面微观结构的孔隙体积；
获取所述图像中土体表面颗粒与孔隙整体体积；
通过获得的所述孔隙体积以及所述土体表面颗粒与孔隙整体体积计算获得孔隙率。
进一步地，所述对所述灰度面积随二值化阈值变化的曲线方程进行积分，获得所述图像中土体表面微观结构的孔隙体积包括：
在0-255的阈值范围内对所述灰度面积随二值化阈值变化的曲线方程进行积分，计算关系式为
其中，Ve为所述孔隙体积，T为阈值，S(T)某确定阈值下得到的所述灰度面积。
进一步地，所述获取所述图像中土体表面颗粒与孔隙整体体积包括：
当阈值为255时，计算所述土体表面颗粒与孔隙整体体积，计算关系式为
V＝255×S(255) (2)
其中，V为所述土体表面颗粒与孔隙整体体积，S(255)为阈值为255时得到的所述灰度面积。
进一步地，所述通过获得的所述孔隙体积以及所述土体表面颗粒与孔隙整体体积计算获得孔隙率包括：
通过获得的所述孔隙体积以及所述土体表面颗粒与孔隙整体体积的比值获得所述孔隙率，计算关系式为
其中，n为所述孔隙率。
进一步地，所述二值化阈值序列为梯度变化的二值化阈值序列，所述二值化阈值的变化梯度不大于10。
进一步地，所述二值化阈值的变化梯度为5或10。
进一步地，所述制备土样包括：
按照预设的土体孔隙率和压实度制备并获得原始土样；
切割所述原始土样，获得具有断面的多个半成土样；
分别清除所述半成土样断面上的浮动颗粒，获得多个成品土样。
进一步地，所述切割所述原始土样，获得具有断面的多个半成土样包括：
通过割土刀在所述原始土样中部切出沟槽；
沿所述沟槽手动掰断所述原始土样；
对所述断面进行打磨，获得所述半成土样。
进一步地，所述分别清除所述半成土样断面上的浮动颗粒包括：
利用胶带粘除所述半成土样断面上的浮动颗粒。
进一步地，所述分别清除所述半成土样断面上的浮动颗粒，获得多个成品土样之后还包括：
对所述成品土样进行干燥处理。
本发明提供的土体三维孔隙率计算方法的有益效果在于：与现有技术相比，本发明土体三维孔隙率计算方法，通过扫描电镜获得土体微观结构的图像，该图像能够精确清晰的呈现土体的微观结构，通过图像的二值化处理，每一个阈值对应的灰度面积(灰度图像实际上是经过二值化处理后，图像中黑色像素点总体的面积)相当于不同深度的孔隙的截面积，在二值化阈值序列范围内对灰度面积进行积分，即可得到孔隙的体积，再通过与土体表面颗粒与孔隙整体体积配合，即可方便的计算出孔隙率。经该方法测试的孔隙率与经过常规土工试验标定的土样的实际孔隙率非常接近，误差限定在了合理的误差范围之内，因此测算结果精确度高，且测算过程简单，有利于提高测算效率。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的土体三维孔隙率计算方法的流程图；
图2为利用本发明实施例提供的土体三维孔隙率计算方法进行测试的土样的扫描电镜图像；
图3为本发明实施例提供的土体三维孔隙率计算方法的原理图；
图4为利用本发明实施例提供的土体三维孔隙率计算方法进行测试的土样的灰度面积随二值化阈值变化的曲线；
图5为对土样的图像进行二值化处理的原理图一；
图6为对土样的图像进行二值化处理的原理图二。
具体实施方式
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请一并参阅图1至图6，现对本发明提供的土体三维孔隙率计算方法进行说明。所述土体三维孔隙率计算方法，包括如下步骤：
制备土样；
利用扫描电镜获取土样断面的图像；
设置二值化阈值序列，通过二值化阈值序列中的二值化阈值分别对图像进行二值化处理，并获取二值化处理后图像中的灰度面积；
根据二值化阈值以及二值化阈值对应的灰度面积拟合灰度面积随二值化阈值变化的曲线；
对灰度面积随二值化阈值变化的曲线方程进行积分，获得图像中土体表面微观结构的孔隙体积；
获取图像中土体表面颗粒与孔隙整体体积；
通过获得的孔隙体积以及土体表面颗粒与孔隙整体体积计算获得孔隙率。
其中，土样的尺寸根据扫描电镜载物台的尺寸制备，并且可以根据测试需求制备不同种类的土样。
需要注意的是，灰度面积也可以用该阈值下显示出的黑色的像素点的个数替代。
扫描电镜(SEM)获得的图像是灰度图像，黑色的灰度为0，白的的灰度为255，接近白色的浅色部分通常表示土体颗粒，接近黑色的深色部分为孔隙。扫描电镜照片的性质适合采用二值化处理方式，二值化是一种图象处理的重要方法，即将图像上像素点设置为0或者255，产生轮廓分明的黑白效果，最终将图像分割，只保留需要关注的部分。
本发明提供的土体三维孔隙率计算方法，与现有技术相比，通过扫描电镜获得土体微观结构的图像，该图像能够精确清晰的呈现土体的微观结构，通过图像的二值化处理，每一个阈值对应的灰度面积(灰度图像实际上是经过二值化处理后，图像中黑色像素点总体的面积)相当于不同深度的孔隙的截面积，在二值化阈值序列范围内对灰度面积进行积分，即可得到孔隙的体积，再通过与土体表面颗粒与孔隙整体体积配合，即可方便的计算出孔隙率。经该方法测试的孔隙率与经过常规土工试验标定的土样的实际孔隙率非常接近，误差限定在了合理的误差范围之内，因此测算结果精确度高，且测算过程简单，有利于提高测算效率。
进一步地，请一并参阅图1及图4，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，以盐渍土土样为例，根据二值化阈值以及二值化阈值对应的灰度面积拟合灰度面积随二值化阈值变化的曲线具体为：
建立以二值化阈值为横坐标、以灰度面积为纵坐标的坐标系，并在阈值为0-255的范围内分别绘制出灰度面积对应的数值，以数值点为基础获得原始曲线。土体扫描后的灰度面积随二值化阈值的变化规律曲线大致分为斜率不同的三段，很明显的发现随着二值化阈值的增大孔隙率增加，经过拟合发现盐渍土土样的曲线分为三段斜率不同的线性关系：
S0＝aiT+bi (7)
其中S0为区间内像素点的灰度面积，ai和bi为三个区间的关系参数。
T＜50时，表示接近黑色区域，由于拍摄限制，远离镜头无法拍摄的区域增多，若将该隐藏区域计算进表面孔隙度会增加误差，应当舍去隐藏部分，同时可以确定最大的计算深度；T≥50时，曲线呈斜率不同的线性关系，将不同斜率曲线方程的积分累加即可。
进一步地，请一并参阅图1至图6，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，对灰度面积随二值化阈值变化的曲线方程进行积分，获得图像中土体表面微观结构的孔隙体积包括：
在0-255的阈值范围内对灰度面积随二值化阈值变化的曲线方程进行积分，计算关系式为
其中，Ve为孔隙体积，T为阈值，S(T)某确定阈值下得到的灰度面积。
如图6所示，(a)表示土样的表面的走势，(b)中的阴影部分表示土体颗粒，空白部分表示孔隙，越是位于底部的部分，其灰度值越小(越接近0)，土体颗粒的顶部灰度值大，接近255。根据关系式
其中，x和y是图像横纵坐标，f(x,y)为未处理图像的灰度值，g(x,y)是二值化后的图像灰度值。当灰度值设定为某一阈值的时候，灰度小于该阈值的像素点显示为黑色，高于该阈值的像素点显示为白色。
如图3所示，该图中模拟了土体颗粒的空间分布侧视结构，由于微观结构下几何形状的足够小，为了便于计算，将真实的土体形状简化为等截面不规则三维立体，不规则孔隙在其内部不均匀分布。由于扫描电镜只能对土体单侧进行拍摄，不能拍摄到土颗粒之间遮藏的孔隙，因此以土颗粒接触点为基准，下部分作等截面立体处理，被遮藏的部分(图3中右侧的两个类似于楔形的阴影部分)忽略不计。图3中的T表示阈值，A表示孔隙。
建立以扫描电镜得到的二维图像为xy像素平面，以阈值T为z轴的三维坐标系，在土体中取出dT的微元体，微元体中阈值下孔隙的体积：
dv＝∫∫Tds (5)
其中，ds为该阈值下电镜扫描后得到的灰度图像的面积。S的值可以由IPP图像处理软件直接扫描得出，因此对(5)式进行整理，变换为对变量T的一元积分，有：
dv＝∫SdT (6)
对扫描得到的灰度图像进行积分，即可得到孔隙的体积Ve为关系式(1)。
进一步地，请参阅图1至图6，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，获取图像中土体表面颗粒与孔隙整体体积包括：
当阈值为255时，计算土体表面颗粒与孔隙整体体积，计算关系式为
V＝255×S(255) (2)
其中，V为土体表面颗粒与孔隙整体体积，S(255)为阈值为255时得到的灰度面积。
在计算过程中，将土体简化为等截面的立体，S(255)为图像定义为全白的时候，扫描得到的断面面积，包括孔隙与土体颗粒的所有面积，由于每一层的面积均一致，因此，可知V的计算公式为关系式(2)。
进一步地，参阅图1，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，通过获得的孔隙体积以及土体表面颗粒与孔隙整体体积计算获得孔隙率包括：
通过获得的孔隙体积以及土体表面颗粒与孔隙整体体积的比值获得孔隙率，计算关系式为
其中，n为所述孔隙率。
土体孔隙率的定义为：土体颗粒之间的孔隙体积占土体总体积的比率。根据该定义，将已经测算得到的与图像对应的孔隙体积除以与图像对应的土体总体积即可得到孔隙率。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，为了保证测算的精确度，二值化阈值序列为梯度变化的二值化阈值序列，二值化阈值的变化梯度不大于10。
优选地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，二值化阈值的变化梯度为5或10。在保证测算精度的前提下，能够节约计算的时间，提高计算效率。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，制备土样包括：
按照预设的土体孔隙率和压实度制备并获得原始土样；
切割原始土样，获得具有断面的多个半成土样；
分别清除半成土样断面上的浮动颗粒，获得多个成品土样。
一个土样通常被截断成两节，进而获得两个断面，在测试的时候，取两个断面上的扫描数据平均值作为本土样的扫描计算孔隙率，基本可以避免截断过程中对原有土样表面造成的不均匀性而导致的测试误差。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，按照预设的土体孔隙率和压实度制备并获得原始土样包括：
对同一种土体制备多个土样。有利于提高计算结果的精确度。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，利用扫描电镜获取土样断面的图像包括：
对同一种土体的多个土样在不同的断面位置进行拍摄，以获取更多的图像。
大量拍摄图像，并根据不同的图像进行计算，可以使数据方差及标准差降低至0.6％，进而有利于提高数据的可靠程度。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，切割原始土样，获得具有断面的多个半成土样包括：
通过割土刀在所述原始土样中部切出沟槽；
沿沟槽手动掰断所述原始土样；
对断面进行打磨，获得半成土样。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，分别清除半成土样断面上的浮动颗粒包括：
利用胶带粘除半成土样断面上的浮动颗粒。
该截断方式为切割后撕皮的方式，操作简单，表面较为光滑，有利于提高图像的清晰度，进而能够提高计算结果的准确性。
进一步地，作为本发明提供的土体三维孔隙率计算方法的一种具体实施方式，分别清除半成土样断面上的浮动颗粒，获得多个成品土样之后还包括：
对成品土样进行干燥处理。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.土体三维孔隙率计算方法，其特征在于，包括如下步骤：
制备土样；
利用扫描电镜获取土样断面的图像；
设置二值化阈值序列，通过所述二值化阈值序列中的二值化阈值分别对所述图像进行二值化处理，并获取二值化处理后图像中的灰度面积；
根据所述二值化阈值以及所述二值化阈值对应的灰度面积拟合灰度面积随二值化阈值变化的曲线具体为：
建立以二值化阈值为横坐标、以灰度面积为纵坐标的坐标系，并在阈值为0-255的范围内分别绘制出灰度面积对应的数值，获得曲线：
以数值点为基础获得所述曲线的三段斜率不同的线性关系：
S0＝aiT+bi (7)
其中S0为区间内像素点的灰度面积，T为阈值，ai和bi为三个区间的关系参数；
当阈值小于某一数值时，该曲线不符合式(7)的线性关系，舍去小于该数值的部分；
对所述灰度面积随二值化阈值变化的曲线方程进行积分，获得所述图像中土体表面微观结构的孔隙体积；
获取所述图像中土体表面颗粒与孔隙整体体积；
通过获得的所述孔隙体积以及所述土体表面颗粒与孔隙整体体积计算获得孔隙率。
2.如权利要求1所述的土体三维孔隙率计算方法，其特征在于，所述对所述灰度面积随二值化阈值变化的曲线方程进行积分，获得所述图像中土体表面微观结构的孔隙体积包括：
在0-255的阈值范围内对所述灰度面积随二值化阈值变化的曲线方程进行积分，计算关系式为
其中，Ve为所述孔隙体积，T为阈值，S(T)为 某确定阈值下得到的所述灰度面积。
3.如权利要求2所述的土体三维孔隙率计算方法，其特征在于，所述获取所述图像中土体表面颗粒与孔隙整体体积包括：
当阈值为255时，计算所述土体表面颗粒与孔隙整体体积，计算关系式为
V＝255×S(255) (2)
其中，V为所述土体表面颗粒与孔隙整体体积，S(255)为阈值为255时得到的所述灰度面积。
4.如权利要求3所述的土体三维孔隙率计算方法，其特征在于，所述通过获得的所述孔隙体积以及所述土体表面颗粒与孔隙整体体积计算获得孔隙率包括：
通过获得的所述孔隙体积以及所述土体表面颗粒与孔隙整体体积的比值获得所述孔隙率，计算关系式为
其中，n为所述孔隙率。
5.如权利要求1-4中任意一项所述的土体三维孔隙率计算方法，其特征在于；所述二值化阈值序列为梯度变化的二值化阈值序列，所述二值化阈值的变化梯度不大于10。
6.如权利要求5所述的土体三维孔隙率计算方法，其特征在于：所述二值化阈值的变化梯度为5或10。
7.如权利要求1所述的土体三维孔隙率计算方法，其特征在于，所述制备土样包括：
按照预设的土体孔隙率和压实度制备并获得原始土样；
切割所述原始土样，获得具有断面的多个半成土样；
分别清除所述半成土样断面上的浮动颗粒，获得多个成品土样。
8.如权利要求7所述的土体三维孔隙率计算方法，其特征在于，所述切割所述原始土样，获得具有断面的多个半成土样包括：
通过割土刀在所述原始土样中部切出沟槽；
沿所述沟槽手动掰断所述原始土样；
对所述断面进行打磨，获得所述半成土样。
9.如权利要求7或8所述的土体三维孔隙率计算方法，其特征在于，所述分别清除所述半成土样断面上的浮动颗粒包括：
利用胶带粘除所述半成土样断面上的浮动颗粒。
10.如权利要求7所述的土体三维孔隙率计算方法，其特征在于，所述分别清除所述半成土样断面上的浮动颗粒，获得多个成品土样之后还包括：
对所述成品土样进行干燥处理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711278075.6A CN108061697B (zh) | 2017-12-06 | 2017-12-06 | 土体三维孔隙率计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711278075.6A CN108061697B (zh) | 2017-12-06 | 2017-12-06 | 土体三维孔隙率计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108061697A CN108061697A (zh) | 2018-05-22 |
| CN108061697B true CN108061697B (zh) | 2020-09-25 |
Family
ID=62136230
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711278075.6A Active CN108061697B (zh) | 2017-12-06 | 2017-12-06 | 土体三维孔隙率计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108061697B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109507084B (zh) * | 2018-12-27 | 2021-06-29 | 辽宁工程技术大学 | 一种散体材料孔隙率智能识别真三轴实验系统及方法 |
| CN111982743B (zh) * | 2020-08-31 | 2022-06-24 | 长春工程学院 | 一种基于物联网的火山岩鉴定方法、系统、终端及介质 |
| CN113523185B (zh) * | 2021-06-24 | 2022-11-18 | 哈尔滨理工大学 | 一种喷墨3d打印砂型/砂芯致密度的检测方法 |
| CN113763351B (zh) * | 2021-09-06 | 2022-09-13 | 深圳大学 | 岩土体孔隙率测量方法、装置及智能终端 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5359194A (en) * | 1992-05-01 | 1994-10-25 | Texaco Inc. | X-ray CT measurement of secondary (vugular) porosity in reservoir core material |
| CN102841220A (zh) * | 2012-07-16 | 2012-12-26 | 天津城市建设学院 | 一种基于孔隙率的黏土扫描电镜照片图像分割方法 |
| CN104535475A (zh) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | 碳酸盐岩微观结构的确定方法及装置 |
| CN105372168A (zh) * | 2015-12-08 | 2016-03-02 | 清华大学 | 一种混凝土孔隙三维结构的构建方法 |
| CN106053316A (zh) * | 2016-06-08 | 2016-10-26 | 燕山大学 | 土体微观结构孔隙分布及分形程序与计算方法 |
| CN106198256A (zh) * | 2016-06-22 | 2016-12-07 | 西安科技大学 | 一种土体剪切带扩展过程观测方法 |
| CN106546521A (zh) * | 2016-10-12 | 2017-03-29 | 北京师范大学 | 一种基于ct扫描技术量化土壤大孔隙空间网络结构的方法 |
| CN106846315A (zh) * | 2017-02-09 | 2017-06-13 | 河海大学 | 一种基于土体图像计算土体裂隙率的方法 |
-
2017
- 2017-12-06 CN CN201711278075.6A patent/CN108061697B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5359194A (en) * | 1992-05-01 | 1994-10-25 | Texaco Inc. | X-ray CT measurement of secondary (vugular) porosity in reservoir core material |
| CN102841220A (zh) * | 2012-07-16 | 2012-12-26 | 天津城市建设学院 | 一种基于孔隙率的黏土扫描电镜照片图像分割方法 |
| CN104535475A (zh) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | 碳酸盐岩微观结构的确定方法及装置 |
| CN105372168A (zh) * | 2015-12-08 | 2016-03-02 | 清华大学 | 一种混凝土孔隙三维结构的构建方法 |
| CN106053316A (zh) * | 2016-06-08 | 2016-10-26 | 燕山大学 | 土体微观结构孔隙分布及分形程序与计算方法 |
| CN106198256A (zh) * | 2016-06-22 | 2016-12-07 | 西安科技大学 | 一种土体剪切带扩展过程观测方法 |
| CN106546521A (zh) * | 2016-10-12 | 2017-03-29 | 北京师范大学 | 一种基于ct扫描技术量化土壤大孔隙空间网络结构的方法 |
| CN106846315A (zh) * | 2017-02-09 | 2017-06-13 | 河海大学 | 一种基于土体图像计算土体裂隙率的方法 |
Non-Patent Citations (4)
| Title |
|---|
| 基于GIS的黏性土SEM图像三维可视化与孔隙度计算;王宝军 等;《岩土力学》;20080131;第29卷(第1期);第251-255页 * |
| 基于SEM和IPP测定软黏土接触面积的试验;徐日庆 等;《浙江大学学报（工学版）》;20150831;第49卷(第8期);第1417-1425页 * |
| 基于SEM图像的软土三维孔隙率计算及影响因素分析;徐日庆 等;《岩石力学与工程学报》;20150731;第34卷(第7期);第1497-1502页 * |
| 阈值选取对黏土微结构参数的影响;李顺群 等;《解放军理工大学学报（自然科学版）》;20110831;第12卷(第4期);第354-360页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108061697A (zh) | 2018-05-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108061697B (zh) | 土体三维孔隙率计算方法 | |
| CN105352873B (zh) | 页岩孔隙结构的表征方法 | |
| CN104535475B (zh) | 碳酸盐岩微观结构的确定方法及装置 | |
| CN110021030B (zh) | 一种岩土体材料数字图像的分割阈值确定方法 | |
| CN104751473B (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
| CN109444015A (zh) | 多视域、多尺度下的页岩储层微孔系统鉴定方法 | |
| CN109584286A (zh) | 一种基于广义回归神经网络的沥青路面构造深度计算方法 | |
| CN110443793B (zh) | 一种沥青混合料空隙分布均匀性评价方法 | |
| CN105479751A (zh) | 一种三维生物打印水凝胶支架的优化控制系统与方法 | |
| CN110223267A (zh) | 基于高度直方图分割的耐火砖深度缺陷的识别方法 | |
| CN105095581B (zh) | 一种铸造缩孔缺陷图像生成方法 | |
| CN110554054B (zh) | 一种机夹式可转位车刀的三维刃口崩缺检测系统 | |
| CN114565658B (zh) | 基于ct技术的孔隙尺寸计算方法及装置 | |
| CN112950659B (zh) | 一种用于压裂支撑剂显微放大的图像处理方法 | |
| CN112116583B (zh) | 一种基于sem图像处理的绝缘纸板老化判别检验方法 | |
| WO2024152526A1 (zh) | 一种基于三维特征分布的路面纹理均匀性评价方法 | |
| CN110264472B (zh) | 一种基于sd-oct的水果近皮下细胞无损成像方法 | |
| CN110490971B (zh) | 一种生物显微镜下细胞动态特征三维图像重建的方法 | |
| CN115078422A (zh) | 一种三维矿物模型建立的方法 | |
| CN112147680A (zh) | 岩石内部微裂缝提取和刻画方法及系统 | |
| CN110599477A (zh) | 一种植草混凝土有效孔隙率的检测方法及系统 | |
| CN110163826A (zh) | 一种黄土微结构图像的处理方法 | |
| CN103424085A (zh) | 一种物体表面三维形貌的测量方法 | |
| CN116452519A (zh) | 一种基于多模态信息融合的轮胎鼓包瑕疵检测方法 | |
| CN114283133A (zh) | 一种基于最邻近分析的沥青混合料均匀性评价方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |