CN116448648B - 一种土体微观孔隙结构的显微观测区间推荐方法及系统 - Google Patents

一种土体微观孔隙结构的显微观测区间推荐方法及系统 Download PDF

Info

Publication number
CN116448648B
CN116448648B CN202310729329.0A CN202310729329A CN116448648B CN 116448648 B CN116448648 B CN 116448648B CN 202310729329 A CN202310729329 A CN 202310729329A CN 116448648 B CN116448648 B CN 116448648B
Authority
CN
China
Prior art keywords
microscopic
observation interval
soil body
pore structure
soil
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
CN202310729329.0A
Other languages
English (en)
Other versions
CN116448648A (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.)
Hainan Research Institute Of Zhejiang University
Original Assignee
Hainan Research Institute Of Zhejiang University
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 Hainan Research Institute Of Zhejiang University filed Critical Hainan Research Institute Of Zhejiang University
Priority to CN202310729329.0A priority Critical patent/CN116448648B/zh
Publication of CN116448648A publication Critical patent/CN116448648A/zh
Application granted granted Critical
Publication of CN116448648B publication Critical patent/CN116448648B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/22—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 measuring secondary emission from the material
    • G01N23/225—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 measuring secondary emission from the material using electron or ion
    • G01N23/2251—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 measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/40—Analysis of texture
    • G06T7/41—Analysis of texture based on statistical description of texture
    • G06T7/48—Analysis of texture based on statistical description of texture using fractals
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Dispersion Chemistry (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明涉及土力学数据量测领域，尤其是一种土体微观孔隙结构的显微观测区间推荐方法及系统。其中，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，包括如下步骤：设定显微倍数尺度；获取对应显微倍数下的待测土体样本图像；利用待测土体样本图像获取对应的分形维数；结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间。本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，通过可表征土体微观表面粗糙度的分形维数在不同显微倍数下的变化趋势，获得了能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的显微观测区间，实现了土体微观孔隙结构的显微观测区间的精准推荐。

Description

一种土体微观孔隙结构的显微观测区间推荐方法及系统
技术领域
本发明涉及土力学数据量测领域，尤其是一种土体微观孔隙结构的显微观测区间推荐方法及系统。
背景技术
土体微观结构是土体宏观存在的基础，这也决定了土体微观结构的研究的必然性和必要性。在土体微观结构中，尤其是膨润土这种整体结构由内部孔隙系统表征的土体，土体的微观孔隙结构更是研究重点。
当前，土体微观结构一般通过电子显微镜所获得的。但是，在电子显微镜不同的放大倍数下，电子显微镜捕获到的土体微观结构界面中涵盖的微观范围面积（体积）不同，所涵盖的土体微观颗粒的数量也不一样。
当显微镜放大倍数小时，显微可视界面内囊括的微观结构范围大，但土体微观结构的表现能力差；当显微镜放大倍数大时，显微可视界面内土体微观结构的表现能力虽然强，但显微可视界面内囊括的微观结构范围小，其捕获的数据不足以表征整体微观结构。
同时，土体微观结构会随着土体的压实密度、土体的含水量等因素变化；且土体微观结构的显微观测区间主要是根据测量工作人员的主观意志和经验选取的，当前并未建立能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的统一的显微观测区间。
发明内容
针对现有技术的不足，为建立能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的统一的显微观测区间，第一方面，本发明提供了一种土体微观孔隙结构的显微观测区间推荐方法，所述土体微观孔隙结构的显微观测区间推荐方法，包括如下步骤：设定显微倍数尺度；获取对应显微倍数下的待测土体样本图像；利用所述待测土体样本图像获取对应的分形维数；结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间。本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，通过可表征土体微观表面粗糙度的分形维数在不同显微倍数下的变化趋势，找到了能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的统一的显微观测区间，实现了土体微观孔隙结构的显微观测区间的精准推荐；同时，本发明的方法适用性强，为多种不同类型的土体提供了统一的可平衡土体微观结构范围大小和土体微观结构的表现能力的显微观测区间的推荐方法。
可选地，所述显微倍数尺度包括800倍、1500倍、2500倍、5000倍、10000倍、25000倍和50000倍。
可选地，所述获取对应显微倍数下的待测土体样本图像，包括如下步骤：提供环境扫描电子显微镜；将待测土体置入所述环境扫描电子显微镜的拍摄视野中；调整并固定所述环境扫描电子显微镜的拍摄视角；根据所述显微倍数尺度，利用所述环境扫描电子显微镜获得不同显微倍数下的待测土体样本图像。环境扫描电子显微镜可以设定不同的温度和水蒸气压力，进而可以容易地控制测量环境的相对湿度，即通过环境扫描电子显微镜可连续获得不同相对湿度下不同显微倍数的待测土体样本图像。
可选地，所述利用所述待测土体样本图像获取对应的分形维数，包括如下步骤：将所述待测土体样本图像二值化，获得黑白待测土体样本图像；设定多个孔隙面积；利用所述孔隙面积，在所述黑白待测土体样本图像中匹配对应的孔隙；分别汇总不同孔隙面积下匹配的孔隙总数量；结合多个孔隙面积以及对应的孔隙总数量，获得所述待测土体样本在对应显微倍数下的分形维数。
可选地，所述待测土体样本在对应显微倍数下的分形维数，满足如下公式：
，
其中，表示分形维数，/>表示孔隙面积，/>表示孔隙面积/>对应的孔隙总数量。
可选地，所述结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间，包括如下步骤：以显微倍数为横坐标，并以分型维数为纵坐标，构建坐标系；选取显微倍数以及对应的分形维数，在所述坐标系内拟合二次函数关系曲线图像；根据所述二次函数关系曲线图像的极值，获得所述极值处对应的显微倍数；利用所述极值处对应的显微倍数，推荐待测土体样本微观孔隙结构的显微观测区间。本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，利用待测土体样本在不同显微倍数下的分形维数，实现了能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的显微观测区间的推荐。
可选地，所述利用所述极值处对应的显微倍数，推荐待测土体样本微观孔隙结构的显微观测区间，包括如下步骤：根据所述极值处对应的显微倍数，构建显微观测区间推荐模型，所述显微观测区间推荐模型满足如下公式：
，/>；
其中，表示显微观测区间的最小显微倍数，/>表示显微观测区间的最大显微倍数，/>表示显微倍数与分型维数的二次函数关系曲线，/>表示显微倍数与分型维数的二次函数关系曲线的极值，/>表示显微观测裕度。本发明所提供的显微观测区间推荐模型通过设定显微观测裕度提升了显微观测区间的容错性，在保证为显微观测区间准确性的同时，进一步地为推荐更加精准的显微观测区间提供了起始显微倍数尺度范围。
可选地，所述结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间，还包括如下步骤：根据分形维数变化的单调性调整所述显微倍数尺度，直至随显微倍数增大，对应的分形维数变化不单调。这保证了用于拟合二次函数关系曲线图像数据的完整性和全面性。
可选地，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，还包括如下步骤：以所述显微观测区间为显微倍数尺度范围；并在所述显微倍数尺度范围设定显微倍数尺度。这进一步缩小了显微观测区间的范围，有助于推荐更加精准的显微观测区间。
第二方面，本发明还提供了一种土体微观孔隙结构的显微观测区间推荐系统，所述土体微观孔隙结构的显微观测区间推荐系统包括输入设备、处理器、存储器和输出设备，所述输入设备、所述处理器、所述存储器和所述输出设备相互连接，其中，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令，执行本发明第一方面所述的土体微观孔隙结构的显微观测区间推荐方法。本发明所述提供土体微观孔隙结构的显微观测区间推荐系统结构紧凑，运行稳定，并且具有极强的拓展性，有助于针对不同待测土体样本进行对应的系统改造，在扩大本发明的适用范围的同时，也保证了土体微观孔隙结构的显微观测区间推荐方法的稳定执行。
附图说明
图1为本发明的土体微观孔隙结构的显微观测区间推荐方法流程图；
图2为本发明实施例中步骤S02获取对应显微倍数下的待测土体样本图像流程图；
图3为本发明实施例中步骤S03利用待测土体样本图像获取对应的分形维数流程图；
图4为本发明实施例中为获取分形维数的待测土体样本图像的处理示意图；
图5为本发明实例中S04推荐待测土体样本微观孔隙结构的显微观测区间流程图；
图6为本发明实例中初始压实度为1.27g/cm3的钙基膨润土的二次函数关系曲线图像的过程示意图；
图7为本发明实例中初始压实度为1.90g/cm3的钙基膨润土的二次函数关系曲线图像的过程示意图；
图8为本发明的土体微观孔隙结构的显微观测区间推荐系统结构示意图。
实施方式
下面将详细描述本发明的具体实施例，应当注意，这里描述的实施例只用于举例说明，并不用于限制本发明。在以下描述中，为了提供对本发明的透彻理解，阐述了大量特定细节。然而，对于本领域普通技术人员显而易见的是：不必采用这些特定细节来实行本发明。在其他实例中，为了避免混淆本发明，未具体描述公知的电路，软件或方法。
在整个说明书中，对“一个实施例”、“实施例”、“一个示例”或“示例”的提及意味着：结合该实施例或示例描述的特定特征、结构或特性被包含在本发明至少一个实施例中。因此，在整个说明书的各个地方出现的短语“在一个实施例中”、“在实施例中”、“一个示例”或“示例”不一定都指同一实施例或示例。此外，可以以任何适当的组合和、或子组合将特定的特征、结构或特性组合在一个或多个实施例或示例中。此外，本领域普通技术人员应当理解，在此提供的示图都是为了说明的目的，并且示图不一定是按比例绘制的。
当前，土体微观结构一般通过电子显微镜所获得的，但是，在电子显微镜不同的放大倍数下，电子显微镜捕获到的土体微观结构界面中涵盖的微观范围面积（体积）不同，所涵盖的土体微观颗粒的数量也不一样。
请参见图1，为建立能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的显微观测区间，本发明提供了一种土体微观孔隙结构的显微观测区间推荐方法，所述土体微观孔隙结构的显微观测区间推荐方法，包括如下步骤：
S01、设定显微倍数尺度。
在一个可选的实施例中，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S01所述的所述显微倍数尺度包括800倍、1500倍、2500倍、5000倍、10000倍、25000倍和50000倍。
S02、获取对应显微倍数下的待测土体样本图像。
请参见图2，在一个可选的实施例中，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S02所述的获取对应显微倍数下的待测土体样本图像，包括如下步骤：
S021、提供环境扫描电子显微镜。
S022、将待测土体置入所述环境扫描电子显微镜的拍摄视野中。
S023、调整并固定所述环境扫描电子显微镜的拍摄视角。
S024、根据所述显微倍数尺度，利用所述环境扫描电子显微镜获得不同显微倍数下的待测土体样本图像。
在本实施例中，提供的是环境扫描电子显微镜(Environmental ScanningElectron Microscope，ESEM)，环境扫描电子显微镜可以设定不同的温度和水蒸气压力，进而可以容易地控制测量环境的相对湿度，即通过环境扫描电子显微镜可连续获得不同相对湿度下不同显微倍数的待测土体样本图像。同时，环境扫描电子显微镜可容易获得在不同显微倍数下孔隙直径范围，在本实施例中，待测土体样本为钙基膨润土，请参见表1，表1为环境扫描电子显微镜的显微倍数对应的孔隙直径：
本发明的通过显微倍数和待测土体样本的孔隙直径相对应，为利用后续获得的不同显微倍数下的分形维数构建立起分形维数与观测孔径尺寸的关系模型提供了数据基础。
S03、利用所述待测土体样本图像获取对应的分形维数。
请参见图3，在一个可选的实施例中，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S03所述的利用所述待测土体样本图像获取对应的分形维数，包括如下步骤：
S031、将所述待测土体样本图像二值化，获得黑白待测土体样本图像。
S032、设定多个孔隙面积。
S033、利用所述孔隙面积，在所述黑白待测土体样本图像中匹配对应的孔隙。
S034、分别汇总不同孔隙面积下匹配的孔隙总数量。
S035、结合多个孔隙面积以及对应的孔隙总数量，获得所述待测土体样本在对应显微倍数下的分形维数。
在本实施例中，请参见图4中的(a)图至(c)图：(a)图表示待测土体样本图像，该待测土体样本为钙基膨润土，该待测土体样本图像中尺度为10μm。(b)图表示黑白待测土体样本图像，即通过选取阈值将待测土体样本图像二值化，使得待测土体样本图像变成黑白待测土体样本图像，其中，图像二值化阈值可根据实际需求设定。(c)图表示利用孔隙面积匹配孔隙的过程示意图，(c)图中的线格表示设定的多个用于匹配孔隙的孔隙面积。
在本实施例中，所述待测土体样本在对应显微倍数下的分形维数，满足如下公式：
，
其中，表示分形维数，/>表示孔隙面积，/>表示孔隙面积/>对应的孔隙总数量。
在又一个可选的实施例中，通过本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S03获得的不同显示倍数下钙基片膨润土的分形维数如下表2所示：
由表2可知，随着显微倍数的增大，可观测到的孔隙直径单调减小，但对应分形维数并非单调变化。
S04、结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间。
请参见图5，在一个可选的实施例中，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S04所述的结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间，包括如下步骤：
S041、以显微倍数为横坐标，并以分型维数为纵坐标，构建坐标系。
S042、选取显微倍数以及对应的分形维数，在所述坐标系内拟合二次函数关系曲线图像。
步骤S042中显微倍数从小到大依次叠加选取的。在一个可选的实施例中，显微倍数为，且/>；其对应的分形维数分别。在本实施例中，首先选取/>与对应的/>拟合对应的二次函数关系曲线图像，若二次函数关系曲线图像中出现极值点，则停止拟合；若二次函数关系曲线图像中未出现极值点，则选取/>与对应的/>再次拟合对应的二次函数关系曲线图像，直至极值出现。
在又一个可选的实施例中，获取的初始压实度为1.27g/cm3的钙基膨润土的不同显示倍数下分形维数，如下表3所示：
根据表3中数据进行的拟合的二次函数关系曲线图像的过程如图6所示，图6中，横坐标为显微倍数，纵坐标为分型维数，菱形片表示拟合数据点。在本实施例中，显微倍数从小到大分别为800、1500、2500、5000、10000、25000以及50000；其对应的分形维数分别1.9004、1.8914、1.8623、1.8318、1.8287、1.8409以及1.8453。
在本实施例中，首先选取800、1500、2500的显微倍数与对应的1.9004、1.8914、1.8623的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中未出现极值点；再次选取800、1500、2500、5000的显微倍数与对应的1.9004、1.8914、1.8623、1.8318的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中未出现极值点；再次选取800、1500、2500、5000、10000的显微倍数与对应的1.9004、1.8914、1.8623、1.8318、1.8287的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中出现极值点，停止二次函数关系曲线图像的拟合。
在又一个可选的实施例中，获取的初始压实度为1.90g/cm3的钙基膨润土的不同显示倍数下分形维数，如下表4所示：
根据表3中数据进行的拟合的二次函数关系曲线图像的过程如图7所示，图7中，横坐标为显微倍数，纵坐标为分型维数，菱形片表示拟合数据点。在本实施例中，显微倍数从小到大分别为800、1500、2500、5000、10000以及50000；其对应的分形维数分别1.9099、1.9025、1.8903、1.8654、1.8277以及1.8330。
在本实施例中，首先选取800、1500、2500的显微倍数与对应的1.9099、1.9025、1.8903的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中未出现极值点；再次选取800、1500、2500、5000的显微倍数与对应的1.9099、1.9025、1.8903、1.8654的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中未出现极值点；再次选取800、1500、2500、5000、10000的显微倍数与对应的1.9099、1.9025、1.8903、1.8654、1.8277的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中未出现极值点；再次选取800、1500、2500、5000、10000、50000的显微倍数与对应的1.9099、1.9025、1.8903、1.8654、1.8277、1.8330的分形维数拟合对应的二次函数关系曲线图像，获得的二次函数关系曲线图像中出现极值点，停止二次函数关系曲线图像的拟合。
S043、根据所述二次函数关系曲线图像的极值，获得所述极值处对应的显微倍数。
在一个可选的实施例中，基于上述实施例中针对初始压实度为1.27g/cm3的钙基膨润土的二次函数关系曲线图像拟合状况，获得二次函数关系曲线图像的极值处对应的显微倍数为7500。
在又一个可选的实施例中，基于上述实施例中针对初始压实度为1.90g/cm3的钙基膨润土的二次函数关系曲线图像拟合状况，获得二次函数关系曲线图像的极值处对应的显微倍数为12500。
S044、利用所述极值处对应的显微倍数，推荐待测土体样本微观孔隙结构的显微观测区间。
在又一个可选的实施例中，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S044所述的利用所述极值处对应的显微倍数，推荐待测土体样本微观孔隙结构的显微观测区间，包括如下步骤：根据所述极值处对应的显微倍数，构建显微观测区间推荐模型，所述显微观测区间推荐模型满足如下公式：
，/>；
其中，表示显微观测区间的最小显微倍数，/>表示显微倍数与分型维数的二次函数关系曲线/>的极值/>，/>表示显微观测裕度，/>表示显微观测区间的最大显微倍数，/>表示显微倍数与分型维数的二次函数关系曲线，/>表示显微倍数与分型维数的二次函数关系曲线的极值。本发明所提供的显微观测区间推荐模型通过设定显微观测裕度提升了显微观测区间的容错性，在保证为显微观测区间准确性的同时，进一步地为推荐更加精准的显微观测区间提供了起始显微倍数尺度范围。
进一步地，在一个可选的实施例中，所述显微观测裕度为1000。因此，基于上述实施例中初始压实度为1.27g/cm3的钙基膨润土，推荐的显微观测区间为。基于上述实施例中初始压实度为1.90g/cm3的钙基膨润土，推荐的显微观测区间为/>。
本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S04通过子步骤S041至步骤S044针对不同的待测土体样本推荐了能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的显微观测区间。
进一步地，为保证用于拟合二次函数关系曲线图像数据的完整和全面性，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法中步骤S04所述的结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间，还包括如下步骤：根据分形维数变化的单调性调整所述显微倍数尺度，直至随显微倍数增大，对应的分形维数变化不单调。
具体地，请参见表2和表3：在表2中，显微倍数从800逐步增加到10000，对应的分形维数逐步减小；显微倍数增加到50000时，对应的分形维数逐步增加；即显微倍数尺度在800到50000时，随显微倍数增大显微，对应的分形维数变化不单调。在表3中，显微倍数从800逐步增加到10000，对应的分形维数逐步减小；显微倍数增加到25000时，对应的分形维数逐步增加；即显微倍数尺度在800到25000时，随显微倍数增大显微，对应的分形维数变化不单调。
应当理解，本发明所提供的所提供的土体微观孔隙结构的显微观测区间推荐方法可通过人工完成，也可通过C语言、C++、Java、VB、PHP等计算机编程语言编写对应的程序执行，其中，利用计算机编程语言编写对应的程序执行本发明的方法可提高本发明方法的执行效率，获得更高精准度的显微观测区间。本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，通过可表征土体微观表面粗糙度的分形维数在不同显微倍数下的变化趋势，找到了能够平衡显微可视界面内囊括的微观结构范围大小和土体微观结构的表现能力的显微观测区间，实现了土体微观孔隙结构的显微观测区间的精准推荐；同时，本发明的方法适用性强，为多种不同类型的土体提供了统一的可平衡土体微观结构范围大小和土体微观结构的表现能力的显微观测区间的推荐方法。
为进一步缩小显微观测区间的范围，以推荐更加精准的显微观测区间，本发明所提供的土体微观孔隙结构的显微观测区间推荐方法，还包括如下步骤：以所述显微观测区间为显微倍数尺度范围；并在所述显微倍数尺度范围设定显微倍数尺度。即通过本发明所获得的显微观测区间作为基础，再次利用本发明的方法进行更加精准的显微观测区间测定。
本发明还提供了一种土体微观孔隙结构的显微观测区间推荐系统，所述土体微观孔隙结构的显微观测区间推荐系统包括输入设备、处理器、存储器和输出设备，所述输入设备、所述处理器、所述存储器和所述输出设备相互连接，其中，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令，执行本发明提供的土体微观孔隙结构的显微观测区间推荐方法。本发明所述提供土体微观孔隙结构的显微观测区间推荐系统结构紧凑，运行稳定，并且具有极强的拓展性，有助于针对不同待测土体样本进行对应的系统改造，在扩大本发明的适用范围的同时，也保证了土体微观孔隙结构的显微观测区间推荐方法的稳定执行。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。

Claims (5)

1.一种土体微观孔隙结构的显微观测区间推荐方法，其特征在于，所述土体微观孔隙结构的显微观测区间推荐方法，包括如下步骤：
设定显微倍数尺度；
获取对应显微倍数下的待测土体样本图像；
利用所述待测土体样本图像获取对应的分形维数；
结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间；
所述利用所述待测土体样本图像获取对应的分形维数，包括如下步骤：
将所述待测土体样本图像二值化，获得黑白待测土体样本图像；
设定多个孔隙面积；
利用所述孔隙面积，在所述黑白待测土体样本图像中匹配对应的孔隙；
分别汇总不同孔隙面积下匹配的孔隙总数量；
结合多个孔隙面积以及对应的孔隙总数量，获得待测土体样本在对应显微倍数下的分形维数，所述分形维数满足如下公式：，其中，/>表示分形维数，/>表示孔隙面积，/>表示孔隙面积/>对应的孔隙总数量；
所述结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间，包括如下步骤：
以显微倍数为横坐标，并以分型维数为纵坐标，构建坐标系；
选取显微倍数以及对应的分形维数，在所述坐标系内拟合二次函数关系曲线图像；
根据所述二次函数关系曲线图像的极值，获得所述极值处对应的显微倍数；
利用所述极值处对应的显微倍数，推荐待测土体样本微观孔隙结构的显微观测区间，所述显微观测区间推荐模型满足如下公式：
，/>；
其中，表示显微观测区间的最小显微倍数，/>表示显微倍数与分型维数的二次函数关系曲线/>的极值/>，/>表示显微观测裕度，/>表示显微观测区间的最大显微倍数；
所述结合显微倍数和对应的分形维数，推荐待测土体样本微观孔隙结构的显微观测区间，还包括如下步骤：
根据分形维数变化的单调性调整所述显微倍数尺度，直至随显微倍数增大，对应的分形维数变化不单调。
2.根据权利要求1所述的土体微观孔隙结构的显微观测区间推荐方法，其特征在于，所述显微倍数尺度包括800倍、1500倍、2500倍、5000倍、10000倍、25000倍和50000倍。
3.根据权利要求1所述的土体微观孔隙结构的显微观测区间推荐方法，其特征在于，所述获取对应显微倍数下的待测土体样本图像，包括如下步骤：
提供环境扫描电子显微镜；
将待测土体置入所述环境扫描电子显微镜的拍摄视野中；
调整并固定所述环境扫描电子显微镜的拍摄视角；
根据所述显微倍数尺度，利用所述环境扫描电子显微镜获得不同显微倍数下的待测土体样本图像。
4.根据权利要求1-3任一所述的土体微观孔隙结构的显微观测区间推荐方法，其特征在于，还包括如下步骤：
以所述显微观测区间为显微倍数尺度范围；
并在所述显微倍数尺度范围设定显微倍数尺度。
5.一种土体微观孔隙结构的显微观测区间推荐系统，其特征在于，所述土体微观孔隙结构的显微观测区间推荐系统包括输入设备、处理器、存储器和输出设备，所述输入设备、所述处理器、所述存储器和所述输出设备相互连接，其中，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令，执行如权利要求1至4任一所述的土体微观孔隙结构的显微观测区间推荐方法。
CN202310729329.0A 2023-06-20 2023-06-20 一种土体微观孔隙结构的显微观测区间推荐方法及系统 Active CN116448648B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310729329.0A CN116448648B (zh) 2023-06-20 2023-06-20 一种土体微观孔隙结构的显微观测区间推荐方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310729329.0A CN116448648B (zh) 2023-06-20 2023-06-20 一种土体微观孔隙结构的显微观测区间推荐方法及系统

Publications (2)

Publication Number Publication Date
CN116448648A CN116448648A (zh) 2023-07-18
CN116448648B true CN116448648B (zh) 2023-08-11

Family

ID=87136042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310729329.0A Active CN116448648B (zh) 2023-06-20 2023-06-20 一种土体微观孔隙结构的显微观测区间推荐方法及系统

Country Status (1)

Country Link
CN (1) CN116448648B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4628468A (en) * 1984-04-13 1986-12-09 Exxon Production Research Co. Method and means for determining physical properties from measurements of microstructure in porous media
CN109444015A (zh) * 2018-10-31 2019-03-08 成都理工大学 多视域、多尺度下的页岩储层微孔系统鉴定方法
CN114018780A (zh) * 2021-10-11 2022-02-08 北京市科学技术研究院分析测试研究所(北京市理化分析测试中心) 一种陶瓷涂层孔隙和分形维数的测定方法及应用
CN115235977A (zh) * 2022-07-22 2022-10-25 中国石油大学(华东) 基于分形特征计算页岩多尺度孔径分布的方法及系统
CN115656481A (zh) * 2022-12-14 2023-01-31 海南浙江大学研究院 一种土体微观结构持水能力表征方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4628468A (en) * 1984-04-13 1986-12-09 Exxon Production Research Co. Method and means for determining physical properties from measurements of microstructure in porous media
CN109444015A (zh) * 2018-10-31 2019-03-08 成都理工大学 多视域、多尺度下的页岩储层微孔系统鉴定方法
CN114018780A (zh) * 2021-10-11 2022-02-08 北京市科学技术研究院分析测试研究所(北京市理化分析测试中心) 一种陶瓷涂层孔隙和分形维数的测定方法及应用
CN115235977A (zh) * 2022-07-22 2022-10-25 中国石油大学(华东) 基于分形特征计算页岩多尺度孔径分布的方法及系统
CN115656481A (zh) * 2022-12-14 2023-01-31 海南浙江大学研究院 一种土体微观结构持水能力表征方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Assessment of multi-scale pore structures and pore connectivity domains of marine shales by fractal dimensions and correlation lengths;Yuchen Fan 等;Fuel;第330卷;1-14 *

Also Published As

Publication number Publication date
CN116448648A (zh) 2023-07-18

Similar Documents

Publication Publication Date Title
Wargo et al. Comparison of focused ion beam versus nano-scale X-ray computed tomography for resolving 3-D microstructures of porous fuel cell materials
WO2019205579A1 (zh) 岩心微米ct孔隙结构重构方法
CA2573732A1 (en) Automated system of processing biological specimens and method
CN109269892B (zh) 一种聚合物弹性模量的测量方法及系统
CN115656481B (zh) 一种土体微观结构持水能力表征方法及系统
CN109523599B (zh) 一种透射电镜高分辨原子图像中矢量的标定方法及系统
CN110441209B (zh) 一种基于致密储层数字岩心计算岩石渗透率的方法
CN111551117A (zh) 显微图像焦点漂移距离的测量方法及系统、计算机设备
CN116448648B (zh) 一种土体微观孔隙结构的显微观测区间推荐方法及系统
Tsugawa et al. Elastic shell theory for plant cell wall stiffness reveals contributions of cell wall elasticity and turgor pressure in AFM measurement
CN109142189A (zh) 一种水泥基材料宏细观孔结构识别与评价方法
Chen et al. Intelligent adaptive sampling guided by Gaussian process inference
CN115146215A (zh) 一种基于数字岩心的微孔径数据多尺度拼接方法及系统
WO1993009396A1 (en) Method for creating a three-dimensional corporeal model from a very small original
CN108320307B (zh) 一种确定储层岩石样品的有效单元体积的方法及装置
CN114518070B (zh) 与扫描探针显微镜联用系统的快速定位方法、系统
CN112734143A (zh) 一种孔隙结构分形表征的方法、装置及系统
CN114778583B (zh) 扫描电镜目标点原位观测快速定位方法、装置及存储介质
Saghi et al. Electron tomography of regularly shaped nanostructures under non‐linear image acquisition
CN115165722A (zh) 表面改性板状元件实际环境应用的换热特性测试方法
CN117095041A (zh) 页岩孔隙结构表征方法、装置、存储介质及电子设备
Abuelfilat et al. Bridging structure and mechanics of three-dimensional porous hydrogel with X-ray ultramicroscopy and atomic force microscopy
CN114487116B (zh) 一种多孔材料声学性能测试系统及方法
Luna et al. Proposed Method for Determining Cross-Sectional Area of Guadua angustifolia Bamboo Fibers
CN116465815B (zh) 一种土体宏微观孔隙结构分形方法及系统

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