1.面向全运行周期的矿区地表形变实时/准实时监测方法,该方法的步骤如下:针对空基/天基/地基对地观测技术获取的面状、线状、点状异构数据融合,应用区域性椭球理论构建矿区地表形变监测的系统框架与转换模型;根据空间尺度理论统一异质数据精度尺度,结合两种典型平差模型的几何意义分析不同几何尺度对模型稳健性、几何尺度不一致对异质数据融合的影响,建立数据处理模型一致性评价手段;应用方差分量估计理论,建立异质数据融合的平差模型;根据矿区地表形变的沉降-失衡-塌陷的基本规律,离散监测区域形成时/空格网,结合不同观测技术获取数据的时/空分辨率互补性对缺失监测区域/监测时段进行插值,应用序列模型对不同时域、不同空域的监测数据进行融合,实现覆盖矿区基建-开采-闭坑全运行周期的地表形变实时/准实时自主化监测;其特征是该方法的具体步骤如下:Step1.构建面向多种对地观测技术获取异构数据融合的矿区参考框架;其中,根据矿区的平均高程面确定监测区域的参考椭球,椭球几何参数选择国家坐标系统椭球参数或者地方坐标系统椭球参考,以经过监测区域中心位置的子午线为坐标投影的中央子午线,投影带度数根据矿区区域确定,其跨度至少覆盖形变监测区域;应用具有大旋转角的Bursa模型建立坐标系统转换模型,构建监测区域统一参考框架;转换模型中的坐标观测值应进行重心化处理,并进行病态性检测;天基/空基获取的遥感影像或者雷达干涉影像干涉点,以及地基三维激光扫描的测站点选择,应优先选用地表已经布设的监测点,若监测点缺乏,则选用矿区不局限于井筒的需要进行形变监测的重要构筑物的特征点,进行地面三维控制测量;控制测量的平面精度不低于《全球定位系统测量规范》(GB/T18314-2009)中B级GPS测量的精度要求,高程精度不低于《国家一、二等水准测量规范》(GB/T12897-2006)中二等水准测量的精度要求;
Step2.Msplit估计理论与S转换模型对基于多种观测技术的异质监测数据一致性检测;其中,异质数据一致性检测主要针对不同观测技术获取的监测数据,包括矿区地表形变常用监测手段获取的监测数据,也包括新兴监测手段;结合矿区地表形变监测的基本规律,定义粗差为一致性检测的指标,将其分为两类情况: (1)同一监测点有两类或两期及以上监测数据,应用Msplit估计理论与S转换模型对监测数据进行粗差的检测-识别-调整;其中,粗差检测的目标函数由Msplit构建,粗差的识别由t分布构建统计模型,粗差的调整由IGG权函数构建模型;(2)同一监测点仅有一类一期监测数据,围绕平差模型的内部可靠性,根据平差模型得到的单位权方差,应用t分布构建统计模型对监测数据的粗差进行检测;矿区地表形变监测是通过多期监测数据的对比分析,揭示其形变过程与形变规律,因此,一类一期监测数据的粗差检测是对原始观测数据质量的初步检核,随着监测数据量的积累,这类数据须纳入到第一类情况处理;其中,常用监测手段不局限于GNSS、导线测量、水准测量;新兴监测手段不局限于雷达干涉测量、三维激光扫描;粗差的识别不局限于t分布构建的统计模型,粗差的调整不局限于IGG权函数构建的模型;
Step3.空间尺度理论定义函数模型几何空间尺度进行模型可靠性检测;其中,针对矿区形变特征,判别、选择合适的数据融合模型,至少对线性与非线性这两类模型进行比较;根据模型的几何意义,定义模型的几何空间尺度,度量不同尺度标准对模型可靠性的影响与可靠性模型的识别;基于观测误差服从正态分布的假设,应用F分布构建统计模型检测不同尺度标准对模型可靠性影响;根据模型几何意义定义的几何空间尺度为指标,度量模型的可靠性,选择适用于所监测矿区的数据处理模型;
Step4.方差分量估计理论建立异质数据融合模型,基于地表形变状态域的监测数据一致性检测;其中,针对不同观测技术获取异质数据具有的异方差性,应用方差分量估计理论对具有异方差的监测数据进行融合;同时,须对方差分量估计构建的平差模型进行病态性分析:(1)当观测数据多于两种,即单位权中误差因子数量大于二时,应用量化的条件数指标对方差分量估计模型进行病态性检测,当方差分量估计模型病态,应用Tikhonov函数对病态模型进行正则化;(2)当方差分量估计模型含有两类观测数据时,不对模型进行病态性检测;应用监测数据构建矿区地表形变时间序列模型,根据序列模型预测形变值,对需要融合的监测数据进行一致性检测;具体方法为,应用预测值与监测值的差值构建统计量,应用t分布构建统计函数,进行融合数据的一致性检测;
Step5.监测区域分别进行时/空格网离散化处理,对缺失监测点的格网单元进行插值;其中,依据矿区已实施监测的时间与空间分辨率,对矿区形变监测区域进行格网处理,格网的空间分辨率根据矿区地表的地貌特征进行划分,总体原则是地貌越复杂格网越精细,并且其分辨率不得小于5m╳5m;时间分辨率根据矿区基建-开采-闭坑的运行周期进行划分,执行标准为《煤矿测量规程》、《工程测量规范》国家规范;格网的空间形状采用三角网/四边网,但至少应确保每个网格单元有一个监测点;若网格缺失监测点形变数据,则通过邻近区域已有监测数据构建形变模型,对形变模型应用克里金插值算法进行插值处理,推估监测点形变数据;
Step6.监测区格网平滑,矿区地表形变四维建模,实现矿区地表形变实时/准实时动态自主监测;其中,以单个格网为单元,应用多面函数构建矿区地表形变四维模型;在形变模型的构建中,单元格网的权值由格网内的监测点一致性检测结果作为先验信息确定;形变模型的参数是空间与时间构成的四维自变量,地表形变值为应变量;根据矿区地表形变的基本特征,自变量中的时间维度(t)以日为计量单位;空间维度采用空间三维直角坐标(X YZ)或者平面坐标加高程(x y H),以米为计量单位,监测点平面坐标由其空间三维直角坐标经定义的矿区区域性椭球投影变换得到;地表形变值以毫米为计量单位。
2.根据权利要求1所述的面向全运行周期的矿区地表形变实时/准实时监测方法,其特征是:Step3中,线性模型不局限于G-M模型,非线性模型不局限于G-H模型;应用不局限于F分布构建的统计模型检测不同尺度标准对模型可靠性影响。