引言

数字高程模型(digital elevation modal)的概念于20世纪60年代被提出，它是伴随着计算机科学、传感器和图像处理等技术发展起来的，例如我国研究人员陈述彭、王之卓、廖克、刘岳等于20世纪70年代末利用计算机辅助制图建立了DEM网格(汤国安，2014). 20世纪七八十年代，由于摄影测量仍处于解析摄影测量阶段，因此航空摄影测量成为当时获取DEM最成熟的手段，也是在此阶段，各国提出了相应的卫星对地观测计划并发射卫星，例如美国的锁眼侦查卫星和俄罗斯的“宇宙”号系列卫星. 80年代以后，法国和印度分别发射了SPOT系列卫星和IRS系列卫星，开启了商用卫星的发展.日本1999年发射的IKONOS被认为是第一颗高分辨率的商业卫星，其全色分辨率优于1 m，自此之后，高分遥感立体测量技术逐渐发展成熟起来(孙刚，2007).在卫星测量发展成熟阶段，一种新的地形测量技术，即激光雷达测距(light detection and ranging，简写为LiDAR)技术悄然兴起，成为目前精度最高的高空间分辨率数字高程模型(very high resolution digital elevation model，简写为VHR DEM)提取技术.

DEM是地表形态的数字化表达，蕴含着丰富的地学应用分析所必需的地形地貌信息(汤国安等，2010). 2000年，美国“奋进”号航天飞机完成了全球DEM的建立，并免费开放了空间分辨率为90 m的数据，极大地促进了宏观地形地貌和地震活动等研究.在构造地貌学领域，大区域的宏观地貌定量化研究已普遍开展(张会平等, 2004, 2006；刘静等，2006；洪顺英等，2007；胡艺，2008；高明星，刘少峰，2008；赵洪壮等，2009；程三友等，2010)；在地震滑坡研究中，许多研究人员利用DEM提取高程、坡度和坡向等重要信息(陈晓利等, 2010, 2013；许冲等，2010；王进等，2012；Niethammer et al，2012；王志恒等，2014；庾露等，2014；刘甲美等，2015；宿星等，2017).近10年来，高分遥感立体测量技术、LiDAR技术以及新兴的数字摄影测量技术已成为获取DEM的主要手段，所提取DEM的空间分辨率达到米级，甚至亚米级，从而使DEM达到高空间分辨率水平，即VHR DEM.因此刘静等(2013)提出了基于VHR DEM的微地形地貌的研究方向.由于地震往往发生在活断层上，并形成一系列特殊的地震地貌(杨景春，李有利，2011)，因此，VHR DEM被快速地应用于地震活断层研究领域.本文参考国内外最新的研究成果，简要介绍了VHR DEM的获取技术，总结了VHR DEM在地震活断层研究中的应用现状，并重点分析了其在活断层识别与几何结构分析、同震位移与累积位移获取及古地震研究等方面的发展.

断层的长期和突发活动是对地表形态的再塑造，将DEM作为基础研究数据是必然的，同时，断层活动造成的地表变形使原始DEM失去了时效性，而有效的DEM数据在抗震救灾和灾后重建等方面至关重要.以上特性决定了DEM在活断层研究中的重要地位，而VHR DEM是地学研究能否达到新的理论水平的关键，DEM的时效性特点也要求DEM获取需从定期转变为动态，这同时也对DEM的构建和分析理论提出了挑战.

1.   VHR DEM获取技术

DEM的概念自Miller于1958年提出(Miller, Laflamme, 1958)后发展至今，已成为地学研究赖以进行的核心数据. DEM根据空间分辨率基本可以分为低(百米)、中(几十米)、高(米、亚米)3个水平.高空间分辨率DEM，即本文所述的VHR DEM，具有高精度、米级至亚米级空间分辨率的特点. VHR DEM提取技术众多，各有其特点和应用范围.本文就现阶段常用的VHR DEM提取技术作简要介绍和比较.

1.1   高分遥感立体测量

立体测量是遥感卫星技术中一个重要的发展和应用方向.从广义上讲，具备立体测图或者高程测量能力的卫星均称为测绘卫星; 从狭义上说，目前一般把从不同视角获取同一地区影像的光学遥感卫星称为测绘卫星(唐新明等，2012).现阶段，具有立体测图能力的高分光学遥感卫星越来越多，表 1列出了部分高分测绘卫星的发射情况.

表  1  高分测绘卫星参数

Table  1.  Parameters of high resolution optical remote sensing mapping satellites

遥感卫星名称	发射 年份	地面分辨率/m		重访周 期/d	幅宽 /km	平面精度/m			高程精度/m
		全色	多光谱			无控制	有控制		无控制	有控制
美国 IKONOS-2	1999	1	4	3	11	12	2		10	3
法国 SPOT-5	2002	2.5/5/10	10	3	60	-	10		10	5
美国 Orbview-3	2003	1	4	3	8	11	-		16	-
美国 Worldview-1	2007	0.41	8	1.7	17.6	7.6	2		-	-
美国 GeoEye-1	2008	0.5/0.41	1.64	3	15.2	4	2		6	3
美国 Worldview-2	2009	0.46	1.84	1	16.4	4	2		3	2
法国 Pleiades	2011	0.5	2	1	20	-	-		1	0.5
法国 SPOT-6	2012	1.5	6	3	60	-	-		-	-
中国资源三号	2012	3.5/2.1/3.7	5.8	5	50	6	1.6		8	1.6
美国 Worldview-3	2014	0.31	1.24	4.5	13.1	2.16	-		1.61	0.62
注：表中各卫星获取DEM的精度仅供参考，因为不同研究人员在进行精度评价时其影像覆盖区域的地形不同，数据处理过程也有差异.表中参数引自Deilami和Hashim (2011)，潘红播等(2013)，朱红等(2014)，兰穹穹等(2015)，Jacobsen和Topan (2015)以及Hu等(2016).

下载: 导出CSV 

| 显示表格

从表 1中可以看出，几乎每年都有高分测绘卫星升空.民用测绘卫星的空间分辨率已提高至0.31 m，测图精度可以满足1:5000地形图的要求.卫星重访周期均低于5天，Worldview-2和Pleiades的重访周期甚至可达1天.另外，立体测量过程中加入地面控制点可以有效提高测量精度，而近几年发射的高分卫星的无控测量精度已与最初发射的高分卫星有控测量精度相当，如Worldview-3的平面无控测量精度达到了IKONOS-2的平面有控测量精度.无地面控制点条件下的几何定位精度是评价测绘卫星性能的重要指标(唐新明等，2012).所以现阶段高分光学遥感立体测量的发展趋势是：高空间分辨率、短重访周期卫星的研发，精密传感器的设计以及卫星轨道和姿态的精密测定，提高无地面控制点测量的精度.

雷达卫星是另一种重要的对地观测卫星.利用合成孔径雷达干涉(interferometric synthetic aperture radar，简写为InSAR)测量技术提取DEM的构想是Rogers和Ingalls (1969)提出的. Graham (1974)首次基于机载雷达制作了InSAR地形图，Goldstein等(1988)又将这项技术应用到星载雷达，从而大大拓宽了InSAR的应用范围. 2000年，美国“奋进”号航天飞机执行雷达地形测绘使命(Shuttle Radar Topography Mission，简写为SRTM)并获取了全球数字高程模型，随后便面向全球开放了空间分辨率为90 m的SRTM DEM，该数据被广泛地应用于宏观地形地貌的研究；2009年，美国航天局与日本经济产业省基于Terra卫星拍摄的130万个立体图像合作推出了先进星载热发射和反射辐射仪全球数字高程模型(advanced spaceborne thermal emission and reflection radiometer global digital elevation model，简写为ASTER GDEM)，其空间分辨率达30 m；2010年，德国航空航天中心(Deutsches Zentrum für Luft-und Raumfahrt, 简写为DLR)发射了TanDEM-X雷达卫星，与之前发射的TeraSAR-X卫星组成了双星SAR系统，构建推出了空间分辨率为12 m的DLR DEM. InSAR DEM的高程精度，在平坦地区约为2 m，在地形起伏较大区域约为5 m(许才军等，2015)，InSAR DEM提取精度除了受地形因素影响外，还受大气、失相干等因素的影响.从目前发展来看，InSAR提取VHR DEM的发展方向是：高分SAR卫星的研发和多源数据融合算法的研究(Bhardwaj et al, 2013; Jiang et al，2014；杜亚男等，2015).

