A discussion on the evaluation method of geological environment carrying capacity of Hebei Province
-
摘要:
随着近几十年河北省经济社会的快速发展、人口增长和城镇化水平的不断提高，产生了大量崩塌、滑坡、泥石流、地面塌陷、地裂缝、地面沉降等地质灾害，影响了河北省经济社会又好又快发展。产生这些地质灾害的主要原因就是人类活动强度超出了当地的地质环境承载能力。为客观分析河北省地质环境承载能力，建立了河北省地质环境承载能力评价指标，并分别进行了河北省地质环境承载能力本底和状态评价。通过对评价结果的研究分析，提出了新的地质环境承载能力状态评价指标和方法，为将来建立多方协同的资源环境承载能力监测预警长效机制提供了参考。
Abstract:With the rapid economic and social development, population growth and urbanization in Hebei Province in the past decades, geological disasters such as collapse, landslides, mudslides, ground subsidence, ground fissures and land subsidence have been increasingly generated, which affects economic development of Hebei. The province's economy and society are developing well and fast.The main reason for the formation of these geological disasters is that the intensity of human activities exceeds the carrying capacity of the local geological environment.To objectively analyze the geological environment carrying capacity of Hebei Province, the authors established the evaluation index of the regional geological environment carrying capacity of Hebei Province, and carried out the geological environment carrying capacity background evaluationand status evaluation of Hebei Province. Throughthe study and analysis of the evaluation results, the authors put forward a new evaluation index and method for the geological environment carrying capacity status evaluation so as to beter conform to the actual situation of the study area, and carried out an analysis and evaluation so as to establish a multi-party coordinated resource and environmental carrying capacit monitoring and early warning long-term effect. The results obtained by the authors provide a reference for these purposes.
-
石英脉型黑钨矿床是重要的钨矿床类型之一，主要分布于南岭成矿带的赣南、湘南、粤北、桂北等地[1-2]。钨矿脉主要分布在碱长花岗岩顶部接触带附近[3]。大多数钨矿脉近垂直产出，垂向延伸可达1000m[2]。石英型黑钨矿床矿脉围岩蚀变呈带状分布在矿脉两侧，围岩蚀变宽度随深度递减，且这种变化趋势不受岩性影响[1, 4]。多位学者认为，形成黑钨矿石英脉的成矿流体不是纯热水溶液，而是富含SiO2和挥发分的岩浆热液过渡性流体[4-7]。这种岩浆热液过渡性流体比水溶液具有更高的粘度和密度。形成钨矿脉的过渡性流体在深部偏向熔体，在浅部偏向热液，因而蚀变宽度上宽下窄。Audétat等[8]的实验表明，只有当熔体的摩尔分数很高时，含水熔体的粘度才会显著高于纯水的粘度。然而，除常海亮等[9]在西华山钨矿脉的绿柱石中发现熔体包裹体外，大多数矿石矿物和脉石矿物的流体包裹体研究结果表明，钨成矿流体是NaCl-H2O±CO2体系[10-22]。因而，形成上述蚀变特征可能另有他因。岩浆热液矿床的形成涉及复杂的化学反应和热液运移过程[23-25]。前人研究多偏重于前者，对热液运移过程的研究较少[26-28]。本文利用有限元方法模拟热液运移和化学组分从裂隙向邻近围岩扩散过程，发现围岩孔隙度随深度增加而递减可能是钨矿床形成上述蚀变特征的有效机制。
1. 石英脉型钨矿床地质和地球化学特征
中国石英脉型钨矿床主要分布于南岭成矿带的赣南、湘南、粤北、桂北等地（图 1）。南岭地区大规模钨成矿作用与燕山早期的花岗岩岩浆活动有关，成矿年龄大多集中在中晚侏罗世(165~ 150Ma) [29]。该类矿床常产在碱长花岗岩顶上带附近[3]。主要的矿石矿物为黑钨矿、白钨矿、锡石等，脉石矿物有石英、长石、萤石、黄铁矿等[1]。
![]()
图 1 南岭成矿带钨矿分布示意图[29]Figure 1. Distribution of tungsten deposits in the Nanling Mountains, South China钨成矿流体主要属于中高温、中低盐度的NaCl-H2O±CO2体系[30]。首先，成矿流体成分以H2O为主，CO2含量低，少量CO、CH4、N2、H2。在主成矿阶段成矿流体主要来自岩浆水，在成矿晚期成矿流体中混有大量大气降水。其次，寄主矿物（石英、黑钨矿、绿柱石、锡石等）捕获的流体包裹体在常温下以富液相的气液两相包裹体为主，纯气相和含CO2的三相包裹体少见，表明溶液不曾沸腾或局部沸腾，压力可达90~160MPa。流体包裹体的均一温度一般为160~390℃，矿石矿物（黑钨矿、锡石等）的均一温度比石英的高40~100℃。最后，成矿流体盐度较低, 成矿早期流体盐度一般为7%~10%NaCl，晚期逐渐降低至1%[1-2, 10-22, 31-35]。
钨矿床围岩类型以寒武系和奥陶系居多，岩性以沉积(变质)碎屑岩为主[2, 36]。蚀变类型包括云英岩化、硅化、黑云母化、电气石化等[2]。含钨石英脉的围岩蚀变有以下特点[1]：①围岩蚀变呈带状分布在矿脉两侧（图 2），蚀变带宽几毫米至几十厘米，带状分布特点不受岩性影响；②围岩蚀变宽度上宽下窄，蚀变宽度变化趋势不受岩性影响；③矿脉以石英为主，但石英脉体两侧的硅化较弱，且上部硅化比下部强。
![]()
图 2 南岭地区石英脉型钨矿典型矿脉a—江西大吉山钨矿深部脉体，标高417m，脉体走向为近东西向，陡倾，向上拍摄；b—江西漂塘钨矿深部大脉，标高268m，脉体走向为近东西向，陡倾，向上拍摄。Qtz—石英；Wol—黑钨矿Figure 2. Typical veins at deeper levels of tungsten deposits in the Nanling Mountains反应性流体（如成矿热液）流经岩石裂隙或孔隙时，会改变邻近围岩的化学成分，形成蚀变晕[37]。蚀变晕的宽度与岩石、流体的物理化学性质、反应速率、反应时间等参数密切相关[37-39]。本文利用热液运移数值模拟研究化学组分从裂隙向邻近围岩扩散过程，为解释石英脉型钨矿床蚀变特征提供新思路。
2. 热液运移数值模拟
2.1 数学理论
多孔介质流体流动和热传输用流体连续性方程、达西定律和热传输方程定量描述，这些偏微分方程组用邢会林教授开发的PANDAS (Parallel Adaptive Nonlinear Deformation Analysis System)来求解。PANDAS是一套基于有限元方法(FEM)和格子Boltzmann方法(LBM)开发的创新性软件平台，已经成功用于模拟地下含裂隙非均质孔隙材料/地层中“应力变形/破坏-流体流动-热传导-化学反应”等多物理场高度非线性耦合的问题，解决从微孔隙尺度到实验室及油气田尺度，甚至全球等多尺度的地球科学和地质资源工程中的主要科学问题和工程技术挑战[40-49]。
NaCl水溶液的密度和粘度与温度、压力和盐度有关。Batzle等提出的盐水密度模型在热液条件下精度较高[50]，因而本文的模拟计算使用了该密度模型。笔者提出一个适于大规模数值计算的盐水粘度模型[43]，用于计算热液的粘度。
某一组分在多孔介质的有效扩散速率De与该组分在水溶液的扩散速率Ds和孔隙度ϕ有关[51]：
D_{e}=D_{s }ϕ ^{2} (1) 扩散速率Ds服从Stokes-Einstein方程[52-53]。
Cathles等[37]提出垂直于脉体的蚀变晕宽度计算公式：
Z（t） = \sqrt {\frac{{ 2D_{e}t}}{{ G}}} (2) 其中，Z为时间t内形成的蚀变晕宽度，G是一个无量纲数，其值为改变1m3岩石所需要的流体体积。本文利用Péclet数（Pe）表征流体传输机制[54]，当Pe < 1时，流体流动的扩散比例较高，流体内的化学反应以反应-扩散过程为主；当Pe > 1时，流体流动的对流比例升高，化学反应以反应-扩散-对流过程为主[55]。
2.2 热液运移模型
本文主要根据以下背景资料建立热液运移模型：①流体包裹体测温研究表明，钨成矿流体的压力较高，可达到90~160MPa，平均值为125MPa；②钨成矿热液来自深部含水岩浆房，岩浆顶上带与围岩接触部位发育大量的裂隙，这些裂隙为热液运移提供构造通道和成矿空间[56-58]。
根据石英脉型钨矿床的地质地球化学特征，本文建立了二维热液运移模型（图 3）。Z轴代表重力方向，X轴为垂直脉带的水平方向。为了减少网格数量和计算时间，本文模型只截取深度3.8~5.0km的部分，模型尺寸为1.2km×1.2km。模型包含裂隙区和围岩2个单元，石英脉型钨矿床脉带宽度从几十米至几百米不等[2]，模型中裂隙区宽度为20m。选取2条剖面用于研究裂隙带及其紧邻围岩的热液运移，其中一条剖面在4.9km，另一条剖面在4.1km处。
![]()
图 3 石英脉型钨矿床热液运移数值模型Figure 3. Numerical model of hydrothermal flow at the tungsten deposits in the Nangling Mountains设定模型底部边界条件为压力125MPa和温度400℃，顶部边界条件为压力38MPa和温度114℃。根据现今地热储库观测数据[59]，设定初始温度梯度为30℃/km，初始压力梯度为10MPa/km，模型中热液盐度为10% NaCl。硅化是南岭石英脉钨矿床较普遍的蚀变类型之一，本文将硅从裂隙到邻近围岩的扩散作为研究对象。硅在25℃稀溶液的扩散系数约为1.0× 10-9m2/s[60]。
本文设计2个数值实验，2个实验中裂隙带比围岩拥有更高的孔隙度和渗透率。在实验1中，围岩和裂隙带的孔隙度和渗透率均一分布，裂隙带孔隙度为10%，渗透率为1.0 × 10-12m2；围岩孔隙度为0.01，渗透率为1.0× 10-16m2。在实验2中，围岩的孔隙度和渗透率与深度呈负指数关系[61]，渗透率在模型截取的深度范围内变化约一个数量级。
3. 模拟结果
3.1 数值实验1
在数值实验1中，围岩的孔隙度和渗透率均一分布。4.9km剖面和4.1km剖面的初始温度分别为147 ℃和123℃。热液运移1a后，裂隙带的温度显著高于邻近围岩的温度，但远离裂隙带的围岩温度变化较小（图 4），在4.9km剖面上，紧邻裂隙带的围岩温度达到340℃；在4.1km剖面上，紧邻裂隙带的围岩温度为320℃。在2条剖面中，裂隙带内的Péclet数为1.7~1.8，围岩的Péclet数远低于1。
![]()
图 4 数值实验1温度(a)和Péclet数(b)的分布Figure 4. The temperature (a) and Péclet number (b) in the first numerical experiment由于4.9km剖面的初始温度高于4.1km剖面的初始温度，前者的初始有效扩散系数高于后者。在紧邻裂隙的围岩内，硅有效扩散系数在热液运移1a后有显著提升：在4.9km剖面最大值达2.3×10-13m2/s，在4.1km剖面最大值为2.1×10-13m2/s（图 5）。热液运移2a后，围岩在4.1km剖面的有效扩散系数仍然低于围岩在4.9km的有效扩散系数。
![]()
图 5 数值实验1热液运移1a后(a)和2a后(b)硅有效扩散系数在裂隙带附近的分布Figure 5. The effective diffusion coefficients of silica in the wallrock adjacent to fracture zones after one year (a) and two years (b) in the first numerical experiment3.2 数值实验2
在数值实验2中，围岩的孔隙度和渗透率随深度增加而按指数规律递减，其他参数与数值实验1一致。在该实验中，围岩孔隙度从3.8km的1.17%降至5.0km的0.58%，对应的渗透率由10-15.8m2降至10-16.7m2。围岩在4.9km的孔隙度为0.64%，渗透率为10-16.6m2。围岩在4.1km的孔隙度为0.99%，渗透率为10-16.0m2。
热液运移1a后，紧邻裂隙带的围岩在4.9km的温度达到315℃，而相应的温度在4.1km为230℃（图 6）。显著低于数值实验1中紧邻裂隙带的围岩温度，裂隙带的Péclet数为1.8~2.0，围岩的Péclet数低于1，说明实验2中热液运移的流体传输机制与实验1相似。
![]()
图 6 数值实验2温度(a)和Péclet数(b)的分布Figure 6. The temperature (a) and Péclet number (b) in the second numerical experiment该实验的初始有效扩散系数与实验1有明显区别。围岩在4.9km的初始有效扩散系数为3.2 × 10-14m2/s，约为4.1km剖面的初始有效扩散系数(6.0 × 10-14m2/s)的一半。热液运移1a后，围岩在4.1km的有效扩散系数最大值为1.4×10-13m2/s，是4.9km剖面对应值的1.7倍（图 7）。热液运移2a后，围岩在4.1km的有效扩散系数最大值升高至2.2× 10-13m2/s，是4.9km剖面对应值的2.4倍。
![]()
图 7 数值实验2中热液运移1a后(a)和2a后(b)硅有效扩散系数在裂隙带邻近围岩的分布Figure 7. The effective diffusion coefficients of silica in the wallrock adjacent to fracture zones after one year (a) and two years (b) in the second numerical experiment4. 讨论
前人在20世纪80年代发现南岭地区石英脉型钨矿床蚀变宽度呈上宽下窄的特征[1]。本文利用方程（2）对比不同深度下的蚀变宽度，讨论蚀变宽度上宽下窄的形成机制。由于成矿流体主要来自深部岩浆，假设方程（2）的参数G是一个常数，在给定时间内垂直脉体的硅化宽度与硅的有效扩散系数呈正比例关系。
在数值实验1中，热液运移在紧邻裂隙带的围岩内以扩散为主，反应-扩散过程在化学反应中占主导地位。由于围岩的孔隙度和渗透率是均一分布的常数，硅在围岩的有效扩散系数只决定于硅在水溶液的扩散系数（见公式1）。围岩在深部的温度高于浅部，因而硅在深部围岩中有更高的有效扩散系数，深部易形成较宽的硅化，这显然与石英脉型钨矿床蚀变特征相反。事实上，前人的数值模拟也得到类似的结论：蚀变宽度随着热液运移远离流体源头而递减[37, 39]。然而，这与南岭地区石英脉型钨矿和某些斑岩型铜矿的蚀变特征不符[38]。通过数值实验1，认为石英脉型钨矿在1km垂直范围内围岩孔隙度可能不是均一分布的。
在数值实验2中，围岩孔隙度和渗透率随深度增大呈指数递减。尽管深部围岩的温度依然高于浅部，由于硅有效扩散系数与孔隙度的平方成正比，随深度递减的孔隙度抵消了温度对于硅有效扩散系数的影响，因而浅部围岩中硅的有效扩散系数高于深部。在这种情况下，裂隙带周围可能形成与南岭石英脉型钨矿床类似的蚀变特征。
本文数值模拟实验的运行时间远小于岩浆热液系统演化的时间尺度，后者可能在数个百万年甚至几十个百万年的时间尺度保持活跃状态[62-63]。热液矿床成矿动力学过程涉及多时空尺度相互作用的热液流动、岩石变形和化学反应。受岩石变形、化学反应等过程的影响，与矿床形成相关的高渗透率具有局部性、短暂性、周期性等特征[64-66]。多方面证据表明，上地壳高渗透会在1~1000a内发生显著衰减[67]。故本文的模拟实验可看作成矿裂隙渗透率近似不变时部分热液充填过程。另外，在本文模型边界条件的控制下，延长模拟实验运行时间会提高围岩整体的温度和硅的扩散速率，但不会改变这2个变量在深部与浅部的相对大小。故数值模拟实验运行时间虽短，但对理解钨矿床蚀变晕形成机制仍具有一定的启发。
石英脉型钨矿床的围岩岩性可能有利于形成随深度递减的孔隙度和渗透率，南岭地区石英脉型钨矿床围岩岩性以沉积(变质)碎屑岩为主[2, 29]。这些围岩较致密，孔隙度也较低。实验表明，由于微小裂隙的闭合，低孔隙度结晶岩石的渗透率随着压力的增加而呈指数规律减小；而高孔隙度沉积岩对压力的敏感性较弱[68-69]。因而，石英脉型钨矿床围岩岩性可能有利于形成随深度递减的孔隙度和渗透率。野外证据表明，石英脉型钨矿床围岩的裂隙发育程度随深度增加而递减[56-57]，这可能是围岩孔隙度和渗透率随深度递减的另一个原因。
5. 结论
通过热液运移和组分扩散数值模拟实验，研究硅从裂隙带向邻近围岩的扩散过程。数值模拟实验使用的流体是NaCl水溶液，其密度和粘度与温度、压力和盐度有关。本文的数值模拟实验为石英脉型钨矿床蚀变宽度垂直分带形成机制提供如下启示。
（1）温度和孔隙度在垂向上的空间分布是影响石英脉型钨矿床蚀变特征的重要变量，高温和高孔隙度有利于硅从裂隙带向邻近围岩扩散，从而形成较宽的蚀变。
（2）在围岩孔隙度均一分布的情况下，由于深部温度高于浅部，深部围岩蚀变比浅部宽，这与石英脉型钨矿床的蚀变特征相反。
（3）围岩孔隙度随深度递减会抵消温度对硅扩散的影响，可使深部围岩形成较窄的蚀变，围岩孔隙度随深度递减可能是形成石英脉型钨矿床蚀变宽度随深度减小的有效机制。相对于前人提出的岩浆热液过渡性流体假说，本文热液运移模拟为理解石英脉型钨矿床的蚀变形成过程提供了新的启示。
-
![]()
图 1 河北省地质环境承载能力本底评价图
Figure 1. Geological environment carrying capacity background evaluation map of Hebei Province
![]()
图 2 河北省地质环境承载能力状态评价
Figure 2. Geological environment bearing capability status evaluation map of Hebei Province
![]()
图 3 河北省地质环境承载能力分级图
Figure 3. Geological environment bearing capability classification map of Hebei Province
![]()
图 4 河北省地质环境承载能力状态评价图(试验方法)
Figure 4. Geological environment bearing capability status evaluation map(test method)of Hebei Province
表 1 河北省地质环境承载能力评价指标体系
Table 1 Evaluation index system of geological environment carrying capacity in Hebei Province
评价因子 本底评价 状态评价 崩塌、滑坡、泥石流 崩塌、滑坡、泥石流易发程度 崩塌、滑坡、泥石流风险性 构造稳定性 断裂活动性、地震动峰值加速度 — 地面塌陷 地面塌陷易发程度 损毁土地程度、地面塌陷风险性 地面沉降 地面沉降累计沉降量、地裂缝发育程度 区域地面沉降速率、沉降中心地面沉降速率 
下载: 导出CSV
表 2 地质环境承载能力评价分级标准
Table 2 Grading criteria for geological environment carrying capacity evaluation
承载能力等级 地质环境评价等级 高 较高 中 较低 低 状态评价等级 好 强 强 较强 中 较弱 中 强 较强 中 较弱 弱 差 较强 中 较弱 弱 弱
下载: 导出CSV
表 3 河北省地质环境承载能力本底评价结果统计
Table 3 Statistical table of geological environment carrying capacity background evaluation results
承载本底等级 县(市)/个 个数占比 总面积/km2 面积占比 高 32 21.3% 24800 13.1% 较高 67 44.7% 67141 35.6% 中 43 28.7% 86101 45.6% 较低 8 5.3% 10747 5.7%
下载: 导出CSV
表 4 河北省地质环境承载能力状态评价计算结果统计
Table 4 Statistical table of calculation results of geological environment carrying capacity status evaluation in Hebei Province
承载状态等级 县(市)/个 个数占比 总面积/km2 面积占比 差 37 24.7% 51393 27.2% 中 33 22.0% 58976 31.2% 好 80 53.3% 78420 41.6%
下载: 导出CSV
表 5 河北省地质环境承载能力超载县统计
Table 5 Statistics of overloaded counties of geological environment carrying capacity in Hebei Province
超载县 个数 状态差因子 备注 赤城、迁西、平山、涉县、承德县、围场、丰宁、隆化、青龙、宽城、阜平、唐县、邢台县、涿鹿 14 地质灾害风险性 宽城地质灾害和地面塌陷均超载 邢台市、张家口、阳原、唐山市、玉田、武安、邯郸县、磁县、宽城、涞水、蔚县、峰峰矿区 12 地面塌陷风险性 蠡县、高阳、衡水市、武邑、阜城、景县、枣强、肃宁、鸡泽、成安、肥乡、广平 12 地面沉降
下载: 导出CSV
表 6 河北省地质环境承载能力状态评价专家会商结果统计
Table 6 Statistical results of expert meeting results of geological environment carrying capacity status evaluation in Hebei Province
承载状态等级 县(市)/个 个数占比 总面积/km2 面积占比 差 28 18.7% 41239 21.8% 中 42 28.0% 77030 40.8% 好 80 53.3% 70520 37.4%
下载: 导出CSV
表 7 河北省地质环境承载能力统计
Table 7 Statistical table of geological environment carrying capacity evaluation in Hebei Province
承载能力等级 县(市)/个 个数占比 总面积/km2 面积占比 强 72 48% 61442 32.5% 较强 32 21% 34081 18.1% 中 25 17% 60561 32.1% 较弱 15 10% 22527 11.9% 弱 6 4% 10178 5.4%
下载: 导出CSV
表 8 河北省地质环境承载能力状态评价
Table 8 Statistical table of results of test methods for geological environment carrying capacity status of Hebei Province
县名(市) 险情变大 险情变小 小计 赤城县 1 5 6 崇礼县 7 5 12 沽源县 11 25 36 康保县 5 1 6 万全县 0 10 10 蔚县 0 6 6 宣化县 1 0 1 阳原县 0 4 4 张北县 0 8 8 张家口市区 0 11 11 涿鹿县 5 5 10 承德市区 2 0 2 丰宁县 1 1 2 宽城县 7 1 8 隆化县 9 28 37 滦平县 2 29 31 平泉县 0 72 72 围场县 2 39 41 兴隆县 0 49 49 玉田县 3 0 3 滦县 1 0 1 迁西县 10 37 47 遵化市 4 0 4 卢龙县 0 2 2 青龙县 1 1 2 阜平县 12 10 22 涞水县 6 9 15 涞源县 15 35 50 满城县 13 0 13 曲阳县 5 1 6 顺平县 0 6 6 唐县 0 2 2 易县 0 6 6 行唐县 1 3 4 井陉县 2 7 9 灵寿县 6 15 21 鹿泉市 0 12 12 平山县 13 24 37 元氏县 9 13 22 赞皇县 3 18 21 临城县 3 12 15 内丘县 1 21 22 邢台市区 3 0 3 邢台县 2 47 49 磁县 6 0 6 涉县 3 0 3 武安县 1 22 23 合计 176 602 778
下载: 导出CSV
-
段永侯, 罗元华, 柳源, 等.中国地质灾害[M].北京:中国建筑工业出版社, 1993:1-25. 李瑞敏, 王祎萍, 王轶, 等.中国主要环境地质问题[M].北京:地质出版社, 2007:1-138. 封志明, 杨艳昭, 闫慧敏, 等.百年来的资源环境承载力研究:从理论到实践[J].资源科学, 2017, 39(3):379-395. http://d.old.wanfangdata.com.cn/Periodical/zykx201703001 殷志强, 李瑞敏, 李小磊, 等.地质资源环境承载能力研究进展与发展方向[J].中国地质, 2018, 45(6):1103-1115 http://d.old.wanfangdata.com.cn/Periodical/zgdizhi201806003 王景明, 倪玉兰, 刘金峰, 等.河北省地裂缝灾害与成因分析[J].中国地质灾害防治学报, 1994, s1:98-102. 邹荟霞, 任建兰.山东省资源环境承载力综合评价与区域差异研究[J].中国环境管理干部学院学报, 2016, 26(4):30-33. http://d.old.wanfangdata.com.cn/Periodical/zghjglgbxyxb201604008 陈兆荣.安徽省资源环境承载力动态趋势与预测研究[J].重庆文理学院学报(社会科学版), 2014, 33(6):104-108. http://d.old.wanfangdata.com.cn/Periodical/yxxyxb-shkxb201406018 杨渺, 甘泉, 叶宏, 等.四川省资源环境承载力预警模型构建[J].四川环境, 2017, 36(1):144-151. doi: 10.3969/j.issn.1001-3644.2017.01.028 陈立鹏.基于PCA和SD模型国土资源环境承载力评价-以浙江省义乌县为例[D].浙江大学硕士学位论文, 2015: 9-19. 张文怡.中部六省资源环境承载力与城镇发展道路研究[D].华中科技大学硕士学位论文, 2015: 19-22. http://cdmd.cnki.com.cn/Article/CDMD-10487-1015907467.htm 余如意.资源环境承载力约束下的建设用地布局优化研究[D].四川师范大学硕士学位论文, 2016: 11-18. http://d.wanfangdata.com.cn/Periodical/jzgcjsysj2017083118 封志明, 杨艳昭, 江东, 等.自然资源资产负债表编制与资源环境承载力评价[J].生态学报, 2016, 36(22):7140-7145. http://d.old.wanfangdata.com.cn/Periodical/stxb201622024 李瑞敏, 孟晖.国土资源环境承载力评价技术要求(试行)(地质部分).北京: 国土资源部.国土资厅函[2016]1213号.2016. 李瑞敏, 孟晖, 李小磊, 等.区域地质资源环境承载能力评价技术要求(征求意见稿). 2018.
计量
- 文章访问数: 3046
- HTML全文浏览量: 296
- PDF下载量: 1946
