联合MT-InSAR技术和MGWR模型的分层含水层组储水系数反演方法技术

本发明专利技术公开一种联合MT

【技术实现步骤摘要】
联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法


[0001]本专利技术涉及一种含水层组储水系数反演的方法，具体而言，涉及一种联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法，特别是针对主要因为地下水开采导致的地面沉降区，联合MT
‑
InSAR技术和空间统计分析方法，量化分层含水层水位变化与地面沉降相关性，根据分层含水层水位变化对地面沉降的贡献度，将沉降数据进行分层，并结合地下水位监测数据，在空间尺度反演不同层位含水层组的储水系数。

技术介绍

[0002]地面沉降是指在一定的地表面积内所发生的地表海拔标高降低的现象，是一种缓变性的地质灾害，主要是由于过量开采地下水，导致地下水位降低，从而产生地面沉降的现象，造成含水层储水能力损失。不均匀的地面沉降会影响和制约当地经济建设可持续发展和社会安定，成为现代城市的重要安全隐患。因此，为了管理地下水资源的可持续利用和控制地面沉降，反演不同层位含水层组储水系数，了解其储水特性至关重要。
[0003]现有的研究估算储水系数的思路大致可以分为两类。第一类是不考虑复杂含水层系统的分层特性，将含水层组作为一个整体来研究。例如Riley等用图形方法，通过绘制应力应变曲线图，根据曲线的反斜率来估计储水系数；胡谢等将InSAR数据和地下水位数据相结合，用谐波序列方法，通过线性回归反演储水系数；白林等使用多通道奇异谱分析方法，分离季节性变形和季节性水头变化，考虑滞后用最小二乘法反演储水系数。然而不同层位的开采情况不一样，含水层组的压缩情况不同，储水系数也不同。第二类是考虑复杂含水层系统的分层特性，利用分层标数据结合水位数据估算不同层位含水层组的储水系数。例如李江涛等依据分层标，用快速独立分量分析(Fast
‑
ICA)与可变预固结水头分解方法相结合的方法，估算了分层含水层的储水系数。雷坤超等利用分层标监测站，根据预固结水头判断并估算北京平原区不同压缩层组储水系数。这类研究需要分层标数据，其实验成本高，监测点位少，仅仅能反映个别点位的储水系数情况，难以大范围推广。
[0004]传统的大地测量技术(如GPS、水准测量等)空间分辨率低、覆盖范围小。多时相合成孔径雷达干涉测量(Multi
‑
temporal Synthetic aperture radar interferometry,MT
‑
InSAR)技术可以同时获取地表对微波的反射强度和相位信息，与传统测量技术相比，具有全天时、全天候、测量范围广等特点，可以在大规模尺度上监测地表变形。MT
‑
InSAR技术利用同一地区的多景SAR影像对时序稳定点(persistent scatterer,PS)进行精确分析，极大地降低了大气延迟等带来的测量误差，使得形变监测精度达到了厘米级到毫米级，在角反射器的辅助下甚至可实现亚毫米级监测精度，相对于传统测量方法，可提取城市区域大范围、高精度的地表三维信息和形变信息。MT
‑
InSAR发展到现在比较有代表性的是以单一影像为主影像的永久散射体干涉(persistent scatterer InSAR，PS
‑
InSAR)方法和以多幅影像为主影像的小基线(small baseline subsets,SBAS)方法。
[0005]在空间统计分析方法中，MGWR模型描述了自变量与因变量之间的关系，还考虑了
数据的空间异质性。将数据的地理位置考虑到回归参数之中，能够体现变量随空间位置的不同而出现差异。允许每个自变量拥有各自不同的空间平滑水平，解决了GWR模型所有变量同一平滑水平的缺陷，降低了估计的偏误。通过MGWR模型可以量化分层含水层地下水位变化与地面沉降的相关性，进而获得各层含水层对地面沉降的贡献率。
[0006]综上所述，以往估算储水系数的方法，大多在少数点位进行，或者将含水层组看作整体计算，将结果视作整个含水层组共用的储水系数。随着全球地面沉降的发育，研究出一种能够在地面沉降区，反演不同层位含水层组区域尺度储水系数，以了解其储水特性的方法至关重要。

技术实现思路

