基于能值分析的矿业城市工业生态效率的动态变化研究——以江西为例

黄晓杏; 王锐

doi:10.13264/j.cnki.ysjskx.2020.03.012

基于能值分析的矿业城市工业生态效率的动态变化研究——以江西为例

黄晓杏^{a, b, ,},
王锐^b

a.
江西理工大学，矿业发展研究中心，江西赣州 341000
b.
江西理工大学，经济管理学院，江西赣州 341000

基金项目:

江西省高校人文社科重点基地项目 JD16082

江西省社科规划青年博士基金项目 17BJ44

详细信息

通讯作者:
黄晓杏（1982—），女，博士, 讲师，主要从事资源环境管理方面的研究。E-mail:50912896@qq.com

中图分类号: F407.32
计量
- 文章访问数: 73
- HTML全文浏览量: 28
- PDF下载量: 6
出版历程
- 收稿日期: 2020-02-19
- 发布日期: 2020-06-29
- 刊出日期: 2020-05-31

Research on the dynamic change of industrial eco-efficiency of mining cities based on emergy analysis——A case study of Jiangxi Province

HUANG Xiaoxing^{a, b, ,},
WANG Rui^b

a.
Mining Development Research Center, Jiangxi University of Science and Technology, Ganzhou 341000, China
b.
School of Economics and Management, Jiangxi University of Science and Technology, Ganzhou 341000, China

摘要

摘要: 矿业城市工业发展在促进经济增长的同时，也带来了诸多环境问题。在此现状下，合理评估矿业城市工业系统的环境影响及可持续发展水平是促进矿业城市工业转型的关键。借助能值分析可衡量和比较不同类别、不同等级能量的真实价值这一特点，对江西3个矿业城市工业系统2009—2018年的能值流量进行了核算，分析了能值效率、生态效率指标变动趋势。结果表明：各矿业城市的工业发展主要依赖本地环境资源，整体系统的对外开放性差；新余的能值自给率最高，萍乡次之，赣州最低，能值投资率的评价结果与之相反，萍乡和赣州二市的净能值产出率显著低于新余；3个矿业城市工业活动对环境的压力均超过了系统可承受范围，导致可持续发展水平均显著偏低，但生态效率却呈现缓慢上升的良好态势。未来3个矿业城市在充分利用本地资源的同时，应加大对外开放程度，大力发展绿色经济和循环经济，减少对环境的压力，促进矿业城市工业系统的可持续发展和生态效率的不断提升。
- 能值分析 /
- 矿业城市 /
- 工业系统 /
- 生态效率 /
- 江西省
Abstract: The industrial development of mining cities inevitably brings about many environmental problems while promoting economic growth. Hence the key to promoting the industrial transformation of mining cities is to reasonably assess the environmental impact and sustainable development level of the cities' industrial systems. By virtue of the characteristic that emergy analysis can be used to measure and compare the real values of different types and levels of energy, the emergy flows of industrial systems of three mining cities (Xinyu, Pingxiang and Ganzhou) in Jiangxi Province from 2009 to 2018 are calculated, and the trends of emergy efficiency and ecological efficiency indexes are analyzed. The results show that the industrial development of the three cities mainly depends on local environmental resources and their whole system are not fully open to the outside. Xinyu has the highest emergy self-sufficiency rate, followed by Pingxiang, and the lowest is Ganzhou. The evaluation of emergy investment rates of the three cities turns out to be a surprise: the net emergy output rates of Pingxiang and Ganzhou are significantly lower than that of Xinyu; the sustainable development levels of the three cities are significantly low as their industrial activities impose too much burden on the environment that are beyond the recovery capacity of the ecosystem. Nevertheless, the ecological efficiency shows a good trend of steady rise. In the future, while making full use of local resources, the three mining cities should be open wider to the outside world and vigorously develop green economy and circular economy to reduce the pressure on the environment and promote the sustainable development of industrial system and the continuous increase of ecological efficiency.
- emergy analysis /
- mining cities /
- industrial system /
- ecological efficiency /
- Jiangxi Province

HTML全文

原地浸出工艺在我国南方离子型稀土矿中推广应用已有二十多年的历史^[1],工艺的不足逐渐显现,主要体现在2个方面：一是资源问题,表现在资源综合回收率偏低,据不完全统计,半数以上企业的资源综合回收率达不到75%^[2-3]；二是环境问题,表现在地质灾害频发,如江西龙南离子型稀土矿山因原地浸矿发生规模小、数量多的滑坡等地质灾害^[4]。不同矿山具有不同的开采技术条件,需要进行专门的浸渗环节设计,但目前国内矿山生产管理粗放,常常简单套用一成不变的方法,甚至生搬硬套其他矿山的工艺和参数^[5-6],是产生上述问题的主要原因。

