第46卷第5期 中国矿业大学学报 Vo1.46 No.5 2017年9月 Journal of China University of Mining&Technology Sep.2017 基于D—InSAR LOS向变形的 开采沉陷预计参数反演方法研究 王 磊 ,张鲜妮 ,陈元非 ,查剑锋 (1.安徽理工大学测绘学院,安徽淮南232001; 2.中国矿业大学国土环境与灾害监测国家测绘地理信息局重点实验室,江苏徐卅I 221116) 摘要:针对基于常规D-InSAR技术成果无法反演全部概率积分参数(特别是水平移动参数)问 题,开展了利用D-InSAR LOS向变形基于模矢法全面反演概率积分参数问题研究:1)依据I)_ InSAR LOS向变形为下沉和南北、东西方向水平移动沿LoS向投影的关系,利用模矢法非线性 规划理论,构建了基于单视线向D—InSAR LOS向变形的概率积分参数反演模型(可反演出全部 概率积分参数),并编制了求参软件;模拟实验表明,q,tan ,b,0求参相对误差均小于2.5 ,求 参拟合误差在一l1~11 mm之间,LOS向变形拟合中误差为m一_-4-5.8 mm,整体求参效果较 好.2)利用建立的参数反演模型,在一定假设条件下,对2012年3月11日9310工作面开采沉 陷预计参数进行了求取,即q===0.170,b一0.12,tan J9—2.12,0—88。,S1一一11,S2—10,S3一 一9O,S4—33,拟舍误差约在一34 ̄38 mm之间,拟合中误差m===±14.67 mm;并认为S3一一90 主要为9310工作面开采导致临近9312采空区边界的悬臂梁或砌体梁失稳“活化”所致.研究成 果对D-InSAR矿山变形监测具有重要参考价值. 关键词:D-InSAR;LOS向移动;模矢法;概率积分参数 中图分类号:TD 32 文献标志码:A 文章编号:1000—1964(2017)05—1159—07 M ethod of mining subsidence prediction parameters inversion based on D—InSAR LOS deformation WANG Lei ,ZHANG Xianni ,CHEN Yuanfei ,ZHA Jianfeng (1.School of Geodesy and Geomatics,Anhui University of Science and Technology,Huainan, Anhui 232001,China;2.Key Laboratory for Land Environment and Disaster Monitoring of SBSM, China University of Mining&Technology,Xuzhou,Jiangsu 221116,China) Abstract:Because traditional D-InSAR technology cannot calculate all parameters(especially horizontal move parameters)in probability integral method parameters inversion,this article studied the method of probability integral method parameters overall inversion based on de— formation in D-InSAR LOS direction as well as pattern search method.According to the theory of D-InSAR technology,in which the deformation in LOS direction reflects the relationship of the proj ection of vertical and horizontal movement on LOS direction,the probability integral method parameters inversion model including all parameters based on single‘—LOS D-InSAR de‘・ formation was established by pattern search method for nonlinear programming theory.The parameter calculation software was also compiled.Resuhs of the simulation experiment show 收稿日期:2016一O6一O5 基金项目:国家自然科学基金项目(41602357,41474026);安徽省博士后基金项目(2014B019);安徽高校自然科学研究项目 (KJ2016A190) 通信作者:王磊(1984一),男,安徽省合肥市人,副教授,工学博士(博士后),从事矿山变形监测及控制方面的研究. E-mail:austwlei@163.corn Tel:18098693642 中国矿业大学学报 第46卷 that the relative errors of q,tan口,b and 0 are all smaller than 2.5 ;the fitting error of pa— rameters calculation iS between一1 1 and 1 1 mm:the fitting error of L0s_direction deformation m一±5.8 mm.The probability integral method parameters on March 1 1 of 2012 of 9310 work face were caculated with this model,and the results are:g一0.170,b一0.12,tan 一2.12, 一 88。,S1一一11,S2—10,S3一一90,S4—33;the fitting error is between一34 and 38 mm;the mean square fitting error 一±14.67 mm.The reason of S3一一90 is that the mining of 9310 work face leads to the activation of cantilever beam or voussoir beam of 9 3 1 2 gob edge.The re— sults of this article have important referential value in the application of D-InSAR mining de— formation monitoring. Key words:D—InSAR;LOS-direction movement;pattern search method;probability integral parameter 开采沉陷预计对建筑物下、铁路下、水体下(简 称“三下”)采煤实践具有重要指导意义.概率积分 法是我国“三下”采煤规程 指定的开采沉陷预计 究表明,基于单视向I)-InsAR技术虽然无法直接监 测开采沉陷三维移动变形,但LOS向移动变形本质 为目标点的下沉、水平移动分量沿L0S向投影之 方法,因而概率积分参数的求取也是矿山变形监测 数据处理的核心内容.D-InSAR由于具有全天时、 和,理论上可根据这种关系模型,结合目标区域地质 采矿条件,利用非线性规划理论反演概率积分参数. 鉴于此,本文针对途径2开展重点研究,结合 采动变形规律和雷达几何成像特点,拟采用模矢法 构建基于单视线向I)_InSAR的开采沉陷概率积分 参数反演模型,并通过模拟实验开展可行性研究, 高精度、似面测量等特点,已经成为当前开采沉陷 监测的研究热点,但当前利用单视线向D-InSAR 技术主要用于监测开采沉陷下沉(低人射角LOS 的D-InSAR监测仅对下沉敏感),无法提取水平移 动.因而,目前基于单视线向D-InSAR监测成果仅 能反演与下沉相关的概率积分参数 (下沉系数 q、最大下沉角 、主要影响角正切tan 、拐点偏移 距s),而并不能求取与水平移动相关的概率积分 参数(水平移动系数6),这严重的制约了D-InSAR 技术在矿山变形监测中的推广及应用. 解决上述难题,可以从两个途径去突破:1)建 立适用于开采沉陷特点的InSAR三维变形监测方 法,然后进行全面求参;2)仍基于单视线向D-In— SAR技术,寻找开采沉陷下沉、水平移动与LOS向 移动变形关系模型,基于非线性理论实施求参.对于 途径1通常解决方案主要有多视线向D-InSAR[。_5], D-InSAR+Offset trackingc 一 ,D-InSAR 4-MAI[ 、 最后将研究成果应用到兖州矿区南屯煤矿9310工 作面概率积分参数反演.研究成果对D—InSAR矿 山变形监测具有重要参考价值. 1 基于InSAR的概率积分参数反演方法 1.1 预计误差函数及求参准则 当卫星回访经过目标沉陷盆地上方成像时刻, 假设移动盆地内任意象元i的累计移动变形值为 (]S(可根据各期LOS向变形累加得到),基于概率 积分法模型和先验参数(经验初值)预计任意象元i 的下沉值 、北方向水平移动 sN和东方向水平移 动UEw.根据 InSAR一维LOS投影原理[15](如图 1所示,其中箭头方向表示正方向),则任意i像元沿 LOS向的预计移动变形值 ◆ 、 融合GPS和D-InSAR法l_9_ 1ol、地表先验模型+I)_ InSARc”。 ,但上述方法主要用于解决地震、滑坡、 火山喷发、冰川移动等特点的移动变形问题,然而上 述对象的移动机理不同于矿山开采沉陷,并不能妥 可表示为 ’ . 善解决矿山开采地表的三维变形监测问题.文献 [13]结合常规开采沉陷模型,利用水平移动与倾斜 的关系,提出了基于单视线向D-InSAR的三维变形 提取方法,但该方法仅适用于水平和缓倾斜煤层开 采沉陷监测,对倾斜与极倾斜煤层、山区等特殊地质 采矿条件并不适用.因而目前依据途径1进行概率 方位 - 东方向 图1 D-InSAR的LOS一维投影原理 Fig.1 LOS one-dimensional projection principle diagram of D-InSAR 积分求参,仍需要进一步开展研究.对于途径2,研 第5期 王 磊等:基于D-InSAR LOS向变形的开采沉陷预计参数反演方法研究 r,l_os一一u sin 01cos( 一—手7c)一 U rwsin 。sin(a 一号丌)+w cos九,(1) 1 为左拐点偏移距;S 为右拐点偏移距;S。为上拐点 偏移距;S 为下拐点偏移距;q为下沉系数;tan卢 为主要影响角正切;0为最大下沉角;b为水平移动 系数. 2)根据地质采矿条件与概率积分参数的经验 W —W(x )一 W。(so )w。( t), (2) UfsN一己,(z , {, sN)一 1 关系模型,选择初始参数近似点B ,并以此为起点 + ( t) 。(M)cos (基点)进行搜索, B1一(S1o,S20,S3o,S40,q0,tan ,0o,bo) 一 U。( )w。(z )sin 酬], (3) (X1o,X2o,X30,X4o,X5o,X6o,X70,X8o) .(8) 【,tEw—U(z , , EW)一 [【,。(z )w。( )co 删+ U。( )w。(z )sin ̄Oirw], (4) 式中: 为入射角 为卫星飞行方向方位角;W 为预计下沉值;U ,U删为北、东方向预计水平移 动;w。为充分采动情况下地表最大下沉值; ,  ̄oirw分别为北、东方向顺时针旋转到预计坐标系z 轴形成的角度;W。( )为倾向方向充分采动时走 向主断面上横坐标为 的点的下沉值;W。( 。)为 走向方向充分采动时倾向主断面上横坐标为 的 点的下沉值;U。( )为倾向方向充分采动时走向主 断面上横坐标为z。的点的水平移动值;U。( )为 走向方向充分采动时倾向主断面上横坐标为 的 点的水平移动值. 定义该像元的LOS向预计移动残差为 —r儿os—r os. (5) 构建求参准则为 ∑ li=1 —min. (6) 式(6)求参准则实质要求寻找一组概率积分参 数,使得参与求参的所有像元LOS向移动预计值 与实测值之差绝对值之和最小,该问题可利用非线 性规划理论对概率积分参数进行寻优. 1.2 融合LOS变形与模矢法的概率积分参数反 演模型 在非线性规划理论中,模矢法已经被广泛地应 用于传统地表移动观测站求参口 .本文在此基础 上,结合上述建立的预计残差模型和求参准则,利 用模矢法进一步构建基于单视线向D-InSAR的概 率积分法求参模型,模型求参的详细思路如下: 1)构筑预计误差函数,根据式(6)令 e(B)一∑I I, (7) i一1 B(S1,S2,S3,S4,q,tan卢,0,6) 式中:e(.B)为预计误差函数;B为概率积分参数;S 3)确定各分量(参数)探索步长,若任意分量 (参数)x ( 一1,2,…,8)步长为厶 ,即 △ 一(O,…,△X ,…,0) . (9) 4)计算初始基点曰 的误差函数值e(B ),比 较B +厶 点,若£(B +厶 )<£(B1),就以点B + △ 作为临时矢点,并记作 . 的第一个下标表 示建立第一个模矢,第二个下标表示x 已被摄 动.若£(B +△ )>e(B ),就试验B 一 点,如果 e(B --zl )<£(B ),就以B 一 点为临时矢点.否 则,仍以B 为临时矢点,上述搜索可用下式进行描 述 r.B +△1 e(B1+厶1)<E(B1), l B △1 e(B1一△1)<£(B1) T11一 <£(Bl+厶1), (10) l B e(B1)<minE ̄(B1一△1), e(Bl+△1)]. 若第一个分量(参数)被摄动,则用Tl 代替原 来的基点B ,对于下一个变量X。进行类似的 摄动,一般情况有 f 』+ +1 e( +厶,+1)<e(T1 ), J一 +1,£( --zlj+1)<e(T1J) ,J+1一{ <e(T1J+ +1), l£( ,)<minE ̄(T1J--zl,+1), l e( +A件1)]. (11) T1o—B1, ( 一0,1,…,7). 上述8个参数(矢量)都被摄动之后,得临时矢 点 。,并令7"1 一B。,由原始基点B 和新基点B 构成第一个模矢. 5)将第一个模矢延长一倍,得第二个模矢的 初始临时矢点 。一B1+2(B2一B1)一2B2一B . (12) 6)在T2。附近进行和上面类似的探索,建立临 时矢点 , ,…, 。,并以 。为第3个基点B。, 中国矿业大学学报 第46卷 此时由B ,B。就确立了第2个模矢.则第3个模矢 的初始临时矢点为 0一B2+2(B3一B2)一2B3一B2. (13) 件. 1.3求参模型应用方法 首先获取覆盖地表移动周期的SAR影像资 料,并按照时间以1,2,…,k进行编号,然后按照时 7)重复上述1)~6)过程,若对于第i个模矢 存在 e( o)<e(.B ). (14) 间顺序分别以1,2,…,k一1影像为主影像,对应的 以2,3,…,k影像为从影像,进行连续差分干涉测 量,并提取各时段D-InSAR LOS向移动变形;然 后对LOS向移动变形进行累加,获取目标区域采 动地表1~k时段累计D-InSAR LOS向移动变 则由 产生不出比B 更好的点,则应退回到 B ,模矢不延长,并在B 附近进行搜索,如果能得 出新的下降点,即可引出新的模矢;否则,将步长缩 小一半,开展更精细的探查.当步长缩小到足够高 的精度时,即可停止迭代,此时得到的概率积分参 数为全局最优. 基于上述求参思想,本文编制了基于D-In— 形;最后利用本文建立的求参模型,结合地质采矿 条件求取开采沉陷预计全部概率积分参数,并利用 求参拟合误差评价求参效果.模型求参的流程图如 2所示. SAR LOS向变形的概率积分参数模矢法反演软 以1为主影像 D—InSAR处理 提取1 ̄2时段 LOS向变形 提取2 ̄3时段 LOS向变形 工作面地质 采矿条件 以2为主影像 D—InSAR处理 LOS向移动 1 f获取l~ 时段f J利用本文求参 变形累加广-,、Los向累积变形广’1 模型求参 SAR影像缸l D InSA R 获取全部 概率积分参数 一 处理l’H撇\ 向变形 LOS 图2模型应用流程 Fig.2 Model application flow chart 伪SAR影像k 2求参模拟实验 2.1 D-InSAR LOS向采动移动变形模拟 D—InSAR变形时,首先在工作面上方地表划分南 北、东西方向间隔为10 m的方格网,模拟分辨率 为10 m(方位、距离向多视处理约1:3)D—InSAR 以淮南矿区煤系地层为背景,模拟工作面上覆 岩层岩性为中硬,煤层倾角为1O。,采高M一3.0 m,采深H一400 m,工作面开采尺寸沿倾向和走向 水平投影尺寸为D ×D。一300 1TI×8OO m,工作面 的变形图.然后根据概率积分法模型和参数真值预 计格网点的下沉值w 、北方向水平移动U 和东 方向水平移动U删,依据式(1)计算目标点沿LOS 方向累计移动量,并以此模拟D—InSAR的LOS向 沿走向布置,其倾向方位角为180。,顶板管理方法 为长壁垮落法.工作面采动程度满足D /H一0.7 实测地表累计移动量,模拟效果图如图3所示.为 了验证求参模型的可行性,在参数真值附近选取任 意初值(工程实践中,参数初值根据地质采矿条件 与概率积分参数的经验关系模型确定),然后基于 模拟的D—InSAR LOS向移动变形,利用本文建立 的融合LOS变形与模矢法的概率积分参数反演模 型进行概率积分参数反演,若重新求得的概率积分 参数和参数真值接近,则认为构建的概率积分参数 <1.2~1.4,D。/H一2>1.2~1.4,工作面倾向非 充分采动,走向超充分采动,总体为非充分采动.该 工作面开采沉陷预计的概率积分参数(参数真值) 如表1所示.假设sAR卫星回访几何成像时(以 TerraSAR为例),分辨率为3 m,雷达入射角为 一44.494。,卫星视线方向方位角a一350.596。. 模拟该工作面开采引起的地表LOS向累计 反演模型是可行的. 第5期 王 磊等:基于D-InSAR LOS向变形的开采沉陷预计参数反演方法研究 下沉/mm g \ 移动/mm 1400 定 椒 U 。rr .链 -.0 200 400 600 800 1000 1200 l400 1600 1800 2000 ...\l200 .极 .1000 控 800 东IIti方向/m (b)模拟南北方向水平移动 LOS向移动/ram l400 最 .0 .200 400 l200 .600 800 1000 1200 1400 1600 妇 .荨1000 链 800 ...—800 2000 2200 2400 2600 2800 3000 东西方向/m (c)模拟东西方向水平移动 东西刀向/m (d)模拟D—InSAR的LOS向移动 红色方框表示工作面,红色圆圈表示提取的求参像元 图3模拟求参数据分布情况 Fig.3 Distribution of simulation parameters 2.2概率积分参数反演模拟分析 L 线上求参像元模拟的LOS向移动变形值,以任 仿照常规煤矿开采地表移动观测站布设方案, 分别沿工作面走向、倾向设计求参模拟观测线L , L ,并以30 m为间隔选取60个离散求参点(像 元)(图3d).根据本文建立的求参模型,联合L , 意初值q一0.750,tan 一2.2O,6—0.25, 一85。,Sl 一15,S:一5,S。一2O,S 一25,进行模矢法求参,求 参结果见表l所示,求参拟合效果如图4所示. 表1 模拟求参实验结果 Table 1 Simulation and experimental results 从表1可以看出,q,tan ,b,0求参相对误差均 小于2.5 ,拐点偏移距S求参误差稍大,但不超过 2o (但对求参的整体效果影响甚微);从图4可以 看出,拟合误差在一ll~1l mm之间,拟合中误差为 ±5.8 mm.模拟求参实验表明,整体求参效果较好, 因而验证了本文建立的求参模型的可行性. 3工程应用研究 3.1 实验区及雷达数据简介 实验区位于兖州煤田南屯煤矿九采区,其上方 有邹济公路和矿用铁路穿过.2011年12月25日 至2012年3月l1日(雷达数据覆盖时间段)期间, 南屯矿正在回3上煤层的931O工作面西段和 l5 9308工作面东段,9310北侧和9308南侧工作面均 暑 10 回采结束形成老采空区且地表沉陷基本稳定.931O 和9308工作面煤层近水平(平均约为8。),工作面 。 5≈Ⅱ 1O 15 沿走向布置,采用综放采煤方法.931o工作面。采 厚4.90 m,工作面采深约588 m,工作面长197 m, 实验期间推进尺寸约220 m;9308工作面,采厚 4.76 m,工作面采深约630 m,工作面长140 m,实 验期间推进尺寸约415 m.9308,9310工作面及部 ..1 5 9 l3 l7 2l 25 29 33 37 4I 45 49 53 57 像元编号 图4拟合效果 Fig.4 Fitting effect 中国矿业大学学报 第46卷 分采空区(仅显示临近易采动“活化”的9312,9306 1 2 3 4 5 6 形梯度大易导致空间失相干特点。实验选取了2011 年12月25日至2012年3月11日期间时间基线为 老采空区)分布和工作面进尺进程如图5所示. 实验数据为Terrasar—X获取的单视复数雷达影 11 d和22 d的7景SAR影像进行实验.Terrasar—X 成像的主要几何参数和干涉方案见表2所示. 像,分辨率为3 1TI,影像大小为16 224像素×10 710 0 }2 2 2 2 }2 2 O 0 O 2 O 2 7 0 0 1 2 2 0 2 2 0 0 2 2 O 1 7 8 像素,实验数据覆盖了整个南屯煤矿.考虑到矿区变 l O 2;2 5 表2 实验区TerrraSAR—X影像几何参数和差分干涉方案 Table 2 Image geometric parameters of TerrraSAR—X and differential interference scheme in the experimental area O 1 5 6 干涉对 主影像 2 辅影像 o6 2口 口 0 时间基线/d O O 2 垂直基线/m 主影像入射角, (。) 方位向,(。) O O O 2 O 1 2 2 2 O O O ;一 l 一5 O O 2 一 7 O 0 O O 2 一 12一 3 1 一7 o0 1 3.2概率积分参数反演 ;l l 06 相等(可按相同考虑,误差可忽略),通过将各干涉对 i 1 l 90 对表2中的6对主从影像分别进行差分干涉测 量,获取了各干涉对沿I.OS向的地表移动变形量. 由于各干涉对主影像I.OS方位向相同、入射角近似 ● 一O 5 6 6 LOS向移动变形量累加,获得了相对与基准时间 (2011年12月25日)的6个时段内工作面开采引起 的LOS向地表累积移动变形如图5所示. 2 5 1 9 0 1 一 4 1 l 9 7 6 4 5 0 O 3 0 2 3 4 4 4 4 4 4 4 5 4 4 4 4 4 4 4 4 4 4 9 7 4 O 0 O 8 9 3 9 8 9 4 3 6 O l 1 3 5 3 5 3 5 3 5 3 5 3 5 O r0 9 6 O 一a 9 6 0 r0 9 6 O 0 r0 9 6 u 9 6 0 0 9 6 图5各时段LOS向累计地表移动变形量 Fig.5 Accumulated amount of surface movement deformat ion of LOS direction in each period 以9310工作面为例,该工作面2012年2月 88。,sl—S2一S 一S 一0,联合c,D线I ()s向移 动变形进行概率积分参数反演. 结果表明,在3月l1日时931O工作面开采沉 28日停采,截止2012年3月11日地表移动盆地 虽然未稳定,但可近似反演止3月】1日移动盆地 动态概率积分参数.由于931O工作面形状为矩形 且煤层为近水平,开采沉陷盆地具有对称性,可沿 走向西侧和倾向南侧两个方向分别取半条求参线 C,D(如图5f所示),并假设2011年12月25日之 前的工作面采动未对C,D线变形产生叠加影响. 陷的动态概率积分预计参数为:q一0.170,b= 0.12,tan口一2.12, =88。,Sl一一1 1,S2一l0,S 一 9O,S ===33,拟合误差约在一34~38 rnm之间, 拟合中误差m一±14.67 mlTl,拟合效果如图6所 示. 以任意初值q:0.200,tan 一2.00,b一0.20, 一 第5期 王 磊等:基于D-InSAR LOS向变形的开采沉陷预计参数反演方法研究 目 ● .∞如0卯∞∞∞ ∞如∞∞ 售 0044 器 S C1 C4 C7 Cl0C13C16C19C22C25C28 D2 D5 D8 Dl1Dl4D17D20D23D26D29 求参点名称 图6求参拟合效果 Fig.6 Fitting effect of simulation results 从求参结果来看,9310上侧拐点偏移距达到 S。一一90(拐点偏移距离出现大幅度外偏),这主要 因为9312工作面开采时在下山方向上存在较大的 悬臂梁或砌体梁存在,当9310工作面开采时引起 了悬臂梁或砌体梁失稳“活化”所致.其它概率积分 参数基本正常. 4 结 论 1)依据LOS向移动实质为下沉、南北和东西 方向水平移动沿L0s向投影的关系,基于模矢法 非线性规划理论,建立了基于单视线向D-InSAR LOS向变形的概率积分参数反演模型(可反演出 全部概率积分参数),编制了求参软件.并通过模拟 实验方式对求参模型的可行性进行了验证,模拟结 果表明,q,tan口,b,0求参相对误差均小于2.5 9/6, 拐点偏移距s求参误差稍大,但不超过20 (但对 整体求参效果影响甚微),求参拟合误差在一11~ ll mm之间,拟合中误差为m一±5.8 mm,整体求 参效果较好. 2)利用本文建立的参数反演模型,在一定的 假设条件下,对2012年3月11日9310工作面开 采沉陷的动态概率积分参数进行了求取,q一 0.170,6—0.12,tan 一2.12, 一88。,S1一一11,S2 —10,S3一一90,S 一33,拟合误差约在一34~38 mm之间,拟合中误差m一±14.67 mm.并认为S。 一一90主要为931O工作面开采导致临近9312采 空区边界的悬臂梁或砌体梁失稳“活化”所致. 3)本文构筑的求参模型是一个8参数非线性 复杂函数,当基于D-InSAR LOS向变形利用模矢 法求参时,防止把局部极值错判为全局最优值,可 根据地质采矿条件与概率积分参数经验关系,在参 数可能的范围内随机设计一个多值初值样本,并以 此为模矢基点进行搜素,若结果样本都近似收 敛于某一点,则可通过对结果样本加权(可以拟合 中误差进行定权)方式获得最优概率积分参数.采 用该策略可进一步增加求参的可靠性. 参考文献: [1] 国家煤炭工业局.建筑物、水体、铁路及主要井巷煤柱 留设与压煤开采规程[M].北京:煤炭工业出版社, 2000:31-47. E2] FAN H D,CHENG D,DENG K Z,et a1.Subsidence monitoring using D-InSAR and probability integral prediction modelling in deep mining areasEJ].Survey Review,2O15,47(345):438-445. [33 GRAY L.Using multiple RADARSAT InSAR pairs tO estimate a full three-dimensional solution for gla— cial ice movementl-J].Geophysical Research Letters, 2011,38(38):132—14O. -143 WRIGHT T J,PARSONS B E,LU Z.Toward map— ping surface deformation in three dimensions using InSAR[J].Geophysical Research Letters,2004,3 1 (1):169-178. -I5] Hu J,DING X L,LI z W,et a1.Kalman-filter based approach for muhisensor,multitrack and muhitempo— ral InSAR I-J 3.IEEE Trans Geosci Remote Sensing, 2013,51(7):4226-4239. -I6] FIALKO Y,sANDwELL D,SIMONS M,et a1. Three-dimensional deformation caused by the Barn.I— ran,earthquake and the origin of shallow slip deficit -lJ].Nature,2005,435:295—299. [77 HU J,LI Z W,ZHU J J,et a1.Inferring threedimen— sional surface displacement field by combining SAR interferometric phase and amplitude information of ascending anddescending orbits[J].Science China Earth Sciences,2O10,53(4):550-560. [8]JUNG H S,LU Z,WON J S,et a1.Mapping threedi— mensional surface deformation by combining multiple— aperture interferometry and conventional interferome— try:Application tO the June 2007 eruption of Kilauea Volcano,Hawaii[J].IEEE Geoscience and Remote Sensing Letters,2011,8(1):34—38. (下转第1180页) 118O 中国矿业大学学报 第46卷 报,2014,9(5):5I5-528. TAN Ying,ZHENG Shaoqiu.Recent advances in fireworks algorithm and application research.Journal of Xi'an Jiaotong University,2015,49(11):82—88. fireworks algorithmEJ].CAAI Transactions O12 In— telligent Systems,2014,9(5):515—528. [14] GAO H,DIAO M.Cultural firework algorithm and its application for digital filters design[J].Interna— tional Journal of Modelling,Identification and Con— trol,2O11。14(4):324—331. [13]李浩,柏鹏,张辉,等.反向烟花算法及其应用 研究口].西安交通大学学报,2015,49(I1):82—88. LI Hao,BAI Peng,ZHANG Hui,et a1.Backward (责任编辑邓群) (上接第1165页) 7DMUNDSS0N S,SIGⅣrUNDSS0N F。CARSTENS— [9] GIlite and in situ ground deformation data integrated by the SISTEM approach:The April 3,2010 earth— EN J M.Three—dimensiona1 surface motion maps esti mated from combined interferometric synthetic aper— quake along the Pernicana fault(Mt.Etna—Italy) ture radar and GPS data[J].Journal of Geophysical Researeh,2002,107(1O):13—14. case study[J].Earth and Planetary Science Letters, 2O11,312(3/4):327-336. Z W,DING X L,et aI.Derivation of 3-D [1o] HU J,LIcoseismic surface displacement fields for the 2011 [13] LI Z W,YANG Z F,ZHU J J,et a1.Retrieving three—・dimensional displacement fields of mining are—I Mw 9.O Tohoku—Oki earthquake from InSAR and as from a single InSAR pair[J].Journal of Geodesy, 2015,89(1):17-32. GPS measurements[J].Geophysical Journal Interna— tional,2O13,192(2):573-585. [14] 吴E15] 胡侃,葛家新,王铃丁,等.开采沉陷预计一体化方 俊.基于现代测量平差的三维形变估计理论与 [11] G0URMELEN N,KIM S W,et a1.Ice velocity de— termined using conventional and multiple—aperture 法[M].徐州:中国矿业大学出版社,1998:37—51. InsAR[J].Earth and Planetary Science Letters, 2O11,307(1/2):156-16O. 方法[D].长沙:中南大学地球科学与信息物理学 院,2O12:43—60. [12] GUGLIELMIN0 F。BIGNAMI C.Analysis of satel一 (责任编辑邓群)