您的当前位置:首页正文

ArcGIS在水源地污染源强度空间分析中的应用

2022-08-04 来源:汇智旅游网
ArcGIS在水源地污染源强度空间分析中的应用

殷绪华;朱亮;陈琳;吴新明;方长风;朱月华;邱云鹏;谭相君

【摘 要】以南京市溧水区方便水库为研究对象,运用ArcGIS对其汇水区域进行子流域划分,同时结合遥感技术对各子流域污染源强度和空间分布进行模拟,确定方便水库各入库河流主要污染源及汇水区域内污染源强度分布情况,然后利用实际监测值对ArcGIS模拟值进行验证.结果表明:ArcGIS技术可直观反映各入库河流汇水面积大小及各污染源强度的空间分布情况;ArcGIS估算值与实测值计算的污染物入库量的误差范围在10%以内.

【期刊名称】《河海大学学报(自然科学版)》 【年(卷),期】2018(046)005 【总页数】7页(P395-401)

【关键词】ArcGIS;水文分析;子流域划分;土地利用;汇水面积;污染源强度;污染物空间分布;污染排放量;方便水库

【作 者】殷绪华;朱亮;陈琳;吴新明;方长风;朱月华;邱云鹏;谭相君

【作者单位】河海大学浅水湖泊综合治理与资源开发教育部重点实验室,江苏南京210098;河海大学环境学院,江苏南京210098;河海大学浅水湖泊综合治理与资源开发教育部重点实验室,江苏南京210098;河海大学环境学院,江苏南京210098;河海大学浅水湖泊综合治理与资源开发教育部重点实验室,江苏南京210098;河海大学环境学院,江苏南京210098;江苏省骆运水利工程管理处,江苏宿迁223800;溧水区水务局,江苏南京211200;溧水区水务局方便水库管理所,江苏南京211212;河海大学浅水湖泊综合治理与资源开发教育部重点实验室,江苏南京210098;河海大学

环境学院,江苏南京210098;河海大学浅水湖泊综合治理与资源开发教育部重点实验室,江苏南京210098;河海大学环境学院,江苏南京210098 【正文语种】中 文 【中图分类】X524

我国各类型地表水源地中,湖库型水源地的比例约为40%,是我国饮用水源的重要类型之一[1]。研究表明,我国湖库型水源地的主要污染物包括化学需氧量(COD)、总氮(TN)和总磷(TP)等,污染来源主要是点源污染、面源污染和内源污染[2-3]。其中,点源污染和面源污染主要是通过入库河流来水携带污染物进入水库[4]。传统的污染源调查方法主要采用现场采样、区域入户调查和水质监测分析等,无论是数据的时效性还是准确度都存在一定的局限,而且消耗大量的人力物力[5-7]。

ArcGIS技术由于其强大的空间制图、空间可视化、空间分析和数据管理功能,在各类规划中应用广泛[8-11]。Khalil等[12]运用数据库和GIS技术评估废弃矿场的环境影响及预测污染形成过程,开发出基于GIS的环境数据库,为其他废弃矿场的污染研究提供依据。周慧平等[13]运用GIS技术强大的可视化功能识别非点源污染的主要来源、迁移途径,预报污染产生负荷、变化趋势及对水体水质的影响。王晓燕等[14]在密云水库石匣小流域非点源污染研究中,运用GIS技术确定流域内非点源污染的分布规律和重点控制区。温兆飞等[15]利用GIS的空间叠加分析法,对三峡库区农田面源污染的关键影响因子进行筛选,获取农田面源污染典型区域,并发现不同典型区分布差异的主导因素是土壤类型。

目前,运用ArcGIS技术模拟污染源分布情况已有了相应的研究,但用其进行污染物入库量的模拟计算却很少。笔者运用ArcGIS技术结合调查数据和遥感技术对方便水库土地利用方式进行分析,同时对入库污染源源强空间分布进行模拟,并用模拟结

果计算各污染物的入库量,然后通过实际监测数据进行验证。将ArcGIS技术应用到污染物入库量的研究中,旨在为河流的污染源源强分布研究提供有效的分析手段。 1 研究区概况及监测布点 1.1 研究区概况

图1 方便水库汇水区域及各入库河流监测点布局Fig.1 Layout of monitoring points around the Fangbian Reservoir and its inflow rivers

方便水库保护区位于南京市溧水区城区东北7km处,库区地处东经119°05′07″~119°07′02″、北纬31°40′~31°43′之间。方便水库在整个汇水区域范围内已进行了保护区的划定:一级保护区范围为校核洪水位(29.15 m)以下的库区,即生态红线范围之内的区域,面积约为12.21 km2;二级保护区范围为生态红线以外所有的汇水区域(生态黄线范围),面积约为64.92 km2。 1.2 监测布点

研究范围为整个方便水库汇水区域,总面积为77.13 km2,其中主要入库河流为青龙桥河、四古凹河、西南村河、西村河、杨家坝河、谢家棚河、王家山河、后吕河这8条河流。为了验证ArcGIS模拟结果的可行性和准确性,在方便水库各入库河流上分别设置采样点(图1),根据实际监测的数据计算各污染物的入库量。监测指标有COD、TN、TP和NH3-N,水质指标分析按照《水和废水监测分析方法》[16]。 2 研 究 方 法 2.1 流量计算

根据方便水库1975—2012年降水资料,采用水文手册中的方法,对陈墩库区入库水量进行计算。频率分析计算表明,1978年为枯水年(保证率P=90%),1992年为平水年(P=50%),2002年为丰水年(P=20%)。根据上述3个代表年全年各月的降雨分配情况,查《江苏省水文手册》秦淮河山丘区月降雨-径流关系曲线图,推算得到陈墩库区枯水年、平水年、丰水年入库流量分别为3.46万m3、32.83万m3、

91.46万m3。2014年总降雨量为1 368.9 mm,为丰水年(P=20%),年入库流量按91.46万m3计算。根据ArcGIS水文分析计算出的各入库河流所占比例,推算得到各条河流的年入库流量。 2.2 子流域划分

通过ArcGIS对数字高程模型(DEM)数据进行水文分析,得到方便水库汇水区域内河流水系分布及各条入库河流的汇水面积[17-18]。分析过程为:DEM数据生成、DEM栅格数据填洼处理、流向和洼地识别处理、计算汇流累积量、计算河网栅格数据、计算河流连接数据、计算子流域栅格、转化矢量数据后合并[19]。分析过程结合流域的实际情况,选择适当的汇流累积量阈值,将整个小流域划分为8个子流域。 2.3 土地利用类型分析

土地利用数据解译是基于谷歌卫星影像遥感信息,结合实地调查和其他辅助数据,采用全数字人机交互作业方法,主要根据对图像光谱、纹理、色调等的认识结合地形图目视解译而成。同时,根据溧水方便水库汇水区域水田为主、坑塘众多的特点,采用由点到面、由易到难的方式进行数据解译。按照全国土地利用分类体系,将土地利用类型分为耕地、林地、草地、水域、城乡工矿居民用地、未利用土地等6个一级类型和24个二级类型。 2.4 ArcGIS模拟污染排放量分布

区域内污染物排放量主要是通过现场调查以及相关资料收集统计。结合土地利用图,将各种污染排放量按土地利用类型进行分类及赋值,区域重合的部分数值进行叠加,以此模拟得到整个汇水区域内各入库河流小流域的污染排放量分布情况。然后,按污染排放量属性对其进行分级,分级按照ArcGIS软件中提供的自然间断点分级法,此方法可自动识别、分组数据,使各级之间差异最大化。最后,运用ArcGIS的矢量剪切功能计算各入库河流汇水面积内的污染排放量,各入库河流小流域污染排放量及贡献率。

2.5 污染排放量分布计算结果验证

由于各个污染源产生点距离库区的距离不能确定,污染物进入库区的途径差别较大,汇入库区过程中自然降解以及削减情况未能有准确、详细的调查数据,因此污染物入库量通过污染物排放量乘以相关入库系数计算得出,最后用实际监测数据计算污染物入库量,作为对计算结果的验证。计算过程中,不考虑入库沿途物理、化学等方面的自然降解作用,从而核算出水源保护区水库最大污染物入库量,以最不利情况考虑水源保护区的水质污染问题。

a. 根据ArcGIS模拟估算污染排放量及排放系数计算污染入库量。各类污染物入库量计算公式如下:

W入库量=(W排放量-θ)β (1)

式中:W排放量——不同来源的污染物年排放量;θ——污水处理厂处理掉的污染物量,取零;β——不同来源生活污染物的入库系数。

