回归分析和灰色理论在地面沉降预测中精度对比.docVIP

回归分析和灰色理论在地面沉降预测中精度对比 　　摘要：本文将灰色系统理论、回归分析引入到因抽汲地下水而导致地面沉降发展趋势预测的研究之中。根据沧州市地面沉降的历史监测数据建立了沧州市地面沉降模型，运用两类不同的模型分别对该区地面沉降发展趋势进行了预测，比较得出灰色系统理论预测精度最大，为该地区今后的地面沉降发展趋势及其防治提供可靠依据。 　　关键词：地面沉降 灰色系统 回归分析 GM（1，1）模型 沧州 　　 　　0 引言 　　地面沉降灾害是环境地质灾害中非常严重的问题之一。主要表现在：地表高程资源损失、地面低洼变形、河道泄洪能力下降、洪涝灾害加重、工程设施破坏、建筑物基础下沉、风暴潮及海岸侵蚀灾害加重等。科学界在努力获得对自然沉降和人为沉降的进一步理解，同时在地面沉降的野外和室内检测、机理、预测方法、防治措施与效果方面也取得了丰硕的成果。目前，地面沉降发展趋势预测是地面沉降研究中的重要内容。 　　1 地面沉降发展趋势预测方法 　　目前，地面沉降预测计算的模型很多，主要有土水模型（多由水位预测模型和土层压密模型构成）和生命旋回模型（直接由沉降量和时间的关系构成）。 　　1.1 回归分析预测法 　　地面沉降灾害的发生与许多因素有着直接与间接的联系，比如地下水（石油、天然气等）的累计开采量、地下水位的埋深、地层岩性结构特征（包括粘性土的岩性、厚度、固结性、渗透性等）、地质构造等。科学实践表明，变量之间的关系可以分为确定性关系和相关关系。在地面沉降发展趋势的预测研究中，累计地面沉降量与上述许多变量之间虽然有非常密切的关系，但是要找出它们之间的确切关系是非常困难的，而且这种确定关系从理论上讲也是不存在的。回归分析就是应用数学方法，根据最小二乘法原理，对这些监测数据与实验数据去粗取精，去伪存真，从而得到反映地面沉降规律的数学预测模型。 　　1.2 灰色理论GM预测法 　　灰色系统理论用于预测多采用GM(1，1)模型，该模型是单变量一维模型。GM(1，1)模型不使用待预测量的原始数据进行处理，它在原始序列的基础上生成一个新的数据序列，使用新数列进行预测，然后对预测数列进行还原处理，得到真正的预测数据序列。灰色系统理论克服了其他理论和方法（回归分析、方差分析）要求大样本量、样本要有较好的分布规律、计算工作量大等方面的缺点，很好的保证了预测的精度。 　　不同因素所诱发的地面沉降的范围、速度及持续时间不同。因此，地面沉降研究的对象是“土层―地质作用―人类活动”的复杂系统。由于表征该系统的参数是灰数、影响地面沉降的因素是灰元、构成系统的各种关系是灰关系，所以该系统是一个不具有物理原型的灰色系统，灰色系统理论在地面沉降研究领域应用是有科学依据的。 　　2 地面沉降发展趋势预测的实例计算 　　2.1 实验数据 　　该区域松散沉积物厚度近1000m，按照水文地质条件与含水层特征，其间被分成五个含水组（层）：第Ⅰ含水层埋深0~50m，通过降雨入渗补给；第Ⅱ含水组埋深30~110m之间，因矿化度高开发成本较高，基本没有利用；第Ⅲ含水层埋深150~350m，主要有粘性土与砂性土组成，该层地下水开采量约占总开采量的90%；第Ⅳ含水层埋深350~450m；第Ⅴ含水层埋深450~1000m。由于沉降层的埋深厚度为160~348m(主要是第Ⅲ含水组)，粘性土约占此厚度的72%，其中以粉质粘土为主，达110m左右；砂性土约占沉降层厚度的28%，其中细砂占砂性土厚度的45%，粉细砂占32%，中细砂占12%，粉砂占11%。因此，本文在选取数据时主要考虑第Ⅲ含水组的情况。 　　观测数据表明，上世纪70年代沧州漏斗中心累计沉降量较小(200mm)，80年代后期沉降量迅速增加，到1990年漏斗中心累计沉降量达到1130mm，地面沉降灾害已成为沧州地区主要地质灾害之一。 　　2.2 预测结果 　　2.2.1 线性回归预测结果 　　回归预测包含线性回归预测与非线性回归预测。为了考虑数据采样的时序性和连续性，在这里选取1978~1983年连续数据进行预测。根据上述原则选取数据，在直角坐标系中对应的做出一系列的点，可得散点图1。 　　根据图1中图形大致走向，我们引用线性回归、乘幂回归、指数回归、对数回归进行预测，部分预测结果见表1和图2（由于计算量大，本文回归预测及图形绘制均使用Excel进行，其中图2进行了合并处理）。 　　通过计算发现对数回归预测模型的相关系数平方值R2= 0.913，精度相对低，而指数回归预测残差太大，所以将其预测结果和预测模型剔除。 　　2.2.2 GM（1，1）模型预测结果 　　①GM(1，1)建模 　　根据沧州地区历史漏斗中心累计地面沉降量实测数据建立GM(1，1)模型，并在此基础上对该地区地面沉降量进行预测