灰色预测在大跨径刚构桥施工线形监控中的应用

1982 年，我国华中科技大学教授邓聚龙先生提出了灰色系统理论［1］，其属管理工程学范畴，是一种研究具有信息不完备特征的系统问题的数学方法，可用于数据控制和预测。该理论以灰色命名，意在突出所研究系统问题初始信息的不完备、半透明状态。与其对应的，内部特征完全已知的系统称为“白色系统”，完全未知称为“黑色系统”［2］。

灰色系统理论以“部分信息已知、部分信息未知”的“小样本、贫信息”的不确定性系统为研究对象，利用对部分已知信息的深度挖掘，提取有价值的信息，获取系统的内在规律，并用量化的方法对系统规律进行描述和建模，实现对不确定系统的运行行为、演化规律的正确描述和有效监控，从而解决数据间的未知关系或原始数据的不确定性问题。

与其他2 种常用的研究不确定性系统的方法（模糊数学、概率统计）相比，三者研究对象的“不确定性”有所区别：模糊数学重点研究“认知不确定”问题，具有“内涵明确、外延不明确”的特点；概率统计重点研究“随机不确定”问题，考察多种可能结果的分布情况；灰色系统理论重点研究“小样本、贫信息”、“部分已知”的问题，具有“外延明确、内涵不明确”的特点。

目前，灰色系统理论在工程控制、经济管理、社会系统、生态系统、农业系统等诸多领域都有成功应用案例［3］。在20世纪90年代，灰色系统理论已开始应用于桥梁工程领域。方志等［4］利用灰色系统理论实现了对斜拉桥的线形和索力的有效控制。张劲泉等［5］研究了灰色系统在悬索桥主缆架设阶段的控制方法。张永水等［6］建立的灰色控制程序在重庆黄花园嘉陵江大桥的施工监控中取得了较好的效果。包仪军等［7］在梁拱组合结构体系桥梁的施工控制中应用了改进型的灰色预测模型。杨正虎等［8］利用灰色理论对成渝客运专线赵家院双线特大桥桥墩沉降进行预测，发现灰色预测模型预测精度较高，有很好的实用性。值得一提的是，王明慧［9］利用灰色系统理论中的灰色关联分析法进行了某跨谷桥桥型方案的选择，是灰色理论在桥梁工程领域除施工控制方面的另类应用范畴。

本文研究了灰色系统理论中灰色预测在一大跨径刚构桥施工线形监控中的应用案例，并对应用过程中的预测精度及影响因素进行了分析。

我把泡好的茶端到杨公子的面前，一不留神手上一滑把茶水洒到他白净的衣服上，当然我没有顺势坐到他身上，那样目的太明显，我只是装着毫不在意地拿出手绢替他擦干，他把我的手轻轻地推开，温雅地说“不要紧”。他抬头看我的时候，我承认自己被电到了，他的眼神很像这个时代的湖水一样，干净清澈。我很没出息地低下了头，剧情明明不是这样写的，难道男主角不该轻易爱上穿越过去的女主角吗？

1 大跨径刚构桥施工线形监控

施工监控是桥梁施工成桥状态达到设计目标的必要手段，连续刚构桥、拱桥、斜拉桥、悬索桥的施工一般都要做施工监控。大跨径桥梁的施工是一项极为复杂的系统工程，受材料性能偏差、施工精度偏差、环境条件变化等各种不确定因素的影响，如果不进行施工监控，很难达到设计预期的成桥线形目标。

大跨径刚构桥的施工监控主要包括墩身线形监控和主梁线形监控2部分内容。其中主梁线形是所有大跨径刚构桥都必须严格控制的要素，因为作为多次超静定体系，合理的成桥线形与受力状态不仅取决于设计，还与科学的施工方法直接相关。线形监控的目的就是通过数据处理、预测分析和实时调整，使大桥实际线形尽可能地与设计理论线形吻合。

目前，大跨径桥梁上部结构施工线形监控的方法主要有2 种：①无应力状态控制法，即基于“各构件无应力构形唯一确定”的原理，利用控制施工过程中各阶段的构件无应力长度来实现成桥线形的最终控制；②自适应控制法，即基于现代控制理论中的自适应原理，对比施工过程中监测数据的实测值和预测值，识别桥梁结构施工中的主要参数，找出实测值与预测值之间产生偏差的原因，从而修正参数，实现对施工线形的监控目标。无应力状态法一般应用于斜拉桥、悬索桥、吊杆拱桥等缆索支承类桥梁的施工线形监控，而大跨径刚构桥的施工线形监控更多采用自适应控制法［10］。