表1 各类污染源入库系数Table 1 Storage coefficient of different pollution sources污染源CODTNTPNH3-N 工业、城镇生活0.8~1.00.8~1.00.8~1.00.8~1.0 农村生活0.1~0.20.1~0.20.1~0.20.1~0.2 农田径流0.1~0.30.1~0.30.1~0.30.1~0.3 禽畜养殖0.2~0.40.2~0.40.2~0.40.2~0.4 水土流失0.1~0.30.1~0.30.1~0.40.1~0.5

根据《全国水环境容量核定技术指南》及其他大中型水库的实际调研结果,水库各类污染源入库系数见表1。

b. 根据河流入库处的实际监测数据及河流入库流量计算污染物入库量。计算公式如下: (2)

式中:C——入库河流入库处污染物实测质量浓度;Q——年入库总流量;α——各入库河流的入库流量占比(表2)。

表2 入库河道汇水情况Table 2 Catchment condition of inflow rivers河道名汇水面积/km2α/%河道名汇水面积/km2α/% 四古凹河1.1771.81谢家棚河11.62517.91西南村河6.59210.16青龙桥河35.61754.87西村河2.7814.28王家山河3.6175.57杨家坝河1.6142.49后吕河1.8832.91

图2 方便水库汇水区河网分布Fig.2 River network distribution of catchment area in the Fangbian reservoir 3 结果与讨论

3.1 河流水系分布及汇水面积

通过ArcGIS对水文DEM数据的模拟分析,得到方便水库汇水区域内河流水系分布(图2)及各条入库河流的汇水面积(表2)。由表2可以看出青龙桥河是汇入方便水库的主要河流,青龙桥河自溧水区共和山区汇流而来,来水量占总入库水量的50%以上。该河系山区性河流,降雨后径流汇聚快,来水较猛,汇流时间短,坡面水流很快汇入河道,并挟带大量泥沙,是主要的入库水沙来源。除青龙桥河外,谢家棚河及西南村河相较于其他河流来水量较大,分别占总来水量的17.91%和10.16%。 3.2 土地利用类型分析

方便水库水源地保护区总面积约77.13 km2,其中水域面积10.25 km2,土地面积66.88 km2,各用地类型面积、所占土地面积比例以及在各级保护区内的分布见表3。从表3可知,方便水库水源地一级保护区内大部分为耕地,耕地比例高达72.4%,其次为林地,以及城乡、工矿、居民用地,所占比例均为11.76%。二级保护区内占比最高的仍为耕地(54.19%),其次为林地(32.78%)。综合整个汇水区域来看,水源保护区内耕地和城乡、工矿、居民用地总的比例过高,而林地和草地所占比例偏低。 表3 方便水库汇水区域土地利用状况Table 3 Land use status in the catchment

area of Fangbian Reservoir土地利用类型 一级保护区二级保护区整个汇水区面积/km2占比/%面积/km2占比/%面积/km2占比/%耕 地

4.2572.4033.0654.1937.3155.79林 地0.6911.7520.0032.7820.6930.94草 地0.244.090.741.210.981.47城乡、工矿、居民用地

0.6911.766.9411.387.6311.41未利用土地000.270.440.270.39 3.3 饮用水源地污染源强度计算与验证 3.3.1 水源地污染排放量

经现场调查以及溧水区水务局提供的相关资料统计,方便水库主要污染源包括工业企业和集约化畜禽养殖造成的点源污染和农村生活污水、农田径流、农村固体废弃物、水土流失以及分散式禽畜养殖造成的面源污染。考虑到方便水库水面以下底泥面积和底泥释放通量等监测困难,故在研究方便水库污染源的过程中,只考虑外源污染。方便水库保护区污染物产生量统计见表4。综合点源和面源2种污染来源,方便水库的主要污染源为面源污染,面源污染的COD、TN、TP、NH3-N排放量占各自总污染排放量的97.60%、99.67%、98.10%、98.71%。其中,农田径流污染是面源污染中的第一大污染源,这与种植面积及施肥量密切相关。

