关于岩土参数随机场空间最优估计精度的分析.docxVIP

关于岩土参数随机场空间最优估计精度的分析 1 空间最优估计模型对于岩土工 在岩石参数的不确定性和离散性研究中，不可避免地会出现岩石参数空间的最佳估计，其精度直接关系到岩石参数的离散度分析的效果以及岩石工程的勘察和设计方案的选择。空间最优估计方法在岩土参数不确定性与离散性研究中的应用,为各种岩土参数的局部估计提供了衡量估计精度的定量指标,即克立格方差。 从分析克立格方差的特点出发,讨论岩土参数随机场空间最优估计的精度问题,并针对岩土介质发育的不均匀性,探讨岩土参数空间最优估计中特异值的处理方法。对完善岩土参数随机场空间变异性分析方法体系具有特殊意义。 2 本地评价精度的定量指标和特点 2.1 空间最优估计体系 岩土参数随机场空间最优估计的成果主要体现在空间最优估计值与局部估计精度定量指标两个方面,二者不可分割,构成了一个完整的空间最优估计体系。其中局部估计精度定量指标——普通克立格方差表达式为: σ2Κ=n∑i=1λ0iγ(hi0)+μ(1)σ2K=∑i=1nλ0iγ(hi0)+μ(1) 式中μ为拉格朗日乘子;λ0i0i为克立格权系数;γ(hi0)代表岩土参数随机场待估点x0与信息点xi之间的变异函数值,二者均可由岩土参数空间结构分析结果获得。 2.2 拉格朗日乘子以及克立格权系数 由局部估计精度表达式(1)可见,普通克立格方差的大小是由参估信息点与待估点之间的变异函数值、拉格朗日乘子以及克立格权系数所决定的,而与参估信息值的大小无关。因此,只要参估信息的数据构型与估计构型一定,便可以在待估点取得实测值之前根据变异函数值以及普通克立格方程组的求解结果,事先计算出勘察场地内待估点的局部估计精度,从而为确定最佳观测、取样、布孔方案及选择设计方案提供依据。 3 最合适的地质参数压缩精度分析 3.1 克立格分析法 岩土参数随机场空间最优估计方法是将具有不确定性的岩土参数视为区域化变量,应用克立格法对其离散性进行最优估计的方法。笔者曾对勘察现场获得的一批离散的岩石透水性指标进行了不同方法空间估计效果的对比分析,其结果表明,该法在所研究的岩土参数离散性的空间估计中确比传统估计方法有着更强的有效性和无偏性。 3.2 参估误差估计精度的限制 为了结合实际分析岩土参数空间最优估计的局部估计精度,本文从勘察现场选择了部分代表性区域化变量进行变异函数计算和空间最优估计,通过对其检验孔的克立格方差(σ2K)及估计误差的分析对比,获得了如下几点认识: (1)区域化变量的空间变异性直接影响着空间最优估计的精度。勘察现场部分区域化变量理论变异函数参数值见表1。不同岩土参数区域化变量空间变异程度存在着明显的差别,这种差别必然会反映在空间最优估计的精度上,空间变异性愈强,其空间估计的总体精度愈低(表2、表3)。 (2)克立格方差(σ2K)是衡量岩土参数随机场空间最优估计精度的指标,它描述了在一定数据构型与估计构型下,利用各种信息数据进行空间估计的估计方差的取值中心,即“期望值”,因而才表现出与参估信息值大小无关并可事先进行精度预测的特点。一般来讲,当估计构型表现为待估点位于周围信息点之间时,其估计精度高(σ2K小);当估计构型表现为待估点位于信息点的边缘或待估点周围信息点较少时,其估计精度低(σ2K大)。另外,当待估点距信息点较远时,估计精度较低(σ2K较大);待估点距信息点较近时估计精度较高(σ2K较小),见表4。 (3)当岩土参数区域化变量Z(x)具有各向同性结构时,在几何形状、位置上与待估点x0有对称性的信息样品在没有其它效应影响的条件下,具有相同的权系数。但自然界中的岩土介质是各向异性的,在这种客观条件下即使信息点位置与待估点有对称性,也不能得出克立格权系数有对称性的结论。 (4)在野外实际工作中,岩土工程勘察取样、测试点的分布通常是不均匀的。在这种情况下,即使空间变异结构是各向同性的,位于与待估点对称位置上的两个信息样品的克立格权系数,也会因受到分布不均匀的其它信息样品的影响而失去对称性。因此,用克立格法对岩土参数区域化变量进行估值时,不会因某个区域的信息样品丛聚在一起而过分增大这一区域信息样品总的权系数,这就是所谓减弱丛聚效应。此外,克立格权系数还具有屏蔽效应,即当块金常数为零或很小时,距离待估点最近一圈样品的权系数明显地要比待估点周围第二圈样品的权系数为大。但这种屏蔽效应随着块金常数值的增大而逐渐减弱。 总之,影响岩土参数随机场克立格估计精度的基本因素主要包括岩土参数空间结构的连续性及变异性γ(h)、信息样品与待估点间的几何关系n∑i=1λ0iγ(hi0)、数据构型特征n∑i=1n∑j=1λ0iλ0jγ(hij)、估计块段时待估点所位于的待估块段本身的几何特征ˉγ(V,V)等4种因素。为提高估计精度,在可能条件下,可由调整后3项因素获得。这也为岩