CN113075106A - 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法 - Google Patents

一种砾石土侵蚀阶段的水力参数及涌水量的确定方法 Download PDF

Info

Publication number
CN113075106A
CN113075106A CN202110326321.0A CN202110326321A CN113075106A CN 113075106 A CN113075106 A CN 113075106A CN 202110326321 A CN202110326321 A CN 202110326321A CN 113075106 A CN113075106 A CN 113075106A
Authority
CN
China
Prior art keywords
soil
particle
erosion
calculating
particle size
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
Application number
CN202110326321.0A
Other languages
English (en)
Other versions
CN113075106B (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.)
Guilin University of Technology
Original Assignee
Guilin University of 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 Guilin University of Technology filed Critical Guilin University of Technology
Priority to CN202110326321.0A priority Critical patent/CN113075106B/zh
Publication of CN113075106A publication Critical patent/CN113075106A/zh
Priority to US17/458,604 priority patent/US20220307964A1/en
Application granted granted Critical
Publication of CN113075106B publication Critical patent/CN113075106B/zh
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
    • G01N15/082—Investigating permeability by forcing a fluid through a sample
    • G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/04—Investigating osmotic effects
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • G01V9/02—Determining existence or flow of underground water

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Engineering & Computer Science (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Fluid Mechanics (AREA)
  • Geophysics (AREA)
  • Revetment (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明涉及一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，该方法包括以下步骤：根据分级侵蚀原理，计算各级粒径的土壤颗粒含量P和土壤孔隙率n，绘制各级粒径的PSD曲线及在每个侵蚀阶段下的PSD曲线簇；根据PSD曲线中的平均粒径Dj计算土壤颗粒的等效直径Dh，根据等效直径Dh计算土壤颗粒的最小等效孔隙直径d0；根据土壤颗粒含量P计算颗粒侵蚀的临界水力梯度icr；根据土壤孔隙率n及PSD曲线簇中的土壤颗粒的几何参数值，计算渗透系数kh；根据渗透系数kh计算渗流流速v及总渗流量Q，本发明通过计算总渗流量，反推出砾石土侵蚀程度和危险状况，以便采取相应措施予以治理和保护，进而避免事故发生。

Description

一种砾石土侵蚀阶段的水力参数及涌水量的确定方法
技术领域
本发明涉及地质岩土工程领域，特别是涉及一种砾石土侵蚀阶段的水力参数及涌水量的确定方法。
背景技术
在影响基坑开挖稳定性和岩溶塌陷安全性的诸多因素中，砾石土的浸水渗透侵蚀受到大家的广泛关注。特别是在岩溶地区分布的中间粒径缺失的双峰型砾石土、或宽级配砾石土，对岩土结构物防水防渗及岩溶塌陷的形成来说是一个重要的影响因素。
现有的技术中已经有很多针对砾石土浸水后因其颗粒分布特征客观上容易引发细颗粒流失、管涌等渗透侵蚀性问题的研究，不同于流土侵蚀，在管涌渗透侵蚀情况下及水力梯度较小的情况下，发生有关土体的渗透侵蚀及其稳定性问题，国外Sherard、Mace和Indraratna和Radampola的研究从渗透侵蚀的物理机理、实验手段到评价方法均有许多研究成果；国内则以刘杰和谢定松、毛昶熙等的研究具有代表性，特别是对于宽级配的砾石土分级评价模型的总结。
由于岩土体渗流量的直接影响基坑或坝体工程的突涌水量，因此，岩土体在不同侵蚀状态下，其水力特征参数的确定是尤为重要的，现有研究和相关规范很难定量评价和界定由渗透侵蚀直至突涌水发生的渗透稳定性问题，即涌水量的问题难以确定，从而导致突发的安全问题，因此，需要提供一种砾石土侵蚀阶段的水力参数及涌水量的确定方法。
发明内容
发明实施例提供了一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，可以解决现有研究和相关规范很难定量评价和界定由渗透侵蚀直至突涌水发生的渗透稳定性问题，即涌水量的问题难以确定，从而导致突发的安全问题。
本发明提供了一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，包括以下步骤：
根据分级侵蚀原理，计算各级粒径的土壤颗粒含量P和土壤孔隙率n，绘制各级粒径与各级粒径的土壤颗粒含量P的PSD曲线及在每个侵蚀阶段下各级粒径与各级粒径的土壤颗粒含量P的PSD曲线簇；
根据PSD曲线中的某两个粒径级区间内的平均粒径Dj计算土壤颗粒的等效直径Dh，根据等效直径Dh计算土壤颗粒的最小等效孔隙直径d0；
根据各级粒径的土壤颗粒含量P计算各级颗粒侵蚀的临界水力梯度icr；
根据土壤孔隙率n及PSD曲线簇中的土壤颗粒的几何参数值，计算渗透系数kh；
根据渗透系数kh计算渗流流速v及总渗流量Q。
优选地，各级粒径的土壤颗粒含量P和土壤孔隙率n的计算方法为：
根据下式(1)计算第j级粒径在(i+1)状态下的土壤颗粒含量Pj (i)及更新至(i+1)状态下的土壤孔隙率n(i+1)：
Figure BDA0002994787270000021
式中，Pj (i)为第j级粒径在(i)状态下的土壤颗粒含量；n i为在(i)状态下的土壤孔隙率；S为渗透侵蚀程度，即被冲走侵蚀的小于某粒径的土壤颗粒质量，其所占原土壤颗粒质量的百分比；Px为细粒含量，即土壤颗粒流失质量占总体土壤质量百分比。
优选地，根据等效直径Dh计算土壤颗粒的最小等效孔隙直径d0的方法为：
根据下式(2)计算土壤颗粒的等效直径Dh
Figure BDA0002994787270000031
根据下式(3)计算最小等效孔隙直径d0
Figure BDA0002994787270000032
式中，Dj是粒径级在j1到第j2之间的平均粒径，ΔSj为第j级粒径的重量与样本总重量之比，n是孔隙度，α是颗粒的形状系数。
优选地，计算各级颗粒侵蚀的临界水力梯度icr的方法为：
根据下式(4)计算临界水力梯度icr
Figure BDA0002994787270000033
式中：(icr)j为第j级颗粒侵蚀临界水力梯度，s为相对密度，即总体土壤密度相对于水体的密度；d85为土壤颗粒流失质量占总体土壤质量的85％的颗粒的粒径；dj是从土壤中被侵蚀的第j级的颗粒粒径，pj第j级粒径的颗粒含量。
优选地，计算渗透系数kh的方法为：
根据下式(5)计算渗透系数kh
Figure BDA0002994787270000034
式中，e为孔隙比，由土壤孔隙率计算得到e＝n/(1-n)，μw是水的动力粘滞系数，γw是水的重度，d10为土壤颗粒流失质量占总体土壤质量的10％的颗粒的粒径，Cu为不均匀系数。
优选地，根据渗透系数kh计算渗流流速v及总渗流量Q的方法为：
根据下式(6)计算土壤颗粒的渗流流速v
v＝Kh·icr (6)
根据下式(7)计算土壤颗粒的总渗流量Q
Q＝n·v·A (7)
式中：v是渗流流速，icr是临界水力梯度，Q是总渗流量，A是面积，n为土壤孔隙率。
相对于现有技术，本发明的优点在于：
本发明的方法是通过砾石土分级侵蚀情况下移动的PSD曲线，计算其动态几何参数及变化的临界水力梯度和渗透系数，再由达西公式计算渗流速度和涌水量，进而得出岩土水力特征参数及突涌水量的方法，使得在出现渗透侵蚀时，通过计算总渗流量，可反推出砾石土侵蚀程度和危险状况，以便采取相应措施予以治理和保护，进而避免事故发生，值得推广。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明方法的流程框图；
图2为本发明的三种不同土类的土壤颗粒粒径分布的PSD曲线；
图3为本发明土壤颗粒的各级粒径与各级粒径的土壤颗粒含量P的PSD曲线；
图4为本发明砾石土S2的在各个侵蚀阶段下的PSD曲线簇；
图5为本发明砾石土S3的在各个侵蚀阶段下的PSD曲线簇；
图6为本发明颗粒分级侵蚀过程中临界水力梯度计算框图；
图7为本发明砾石土临界水力梯度与各级土颗粒粒径之间的关系图；
图8是本发明砾石土临界水力梯度与颗粒百分含量之间的关系图；
图9是本发明不同侵蚀度下砾石土S2的几何参数及渗透系数的关系图；
图10是本发明不同侵蚀度下砾石土S3的几何参数及渗透系数的关系图；
图11是本发明在不同阶段的临界水力梯度icr低值区对渗流流速v和涌水量Q变化图；
图12是本发明在不同阶段的临界水力梯度icr高值区对渗流流速v和涌水量Q变化图。
具体实施方式
下面将结合本发明实施例中的附图1到附图11，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供了一种砾石土侵蚀阶段的水力参数及涌水量的确定方法：
在实施该方法前，需要根据含砾石土的主要建筑场地勘察资料，绘制出如图2所示的三种不同土类的土壤颗粒粒径分布的PSD曲线，从PSD曲线中获取各类砾石土的几何参数，几何参数包括d10、d30、d60、d70、d85、dq，d10表示为土壤颗粒流失质量占总体土壤质量的10％的颗粒的粒径；d30、d60、d70、d85与d10的意义相同，dq为两种砾石土类的粗细颗粒的分界直径，将几何参数代入
Figure BDA0002994787270000051
计算砾石土类的不均匀系数Cu，代入
Figure BDA0002994787270000052
计算曲率系数Cc，代入
Figure BDA0002994787270000053
计算dq，确定粗细颗粒分界直径dq，获得图2中对应的细粒含量Px。
表1为基于PSD曲线三种不同土类的几何参数
Figure BDA0002994787270000061
根据分级侵蚀原理，计算各级粒径的土壤颗粒含量P和土壤孔隙率n如图4和图5所示，根据每个状态下的土壤颗粒含量P绘制各级粒径分布的PSD曲线，两种砾石土分级侵蚀情况下右移动的PSD曲线簇，每条曲线对应的是细颗粒逐级被侵蚀冲走，土壤颗粒组成逐渐变得简化和单一，最终接近级配为卵砾石的新土类；比较两种砾石土的级配曲线(如图3)和移动的PSD曲线簇(图如图4和图5)可以看出，砾石土S2的初始状态(0％侵蚀度PSD曲线)相当于S3的第五级侵蚀状态(25％侵蚀度PSD曲线)情况。由此可以看出两种砾石土具有相似的沉积物质来源，前者的组分结构是后者在长期水力侵蚀作用之后所形成的阶段产物。
其中，根据下式(1)计算各级粒径的土壤颗粒含量P和土壤孔隙率n
Figure BDA0002994787270000062
如图3所示，根据每个状态下的土壤颗粒含量P绘制各级粒径分布的PSD曲线；从图中获取粒径级从j1到第j2之间的平均粒径Dj，根据平均粒径Dj确定土壤颗粒的等效直径Dh，根据等效直径Dh计算最小等效孔隙直径d0，
其中，根据下式(2)计算土壤颗粒的等效直径Dh
Figure BDA0002994787270000063
根据下式(3)计算最小等效孔隙直径d0
Figure BDA0002994787270000071
根据土壤颗粒含量P，计算各级颗粒侵蚀的临界水力梯度icr
其中，根据下式(4)计算各级颗粒侵蚀的临界水力梯度icr：
Figure BDA0002994787270000072
如图6所示，通过Pj(i)与Pj(i+1)之间交替循环迭代，同时判定每次迭代后是否满足继续侵蚀条件，从而实现不同侵蚀度下，一系列更新的PSD曲线簇和各级颗粒侵蚀临界水力梯度icr。
根据各级颗粒侵蚀的临界水力梯度icr绘制出如图7所示的两种砾石土临界水力梯度与各级土颗粒粒径之间的关系图，图7中，砾石土S2的临界水力梯度icr大于S3的临界水力梯度icr，这是由于其粒径范围小且粒径大，土颗粒的移动需要更大的渗透力，且砾石土S3在25％侵蚀度下，将其PSD曲线作为新土类计算获得的临界水力梯度icr随颗粒粒径变化情况，它与砾石土S2的初始临界水力梯度icr变化极为相似，Pearson相关系数达到90.4％。由此可进一步证实两种砾石土沉积物质具有同源性，前者是后者长期侵蚀下的阶段产物。
根据各级颗粒侵蚀的临界水力梯度icr绘制出如图8所示绘制出两种砾石土临界水力梯度与颗粒百分含量之间的关系图，图中砾石土S2中以砂粒作为可冲蚀细颗粒的主要成分；而S3中则以粉粒和少量粘粒作为可冲蚀细颗粒的主要成分，砂粒可冲蚀量偏少。
根据各土类几何参数值及所述最小等效孔隙直径d0，计算渗透系数kh。
通过式(2)和式(3)进行计算得到的最小等效孔隙直径d0及从PSD曲线直接读取方式获得几何参数，再代入式(5)获得砾石土的渗透系数kh。
Figure BDA0002994787270000073
如图9所示，不同侵蚀度下砾石土S2的几何参数及渗透系数图，其中的数据从移动的PSD曲线读取或由公式计算获得的几何参数取值，如粒径d10、d15可直接读取，而孔隙直径d0需通过式(3)计算获得，从图9中可见，,考虑孔隙度n值(0.20～0.47)影响，以充填区域表示孔隙度n值影响下几何参数和渗透系数的取值范围，其随着侵蚀度S的增加，S＝20％的时候，几何参数或渗透系数取值将趋于一致。
如图10所示，不同蚀度下砾石土S3的几何参数及渗透系数关系图，无论是取有效粒径d10还是d15，其渗透系数取值在侵蚀度S<15％情况下，计算结果与经典水力学公式的结果吻合，考虑了式(5)中不均匀系数Cu的影响，渗透性在数量级上变化范围较大，细颗粒的影响得以突显，如砾石土S3从10-6～10cm/s，能够涵盖从粉土到砾石的渗透系数范围。
根据达西定律和土壤渗流流量计算公式，计算渗流流速v及总渗流量Q；
其中，根据下式(6)计算土壤颗粒的渗流流速v
v＝Kh·icr (6)
根据下式(7)计算土壤颗粒的总渗流量Q
Q＝n·v·A (7)
根据图11和图12，得出如表2所示，在各个侵蚀阶段下砾石土的水力侵蚀性质；
在渗滤阶段(S＜5％)，临界水力梯度为icr＜0.01，渗透流速v小于5.0×10-4cm/s，可侵蚀的颗粒直径为0.01mm以下的粉细颗粒；
在细粒侵蚀阶段(5％≤S＜30％)，临界水力梯度icr为0.01～0.13，渗透流速v＝5.0×10-4～4.0×10-1cm/s；可侵蚀的颗粒直径为0.01～0.1mm的粉粒至细砂颗粒；
在粗粒侵蚀阶段(30％≤S≤40％)，临界水力梯度icr为0.13～0.50，渗透流速v＝0.4～2.8cm/s；可侵蚀的颗粒直径为0.1～0.57mm的细至中粗砂颗粒；
在突涌水或水淹阶段(S不考虑)，临界水力梯度icr为0.50～0.89，渗透流速v>2.8cm/s。
表2砾石土侵蚀各阶段参数及取值范围表
Figure BDA0002994787270000091
注：1：孔隙度n取0.47，即砾石土处于松散状态；A取1.0m2。
2：参考国内矿井突水流量大小，当流量小于50.0m3/h，即为小规模突水点。
综上，本发明的计算方法，通过砾石土分级侵蚀情况下移动的PSD曲线，计算其动态几何参数及变化的临界水力梯度和渗透系数，再由达西公式计算渗流速度和涌水量，进而得出岩土水力特征参数及突涌水量的方法，使得在出现渗透侵蚀时，通过计算总渗流量，可反推出砾石土侵蚀程度和危险状况，以便采取相应措施予以治理和保护，进而避免事故发生，值得推广。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (6)

1.一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，其特征在于，包括以下步骤：
根据分级侵蚀原理，计算各级粒径的土壤颗粒含量P和土壤孔隙率n，绘制各级粒径与各级粒径的土壤颗粒含量P的PSD曲线及在每个侵蚀阶段下各级粒径与各级粒径的土壤颗粒含量P的PSD曲线簇；
根据PSD曲线中的某两个粒径级区间内的平均粒径Dj计算土壤颗粒的等效直径Dh，根据等效直径Dh计算土壤颗粒的最小等效孔隙直径d0；
根据各级粒径的土壤颗粒含量P计算各级颗粒侵蚀的临界水力梯度icr；
根据土壤孔隙率n及PSD曲线簇中的土壤颗粒的几何参数值，计算渗透系数kh；
根据渗透系数kh计算渗流流速v及总渗流量Q。
2.如权利要求1所述的一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，其特征在于，计算各级粒径的土壤颗粒含量P和土壤孔隙率n的方法为：
根据下式(1)计算第j级粒径在(i+1)状态下的土壤颗粒含量Pj (i)及更新至(i+1)状态下的土壤孔隙率n(i+1)：
Figure FDA0002994787260000011
式中，Pj (i)为第j级粒径在(i)状态下的土壤颗粒含量；ni为在(i)状态下的土壤孔隙率；S为渗透侵蚀程度，即被冲走侵蚀的小于某粒径的土壤颗粒质量，其所占原土壤颗粒质量的百分比；Px为细粒含量，即土壤颗粒流失质量占总体土壤质量百分比。
3.如权利要求1所述的一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，其特征在于，根据等效直径Dh计算土壤颗粒的最小等效孔隙直径d0的方法为：
根据下式(2)计算土壤颗粒的等效直径Dh
Figure FDA0002994787260000021
根据下式(3)计算最小等效孔隙直径d0
Figure FDA0002994787260000022
式中，Dj是土壤颗粒粒径级在j1到第j2之间的平均粒径，ΔSj为第j级粒径的重量与样本总重量之比，n是孔隙度，α是颗粒的形状系数。
4.如权利要求1所述的一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，其特征在于，计算各级颗粒侵蚀的临界水力梯度icr的方法为：
根据下式(4)计算第j级颗粒侵蚀临界水力梯度(icr)j
Figure FDA0002994787260000023
式中：(icr)j为第j级颗粒侵蚀临界水力梯度，s为相对密度，即总体土壤密度相对于水体的密度；d85为土壤颗粒流失质量占总体土壤质量的85％的颗粒的粒径；dj是从土壤中被侵蚀的第j级的颗粒粒径，pj第j级粒径的颗粒含量。
5.如权利要求1所述的一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，其特征在于，计算渗透系数kh的方法为：
根据下式(5)计算渗透系数kh
Figure FDA0002994787260000024
式中，e为孔隙比，由土壤孔隙率计算得到e＝n/(1-n)，μw是水的动力粘滞系数，γw是水的重度，d10为土壤颗粒流失质量占总体土壤质量的10％的颗粒的粒径，Cu为不均匀系数。
6.如权利要求1所述的一种砾石土侵蚀阶段的水力参数及涌水量的确定方法，其特征在于，根据渗透系数kh计算渗流流速v及总渗流量Q的计算方法为：
根据下式(6)计算土壤颗粒的渗流流速v
v＝Kh·icr (6)
根据下式(7)计算土壤颗粒的总渗流量Q
Q＝n·v·A (7)
式中：v是渗流流速，icr是临界水力梯度，Q是总渗流量，A是面积，n为土壤孔隙率。
CN202110326321.0A 2021-03-26 2021-03-26 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法 Active CN113075106B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110326321.0A CN113075106B (zh) 2021-03-26 2021-03-26 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法
US17/458,604 US20220307964A1 (en) 2021-03-26 2021-08-27 Method for determining hydraulic parameters and water inflow in erosion stage of gravel soil

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110326321.0A CN113075106B (zh) 2021-03-26 2021-03-26 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法

Publications (2)

Publication Number Publication Date
CN113075106A true CN113075106A (zh) 2021-07-06
CN113075106B CN113075106B (zh) 2022-11-11

Family

ID=76610697

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110326321.0A Active CN113075106B (zh) 2021-03-26 2021-03-26 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法

Country Status (2)

Country Link
US (1) US20220307964A1 (zh)
CN (1) CN113075106B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116732979A (zh) * 2023-08-11 2023-09-12 中冶成都勘察研究总院有限公司 一种利用注浆钢管桩处理深厚红层回填土地基的方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115935710B (zh) * 2023-01-13 2023-07-07 水利部交通运输部国家能源局南京水利科学研究院 一种计算心墙堆石坝掺砾黏土心墙料渗透系数、评估渗流安全的方法
CN116930036B (zh) * 2023-07-24 2024-02-02 中国水利水电科学研究院 一种联合室内外确定原位临界水力坡降的方法
CN117607398B (zh) * 2024-01-23 2024-04-12 昆明理工大学 一种砾石土边坡失稳临界含水量的预测方法

Citations (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5269469A (en) * 1990-11-13 1993-12-14 Buhler Ag Method for measuring the fineness or bulk density, apparatus for carrying out the method and control system with such an apparatus
JP2006226000A (ja) * 2005-02-18 2006-08-31 Jdc Corp 改良土における固化材の配合方法および改良土における固化材と助剤の配合方法
CN102109447A (zh) * 2009-12-24 2011-06-29 上海张江中药现代制剂技术工程研究中心 快速判定中药提取物粉末吸湿性的方法
CN102608013A (zh) * 2012-03-02 2012-07-25 河海大学 一种测定管涌发展过程中孔隙率的方法
CN102930148A (zh) * 2012-10-22 2013-02-13 河海大学 一种基于随机起动的管涌渗透系数确定方法
CN103018424A (zh) * 2012-12-11 2013-04-03 重庆交通大学 管涌临界水力梯度与颗粒流失率室内联合测定装置及方法
WO2014078624A2 (en) * 2012-11-16 2014-05-22 Chevron U.S.A. Inc. Methods and systems for determining minimum porosity for presence of clathrates in sediment
CN104007045A (zh) * 2014-05-12 2014-08-27 河海大学 一种泥水盾构机泥浆成膜数值模拟方法
CN104021280A (zh) * 2014-05-19 2014-09-03 中冶集团武汉勘察研究院有限公司 一种适用于尾粉砂发生管涌的临界水力梯度的计算方法
CN104021277A (zh) * 2014-05-14 2014-09-03 河海大学 一种管涌现象的数值分析方法
CN104931401A (zh) * 2015-06-02 2015-09-23 中国科学院力学研究所 一种砂砾石土管涌侵蚀过程中渗透系数动态变化模型
CN105178259A (zh) * 2015-08-24 2015-12-23 中国科学院力学研究所 一种面板堆石坝垫层与过渡层的联合设计方法
CN107680131A (zh) * 2017-09-08 2018-02-09 北京理工大学 一种快速确定多孔介质表征单元体积尺寸的方法
CN107917865A (zh) * 2016-10-11 2018-04-17 中国石油化工股份有限公司 一种致密砂岩储层多参数渗透率预测方法
CN108008114A (zh) * 2017-11-30 2018-05-08 西南交通大学 一种粗粒土土体内部稳定性的判定方法
CN108693333A (zh) * 2018-06-14 2018-10-23 中铁二院成都勘察设计研究院有限责任公司 一种确定粗颗粒硫酸钠盐渍土盐胀系数的方法
CN109061105A (zh) * 2018-08-02 2018-12-21 中国水利水电科学研究院 一种土壤盐渍化临界地下水埋深的计算方法
CN110516322A (zh) * 2019-08-06 2019-11-29 湖北工业大学 一种不同水力梯度下黏土饱和非线性渗透系数的预测方法
CN110579427A (zh) * 2019-10-22 2019-12-17 桂林理工大学 一种裂隙-孔隙双渗透介质优势流模拟装置及实验方法
CN110598323A (zh) * 2019-09-12 2019-12-20 中国水利水电科学研究院 一种渗透破坏离散元模拟方法
CN111579454A (zh) * 2020-05-14 2020-08-25 上海交通大学 模拟砂土内细颗粒水平渗流侵蚀的试验装置及试验方法
CN111709148A (zh) * 2020-06-22 2020-09-25 河北工业大学 一种粘性砂土水力侵蚀破坏的离散元数值模拟方法
CN112098295A (zh) * 2020-09-11 2020-12-18 浙大城市学院 一种基于核磁共振技术测定无粘性土渗透系数的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110411916B (zh) * 2019-08-01 2021-07-20 国网四川省电力公司 一种巨粒土的颗粒级配测试方法

Patent Citations (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5269469A (en) * 1990-11-13 1993-12-14 Buhler Ag Method for measuring the fineness or bulk density, apparatus for carrying out the method and control system with such an apparatus
JP2006226000A (ja) * 2005-02-18 2006-08-31 Jdc Corp 改良土における固化材の配合方法および改良土における固化材と助剤の配合方法
CN102109447A (zh) * 2009-12-24 2011-06-29 上海张江中药现代制剂技术工程研究中心 快速判定中药提取物粉末吸湿性的方法
CN102608013A (zh) * 2012-03-02 2012-07-25 河海大学 一种测定管涌发展过程中孔隙率的方法
CN102930148A (zh) * 2012-10-22 2013-02-13 河海大学 一种基于随机起动的管涌渗透系数确定方法
WO2014078624A2 (en) * 2012-11-16 2014-05-22 Chevron U.S.A. Inc. Methods and systems for determining minimum porosity for presence of clathrates in sediment
CN103018424A (zh) * 2012-12-11 2013-04-03 重庆交通大学 管涌临界水力梯度与颗粒流失率室内联合测定装置及方法
CN104007045A (zh) * 2014-05-12 2014-08-27 河海大学 一种泥水盾构机泥浆成膜数值模拟方法
CN104021277A (zh) * 2014-05-14 2014-09-03 河海大学 一种管涌现象的数值分析方法
CN104021280A (zh) * 2014-05-19 2014-09-03 中冶集团武汉勘察研究院有限公司 一种适用于尾粉砂发生管涌的临界水力梯度的计算方法
CN104931401A (zh) * 2015-06-02 2015-09-23 中国科学院力学研究所 一种砂砾石土管涌侵蚀过程中渗透系数动态变化模型
CN105178259A (zh) * 2015-08-24 2015-12-23 中国科学院力学研究所 一种面板堆石坝垫层与过渡层的联合设计方法
CN107917865A (zh) * 2016-10-11 2018-04-17 中国石油化工股份有限公司 一种致密砂岩储层多参数渗透率预测方法
CN107680131A (zh) * 2017-09-08 2018-02-09 北京理工大学 一种快速确定多孔介质表征单元体积尺寸的方法
CN108008114A (zh) * 2017-11-30 2018-05-08 西南交通大学 一种粗粒土土体内部稳定性的判定方法
CN108693333A (zh) * 2018-06-14 2018-10-23 中铁二院成都勘察设计研究院有限责任公司 一种确定粗颗粒硫酸钠盐渍土盐胀系数的方法
CN109061105A (zh) * 2018-08-02 2018-12-21 中国水利水电科学研究院 一种土壤盐渍化临界地下水埋深的计算方法
CN110516322A (zh) * 2019-08-06 2019-11-29 湖北工业大学 一种不同水力梯度下黏土饱和非线性渗透系数的预测方法
CN110598323A (zh) * 2019-09-12 2019-12-20 中国水利水电科学研究院 一种渗透破坏离散元模拟方法
CN110579427A (zh) * 2019-10-22 2019-12-17 桂林理工大学 一种裂隙-孔隙双渗透介质优势流模拟装置及实验方法
CN111579454A (zh) * 2020-05-14 2020-08-25 上海交通大学 模拟砂土内细颗粒水平渗流侵蚀的试验装置及试验方法
CN111709148A (zh) * 2020-06-22 2020-09-25 河北工业大学 一种粘性砂土水力侵蚀破坏的离散元数值模拟方法
CN112098295A (zh) * 2020-09-11 2020-12-18 浙大城市学院 一种基于核磁共振技术测定无粘性土渗透系数的方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
BAIYU ZHU 等: "Coupled CFD-DEM simulation of granular LCM bridging in a fracture", 《PARTICULATE SCIENCE AND TECHNOLOGY》 *
BIN ZHU 等: "Analysis on water-inrush the process of deep excavation in karst area caused by soil internal erosion", 《IOP CONFERENCE SERIES: EARTH AND ENVIRONMENTAL SCIENCE》 *
DAN MA 等: "Experimental Investigation on Hydraulic Properties of Granular Sandstone and Mudstone Mixtures", 《GEOFLUIDS》 *
何修仁 等: "《注浆加固与堵水》", 31 March 1990, 东北工业学院出版社 *
俞茂宏: "《土力学新论》", 30 November 2019, 浙江大学出版社 *
刘汉东: "《岩土力学》", 30 November 2003, 中国广播电视大学出版社 *
徐萍: "鲁中山区小流域不同土地利用类型的土壤分形及水分入渗特征", 《中国水土保持科学》 *
房营光 等: "基于有效比表面积修正的Kozeny-Carman方程在黏土渗透中的适用性研究", 《岩土力学》 *
朱斌 等: "深部薄层煤岩体裂隙-孔隙双渗流模拟研究", 《中国矿业大学学报》 *
朱斌 等: "矿井断裂构造带滞后突水的流–固耦合模拟方法分析与滞后时间确定", 《岩石力学与工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116732979A (zh) * 2023-08-11 2023-09-12 中冶成都勘察研究总院有限公司 一种利用注浆钢管桩处理深厚红层回填土地基的方法
CN116732979B (zh) * 2023-08-11 2023-10-20 中冶成都勘察研究总院有限公司 一种利用注浆钢管桩处理深厚红层回填土地基的方法

Also Published As

Publication number Publication date
CN113075106B (zh) 2022-11-11
US20220307964A1 (en) 2022-09-29

Similar Documents

Publication Publication Date Title
CN113075106B (zh) 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法
Ahlinhan et al. Experimental investigation of critical hydraulic gradients for unstable soils
Cividini et al. Seepage-induced erosion in granular soil and consequent settlements
Estabragh et al. Models for predicting the seepage velocity and seepage force in a fiber reinforced silty soil
Jungkyu et al. Numerical modeling of sediment flushing from Lewis and Clark Lake
Chen et al. Runout analysis of slurry flows with Bingham model
Van Maren et al. Modelling hyperconcentrated flow in the Yellow River
Noseda et al. Upstream erosion and sediment passage at piano key weirs
CN104931401B (zh) 一种砂砾石土管涌侵蚀过程中渗透系数动态变化模型
Mahmood et al. The effect of soil type on matric suction and stability of unsaturated slope under uniform rainfall
CN106960105B (zh) 一种层理淤积河床上桥墩局部冲刷深度计算方法
CN102116014A (zh) 一种滑坡泥石流坝溃口特征测算方法及其应用
Cao et al. A large-strain vacuum-assisted radial consolidation model for dredged sludge considering lateral deformation
Ferrer‐Boix et al. Sorting of a sand–gravel mixture in a Gilbert‐type delta
CN109425561B (zh) 一种基于拦砂坝前淤积土壤剖面颗粒组成的已治理崩岗稳定性诊断技术
Sisinggih et al. Sediment transport functions in HEC-RAS 4.0 and their evaluation using data from sediment flushing of Wlingi reservoir-Indonesia
Chen et al. Prediction model investigations of the reinforced fine-grained tailings dam overtopping failure
Murray Erodibility of coarse sand-clayey silt mixtures
Arabi et al. Effects of fines content and relative density on erosion and recession of predominantly sandy beach–bluff system
Poulos et al. Density measurements in a hydraulic fill
Addis et al. Comparison of beaching slopes from two centrally discharging tailings storage facilities
Ellithy et al. Determination of erosion parameters of coarse-grained materials using a small flume
Mi et al. Roadway Settlement Characteristics on Mekong Delta for Various Landforms
Ellithy et al. Erosion Rate Equations for Coarse-Grained Materials Using a Small Flume Testing
Nordin et al. Numerical modelling of seepage analysis using SEEP/W: A case study for the Kerian River Flood Mitigation Project (Phase 3) in Bandar Baharu, Kedah

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