表4 方便水库保护区不同来源污染物排放量及贡献率Table 4 Emission amount and contribution rate of pollutants from different sources in the Fangbian reservoir污染类别 CODTNTPNH3-N排放量/(t·a-1)贡献率/%排放量/(t·a-1)贡献率/%排放量/(t·a-1)贡献率/%排放量/(t·a-1)贡献率/%工业污染源1.810.080.180.020.050.030.140.03集约化禽畜养殖36.161.602.730.313.031.875.831.26生活污水222.839.8662.377.086.674.1143.629.43农田径流

1830.1280.98694.3678.82134.6883.01385.1483.27农村固体废弃物

59.736.787.014.32水土流失165.207.3160.796.9010.406.4127.475.94分散畜

禽养殖3.840.170.790.090.410.250.320.07

注:工业污染源、集约化禽畜养殖为点源污染,其他为面源污染。 3.3.2 ArcGIS模拟计算汇水区域内污染物排放量分布

结合RS技术对区域内土地利用方式的调查及各类型土地利用方式污染排放系数,同时结合表4,模拟计算整个汇水区域内各入库河流小流域的污染物排放量分布情况(图3)。从图3可以直观地看出不同入库河流小流域污染物排放量存在明显差异。图3中红色区域分别表示COD排放量为175.0~325.4 t/a、TN排放量为61.2~121.2 t/a、TP排放量为12.6~23.4 t/a和NH3-N排放量为35.4~65.8 t/a,均分布在青龙桥河流域和谢家棚河流域内;其次为橙红色区域,分别表示COD排放量为90.6~175.0 t/a、TN排放量为33.7~66.2 t/a、TP排放量为6.5~12.6 t/a和NH3-N排放量为18.3~34.5 t/a,主要分布在青龙桥河流域、西村河流域、西南村河流域以及四古凹河流域小部分区域内。后吕河流域和王家山河流域内多为深绿色和浅绿色,可见这2个流域污染负荷年排放量较小。

图3 方便水库保护区污染物年排放量分布Fig.3 Annual discharge distribution of pollution in the Fangbian Reservoir

运用ArcGIS的矢量剪切功能计算各入库河流汇水面积内的污染物排放量,各入库河流小流域污染物排放量及贡献率见表5。由表5可知,青龙桥河流域及谢家棚河流域污染负荷产生量较大,主要原因:(a)二者的流域汇水面积较大,青龙桥河、谢家棚河流域汇水面积各占水库总汇水面积的46.20%和15.08%;(b)处在这2条河流汇水区域内的农田及人口较多,且苗木种植也多在这2个流域范围内,故农田径流污染及农村生活、农村固废污染等面源污染的排放量大。

表5 方便水库保护区入库河流小流域污染排放量及贡献率Table 5 Pollution emission amount and contribution rate of sub-watershed in the Fangbian Reservoir编号流域CODTNTPNH3-N排放量/(t·a-1)贡献率/%排放量/(t·a-1)贡

献率/%排放量/(t·a-1)贡献率/%排放量/(t·a-1)贡献率/%1青龙桥河流域1330.7758.88508.4157.7195.7758.88271.2858.40 2谢家棚河流域393.2117.40159.0818.0628.4917.5281.8617.62 3西南村河流域224.979.9589.1410.1216.189.9546.6210.04 4西村河流域113.205.0144.925.108.114.9923.525.06 5杨家坝河流域70.583.1227.313.104.943.0414.693.16 6王家山河流域66.412.9426.152.974.702.8913.852.98 7后吕河流域32.161.4214.231.622.371.456.891.48 8四古凹河流域28.671.2811.731.322.101.285.811.26 3.3.3 污染排放量分布计算结果的验证

按照不同污染物来源及相应的入库系数,计算出不同河流小流域污染入库量及贡献率。结果显示按照排放量及入库系数计算出来的COD、TN、TP、NH3-N年入库量分别为:903.99 t、88.10 t、14.60 t、41.60 t。青龙桥河流域对COD、TN、TP、NH3-N贡献率最大,分别达到了58.88%、57.71%、58.88%、58.40%;其次为谢家棚河,其4种污染物的贡献率分别为17.40%、18.06%、17.52%、17.62%。同时,按照实际监测数据及入库水量计算出来的COD、TN、TP、NH3-N年入库量分别为:1 000.33 t、95.63 t、14.15 t、45.41 t。此种计算方法下仍然是青龙桥河流域对COD、TN、TP、NH3-N贡献率最大,分别达到了57.40%、51.74%、54.35%、50.21%;其次为谢家棚河,其4种污染物的贡献率分别为18.23%、 16.49%、13.55%、20.13%。

