CN112816387B - 一种渗透系数的确定方法、装置及存储介质 - Google Patents
一种渗透系数的确定方法、装置及存储介质 Download PDFInfo
- Publication number
- CN112816387B CN112816387B CN202011625598.5A CN202011625598A CN112816387B CN 112816387 B CN112816387 B CN 112816387B CN 202011625598 A CN202011625598 A CN 202011625598A CN 112816387 B CN112816387 B CN 112816387B
- Authority
- CN
- China
- Prior art keywords
- soil sample
- particle
- test soil
- particle size
- determining
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000003860 storage Methods 0.000 title claims abstract description 17
- 239000002245 particle Substances 0.000 claims abstract description 331
- 239000002689 soil Substances 0.000 claims abstract description 255
- 238000012360 testing method Methods 0.000 claims abstract description 187
- 238000004458 analytical method Methods 0.000 claims abstract description 73
- 239000010419 fine particle Substances 0.000 claims abstract description 31
- 239000011362 coarse particle Substances 0.000 claims abstract description 18
- 230000035515 penetration Effects 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 25
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000005553 drilling Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005029 sieve analysis Methods 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000002349 well water Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 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
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)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明实施例提供一种坡洪积混合土渗透系数的确定方法、装置及存储介质,属于岩土工程领域。所述确定方法包括:获取所述坡洪积混合土的测试土样;基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径;对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径;以及根据所述有效粒径与所述渗透参数之间的关联关系,确定所述测试土样的渗透系数。本发明实施例坡洪积渗透系数确定方法简单易行,更适用于深层地表的坡洪积混合土渗透系数的确定。
Description
技术领域
本发明涉及岩土工程领域,具体涉及一种渗透系数的确定方法、装置及存储介质。
背景技术
混合土是由细粒土和粗粒土混杂形成且缺乏中间粒径的土,其主要由黏粒、粉粒、砾粒和漂砾粒组成,成因主要为洪积、坡积和残积。目前,普遍认为有效粒径对混合土渗透性有较大影响,但对于不均匀系数、平均粒径等参数对混合土渗透性的影响研究较少。
混合土颗粒组成具有较大的分散性特点,其渗透系数受到多方面因素影响,其中颗粒级配是影响渗透系数的重要因素。细粒的存在改变了土的孔隙特性,而孔隙特性是影响土渗透性能的最重要因素。显然,分析细粒对混合土的孔隙特性的影响,建立细粒等状态参量和混合土渗透系数之间的规律性关系,进而揭示混合土渗透性能的细粒效应意义十分明确。
目前对混合土渗透性进行计算的方法中,通常采用现场试验的方法。例如,浅地表的混合土一般采用单环注水试验,地下深处的混合土一般采用单井注水试验或抽水试验来确定,根据现场获得的数据得到渗透系数。但是坡洪积混合土夹杂大粒径块石,存在钻井困难、清水钻进容易塌孔、钻孔进行泥浆护壁后会形成泥皮等问题,故通过现场试验很难获得准确地的试验数据。
发明内容
本公开实施例提供一种坡洪积混合土渗透系数的确定方法、装置及存储介质,用于解决存在的相关技术问题。
根据本发明实施例的第一方面,提供一种坡洪积混合土渗透系数的确定方法,所述确定方法包括:获取所述坡洪积混合土的测试土样;基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径;对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径;以及根据所述有效粒径与所述渗透参数之间的关联关系,确定所述测试土样的渗透系数。
优选的,所述测试土样包括深层地表土样和/或浅层地表土样。
优选的,所述基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径,包括:获取所述测试土样的颗粒分布,以得到所述测试土样的第一颗粒分析曲线;基于所述第一颗粒分析曲线,确定用于所述测试土样中的粗颗粒和细颗粒的区分粒径。
优选的,所述获取所述测试土样的颗粒分布,以得到所述测试土样的第一颗粒分析曲线,包括:基于预设的多个颗分粒径,对所述测试土样进行颗分试验,得到用于表征所述测试土样颗粒分布的第一颗粒数据;基于所述第一颗粒数据,得到所述测试土样的第一颗粒分析曲线。
优选的,所述基于所述第一颗粒分析曲线,确定所述测试土样的区分粒径,包括:分析所述第一颗粒分析曲线的平缓曲线段,以得到所述平缓曲线段在所述第一颗粒分析曲线中与相邻曲线段的第一临界点和第二临界点;获取所述第一临界点和第二临界点所分别对应的第一粒径和第二粒径,所述区分粒径大于等于所述第一粒径且小于等于所述第二粒径。
优选的,所述对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径,包括:基于所述区分粒径,分析所述测试土样中小于或者等于所述区分粒径下的颗粒分布,以得到所述测试土样的第二颗粒分析曲线;基于所述第二颗粒分析曲线,确定所述测试土样的有效粒径。
优选的,所述基于所述区分粒径,分析所述测试土样中小于或者等于所述区分粒径下的颗粒分布,以得到所述测试土样的第二颗粒分析曲线,包括:基于所述区分粒径,将所述测试土样中小于或者等于所述区分粒径的颗粒保留,得到所述测试土样的第二颗粒数据;基于所述第二颗粒数据,得到所述测试土样的第二颗粒分析曲线。
优选的,基于所述有效粒径,采用如下公式计算得到所述测试土样的渗透系数:
k=Cd10(a+bt)
其中,k为渗透系数,C、a、b为常数系数,d10为有效粒径,t为水温。
根据本发明实施例的第二方面,提供一种坡洪积混合土渗透系数的确定装置,所述确定装置包括:获取单元,获取所述坡洪积混合土的测试土样;第一确定单元,基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径;第二确定单元,对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径;以及第三确定单元,根据所述有效粒径与所述渗透参数之间的关联关系,确定所述测试土样的渗透系数。
根据本发明实施例的第三方面,提供一种电子设备,该电子设备包括:至少一个处理器;存储器;其中,所述存储器存储有能被所述至少一个处理器执行的指令,所述至少一个处理器通过执行所述存储器存储的指令实现本发明实施例第一方面或者第一方面任意一项所述的确定方法。
根据本发明实施例的第四方面,提供一种机器可读存储介质,该机器可读存储介质上存储有指令,该指令用于使得机器执行第一方面或者第一方面任意一项所述的确定方法。
通过上述技术方案,获取坡洪积混合土的测试土样,确定测试土样的区分粒径,基于区分粒径确定有效粒径,从而最终得到渗透系数。该方法相比现有技术方案简便易行,尤其对于深层地表的混合土渗透系数计算更为准确。
本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施例的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明实施例,但并不构成对本发明实施例的限制。在附图中:
图1是根据一示例性实施例示出的一种坡洪积混合土渗透数的确定方法流程图;
图2是根据一示例性实施例示出的一种区分粒径确定方法流程图;
图3是根据一示例性实施例示出的一种获取测试土样的颗粒分布方法流程图;
图4是根据一示例性实施例示出的一种第一颗粒分析曲线示意图;
图5是根据一示例性实施例示出的一种第一颗粒曲线分析方法流程图;
图6是根据一示例性实施例示出的确定有效粒径方法流程图;
图7是根据一示例性实施例示出的得到第二颗粒分析曲线方法流程图;
图8是根据一示例性实施例示出的一种第二颗粒分析曲线示意图;
图9是根据一示例性实施例示出的一种坡洪积混合土渗透系数的确定装置框图。
具体实施方式
以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明实施例,并不用于限制本发明实施例。
渗透系数是反映土壤渗透能力的一个重要指标,尤其渗透系数的确定对工程现场施工有着非常重要的意义。而目前,对于坡洪积混合土的渗透系数确定,通常采用现场试验方法,其中深层地表的坡洪积混合土一般采用单孔注水试验或者抽水试验来确定其渗透系数,但是由于坡洪积混合土存在较大粒径块石,造成在现场试验中钻进困难,容易塌孔,而钻孔进行泥浆护壁厚形成泥皮,难以获得准确的试验数据。
鉴于上述存在的技术问题,本公开提供了一种坡洪积混合土渗透系数的确定方法,如图1所示,确定方法分为步骤S11至步骤S14:
步骤S11,获取坡洪积混合土的测试土样。
本公开优选实施例采用室内试验方法,对于坡洪积混合土渗透系数的确定,需要先获取该测试土样,从而对获取到的测试土样进行在室内进行相应的试验分析。
步骤S12,基于测试土样,确定用于区分测试土样中的粗颗粒和细颗粒的区分粒径。
坡洪积混合土的颗粒多由粗颗粒和细颗粒组成,细颗粒的存在改变了坡洪积混合土的土孔隙特性,而土孔隙特性是影响坡洪积混合土渗透性能的最重要因素。本公开实施例中采用的确定方法,考虑细颗粒对混合土土孔隙特性的影响,建立细颗粒含量和坡洪积混合土之间的关联关系。因此,需要基于获取到的测试土样,确定测试土样的区分粒径,该区分粒径介于粗颗粒的粒径和细颗粒的粒径之间,通过确定的区分粒径可以将测试土样中的粗颗粒和细颗粒加以区分,从而可以针对测试土样的细颗粒进行进一步的分析坡洪积混合土中细颗粒含量对于渗透系数的影响。
步骤S13,对测试土样中粒径小于区分粒径的颗粒进行分析,并根据分析结果确定与测试土样的渗透参数相关联的有效粒径。
本公开实施例考虑到细颗粒对于渗透系数的影响,当步骤S12确定出区分粗颗粒和细颗粒的区分粒径,进而根据该区分粒径,分析小于区分粒径的颗粒,即对于测试土样中的细颗粒的颗粒含量进行进一步分析,进而得到与测试土样的渗透参数相关联的有效粒径。
步骤S14,根据有效粒径与渗透参数之间的关联关系,确定测试土样的渗透系数。
通过步骤S13得到的有效粒径,根据有效粒径与渗透参数之间的关联关系建立的计算公式,得到测试土样的渗透系数。在此,可选的计算公式将在下文结合示例描述。
在步骤S11至S14中,获取坡洪积混合土的测试土样,确定测试土样的区分粒径,基于区分粒径确定有效粒径,从而最终得到渗透系数。该方法相比现有技术方案简便易行,尤其适用于深层地表的坡洪积混合土的渗透系数的确定。
下面针对上述确定方法的确定过程进行进一步详细的说明。
测试土样包括深层地表土样和/或浅层地表土样。本公开实施例中,浅层地表土样可以采用挖探方法坑取大块土获得,而深层地表土样可以采用钻探方法钻取深层土壤获得。本公开实施例的确定方法适用于不同地表深度的坡洪积混合土。
图2是根据一示例性示出的关于上述的区分粒径的确定方法的流程图。如图2所示,即对应于上述的步骤S12,区分粒径的确定方法包括步骤S121-S122:
步骤S121,获取测试土样的颗粒分布,以得到测试土样的第一颗粒分析曲线。
在对测试土样渗透系数的确定过程中,首先需要获取测试土样的颗粒分布,即测试土样中不同粒径的颗粒占总测试土样总质量的百分比,也称颗粒级配。从而可以基于测试土样的颗粒分布,得到测试土样的第一颗粒分析曲线,通过第一颗粒分析曲线可以直观了解测试土样的粗细程度等。
步骤S122,基于第一颗粒分析曲线,确定用于区分测试土样中的粗颗粒和细颗粒的区分粒径。
当得到测试土样的第一颗粒分析曲线,对第一颗粒分析曲线的曲线特性进行分析,确定出用于区分测试土样中的粗颗粒和细颗粒的区分粒径。
举例而言,获取测试土样的颗粒分布,通常可以采用筛析法、密度计法等,本公开实施例采用筛析法来确定颗粒分布,利用筛析装置对测试土样的颗粒逐级分离,从而得到测试土样的颗粒分布情况,以得到该测试土样的第一颗粒分析曲线,如图3所示,具体方法包括步骤S1211至S1212。
步骤S1211,基于预设的多个颗分粒径,对测试土样进行颗分试验,得到用于表征测试土样颗粒分布的第一颗粒数据。
筛析装置预先设置多个不同的颗分粒径,例如设置颗分粒径为200mm,60mm,20mm、5mm、0.5mm等。进一步的,将测试土样按照颗分粒径从大到小的顺序,逐级经过具有不同颗分粒径的筛子,可以得到对应每个颗分粒径的筛子上面的颗粒质量。将每个颗分粒径的筛子上面的颗粒质量进行统计,最终得到测试土样的第一颗粒数据。
在步骤S1212中,基于第一颗粒数据,得到测试土样的第一颗粒分析曲线。
基于第一颗粒数据,统计小于每个颗分粒径下的颗粒质量占测试土样总质量的百分比。以小于某粒径的颗粒质量占测试土样总质量的百分比为纵坐标,颗粒粒径为横坐标,在单对数坐标轴上绘制测试土样的第一颗粒分析曲线。因此,得到测试土样的第一颗粒分析曲线。
下面基于上述方法,以某个工程现场的混合土的具体数据为例来进一步说明。现场获取到四种不同地表的混合土测试土样,四种测试土样取样编号分别为TK2、TK4、TK7和TK9。预设的颗分粒径为400mm、200mm、60mm、20mm、5mm、2mm、0.5mm、0.25mm、0.075mm以及0.005mm。基于预设的颗分粒径,对四种测试土样进行颗分试验,分别得到四种测试土样的第一颗粒数据,如表1所示:
表1
在表1中,例如200mm~60mm所对应的数值10.4,表示该测试土样中颗粒粒径大于60mm且小于200mm的颗粒质量占测试土样总质量的百分比为10.4%。通过表1中的数据可以得到四种测试土样的颗粒分布情况。根据表1的第一颗粒数据,分别计算小于400mm、200mm、60mm、20mm、5mm、2mm、0.5mm、0.25mm、0.075mm以及0.005mm时的颗粒质量占测试土样总质量的百分比,得到下表中的统计数据:
表2
以表3中测试土样TK2的2mm粒径对应的数据为例,由上面表1中的数据可以看出,小于2mm的颗粒质量占测试土样总质量的百分比为表1中2~0.5mm、0.5~0.25mm、0.25~0.075mm、0.075~0.005mm所对应的颗粒数据之和,即11.2+5.6+5.6+1.7=24.1。根据该计算方法,分别计算小于其他颗分粒径所对应的百分比数据,从而得到表3中的统计数据。
进一步的,根据表3中的统计数据,以各个颗分粒径为横坐标,每个颗分粒径所对应的百分比数据为纵坐标,在单对数坐标轴上,绘制得到四种土样的第一颗粒分析曲线。如图4所示,图4中分别包含四种测试土样TK2、TK4、TK7、TK9的第一颗粒分析曲线。
通过图4中的第一颗粒分析曲线,不仅可以直观得到测试土样的颗粒分布情况,还可以获取到用于表征测试土样属性的相关参数,比如测试土样的不均匀系数、曲率系数、测试土样的土壤类别等,如表3中所示:
表3
如表2所示,可以得到四种测试土样的土壤类别包括砾石混土、卵石混土和碎石混土。而TK7和TK9虽然土壤类别同为卵石混混土,但是两种土样的不均匀系数、曲率系数以及不同界线粒径值都不相同,通过对第一颗粒曲线进行分析,可以有效的得到土壤的属性相关参数。
进一步的,在得到测试土样的第一颗粒分析曲线后,分析该曲线的曲线特征,确定测试土样的区分粒径。具体方法如图5所示,包括步骤S1221至S1222。
步骤S1221中,分析第一颗粒分析曲线的平缓曲线段,以得到平缓曲线段的前端点和后端点,其中,前端点所对应的第一粒径大于后端点所对应的第二粒径。
由于坡洪积混合土是由细粒土和粗粒土组成而缺乏中间粒径的土,故而,在对测试土样分析中,小于某些中间粒径的颗粒质量占测试土样的百分比数据之间变化很小。对应到第一颗粒分析曲线上,会相应出现一条平缓曲线段,平缓曲线上的各个点曲率相差很小,说明该测试土样缺乏平缓曲线上各个点所对应的颗粒粒径。进一步的,基于第一颗粒分析曲线的平缓曲线段,得到平缓曲线段前端点和后端点。如图4所示,图4中,以TK2所对应的颗粒分析曲线为例,可以看出平缓曲线段S1,相对其他曲线段相对平缓,那么可以得到曲线段S1的前端点a1和后端点a2。可以看出a1和a2之间的点曲率相差很小,从而形成一段平缓的曲线段S1。另外,从图4中可以看出,前端点a1所对应的第一粒径大于后端点a2所对应的第二粒径。
步骤S1222中,基于前端点所对应的第一粒径和后端点所对应的第二粒径,确定区分粒径,区分粒径大于等于第二粒径且小于等于第一粒径。
如图4所示,基于前端点a1和后端点a2画一条垂直直线,所对应横坐标中的点r1和r2即为第一粒径和第二粒径。测试土样的区分粒径介于第一粒径和第二粒径之间。通过大量实验得出,通常将平缓曲线段中后端点所对应的第二粒径作为区分粒径。
进一步的,当得到测试土样的区分粒径,基于该区分粒径对测试土样中小于该区分粒径的颗粒进行分析,得到有效粒径,如图6所示,具体方法包括S131至S132。
步骤S131,基于区分粒径,分析测试土样中小于或者等于区分粒径下的颗粒分布,以得到测试土样的第二颗粒分析曲线。
本公开实施例中,渗透系数的确定主要考虑到细颗粒含量的影响,区分粒径用于区分测试土样的粗颗粒和细颗粒,从而根据得到的区分粒径,分析小于或者等于区分粒径下的颗粒分布,可以理解为,将测试土样的粗颗粒去除,保留小于或者等于区分粒径的颗粒,进一步具体分析测试土样中的细颗粒的分布情况,以相应得到测试土样的第二颗粒分析曲线。
在步骤S132中,基于第二颗粒分析曲线,确定测试土样的有效粒径。
由于第二颗粒分析曲线考虑到细颗粒对渗透系数的影响,因此,第二颗粒分析曲线是基于得到的区分粒径,对前述的第一颗粒分析曲线的修正曲线。可以看出,基于第二颗粒分析曲线,确定出的测试土样的有效粒径较之第一颗粒分析曲线确定出的有效粒径更为精准。
具体的,如图7所示,第二颗粒分析曲线得到过程包括步骤S1311-S1312:
步骤S1311,基于区分粒径,将测试土样中小于或者等于区分粒径的颗粒保留,得到测试土样的第二颗粒数据。
步骤S1312,基于第二颗粒数据,得到测试土样的第二颗粒分析曲线。
具体的,还以上述实验数据为例,以第二临界点对应的第二粒径作为区分粒径,图4中,TK2对应的区分粒径为r2,r2介于2mm和3mm之间,因此将表1中小于r2的颗粒数据保留,重新计算小于或者等于该区分粒径下的颗粒质量占保留后额颗粒总质量的百分比。从而,TK4、TK7以及TK9三种土样基于与TK2相同的方法,得到表4中的第二颗粒数据:
表4
进一步,根据表4中的数据,分别统计计算小于2mm、0.5mm、0.25mm、0.075mm以及0.005mm时的颗粒质量占测试土样总质量的百分比,得到表5中的统计数据。
表5
根据表5中的统计数据,依据第一颗粒分析曲线同样的绘制方法,相应得到测试土样的第二颗粒分析曲线,如图8所示。
本公开实施例中,以第二颗粒分析曲线上小于某粒径的颗粒质量占测试土样总质量的10%的粒径作为测试土样的有效粒径。由图8中的第二颗粒分析曲线可以看出,以TK2的曲线为例,TK2的曲线上点c1对应纵轴的数值为10,而c1所对应的到横轴的数值即为测试土样的有效粒径,即图中的0.1mm。因此,从图8中的第二颗粒分析曲线,可以分别得到四种测试土样TK2、TK4、TK7以及TK9的有效粒径。
进一步的,基于上述得到的有效粒径,采用相关计算公式,得到测试土样的渗透系数,采用的计算公式如下:
k=Cd10(a+bt)
在上述公式中,k为渗透系数,C、a、b为常数系数,d10为有效粒径,t为水温。常数C通常的取值范围为700~1000,本公开实施例中的计算取值为1000,常数a通常取值为0.7,常数b取值为0.03,水温t取值15℃,但这些计算取值是示例性,本领域技术人员可根据需要进行调整。根据每个公式中每个参数的取值,将四种测试土样的有效粒径分别代入到上述渗透系数计算公式中,分别得到四种测试土样的渗透系数。见下表中所示:
表6
由表6中可以看到四种测试土样的基于第一颗粒曲线得到的有效粒径和基于第二颗粒曲线得到的有效粒径,以及基于修正后混合土有效粒径计算的渗透系数。
根据前文所述,本公开实施例考虑到细颗粒对于混合土渗透系数的影响,因此在对渗透系数的确定过程中,采用基于确定的区分粒径,在对原级配混合土测试土样的第一颗粒分析曲线修正得到的第二颗粒分析曲线上获取有效粒径,为了验证本公开实施例得到的渗透系数的准确性,下表给出了基于第一颗粒级配的有效粒径,第二颗粒级配的有效粒径,以及进行现场单环注水试验得到的三种情况的渗透系数进行对比,见下表7:
表7
由表7中的数据对比可以看出,按照第一颗粒分析曲线得到的有效粒径计算的渗透系数是现场试验获得渗透系数的4.8~209.3倍,而按照第二颗粒分析曲线得到的有效粒径计算的渗透系数是现场试验获得渗透系数的0.8~1.1倍,相比原级配中获得的渗透系数更准确。由此可以得出,采用本公开实施例的方法确定的渗透系数更为准确。
综上所述,本公开实施例的渗透系数确定方法,通过现场获取测试土样,对测试土样进行室内分析实验,确定出区分测试土样粗细粒径的区分粒径,并基于该区分粒径,得到测试土样的有效粒径,从而计算出相应的渗透系数,该方法考虑到坡洪积混合土的颗粒分布特性,得到的渗透系数更为准确,同时对于深层地表的坡洪积混合土相比现场试验方法更为简单易行。
相应的,基于相同的发明构思,图9示出一种坡洪积混合土渗透系数的确定装置100,该确定装置100包括获取单元110,用于获取坡洪积混合土的测试土样。第一确定单元120,用于基于测试土样,确定用于区分测试土样中的粗颗粒和细颗粒的区分粒径。第二确定单元130,用于对测试土样中粒径小于区分粒径的颗粒进行分析,并根据分析结果确定与测试土样的渗透参数相关联的有效粒径;以及第三确定单元140,用于根据有效粒径与渗透参数之间的关联关系,确定测试土样的渗透系数。
本发明实施例的坡洪积混合土渗透系数确定装置与上述坡洪积混合土渗透系数确定方法的实施例的具体实施细节及效果相同,在此则不再赘述。
相应的,本发明实施例提供了一种电子设备,该电子设备包括:至少一个处理器;存储器;其中,所述存储器存储有能被所述至少一个处理器执行的指令,所述至少一个处理器通过执行所述存储器存储的指令实现上述所述的坡洪积混合土渗透系数确定方法。
本发明实施例还提供了一种存储介质,其上存储有程序,该程序被处理器执行时实现所述坡洪积混合土渗透系数确定方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、装置、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、装置、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种坡洪积混合土渗透系数的确定方法,其特征在于,所述确定方法包括:
获取所述坡洪积混合土的测试土样;
基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径;
对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径;以及
根据所述有效粒径与所述渗透参数之间的关联关系,采用如下公式确定所述测试土样的渗透系数:
k=Cd10(a+bt)
其中,k为所述渗透系数,C、a、b为常数系数,d10为所述有效粒径,t为水温。
2.根据权利要求1所述的确定方法,其特征在于,所述测试土样包括深层地表土样和/或浅层地表土样。
3.根据权利要求1所述的确定方法,其特征在于,所述基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径,包括:
获取所述测试土样的颗粒分布,以得到所述测试土样的第一颗粒分析曲线;以及
基于所述第一颗粒分析曲线,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径。
4.根据权利要求3所述的确定方法,其特征在于,所述获取所述测试土样的颗粒分布,以得到所述测试土样的第一颗粒分析曲线,包括:
基于预设的多个颗分粒径,对所述测试土样进行颗分试验,得到用于表征所述测试土样颗粒分布的第一颗粒数据;
基于所述第一颗粒数据,得到所述测试土样的第一颗粒分析曲线。
5.根据权利要求3所述的确定方法,其特征在于,所述基于所述第一颗粒分析曲线,确定所述测试土样的区分粒径,包括:
分析所述第一颗粒分析曲线的平缓曲线段,以得到所述平缓曲线段的前端点和后端点,其中,所述前端点所对应的第一粒径大于所述后端点所对应的第二粒径;以及
基于所述前端点所对应的第一粒径和所述后端点所对应的第二粒径,确定所述区分粒径,所述区分粒径大于等于所述第二粒径且小于等于所述第一粒径。
6.根据权利要求1所述的确定方法,其特征在于,所述对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径,包括:
基于所述区分粒径,分析所述测试土样中小于或者等于所述区分粒径下的颗粒分布,以得到所述测试土样的第二颗粒分析曲线;以及
基于所述第二颗粒分析曲线,确定所述测试土样的有效粒径。
7.根据权利要求6所述的确定方法,其特征在于,所述基于所述区分粒径,分析所述测试土样中小于或者等于所述区分粒径下的颗粒分布,以得到所述测试土样的第二颗粒分析曲线,包括:
基于所述区分粒径,将所述测试土样中小于或者等于所述区分粒径的颗粒保留,得到所述测试土样的第二颗粒数据;以及
基于所述第二颗粒数据,得到所述测试土样的第二颗粒分析曲线。
8.一种坡洪积混合土渗透系数的确定装置,其特征在于,所述确定装置包括:
获取单元,用于获取所述坡洪积混合土的测试土样;
第一确定单元,用于基于所述测试土样,确定用于区分所述测试土样中的粗颗粒和细颗粒的区分粒径;
第二确定单元,用于对所述测试土样中粒径小于所述区分粒径的颗粒进行分析,并根据分析结果确定与所述测试土样的渗透参数相关联的有效粒径;以及
第三确定单元,用于根据所述有效粒径与所述渗透参数之间的关联关系,采用如下公式确定所述测试土样的渗透系数:
k=Cd10(a+bt)
其中,k为所述渗透系数,C、a、b为常数系数,d10为所述有效粒径,t为水温。
9.一种电子设备,其特征在于,该电子设备包括:至少一个处理器;存储器;其中,所述存储器存储有能被所述至少一个处理器执行的指令,所述至少一个处理器通过执行所述存储器存储的指令实现权利要求1至7中任意一项权利要求所述的确定方法。
10.一种机器可读存储介质,该机器可读存储介质上存储有指令,该指令用于使得机器执行权利要求1至7中任意一项所述的确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011625598.5A CN112816387B (zh) | 2020-12-31 | 2020-12-31 | 一种渗透系数的确定方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011625598.5A CN112816387B (zh) | 2020-12-31 | 2020-12-31 | 一种渗透系数的确定方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112816387A CN112816387A (zh) | 2021-05-18 |
CN112816387B true CN112816387B (zh) | 2022-07-12 |
Family
ID=75856306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011625598.5A Active CN112816387B (zh) | 2020-12-31 | 2020-12-31 | 一种渗透系数的确定方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112816387B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116930042B (zh) * | 2023-09-19 | 2023-12-01 | 常州建昊建筑鉴定检测有限公司 | 建筑防水材料性能检测设备及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2438134A1 (en) * | 2002-09-30 | 2004-03-30 | William Daniel Gunter | Process for predicting porosity and permeability of a coal bed |
CN2630839Y (zh) * | 2003-07-22 | 2004-08-04 | 河海大学岩土工程科学研究所 | 粗颗粒土渗透系数及土体渗透变形测试仪 |
CN107228817A (zh) * | 2017-06-12 | 2017-10-03 | 山东大学 | 地质构造内部充填介质/土体渗透稳定性综合判别方法 |
CN108956415A (zh) * | 2018-05-31 | 2018-12-07 | 中国科学院力学研究所 | 一种针对非常规储层样品的相对渗透率曲线计算方法 |
CN112067514A (zh) * | 2020-08-12 | 2020-12-11 | 中铁十二局集团有限公司 | 基于土工筛分试验的土壤粒径检测方法、系统及介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0983501A1 (en) * | 1996-11-08 | 2000-03-08 | Purdue Research Foundation | Particle analysis system and method |
JP3723779B2 (ja) * | 2002-02-15 | 2005-12-07 | 日本国土開発株式会社 | 多層覆土の遮水性能評価装置 |
US20100089124A1 (en) * | 2008-09-26 | 2010-04-15 | North Dakota State University | Integrated porous rigid wall and flexible wall permeability test device for soils |
-
2020
- 2020-12-31 CN CN202011625598.5A patent/CN112816387B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2438134A1 (en) * | 2002-09-30 | 2004-03-30 | William Daniel Gunter | Process for predicting porosity and permeability of a coal bed |
CN2630839Y (zh) * | 2003-07-22 | 2004-08-04 | 河海大学岩土工程科学研究所 | 粗颗粒土渗透系数及土体渗透变形测试仪 |
CN107228817A (zh) * | 2017-06-12 | 2017-10-03 | 山东大学 | 地质构造内部充填介质/土体渗透稳定性综合判别方法 |
CN108956415A (zh) * | 2018-05-31 | 2018-12-07 | 中国科学院力学研究所 | 一种针对非常规储层样品的相对渗透率曲线计算方法 |
CN112067514A (zh) * | 2020-08-12 | 2020-12-11 | 中铁十二局集团有限公司 | 基于土工筛分试验的土壤粒径检测方法、系统及介质 |
Non-Patent Citations (2)
Title |
---|
大孤山水电站坝基覆盖层渗透特性及渗流控制参数的确定;张建红;《甘肃水利水电技术》;20100430;第46卷(第4期);第2-3节 * |
级配不连续碎砾石粗细颗粒区分粒径比较分析;齐俊修等;《水利与建筑工程学报》;20160229;第14卷(第1期);第6节 * |
Also Published As
Publication number | Publication date |
---|---|
CN112816387A (zh) | 2021-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10641088B2 (en) | Method and device for determining karst development degree of reservoir, computer readable storage medium and device | |
Muszynski et al. | Particle shape estimates of uniform sands: visual and automated methods comparison | |
CN111323844A (zh) | 一种基于曲线重构的复杂砂砾岩体的岩性识别方法及系统 | |
CN110411916B (zh) | 一种巨粒土的颗粒级配测试方法 | |
CN112816387B (zh) | 一种渗透系数的确定方法、装置及存储介质 | |
CN107941670B (zh) | 一种岩屑孔隙度测定方法 | |
Avşar | Contribution of fractal dimension theory into the uniaxial compressive strength prediction of a volcanic welded bimrock | |
CN109389154B (zh) | 洪积扇砂砾岩储层岩相识别方法及装置 | |
CN105089662B (zh) | 碳酸盐岩地层录井综合柱状图的校正方法及系统 | |
CN109298464B (zh) | 致密砂岩储层成岩相测井识别方法及装置 | |
Vakili et al. | Filtration of broadly graded cohesive dispersive base soils | |
CN103149130B (zh) | 砾岩岩心颗粒结构中粒度的分析方法 | |
CN108153947B (zh) | 一种基于软岩现场状态的崩解特征的描述方法 | |
CN115390155A (zh) | 针对非均质储层的测井解释方法、装置、电子设备及介质 | |
Sahu et al. | Empirical methods and estimation of hydraulic conductivity of fluvial aquifers | |
Yu et al. | Model identification and control of development of deeply buried paleokarst reservoir in the central Tarim Basin, northwest China | |
CN106094053A (zh) | 一种烃源岩生排烃强度计算方法及装置 | |
CN110399649A (zh) | 一种基于成岩指示元素的变质岩量化识别方法 | |
CN117388945A (zh) | 一种断缝体油藏储层类型的识别方法 | |
CN112989257B (zh) | 海相页岩油气储层的产气量测定方法 | |
CN110489769B (zh) | 一种储油层泥质含量计算方法及电子设备 | |
CN112987126A (zh) | 一种岩石物理相的地球物理预测方法 | |
Melzer et al. | Retention of high permeability during shallow burial (300 to 500 m) of carbonate grainstones | |
CN112348330A (zh) | 确定蚀变火山碎屑砂砾岩储层定量分类的方法 | |
Lipiński et al. | A study of granulation distribution in a beach of tailings dam |
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 |