1.2   机载LiDAR技术

LiDAR技术出现于1970年(Kaula et al, 1974). 20世纪八九十年代，大多数研究主要集中于激光测高传感器的研制，而全球定位系统(global positioning system，简写为GPS)和惯性导航系统(inertial navigation system，简写为INS)的出现彻底解决了机载LiDAR的定位和姿态测量问题(刘经南，张小红，2003). 1988年德国斯图加特大学利用GPS接收机、INS及激光扫描仪组成地形断面测量系统，被视为现代机载LiDAR系统的原型(Ackermann, 1999).中国科学院和武汉大学最早开始自主研制我国机载LiDAR系统(李清泉等，2000；李树楷等，2000).目前，机载LiDAR系统的硬件和软件由各大生产厂商提供，如Leica，Trimble公司，其扫描频率高达250 kHz，精度为厘米级，点云密度达到20—50个/m²，可直接生成高精度的DEM数据(陈涛，2014).未来机载LiDAR在DEM提取方面的发展趋势是：多源数据的融合处理，例如融合影像数据、多光谱数据、地面已知GIS数据等；点云滤波和分类算法的研究，需注意的是，该方法目前在高山和峡谷等地区容易出现分类错误(张小红，刘经南，2004；周晓明，2011).

与传统遥感手段相比，机载LiDAR的主要技术优势为(Baltsavias，1999；马洪超，2011)：激光雷达能部分穿透植被，获得真实的地表高程信息，利用高密度点云生成高精度DEM；受天气影响比光学遥感卫星测量要小；可在地震、火山等自然灾害情况下快速开展应急.面临的主要问题是开展LiDAR扫描工作的成本太高，在高山、峡谷等坡度较陡地区，实现点云滤波处理困难较大.

1.3   SfM技术

运动重建(structure from motion，简写为SfM)技术不是一种严格意义上的数字摄影测量技术，在地形地貌研究上的运用是近几年才发展起来的.摄影测量技术的发展经历了模拟摄影测量、解析摄影测量和数字摄影测量.模拟摄影测量已经基本淘汰(张祖勋，2007).解析摄影测量有专门的解析测图仪，但现阶段已很少使用解析测图仪进行DEM提取，而是使用在普通计算机上运行的专业软件，如EREAS，完成DEM的解析计算.数字摄影测量是目前研究的热点，主要有两个研究方向：一是数据采集端，二是数据后处理系统(刘先林等，2007)，基于SfM的新型摄影测量技术就是在后者上发展起来的. SfM开发的初衷是为了自动生成来自无序图像集合的3D模型(Brown, Lowe，2005)，与摄影测量计算像点点位的目的相同.该技术在运用于地表模型恢复之前，已经开发了众多的成熟软件，例如Agisoft Photoscan，Bunlder Photogrammetry Package，PhotoModeler等.近几年，一些研究人员利用航空摄影测量平台拍摄了大量航片，使用SfM处理软件得到地表的点云数据(James，Robson，2012；Westoby et al，2012).由于得到的点云数据是在任意坐标系下，因此需要结合地面控制点将其转化至地理坐标系下. Westoby等(2012)在陡峭的沿海山坡、冰川碛和基岩山脊处拍摄了大量照片，并利用SfM技术构建了地面模型，该模型与LiDAR DEM的残差小于0.5 m. James和Robson(2012)也用近距离(20 m)拍摄的地面照片，制作了7个时间序列的SfM模型，与LiDAR DEM的残差仅几厘米.以上实验研究表明，使用SfM技术构建的VHR DEM具有很高的精度; 但SfM技术与LiDAR技术一样，均需对地面点云数据作后处理，不同的是SfM生成的点云数据除了具有位置属性外，还具有色彩属性(RGB值)(Fonstad et al, 2013; Johnson et al, 2014).

以上技术在提取VHR DEM时各有其优势.从研究的空间尺度上来看：LiDAR和SfM均可实施地面测量，通过将其搭载于航空器(飞机、氢气球)上，对地面物体进行三维提取，建立某一区域的地面模型；而高分遥感卫星的覆盖范围更广，系统更加稳定可靠，可以完成大范围的地形测绘任务.从对地表观测的精度上来看：LiDAR可以获得更真实的地形，这是其它技术无法比拟的，其生成的VHR DEM精度最高，常用其评价其它技术提取的DEM的精度；SfM技术是一种广义上的摄影测量技术，比大多数测绘卫星的立体测量精度要高.从资金成本上来看：LiDAR的成本最高，SfM比LiDAR低得多，在能够满足精度要求的情况下，建议优先选择高分遥感立体测量技术.