在自适应控制法中又可利用多种方法预测主梁各节段的高程，包括卡尔曼滤波法、最小二乘法、人工神经网络法以及灰色系统理论法。本文研究利用灰色系统理论中的灰色预测方法控制大跨径刚构桥悬浇施工线形的技术，达到精准预测、跟踪控制、随时纠偏的目标，进而最终实现设计理论线形。

2 工程概况

大桥全长960 m，为6跨变截面预应力混凝土连续刚构桥，采用（95+170+3×190+125）m 的跨径布置方案。桥面宽23.4 m，双向4 车道，中央设宽4 m 人行道（含护栏及管线区），两侧各设宽2 m检修道（含护栏及管线区）。桥墩为箱形墩，其中1#墩与5#墩为双肢墩，2#，3#，4#墩下半部分为整体墩，上半部分为双肢墩。各墩基础均为扩大基础。

大桥主梁采用单箱双室箱形截面。顶板宽度为22.7 m，厚0.3 m，悬臂长度为3.5 m；斜腹板构造，从跨中到墩顶腹板斜率保持不变，中墩墩顶位置腹板厚0.85 m，跨中位置腹板厚0.45 m；箱梁高度及底板厚度均按1.8次抛物线变化，中墩墩顶箱梁高11.5 m，底板厚度1.2 m，宽9 m，跨中位置底板箱梁高3.8 m，底板厚0.3 m，宽14 m。桥型布置及典型横断面见图1。

该桥上部结构主梁采用挂篮悬臂浇筑法施工，其中长度为94 m 的悬臂共划分为25 个节段，长度为74 m 的悬臂共划分为20 个节段，95 m 边跨设长19 m的支架现浇段，125 m边跨设长29 m的支架现浇段，各跨合龙段长度皆为2 m。

图1 桥型布置及典型横断面（单位：cm）

大桥整体施工流程如下：①施工各墩基础及墩身；②以桥墩为中心，挂篮悬臂浇筑箱梁节段，形成T构；③支架现浇边跨合龙段并合龙边跨；④合龙中跨；⑤合龙次边跨；⑥施工桥面系。

3 有限元模型

利用MIDAS/Civil 建立该桥有限元模型。主梁、桥墩、扩大基础均采用三维梁单元模拟，墩顶单元与主梁0 号段单元之间的连接采用刚性连接模拟，边墩支座采用弹簧单元模拟，扩大基础单元底部按固结约束模拟。挂篮和施工过程中的其他构件均采用等效荷载模拟。全桥共划分为613 个节点，543 个单元，其中1 ～344 为主梁单元，从基础施工至成桥通车共设113个施工阶段。有限元模型见图2。

回家后，一杭不放心，给雪萤打了一个电话，电话里传来的却是范坚强的声音。一杭大吃一惊，说：“这件事情与雪萤无关，你不要为难她，想报仇直接冲我来。”范坚强皮笑肉不笑地说：“也可以，你把记事本的原件拿来，不然——”电话里突然传来雪萤的尖叫，一杭心里一紧，范坚强继续说：“听到了吧？明天，上午9点，

图2 有限元模型

4 灰色预测在大跨径刚构桥施工线形监控中的应用研究

4.1 灰色预测

灰色预测是灰色系统理论的主要内容之一，其基本原理是利用灰生成对原始数据序列进行加工处理，挖掘数据序列关系，建立灰色模型，预测数据变化趋势。

2.1 两组患者疗效比较 25.0%(10/40)的手术组患者术后出现声音嘶哑、肌痛、肠胃不适及高血压等并发症，显著高于研究组的7.5%(3/40)，差异有统计学意义(P<0.05)。研究组的治疗总有效率为97.5%(39/40)，显著高于手术组的90.0%(36/40)，差异有统计学意义(P<0.05)。见表1。

她说自己叫韩贝，捡到了阿巴的默写本，特地来还给他。韩贝看到默写本上那么多红叉叉时，皱起了眉头：“我可不想我的新朋友是个差生。如果你明天默写得满分，我就带你上船、出岛。”

灰色模型即GM（n，h）模型，代表n 阶h 个变量的灰 方程，G 即gray，M 即model。常用的灰色模型有GM（1，1）模型、GM（1，N）模型、GM（0，N）模型、GM（2，1）模型，其中GM（1，1）模型是灰色预测的核心。

