基于条件数学期望花岗岩孤石方量估算方法.docVIP

基于条件数学期望花岗岩孤石方量估算方法 　　摘要：在花岗岩地区进行场地平整设计，一般将孤石并入剥离层计算，因此最终石方量与计算值存在一定的偏差。现行规程规范并无估算孤石方量的方法，因此研究孤石方量的估算方法具有现实意义。本文根据工程实例采用概率统计方法中的条件数学期望进行统计分析，探讨孤石方量的估算方法。 　　关键词：花岗岩孤石方量估算概率统计条件数学期望 　　 　　 　　1前言 　　 在花岗岩地区进行场地平整设计，一般将孤石并入剥离层计算，这种简化的做法造成了石方量计算的偏差。本文所用工程案例位于广州市科学城，工程场地位于覆盖层厚度大且花岗岩孤石较发育的剥蚀残丘，虽然场平设计时钻孔间距达到了20m（图1-1），但最终石方占总方量的比例（13.5%）与计算值（6.5%）仍存在一定的偏差。 　　 现行规程规范并无估算孤石方量的方法，研究花岗岩孤石的估算方法具有现实意义。地表孤石的分布、块径、形状等和地下的孤石有区别，本文仅讨论地下孤石方量估算方法。 　　 本文以钻孔揭露的水平方向遇孤石率，垂直方向线孤石率构造二维随机向量，通过条件数学期望获得二者的函数关系，并进行可靠度估值。最后通过与某工程的平均厚度法计算结果和施工记录进行对比，讨论计算方法的合理性。 　　2借鉴储量计算方法的孤石方量估算 　　 场平设计时常将孤石视为夹层，借鉴储量计算常用的平行断面法、三角法、平均厚度法粗略估算孤石方量。 　　 根据工程经验，平行断面法计算的储量一般＜平均厚度法计算结果；三角法属于一种精度更高的平均厚度法，但在钻孔较多时效率太低，一般不采用，以下只讨论平均厚度法。 　　 平均厚度法中，剥离层、夹层、有用层相当于同底（S）不同高的柱体，该法不考虑孤石的具体形状，因此计算结果较为保守。设开挖面以上，全部钻孔的孤石平均厚度为，基岩面以上风化带厚度（含孤石）平均值为，则孤石方量，剥离层总方量。 　　 　　 　　 本工程采用一种简化的计算方法，通过对各孔孤石线性比例进行加权平均来获得v/V，，式中yi为各孔孤石总厚度Hi与基岩面以上风化带厚度Ti（???孤石）的比值，si为各孔基岩面以上风化带厚度，n为钻孔总数。可见平均厚度法和加权平均法是等价的。 　　 由于钻孔不可能把孤石全部揭露，平均厚度法一方面会遗漏未揭露的石方量，另一方面把孤石作为成层土处理又夸大的石方量，因此其计算精度难以评估。为反映孤石分布的不确定性，下文把概率统计理论应用到方量估算中来。 　　3基于概率统计的估算方法 　　3.1 计算参数 　　 由于孤石之间是相互独立的，水平方向的“遇孤石率”，垂直方向的“线孤石率”具有一定的统计意义。 　　 （1）水平方向计算参数――遇孤石率p 　　 单个钻孔是否遇孤石的随机事件可看作贝努利试验，该事件服从概率为p的两点分布，而n个钻孔相当于n重贝努利试验，因此遇孤石的孔数x服从概率为p的二项分布，即，其数学期望E(x)=np。 　　 　　 　　 对遇孤石率的分区统计见表3-1。 　　 （2）垂直方向计算参数――线孤石率y 　　 把线孤石率y定义为孤石总厚度H与基岩面以上风化带厚度T（含孤石）的比值，y可直观地表示孤石在垂向所占比例，根据钻探资料，对场地线孤石率y进行分区统计如表3-1。 　　 　　 　　 若以y的均值和样本标准差为正态分布参数，对y作Kolmogorov-Smirnov双侧检验，其方法是根据样本经验分布函数Fn(y)和正态分布函数N(y)构造检验统计量，在显著水平0.05下，求得临界值KS0.05=0.2873＞KS=0.271，概率p=0.0747，虽然非常接近拒绝域，但仍接受了正态分布的假设。 　　 y近似服从正态分布是下文计算和讨论的重要前提。 　　3.2 基于条件数学期望的计算方法 　　 对x进行标准化变换，则服从正态分布，因y近似服从正态分布，随机向量(、y)服从二元正态分布，设的均值和方差分别为、，y的均值和方差分别为、。当已确定，则可通过条件数学期望和描述(、y)的统计学特征。 　　 对场地各分区的、y及其相关系数r统计成果见表3-3。 　　 　　 方量估算与平均厚度法采用同样的假设，全场地揭露孤石的面积为p*S，孤石所占厚度为y*，为基岩面以上风化土和孤石的总厚度平均值，则孤石体积为v=φpSy，基岩面以上的总方量为V=S。式中φ为形状系数，取值0.526（椭球体）～0.67（圆柱）～1（六面体），因此孤石体积占基岩面以上总方量的比例为v/V=φpy。为方便与平均厚度法比较，取φ=1。 　　 因，、n、x确定时，遇孤石率p也是定值，。另有，式中、、、、、已知，因此E(v/V)也是的线性函数，可以用均值的置信区间描述E(v/V)的置信区间。