2. 北京市城市环境过程与数字模拟国家重点实验室培育基地, 北京 100048;
3. 首都师范大学资源环境与旅游学院, 北京 100048;
4. 北京市育英学校, 北京 100072
2. State Key Laboratory Breeding Base of Process of Urban Environment and Digital Simulation, Beijing 100048, China;
3. College of Resource Environment and Tourism, Capital Normal University, Beijing 100048, China;
4. Beijing Yuying School, Beijing 100072, China
0 引言
地面沉降是在自然和人为双重因素的影响下,由于地壳表层土体的压缩,最终导致地面标高缓慢降低的环境地质灾害[1]。这种现象通常不易被人们察觉,但当它以房屋开裂、管线折断、洪水淤积等形式体现出来时,往往已造成不可补偿的永久性经济和资源损失。目前国内外学者主要针对地面沉降监测[2-5]、成因机理[6-8]、数值模拟[9-10]及评估等方面进行了诸多的研究与讨论。在地面沉降风险评估的研究中:Ki-Dong Kim等[11]根据图像数据库构建的地形图、地质图以及挖掘隧道剖面图,对地面沉降的相关因素进行了评价;宫辉力等[12]利用16景ENVISAT ASAR (advanced synthetic aperture radar)数据,结合水文地质学与地理信息系统学,对北京平原区进行了四维多尺度的地面沉降监测,同时对引起沉降的动静载荷进行了探讨;胡蓓蓓等[13]对累计地面沉降量及近几年地面沉降速率进行了分析与叠加评价,对天津市滨海地区制定了地面沉降危险性分级标准。
尽管国内外学者对地面沉降风险评估方法的研究取得了一定的进展,但是在分析沉降数据时,由于数据量的不确定性,其数据的时效性和复杂性远远超过了当前的信息处理能力,因此,数据挖掘技术应运而生。本文应用数据挖掘中的数据场模型,根据风险系统的定义以及沉降数据的特点,以地铁站点和道路节点为指标进行地面沉降交通载荷程度评价,获得地铁站点和道路节点与沉降数据相互作用的数据势场,通过与地面沉降速率的叠加聚类分析,得到区域地面沉降交通载荷程度等级划分,以期对地面沉降的风险管理与城市的恢复建设提供有益的参考。
1 研究区概况研究区位于北京东南部平原地区,包括东城区、朝阳区、大兴区东部以及通州区西部, 地处燕山南麓,华北平原北端,位置在北京冲洪积扇中上部,属潮白河冲洪积扇和温榆河冲洪积扇交汇地带,平原区海拔范围为25~45 m。研究区地理位置为116°21′00″E—116°38′50″E、39°36′25″N—40°05′27″N,面积约为746 km2。区内流动人口约为330.27万人,2010年底辖37个街道和29个乡镇。研究区交通条件便利,多条地铁线路贯穿全区;公路主要包括首都机场高速公路、京哈高速公路、京津高速公路、六环路以及多条省道,并与乡、镇、街道公路连接成网,四通八达,如图 1所示。
2 研究方法 2.1 地面沉降监测方法永久散射体干涉测量(PS-InSAR)技术是在传统差分干涉测量技术(D-InSAR)基础上发展起来的,它能够解决传统D-InSAR的时间、空间失相关和大气延迟等影响,通过获取永久散射体(PS)点,并对PS点进行相位分析,最终获得地表形变信息[14]。PS-InSAR处理流程如图 2所示。
PS-InSAR技术的优点主要包括3个方面[15-16]:第一,可以去除由于时空不一致而引起的误差;第二,可以通过多幅SAR (synthenic aperture radar)影像进行地形误差的去除,并进一步获取区域尺度长时间序列的地面沉降信息[17];第三,可以利用多幅影像、采用不同算法进行PS点的识别,进而最大程度地减少时间和空间失相关现象的发生。
本文充分考虑了研究区多幅SAR的时空相干性,并结合后续进一步的PS处理,选取2012年7月5日的TerraSAR-X数据作为主影像、其余29幅影像为辅影像进行差分干涉和PS处理,空间、时间基线情况如图 3所示。进一步解算时间序列的形变结果,得到地面年均沉降速率,如图 4所示。
2.2 数据场原理数据场就是数据通过辐射将其能量从样本空间辐射到整个母体空间,接受数据能量并被数据辐射所覆盖的空间[18]。它描述并计算每个数据点对整个空间域作用的虚拟空间场,并能够有效地处理识别其不确定性[19]。数据场具有独立性。每个目标数据都以其数据本身为中心,向外独立释放能量,其他数据不会影响此数据的特性。如同本次研究中利用PS-InSAR技术获得的地面沉降监测点(PS点),都独立地通过沉降速率或沉降量反映很小范围内的地面形变趋势。
数据场的场强函数可以描述数据在空间中辐射能量的分布规律。度量数据场中不同位置的场强,重点讨论单个数据对概念的确定度随距离的变化。对于地面沉降等一系列大量空间数据的数据挖掘研究,更多关心的是多个空间对象在一个或一个以上属性数据值中所体现出来的空间特性。即一个数据点对空间中任何一点都有作用,每个点的数据场都是它接受所有数据的数据场的叠加,任何一点的数据场的场强都是这些数据场在该点的代数和,即在数域空间中的一点所引起的数据能量之和,称为此处数据场的势。
势函数φ(x)定义为:设在数域空间Ω ⊆Rp中,数据集D={x1, x2, …, xn}中的各数据对象为具有一定辐射亮度的场源点,且按照拟核场辐射的方式向空间中其他位置x(x ⊆Ω)辐射能量,则数据场中该位置处的势值为所有场源点的作用之和,可表示为
式中:n1是场源点数;K (x)为单位势函数;δ为影响因子,它决定对象间相互作用的有效范围;mi为对象xi的质量,满足归一化条件,即
利用PS-InSAR技术监测到的地表形变信息,基于ArcGIS空间分析平台,选取空间地统计Kriging方法对PS点的形变速率进行空间插值,得到研究区年地面沉降速率分布图(图 5a)。
由图 5a可知,2010 — 2013年间,北京地面沉降速率空间分布差异性很大,沉降不均性较强,区域内形成多个沉降漏斗,最大年沉降速率达到77.69 mm/a。在研究地面沉降速率空间演化趋势基础上,采用自然分类法,将地面沉降速率分成5类(图 5b),并统计得到各类沉降速率的面积比例和在此沉降区间内的PS点数(表 1)。
沉降速率类别 | 沉降速率/(mm/a) | 面积/km2 | 面积比例/% | PS点数 |
第一类 | 77.69~50.04 | 76.14 | 8.89% | 66 652 |
第二类 | 50.04~34.13 | 154.64 | 18.05% | 130 176 |
第三类 | 34.13~20.92 | 125.58 | 14.66% | 124 035 |
第四类 | 20.92~10.42 | 325.77 | 38.04% | 250 945 |
第五类 | 10.42~0.39 | 174.24 | 20.34% | 259 764 |
从表 1可以发现:第四类地面沉降速率(20.92~10.42 mm/a)的区域面积最大,达到325.77 km2,面积所占比例为38.04 %,而地面沉降速率区间为77.69~50.04 mm/a的第一类地面沉降速率区域面积最小,占总区域面积的8.89 %,区域面积为76.14 km2;PS点数占比例最大的为第五类,地面沉降速率处于10.42~0.39 mm/a,PS点数达到259 764个,而在第一类地面沉降速率(77.69~50.04 mm/a)内的PS点数为66 652个。整体而言,地面沉降速率的分布不均匀性较大,即不均匀沉降演化趋势明显。
3.2 研究区数据场势值分析地面沉降灾害风险包括沉降程度的大小及其对社会经济与自然环境系统造成的影响[20]。本文采用基于数据场的方法对北京地区66个乡镇街道的地面沉降交通载荷程度进行初步获取,选取地铁站信息和道路节点信息作为交通载荷程度获取的指标,其中地铁站信息为整个研究区的所有地铁站地理坐标,道路节点信息为研究区高速公路和城市一级公路交叉点的地理坐标。结合沉降监测所得到的PS点数据,应用Matlab软件对地铁站点与道路节点数据信息进行扩散,得到地铁站点以及道路节点的数据场势图,结果如图 6所示。
3.3 利用因子贡献权重法划分地面沉降交通载荷程度等级在获取地铁站点和道路节点与地面沉降信息势值的基础上,采用因子贡献权重方法计算地铁站点和道路节点指标权重,其中应用贡献率的概念量化每种指标的影响程度[21]。该方法可以定量表述各指标贡献程度的大小,并据此对地面沉降交通载荷程度等级进行划分。
3.3.1 各指标归一化处理统计地铁站点和道路节点指标的数据信息,以乡镇街道为单位区域平均联立得到各指标数值。为解决数据的可比统一性,采用归一化计算方法,将各指标数据化为无量纲表达,公式如下:
式中:Ai为各级指标的归一化数值;Xi为各指标值;n2为指标的数量。
3.3.2 各指标权重确定首先按照等差通项公式将各指标归一化值划分出5个级别:高级、较高级、中级、较低级、低级,求得各指标对应等差数列的公差,公式为
式中:D为指标数列公差; Amax为指标中的最大值; Amin为指标中的最小值。由此得出5级划分区间:低贡献率区间X1=[Amin, Amin+D];较低贡献率区间X2=[Amin+D, Amin+2D];中贡献率区间X3=[Amin+2D, Amin+3D];较高贡献率区间X4=[Amin+3D, Amin+4D];高贡献率区间X5=[Amin+4D, Amax]。各指标贡献率划分结果如表 2所示。
贡献率区间 | 沉降速率 | 地铁站点势值 | 道路节点势值 |
X1 | 0.000 3 | 0.000 0 | 0.000 0 |
0.000 4 | 0.000 0 | 0.000 0 | |
0.001 2 | 0.000 0 | 0.000 0 | |
0.001 4 | 0.000 0 | 0.000 0 | |
0.001 6 | 0.000 0 | 0.000 0 | |
0.001 8 | 0.000 0 | 0.000 1 | |
0.001 9 | 0.000 0 | 0.000 2 | |
0.002 1 | 0.000 0 | 0.000 2 | |
0.002 1 | 0.000 1 | 0.000 7 | |
0.002 2 | 0.000 2 | 0.000 8 | |
0.002 5 | 0.000 6 | 0.000 9 | |
0.002 6 | 0.001 0 | 0.000 9 | |
0.002 7 | 0.001 2 | 0.001 3 | |
0.003 1 | 0.001 3 | 0.001 8 | |
0.003 3 | 0.001 8 | 0.002 1 | |
0.004 0 | 0.002 6 | 0.002 2 | |
0.004 1 | 0.002 7 | 0.002 5 | |
0.004 1 | 0.002 8 | 0.002 8 | |
0.004 1 | 0.003 1 | 0.003 0 | |
0.004 3 | 0.003 9 | 0.003 4 | |
0.004 9 | 0.004 1 | 0.003 5 | |
0.005 5 | 0.004 9 | 0.003 6 | |
0.005 5 | 0.005 1 | 0.004 1 | |
0.006 2 | 0.005 2 | 0.004 7 | |
0.006 3 | 0.005 7 | 0.004 9 | |
0.006 4 | 0.006 5 | 0.005 2 | |
0.006 7 | 0.007 0 | 0.007 0 | |
0.006 7 | |||
0.006 8 | |||
0.007 2 | |||
0.007 5 | |||
0.008 2 | |||
0.008 3 | |||
0.008 3 | |||
0.008 5 | |||
0.009 0 | |||
0.009 0 | |||
0.009 2 | |||
0.010 3 | |||
0.010 9 | |||
X2 | 0.011 6 | 0.008 5 | 0.008 7 |
0.012 0 | 0.011 7 | 0.008 9 | |
0.012 7 | 0.012 1 | 0.010 4 | |
0.013 2 | 0.014 2 | 0.012 4 | |
0.013 3 | 0.014 3 | 0.012 8 | |
0.018 8 | 0.014 4 | 0.013 0 | |
0.018 9 | 0.015 0 | 0.013 1 | |
0.021 5 | 0.015 7 | 0.013 3 | |
0.022 3 | 0.015 8 | 0.013 6 | |
0.015 9 | 0.013 7 | ||
0.016 0 | 0.014 9 | ||
0.015 2 | |||
0.016 9 | |||
X3 | 0.023 2 | 0.016 5 | 0.018 5 |
0.026 4 | 0.017 1 | 0.020 7 | |
0.027 4 | 0.017 1 | 0.020 7 | |
0.032 9 | 0.018 2 | 0.023 0 | |
0.032 9 | 0.019 0 | 0.024 0 | |
0.033 1 | 0.021 9 | 0.024 4 | |
0.022 1 | 0.024 5 | ||
0.024 9 | |||
0.025 9 | |||
X4 | 0.033 9 | 0.024 4 | 0.026 1 |
0.035 6 | 0.024 4 | 0.026 2 | |
0.037 4 | 0.024 7 | 0.026 7 | |
0.038 4 | 0.025 0 | 0.027 7 | |
0.038 9 | 0.026 4 | 0.028 1 | |
0.043 6 | 0.027 0 | 0.029 7 | |
0.028 7 | 0.030 2 | ||
0.029 1 | 0.030 9 | ||
0.031 1 | 0.031 0 | ||
0.031 1 | 0.032 2 | ||
X5 | 0.045 7 | 0.032 0 | 0.034 9 |
0.046 4 | 0.032 5 | 0.036 2 | |
0.049 2 | 0.032 5 | 0.037 8 | |
0.053 9 | 0.032 9 | 0.042 8 | |
0.055 6 | 0.033 2 | 0.043 2 | |
0.033 6 | 0.043 4 | ||
0.033 8 | 0.043 4 | ||
0.036 9 | |||
0.037 0 | |||
0.038 2 | |||
0.040 0 |
然后对每一个指标的贡献率按5级划分后进行贡献率均值化,通过对不同区间的所有值进行加权平均,得到高、较高、中、较低、低贡献率均值:
式中:Ah, Au, Am, As, Al分别为高、较高、中、较低、低贡献率各区间对应指标和的平均值;Xi为各区间对应指标的平均值;n3为高、较高、中、较低、低贡献率区间中各指标数。据此得到各指标各级贡献率的均值数据(表 3)。
沉降速率 | 地铁势值 | 道路势值 | |
Ah | 0.050 1 | 0.034 7 | 0.040 2 |
Au | 0.037 9 | 0.027 1 | 0.028 8 |
Am | 0.029 3 | 0.018 8 | 0.022 9 |
As | 0.016 0 | 0.014 0 | 0.012 8 |
Al | 0.005 0 | 0.002 2 | 0.002 0 |
最后,计算地铁站点和道路节点指标对地面沉降交通载荷的贡献比重,用贡献比重来代表地面沉降交通载荷影响指标的权重。表达式为
表达式中:ws为地铁站点指标的权重;wr为道路节点指标的权重;Ave (s)是5个区间地铁势值的平均值;Ave (v)是5个区间沉降速率的平均值;Ave (r)是5个区间道路势值的平均值。计算得到地铁势值指标和道路势值指标的权重分别为0.411 9和0.435 9。
应用ArcGIS分别对研究区地面沉降交通载荷指标进行空间分析,采用叠加分析对各指标重分类的栅格图层赋以分级权重,得到地面沉降交通(图 7)和道路节点(图 8)载荷程度区划结果。
从图 7可知:研究区地铁站高载荷程度区主要位于平房乡、楼梓庄乡、管庄乡、崔各庄乡和常营乡等;较高载荷程度区主要位于东坝乡、王四营乡、东风乡和高碑店乡等。这些地区主要有地铁1号线、6号线、10号线以及机场地铁专线贯穿其中,地铁站相对密集,人口流动性较强。其中地面沉降灾害对管庄、双桥、传媒大学、四惠东站、国贸以及青年路站作用程度较严重。
从图 8可知:研究区道路节点高载荷程度区主要位于崔各庄乡、楼梓庄乡、常营乡和三间房乡等;较高载荷程度区主要位于东四街道、建国门街道和朝阳门街道等。这些区域均属于北京市区交通要道处,车流量较大,其中来广营桥、四惠桥、永定门桥等主路段受沉降灾害影响较大。
综上,地面沉降地铁站点和道路节点高载荷区域主要分布在北京市朝阳区北部与中部区域,而低载荷区域主要位于远离城市重要交通道路和地铁线路不发达乡镇。
统计地铁站点和道路节点载荷程度结果如表 4所示。从表 4可以看到:研究区内地铁站高载荷程度区、较高程度区和中等程度乡镇分别有7个、15个和23个,占总乡镇的68.19%;而道路节点高载荷程度区、较高程度区和中等程度乡镇分别有5个、8个和15个,占总乡镇的42.43%。这在一定程度上说明,地面沉降交通载荷对研究区地铁站的影响程度比对交通道路的程度大。
地铁站点载荷程度 | 乡镇数 | 占总乡镇比例/% |
高 | 7 | 10.61 |
较高 | 15 | 22.73 |
中等 | 23 | 34.85 |
较低 | 14 | 21.21 |
低 | 7 | 10.61 |
高 | 5 | 7.58 |
较高 | 8 | 12.12 |
中等 | 15 | 22.73 |
较低 | 21 | 31.82 |
低 | 17 | 25.76 |
1) 研究区地面沉降速率空间分布差异性很大,地面沉降呈现不均匀,且不均性较强,研究区内形成多个沉降漏斗,最大年沉降速率达到了77.69 mm/a。
2) 在66个乡镇街道中,地铁站交通载荷程度3级(中等载荷程度)以上的乡镇街道就有45个,占68.19%;道路节点交通载荷程度3级(中等载荷程度)以上的乡镇街道28个,占42.43%。
3) 研究区内高等载荷程度主要分布在朝阳区北部与中部区域,这些乡镇街道区域经济发达,铺设地铁线路集中,人口密集,单位路段车辆承载量较大,而低等交通载荷程度的地区位于远离重要交通道路与地铁线路的经济欠发达乡镇。
本文侧重利用数据场模型建立地面沉降交通载荷程度的评价指标,通过因子贡献权重方法分析方法获取研究区地面沉降交通载荷程度分区情况。但本次研究中,选择地铁站点和道路节点作为地面沉降交通载荷的评价指标略显单一,在后期研究中可加入更多指标对地面沉降交通载荷程度进行评价。
[1] | Deng Zeng, Ke Yinghai, Gong Huili, et al. Land Subsidence Prediction in Beijing Based on PS-InSAR Technique and Improved Grey-Markov Model[J]. GIScience & Remote Sensing, 2017: 1-22. |
[2] |
陈蓓蓓, 宫辉力, 李小娟, 等. 北京地下水系统演化与地面沉降过程[J].
吉林大学学报(地球科学版), 2012, 42(1): 373-379.
Chen Beibei, Gong Huili, Li Xiaojuan, et al. Groundwater System Evolution and Land Subsidence Process in Beijing[J]. Journal of Jilin University (Earth Science Edition), 2012, 42(1): 373-379. |
[3] | Lesniak A, Porzycka S. Environment Monitoring Using Satellite Radar Interferometry Technique (PSInSAR)[J]. Polish Journal of Environmental Studies, 2008, 17(3A): 382-387. |
[4] | Kim S W, Wdowinski S, Dixon T H, et al. Measurements and Predictions of Subsidence Induced by Soil Consolidation Using Persistent Sscatterer InSAR and a Hyperbolic Model[J]. Geophysical Research Letters, 2010, 37(5): 87-105. |
[5] | Ciampalini A, Raspini F, Lagomarsino D, et al. Landslide Susceptibility Map Refinement Using PSInSAR Data[J]. Remote Sensing of Environment, 2016, 184: 302-315. DOI:10.1016/j.rse.2016.07.018 |
[6] | Budhu M, Adiyaman I B. Mechanics of Land Subsi-dence Due to Groundwater Pumping[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2010, 34(14): 1459-1478. |
[7] | Bakr M. Influence of Groundwater Management on Land Subsidence in Deltas[J]. Water Resources Management, 2015, 29(5): 1541-1555. DOI:10.1007/s11269-014-0893-7 |
[8] |
贾三满, 王海刚, 赵守生, 等. 北京地面沉降机理研究初探[J].
城市地质, 2007, 2(1): 20-26.
Jia Sanman, Wang Haigang, Zhao Shousheng, et al. A Tentative Study of the Mechanism of Land Subsidence in Beijing[J]. City Geology, 2007, 2(1): 20-26. |
[9] |
付延玲, 骆祖江, 廖翔, 等. 高层建筑引发地面沉降模拟预测三维流固全耦合模型[J].
吉林大学学报(地球科学版), 2016, 46(6): 1781-1789.
Fu Yanling, Luo Zujiang, Liao Xiang, et al. A Three-Dimensional Full Coupling Model to Simulate and Predict Land Subsidence Caused by High-Rise Building[J]. Journal of Jilin University (Earth Science Edition), 2016, 46(6): 1781-1789. |
[10] |
骆祖江, 刘金宝, 李朗. 第四纪松散沉积层地下水疏降与地面沉降三维全耦合数值模拟[J].
岩土工程学报, 2008, 30(2): 193-198.
Luo Zujiang, Liu Jinbao, Li Lang. Three-Dimensional Full Coupling Numerical Simulation of Groundwater Dewatering and Land-Subsidence in Quaternary Loose Sediments[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(2): 193-198. |
[11] | Kim K, Lee S, Oh H, et al. Assessment of Ground Subsidence Hazard near an Abandoned Underground Coal Mine Using GIS[J]. Environmental Geology, 2006, 50(8): 1183-1191. DOI:10.1007/s00254-006-0290-5 |
[12] |
宫辉力, 张有全, 李小娟, 等. 基于永久散射体雷达干涉测量技术的北京市地面沉降研究[J].
自然科学进展, 2009, 19(11): 1261-1266.
Gong Huili, Zhang Youquan, Li Xiaojuan, et al. The Research of Land Subsidence in Beijing Based on Permanent Scatterers Interferometric Synthetic Aperture Radar(PS-InSAR) Technique[J]. Progress in Natural Science, 2009, 19(11): 1261-1266. DOI:10.3321/j.issn:1002-008X.2009.11.015 |
[13] |
胡蓓蓓, 姜衍祥, 周俊, 等. 天津市滨海地区地面沉降灾害风险评估与区划[J].
地理科学, 2008, 28(5): 693-697.
Hu Beibei, Jiang Yanxiang, Zhou Jun, et al. Assessment and Zonation of Land Subsidence Disaster Risk of Tianjin Binhai Area[J]. Scientia Geographica Sinica, 2008, 28(5): 693-697. |
[14] |
万小琴, 胡波, 马清林. 基于PS-InSAR技术的澳门地表沉降研究[J].
测绘工程, 2012, 21(3): 39-43.
Wan Xiaoqin, Hu Bo, Ma Qinglin. Monitoring Ground Subsidence in Macau with Persistent Scatters Interferometry[J]. Engineering of Surveying and Mapping, 2012, 21(3): 39-43. |
[15] | Colesanti C, Ferretti A, Prati C, et al. Monitoring Landslides and Tectonic Motions with the Permanent Scatterers Technique[J]. Engineering Geology, 2003, 68(1): 3-14. |
[16] | Teatini P, Tosi L, Strozzi T, et al. Mapping Regional Land Displacements in the Venice Coastland by an Integrated Monitoring System[J]. Remote Sensing of Environment, 2005, 98(4): 403-413. DOI:10.1016/j.rse.2005.08.002 |
[17] |
雷坤超, 贾三满, 陈蓓蓓, 等. 基于PS-InSAR技术的廊坊市地面沉降监测研究[J].
遥感技术与应用, 2013(6): 1114-1119.
Lei Kunchao, Jia Sanman, Chen Beibei. Land Subsidence Detection Based on PS-InSAR in Langfang[J]. Remote Sensing Technology and Application, 2013(6): 1114-1119. |
[18] |
李德毅, 刘常昱, 杜鹢, 等. 不确定性人工智能[J].
软件学报, 2004, 15(11): 1583-1594.
Li Deyi, Liu Changyu, Du Yi, et al. Artificial Intelligence with Uncertainty[J]. Journal of Software, 2004, 15(11): 1583-1594. |
[19] |
王树良, 邹珊珊, 操保华, 等. 利用数据场的表情脸识别方法[J].
武汉大学学报(信息科学版), 2010, 35(6): 738-742.
Wang Shuliang, Zou Shanshan, Cao Baohua, et al. Facial Expression Recognition Based on Data Field[J]. Geomatics and Information Science of Wuhan University, 2010, 35(6): 738-742. |
[20] |
蒋新宇, 范久波, 张继权, 等. 基GIS的松花江干流暴雨洪涝灾害风险评估[J].
灾害学, 2009, 24(3): 51-56.
Jiang Xinyu, Fan Jiubo, Zhang Jiquan, et al. GIS-Based Risk Assessment on Rain and Flood Disasters of Songhua River[J]. Journal of Catastrophology, 2009, 24(3): 51-56. |
[21] |
鲍新华, 张宇, 李野, 等. 松辽盆地增强型地热系统开发选区评价[J].
吉林大学学报(地球科学版), 2017, 47(2): 564-572.
Bao Xinhua, Zhang Yu, Li Ye, et al. Evaluation of Development Selection for Enhanced Geothermal System in Songliao Basin[J]. Journal of Jilin University (Earth Science Edition), 2017, 47(2): 564-572. |