1.3 评价指标 ①痛风相关知识了解情况：出院时利用自制痛风相关知识调查问卷，对患者通风知识了解程度进行判定，包括痛风危险因素、病因、并发症、治疗方案及依从性等多方面，每方面分为了解(1分)、部分了解(0.5分)及不了解(0分),满分10分，得分越高代表患者对相关知识越了解。②疼痛及疼痛缓解情况：利用疼痛视觉模拟评分方法(VAS)，对患者治疗前及出院时疼痛情况进行评估，并记录患者疼痛缓解时间。③生活质量：依据生活质量综合问卷(GQOLI-74)进行评价，包括躯体功能、心理功能、社会功能以及物质生活状态等方面，共20各问题74个项目，每个项目1～5分，分数增高表示患者生活质量提高。

4.2 GM（1，1）模型

建立GM（1，1）模型的基本步骤如下。

1）设有原始数据序列X（0）

2）X（1）的1-AGO序列，即一次累加数列为

3）X（1）的紧邻均值生成序列为

4）X（0），Z（1）中各时刻的数据满足关系：

式中：a为发展系数，反映的发展态势；b为灰色作用量，反映数据变化的关系。

式（3）为GM（1，1）模型的方程表达式。

4.2.1 预测模型

式（3）的矩阵形式表达为

由表3可以看出，不同退耕还林地间作黑麦草后均提高了林木青海云杉的生长量，不同退耕还林地林木青海云杉的树高、冠幅和新梢长度平均值分别为1.73m、1.64m和16.44cm，较对照(CK)平均值分别提高了14.57%、22.39%和19.22%。

式中：

墙夼水库位于潍河上游、山东省诸城市枳沟镇，1960年8月建成蓄水，是一座以灌溉、防洪为主，兼有发电、养殖、供水等功能的大型水库。水库由既独立又相互连通的东西两座水库组成，控制流域面积656km2。其中东库拦截潍河上游两条主要支流高泽河和仲崮河，控制流域面积386 km2；西库拦截潍河主河道，控制流域面积270 km2。东库流域长度29.6 km,干流平均坡度0.002 62；西库流域长度30.1 km，干流平均坡度0.002 2。流域形状呈双叶形。地势东南高，西北低。

灰微分方程的最小二乘参数列满足：

式（6）为式（3）的白化方程，又称影子方程。求解式（6），得到时间相应函数为

则GM（1，1）模型的时间响应序列为

因故有

经累减还原后得到原始数据的响应为

当GM（1，1）模型的精度不符合要求时，可用残差序列建立GM（1，1）模型对原来的模型进行修正，以提高精度。

1) 舒适度。中国船级社(China Classification Society, CCS)在国际海事组织(International Maritime Organization, IMO)SOLAS修正案决议的基础上制定了“舒适度附加标志”的细则，包括COMF(NOISE N)噪声舒适性附加标志和COMF(VIB N)振动舒适性附加标志,N为舒适度等级1，2，3，其中：1表示舒适度最高等级；3表示可接受舒适度等级。由此可看出，船级社在随着航运发展的变革,逐步将客滚船舒适度要求标准化、具体化、规范化。

4.2.2 残差修正

令

◎一旦发生高热惊厥，要冰水敷额、使用退热药物、侧卧位防止呕吐堵塞呼吸道、保持呼吸道畅通，同时立即就医。

目前rGERD的治疗尚无肯定的长期疗效，而导致rGERD治疗失败的原因是多因素的，如患者依从性差、PPI快代谢、抑制胃酸不足、不正确的用药时间、解剖异常、精神心理异常[10，11]，为此各国研究者努力探索安全有效、不良反应少的治疗方法，在药物、内镜和外科手术等方面制定了多种治疗策略，目前各种疗法尚缺乏长期有效有证据。

则的残差序 列。若∃k0，满 足：①∀k ≥k0，ε( 0 )(k)的符号一致；，则称为可建模残差尾段，仍记作ε( 0)。该序列的1-AGO序列的GM（1，1）模型的时间相应为残差尾段的模拟序列其中：

旅客快走完了，谁也没停下来。咦？难道没来？掌声、锣鼓声稀拉下来。就在大家快要失望的时候，一个姑娘微笑着款款地向我们走来。我顿觉眼前一亮，她身穿紧身连衣裙，明显勾勒出健美的曲线，乌黑的披肩发裹着一张酷似桃花的脸，毛嘟嘟的睫毛包着一双神采飞扬的大眼睛——真乃顾盼生姿，窈窕清秀。这是一张多么熟悉的脸，好像在哪见过？呵！我终于想起来了。与此同时，巴克夏也惊喜地俯在我的耳边说：“是她？”

用，修正后的时间响应式为

式（14）称为残差修正GM（1，1）模型，简称残差GM（1，1），其中残差修正项的符号与残差尾段的符号一致。若则相应的残差修正时间响应式为

虽然各种发酵剂菌株在发酵香肠中已有广泛的应用，但鲜有将鲁氏酵母加入发酵香肠中，并研究其对发酵香肠中品质影响的相关报道。因此，将鲁氏酵母菌加入发酵香肠中，探究其是否能够提高发酵香肠的风味，是否具有增鲜增香的作用，将对鲁氏酵母的开发提供良好的基础，为生产新型发酵食品提供依据。本研究通过菌种的耐酸性、耐食盐能力、耐亚硝酸盐能力、生长温度和最适生长温度等试验，对其发酵生产性能进行测定，同时，通过吲哚试验、溶血试验及药敏试验，确定菌种的安全性能。判断其是否符合应用于发酵肉制品的要求，以便生产出新型发酵肉制品。

式（15）被称为累减还原式的残差修正模型。

当预测结果与实际数据出现很不相符的现象，可以考虑用建立残差GM（1，1）模型对预测值进行修正，提高预测精度，但需先满足残差模型的假设条件。

4.3 灰色预测在施工线形监控中的应用

大跨径刚构桥采用悬臂浇筑施工时，悬臂段每节段施工过程可分为3道工序：挂篮立模、浇筑混凝土和张拉预应力，其中最关键的工作是挂篮立模时给出准确合理的抛高值。通过有限元模型仿真分析，可以计算出基于设计理论状态的抛高值，但受到多种因素的影响，实际施工中节段完成后实测抛高值与该理论值往往存在偏差。因此可将实测抛高值与理论值的比值作为GM（1，1）模型的原始数据列。通过有限元分析得到的理论抛高值H =(h1，h2，…，hn)；节段完成后实际测得的抛高值模型原始数据列为λ=（λ1，λ2，…，λn），其中λi = h′i/hi，i =1，2，…，n。再利用上述GM（1，1）模型的灰色预测理论推出下一节段实测抛高与理论预抛高的比值λn + 1，进而得到经预测优化后的抛高为

以大桥2#墩悬臂浇筑施工过程为例，其单侧悬臂浇筑节段划分见图3。

图3 2#墩单侧悬臂浇筑节段划分（单位：cm）

部分节段有限元模型仿真分析结果（理论抛高值）与实测抛高值数据见表1。

表1 2#墩7—16号节段抛高值与λi的实测计算值

节段号λi λi 7 8 9 10 11 H/m 0.071 0.081 0.090 0.101 0.117 H′/m 0.077 0.085 0.095 0.101 0.117 1.071 1.055 1.050 0.994 1.000节段号12 13 14 15 16 H/m 0.133 0.151 0.169 0.187 0.207 H′/m 0.125 0.138 0.166 0.175 0.216 0.938 0.913 0.985 0.937 1.046

以7—10号节段理论抛高值与实测抛高值来预测11号节段的抛高值为例，过程为：以7—10号节段的λi为原始数据序列，其1-AGO 生成分别为1.071，2.126，3.176，4.170；1-AGO生成序列的紧邻均值生成分别为1.598 5，2.651 0，3.673 0；计算得a=0.029，b=1.110。计算过程的模拟值及误差见表2，计算平均模拟相对误差为1.111%。预测11号节段λi为0.974，反算其预测抛高值为0.114 m。

表2 λ11预测过程模拟值及误差

节段号8 9 10实测计算值1.055 1.050 0.994模拟值1.063 1.033 1.003误差-0.008 0.017-0.009相对模拟误差/%0.789 1.653 0.889

11 号节段施工完成后实测抛高值为0.117 m，误差仅为3 mm。

4.4 抛高值预测精度影响因素分析

4.4.1 样本数量对抛高值预测精度的影响

根据建模计算原理，当数据样本数量多于2 个时即可进行预测。对于大跨径刚构桥悬臂浇筑施工线形监控工作，因桥梁跨径大、悬臂长、节段多，越接近跨中的节段，其抛高值预测时可利用的数据样本越多，所以有必要分析样本数量对抛高值预测精度的影响。预测精度通过偏差率来反映：

可知偏差率越大，预测精度越低，反之则越高。

分别利用表1 中2 个样本数量（11—12 号节段λi值）、3 个样本数量（10—12 号节段λi 值）、4 个样本数量（9—12号节段λi值）、5个样本数量（8—12号节段λi值）、6 个样本数量（7—12 号节段λi 值）预测13 号节段的λ13值，分析不同样本数量的预测精度，结果见表3。

表3 不同样本数量下λ13的预测精度

样本数量实测计算值2 3 4 5 6预测值2.462 0.880 0.923 0.916 0.925 0.913偏差率/%169.660 3.614 1.095 0.329 1.314

从表3 可以看出：当样本数量≥4 个时，预测的λ13的偏差率基本稳定且预测精度很高，其值明显小于3个样本的预测值；偏差率不会随样本量的增加而减小，而是呈波动状态；2个样本的预测效果偏差过大。

4.4.2 预测步数对抛高值预测精度的影响

根据建模计算原理，GM（1，1）模型还可进行多步预测。对于大跨径刚构桥悬臂浇筑施工线形监控，应用GM（1，1）模型不仅能预测下一节段的抛高值，还可以预测之后任意节段的抛高值。

将实验结果与现有方法的结果进行对比,根据文献[5-6，8]中的数据,文献[5]采用单目视觉和2组点结构光的方法实现了纸箱三维尺寸的自动化测量,对真实纸箱的测量相对误差在0.20%到11.20%之间;文献[6]采用精确标定的单相机对长方体标准块的测量相对误差在2%以内,但要求手工确定特征点,没有达到本文自动化处理的水平;文献[8] 采用双目视觉并手工选择特征点的方法,对纸箱的测量相对误差在1.5%到11.7%之间。综合以上对比结果,本文方法的测量精度不低于现有最好的方法,在测量系统结构的简单性、自动化程度、鲁棒性等方面都具有优势。

利用表1 中7—12 号节段的λi 经过1～5 步分别预测13—17 号节段的λi，并与实测计算值进行比较，利用λi的偏差率判断预测精度，结果见表4。

表4 λi多步预测精度分析

λi λ13 λ14 λ15 λ16 λ17预测值0.925 0.900 0.875 0.851 0.827实测计算值0.913 0.985 0.937 1.046 1.054偏差率/%1.314 8.629 6.617 18.642 21.537

从表4可以看出，抛高值多步预测的效果较差，且预测精度随预测步数的增加基本呈降低趋势。对于大跨径刚构桥悬臂浇筑施工线形监控，当节段预应力张拉完成后便可获得当前节段的实测抛高值，因此实际施工时仅进行下一节段抛高值的预测即可。

5 结论

以大跨径连续刚构桥为研究对象，研究了GM（1，1）模型灰色预测模型在大跨径刚构桥悬臂浇筑施工线形监控中的应用技术。得到以下结论：

1）应用GM（1，1）模型进行大跨径刚构桥悬臂浇筑施工线形监控可以通过对抛高值的预测调整来实现，即结合各节段的理论抛高值，利用已完成节段的实测抛高值预测待进行节段的抛高值，并在挂篮立模时进行调整，且精度很高，能够满足施工监控的要求。

2）根据本文实例，采用灰色预测模型进行大跨径刚构桥悬臂浇筑施工线形监控时，按本文采取的主梁节段抛高值预测方法，当使用的样本数量≥4 个时，即能达到线形施工监控的精度要求，因此可采用已完成的4个邻近节段的实测抛高值预测待进行节段的抛高值，且一步预测精度较高。

参 考 文 献

［1］邓聚龙.灰色控制系统［J］.华中工学院学报，1982，10（3）：9-18.

［2］刘思峰，杨英杰，吴利丰，等.灰色系统理论及其应用［M］.北京：科学出版社，2014.

［3］刘震.灰色关联与预测建模技术及其应用研究［D］.南京：南京航空航天大学，2016.

［4］方志，刘光栋，王光炯.斜拉桥施工的灰色预测控制系统［J］.湖南大学学报（自然科学版），1997，24（3）：74-81.

［5］张劲泉，徐岳，鲜正洪.悬索桥主缆架设阶段灰色控制系统的研究［J］.西安公路交通大学学报，1997，17（4）：51-55.

［6］张永水，顾安邦.灰色系统理论在连续刚构桥施工控制中的应用［J］.公路，2001（6）：42-45．

［7］包仪军，王常峰，赵继康.基于改进灰色预测GM（1，1）模型的大跨度桥梁施工控制［J］.铁道建筑，2016，56（2）：18-22.

［8］杨正虎，尹紫红，彭文件，等.客运专线群桩基础桥梁沉降监测与预测方法研究［J］.铁道建筑，2012，52（9）：1-4.

［9］王明慧.用灰色关联分析法选择桥型设计方案［J］.铁道建筑，2012，52（9）：23-25.

［10］杨昀，周列茅，周勇军.弯桥与高墩［M］.北京：人民交通出版社，2011.