研究表明,上述“资源-环境”问题均主要与矿山“注液-收液”的浸渗过程密切相关。浸渗过程的资源损失主要有浸矿盲区、母液渗漏和反吸附3部分^[7],其中,前２项损失主要是注液工程、收液工程布置不合理及注液强度不当引起的。浸矿液流经矿体会对矿体的抗剪强度产生弱化作用^[8-9],尤其是当注液强度大于矿山的出渗能力时,矿山浸润面不断上升,当浸润面达到临界水位面以上时,即可在瞬间发生滑坡^[10]。而浸润面的高低形状则取决于矿山的水文地质条件、渗透性能、注液强度与出渗能力等因素。

对“浸渗环节”,很多学者基于室内试验,在机理层面研究了浸矿过程的浸出动力学^[11-13]与传质过程^[14-17],以筛选浸矿剂^[18-20]、优化浸矿工艺；部分学者基于矿山生产现场试验,研究了原地浸矿过程中矿山内浸润面、土水压力及边坡表面裂缝的变化规律^[4,21],提出矿山滑坡的防控措施^[22]。这些研究成果为浸渗环节“事中观测,事后总结”阶段提供了有力的帮助,但难以用于“事先预测（设计）”阶段。因此,池汝安等^[23]指出,从矿山宏观水动力学角度研究浸取剂溶液在矿体内部的渗流规律,可以更好地指导离子型稀土矿的原地浸出。

本文选取浑圆形裸脚式稀土矿山为研究对象,建立原地浸矿“注液-收液”轴对称简化模型,基于地下水动力学和Dupuit假定,在已知矿山地形条件及注液条件下推导矿山稳定浸润线方程,为原地浸矿浸渗环节的“事先预测（设计）”提供理论依据。

1 浑圆形裸脚式矿山原地浸矿稳定浸润线理论分析

1.1 矿山“注液-收液”轴对称简化模型

我国南方离子型稀土矿区多为低山丘陵地形,山头多为不规则的浑圆形^[24]。离子型稀土矿原地浸出工艺注液工程一般为梅花形布置的注液孔,沿山顶往山腰的一定范围内布置,注液孔直径一般为Φ=0.15～0.30 m,深度一般为见矿0.5～1.0 m。收液工程则根据矿山具体的工程地质条件确定,其中,裸脚式矿山具有较完整的天然基岩底板,且在山脚或山谷处有基岩出露,常在山脚开挖集液沟收集母液^[25],其工艺如图1所示,图2所示为江西某浑圆形稀土矿山照片。

图 1 裸脚式稀土矿山原地浸出工艺示意

Figure 1. Schematic diagram of the in-situ leaching process of a bare-foot rare earth mine

下载: 全尺寸图片幻灯片

图 2 江西某浑圆形稀土矿山

Figure 2. A rounded rare earth mine in Jiangxi Province

下载: 全尺寸图片幻灯片

对于浑圆形裸脚式矿山,沿过中心轴线截取剖面,可得到“注液-收液”轴对称简化模型,如图3所示,对该模型作如下主要假设：

图 3 浑圆形裸脚式矿山“注液-收液”轴对称简化模型

Figure 3. The “liquid injection- liquid collection” simplified axisymmetric model of a rounded barefoot type mine

下载: 全尺寸图片幻灯片

1）假设地表径流与地下径流分水岭在平面上重合处可作为中心轴线（r = 0）,过中心轴线作垂直的H轴,取水平r轴与集液沟底板山脚出露点（图3中点C,坐标为r₃,单位为m）相交,建立坐标系。

2）原地浸矿工艺要求集液沟能及时排除出渗母液,因此,可取点C为水头基准点（H_c= 0）。

3）基岩底板坡度一般较小,地形一般会有起伏,模型中简化为直线坡,按平均值β取值（单位：（°））。

4）修建集液沟时会挖除坡脚表层黏土直至矿层（见图1）,因此,山坡出渗点B（坐标为r₂,单位为m）处的坡度α（单位（°））一般较大。

5）图3中地下水的运动分为2个部分：①浸润线之上的注液区为非饱和渗流,属于土壤水分运动范畴,此处水流主要为在重力作用下的垂直向下渗流。单孔注液时,溶液从注液孔往周围非饱和矿体入渗,距离孔越远,矿体饱和度越低,湿润体的形状类似为椭球体^[26-28]；孔网注液时,由于孔与孔之间的相互影响,使得溶液在孔底以下一定范围外的非饱和矿体中的入渗变得均匀^[25]。因此,图3中注液区（注液范围为r₁,单位为m）的边界条件简化为均匀的入渗强度W（单位水平面积、单位时间内入渗补给到地下水的水量,单位为m/d）。②浸润线之下基岩底板之上的区域为饱和渗流,属于地下水动力学范畴,此处水流主要为在重力、基岩坡度及液面比降作用下的水平流动。本文浸润线方程主要由第二部分水流在地下水动力学原理下进行推导,因此,当基岩底板坡度比较平缓、液面比降不太大的情况下,可以忽略垂向分流速而采用Dupuit假定。