从各条入库河流的贡献率结果比较中可以看出(图4),与实际监测数据计算结果相比,用污染排放总量计算出的入库河流污染负荷入库量与其较为吻合,只有贡献率较小的几条河流结果稍有不同,但是相差不大,用污染排放量总量计算出的COD、TN、TP、NH3-N总入库量分别相差10%、8%、3%和8%。

图4 不同算法下污染入库量比较Fig.4 Comparison of pollution storage under different algorithms 4 结 论

a. 在方便水库8条主要入库河流中,青龙桥河的汇水面积最大,其次为谢家棚河和西南村河。方便水库各类污染源中,面源污染所占比例最大,其对COD、TN、TP、NH3-N污染负荷的贡献率分别达97.60%、99.67%、98.10%、98.71%。 b. ArcGIS模拟的污染排放量分布结果显示:污染负荷排放量最大的区域主要集中在青龙桥河流域和谢家棚河流域,其次为西村河流域和西南村河流域。位于水库南边的青龙桥河流域对COD、TN、TP、NH3-N贡献率最大,分别达到了58.88%、57.71%、58.88%、58.40%。污染物入库量的估算值和实测值的误差范围在10%以内,进而说明用ArcGIS模拟计算污染源强空间分布具有一定的可行性和准确度。

【相关文献】

[ 1 ] 冉飞. 湖库型饮用水水源地水环境生态管理和保护[J]. 山西建筑,2011,37(34): 184-186.(RAN Fei. The water environment ecosystem management and protection to the drinking water source areas of lakes and reservoirs type[J]. Shanxi Architecture,2011,37(34): 184-186.(in Chinese))

[ 2 ] 边博,朱伟,李冰,等. 太湖流域西部地区面源污染特征及其控制技术[J]. 水资源保护,2015,31(1): 48-55. (BIAN Bo, ZHU Wei, LI Bing, et al.Characteristics and control techniques of non-point pollution in the western region of Taihu Basin[J]. Water Resources Protection,2015,31(1): 48-55.(in Chinese))

[ 3 ] 郑丙辉,付青,刘琰. 中国城市饮用水源地环境问题与对策[J]. 环境保护,2007(19): 59-61. (ZHENG Binghui, FU Qing,LIU Yan. Environmental problems and countermeasures of urban drinking water sources in China[J]. Environmental Protection,2007(19): 59-61.(in Chinese))

[ 4 ] 唐访良,张明,徐建芬,等. 千岛湖库区及其主要入库河流水中有机氯农药残留污染特征及健康风险评价[J]. 环境科学,2014,35(5): 1735-1741.(TANG Fangliang, ZHANG Ming, XU Jianfen, et al.

Pollution characteristics and health risk assessment of organochlorine pesticides (ocps) in the water of Lake Qiandao and its major input rivers[J]. Environmental Science,2014,35(5): 1735-1741.(in Chinese))

[ 5 ] 李晓燕,李壁成,安韶山. 基于GIS的小流域坡耕地调查方法探讨[J]. 水土保持通报, 2001,21(5): 37-40.(LI Xiaoyan, LI Bicheng, AN Shaoshan. Sloping farmland investigation method in small watershed based on GIS[J].Bulletin of Soil and Water Consercation,2001,21(5): 37-40.(in Chinese))

[ 6 ] 周超,高乃云,赵世嘏,等. 上海青草沙水库水质调查与评价[J]. 同济大学学报(自然科学版), 2012,40(6): 894-899. (ZHOU Chao, GAO Naiyun, ZHAO Shijia, et al.Evaluation of water quality in Qingcaosha Reservoir of Shanghai[J]. Journal of Tongji University(Natural Science), 2012,40(6): 894-899.(in Chinese))

[ 7 ] 卫丹,万梅,刘锐,等. 嘉兴市规模化养猪场沼液水质调查研究[J]. 环境科学,2014,35(7): 2650-2657. (WEI Dan, WAN Mei, LIU Rui, et al. Study on the quality of digested piggery wastewater in large-scale farms in Jiaxing[J]. Environmental Science,2014,35(7): 2650-2657.(in Chinese))