[0007]为解决上述问题，本专利技术的目的在于提供一种联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法，适用于空间尺度且不同层位的储水系数反演。该方法联合MT
‑
InSAR技术和MGWR模型，量化分层含水层水位变化与地面沉降相关性，获得各层含水层对地面沉降的贡献率，根据分层含水层水位变化对地面沉降的贡献度，在空间尺度反演得到不同层位含水层组的储水系数。
[0008]为达到上述目的，本专利技术提供了一种联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法，具体包括以下步骤：
[0009]步骤S1：利用MT
‑
InSAR处理SAR影像，获取地面沉降信息，所述地面沉降信息包括地面沉降数据和地下水位数据；
[0010]步骤S2：利用ArcGIS软件的克里金插值法将需要输入MGWR模型的沉降数据和地下水位数据进行预处理；
[0011]步骤S3：构建MGWR模型，并将预处理后的数据输入模型，获取分层含水层地下水位变化与地面沉降的相关性；
[0012]步骤S4：基于分层含水层地下水位变化与地面沉降的相关性，计算分层含水层对地面沉降的贡献率；
[0013]步骤S5：根据得到的贡献率将MT
‑
InSAR获取的地面沉降信息分层，具体为将步骤S1中获取的地面沉降信息分别乘步骤S4计算得到的分层含水层对地面沉降的贡献率，得到不同层位含水层组的地面沉降信息；
[0014]步骤S6：将分层的地面沉降信息结合各层位地下水位变化数据在空间尺度反演得到不同层位含水层组的储水系数。
[0015]在本专利技术一实施例中，其中，步骤S1具体包括：
[0016]步骤S101：设SAR影像共有N+1幅，选取其中一幅作为主影像，将剩余的N幅影像作为辅影像与主影像进行配准；
[0017]步骤S102：利用外部获取的DEM数据模拟地形相位，然后将辅影像与主影像进行差分干涉以去除地形相位和平地相位；确定基线对，并采用确定的基线组合再次进行差分干涉处理；
[0018]步骤S103：通过相关系数阈值法或相位离差阈值法选取永久散射体；
[0019]步骤S104：进行相位解缠和噪声去除处理，其中，所述噪声包括大气误差及DEM引入的地形误差；
[0020]步骤S105：去除轨道误差并生成时间序列形变信息，其中轨道误差是指垂直轨道分量；
[0021]其中，在步骤S102进行干涉处理后能够得到N幅干涉相位图，然后通过式(1)对干涉相位进行分解：
[0022][0023]其中，为干涉相位；为平地相位，其通过在读取SAR数据时利用卫星精确轨道数据去除；为地形相位，其是本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法，其特征在于，包括以下步骤：步骤S1：利用MT
‑
InSAR处理SAR影像，获取地面沉降信息，所述地面沉降信息包括地面沉降数据和地下水位数据；步骤S2：利用ArcGIS软件的克里金插值法将需要输入MGWR模型的沉降数据和地下水位数据进行预处理；步骤S3：构建MGWR模型，并将预处理后的数据输入模型，获取分层含水层地下水位变化与地面沉降的相关性；步骤S4：基于分层含水层地下水位变化与地面沉降的相关性，计算分层含水层对地面沉降的贡献率；步骤S5：根据得到的贡献率将MT
‑
InSAR获取的地面沉降信息分层，具体为将步骤S1中获取的地面沉降信息分别乘步骤S4计算得到的分层含水层对地面沉降的贡献率，得到不同层位含水层组的地面沉降信息；步骤S6：将分层的地面沉降信息结合各层位地下水位变化数据在空间尺度反演得到不同层位含水层组的储水系数。2.根据权利要求1所述的联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法，其特征在于，步骤S1具体包括：步骤S101：设SAR影像共有N+1幅，选取其中一幅作为主影像，将剩余的N幅影像作为辅影像与主影像进行配准；步骤S102：利用外部获取的DEM数据模拟地形相位，然后将辅影像与主影像进行差分干涉以去除地形相位和平地相位；确定基线对，并采用确定的基线组合再次进行差分干涉处理；步骤S103：通过相关系数阈值法或相位离差阈值法选取永久散射体；步骤S104：进行相位解缠和噪声去除处理，其中，所述噪声包括大气误差及DEM引入的地形误差；步骤S105：去除轨道误差并生成时间序列形变信息，其中轨道误差是指垂直轨道分量；其中，在步骤S102进行干涉处理后能够得到N幅干涉相位图，然后通过式(1)对干涉相位进行分解：其中，为干涉相位；为平地相位，其通过在读取SAR数据时利用卫星精确轨道数据去除；为地形相位，其是利用从美国SRTM获取的DEM数据去除；为变形相位；为大气相位，其通过SARPROZ软件中的APS处理消除；为APS处理过程中的热噪声和配准误差，其通过线性模型消除；在消除和之后，得到雷达视距方向的变形相位然后将雷达视距方向的变形相位通过式(2)转化为垂直位移d
v
，即得到步骤S105的垂直轨道分量：d
v
＝d
los
/cosθ
ꢀꢀꢀꢀ
(2)
其中，d
los
为雷达视距方向的变形相位，θ为入射角。3.根据权利要求1所述的联合MT
‑
InSAR技术和MGWR模型的分层含水层组储水系数反演方法，其特征在于，步骤S2是利用ArcGIS软件的克里金插值法将需要获取到的基于点的水位数据插值成面数据，需要计算每一层地下水位数据的变化量，具体为：利用ArcGIS软件中的创建渔网工具建立1km
×
1km的格网，然后将各层地下水位变化数据和沉降数据提取均值到每个格网中，以作为...

【专利技术属性】
技术研发人员：陈蓓蓓，宫辉力，赵单妮，周超凡，雷坤超，胡锦铭，
申请(专利权)人：首都师范大学，
类型：发明
国别省市：