VHR DEM主要是对微地形地貌的表达，因此对DEM的精度要求很高.所以，LiDAR技术在地震活断层的研究中发展得最快、应用得最广，由此产生的关于活断层的新研究方法开拓了其它技术的应用方向. SfM技术是近几年才被应用于地学研究，因此研究成果少，但发展潜力大.随着近几年民用高分遥感卫星的发射，开始有研究者使用高分立体像对提取高精度的VHR DEM，并将其应用于活断层研究.在开展活断层研究时，应以研究内容为主，根据研究区地形地貌特点选择合适的测量技术.

2.   VHR DEM在活断层研究中的应用

活断层的研究手段主要包括传统的野外调查、航片卫片解译、GPS、InSAR空间对地观测以及数值模拟等.技术手段的变革深刻地影响着不同学科对活断层的研究方式，VHR DEM极高的空间分辨率使我们从形态特征、运动学特征等方面对断层有了更深刻的认识.下文将主要基于相关研究者的最新成果和前人的总结(马洪超，2011；Meigs, 2013；任治坤等，2014)，从活断层识别及几何结构分析、同震位移与累积位移及古地震等方面论述VHR DEM在活断层研究中的独到优势及其所带来的变革.

2.1   活断层识别及几何结构分析

断层位置的准确识别和提取是活断层研究的基础，亦是评价和研究断层体系的关键.最初将LiDAR技术运用到断层研究中源于一次偶然的发现，1997年美国地质学家在利用LiDAR扫描地形的过程中发现了一个高达5 m的断裂陡坎，陡坎周围森林密布，此前的野外调查和航空解译均未发现此陡坎.相关研究者敏锐地洞察到了这一技术优势(Haugerud et al, 2003)，随后利用LiDAR技术的大型断裂带精细扫描项目陆续开展，例如美国B4项目对圣安德烈斯断裂南段和圣哈辛托断裂的扫描(Bevis et al, 2005)，该项目的开展对活断层研究产生了深远的影响(Arrowsmith, Zielke, 2009; Zielke et al, 2010, 2012；Zielke，Arrowsmith，2012).我国自汶川地震后开始应用该技术进行相关研究(李峰等，2008；马洪超等，2008；袁小祥等，2012)，并在2011年对海原断裂带开展大范围机载LiDAR扫描，获得丰硕成果(刘静等，2013；Chen et al，2014；李占飞等，2016).

在植被覆盖茂密的山区和城市地区，断层填图一直是难点.另外，一些破坏性大的地震常常发生在研究较少的分支断层或未识别断层上，而我们对这些断层却知之甚少，一个重要的原因就是断层信息难于获取.多位学者先后在研究过程中分析比较了VHR DEM在活断层识别中的优势. Lin等(2013)通过对比分析分辨率为10 m，2 m，0.5 m，0.25 m的DEM数据，认为0.5 m的DEM更适于识别微小构造断裂，结合赤色图(red relief image map，简写为RRIM)可视化方法，可识别出许多茂密植被覆盖下的微小陡坎(Chiba et al 2007, 2008). Arrowsmith和Zielke(2009)采集LiDAR点云数据，利用反距离加权法对点云数据进行插值运算，制作分辨率为0.25 m和0.5 m的DEM，仅依靠此VHR DEM进行了圣安德烈斯断裂局部活断层的构造地貌填图，并与航片和野外观测相结合的经典断层提取方法相比较，发现两者精度基本一致. 魏占玉等(2014)也进行了相似研究，其结果表明LiDAR DEM在宏观上足以获得与航片解译和野外调查一致的断裂分布特征，但在微观上具有更高的精细度. Johnson等(2014)首次利用SfM低空摄影测量绘制了断裂带地形，拍摄了几百张研究区照片，生成每平米500点以上的点云数据，SfM DEM分辨率达厘米级.与LiDAR DEM相比，SfM DEM分辨率更高，且能更清晰地展示微小陡坎等断错地貌.以上相关实验研究表明，VHR DEM数据在断层识别中充分体现了其高空间分辨率的优势，可以准确识别、测量几十厘米的微小陡坎、断错冲沟和冲积扇等地貌(Nissen et al, 2014；Scharer et al, 2014)，便于补充完善早期断裂地质图，另外可在研究人员难以到达的地区提取断层，极大地减少野外工作.