6）随着注液的持续进行,矿体内部逐渐从基岩底板开始往上饱和,浸润线也越来越高,直至趋于稳定。一般情况下,原地浸矿生产周期6~8个月,浸矿剂溶液在前40天左右注完^[4],而后注入上清液,此时矿体内的浸润线已接近于稳定浸润线（图3中虚线所示）。因此,可以用稳定浸润线近似代替实际最高浸润线。

7）本文不考虑矿山土体渗透系数的空间变异性及各向同性。

1.2 稳定浸润线方程的推导

将图3视为轴对称问题,即假设在环向没有水流,过水断面为以点O为圆心的圆柱面,水头线为上凸的曲线。分析图3的流线可知,水流有径向和垂向2个方向的分流速,属于轴对称二维流。若允许引进Dupuit假定,忽略垂向分流速,则水流可简化为径向的一维流动计算^[29]。

对于裸脚式稀土矿山,基岩底板一般从分水岭向集液坑方向（流线方向）呈顺坡倾斜（本文不考虑角度的正负值,取图3中α＞β＞0 °）,当β＜20 °时,渗流长度可以用径向距离来近似表示^[29],水力坡度表示为-dH/dr。此时,H（r）=h（r）+tanβ·（r₃-r）,dH/dr= dh/dr-tanβ。依据水均衡原理,利用分段法,可得任意渗流断面处的总流量公式为：

