【摘要】回顾了随机变量的几个重要概率分布,探索了岩土参数标准值的本质,提出岩土参数是随机变量又非随机变量,对岩土参数的统计分析和最小样本数量提出见解。
1、前言
对岩土参数进行分析评价,是岩土工程师的一项重要任务。参数评价之前要进行数理统计,数理统计的基本任务是要给出该参数的平均值、变异系数和标准值。之后,由岩土参数经综合考虑升华为设计参数,再进行岩土工程设计。
本文对岩土参数标准值的内涵实质进行分析,进而提出设计取值建议。
2、随机变量的概率分布回顾
(1)正态分布
设随机变量 X 服从正态分布,其概率密度函数为:
,
累积概率分布函数(随机变量X
其中μ、σ² > 0 为常数,μ 表示均值(又称数学期望值或简称期望值),σ² 表示方差,σ 表示标准差(又称均方差)。
随机变量 X 服从参数为μ、σ² 的正态分布,记为:
对于正态分布,期望值 μ(又称均值)决定了其位置,标准差 σ 决定了分布的幅度。随机变量的几个重要取值区间与其概率面积的关系见下图:
当 μ = 0, σ = 1 时的正态分布称为标准正态分布,其概率密度函数为:
不难导出:
总体的方差计算公式:
其中,为总体方差,
为变量,
为总体均值,
为总体例数。
实际工作中,总体参数难以得到时,可用样本统计量代替总体参数,经校正后,样本方差计算公式为:
式中,S²为样本方差(S 表示样本的标准差,又称均方差),X 为变量, 为样本均值,n 为样本例数。
(2)分布
χ²分布简称卡方分布。设X1,X2,......Xn相互独立,均服从标准正态分布N(0,1), 则称随机变量
的分布为自由度为n 的χ²分布。卡方分布的概率密度函数如下:
χ²分布的期望E(χ2)=n,方差D(χ2)=2n。
(3)t分布
假设X服从标准正态分布即X~N(0,1),Y服从自由度n的卡方分布即Y~χ²(n),且X与Y是相互独立的,那么Z=X/sqrt(Y/n)的分布成为自由的为n的t分布,记为Z~t(n)。
t分布常常用于根据小样本来估计呈正态分布且方差值为知的样本的均值。(如果总体的方差已知的话,则应该用正态分布来估计总体的均值。)
t(n)分布的密度函数表达式、图示如下。
其中表示伽玛函数。
对于Z~t(n),其数学期望E(Z) = 0,n>1;方差D(Z)=n/(n-2) , n>2 。
3、岩土参数的估计
岩土参数的获取首先要进行岩土工程勘察,从数学上看,勘察就是抽样。工程上岩土体的参数平均值是未知的,但存在,要靠勘察得到的参数小样本进行计算估计。
根据高大钊《土力学可靠性原理》,
其中tα为失效概率为α时的 t 分布值,表列如下:
注意表格最后一行的值,自由度n为无限大时的t-分布和正态分布等价。
4、岩土参数标准值的本质
岩土参数标准值,是根据抽样子样本,按照一定失效概率进行对总体样本(假定符合正态分布)平均值的估计。从计算公式看,如取α=5%,当n趋于无穷大时,tα=1.645 ,标准值的计算结果就是平均值。
由此看来,岩土参数具有变异性和随机性双重属性,采用t分布理论进行参数置信估计,当勘察样本数量较小时,就应该运用上述公式,计算出相应于一定置信概率的标准值来对总体样本的平均值进行估计。但当勘察得到样本较大时,可直接采用大容量样本的平均值来代表总体的平均值。
另外,岩土参数标准值的本质,是总体参数平均值的估计值。
从上述表列数值中,摘取n=1,2,3,......,同时查取相应的 tα,用n做横坐标, tα做纵坐标,绘图如下(取α=0.05) 。
从图中可以看出,岩土参数标准值计算的最小数据量可取 5 。当前《岩土工程勘察规范》规定最小取值为6 。
5、结束语
本文就岩土参数的估计进行回顾和分析,对勘察抽取岩土参数样本的标准值,进行概念分析,得出结论,样本参数标准值的本质,实际就是基于一定置信概率的总体参数平均值的估计值。当勘察得到样本较大时,可以采用大样本的平均值来直接代表总体的参数平均值。
对于勘察抽取样本的最少数量可采用 5 。
参考文献:
1、GB/T 15932-1995非中心t分布分位数表,国家技术监督局 ,1996-06-01
2、GB 4086.3-1983统计分布数值表t分布,国家标准化管理委员,1984-10-01
3、高大钊. 土力学可靠性原理. 北京:中国建筑工业出版社,1989
延伸阅读:
关于岩土参数抽样统计的代表性
岩土参数的确定是四维空间问题
压缩模量Es并不是土的基本参数
小程序:
不土不木规范链接