CN107389519A - 一种基于压汞法的土‑水特征曲线两特征参数预测方法 - Google Patents
一种基于压汞法的土‑水特征曲线两特征参数预测方法 Download PDFInfo
- Publication number
- CN107389519A CN107389519A CN201710418707.8A CN201710418707A CN107389519A CN 107389519 A CN107389519 A CN 107389519A CN 201710418707 A CN201710418707 A CN 201710418707A CN 107389519 A CN107389519 A CN 107389519A
- Authority
- CN
- China
- Prior art keywords
- aperture
- soil
- characteristic
- mrow
- msub
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 61
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 title claims abstract description 37
- 229910052753 mercury Inorganic materials 0.000 title claims abstract description 37
- 238000002347 injection Methods 0.000 title claims abstract description 27
- 239000007924 injection Substances 0.000 title claims abstract description 27
- 239000002689 soil Substances 0.000 claims abstract description 44
- 239000011148 porous material Substances 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims abstract description 18
- 238000009826 distribution Methods 0.000 claims abstract description 7
- 230000003993 interaction Effects 0.000 claims description 32
- 239000000243 solution Substances 0.000 claims description 9
- 239000007787 solid Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 239000011800 void material Substances 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 5
- 238000000643 oven drying Methods 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000002474 experimental method Methods 0.000 description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- HPNSNYBUADCFDR-UHFFFAOYSA-N chromafenozide Chemical compound CC1=CC(C)=CC(C(=O)N(NC(=O)C=2C(=C3CCCOC3=CC=2)C)C(C)(C)C)=C1 HPNSNYBUADCFDR-UHFFFAOYSA-N 0.000 description 6
- 230000009471 action Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 229920006300 shrink film Polymers 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 235000012431 wafers Nutrition 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004821 distillation Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
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
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)
- Sampling And Sample Adjustment (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明公开了一种基于压汞法的土‑水特征曲线两特征参数预测方法，其特征在于：根据孔隙分布数据求解两参数，分别是分维数和进气值，再基于分形模型预测土‑水特征曲线。在分维数预测时，本发明提出了“临界孔径”的概念，选取小于“临界孔径”的孔隙数据进行拟合可得出最优分维数。在进气值预测时，本发明提出了“最大特征孔径”的概念，基于孔隙率分形模型预测土体的“最大特征孔径”，然后利用“最大特征孔径”反算进气值。试验结果表明，通过两特征参数预测的土‑水特征曲线均与实测值吻合较好。本发明提供的基于压汞技术的土‑水特征曲线分形预测方法，对于非饱和土水力特性研究具有极其重要的意义。
Description
技术领域
本发明涉及到非饱和土水力特性研究技术领域，更加具体来说是 一种基于压汞法的土-水特征曲线两特征参数预测方法。
背景技术
土-水特征曲线反映的是土体基质吸力与含水量之间的关系，在 非饱和土的孔隙中不但充填有水，而且还有空气，水—气分界面(收 缩膜)具有表面张力，在非饱和土中，孔隙气压力与孔隙水压力不相 等，并且孔隙气压力大于孔隙水压力，收缩膜承受着大于水压力的空 气压力，这个压力差值称为基质吸力。
直接测量的方法只能给出特定状态下(如特定的孔隙结构)的土- 水特征曲线，而非饱和土土力学理论需要获得任意变形条件下的土- 水特征曲线。此外，直接测定土-水特征曲线试验耗时太长、工作量 大，特别是在区域尺度上进行实际问题研究时，这一方法更是受到限 制。因此，土-水特征曲线预测研究显得尤为重要。
土-水特征曲线受孔隙结构、矿物成分、颗粒结构、密实程度和 土的应力历史等众多因素影响。对于同一种土(假设矿物成分一定)， 土-水特征曲线变化主要由孔隙结构所控制，密实程度和应力历史等 是通过影响孔隙结构来影响土-水特征曲线。近些年来，从孔隙结构 预测土-水特征曲线已成为许多学者研究的焦点。
已有专利：一种基于轴平移测试的非饱和土土水特征曲线的修正 方法(201510786347.8)
该方法主要应用于基于轴平移技术测试的土水特征曲线，为了提 高试验结果的准确性，解决试验的某些误差。只是修正土水特征曲线， 并不能对土水特征曲线进行预测。
已有专利：一种基于亚甲蓝的细集料土水特征曲线预估方法(申 请号：201610126859.6)
该方法主要基于亚甲蓝的细集料土水特征曲线预估方法，所述方 法涉及部分经验模型，并没有相关理论分析。此方法仅能说明针对细 集料这种方法有可行性，针对土壤此方法是否具有可行性不得而知。
发明内容
本发明的目的在于克服上述背景技术的不足之处，而提出一种基 于压汞法的土-水特征曲线两特征参数预测方法。
本发明的目的是通过如下技术方案来实施的：一种基于压汞法的 土-水特征曲线两特征参数预测方法，包括如下步骤：
①、基于压汞法试验，获取微观孔隙分布数据；
②、利用分维数求解公式对压汞法试验的数据进行拟合处理得出 分维数，分维数求解公式如下：
(3-D)ln di∝ln(1-φ+V(≤di)) I
式中：V(≤di)表示单位土样总体积内包含的孔径小于等di的孔隙 总体积；φ为土样总孔隙率，单位为％；D为分维数；
③、基于孔隙率分形模型预测土体的“最大特征孔径”，所用模型表 达式如下：
φ＝1-(dmin/dmax)3-D； II
式中，φ表示土样总孔隙率，单位％；dmin为“最小特征孔径”； dmax为“最大特征孔径”；
④、根据预测的“最大特征孔径”计算进气值，计算公式如下：
式中：ψa为进气值，无量纲；Ts为表面张力，单位为N；α为接 触角，单位为°；dmax为最大特征孔径，单位为mm；参数λ<1；
⑤、利用分维数和进气值两特征参数预测结果，根据土水特征曲线模 型预测土-水特征曲线，模型表达式如下：
式中，wi为质量含水率；Gs表示土粒相对密度；ψa为进气值；ψa表示相应于孔径di的基质吸力；e分为孔隙比。
在上述技术方案中：在所述步骤2中；计算分维数D的步骤包括：
2.1)、利用维数求解公式I对压汞法所得数据进行拟合，以 ln(1-φ+V(≤di))纵坐标，ln di横坐标，绘制散点图；
2.2)、在散点图中，存在一个特殊孔径正处在“拐点”位置，在 孔径小于临界孔径之前的数据点，几乎呈直线分布；在孔径大于临界 孔径之后的数据点分布明显“上凸”，该特殊孔径称为“临界孔径”； 选取小于“临界孔径”的压汞法试验数据同步骤①进行拟合，得出分 维数D。
在上述技术方案中：在步骤③中；利用孔隙率分形模型预测土体 的“最大特征孔径”dmax，通过土的基本物理试验“烘干法”可求得 总孔隙率，取“最小特征孔径”dmin为0.2nm，通过总孔隙率利用公式 II反算出“最大特征孔径”dmax，再根据预测的“最大特征孔径”dmax利用公式III计算进气值ψa。
在上述技术方案中：在步骤⑤中；根据小于“临界孔径”拟合得 出的分维数D和取“最小特征孔径”dmin为0.2nm时算出的“最大特 征孔径”dmax最后得出的进气值ψa，通过两特征参数的预测值，利用 公式IV可以绘制出土-水特征曲线。
本发明包括如下优点：本发明可以获得任意变形条件下的土-水 特征曲线，为土-水特征曲线研究提供了一种新的预测方法。
附图说明
图1是本发明一具体实施中对不同干密度黏性土土-水特征曲线 试验结果。
图2是本发明一具体实施中对不同干密度黏性土压汞试验结果。
图3是本发明一具体实施中对全部压汞试验数据直接进行直线 拟合。
图4是本发明一具体实施中对压汞试验数据小于“临界孔径” 的孔径进行直线拟合。
图5对本发明中全部压汞试验数据拟合分维数的计算结果。
图6对小于“临界孔径”的数据拟合分维数的计算结果。
图7是本发明一具体实施中不同干密度下根据dmin取值为0.1nm、 0.2nm、0.3nm得出的“最大特征孔径”dmax及进气值ψa。
图8为当dmin＝0.1nm时，根据分维数D和预测进气值ψa得出的 土-水特征曲线预测曲线与实测图。
图9为当dmin＝0.2nm时，根据分维数D和预测进气值ψa得出的 土-水特征曲线预测曲线与实测图。
图10为当dmin＝0.3nm时，根据分维数D和预测进气值ψa得出的 土-水特征曲线预测曲线与实测图。
具体实施方式
下面结合附图详细说明本发明的实施情况，但它们并不构成对本 发明的限定，仅作举例而已；同时通过说明对本发明的优点将变得更 加清楚和容易理解。
参照图1-10所示：本发明实现通过分维数和进气值的预测，更 准确的预测土-水特征曲线，预测了一种预测值和实测值几乎重合的 曲线。
本发明预测土-水特征曲线方法分为3种，如图8-10所示，对土 -水特征曲线的预测值和实测值进行比较，得出最优预测土-水特征曲 线方法。
下面详述本发明得具体实施过程：
本发明所涉及实验的过程如下，土样为武汉地区黏性土。土样晾 晒后碾碎，并过2mm筛，利用液压千斤顶制备2组平行试样，每组 试样包含干密度ρd为1.30、1.35、1.40、1.45、1.50、1.60、1.71g/cm3的7个试样。2组试样均抽真空饱和，饱和后的试样一组用于压力板 仪试验，得出试样实测的土-水特征曲线，如图1所示；另一组以备 压汞法试验之用，为获取压汞数据预测土-水特征曲线，用于预测土 水特征曲线，如图2所示。
压力板仪试验
本次发明所使用的压力板仪试验为选用美国Soilmoisture公司 生产的压力板仪，其主要由压力室、高进气值陶瓷板、压力表和氮气 源等组成。试验前先将所制环刀土样抽气饱和，放置于饱和HAE陶瓷 板上，使两者紧密贴合。由氮气源提供气压，通过调压阀调节到设定 气压，容器气压增加到某一预定值，孔隙水可逐渐排出土样，直至达 到平衡状态。
试验过程中，要时刻监测并记录水的流出情况(即称水瓶的质量)， 当水的流出不在变化是，当存水瓶中水的质量不再变化后，认为达到 平衡状态，即可以结束此压力下的试验。然后打开容器，测量试样的 含水率。最后根据某一压力平衡下试样的质量含水率，压力即为基质 吸力作为横坐标，质量含水量作为纵坐标拟合绘制土-水特征曲线。 压汞法试验
本次发明所使用的压力板仪试验采用的仪器为美国康塔公司生 产的PoreMaster33压汞仪，试样真空饱和达到要求后，将试样进行 液氮极速冷冻，使土中的液体成为不具有膨胀性的非结晶态冰，然后 利用冷干机抽真空，使土中非结晶态冰冷升华。
用推土器将预处理过的试样从环刀里小心的取出，使用细钢丝锯将试 样切取1cm3的小方块，进行压汞试验.实际操作过程应尽量避免对 试样的破坏扰动，保证试验结果的可靠性.控制进汞压力和进汞速度， 将汞压入试样孔隙中，使汞充充满孔隙，记录每一级进汞压力时的进 汞体积，利用压力和孔半径的关系，将进汞压力换算成孔隙孔径，得 到试样的孔隙孔径分布数据。
预测分维数推导公式：
对获取的压汞法试验利用下列公式直接进行直线拟合，以 ln(1-φ+V(≤di))纵坐标，ln di横坐标，绘制散点图，如图3所示。
根据文献(表征孔隙及颗粒体积与尺度分布的两类岩土体分形模 型)建议的表征孔隙孔径分形特性的数学表达式可用V式表示：
式中V(>di)为孔径大于di的孔隙总体积，Va表示土样总体积，L2表示 研究区域尺寸，D为分维数。
根据公式V，有：(di)3-D∝Va-V(>di) VI
现假设分析的土样总体积为单位1，即令Va＝1，若土样总孔隙率 为φ，则有V(>di)＝φ-V(≤di)，其中V(≤di)表示单位土样总体积内包 含的孔径小于等于di的孔隙总体积。那么公式VI可变为：
(di)3-D∝1-φ+V(≤di) VII
对公式VII两边同时取对数，则有：
(3-D)ln di∝ln(1-φ+V(≤di)) I
式中：V(≤di)表示单位土样总体积内包含的孔径小于等di的孔隙 总体积；φ为土样总孔隙率，单位为％；D为分维数；
根据公式I便可以求出分维数D。
对以上直线拟合图分析发现途中存在一个“特殊孔径”位于“拐 点位置”，称为“临界孔径”，同样对小于“临界孔径”的孔隙进行直 线拟合，获取相关拟合效果图，如图4所示。
根据对小于“临界孔径”的孔隙进行直线拟合得出的图形，拟合 图中相关系数更高，直线效果更加明显，分形特征更加突出，根据公 式I进行计算，得出分维数。
进气值推导公式：
利用Young-Laplace理论认为基质吸力与曲率半径之间存在反比 的关系，可以提出了一种基于孔隙分布预测土-水特征曲线的方法， 可用公式VIII表示：
式中：di表示孔隙等效直径，wi表示相应于孔径di的基质吸力， Ts为表面张力，α为接触角，温度一定时4Tscosα为常数。
引入参数λ来考虑该尺度效应，即压汞法试验获取的某等效连 通孔径为di，则相应于土-水特征曲线试样的等效连通孔径为λdi， 这里λ<1。
则提出了一种新的预测方法，如公式III式所示：
式中：ψa为进气值，无量纲；Ts为表面张力，单位为N；α为接 触角，单位为°；dmax为最大特征孔径，单位为mm；参数λ<1；
土水特征曲线模型公式推导：
根据文献(表征孔隙孔径分布的岩土体孔隙率模型及其应用)建 议了一种孔隙率模型，该模型表达式如下：
φ(>di)＝1-(di/dmax)3-D IX
其中，D是分维数，dmax为最大特征孔径，φ(>di)是孔径大于di的孔隙率。用dmin为最小特征孔径，取di＝dmin，则有：
φ＝1-(dmin/dmax)3-D II
在公式II中，假设非饱和土中孔径小于等于di的孔隙充满了水， 其体积含水率为θ，总孔隙率为φ，，则θi＝φ(≤di)＝φφ(>di)， 代入到公式IX，则有：
θi＝φ-1+(di/dmax)3-D X
结合公式III式和公式X可推导出：
θi＝φ-1+(ψa/ψi)3-D XI
式中ψa为进气值，其值按公式III计算：
其中公式III是体积含水量表示的土-水特征曲线模型，根据体 积含水量与质量含水量之间的关系θ＝wi·ρd，干密度其中 质量含水量wi表示的数学模型公式为IV表示：
式中，wi为质量含水率；Gs表示土粒相对密度；ψa为进气值；ψa表示相应于孔径di的基质吸力；e分为孔隙比。
式中Gs、e均为已知值，分别表示土粒相对密度及孔隙比均为已 知值。值得说明的时，公式IX成立的条件是ψ≤ψa，当ψi≤ψa时，统一近似取值饱和质量含水量，wi＝e/Gs。
根据最优分维数和进气值两特征参数所预测的土-水特征曲线， 以及实测的土-水特征曲线拟合对比，(如图9所示)可以很明显得出 的土-水特征曲线预测值和实测值几乎重合。
最大特征孔径计算方法：
对进气值ψa进行预测，关键求出“最大特征孔径”dmax即所有孔 隙中最大的孔隙直径，对“最大特征孔径”的进行计算，通过总孔隙 率预测“最大特征孔径”的求法，试样的总孔隙率是已知的，利用“最 小特征孔径”dmin即烘干法中所能烘干的最小孔隙直径，取为0.1、0.2、 0.3nm。利用公式如下：
φ＝1-(dmin/dmax)3-D； II
式中，φ表示土样总孔隙率，单位％；dmin为“最小特征孔径”；dmax为“最大特征孔径”。
最小特征孔径取值问题：
总孔隙率可以通过“烘干法”实验得到，“烘干法”是在100-105℃ 烘干至恒温。也就是说，最小特征孔径dmin应该考虑对应100-105℃ 所能烘干的最小孔隙直径。土中水分为液态水、固态水及气态水。
气态水质量在这里忽略不计，所以土中水的质量这里主要考虑液 态水和固态水。
液态水包括：自由水和结合水。其中，自由水包括重力水和毛细 水，重力水在本身重力作用下可以运动，毛细水受表面张力作用，力 的作用范围1nm-104nm；结合水受分子引力作用，力的作用范围约为 0.3nm-102nm。在100-105℃恒温烘干时，液态水(自由水和结合水) 几乎会全部挥发，则最小特征孔径dmin≤0.3nm。对于固态水认为受 氢键和共价键作用，力的作用范围分别为0.2nm-0.3nm，0.1nm-0.2nm。 对于固态水，100-105℃恒温可能使部分水分挥发，但仍有部分水分 残留。也就是说，所考虑的最小特征孔径范0.1nm≤dmin≤0.3nm。预 测最大特征孔径时，发明中采用dmin＝0.1nm、dmin＝0.2nm、dmin＝0.3nm 分别进行试算。
基于小于“临界孔径”的数据拟合分维数的计算结果得出的分 维数，根据dmin取值为0.1nm、0.2nm、0.3nm，将分为f、g、h三种 方法预测的土-水特征曲线，预测和实测图形如下所示(参照图8-10 所示)。分析不同干密度试样质量含水量预测值与实测值的关系，发现相应于三种方法的平均相关系数分别为0.972、0.988、0.987，平 均绝对误差分别为0.017、0.008、0.009。可见三种方法中，g法相 关系数最高，达到0.988；平均绝对误差最小，仅为0.008。因此，g 法预测效果最佳，建议采用g法(dmin取值为0.2nm)进行预测。
上述未详细说明的部分均为现有技术。
Claims (4)
1.一种基于压汞法的土-水特征曲线两特征参数预测方法，其特征在于：包括如下步骤：
①、基于压汞法试验，获取微观孔隙分布数据；
②、利用分维数求解公式对压汞法试验的数据进行拟合处理得出分维数，分维数求解公式如下：
(3-D)lndi∝ln(1-φ+V(≤di))I
式中：V(≤di)表示单位土样总体积内包含的孔径小于等di的孔隙总体积；φ为土样总孔隙率，单位为％；D为分维数；
③、基于孔隙率分形模型预测土体的“最大特征孔径”，所用模型表达式如下：
φ＝1-(dmin/dmax)3-D； II
式中，φ表示土样总孔隙率，单位％；dmin为“最小特征孔径”；dmax为“最大特征孔径”；
④、根据预测的“最大特征孔径”计算进气值，计算公式如下：
<mrow>
<msub>
<mi>&psi;</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mrow>
<msub>
<mi>&lambda;d</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mi>I</mi>
<mi>I</mi>
<mi>I</mi>
</mrow>
式中：ψa为进气值，无量纲；Ts为表面张力，单位为N；α为接触角，单位为°；dmax为最大特征孔径，单位为mm；参数λ<1；
⑤、利用分维数和进气值两特征参数预测结果，根据土水特征曲线模型预测土-水特征曲线，模型表达式如下：
<mrow>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>e</mi>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>a</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&psi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>-</mo>
<mi>D</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>G</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mi>I</mi>
<mi>V</mi>
</mrow>
式中，wi为质量含水率；Gs表示土粒相对密度；ψa为进气值；ψa表示相应于孔径di的基质吸力；e分为孔隙比。
2.根据权利要求1所述的一种基于压汞法的土-水特征曲线两特征参数预测方法，其特征在于：在所述步骤2中；计算分维数D的步骤包括：
2.1)、利用分维数求解公式I对压汞法所得数据进行拟合，以ln(1-φ+V(≤di))为纵坐标，lndi为横坐标，绘制散点图；
2.2)、在散点图中，存在一个特殊孔径正处在“拐点”位置，在孔径小于临界孔径之前的数据点，几乎呈直线分布；在孔径大于临界孔径之后的数据点分布明显“上凸”，该特殊孔径称为“临界孔径”；选取小于“临界孔径”的压汞法试验数据同步骤①进行拟合，得出分维数D。
3.根据权利要求1或2所述的一种基于压汞法的土-水特征曲线两特征参数预测方法，其特征在于：在步骤③中；利用孔隙率分形模型预测土体的“最大特征孔径”dmax，通过土的基本物理试验“烘干法”可求得总孔隙率，取“最小特征孔径”dmin取0.2nm，通过总孔隙率利用公式II反算出“最大特征孔径”dmax。
4.根据权利要求3所述的一种基于压汞法的土-水特征曲线两特征参数预测方法，其特征在于：在步骤⑤中；根据小于“临界孔径”拟合得出的分维数D和取“最小特征孔径”dmin为0.2nm时算出的“最大特征孔径”dmax最后得出的进气值ψa，通过两特征参数的预测值，利用公式IV可以绘制出土-水特征曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710418707.8A CN107389519B (zh) | 2017-06-06 | 2017-06-06 | 一种基于压汞法的土-水特征曲线两特征参数预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710418707.8A CN107389519B (zh) | 2017-06-06 | 2017-06-06 | 一种基于压汞法的土-水特征曲线两特征参数预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107389519A true CN107389519A (zh) | 2017-11-24 |
| CN107389519B CN107389519B (zh) | 2019-11-08 |
Family
ID=60333199
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710418707.8A Expired - Fee Related CN107389519B (zh) | 2017-06-06 | 2017-06-06 | 一种基于压汞法的土-水特征曲线两特征参数预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107389519B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108426820A (zh) * | 2018-04-12 | 2018-08-21 | 东南大学 | 一种基于微观孔隙结构的多孔介质的渗透系数预测方法 |
| CN108444895A (zh) * | 2018-06-14 | 2018-08-24 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
| CN109238911A (zh) * | 2018-09-10 | 2019-01-18 | 湖北工业大学 | 一种基于土中水分蒸发曲线预测土-水特征曲线vg模型的方法 |
| CN111141650A (zh) * | 2019-12-11 | 2020-05-12 | 中国地质大学（武汉） | 一种利用冷冻高压压汞表征页岩中水的赋存状态的方法 |
| CN114002409A (zh) * | 2021-11-02 | 2022-02-01 | 宁波大学 | 一种非饱和压实土过渡区土水特征曲线的快速测定方法 |
| CN114329700A (zh) * | 2021-12-17 | 2022-04-12 | 湖北工业大学 | 一种新型的预应力钢束实验参数计算方法 |
| CN114813441A (zh) * | 2022-04-21 | 2022-07-29 | 兰州大学 | 一种基于轴平移的高精度土水特征曲线测量方法 |
| CN116482001A (zh) * | 2023-06-20 | 2023-07-25 | 浙江大学海南研究院 | 一种土体宏微观孔隙分界孔径的确定及影响评估方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104697915A (zh) * | 2015-03-20 | 2015-06-10 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
| CN105136648A (zh) * | 2015-10-19 | 2015-12-09 | 江西理工大学 | 土壤有效孔径及其分布参数的测试方法 |
| CN106769688A (zh) * | 2017-01-04 | 2017-05-31 | 中国石油天然气股份有限公司 | 一种预测地质条件下致密砂岩中天然气扩散系数的方法 |
-
2017
- 2017-06-06 CN CN201710418707.8A patent/CN107389519B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104697915A (zh) * | 2015-03-20 | 2015-06-10 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
| CN105136648A (zh) * | 2015-10-19 | 2015-12-09 | 江西理工大学 | 土壤有效孔径及其分布参数的测试方法 |
| CN106769688A (zh) * | 2017-01-04 | 2017-05-31 | 中国石油天然气股份有限公司 | 一种预测地质条件下致密砂岩中天然气扩散系数的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 庄心善等: "压汞试验技术研究黏性土微观孔隙分布特性", 《湖北工业大学学报》 * |
| 陶高梁等: "土-水特征曲线的分形特性及其分析拟合", 《岩土力学》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108426820A (zh) * | 2018-04-12 | 2018-08-21 | 东南大学 | 一种基于微观孔隙结构的多孔介质的渗透系数预测方法 |
| CN108444895A (zh) * | 2018-06-14 | 2018-08-24 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
| CN108444895B (zh) * | 2018-06-14 | 2020-08-25 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
| CN109238911A (zh) * | 2018-09-10 | 2019-01-18 | 湖北工业大学 | 一种基于土中水分蒸发曲线预测土-水特征曲线vg模型的方法 |
| CN111141650A (zh) * | 2019-12-11 | 2020-05-12 | 中国地质大学（武汉） | 一种利用冷冻高压压汞表征页岩中水的赋存状态的方法 |
| CN114002409A (zh) * | 2021-11-02 | 2022-02-01 | 宁波大学 | 一种非饱和压实土过渡区土水特征曲线的快速测定方法 |
| CN114329700A (zh) * | 2021-12-17 | 2022-04-12 | 湖北工业大学 | 一种新型的预应力钢束实验参数计算方法 |
| CN114329700B (zh) * | 2021-12-17 | 2024-06-07 | 湖北工业大学 | 一种预应力钢束实验参数计算方法 |
| CN114813441A (zh) * | 2022-04-21 | 2022-07-29 | 兰州大学 | 一种基于轴平移的高精度土水特征曲线测量方法 |
| CN114813441B (zh) * | 2022-04-21 | 2024-05-10 | 兰州大学 | 一种基于轴平移的高精度土水特征曲线测量方法 |
| CN116482001A (zh) * | 2023-06-20 | 2023-07-25 | 浙江大学海南研究院 | 一种土体宏微观孔隙分界孔径的确定及影响评估方法 |
| CN116482001B (zh) * | 2023-06-20 | 2023-09-01 | 浙江大学海南研究院 | 一种土体宏微观孔隙分界孔径的确定及影响评估方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107389519B (zh) | 2019-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107389519A (zh) | 一种基于压汞法的土‑水特征曲线两特征参数预测方法 | |
| CN107402175B (zh) | 一种变形条件下土-水特征曲线的分形预测方法 | |
| Flint et al. | 2.2 Particle Density | |
| Campbell et al. | Application of a dew point method to obtain the soil water characteristic | |
| Likos et al. | Pore-scale model for water retention and fluid partitioning of partially saturated granular soil | |
| Yan et al. | Some issues on nanoindentation method to measure the elastic modulus of particles in composites | |
| CN107709964B (zh) | 用于评价储集岩润湿性的核磁共振气体等温线技术 | |
| Salager et al. | An innovative device for determining the soil water retention curve under high suction at different temperatures | |
| US8256268B2 (en) | System and method for measuring porosity of high strength and high performance concrete using a vacuum-pressure saturation method | |
| Vanstreels et al. | Effect of pore structure of nanometer scale porous films on the measured elastic modulus | |
| Khoury et al. | Effect of cyclic suction history (hydraulic hysteresis) on resilient modulus of unsaturated fine-grained soil | |
| Fratta et al. | Combined TDR and P-wave velocity measurements for the determination of in situ soil density—experimental study | |
| EP2954306A1 (en) | Apparatus and methodology for measuring properties of microporous material at multiple scales | |
| Li et al. | Calculation of capillary rise height of soils by SWCC model | |
| Liu et al. | Capillary rise method for the measurement of the contact angle of soils | |
| CN107389241B (zh) | 一种用激光冲击强化孔壁的应力峰值检测方法 | |
| CN109815628B (zh) | 快速确定土水特征曲线和动态回弹模量模型参数的方法 | |
| Azmi et al. | Estimation of Soil Water Characteristic Curves (SWCC) of mining sand using soil suction modelling | |
| CN102735584A (zh) | 液体浮力称量法测定钢砂铝中铝铁含量的方法 | |
| Chen et al. | Influence of degree of compaction on unsaturated hydraulic properties of a compacted completely decomposed granite | |
| Wiese et al. | Early-stage interaction between settlement and temperature-gradient metamorphism | |
| Darnbrough et al. | Elastic and plastic mechanical properties of lithium measured by nanoindentation | |
| Liu et al. | Influence of contact angle on soil–water characteristic curve with modified capillary rise method | |
| Alekseev et al. | Aerodynamic method for obtaining the soil water retention curve | |
| Wang et al. | Investigating icing behavior in cementitious material during freeze–thaw using low-temperature low-field NMR |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191108 |