CN105938193B - 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 - Google Patents
一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 Download PDFInfo
- Publication number
- CN105938193B CN105938193B CN201610554394.4A CN201610554394A CN105938193B CN 105938193 B CN105938193 B CN 105938193B CN 201610554394 A CN201610554394 A CN 201610554394A CN 105938193 B CN105938193 B CN 105938193B
- Authority
- CN
- China
- Prior art keywords
- mrow
- deformation
- rail
- insar
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法，首先利用升轨和降轨InSAR数据获取地表在两个不同斜距向上的相对地表形变测量值；然后融合升轨和降轨InSAR相对形变场估计地表在垂直向和东西向上的相对形变；进而利用东西向相对形变量级小于一定阈值的地面点计算升轨和降轨InSAR相对形变测量值的绝对偏移量；最后利用经绝对偏移量改正后的升轨和降轨InSAR形变值估计地表在垂直向和东西向上的绝对形变。突破了传统InSAR只能获取相对地表形变的技术瓶颈，无需布设地面控制点或者假设远场形变为零就可以得到地表在垂直向和东西向上的绝对形变结果。
Description
技术领域
本发明属于基于遥感影像的大地测量领域，尤其涉及一种无需地面辅助的升降轨InSAR 监测沉降区绝对地表形变的方法。
背景技术
合成孔径雷达干涉测量(InSAR)技术是近几十年发展起来的一种空间大地测量和遥感手段，具有监测范围大、测量精度高、空间近似连续、全天时全天候等优势，已经被广泛地用于地震、火山喷发、冰川漂移、滑坡、泥石流等地质灾害引起的地表形变监测。特别是随着多时相InSAR技术(如永久散射体和短基线集技术)的日益成熟，可以对同一轨道获取的时序SAR影像展开联合分析，从而更好的抑制了InSAR干涉图中的地形残差、大气延迟和失相关噪声等，进一步提高了InSAR形变监测的精度和可靠性，使得InSAR技术在监测由于地下流体和矿产资源开采、冻土冻胀和融沉等导致的缓慢长期的地面沉降时也展现出了极大的潜力，详见参考文献[1]和[2]。但是，由于InSAR技术采用的是相位测量方式，在形变解算过程中需要已知或假设其中一个监测点的形变量才能够解算其他监测点上的形变结果，因此 InSAR技术得到的地表形变场是相对的，和绝对形变场之间存在一个常数差。
如何将InSAR技术监测得到的地表相对形变转化为绝对形变，目前国际上的主流方法有以下两种：(1)地面控制点法，即在监测区域布设一个或多个地面观测站，利用水准、GPS 等手段获取这些地面离散点上的形变值，然后利用这些点上的实测形变值和InSAR形变测量值之差将InSAR的相对形变转化为绝对形变。但是对于没有布设地面控制点的地区，该方法无法适用。此外，InSAR所监测的地面点和实地测量的地面点也往往难以统一，影响了形变转化的精度。(2)远场形变假设法，即在远离地表形变区域的区域寻找稳定的地面点，并假设稳定点上的地表形变为零，从而得到InSAR形变场中其他地面点相对于该稳定点的地表形变结果。很明显，对于形变影响区域较大而导致无法找到稳定点的情况，该方法无法得到地表的绝对形变结果，详见参考文献[3]。
通过上述分析可以看出，对于没有地面控制点的地区，而地表形变的影响范围又超过了整个SAR影像的覆盖范围，那么现有的方法都无法得到地表的绝对形变结果。而在实际中，上述情况是非常常见的，例如在面对地下流体和矿产资源开采导致的大范围地面沉降时，以及冻土由于冻胀和融沉产生季节性垂直形变时，往往会出现难以找到可用的地面控制点或稳定点的情况，导致所监测到的形变结果不能真实反映地面沉降。特别是随着近年来高分SAR 数据的普及，InSAR形变场的分辨率越来越高，但覆盖范围却越来越小，使得InSAR技术的这个局限性愈加凸显。
发明内容
本发明的目的在于，克服现有InSAR技术难以监测沉降区绝对地表形变的局限性，提供一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法。
一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法，首先利用升轨和降轨InSAR数据获取沉降区地表在两个不同斜距向上的相对地表形变测量值；然后融合升轨和降轨InSAR相对形变场估计地表在垂直向和东西向上的相对形变；接着，利用东西向相对形变量级小于设定阈值的地面点计算升轨和降轨InSAR相对形变测量值的绝对偏移量；最后利用经绝对偏移量修正后的升轨和降轨InSAR形变值估计地表在垂直向和东西向上的绝对形变。
沉降区地表在两个不同斜距向上的相对地表形变测量值是通过收集所监测沉降区的升轨和降轨SAR数据，利用差分InSAR或多时相InSAR技术获取该区域地表在两个不同雷达视线方向上的相对形变测量值，将它们编码至统一的地理坐标系和格网下，并将它们的时间尺度和地面参考点统一；
所述沉降区即地壳运动以垂直向形变为主要特征的区域；
所述雷达视线方向即为斜距向；
所述升轨和降轨InSAR相对形变测量值的绝对偏移量按照以下公式计算：
其中，yi表示第i个挑选出来的地面观测点，i的取值为1-N，且地表观测点yi的东西向相对形变不大于设定的形变阈值即ε为事先设定的形变阈值；
代表第i个挑选出来的地面观测点在垂直向上的绝对形变；
Kas和Kdes分别代表升轨和降轨InSAR斜距向相对形变与绝对形变测量值之间的常数差，即升轨和降轨InSAR绝对偏移量；
和分别代表地面观测点yi上升轨和降轨InSAR斜距向相对形变测量值；
和分别为地表观测点yi的垂直向地表形变在升轨和降轨InSAR斜距向上的投影系数，θas和θdes分别表示地表观测点yi上的升轨和降轨雷达局部入射角。
利用SAR卫星的极轨特性，解算地表在垂直向和东西向上的相对形变。由于现有SAR 卫星都是在近南北的极轨轨道上飞行，因此无论是升轨还是降轨InSAR斜距向测量值都对南北向形变极其不敏感。
为了能够利用升降轨InSAR斜距向形变测量值估计垂直向和东西向形变，通常忽略南北向地表形变对InSAR斜距向形变测量值的贡献：
利用稀疏最小二乘算法解算上述方程，得到升降轨InSAR绝对偏移量，以及所挑选出来的N个地面观测点在垂直向上的绝对形变。
建立地面观测点xi上升降轨InSAR斜距向绝对形变测量值与垂直向和东西向绝对形变之间的函数关系，利用最小二乘平差解算出地面观测点在垂直向和东西向上的绝对形变和
其中，和分别代表地面观测点xi上升轨和降轨InSAR斜距向相对形变测量值；分别为地表观测点xi的垂直向、东西向地表形变在升轨InSAR斜距向上的投影系数，θas(xi)和αas(xi)则分别为地表观测点xi上的升轨雷达局部入射角和卫星飞行方向角；
为升轨和降轨InSAR斜距向相对形变测量值的误差；
分别为地表观测点xi的垂直向、东西向和南北向地表形变在降轨InSAR 斜距向上的投影系数；
θdes(xi)和αdes(xi)则分别为地表观测点xi上的升轨雷达局部入射角和卫星飞行方向角。
所述的形变阈值ε的取值小于或等于3mm，以保证忽略东西向形变的假设具有合理性。
所述沉降区表示地壳运动以垂直向形变为主要特征的区域，包括由于地下水开采、地下油气开采、地下采矿、地热开采、地下工程、冻土冻胀和融沉等引起地表形变的区域。
所述的升轨和降轨InSAR斜距向相对形变测量值由单一SAR卫星数据或联合不同SAR 卫星数据获取。
有益效果
本发明提供了一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法，1)利用InSAR技术获取沉降区地理编码后的升轨和降轨斜距向地表形变场的相对测量值，并统一它们的格网、时间尺度和地面参考点；2)根据SAR卫星的成像几何建立升降轨InSAR斜距向相对形变测量值与三维相对形变之间的函数关系；3)根据SAR卫星的极轨特性忽略南北向形变对InSAR斜距向测量值的贡献，估计垂直向和东西向上的相对形变；4)针对东西向相对形变小于一定阈值的地面观测点估计升轨和降轨InSAR斜距向相对形变测量值与绝对形变测量值之间的绝对偏移量；5)利用绝对偏移量将升降轨InSAR相对形变测量值转化为绝对形变测量值，并在此基础上估计垂直向和东西向上的绝对形变结果。该方法实现简单，无需在地面布设控制点或者寻找稳定区域，对地下开采、冻土活动等导致的绝对地表形变监测而言是一种高效率、低成本、高精度和大范围的方法，突破了InSAR技术本身只能监测相对地表形变的局限，进一步推动了InSAR大地测量技术的市场化和工程化进程。
附图说明
图1是升轨和降轨SAR卫星的成像几何图；
图2是本发明所述方法的流程图；
图3是模拟的沉降区地表形变，其中，(a)为垂直向形变；(b)为东西向形变；(c)为南北向形变；单位：cm/yr；
图4是含噪的InSAR斜距向相对形变测量值，其中，(a)为升轨测量值；(b)为降轨测量值；单位：cm/yr；
图5是垂直向和东西向上的相对形变结果；单位：cm/yr；
图6是东西向相对形变小于阈值的地面观测点分布；
图7是本发明得到的垂直向和东西向上的绝对形变结果；其中，(a)是垂直向绝对地表形变，(b)是东西向绝对地表形变，单位：cm/yr；
图8为应用本发明得到的绝对形变与模拟形变之间的差值；(a)是垂直向绝对和模拟的地表形变之间的差值，(b)为是东西向绝对和模拟的地表形变之间的差值；单位：cm/yr。
具体实施方式
下面将结合附图和实施例对本发明做进一步的说明。
为了便于理解本发明，首先提供本发明的理论基础：
众所周知，地质灾害的直接后果是地表发生形变。而对于InSAR技术而言，其只能监测地面观测点xi在斜距向上相对于某个参考点的相对形变，即
Drat(xi)＝Dreal(xi)+K (1)
其中，Drat(xi)为地面观测点xi上的InSAR斜距向相对形变测量值；Dreal(xi)为地面观测点xi上的斜距向绝对形变；K为InSAR斜距向相对形变测量值与斜距向绝对形变之间的常数差，即绝对偏移量。
然而在实际中，地表形变通常是发生在三维框架中的，即会同时发生垂直向、东西向和南北向形变。因此，斜距向绝对形变只是地表真实三维形变在斜距向上的投影，两者之间的关系可由下式表达：
其中，和分别为地表真实三维形变；Su(xi)、Se(xi)和Sn(xi)分别为地表观测点xi的垂直向、东西向和南北向地表形变在InSAR斜距向上的投影系数：
其中，θ(xi)和α(xi)则分别为地表观测点xi上的雷达局部入射角和卫星飞行方向角(以北方向为起始顺时针旋转)。升轨和降轨SAR卫星的成像几何图如图1所示。由于目前的SAR卫星均在太阳同步轨道(即极轨)上飞行，因此α(xi)通常在350度(升轨)和190度(降轨)左右。而作为侧视成像雷达，SAR的雷达局部入射角一般在20度到50度之间变化。因此，相应的投影系数变化范围分别大约为：0.7≤Su(xi)≤0.9、0.3≤Se(xi)≤0.6、0.08≤Sn(xi)≤ 0.10，从而导致InSAR的斜距向形变测量值对垂直向地表形变最为敏感，东西向地表形变次之，而南北向地表形变最差，详见参考文献[4]和[5]。
如图2所示，一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法，包括以下步骤：
(1)利用差分InSAR或多时相InSAR技术获取待监测沉降区地表分别在升轨和降轨雷达视线方向(即斜距向)上的相对形变测量值，并对它们进行地理编码，并统一它们的时间尺度和地面参考点；
(2)利用SAR影像头文件中包含的雷达局部入射角、卫星飞行方位角，按照公式(3)计算每个地面观测点上的投影系数Su(xi)、Se(xi)和Sn(xi)。然后根据公式(2)就可以构建升轨和降轨InSAR斜距向相对形变测量值与三维相对地表形变之间的函数模型：
其中为升轨和降轨InSAR斜距向相对形变测量值的误差；
(3)忽略南北向地表形变对InSAR斜距向形变测量值的贡献，则公式(4)可以写成如下形式：
根据最小二乘平差原理就可以解算出垂直向和东西向上的相对形变：
(4)挑选的地面观测点，其中ε为事先设定的形变阈值，一般不超过3mm/yr。对于地下开采、冻土活动等导致的地表形变而言，垂直形变较大的地方往往水平形变较小(如沉降漏斗中心)，而垂直形变较小的地方往往水平形变较大(如沉降漏斗边缘)。因此，可以假设挑选出来的N个地面观测点上的升轨和降轨InSAR斜距向相对形变测量值主要由垂直向形变所贡献，因此忽略东西向形变对其的贡献，从而建立以下函数模型：
其中，为挑选出来的N个地面观测点上升轨和降轨InSAR斜距向相对形变测量值，为相应的观测误差； Kas和Kdes分别代表升轨和降轨InSAR斜距向相对形变与绝对形变测量值之间的常数差，即升轨和降轨InSAR绝对偏移量。
上式中观测量的个数为2N个，而待求参数的个数为N+2个，因此观测量的个数大于待求参数的个数，可以通过最小二乘平差对其进行解算。但是由于设计矩阵为一个大型的稀疏矩阵，因此需要利用稀疏最小二乘对其进行求解，从而得到升轨和降轨InSAR绝对偏移量，以及所挑选出来的N个地面观测点在垂直向上的绝对形变。
(5)根据公式(1)将升降轨InSAR斜距向相对形变测量值转化为绝对形变测量值：
再根据公式(2)，同样忽略南北向形变对InSAR斜距向形变测量值的贡献，构建地面观测点xi上升降轨InSAR斜距向绝对形变测量值与垂直向和东西向绝对形变之间的函数模型：
其中为升轨和降轨InSAR斜距向绝对形变测量值的误差。同样利用最小二乘平差原理就可以解算出垂直向和东西向上的绝对形变：
在400×450的规则格网中模拟沉降区的三维绝对地表形变，格网尺寸为10m×10m，其中垂直向、东西向和南北向的绝对形变分别如图3(a)、3(b)和3(c)所示。然后利用公式(2) 模拟出升轨和降轨InSAR斜距向相对形变测量值，结果分别如图4(a)和4(b)所示，其中的三角形代表地面参考点。为了让模拟数据具有真实性，升轨数据的局部雷达入射角和卫星飞行方位角采用ALOS/PALSAR卫星影像头文件中提供的参数，并将均值为零、标准偏差为2mm 的高斯白噪声加入到升轨InSAR斜距向形变测量值；降轨数据的局部雷达入射角和卫星飞行方位角则采用ENVISAT/ASAR卫星影像头文件中提供的参数，并将均值为零、标准偏差为 1.5mm的高斯白噪声加入到降轨InSAR斜距向形变测量值。由于升轨和降轨模拟数据是直接模拟的同一地理坐标系下的年均形变速率结果，并且采用同一个参考点，因此在这次试验中无需时间尺度统一、地理编码、参考点统一等步骤。
通过本发明所提出的方法处理，就可以利用上述模拟的含噪升降轨InSAR斜距向相对形变测量值估计出垂直向和东西向上的绝对地表形变。图5(a)和5(b)分别是利用公式(6)解算出的地表在垂直向和东西向上的相对地表形变。可以看出，虽然在东西向上该结果与原始模拟的东西向绝对地表形变比较一致，但在垂直向上该结果的变化区间与原始模拟的垂直向绝对地表形变相差较大，这也说明传统的InSAR方法难以得到可靠、真实的地面沉降结果。图 6是本次试验挑选出来的东西向相对形变小于3mm的地面观测点，总共有48833个。图7(a) 和7(b)分别是利用公式(10)解算出的地表在垂直向和东西向上的绝对地表形变，总体而言和原始模拟的垂直向和东西向形变非常一致。图8(a)和8(b)给出的分别是该方法解算出来的垂直向和东西向形变场与模拟的垂直向和东西向形变场之间的差值。为了定量验证本发明的效果，实施例中分别计算垂直向和东西向形变的均方根误差，分别为2.1mm和2.6mm，与 InSAR斜距测量值中噪声的标准差相当，从而说明本发明是可行的，能够得到沉降区可靠的地表垂直向和东西向绝对形变结果。
参考文献：
[1]Berardino,P.,Fornaro,G.,Lanari,R.,Sansosti,E.,2002.A new algorithmfor surface deformation monitoring based on small baseline differential SARinterferograms.IEEE Transactions on Geoscience and Remote Sensing 40(11),2375-2383；
[2]Ferretti,A.,Prati,C.,Rocca,F.,2001.Permanent scatterers in SARinterferometry.IEEE Transactions on Geoscience and Remote Sensing 39(1),8-20；
[3]Liu,L.,Zhang,T.J.,Wahr,J.,2010.InSAR measurements of surfacedeformation over permafrost on the North Slope of Alaska.Journal ofGeophysical Research115,F03023；
[4]Fialko,Y.,Sandwell,D.,Simons,M.,Rosen,P.,2005.Three-dimensionaldeformation caused by the Bam,Iran,earthquake and the origin of shallow slipdeficit.Nature 435,295-299；
[5]Hu,J.,Li,Z.W.,Ding,X.L.,Zhu,J.J.,Zhang,L.,Sun,Q.,2014.ResolvingThree-Dimensional Surface Displacements from InSAR Measurements:AReview.Earth-Science Reviews 133, 1-17.
Claims (4)
1.一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法，其特征在于，首先利用升轨和降轨InSAR数据获取沉降区地表在两个不同斜距向上的相对地表形变测量值；然后融合升轨和降轨InSAR相对形变场估计地表在垂直向和东西向上的相对形变；接着，利用东西向相对形变量级小于设定阈值的地面点计算升轨和降轨InSAR相对形变测量值的绝对偏移量；最后利用经绝对偏移量修正后的升轨和降轨InSAR形变值估计地表在垂直向和东西向上的绝对形变；
所述升轨和降轨InSAR相对形变测量值的绝对偏移量按照以下公式计算：
其中，yi表示第i个挑选出来的地面观测点，i的取值为1-N，且地表观测点yi的东西向相对形变不大于设定的形变阈值即ε为事先设定的形变阈值；
代表第i个挑选出来的地面观测点在垂直向上的绝对形变；
Kas和Kdes分别代表升轨和降轨InSAR斜距向相对形变与绝对形变测量值之间的常数差，即升轨和降轨InSAR绝对偏移量；
和分别代表地面观测点yi上升轨和降轨InSAR斜距向相对形变测量值，为地面观测点上升轨和降轨InSAR斜距向相对形变测量值的观测误差；
和分别为地表观测点yi的垂直向地表形变在升轨和降轨InSAR斜距向上的投影系数，θas和θdes分别表示地表观测点yi上的升轨和降轨雷达局部入射角；
利用稀疏最小二乘算法解算上述方程，得到升降轨InSAR绝对偏移量，以及所挑选出来的N个地面观测点在垂直向上的绝对形变；
建立地面观测点xi上升降轨InSAR斜距向绝对形变测量值与垂直向和东西向绝对形变之间的函数关系，利用最小二乘平差解算出地面观测点在垂直向和东西向上的绝对形变和
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>a</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>a</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>a</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>a</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，和分别代表地面观测点xi上升轨和降轨InSAR斜距向相对形变测量值；分别为地表观测点xi的垂直向、东西向地表形变在升轨InSAR斜距向上的投影系数，θas(xi)和αas(xi)则分别为地表观测点xi上的升轨雷达局部入射角和卫星飞行方向角；
为升轨和降轨InSAR斜距向相对形变测量值的误差；
分别为地表观测点xi的垂直向、东西向和南北向地表形变在降轨InSAR斜距向上的投影系数；
θdes(xi)和αdes(xi)则分别为地表观测点xi上的升轨雷达局部入射角和卫星飞行方向角。
2.根据权利要求1所述的方法，其特征在于，所述的形变阈值ε的取值小于或等于3mm。
3.根据权利要求1所述的方法，其特征在于，所述沉降区表示地壳运动以垂直向形变为主要特征的区域，包括由于地下水开采、地下油气开采、地下采矿、地热开采、地下工程、冻土冻胀和融沉引起地表形变的区域。
4.根据权利要求1所述的方法，其特征在于，所述的升轨和降轨InSAR斜距向相对形变测量值由单一SAR卫星数据或联合不同SAR卫星数据获取。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610554394.4A CN105938193B (zh) | 2016-07-14 | 2016-07-14 | 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610554394.4A CN105938193B (zh) | 2016-07-14 | 2016-07-14 | 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105938193A CN105938193A (zh) | 2016-09-14 |
| CN105938193B true CN105938193B (zh) | 2018-04-06 |
Family
ID=56873466
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610554394.4A Active CN105938193B (zh) | 2016-07-14 | 2016-07-14 | 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105938193B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106556834B (zh) * | 2016-11-24 | 2018-11-16 | 首都师范大学 | 一种从两平行轨道sar数据集中精确提取地面垂直形变方法 |
| CN106767380B (zh) * | 2017-01-19 | 2019-04-19 | 中南大学 | 一种基于两景sar强度影像的矿区地表大量级三维形变估计方法 |
| CN107102332B (zh) * | 2017-05-11 | 2019-12-06 | 中南大学 | 基于方差分量估计与应力应变模型的InSAR三维地表形变监测方法 |
| CN110244298B (zh) * | 2019-07-26 | 2021-09-10 | 北京东方至远科技股份有限公司 | 一种InSAR数据升降轨联合滑坡分析方法 |
| CN110568440B (zh) * | 2019-09-10 | 2020-09-15 | 四川省地质工程勘察院集团有限公司 | 一种基于DS-InSAR技术监测复杂山区形变的方法 |
| CN111812645A (zh) * | 2020-06-10 | 2020-10-23 | 西南交通大学 | 冻土季节形变的卫星干涉测量方法 |
| CN112731397B (zh) * | 2020-12-28 | 2022-11-11 | 上海卫星工程研究所 | 一种无地面控制的InSAR绝对相位确定方法及系统 |
| BR102021005841A2 (pt) * | 2021-03-25 | 2022-09-27 | Radaz Industria E Comercio De Produtos Eletronicos Ltda | Método de localização de alterações no subsolo |
| CN113091598B (zh) * | 2021-04-06 | 2022-02-08 | 中国矿业大学 | 一种InSAR划定采空区建筑场地稳定性等级范围的方法 |
| CN113487105B (zh) * | 2021-07-27 | 2023-04-07 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
| CN115201822B (zh) * | 2022-07-07 | 2023-03-14 | 长沙理工大学 | 一种钻井水溶岩盐矿区采卤量估计方法 |
| CN116522517B (zh) * | 2022-12-05 | 2023-11-21 | 首都师范大学 | 一种量化地面沉降不均匀程度及沉降漏斗稳定性的方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706577A (zh) * | 2009-12-01 | 2010-05-12 | 中南大学 | InSAR监测高速公路路面沉降方法 |
| CN102927934A (zh) * | 2012-11-07 | 2013-02-13 | 中南大学 | 一种利用单个InSAR干涉对获取矿区地表三维形变场的方法 |
| CN103091675A (zh) * | 2013-01-11 | 2013-05-08 | 中南大学 | 一种基于InSAR技术的矿区开采监测方法 |
| CN104062660A (zh) * | 2014-07-14 | 2014-09-24 | 中南大学 | 一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法 |
| CN105158760A (zh) * | 2015-08-10 | 2015-12-16 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
-
2016
- 2016-07-14 CN CN201610554394.4A patent/CN105938193B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706577A (zh) * | 2009-12-01 | 2010-05-12 | 中南大学 | InSAR监测高速公路路面沉降方法 |
| CN102927934A (zh) * | 2012-11-07 | 2013-02-13 | 中南大学 | 一种利用单个InSAR干涉对获取矿区地表三维形变场的方法 |
| CN103091675A (zh) * | 2013-01-11 | 2013-05-08 | 中南大学 | 一种基于InSAR技术的矿区开采监测方法 |
| CN104062660A (zh) * | 2014-07-14 | 2014-09-24 | 中南大学 | 一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法 |
| CN105158760A (zh) * | 2015-08-10 | 2015-12-16 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 利用InSAR资料反演缅甸Mw6.8地震断层滑动分布;周辉 等;《地球物理学报》;20130930;第56卷(第9期);第3011-3021页 * |
| 升降轨道ASAR雷达干涉揭示的巴姆地震（Mw6.5）3D同震形变场;孙建宝 等;《遥感学报》;20060731;第10卷(第4期);第489-496页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105938193A (zh) | 2016-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105938193B (zh) | 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 | |
| Zhu et al. | Land subsidence due to groundwater withdrawal in the northern Beijing plain, China | |
| Dehghani et al. | Hybrid conventional and persistent scatterer SAR interferometry for land subsidence monitoring in the Tehran Basin, Iran | |
| Tong et al. | High‐resolution interseismic velocity data along the San Andreas Fault from GPS and InSAR | |
| Catalão et al. | Merging GPS and atmospherically corrected InSAR data to map 3-D terrain displacement velocity | |
| Tung et al. | Assessments of serious anthropogenic land subsidence in Yunlin County of central Taiwan from 1996 to 1999 by Persistent Scatterers InSAR | |
| CN108983232B (zh) | 一种基于邻轨数据的InSAR二维地表形变监测方法 | |
| CN105158760B (zh) | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 | |
| CN106556834B (zh) | 一种从两平行轨道sar数据集中精确提取地面垂直形变方法 | |
| Evain et al. | Structure of the Lesser Antilles subduction forearc and backstop from 3D seismic refraction tomography | |
| Farolfi et al. | A procedure to use GNSS data to calibrate satellite PSI data for the study of subsidence: an example from the north-western Adriatic coast (Italy) | |
| Kamiński et al. | Electrical resistivity tomography and digital aerial photogrammetry in the research of the “Bachledzki Hill” active landslide–in Podhale (Poland) | |
| Falorni et al. | Advanced InSAR techniques for geothermal exploration and production | |
| Oja et al. | Gravity anomaly field over Estonia | |
| Berrocoso et al. | Determination of geomorphological and volumetric variations in the 1970 land volcanic craters area (Deception Island, Antarctica) from 1968 using historical and current maps, remote sensing and GNSS | |
| Xing et al. | Time series ground subsidence inversion in mining area based on CRInSAR and PSInSAR integration | |
| Su et al. | Monitoring interseismic activity on the Ilan Plain (NE Taiwan) using Small Baseline PS-InSAR, GPS and leveling measurements: partitioning from arc-continent collision and backarc extension | |
| Xiao et al. | GPS and InSAR time series analysis: deformation monitoring application in a hydraulic engineering resettlement zone, Southwest China | |
| Artese et al. | Monitoring of land subsidence in Ravenna municipality using integrated SAR-GPS techniques: Description and first results | |
| Morelli et al. | Iso-Kinematic Maps from statistical analysis of PS-InSAR data of Piemonte, NW Italy: Comparison with geological kinematic trends | |
| Ahmad et al. | Characterization of Upper Indus Basin based on DEM hypsometric analysis | |
| Walstra et al. | Extracting landslide movements from historical aerial photographs | |
| Densmore et al. | Evaluation of land subsidence and ground failures at Bicycle Basin, Fort Irwin National Training Center, California, 1992–2017 | |
| Liu et al. | Kilometer-resolution three-dimensional crustal deformation of Tibetan Plateau from InSAR and GNSS | |
| Fabris et al. | Monitoring Land Subsidence Using Remote Sensing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |