CN115598023B - 一种确定松散破碎岩土体三维迂曲度方法 - Google Patents

一种确定松散破碎岩土体三维迂曲度方法 Download PDF

Info

Publication number
CN115598023B
CN115598023B CN202211389060.8A CN202211389060A CN115598023B CN 115598023 B CN115598023 B CN 115598023B CN 202211389060 A CN202211389060 A CN 202211389060A CN 115598023 B CN115598023 B CN 115598023B
Authority
CN
China
Prior art keywords
loose
broken rock
soil body
tortuosity
rock
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
CN202211389060.8A
Other languages
English (en)
Other versions
CN115598023A (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.)
Sichuan Zhongde Lujiang Equipment Technology Co ltd
Kunming University of Science and Technology
Original Assignee
Sichuan Zhongde Lujiang Equipment Technology Co ltd
Kunming University of Science and Technology
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 Sichuan Zhongde Lujiang Equipment Technology Co ltd, Kunming University of Science and Technology filed Critical Sichuan Zhongde Lujiang Equipment Technology Co ltd
Priority to CN202211389060.8A priority Critical patent/CN115598023B/zh
Publication of CN115598023A publication Critical patent/CN115598023A/zh
Application granted granted Critical
Publication of CN115598023B publication Critical patent/CN115598023B/zh
Priority to US18/499,251 priority patent/US11946845B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/02—Investigating particle size or size distribution
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Dispersion Chemistry (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种确定松散破碎岩土体三维迂曲度方法，包括：利用颗粒分析试验获得松散破碎岩土体的颗粒级配曲线，并求得等效粒径和平均粒径；利用含水率试验、密度试验与比重试验得到松散破碎岩土体的孔隙率；利用等效粒径、平均粒径和孔隙率，求得松散破碎岩土体三维迂曲度。通过上述方案，本发明具有逻辑简单、准确可靠、参数确定简单快捷等优点，在环境保护与生态修复技术领域具有很高的实用价值和推广价值。

Description

一种确定松散破碎岩土体三维迂曲度方法
技术领域
本发明涉及环境保护与生态修复技术领域，尤其是一种确定松散破碎岩土体三维迂曲度方法。
背景技术
迂曲度作为表征流体流经多孔介质孔隙通道的过程参量，可反映多孔介质孔隙通道的分布发育程度。目前，一些学者对多孔介质迂曲度开展了一定的研究，例如，郁伯铭等通过假定土体颗粒为正方形形状推导了迂曲度和孔隙率间的几何关系式，再如LIJianhua、袁培等基于Sierpinski carpet的自相似性推导了迂曲度与孔隙率呈分形维数的函数关系，MATYKA基于格子玻尔兹曼方法研究了多孔介质迂曲度-孔隙度关系，Lanfrey等在假定孔隙通道面积不变的情况下提出了一个关于孔隙率的迂曲度模型，李滔分析了微多孔介质迂曲度与孔隙结构的关系。颜瀚、戴邵衡、甘进等在假设土体颗粒为正方形颗粒及大小相等的条件下研究了考虑孔隙率和颗粒排列影响的土体迂曲度模型，黎章等在假设土体颗粒为圆形的情况下建立了考虑土体颗粒级配的几何迂曲度模型。薛东杰基于毛细管模型有与Hagen-Poiseuille 公式及达西定律建立微观的迂曲度表达式，饶登宇基于光滑粒子流体动力学（SPH）方法对重构三维土柱进行仿真试验及采用曲线拟合获取其扩散迂曲度值，AHMADI等通过考虑球形颗粒在三维空间内正方体和四面体排列方式与采用平均体积法建立了等粒径球体迂曲度模型。
另外，在专利技术中，如专利公开号为“CN109993786A”、名称为“迂曲度获取方法、装置、设备以及存储介质”的中国发明专利，其通过获取储集岩的孔隙网络模型，根据第一坐标数据和第二坐标数据确定第一孔隙和第二孔隙之间的直线距离；根据分布数据确定第一通道集合，并从第一通道集合筛选距离最短的第一通道作为第二通道，根据直线距离和第二通道的距离确定第一孔隙和第二孔隙之间的第一路径迂曲度，根据第一路径迂曲度集合确定储集岩的迂曲度。该技术以取储集岩为对象开展研究，通过图像处理得到其孔隙网络模型，在利用取储集岩的孔隙网络模型获得它的迂曲度。
再如专利公开号为“CN113297815A”、名称为“一种考虑多孔松散介质迂曲度的宾汉姆型水泥浆液渗透注浆扩散半径计算方法”的中国发明专利，又如专利公开号为“CN114720331A”、名称为“岩土体迂曲度幂律浆液时变性的柱半球渗透半径确定方法”的中国发明专利。上述两个发明专利技术仅考虑了多孔介质迂曲度影响的渗透注浆扩散半径确定方法且将迂曲度均认为一定值，而未提出迂曲度的具体计算方法。
综上所述，目前确定岩土体迂曲度的技术方法主要存在以下两个方面的问题：
第一，假设岩土体为正方形或圆形（球形）颗粒，这不能反应岩土体实际的颗粒分布特征；
第二，当前的迂曲度计算模型，要么包含的岩土体性质指标单一，如目前的绝大多数迂曲度计算模型仅考虑了岩土体孔隙率的影响，由此计算得到的迂曲度理论值与岩土体实际的迂曲度指标相差大；要么复杂包含的参数多且这些参数确定极其困难，不能在实践中应用来确定岩土体的迂曲度。
由此可见，目前的岩土体迂曲度计算方法难以满足工程实践的需求，不能合理有效地确定松散破碎岩土体实际的三维迂曲度指标。
因此，急需要提出一种逻辑简单、准确可靠的确定松散破碎岩土体三维迂曲度方法。
发明内容
针对上述问题，本发明的目的在于提供一种确定松散破碎岩土体三维迂曲度方法，本发明采用的技术方案如下：
一种确定松散破碎岩土体三维迂曲度方法，其包括以下步骤：
利用颗粒分析试验获得松散破碎岩土体的颗粒级配曲线，并求得等效粒径和平均粒径；
利用含水率试验、密度试验与比重试验得到松散破碎岩土体的孔隙率；
利用等效粒径、平均粒径和孔隙率，求得松散破碎岩土体三维迂曲度，其表达式为：
其中，为水的密度；为重力加速度；为水的黏度；表示以10为底的对数函数。
进一步地，所述确定松散破碎岩土体三维迂曲度方法，还包括：在松散破碎岩土体的颗粒级配曲线的纵坐标上选取20%、50%所在的点对应的横坐标值，分别记为等效粒径和平均粒径。
进一步地，联合采用筛析法与密度计法或移液管法测量获得松散破碎岩土体的颗粒级配曲线。
进一步地，所述松散破碎岩土体的孔隙率的表达式为：
其中，为4℃时纯蒸馏水密度；表示松散破碎岩土体密度；表示松散破碎岩土体质量含水量；表示松散破碎岩土体比重。
优选地，所述松散破碎岩土体质量含水量采用烘干法测量获得。
优选地，所述松散破碎岩土体比重采用比重瓶法与虹吸筒法结合测量获得。
与现有技术相比，本发明具有以下有益效果：
（1）本发明针对目前确定岩土体迂曲度方法的不足与实际工程中面临的技术难题，提出的一种确定松散破碎岩土体三维迂曲度的计算方法更能反映松散破碎岩土体实际的颗粒级配特性及其孔隙通道分布发育规律。
（2）本发明结合等效粒径、平均粒径和孔隙率求得的松散破碎岩土体三维迂曲度较目前的技术得到的迂曲度指标更符合实际情况，进而可为松散破碎岩土体的环境保护与生态修复工程设计提供理论指导与技术参考。
综上所述，本发明具有逻辑简单、准确可靠、参数确定简单快捷等优点，在环境保护与生态修复技术领域具有很高的实用价值和推广价值。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需使用的附图作简单介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对保护范围的限定，对于本领域技术人员来说，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本发明的逻辑流程图。
图2为本发明中实施例2-4的松散破碎岩土体颗粒级配曲线图。
图3为本发明中实施例5的松散破碎岩土体颗粒级配曲线图。
图4为本发明中实施例6的松散破碎岩土体颗粒级配曲线图。
图5为本发明中实施例7的松散破碎岩土体颗粒级配曲线图。
图6为本发明的松散破碎岩土体三维迂曲度与孔隙率的变化规律图。
图7为本发明的松散破碎岩土体三维迂曲度与等效粒径的变化规律图。
图8为本发明的松散破碎岩土体三维迂曲度与平均粒径的变化规律图。
具体实施方式
为使本申请的目的、技术方案和优点更为清楚，下面结合附图和实施例对本发明作进一步说明，本发明的实施方式包括但不限于下列实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。
本本实施例中，术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。本实施例的A、B、C、D和E无实际含义，由于松散破碎岩土体三维迂曲度的表达式较大，其不便于显示，因此，利用A、B、C、D和E进行拆分。
本实施例的说明书和权利要求书中的术语“第一”和“第二”等是用于区别不同的对象，而不是用于描述对象的特定顺序。例如，第一目标对象和第二目标对象等是用于区别不同的目标对象，而不是用于描述目标对象的特定顺序。
在本申请实施例中，“示例性的”或者“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性的”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言，使用“示例性的”或者“例如”等词旨在以具体方式呈现相关概念。
在本申请实施例的描述中，除非另有说明，“多个”的含义是指两个或两个以上。例如，多个处理单元是指两个或两个以上的处理单元；多个系统是指两个或两个以上的系统。
实施例1
如图1所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，其包括以下步骤：
第一步，联合采用筛析法与密度计法或移液管法测量获得获得松散破碎岩土体的颗粒级配曲线，并求得等效粒径和平均粒径，具体来说，在松散破碎岩土体的颗粒级配曲线的纵坐标上选取20%、50%所在的点对应的横坐标值，分别记为等效粒径和平均粒径。
第二步，利用含水率试验、密度试验与比重试验得到松散破碎岩土体的孔隙率。其中，松散破碎岩土体的孔隙率的表达式为：
其中，为4℃时纯蒸馏水密度；表示松散破碎岩土体密度；表示松散破碎岩土体质量含水量；表示松散破碎岩土体比重。松散破碎岩土体质量含水量采用烘干法测量获得。另外，松散破碎岩土体比重采用比重瓶法与虹吸筒法结合测量获得。
第三步，利用等效粒径、平均粒径和孔隙率，求得松散破碎岩土体三维迂曲度，其表达式为：
其中，为水的密度；为重力加速度；为水的黏度；表示以10为底的对数函数。
实施例2
如图2所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，具体步骤如下：
1）利用颗粒分析试验获得松散破碎岩土体颗粒级配曲线，并以此获得其等效粒径=0.97mm、平均粒径=6.83mm。
2）利用含水率试验、密度试验与比重试验得到松散破碎岩土体孔隙率=36.04%；
求得松散破碎岩土体的孔隙率；其中，取1000kg/m3；另外，分别由密度试验、含水率试验与比重试验测得的松散破碎岩土体密度=1830 kg/m3、质量含水量=6.37%、比重=2.69。
3）利用等效粒径、平均粒径和孔隙率，求得松散破碎岩土体三维迂曲度，取值为1000kg/m3；取9.8m/s2；水的黏度取1.52×10-3 Pa·s。通过上述计算，求得本实施例的松散破碎岩土体三维迂曲度=1.65。
实施例3
如图2所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，具体步骤如下：
1）利用颗粒分析试验获得松散破碎岩土体颗粒级配曲线，并以此获得其等效粒径=0.97mm、平均粒径=6.83mm；
2）利用含水率试验、密度试验与比重试验得到松散破碎岩土体孔隙率=43.13%；
松散破碎岩土体孔隙率中参数如下：松散破碎岩土体密度=1628 kg/m3、质量含水量=6.42%、比重=2.69。
3）求得松散破碎岩土体三维迂曲度，确定得到本实施例的松散破碎岩土体三维迂曲度=1.39。
实施例4
如图2所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，具体步骤如下：
1）利用颗粒分析试验获得松散破碎岩土体颗粒级配曲线，并以此获得其等效粒径=0.97mm、平均粒径=6.83mm；
2）利用含水率试验、密度试验与比重试验得到松散破碎岩土体孔隙率=32.73%；
松散破碎岩土体孔隙率中参数如下：松散破碎岩土体密度=1924 kg/m3、质量含水量=6.32%、比重=2.69。
3）确定得到本实施例的松散破碎岩土体三维迂曲度=1.78。
实施例5
如图3所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，具体步骤如下：
1）利用颗粒分析试验获得松散破碎岩土体颗粒级配曲线，并以此获得其等效粒径=0.83mm、平均粒径=3.68mm；
2）利用含水率试验、密度试验与比重试验得到松散破碎岩土体孔隙率=33.95%；
松散破碎岩土体孔隙率中参数如下：松散破碎岩土体密度=1770 kg/m3、质量含水量=2.67%、比重=2.61。
3）利用等效粒径、平均粒径和孔隙率，求得松散破碎岩土体三维迂曲度，取值为1000kg/m3；取9.8m/s2；水的黏度取1.01×10-3 Pa·s。通过上述计算，求得本实施例的松散破碎岩土体三维迂曲度=2.50。
实施例6
如图4所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，具体步骤如下：
1）利用颗粒分析试验获得松散破碎岩土体颗粒级配曲线，并以此获得其等效粒径=1.23mm、平均粒径=4.77mm；
2）利用含水率试验、密度试验与比重试验得到松散破碎岩土体孔隙率=33.71%；
松散破碎岩土体孔隙率中参数如下：松散破碎岩土体密度=1780 kg/m3、质量含水量=2.48%、比重=2.62。
3）求得本实施例的松散破碎岩土体三维迂曲度=1.69。
实施例7
如图5所示，本实施例提供了一种确定松散破碎岩土体三维迂曲度方法，具体步骤如下：
1）利用颗粒分析试验获得松散破碎岩土体颗粒级配曲线，并以此获得其等效粒径=1.53mm、平均粒径=5.33mm；
2）利用含水率试验、密度试验与比重试验得到松散破碎岩土体孔隙率=33.81%；
松散破碎岩土体孔隙率中参数如下：松散破碎岩土体密度=1780 kg/m3、质量含水量=2.64%、比重=2.62。
3）求得本实施例的松散破碎岩土体三维迂曲度=1.35。
从图6至图8中可以看出，松散破碎岩土体三维迂曲度与孔隙率、等效粒径和平均粒径之间均呈反向的变化关系，即松散破碎岩土体的孔隙率、等效粒径和平均粒径越大，其三维迂曲度越小。
上述实施例仅为本发明的优选实施例，并非对本发明保护范围的限制，但凡采用本发明的设计原理，以及在此基础上进行非创造性劳动而作出的变化，均应属于本发明的保护范围之内。

Claims (5)

1.一种确定松散破碎岩土体三维迂曲度方法，其特征在于，包括以下步骤：
利用颗粒分析试验获得松散破碎岩土体的颗粒级配曲线，并求得等效粒径d20和平均粒径d50；即在松散破碎岩土体的颗粒级配曲线的纵坐标上选取20％、50％所在的点对应的横坐标值，分别记为等效粒径d20和平均粒径d50；
利用含水率试验、密度试验与比重试验得到松散破碎岩土体的孔隙率φ；
利用等效粒径d20、平均粒径d50和孔隙率φ，求得松散破碎岩土体三维迂曲度τ，其表达式为：
Figure FDA0004057945040000011
Figure FDA0004057945040000012
Figure FDA0004057945040000013
Figure FDA0004057945040000014
Figure FDA0004057945040000015
Figure FDA0004057945040000016
其中，ρw为水的密度；g为重力加速度；μw为水的黏度；ln表示以10为底的对数函数。
2.根据权利要求1所述的一种确定松散破碎岩土体三维迂曲度方法，其特征在于，联合采用筛析法与密度计法或移液管法测量获得松散破碎岩土体的颗粒级配曲线。
3.根据权利要求1所述的一种确定松散破碎岩土体三维迂曲度方法，其特征在于，所述松散破碎岩土体的孔隙率φ的表达式为：
Figure FDA0004057945040000017
其中，
Figure FDA0004057945040000021
为4℃时纯蒸馏水密度；ρ表示松散破碎岩土体密度；ω表示松散破碎岩土体质量含水量；GS表示松散破碎岩土体比重。
4.根据权利要求3所述的一种确定松散破碎岩土体三维迂曲度方法，其特征在于，所述松散破碎岩土体质量含水量ω采用烘干法测量获得。
5.根据权利要求3所述的一种确定松散破碎岩土体三维迂曲度方法，其特征在于，所述松散破碎岩土体比重GS采用比重瓶法与虹吸筒法结合测量获得。
CN202211389060.8A 2022-11-08 2022-11-08 一种确定松散破碎岩土体三维迂曲度方法 Active CN115598023B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211389060.8A CN115598023B (zh) 2022-11-08 2022-11-08 一种确定松散破碎岩土体三维迂曲度方法
US18/499,251 US11946845B1 (en) 2022-11-08 2023-11-01 Method for determining three-dimensional tortuosity of loose and broken rock-soil mass

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211389060.8A CN115598023B (zh) 2022-11-08 2022-11-08 一种确定松散破碎岩土体三维迂曲度方法

Publications (2)

Publication Number Publication Date
CN115598023A CN115598023A (zh) 2023-01-13
CN115598023B true CN115598023B (zh) 2023-04-07

Family

ID=84852672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211389060.8A Active CN115598023B (zh) 2022-11-08 2022-11-08 一种确定松散破碎岩土体三维迂曲度方法

Country Status (2)

Country Link
US (1) US11946845B1 (zh)
CN (1) CN115598023B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109490166A (zh) * 2018-11-06 2019-03-19 北京理工大学 一种多孔介质任意横截面形状的通道渗透率确定方法
JP2020054985A (ja) * 2018-10-01 2020-04-09 株式会社デンソー 排ガス浄化フィルタ及び排ガス浄化フィルタの製造方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8881587B2 (en) * 2011-01-27 2014-11-11 Schlumberger Technology Corporation Gas sorption analysis of unconventional rock samples
WO2016007170A1 (en) * 2014-07-11 2016-01-14 Halliburton Energy Services, Inc. Imaging a porous rock sample using a nanoparticle suspension
KR101862807B1 (ko) * 2016-09-30 2018-05-31 한국가스공사 다공성 매질의 비틀림 수력 직경의 산출 방법 및 이를 이용한 다공성 매질 내의 유동 해석 방법
CN108008114A (zh) * 2017-11-30 2018-05-08 西南交通大学 一种粗粒土土体内部稳定性的判定方法
CN113189122B (zh) * 2021-05-13 2023-12-22 中海石油(中国)有限公司海南分公司 一种射孔损伤室内综合评价方法
CN113297815B (zh) * 2021-05-27 2022-12-30 昆明理工大学 松散介质迂曲度的宾汉浆液渗透注浆扩散半径计算方法
CN113984589B (zh) * 2021-11-01 2023-05-12 大连理工大学 一种计算岩石迂曲度和气体扩散系数的方法
CN114720331B (zh) * 2022-05-31 2022-09-02 成都理工大学 岩土体迂曲度幂律浆液时变性的柱半球渗透半径确定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020054985A (ja) * 2018-10-01 2020-04-09 株式会社デンソー 排ガス浄化フィルタ及び排ガス浄化フィルタの製造方法
CN109490166A (zh) * 2018-11-06 2019-03-19 北京理工大学 一种多孔介质任意横截面形状的通道渗透率确定方法

Also Published As

Publication number Publication date
CN115598023A (zh) 2023-01-13
US11946845B1 (en) 2024-04-02

Similar Documents

Publication Publication Date Title
CN111521532B (zh) 一种致密储层含气饱和度的测定方法
Vincens et al. Approaches to determine the constriction size distribution for understanding filtration phenomena in granular materials
CN111624147B (zh) 岩心的相对渗透率测定方法及装置
Xie et al. Fractal and multifractal analysis of carbonate pore-scale digital images of petroleum reservoirs
Gribble et al. Porometry, porosimetry, image analysis and void network modelling in the study of the pore-level properties of filters
CN104729972B (zh) 确定碳酸盐岩气水相对渗透率的方法及装置
Clerke Permeability, Relative Permeability, Microscopic Displacement Efficiency and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone
Dou et al. Pore structure, fractal characteristics and permeability prediction of tight sandstones: A case study from Yanchang Formation, Ordos Basin, China
CN107132171B (zh) 一种基于压汞-氮吸附联测数据确定致密储层孔径分布的方法
CN113051841B (zh) 一种强非均质页岩油储层表观渗透率确定方法和装置
AU2013289017A1 (en) Digital rock analysis systems and methods with reliable multiphase permeability determination
CN110927035A (zh) 一种低渗致密砂岩束缚水饱和度计算方法
CN114609010B (zh) 一种页岩储层油水相对渗透率的测定方法及装置
CN110702580B (zh) 一种基于信息熵的致密砂岩储层孔喉非均质性表征方法
CN109632604B (zh) 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法
Liu et al. Coupling lattice Boltzmann and finite volume methods for simulation of reactive transport
Mahoney et al. The effect of rank and lithotype on coal wettability and its application to coal relative permeability models
CN115618643A (zh) 宾汉浆液时水效应岩土体迂曲度柱半球渗透半径确定方法
CN115598023B (zh) 一种确定松散破碎岩土体三维迂曲度方法
Gerke et al. Studying structure and determining permeability of materials based on X-Ray microtomography data (using porous ceramics as an example)
Chu et al. Floc model and intrafloc flow
CN115526031A (zh) 一种反映矿物组成的多尺度孔隙网络模型构建及求解方法
Rahmati et al. Prediction of soil hydraulic conductivity at saturation using air permeability at any individual soil water content
Wulong et al. Pore-scale lattice Boltzmann modeling of soil water Distribution
Zhang et al. Characterization of Micro-Nano Pore Structures in Shale Oil Reservoirs and FIB-SEM 3D Reconstruction: A Case Study

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
CB02 Change of applicant information

Address after: 650000 No. 68, Wenchang Road, 121 Avenue, Kunming, Yunnan

Applicant after: Kunming University of Science and Technology

Applicant after: Sichuan Zhongde Lujiang Equipment Technology Co.,Ltd.

Address before: No. 1888, Jinggong East 1st Road, Xinxing street, Tianfu New District, Chengdu, Sichuan 610000

Applicant before: Sichuan Zhongde Lujiang Equipment Technology Co.,Ltd.

Applicant before: Kunming University of Science and Technology

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant