CN108801872A - 一种岩土材料渗流系数相关偏度确定方法 - Google Patents
一种岩土材料渗流系数相关偏度确定方法 Download PDFInfo
- Publication number
- CN108801872A CN108801872A CN201810350260.XA CN201810350260A CN108801872A CN 108801872 A CN108801872 A CN 108801872A CN 201810350260 A CN201810350260 A CN 201810350260A CN 108801872 A CN108801872 A CN 108801872A
- Authority
- CN
- China
- Prior art keywords
- rock
- seepage coefficient
- bias
- soil
- correlation degree
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 50
- 239000000463 material Substances 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000000528 statistical test Methods 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims description 16
- 238000002474 experimental method Methods 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 3
- 230000035515 penetration Effects 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 5
- 238000001914 filtration Methods 0.000 abstract 1
- 239000011435 rock Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/082—Investigating permeability by forcing a fluid through a sample
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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种岩土材料渗流系数相关偏度确定方法,包括步骤:1)岩土层原状土现场取样;2)获取不同层位岩土材料渗流系数;3)统计测试渗流系数数字特征;4)分析局部平均渗流系数数字特征;5)计算局部平均条件下的方差折减系数Г2(L);6)绘制(1+L2)/Г2(L)‑L3关系曲线图;7)确定岩土材料渗流系数相关偏度。通过获得渗流系数相关偏度的定量值,实现不同复杂地质条件下形成的岩土层渗流系数随机性评价,为岩土结构随机渗流场的计算提供有效依据。该方法实施简单,操作方便,计算准确,实用性强,可广泛应用于不同岩土地层渗流系数随机性评价分析中。
Description
技术领域
本发明涉及一种渗流系数相关偏度确定方法,尤其是一种适用于岩土材料垂直方向渗流系数相关偏度确定方法。
背景技术
岩土体各点参数是确定的(不考虑时间效应),然而人们无法确切知道参数在每点的数值,只能通过有限样本点去估计未知点。同时,岩土体往往是成层沉积而成,成层状分布。由于在沉积过程中扁平颗粒的取向关系,使得土体在水平面方向和竖直方向的性质存在较大的差异。如何依据已有的岩土资料及工程情况对岩土工程可靠性做出合理的评价是当前工程设计及可靠性评价中遇到的难题,解决该问题的关键在于岩土参数空间变异性的合理评价。Vanmarcke提出的土性参数随机场模型,其实质是用齐次正态随机场模拟土性剖面参数,用自相关函数刻画岩土材料的自相关性,确立了由试验数据求得的点特性过渡到空间平均特性的方差计算方法,能够有效反应空间变异性及相关性特征,建立在该模型上的随机有限元法正成为岩土工程可靠性分析强有力的工具。
然而,在渗流条件下,要想准确描述岩土材料渗流系数的空间变异性,其问题之一是如何确定岩土材料垂直方向渗流系数相关偏度。目前,可准确确定岩土材料渗流系数相关偏度的方法尚不多见,特别是可供一般工程人员确定岩土材料垂直方向渗流系数相关偏度的简便实用方法鲜见报道,因此,亟需一种可以有效确定岩土材料渗流系数相关偏度的分析方法,进而为岩土工程随机渗流场的计算提供有效依据。
发明内容
技术问题:本发明的目的是针对现有技术中存在的问题,提供一种计算可靠、精度较高、有利于工程设计人员使用的岩土材料渗流系数相关偏度确定方法,以解决现有岩土材料渗流系数相关偏度无法有效确定的问题。
技术方案:为实现上述目的,本发明的岩土材料渗流系数相关偏度确定方法,包括以下步骤:
步骤1:岩土层原状土现场取样
采用钻孔法在岩土层垂直深度范围内进行间隔距离Δr的多个原状岩土现场取样;
步骤2:获取不同层位岩土材料渗流系数
将不同层位获得的多组原状岩土分别进行变水头渗透试验,根据公式(1)确定试验原状岩土的渗流系数:
式中:a为玻璃管断面积,l为试样长度,A为试样断面积,Δh1为起始水头差,Δh2为终了水头差,t为渗透时间。
步骤3:统计测试渗流系数数字特征
依据数理统计方法,对试验获取的渗流系数测试值进行样本点统计分析,获取点均值μ及点方差σ2;
步骤4:分析局部平均渗流系数数字特征
沿取样深度范围内选定多个取值间距L,在每个间距L内取样点测试值均值构成一组数据,依据数理统计方法,求出该组数据的局部平均均值μ(L)及局部平均方差D2(L);
步骤5:计算局部平均条件下的方差折减系数Г2(L)
根据公式(2)确定与不同取值间距L相对应的渗流系数方差折减系数Г2(L);
式中:L为局部平均范围,D2(L)为局部平均方差,σ2为测试值点方差。
步骤6:绘制(1+L2)/Г2(L)-L3关系曲线图;
以(1+L2)/Г2(L)为横坐标,L3为纵坐标,将渗流系数离散点绘制在该坐标系中;
步骤7:确定岩土材料渗流系数相关偏度;
以一条直线拟合(1+L2)/Г2(L)-L3坐标系中的离散点,直线的斜率即为该场地竖直方向渗流系数相关偏度。
所述的步骤1中,岩土层垂直深度范围为5~15m,间隔距离Δr为0.8~1.2m。
所述的步骤4中,选定取值间距L为实际取样间距的整数倍,即L=nΔr。
所述的步骤4中,每一组数据由相同取值间距L对应的均值组成。
有益效果:由于采用上述技术方案,本发明计算可靠、精度较高、有利于工程设计人员使用的岩土材料渗流系数相关偏度确定方法,解决了现有岩土材料渗流系数相关偏度无法有效确定的问题。本发明与现有技术相比具有如下优点:
1.本发明通过简单的现场取样、渗透试验、统计分析及直线拟合即可有效确定岩土材料渗流系数相关偏度,整个步骤清晰,物理意义明确,操作方便,实用性强,得到的渗流系数相关偏度可直接应用于岩土结构随机渗流场的计算;
2.本发明的渗流系数相关偏度确定方法不仅仅局限于渗流系数,对于具有空间变异性及相关性特征的其他水力学参数(如扩散系数、比水容)随机场分析同样适用,通用性强。
附图说明
图1为本发明的渗流系数相关偏度确定方法的流程图。
图2为本发明具体实施例的直线拟合图。
具体实施方式
下面结合附图对本发明的实施例作进一步的描述:
本发明的岩土材料渗流系数相关偏度确定方法包括以下步骤:
步骤1:岩土层原状土现场取样
采用钻孔法在5~15m岩土层垂直深度范围内进行间隔Δr=0.8~1.2m的多个原状岩土现场取样;
步骤2:获取不同层位岩土材料渗流系数
将不同层位获得的多个原状岩土分别进行变水头渗透试验,根据公式(1)确定试验原状岩土各自的渗流系数:
式中:a为玻璃管断面积,l为试样长度,A为试样断面积,Δh1为起始水头差,Δh2为终了水头差,t为渗透时间;
步骤3:统计测试渗流系数数字特征
依据数理统计方法,对试验获取的原状岩土渗流系数测试值进行样本点统计分析,获渗流系数的点均值μ及点方差σ2;
步骤4:分析局部平均渗流系数数字特征
沿取样深度范围内选定多个取值间距L,在每个间距L内取样点测试值均值构成一组数据,依据数理统计方法,求出该组数据的局部平均均值μ(L)及局部平均方差D2(L);选定取值间距L为实际取样间距的整数倍,即L=nΔr;每一组数据由相同取值间距L对应的均值组成;
步骤5:计算局部平均条件下的方差折减系数Г2(L)
根据公式(2)确定与不同取值间距L相对应的渗流系数方差折减系数Г2(L);
式中:L为局部平均范围,D2(L)为局部平均方差,σ2为测试值点方差;
步骤6:绘制(1+L2)/Г2(L)-L3关系曲线图;
以(1+L2)/Г2(L)为横坐标,L3为纵坐标,将渗流系数离散点绘制在该坐标系中;
步骤7:确定岩土材料渗流系数相关偏度
以一条直线拟合(1+L2)/Г2(L)-L3坐标系中的离散点,直线的斜率即为该场地竖直方向渗流系数相关偏度。
具体的:
某岩土层原状土的10组渗流系数离散点绘制在(1+L2)/Г2(L)-L3坐标系中,结果如图2所示,渗流系数的拟合直线斜率为1.46,由此可确定该岩土材料竖直方向渗流系数相关偏度为1.46m。
Claims (4)
1.一种岩土材料渗流系数相关偏度确定方法,其特征是:包括以下步骤:
步骤1:岩土层原状土现场取样
采用钻孔法在岩土层垂直深度范围内进行间隔距离Δr的多个原状岩土现场取样;
步骤2:获取不同层位岩土材料渗流系数
将不同层位获得的多组原状岩土分别进行变水头渗透试验,根据公式(1)确定试验原状岩土的渗流系数:
式中:a为玻璃管断面积,l为试样长度,A为试样断面积,Δh1为起始水头差,Δh2为终了水头差,t为渗透时间;
步骤3:统计测试渗流系数数字特征
依据数理统计方法,对试验获取的渗流系数测试值进行样本点统计分析,获取点均值μ及点方差σ2;
步骤4:分析局部平均渗流系数数字特征
沿取样深度范围内选定多个取值间距L,在每个间距L内取样点测试值均值构成一组数据,依据数理统计方法,求出该组数据的局部平均均值μ(L)及局部平均方差D2(L);
步骤5:计算局部平均条件下的方差折减系数Г2(L)
根据公式(2)确定与不同取值间距L相对应的渗流系数方差折减系数Г2(L);
式中:L为局部平均范围,D2(L)为局部平均方差,σ2为测试值点方差;。
步骤6:绘制(1+L2)/Г2(L)-L3关系曲线图
以(1+L2)/Г2(L)为横坐标,L3为纵坐标,将渗流系数离散点绘制在该坐标系中;
步骤7:确定岩土材料渗流系数相关偏度
以一条直线拟合(1+L2)/Г2(L)-L3坐标系中的离散点,直线的斜率即为该场地岩土竖直方向渗流系数相关偏度。
2.根据权利要求1所述的一种岩土材料渗流系数相关偏度确定方法,其特征在于:所述的步骤1中,岩土层垂直深度范围为5~15m,间隔距离Δr为0.8~1.2m。
3.根据权利要求1所述的一种岩土材料渗流系数相关偏度确定方法,其特征在于:所述的步骤4中,选定取值间距L为实际取样间距的整数倍,即L=nΔr。
4.根据权利要求1所述的一种岩土材料渗流系数相关偏度确定方法,其特征在于:所述的步骤4中,每一组数据由相同取值间距L对应的均值组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810350260.XA CN108801872B (zh) | 2018-04-18 | 2018-04-18 | 一种岩土材料渗流系数相关偏度确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810350260.XA CN108801872B (zh) | 2018-04-18 | 2018-04-18 | 一种岩土材料渗流系数相关偏度确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108801872A true CN108801872A (zh) | 2018-11-13 |
CN108801872B CN108801872B (zh) | 2021-02-12 |
Family
ID=64094392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810350260.XA Active CN108801872B (zh) | 2018-04-18 | 2018-04-18 | 一种岩土材料渗流系数相关偏度确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108801872B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110044746A (zh) * | 2019-03-22 | 2019-07-23 | 驿涛项目管理有限公司 | 一种岩土工程勘察用测定装置及测定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011089367A1 (fr) * | 2010-01-22 | 2011-07-28 | Total Sa | Mesure de parametres lies a l'ecoulement de fluides dans un materiau poreux |
CN105628586A (zh) * | 2016-02-19 | 2016-06-01 | 北京理工大学 | 一种多孔介质渗透率的确定方法 |
CN105701274A (zh) * | 2016-01-04 | 2016-06-22 | 中国矿业大学 | 一种岩土参数三维局部平均随机场样本的生成方法 |
CN107063972A (zh) * | 2017-05-23 | 2017-08-18 | 山东大学 | 一种常‑变水头一体化渗透系数测定装置及其试验方法 |
CN107832513A (zh) * | 2017-11-01 | 2018-03-23 | 中科吉奥(北京)地质工程咨询有限责任公司 | 基于统计岩体力学的岩体工程参数计算方法及系统 |
-
2018
- 2018-04-18 CN CN201810350260.XA patent/CN108801872B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011089367A1 (fr) * | 2010-01-22 | 2011-07-28 | Total Sa | Mesure de parametres lies a l'ecoulement de fluides dans un materiau poreux |
CN105701274A (zh) * | 2016-01-04 | 2016-06-22 | 中国矿业大学 | 一种岩土参数三维局部平均随机场样本的生成方法 |
CN105628586A (zh) * | 2016-02-19 | 2016-06-01 | 北京理工大学 | 一种多孔介质渗透率的确定方法 |
CN107063972A (zh) * | 2017-05-23 | 2017-08-18 | 山东大学 | 一种常‑变水头一体化渗透系数测定装置及其试验方法 |
CN107832513A (zh) * | 2017-11-01 | 2018-03-23 | 中科吉奥(北京)地质工程咨询有限责任公司 | 基于统计岩体力学的岩体工程参数计算方法及系统 |
Non-Patent Citations (4)
Title |
---|
TAO WANG 等: "Assessment for the spatial variation characteristics of uncertain thermal parameters for warm frozen soil", 《APPLIED THERMAL ENGINEERING》 * |
刘玲玲 等: "不同地质统计方法在确定渗透系数场中的对比研究", 《水文地质工程地》 * |
张小筱: "泾阳南塬黄土边坡饱和渗透系数特性研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
程演: "基于随机场的非饱和土渗流变形耦合分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110044746A (zh) * | 2019-03-22 | 2019-07-23 | 驿涛项目管理有限公司 | 一种岩土工程勘察用测定装置及测定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108801872B (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
CN106875471B (zh) | 煤系含或隔水层三维可视化建模方法 | |
Painter et al. | Fractional Lévy motion as a model for spatial variability in sedimentary rock | |
CN103628871A (zh) | 一种基于阿尔奇公式的电阻率侵入校正的新方法 | |
CN107678064B (zh) | 一种声波时差实时提取方法 | |
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
CN109187744B (zh) | 一种基于静力触探试验的剪切波速估算方法 | |
Zheng et al. | Estimation of the hydraulic parameters of leaky aquifers based on pumping tests and coupled simulation/optimization: Verification using a layered aquifer in Tianjin, China | |
CN107609759A (zh) | 一种海洋天然气水合物开采的海底工程地质影响评价方法 | |
CN103061320B (zh) | 基于孔压静力触探确定土体渗透系数的方法 | |
CN108801872A (zh) | 一种岩土材料渗流系数相关偏度确定方法 | |
CN107169224B (zh) | 基于cptu测试进行大直径超长管桩可打入性分析方法 | |
CN113409463A (zh) | 一种包括尖灭处理的三维地质模型构建方法及装置 | |
CN103074880A (zh) | 一种有效识别土层界面的微型孔压静力触探探头 | |
CN117131971A (zh) | 基于Xgboost算法的储层优势渗流通道预测方法 | |
CN102134994A (zh) | 一种基于基岩油藏油水层电阻率的地层数据处理方法 | |
CN113534263B (zh) | 一种不依赖测井资料的含油气饱和度预测方法 | |
CN115577616A (zh) | 基于深度学习的碳酸岩缝洞地震刻画方法及装置 | |
CN108535122A (zh) | 一种粘土介质强度参数波动范围确定方法 | |
CN105373643B (zh) | 一种注入聚合物推进的前缘模拟监测方法 | |
Su et al. | Mechanics of aquitard drainage by aquifer-system compaction and its implications for water-management in the North China Plain | |
CN104963328A (zh) | 确定大埋深下粗粒土内摩擦角的方法 | |
Zou et al. | Geostatistical modeling resistivity of cohesionless soil using RCPTU data | |
Kenmogne et al. | Correlation studies between SPT and Pressuremeter tests | |
CN117724167B (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 |