Q (r) = \{\begin{array}{l} W π r^{2} 0 \leq r \leq r_{1} \\ W π r_{1}^{2} r_{1} \leq r \leq r_{2} \end{array}

(1)

式（1）中：Q（r）为任意渗流断面处的总流量数值,单位m³/d；W为注液入渗强度数值,单位m/d；r为任意渗流断面径向坐标数值,单位m；r₁为注液范围数值,单位m；r₂为出渗点B坐标数值,单位m。

根据Dupuit微分方程

\begin{array}{l} Q (r) = 2 π r h \cdot (- k_{s} \frac{d H}{d r}) = \\ - 2 π r k_{s} [H - t a n β \cdot (r_{3} - r)] \frac{d H}{d r} \end{array}

(2)

式（2）中：β为基岩坡度数值,单位（°）。

对于注液区域（0≤r≤r₁）,可得

- 2 π r k_{s} [H - t a n β \cdot (r_{3} - r)] \frac{d H}{d r} = W π r^{2}

(3)

分离变量后,令λ = W/k_s,作0→r的不定积分

\int_{H_{0}}^{H} H d H - t a n β \cdot r_{3} \int_{H_{0}}^{H} d H + t a n β \cdot \int_{H_{0}}^{H} r d H = - \frac{λ}{2} \int_{0}^{r} r d r

(4)

由于H与r之间的函数关系未知,式（4）左边第三项无法积分,但可根据积分中值定理近似求解。令 $\int_{H_{0}}^{H} r d H = r_{m} (H - H_{0})$ ,r_m介于0→r之间,由于浸润线方程光滑,可近似取r_m= r/2。因此,得到倾斜底板注液区域（0≤r≤r₁）稳定浸润线方程

H^{2} - t a n β \cdot (H - H_{0}) \cdot (2 r_{3} - r) + \frac{λ}{2} r^{2} = H_{0}^{2} 0 \leq r \leq r_{1}

(5)

对于非注液区域（r₁≤r≤r₂）,可得

- 2 π r K_{s} [H - t a n β \cdot (r_{3} - r)] \frac{d H}{d r} = W π r_{1}^{2}

(6)

分离变量后,作r₁→r的不定积分

\int_{H_{1}}^{H} H d H - t a n β \cdot r_{3} \int_{H_{1}}^{H} d H + t a n β \cdot \int_{H_{1}}^{H} r d H = - \frac{λ r_{1}^{2}}{2} \int_{r_{1}}^{r} \frac{d r}{r}

(7)

同理,根据积分中值定理近似求解式（7）左边第三项。令 $\int_{H_{1}}^{H} r d H = r_{d} (H - H_{1})$ ,r_d介于r₁→r之间,由于浸润线方程光滑,可近似取r_d=（r₁+r）/2。因此,得到倾斜底板非注液区域（r₁≤r≤r₂）稳定浸润线方程

\begin{array}{l} H^{2} - t a n β \cdot (H - H_{1}) \cdot (2 r_{3} - r_{1} - r) + λ r_{1}^{2} l n (\frac{r}{r_{1}}) = H_{1}^{2} \\ r_{1} \leq r \leq r_{2} \end{array}

(8)

由于浸润线在渗出点B处与下游坡面相切^[30],即B点应该满足以下2个条件

H_{2} = t a n α \cdot (r_{3} - r_{2})

(9)

J_{B} = - \frac{d H}{d r} |_{r = r_{2}} = s i n α

(10)

式（9）、式（10）中：H₂为出渗点B（r = r₂）处水头数值,单位m；r₃为C点坐标数值,单位m；α为出渗点B处的坡度数值,单位（°）；J_B为出渗点B处水力梯度。

将式（10）及点（r₂, H₂）代入式（8）,并联合式（9）,可得

r_{2} = \frac{1}{2} \{r_{3} + \sqrt[]{r_{3}^{2} - \frac{2 λ r_{1}^{2}}{s i n α \cdot [t a n α - t a n β]}}\}

(11)

H_{2} = \frac{t a n α}{2} \cdot \{r_{3} - \sqrt[]{r_{3}^{2} - \frac{2 λ r_{1}^{2}}{s i n α \cdot [t a n α - t a n β]}}\}

(12)

由于H₁＞H₂＞0,定f₁=tanβ·（2r₃-r₁-r₂）,将点（r₂, H₂）代入式（11）,可得

H_{1} = \frac{1}{2} \{f_{1} + \sqrt[]{f_{1}^{2} + 4 [H_{2}^{2} + 2 λ r_{1}^{2} l n (\frac{r_{2}}{r_{1}}) - f_{1} H_{2}]}\}

(13)

由于H₀＞H₁,定f₂=tanβ·（2r₃-r₁）,将点（r₁, H₁）代入式（6）,可得

H_{0} = \frac{1}{2} [f_{2} + \sqrt[]{f_{2}^{2} + 4 (H_{1}^{2} + \frac{λ r_{1}^{2}}{2} - f_{2} H_{1})}]

(14)

综合式（5）、式（8）—式（14）可得到倾斜底板稳定浸润线方程

\{\begin{array}{l} H^{2} - t a n β \cdot (H - H_{0}) \cdot (2 r_{3} - r) + \frac{λ}{2} r^{2} = H_{0}^{2} 0 \leq r \leq r_{1} \\ H^{2} - t a n β \cdot (H - H_{1}) \cdot (2 r_{3} - r_{1} - r) + λ r_{1}^{2} l n (\frac{r}{r_{1}}) = H_{1}^{2} r_{1} \leq r \leq r_{2} \\ r_{2} = \frac{1}{2} \{r_{3} + \sqrt[]{r_{3}^{2} - \frac{2 λ r_{1}^{2}}{s i n α \cdot [t a n α - t a n β]}}\} \\ H_{2} = \frac{t a n α}{2} \cdot \{r_{3} - \sqrt[]{r_{3}^{2} - \frac{2 λ r_{1}^{2}}{s i n α \cdot [t a n α - t a n β]}}\} \\ H_{1} = \frac{1}{2} \{f_{1} + \sqrt[]{f_{1}^{2} + 4 [H_{2}^{2} + λ r_{1}^{2} l n (\frac{r_{2}}{r_{1}}) - f_{1} H_{2}]}\} \\ H_{0} = \frac{1}{2} [f_{2} + \sqrt[]{f_{2}^{2} + 4 (H_{1}^{2} + \frac{λ r_{1}^{2}}{2} - f_{2} H_{1})}] \end{array}

(15)

式（15）为一般情况下矿山稳定浸润线方程,对于基岩底板水平的情况,令β = 0°,可将式（15）简化为

\{\begin{array}{l} H^{2} + \frac{λ}{2} r^{2} = H_{0}^{2} 0 \leq r \leq r_{1} \\ H^{2} + λ r_{1}^{2} \cdot l n (\frac{r}{r_{1}}) = H_{1}^{2} r_{1} \leq r \leq r_{2} \\ r_{2} = \frac{1}{2} (r_{3} + \sqrt[]{r_{3}^{2} - \frac{2 λ r_{1}^{2}}{s i n α \cdot t a n α}}) \\ H_{2} = \frac{t a n α}{2} \cdot [r_{3} - \sqrt[]{r_{3}^{2} - \frac{2 λ r_{1}^{2}}{s i n α \cdot t a n α}}] \\ H_{1} = \sqrt[]{H_{2}^{2} + λ r_{1}^{2} \cdot l n (\frac{r_{2}}{r_{1}})} \\ H_{0} = \sqrt[]{H_{1}^{2} + \frac{λ r_{1}^{2}}{2}} \end{array}

(16)

根据式（15）,均匀稳定注液入渗（λ＞0）的条件下,注液区的浸润线是椭圆线的上半支,非注液区的浸润线是抛物线与对数曲线组成的复杂形状。令式（15）中的λ = 0、r₁=0（即不注液）,可求得H₀=H₁=H₂=tanβ·r₃,过水断面高度h₀=h₁=h₂=0,与实际情况相符。经验表明,当注液入渗量不足时,基岩面难以形成完整的浸润线。因此,如果得到完整的浸润线,应该满足条件H₀≥tan β·r₃,将该条件代入式（15）,即可求得给定注液范围下的最低注液强度或给定注液强度下的最小注液范围。

2 浸润线计算误差的参数分析

2.1 有限元模型与工况

为了进一步分析浸润线方程式（15）及式（16）的计算误差,根据矿山的地形条件和注液条件,提出4个参数：出渗坡度（α）、基岩坡度（β）、注液范围（r₁²/r₃²）和相对渗透系数（λ）,按图3简化模型建立有限元模型,将有限元计算结果作为对照值,分析上述4个参数对浸润线方程计算误差的影响规律。其中,α可根据山脚集液沟设计和施工情况确定；β应根据矿山工程勘探资料给出的地质剖面图确定；r₁²/r₃²根据矿山注液方案中矿体分布情况确定,同时避开边坡坡度陡峭处；λ根据矿山注液方案和土体渗透系数k_s确定,注液方案给出了矿山注液范围的面积A（单位：m²）和注液总强度Q（单位：m³/d）,k_s通过室内、外土体渗透性实验确定,则相对渗透系数λ=Q/（Ak_s）。

矿山建模高度为30 m,长度r₃= 50 m,计算工况见表1。

表 1 计算工况汇总

Table 1. Summary of calculation conditions

参数	取值										工况数
α /（°）	30.96	34.29	36.87	41.42	45.00	51.34	56.31	63.43^*	71.57	82.41	10
β /（°）	0	2.29	4.57	7.97^*	10.20	12.41	14.57	—	—	—	7
r₁²/r₃²	0.01	0.04	0.09	0.16	0.25^*	0.36	0.49	—	—	—	7
λ	0.04	0.08	0.12	0.16	0.20^*	0.24	0.28	0.32	0.36	—	9
注：“—”现场资料未提供数据；标有“*”的为模型参数的基本值,即：进行某一个参数的敏感性分析时,该参数按表1中数据从小到大取值,另外3个参数则取基本值。

下载: 导出CSV

| 显示表格

2.2 计算误差参数分析

典型工况（各参数均取基本值）下,浸润线对比分析如图4所示。由图4可知,方程计算浸润线基本位于数值计算浸润线之上,但两者非常接近,其计算误差沿流场方向逐渐增大,即注液段误差最小,中间段误差次之,出渗段（出渗点附近）误差最大；但在出渗点处,方程计算浸润线却位于数值计算浸润线之下。这是因为：本文浸润线方程的推导引入了Dupuit假定,该假定在地下分水岭附近及出渗点附近并不满足^[29],分水岭附近,垂向流速较大,铅垂面十分接近流面,不能假定为等水头面,Dupuit假定忽略了垂向流速,因此,公式计算值一般要比实际值偏大；出渗段等水头面往往比较弯曲,水力梯度逐渐变大,以利地下水渗出成为地表水,并产生“水跃”现象,而Dupuit假定不能考虑“水跃”现象,且本文浸润线方程在出渗点处满足浸润线与出渗边坡相切的条件,造成浸润线在该处“弯曲”严重,因此,公式计算值一般要比实际值更小。

图 4 典型工况下浸润线对比

Figure 4. Comparison chart of steady seepage lines under typical working conditions

下载: 全尺寸图片幻灯片

为进一步分析各参数对不同部位浸润线方程计算结果的影响,下面以分水岭（r = 0 m）、中间点（r₁²/r₃² = 0.5）与出渗点（r = r₂）分别分析4个参数对浸润线方程计算误差的影响规律。

浸润线计算误差随出渗坡度的变化曲线如图5所示,由图5可知：浸润线计算值在分水岭及中间点较数值解偏大,且随出渗坡度的增大而小幅增大,误差最大值为11.42%；浸润线计算值在出渗点较数值解偏小,且误差的绝对值随出渗坡度的增大呈现先增大后减小的变化趋势,当出渗坡度α=41.42°时,误差最小值为-13.51%。由于修建集液沟时会挖除坡脚表层黏土,出渗坡度α一般较大,因此,出渗坡度α对浸润线计算误差较小。

图 5 浸润线计算误差随出渗坡度的变化曲线

Figure 5. The changing curve of calculation error with the seepage gradient

下载: 全尺寸图片幻灯片

浸润线计算误差随基岩坡度的变化曲线如图6所示,由图6可知：浸润线计算值在分水岭及中间点较数值解偏大,在出渗点较数值解偏小,且误差的绝对值随基岩坡度的增大近似线性增大,该误差的产生是由于本文浸润线方程中水力梯度的计算忽略了基岩坡度的影响,而采用水平距离近似表示渗流长度,因此基岩坡度越大,误差越大。由于原地浸矿基岩坡度一般较小（β＜15°）,因此,出渗坡度对浸润线计算误差的影响较小。

图 6 浸润线计算误差随基岩坡度的变化曲线

Figure 6. The changing curve of calculation error with the gradient of bed rock

下载: 全尺寸图片幻灯片

浸润线计算误差随注液范围和相对渗透系数的变化曲线分别如图7和图8所示,由图7和图8可知：浸润线计算误差随注液范围和相对渗透系数具有类似的变化规律,即误差的绝对值随着注液范围或相对渗透系数的增大而减小。定义矿山相对注液强度W_r为矿山实际注液流量与矿山总面积下饱和注液流量之比,即W_r= λ·r₁²/r₃²,可知,图7和图8均反映的是浸润线方程计算误差随矿山相对注液强度W_r的变化规律,当相对注液强度W_r很小时,矿山基岩上甚至不能形成完整的浸润线,此时浸润线计算误差较大,这种情况只发生在原地浸矿注液初期；在注液稳定期,矿山相对注液强度W_r较大,矿山基岩上形成了较高的浸润线,浸润线方程计算误差较小。

图 7 浸润线计算误差随注液范围的变化曲线

Figure 7. The changing curve of calculation error with the range of liquid injection

下载: 全尺寸图片幻灯片

图 8 浸润线计算误差随相对渗透系数的变化曲线

Figure 8. The changing curve of calculation error with the relative permeability coefficient

下载: 全尺寸图片幻灯片

3 结论

1）选取浑圆形裸脚式稀土矿山为研究对象,建立原地浸矿“注液-收液”轴对称简化模型,基于地下水动力学和Dupuit假定,推导得到矿山稳定浸润线方程,方程以分段函数表示,其中注液区的浸润线是椭圆线的上半支,非注液区的浸润线是抛物线与对数曲线组成的复杂形状。

2）以有限元计算结果为对照值,分析出渗坡度、基岩坡度、注液范围和相对渗透系数等参数对浸润线方程计算误差的影响,发现方程计算浸润线基本位于数值计算浸润线之上,但两者非常接近,其计算误差沿流场方向逐渐增大,但对于原地浸矿而言,各参数取值对方程计算误差的影响较小。

3）原地浸矿浸润线的高低与形状是矿山“注液-收液”渗流场形成与发展的结果,同时浸润线的高低与形状又跟矿山边坡稳定性、矿体浸泡体积百分比及资源浸取率等密切相关,因此浸润线方程可为原地浸矿“注液-收液”浸渗环节的“事先预测（设计）”提供理论依据。

4）浸润线方程的求解需要确定地形条件（α、β和r₃）和注液条件（r₁和λ）,模型参数比较多,某些参数不易获取,且模型和参数均作了一定程度上的简化,计算结果应充分评估后使用。

图 1 矿业城市工业系统可更新资源能值（R）变化趋势

Fig 1. Trend of renewable resource emergy industrial system of mining cities

下载: 全尺寸图片幻灯片

图 2 矿业城市工业系统不可更新资源能值（N）变化趋势

Fig 2. Trend of non-renewable resource emergy industrial system of mining cities

下载: 全尺寸图片幻灯片

图 3 矿业城市工业系统购买资源能值（I）变化趋势

Fig 3. Trend of input resource emergy industrial system of mining cities

下载: 全尺寸图片幻灯片

图 4 矿业城市工业系统废弃物能值（W）变化趋势

Fig 4. Trend of waste emergy industrial system of mining cities

下载: 全尺寸图片幻灯片

图 5 矿业城市工业系统输出能值（EXP）变化趋势

Fig 5. Trend of export emergy industrial system of mining cities

下载: 全尺寸图片幻灯片

图 6 矿业城市工业系统能值自给率（ESR）

Fig 6. Emergy self-sufficiency rate of industrial system of mining cities

下载: 全尺寸图片幻灯片

图 7 矿业城市工业系统能值投资率（EIR）

Fig 7. Emergy investment rate of industrial system of mining cities

下载: 全尺寸图片幻灯片

图 8 矿业城市工业系统净能值产出率（EYR）

Fig 8. Net emergy yield rate of industrial system of mining cities

下载: 全尺寸图片幻灯片

图 9 矿业城市工业系统环境负载率（ELR）

Fig 9. Environmental load ratio of industrial system of mining cities

下载: 全尺寸图片幻灯片

图 10 矿业城市工业系统可持续发展指数（ESI）

Fig 10. Environment sustainable development index of industrial system of mining cities

下载: 全尺寸图片幻灯片

图 11 矿业城市工业系统生态效率指数（EEI）

Fig 11. Ecological efficiency index of industrial system of mining cities

下载: 全尺寸图片幻灯片

表 1 矿业城市工业系统能值流量评价指标

Table 1 Emergy flow evaluation indicators of industrial system of mining cities

下载: 导出CSV

表 2 矿业城市工业系统能值效率评价指标

Table 2 Emergy efficiency evaluation indicators of industrial system of mining cities

下载: 导出CSV

表 3 矿业城市工业生态系统生态效率评价指标

Table 3 Ecological efficiency evaluation indicators of industrial system of mining cities

下载: 导出CSV

表 4 矿业城市工业系统能值数据表

Table 4 Emergy data of industrial system of mining cities in Jiangxi Province

下载: 导出CSV

表 5 江西矿业城市工业系统能值流量表

Table 5 Emergy flow of industrial system of mining cities in Jiangxi Province

下载: 导出CSV

参考文献(29)

[1]	陈勇, 黄冉冉, 唐荣彬, 等.基于DPSIR和云模型的矿业城市生态安全评价[J].矿业研究与开发, 2017, 37(12):32-38. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kyyjykf201712008
[2]	尹科, 王如松, 周传斌, 等.国内外生态效率核算方法及其应用研究述评[J].生态学报, 2012(11):3595-3605. http://d.old.wanfangdata.com.cn/Periodical/stxb201211032
[3]	栾桂林.利用生态效率指标体系对青岛循环经济发展的评价与对策研究[D].青岛: 青岛大学, 2009. http://cdmd.cnki.com.cn/Article/CDMD-11906-2009128923.htm
[4]	周国梅, 彭昊, 曹凤中.循环经济和工业生态效率指标体系[J].城市环境与城市生态, 2003, 16(6):201-203. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cshjycsst200306080
[5]	NIEMINEN E, LINK M, TOBLER M, et al. EU COST Action 628: Life cycle assessment (LCA) of textile products, eco-efficiency and definition of best available technology (BAT) of textile processing[J]. Journal of Cleaner Production, 2007, 15(13/14): 1259-1270. http://cn.bing.com/academic/profile?id=88abccdad9a37047dcd09dfb6e584a40&encoded=0&v=paper_preview&mkt=zh-cn
[6]	胡彪, 付业腾.中国生态效率测度与空间差异实证—基于SBM模型与空间自相关性的分析[J].干旱区资源与环境, 2016(6):6-12. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ghqzyyhj201606002
[7]	潘兴侠, 何宜庆.工业生态效率评价及其影响因素研究—基于中国中东部省域面板数据[J].华东经济管理, 2014, 28(3):33-38. http://d.old.wanfangdata.com.cn/Periodical/hdjjgl201403007
[8]	王雪, 施晓清.基于投入产出表的北京市产业生态效率[J].生态学报, 2018(8):2660-2671. http://d.old.wanfangdata.com.cn/Periodical/stxb201808006
[9]	王震, 石磊, 刘晶茹, 等.区域工业生态效率的测算方法及应用[J].中国人口·资源与环境, 2008, 18(6):121-126. http://d.old.wanfangdata.com.cn/Periodical/zgrkzyyhj200806022
[10]	汪东, 朱坦.基于数据包络分析理论的中国区域工业生态效率研究[J].生态经济, 2011(4):24-28. http://d.old.wanfangdata.com.cn/Periodical/stjj201104004
[11]	汪克亮, 黄晴晴, 孟祥瑞.基于环境压力的矿业城市工业生态效率[J].系统工程, 2017, 35(2):36-44. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xtgc201702006
[12]	孙欣, 赵鑫, 宋马林.长江经济带生态效率评价及收敛性分析[J].华南农业大学学报(社会科学版), 2016, 15(5):1-10. http://d.old.wanfangdata.com.cn/Periodical/hnnydxxb-shkx201605001
[13]	李美娟, 陈国宏.数据包络分析法(DEA)的研究与应用[J].中国工程科学, 2003, 5(6):88 －94. http://d.old.wanfangdata.com.cn/Periodical/zggckx200306016
[14]	宋改凤, 刘艳中, 朱晓南, 等.基于能值分析的武汉市生态效率动态变化研究[J].生态经济, 2019, 35(7):103-109. http://www.cqvip.com/QK/96795X/20197/7002359767.html
[15]	蓝盛芳, 钦佩, 陆宏芳.生态经济系统能值分析[M].北京: 化学工业出版社, 2002.
[16]	朱玉林, 李明杰.湖南省农业生态系统能值演变与趋势[J].应用生态学报, 2012, 23(2):499－505. http://d.old.wanfangdata.com.cn/Periodical/yystxb201202029
[17]	聂艳, 乌云嘎, 于婧, 等.基于能值分析的武穴市耕地利用效益时空特征分析[J].长江流域资源与环境, 2015, 24(6):987-993. http://d.old.wanfangdata.com.cn/Periodical/cjlyzyyhj201506013
[18]	汤萃文, 杨莎莎, 刘丽娟, 等.基于能值理论的东祁连山森林生态系统服务功能价值评价[J].生态学杂志, 2012, 31(2):433－439. http://d.old.wanfangdata.com.cn/Periodical/stxzz201202029
[19]	李凯, 崔丽娟, 李伟, 等.基于能值代数的湿地生态系统服务评价去重复性计算[J].生态学杂志, 2016, 35(4): 1108－1116. http://d.old.wanfangdata.com.cn/Periodical/stxzz201604036
[20]	李丽锋, 惠淑荣, 宋红丽, 等.盘锦双台河口湿地生态系统服务功能能值价值评价[J].中国环境科学, 2013, 33(8): 1454－1458. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zghjkx201308016
[21]	赵晟, 李梦娜, 吴常文.舟山海域生态系统服务能值价值评估[J].生态学报, 2015, 35(3):1－11. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stxb201503010
[22]	管新建, 齐雪艳, 吴泽宁, 等.东居延海生态系统服务功能价值的能值分析[J].水土保持研究, 2012, 19(5):253－261. http://d.old.wanfangdata.com.cn/Periodical/stbcyj201205052
[23]	郭丽英, 雷敏, 刘晓琼.基于能值分析法的绿色GDP核算研究—以陕西省商洛市为例[J].自然资源学报, 2015, 30(9):1523－1533. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zrzyxb201509009
[24]	康文星, 王东, 邹金伶, 等.基于能值分析法核算的怀化市绿色GDP[J].生态学报, 2010, 30(8):2151-2158. http://d.old.wanfangdata.com.cn/Periodical/stxb201008024
[25]	孙玉峰, 郭全营.基于能值分析法的矿区循环经济系统生态效率分析[J].生态学报, 2014, 34(3):710-717. http://www.cnki.com.cn/Article/CJFDTotal-STXB201403022.htm
[26]	ODUM H T. Environment accounting: emergy and environmental decision making [M]. New York: John Wiley & Soon, 1996: 163-168.
[27]	BROWN M T, PROTANO G, ULGTATI S. Assessing geobiosphere work of generating global reserves of coal, crude oil and natural gas[J]. Ecological Modelling, 2011, 222: 879－887. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f50193147f996bd5606da94db424d51b
[28]	GIANNETTI B F, OGURA Y, BONILLA S H, et al. Emergy assessment of a coffee farm in Brazilian Cerrado considering in abroad form the environmental services, negative externalities and fair price[J]. Agricultural Systems, 2011, 104(9): 679－688. http://www.researchgate.net/publication/251575254_Emergy_assessment_of_a_coffee_farm_in_Brazilian_Cerrado_considering_in_a_broad_form_the_environmental_services_negative_externalities_and_fair_price?ev=prf_cit
[29]	吴婧慈.基于能值分析的海岛生态系统服务价值研究[D].舟山: 浙江海洋大学, 2018. http://cdmd.cnki.com.cn/Article/CDMD-10340-1018184394.htm

施引文献(3)

期刊类型引用(2)

1.	刘德峰，张世佳，张臻悦，郭文达，池汝安. 风化壳淋积型稀土矿浸润线空间分布特征及影响因素敏感性分析. 稀土. 2025(01): 1-13 . 百度学术
2.	王忠锋，冯羽生，黄伟玲. 不同结构组态对导流筒搅拌槽混合效率影响分析. 有色金属科学与工程. 2024(06): 814-821+855 . 本站查看

其他类型引用(1)

资源附件(0)

图(11) / 表(5)

计量

文章访问数: 73
HTML全文浏览量: 28
PDF下载量: 6
被引次数: 3

1 浑圆形裸脚式矿山原地浸矿稳定浸润线理论分析
1.1 矿山“注液-收液”轴对称简化模型
1.2 稳定浸润线方程的推导
2 浸润线计算误差的参数分析
2.1 有限元模型与工况
2.2 计算误差参数分析
3 结论

基于能值分析的矿业城市工业生态效率的动态变化研究——以江西为例

通讯作者: 黄晓杏（1982—），女，博士, 讲师，主要从事资源环境管理方面的研究。E-mail:50912896@qq.com

计量

出版历程