[ 8 ] 邓必平,严恩萍,洪奕丰,等. 基于GIS和DEM的东江湖流域水文特征分析[J]. 湖北农业科学,2013,52(15): 3531-3536. (DENG Biping, YAN Enping, HONG Yifeng, et al. Hydrological characteristic analysis of Dongjiang Lake Watershed Based on GIS and DEM[J]. Hubei Agricultural Science,2013,52(15): 3531-3536.(in Chinese))

[ 9 ] 仉天宇,周成虎,邵全琴. 海洋GIS数据模型与结构[J]. 地球信息科学, 2003,5(4): 25-29.(ZHANG Tianyu, ZHOU Chenghu,SHAO Quanqin. Ocean GIS data model and structure[J]. Journal of Geo-Information Science,2003,5(4): 25-29.(in Chinese))

[10] 王云鹏,闵育顺,傅家谟,等. 水体污染的遥感方法及在珠江广州河段水污染监测中的应用[J]. 遥感学报,2001,5(6): 460-465.(WANG Yunpeng, MIN Yushun, FU Jiamo, et al. Remote sensing method of water pollution and application on water pollution monitoring in Guangzhou Section of Pearl River[J]. Journal of Remote Sensing,2001,5(6): 460-465.(in Chinese)) [11] 王光明,梁秀娟,肖长来,等. GIS技术在水文水资源领域中的应用现状与发展趋势[J]. 吉林水利,2009(6): 1-5. (WANG Guangming, LIANG Xiujuan, XIAO Changlai, et al. The status of the application of GIS technology in the field of hydrology and water resources and development trend[J]. Jilin Water Resources,2009(6): 1-5.(in Chinese))

[12] KHALIL A, HANICH L, HAKKOU R, et al. GIS-based environmental database for assessing the mine pollution: a case study of an abandoned mine site in Morocco[J]. Journal of Geochemical Exploration,2014, 144: 468-477.

[13] 周慧平,葛小平,许有鹏,等. GIS在非点源污染评价中的应用[J]. 水科学进展, 2004,15(4): 441-444.(ZHOU Huiping, GE Xiaoping, XU Youpeng, et al. Application of GIS to nonpoint source pollution evaluation[J]. Advances in Water Science,2004,15(4): 441-444.(in Chinese)) [14] 王晓燕,王振刚,王晓峰. GIS 支持下密云水库石匣小流域非点源污染[J]. 城市环境与城市生态,

2003 (增刊1):29-31.(WANG Xiaoyan, WANG Zhengang, WANG Xiaofeng. A GIS-based study on loss of non-point sourcepollutions from Shixia Small Watershed, Miyun Reservoir[J]. Urban Environment & Urban Ecology, 2003 (Sup1):29-31.(in Chinese)) [15] 温兆飞,吴胜军,陈吉龙,等. 三峡库区农田面源污染典型区域制图及其研究现状评价[J]. 长江流域资源与环境, 2014, 23(12): 1684-1692. (WEN Zhaofei, WU Shengjun, CHEN Jilong, et al. Typical region mapping of agricultural non-point source pollution in the Three Gorges Reservoir Area: pattern analysis and status evaluation[J]. Resources and Environment in The Yangtze Basin,2014, 23(12): 1684-1692.(in Chinese))

[16] 国家环境保护总局《水和废水监测分析方法》编委会.水和废水监测分析方法[M].4版.北京:中国环境科学出版社,2002.

[17] 李晶,张征,朱建刚,等. 基于DEM的太湖流域水文特征提取[J]. 环境科学与管理,2009,34(5): 138-142. (LI Jing, ZHANG Zheng, ZHU Jiangang, et al. Hydrological feature extraction of Taihu Lake Basin based on DEM[J]. Environmental Science and Management,2009,34(5): 138-142.(in Chinese))

[18] 顾用红,舒振文,张火青. DEM在流域水文特征分析中的应用[J]. 人民珠江,2001,22(4): 5-6.(GU Yonghong, SHU Zhenwen, ZHANG Huoqing. Application of DEM method to analysis of hydrological features in the water basin[J]. Pearl River,2001,22(4): 5-6.(in Chinese)) [19] OMRAN A, DIETRICH S, ABOUELMAGD A, et al. New ArcGIS tools developed for stream network extraction and basin delineations using Python and java script[J]. Computers & Geosciences, 2016, 94: 140-149.

因篇幅问题不能全部显示,请点此查看更多更全内容