工程爆破二元回归分析参数求解的方法.pdfVIP

2011年第 9期 西部探矿工程 9 工程爆破二元 回归分析参数求解的方法 占学军 ，任贤斌 (长江水利委员会长江科学院，湖北 武汉 430010) 摘 要：给出了工程爆破一元回归及二元回归的回归参数详细求解过程以及计算实例，供技术人员回 归分析时参考。 关键词 ：二元回归；工程爆破 ；质点振动速度 中图分类号 ：TU751．9 文献标识码 ：A 文章编号 ：1O04—5716(2O11)O9一O0O9—03 1 二元 回归分析提 出 lnv=lnK+aln(Q ／。／R) (4) 工程爆破中的爆破质点振动速度衰减规律，通常是 令 y~lnv，口一lnK，b=a，—In(Q ／R) 采用最小二乘法回归分析实测振速 、单响药量Q、测 则式 (4)化为： 点距爆心的距离R的相互关系。最常用的公式是萨道 y=a+ 6z 夫斯基公式，即： 应用数理统计知识，可以求出： — K(Q ／。／R) (1) 6一 )-jxiyi--nxy ，口一 一 6 式中：K、a——待求的与场地、爆破方式等有关的参数 ， 通过统计学上的一元回归分析可求出这两个值。 其中 一 ，一 ，为样本点个数。 这一公式在工程爆破中得到广泛的应用，但是在大 型水电站高边坡、主厂房等部位，由于高边坡高差的放 求出 b之后，就可以求出爆破一元回归系数K— 大效应，导致沿高程增大的方 向振速衰减很慢，甚至出 ea，a=b。 现高程较大的测点的振速大于高程较小的测点，使得该 相关系数： 公式不能适用。在这种情况下，应采用考虑高差的回归 ∑( 一 )( 一 ) 公式： y一 ~／∑(z--x)∑( -y)。 一 K(Q ／／R)epH (2) F检验值 ： 或者是 ： )，2(一2) 一 一 K(Q 。／R)(Q ／。／H)一 (3) 1一y2 式中：H——高差 ； 当F≥F1一 (1，一2)时 (其 中 为显著性水平 ，一 K、a、 待求的参数，卢表示高程影响因子， 般取 co=O．05或0．O1)则认为回归效果显著；否则认为 0为放大，g0为衰减，p=0时表示