研究活断层与地震的关系，除了需要准确识别断层位置外，更重要的是获取断层的长度、走向、倾向、倾角、组合形式等几何要素，例如较大断层之间的阶区分布和断层本身弯曲变化往往抑制地表破裂的传播(Elliott，2014).活断层的几何结构关系到地震灾害的破坏程度，也是地震预测研究的重要内容.目前常用的研究方法有构造地质法、地震剖面及数值模拟推断法等，VHR DEM在断层几何结构分析方面的应用将是重要的发展方向之一，但现阶段该方面的研究成果比较少. Zhou等(2016)利用高分立体像对提取的VHR DEM计算了2013年9月24日巴基斯坦M_W7.7地震发震断层的倾角.他们假设研究区断层面为平面且破裂面穿透地表，建立如式(1) 所示的断层面模型，即

式中，A，B，C，D为平面参数，x, y, z为破裂面与地面交线上的点，则断层倾角的计算公式为

计算结果显示Hoshab断层西部端点处的倾向发生了变化，进一步通过遥感解译发现断层在此处分成东、西两部分，存在抑制地表破裂发育的阶区，而在此前的研究中均未发现该阶区.需要注意的是，Zhou等(2016)的模型计算的是地表处的断层面倾角，即破裂面穿透地面的有限部分.

2.2   同震位移与累积位移

理论上数字高程模型可以直接捕捉地表形变信息，但中等分辨率的DEM精度较低，即使在较大的地震中也难以发挥作用.随着VHR DEM数据的使用，获取微小地表形变成为可能.此外，基于VHR DEM还可以获得垂直方向上的位移.利用VHR DEM获取同震位移的方法大致有两种：断错地貌标志法和DEM差分法，其中传统的断错地貌标志法虽然只需震后DEM数据，但需研究人员主观识别断错山脊、冲沟、阶地等(Arrowsmith, Zielke, 2009; Chen et al，2014)，而DEM差分需要地震前后的DEM数据，数据处理过程自动化程度高且结果可靠，但该方法对DEM数据的空间分辨率和精度要求较高.

Hudnut等(2002)首次运用LiDAR技术获取了1999年10月16日加利福尼亚M_W7.1地震的同震位移.对发震断层进行全覆盖激光扫描，提取断层两侧地形剖面并进行比较，获得同震水平位移约为4.2 m，垂直位移约为0.9 m，与前人野外测量和InSAR观测结果相近，从而证明了该方法的可行性和可靠性. Zielke等(2010, 2012)及Zielke和Arrowsmith (2010)对1857年1月9日圣安德烈斯断层上发生的M_W7.9地震基于VHR DEM开展同震位移研究，在测量了大量的断错冲沟、洪积扇、阶地等后重新评定此次地震的同震水平位移约为6 m，而不是之前的野外测量结果(约9 m)，由此对圣安德烈斯断层破裂模型提出了质疑. Zielke和Arrowsmith (2012)提出的位错测量原理及开发的LaDiCaoz软件，得到了相关研究者的认同，如Chen等(2014)基于该方法对海原断裂1920年12月16日M_W8.3地震哨马饮区段进行了研究. 王朋涛等(2016)将数码相机搭载于无人机上，对海原断裂上骟马沟及其西侧阶地断错地貌进行了图像采集，将采集到的106幅照片利用SfM技术处理生成了0.1 m空间分辨率的SfM DEM，更加真实地反映了微地形地貌. 魏占玉等(2015)在独山子逆断裂-背斜带的机载LiDAR测量过程中，借助飞行平台用数码相机拍摄了大量照片，分别生成了SfM和LiDAR点云，SfM的点云密度为25.5个/m²，生成分辨率为0.2 m和0.5 m的SfM DEM均可以清楚地表示小冲沟，而LiDAR的点云密度仅为6.9个/m²，生成的0.5 m分辨率的LiDAR DEM表示不出这些小冲沟，并且LiDAR点云与SfM点云相比，垂直偏差小于0.1 m的比例为58.3%，垂直偏差小于0.2 m的比例为88.3%，这说明SfM DEM的精度与LiDAR DEM相当.

