1. 引言
为确保施工安全及施工质量,施工期沉降观变得越来越重要。其一方面可验证设计成果,确保施工沉降在允许范围内,另一方面又可指导施工,为优化设计提供参考。特大桥梁桩基沉降是影响工程安全的关键因素,国内外学者已进行了大量研究。林智勇等(2009)建立了与实际工程中相近的低承台群桩模型,用数值模拟手段揭示了在上部荷载作用下群桩的荷载-沉降特性,并对群桩基础荷载-沉降曲线的影响因素分别进行了分析与探讨。陈茜等(2010)提出了在非线性土体中群桩中桩和桩之间的相互影响因子,及预测群桩整体沉降、桩顶刚度的计算方法。邹力等(2010)分析了注浆与未注浆群桩基础的长期沉降监测数据,探讨了后注浆对群桩沉降及其发展规律的影响机制。占川等(2013)提出有效桩长计算方法并对超长桩的受力机理及群桩效应问题进行了探讨,对当前群桩沉降计算方法进行了评价,为超长桩群桩的沉降计算提供一些参考。张乾青等(2013)采用双折线硬化模型模拟桩端位移与单位端阻间的关系,提出单桩受力性状的迭代算法,结合等代墩法建立群桩沉降的简化算法。茅燕兵等(2013)建立了单桩桩顶位移与等效弹性模量的函数关系,并引入一个修正弹性理论误差的经验公式来计算群桩的沉降。陈洪运等(2014)以某一软土地基高速铁路桥梁墩台为监测对象,采用液位-单点联合自动监测法对天津某特大桥沉降进行了长期观测。秋仁东(2015)对长群桩基础的沉降性状进行了较为系统的研究,结果表明:工作荷载条件下,群桩基础中基桩的桩身压缩比例要小于单桩基础。杨军等(2017)针对群桩解析算法用于大型群桩非线性分析时计算量过大的问题,建立一种分层地基中的群桩沉降简化计算模型,避免对大规模群桩满秩矩阵进行求逆运算,有效提高了算法的计算效率 [4] - [14] 。
本文在前人研究的基础上,以灰岩区特大桥梁桩基沉降为例,分析了由于岩溶发育,极有可能产生不均匀沉降或沉降过大等问题,主要探讨了灰岩地区长短桩群众的沉降计算及基于实测数据的沉降预测方法。
2. GM(1,1)模型
灰色简单模型GM(1,1)表示一阶的,一个变量的微分预测模型,用于时间序列预测的是其离散形式的模型 [15] 。
设原始系列为
其1-AGO系列X(1)为:
对一阶生成数列x(1),建立GM(1,1)模型
(1)
由导数的定义有
,因为一般预测的系列都在时间上是离散的,所以我们以离散的形式表示,则有:
又由于离散的关系,我们取
为
和
的均值。于是我们可以把公式(1)的微分方程表示为如下的离散方程:
(2)
其中:
。称(2)式为GM(1,1)模型,而把公式(1)的方程称为(2)的白化方程。这样可以得到如下的方程组:
(3)
3. 模型参数求解
式(3)中,a、u为待估参数,-a为发展灰数,反映了
及
的发展态势;u为灰色作用量是从背景值挖掘出来的数据,它反映数据变化的关系,其确切内涵是灰的。灰色作用量是内涵外延化的具体体现,它的存在世区别灰色建模与一般输入输出建模的分水岭 [16] 。将两个待估参数表示为向量形式:
。
对于(1)的方程组,用最小二乘法求解,和一元线性回归的参数估计方法相同,可得:
(4)
式中:
,
4. GM(1,1)预测
式(1)即GM(1,1)的白化方程的解为:
(5)
它的离散形式
(6)
即为GM(1,1)方程的时间响应系列。一般有
,则:
(7)
还原得:
(8)
这就是灰色简单模型的预测公式。
5. 残差GM(1,1)模型
GM(1,1)模型只能精心粗略的模拟,精度不高,平均相对误差可达30%,相对精度不足70%可靠度相对较低。当GM(1,1)模型精度不符合要求时,可使用残差序列建立GM(1,1)模型,对原来模型进行修正,以提高精度 [17] 。
按GM(1,1)模型求解得到
,即
的预测值,计算
的残差序列:
判断可建模残差尾段:
若存在
,
,
的符号一致,
,则称:
为可建模残差尾段。
通过计算可建模残差尾段的一次累加序列,按GM(1,1)模型计算可建模残差尾段的时间响应式。计算残差尾段
的模拟序列:
,这里,
为导数还原值
即:
用
修正
(用一次累加序列的残差修正一次累加序列预测值),称修正后的时间响应式:
(9)
其中残差修正值
的符号应与残差尾段
的符号保持一致。
用
修正
(用原始序列的残差修正原始序列预测值),根据由
到的不同还原方式,得到不同的残差修正时间响应式。
若
则相应的残差修正时间响应式为:
(10)
若
,则相应的残差修正时间响应式为:
(11)
对从式(9)到式(11)中的
还原,也可以采用累减还原式,即:
(12)
6. 新陈代谢GM(1,1)模型
在用等时距GM(1,1)模型对桥梁群桩基础沉降情况进行预测时,随着施工的进行,预测时间逐步推移,原有的预测模型会老化,预测的精度会越来越差,这主要是因为GM(1,1)模型是以灰色模块为基础,其灰平面呈喇叭型展开,未来时刻越远,预测的灰区间越大,因此需要在此基础上做出了改进 [18] [19] [20] [21] 。
由原始序列
建立GM(1,1)模型求得预测值,将最新信息
加入序列,并去掉老信息
,用
,
建立的模型称为新陈代谢GM(1,1)模型。新陈代谢模型在不断补充新信息的同时,及时去掉老化信息,建模序列更能反映系统目前的特征,能更好的揭示系统的发展趋势。
由于施工沉降量测是定期进行的,会不断有新的实测资料可以充实时间序列数据,这种补充了新数据后建立的模型称为新陈代谢模型。新陈代谢GM(1,1)模型是根据已有序列建立GM(1,1)模型预测下一个值,然后把这个预测值补充到已知序列之后,同时去掉第一个数据,保持数列等维,接着再建立GM(1,1)模型预测下一个值,如此逐步预测。考虑到现场条件的不断变化,虽然前期量测的数据对研究后面的沉降特性有一定的作用,但与新得到的量测数据相比,显然后者对预测模型更能反映沉降的最新趋势和变化规律。因此,为了不增加数据序列的长度,可以在增加新数据的同时去掉其第一个已知数据,从而保持数据序列等维。通过大量的建模预测证明,等维新陈代谢模型的预测精度比全数列模型的预测精度高。
7. 某斜拉桥计算实例
该斜拉桥位于湖南宜章县赤石乡境内,桥位跨青头江河道。桥位区内基岩为下石炭统灰岩、炭质灰岩,两者呈互层关系接触。区内岩溶发育,工程地质条件极其复杂。主桥为165 + 3 × 380 + 165 m。四塔预应力混凝土斜拉桥,主塔塔高286 m。主塔7#墩为钻孔灌注桩群桩承台基础,每墩30根大直径钻孔灌注嵌岩桩(30根Ф2.8 m)。最长桩长达到96 m,成孔深度达到102 m。桩基相关参数如表1所示,各承台桩基平面布置如图1所示。
相应的测量(图2)及数据见表2。
预测结果见表3和图3,表3中加粗的数据为依据前期的测量数据的预测值。对比结果发现新陈代谢模型误差较小,预测效果较好。
8. 结论
将灰色预测理论及一般GM(1,1)模型、残差模型和改进的新陈代谢模型对沉降进行预测,通过与湖南某斜拉桥7#墩的沉降观测资料的对比分析,分析了各种预测方法的预测精度。
Table 1. Design parameters for pile of cable-stayed bridge
表1. 斜拉桥桩基设计参数表
(a) (b)
Figure 1. Pile schematic of 7# cable support tower; (a) Plan; (b) Section
图1. 7号索塔桩基布置示意;(a) 平面图;(b) 剖面图
Figure 3. Comparation of measured curve and predicted curve o f 7# pier
图3. 7#墩沉降实测曲线及预测曲线比较
Table 3. Comparation of measured settlement and predicted settlement o f 7# pier
表3. 7#墩群桩基础沉降预测与实测值比较表
在采用灰色理论对斜拉桥群桩基础沉降预测过程中,基于前期实测沉降值按一般GM(1,1)模型、残差模型、改进的新陈代谢模型拟合的结果来看,等时距曲线预测误差较大,随着时间的推移偏离加大,在后期已不再适用。残差模型在前期与实测结果偏离较大,但随着时间的推移预测值与实测值吻合较好;新陈代谢预测的沉降量与实测值拟合较理想,通过补充新的观测据构建动态预测模型,预测精度更高,反映了沉降的变化,更接近实际情况。由此可知新陈代谢GM(1,1)模型在群桩沉降预测上具有较好的实用性。新陈代谢GM(1,1)模型能够很好预测群桩的施工沉降,从而验证高速公路设计并为施工提供参考,指导施工以达到缩短工期、降低造价的目的。
基金项目
国家自然科学基金资助(41174061)。