论文部分内容阅读
随着大规模的城市建设及过量的地下水抽取,导致了城市地下水位的下降,从而引起严重的地面沉降与地表塌陷等一系列地质灾害的发生。加强城市地表沉陷监测和建立地下水-地面沉降预测模型是预防地表沉陷灾害发生有效手段。利用InSAR(Interferometric Synthetic Aperture Radar,InSAR)技术对大尺度的城市地面沉降进行监测,可以得到大面积、高分辨率和高精度的持续的时序沉降监测结果,为地面沉降预测模型参数的反演提供大量的实测数据。本文系统分析了常规InSAR与时序InSAR技术原理、处理流程以及主要误差源与敏感度。针对时序InSAR地面沉降监测技术中的大气信号分离问题,研究了常规时空滤波与平滑样条滤波两种大气信号分离算法,分别采用常规时空滤波方法与平滑样条滤波方法进行大气滤波,实现大气信号与形变信息的分离,并对两种滤波方法得到的大气与形变结果进行对比与分析,通过对滤波算法的实验评价两种滤波对大气与形变信号分离的精度与抗差程度。针对地下水-地面沉降预测模型参数反演问题,对模型参数与沉降值的相关性进行了分析,采用InSAR地面沉降监测数据约束下的LS(Least Square,LS)反演方法得到模型中的固结系数;采用连续小波变换方法分离出地面沉降量中因含水层压缩变形引起的地表弹性与非弹性沉降量,并结合地面沉降预测模型反演得到弹性释水系数和非弹性释水系数;利用InSAR地面沉降监测值,根据Biot固结理论中的自由应力下体积应变公式反演得到模型的地下水体积变化量参数。以天津市InSAR地面沉降监测为例,采用短基线集InSAR(Small Baseline InSAR,SBAS-InSAR)方法,获取了地面沉降变形值,并基于PS点的地面沉降值解译了地下水位、地下水系统、基地构造等相关水文地质信息,通过InSAR地面沉降值反演了地面沉降模型参数,建立了天津市地下水-地面沉降预测模型,采用Modflow地面沉降预测数值模拟软件模拟计算,将模拟计算结果与地面沉降测量值进行比较,验证了模型参数反演方法的正确性。实验表明地面沉降误差小于10mm,最大残差15mm。主要研究成果如下:(1)研究了InSAR时序处理过程中常规时空滤波与平滑样条滤波两种大气信号分离算法,分别采用常规时空滤波方法与平滑样条滤波方法进行大气滤波,实现大气信号与形变信息的分离,并对两种滤波方法得到的大气与形变结果进行对比与分析,通过实验得出:在非均匀采样状态下常规滤波方法滤波参数的固定使滤波结果存在着不确定性,平滑样条滤波也可以有效的分离大气与沉降信号,该滤波方法通过选取最佳平滑系数,能够较好的逼近形变相位分量。(2)通过研究地面沉降模型建立与参数反演理论,分析地下水三维渗流-地面沉降预测模型的建立过程,探讨了渗流偏微分方程定解条件及耦合后的收敛性能,并对模型中的几个重要模型参数如地下水位、弹性释水系数、非弹性释水系数、固结系数、开采量与沉降的关系进行参数相关性分析,解释每个参数对沉降模型的主要贡献度与灵敏度。(3)提出了一种含水层弹性、非弹性释水系数的InSAR反演方法。研究承压含水层形变机理,结合垂直时序InSAR形变分量和水位变化因素,采用连续小波变换(Continuous Wavelet Transform,CWT)方法分离周期形变信号分量和长期趋势,分别得到含水层弹性和非弹性时序形变信息,孤立的信号分量基于地面沉降模型用于估计弹性储存系数、非弹性骨架释水系数和压实时间常数,从而可以得到精确的弹性、非弹性释水系数,解决了以前靠InSAR只能反演平均释水系数一个参数的困难,为后续的水文地质研究提供更准确的参数依据。采用最小二乘方法,建立了InSAR时序数据源下的固结系数约束方程,提出基于InSAR时序技术的含水层固结系数反演方法。结合Biot固结理论中的自由应力下体积应变模型,讨论不同水文地质概念模型下InSAR对地下水探测的敏感性,提出采用卫星视角分量矢量化的方法将其垂直分量融入Biot模型中,从而达到体积应变积分方程求解收敛、地下水体积变化率反演的目的。(4)为了验证本文提出的基于InSAR技术的地下水-地面沉降预测模型参数反演方法的准确性与实用性,以天津市沉降区为例,采用SBAS-InSAR技术对主要沉降大区域与单个漏斗区进行监测,并采用定性与定量的方法作了基于PS点的地面沉降特征、地下水位、地下水系统、基地构造等相关水位地质解译。基于上述提出的反演方法反演了地下水-地面沉降预测模型所需的弹性释水系数、非弹性释水系数、开采量变化和固结系数,并利用Modflow三维地下水流软件,Interbed-Storage程序包模拟地面沉降,先用初始水位地质实测单点参数代入模型作为初始参数模拟得到一个与InSAR时间内对应的沉降值序列,再利用InSAR反演得到的参数代入模型得到一个与InSAR时间内对应的沉降值序列并与实测水准、InSAR进行精度对比验证地下水-地面沉降预测模型的可靠性,最后提出基于赤池信息准则(Akaike’s Information Criterion,AIC)下的地面沉降模型参数、与沉降值精度综合评价方法。