Oskin等(2012)对2010年4月4日EI Mayor-Cucapah M_W7.2地震前后的LiDAR DEM进行差分，获取了该地震的近场形变.虽然该地震发生在多级断层组成的断裂带上，构造背景十分复杂，但VHR DEM差分图清晰、直观地显示了120 km长的地表破裂带，与三维弹性模型的模拟结果吻合. Zhou等(2015a)首次评定了Pleiades-1A/B卫星提取的VHR DEM在获取EI Mayor-Cucapah地震同震垂直位移中的稳定性，并从点云密度、地形坡度等方面评定了VHR DEM精度，认为在坡度较小、植被稀疏、特征地物丰富的地区所提取的VHR DEM具有较高精度，通过与震前LiDAR DEM差分获得同震水平位移，且其结果与Oskin等(2012)所得结果相近.

此外，构造地貌法常常根据断错的地貌标志来研究活断层的长期运动规律，例如通过测量河流阶地的累积位移研究断层的第四纪活动速率. Cowgill(2007)及Cowgill等(2009)在此方面作了大量的研究，但大多基于航片、卫片及野外观测.在2010年1月12日海地M_W7.0地震中，Cowgill等(2012)基于震后快速获取的LiDAR数据，提出了可交互的沉浸式虚拟现实技术在快速开展地震科学研究中的应用，LiDAR DEM可通过3D形式再现断错地貌标志，达到与野外观测相似的位错测量结果，尤其是垂直位错. Zhou等(2015b)在无野外观测的情况下，利用2013年9月24日巴基斯坦M_W7.7地震后Pleiades立体像对提取了VHR DEM，并测量了大量断错的累积位移和同震位移，之后分别计算位移比值R，即

式中，L为水平位移，V为垂直位移.结果表明同震位移与累积位移其R值相似，由此判断该断层的运动方式并未发生变化.由此可见，VHR DEM给基于构造地貌法研究断层长期活动规律的传统方法带来了新的启示.

2.3   古地震

古地震研究对认识断层活动的时空规律及评价地震危险性均具有重要意义.古地震的研究方法非常广泛，包括地貌标志法、地层分析法、树木年轮变异法等.确定断层的活动幅度、年代、次数、时间间隔和活动速率等是古地震研究的主要内容.例如利用断错地貌标志测定累积位移，再结合相应的测年数据就可计算出断层的长期活动速率(Cowgill，2007；张培震等，2008；Cowgill et al, 2009).

反映古地震的地质地貌现象常因被埋藏在地下而未受到侵蚀破坏，所以开挖探槽是进行古地震研究的常用手段(杨景春，李有利，2011). Haddad等(2012)利用地基LiDAR扫描古地震探槽，发现LiDAR技术比传统研究手段更具优势：传统人工手绘作业需要大量的作业时间，且易遗漏细节信息，对探槽壁拍照后需要对几百张照片进行室内拼接，不同拍摄角度和探槽壁粗糙度会加大拼接变形，后期校正处理工作量大；而LiDAR扫描的点云具有空间坐标，能够完整地保存和记录探槽信息，避免影像拼接造成的影响，更可避免塌方等危险对工作人员的伤害(郑文俊等，2015).此外，可将采集到的古地震探槽信息统一到同一地理框架下进行集成化管理.在无法开挖探槽的基岩地区，依据保存较好的出露断层面形貌特征进行古地震识别已得到相关研究人员的关注(Giaccio et al, 2003；Bubeck et al, 2015). 何宏林等(2015)利用地基LiDAR扫描了霍山山前断裂的基岩断层面，生成了分辨率为2 mm的断面DEM，精确地反映了断层面的粗糙程度，进而判断断层面出露时间，之后又基于分形理论计算断层面2D分维值并分析其在垂直方向上的变化，最终识别出3次古地震事件.另外，断层的垂直运动会造成河流的溯源侵蚀并形成河流裂点，裂点的分布可以揭示古地震的活动特性. 毕丽思等(2011)同样以霍山山前断裂为试验区，从立体像对中提取DEM数据，获取横穿断裂的23条冲沟的地形剖面，最终揭示了霍山山前断裂的3次古地震事件，并确定强震重复周期为1500—2600年.可见，VHR DEM除了在空间分辨率上具有无可比拟的优势外，获取方式的灵活性也使其获得了广泛的应用，给传统古地震研究带来了新的启示和研究策略.

3.   讨论与结论

本文主要介绍了现阶段VHR DEM的提取技术，并结合最新研究成果概述了VHR DEM在活断层研究中的应用，得到如下主要结论：

1) 高分遥感立体测量、InSAR技术、LiDAR技术及SfM技术是现阶段获取VHR DEM的有效手段.机载LiDAR技术作为目前的研究热点，能够去除植被影响并获得真实的地面高程信息，高密度点云生成的VHR DEM精度高，因此对通过其它手段提取的VHR DEM进行精度评价时，常以LiDAR DEM为依据.机载LiDAR在地震活断层研究中的应用最广泛，最具发展前景，但高成本对其也有所制约；SfM和高分遥感立体测量是近几年才被应用于地震活断层的研究，因此研究成果较少，但其成本比LiDAR要低得多，只是受天气、植被、地形等影响较大；InSAR技术在获取DEM方面一直备受关注，尤其是TeraSAR-X/TanDEM-X双星系统的建成，将给全球高精度DEM的构建及应用研究带来更多突破.目前，高精度真实地面模型的建立主要依靠LiDAR技术和InSAR技术.

从发展前景来看：LiDAR技术仍然保持着迅猛的发展速度，其数据后处理算法是重要的研究方向，系统运行成本的降低也是该领域研究人员所期望的；SfM技术具有广阔的发展和应用空间，因其数据采集端基于航空摄影测量平台，所以系统构建比较容易，主要的发展方向仍在数据点云的后处理上；高分遥感立体测量基于卫星平台，覆盖范围广，系统运行稳定，高空间分辨率、高精度、短重访周期测绘卫星的研制和无控制点立体测量算法的研究是未来主要的发展方向.

2) VHR DEM亚米级的高空间分辨率在活断层研究中已发挥了独到的优势，逐渐成为一种新的研究手段，使基于地形地貌的活断层研究真正从宏观尺度进入到精细化微观尺度.相关研究人员通过实验得出VHR DEM能精确地识别并测量微小陡坎、冲沟等断错地貌；在茂密植被覆盖下的山区、建筑物密集的城市地区，利用LiDAR DEM能准确地发现“盲断层”，可用于精细地质断层的填图工作，有利于深入理解活断层的几何学、运动学及动力学等特征；极大地减少无人区、地形复杂山区等地的野外观测工作，避免不必要的危险，提高工作效率.获取同震位移场是VHR DEM差分技术的一个新的应用方向，未来大量VHR DEM数据的采集将使VHR DEM差分技术在地震形变研究中得到广泛应用，但该技术主要受VHR DEM精度的影响.

高精度VHR DEM在活断层几何结构定量分析方面的应用可作为未来研究的重要方向.现阶段，基于构造地貌理论的活断层定量化研究，大多是从高精度数字高程模型上进行的“类实地测量”，需要经验丰富的地质工作者主观识别地貌标志，其实质是基于VHR DEM的人工精细测量及相关统计分析，VHR DEM在这方面所体现的优势是前所未有的高空间分辨率、垂直方向测量和提供“大量”可供测量的断错地貌.随着VHR DEM提取技术和精度的提高，可进行适宜尺度的活断层自动识别和参数自动提取，如断层走向、倾向、倾角等，研究断层分级、组合形式等，当然这需要相关研究理论的支持和实验验证，虽然现阶段这方面的研究较少.地震形变方面的研究对数据质量要求较高，往往需要地震发生前的VHR DEM数据，因此有必要开展大范围主要活动断裂带的VHR DEM数据采集工作，这不仅会推动相关科学研究的进步，更重要的是为地震等地质灾害应急和监测预报提供宝贵资料.现阶段，VHR DEM获取技术众多，部分已比较成熟，在实际应用时需综合考虑研究区本身的地形特点、数据处理的难度及精度，选择最有利的方法. VHR DEM不仅为活断层研究带来了新的启示和创新，也将在其它地球科学研究中发挥重要作用.