基于地质-环境对露天煤矿开采计划的优化方法及装置

1.本技术涉及露天开采技术领域，尤其是涉及到一种基于地质-环境对露天煤矿开采计划的优化方法及装置。


背景技术：

2.露天矿生产过程中不可避免地造成对生态环境的冲击，所产生的冲击量及其随时间的变化与开采计划密切相关。相关技术中，露天矿开采境界的确定均是以经济合理开采比等生产成本为基本准则，来测定最优露天矿开采境界，并未将环境成本像生产成本一样纳入到开采计划之中。


技术实现要素：

3.有鉴于此，本技术提供了一种基于地质-环境对露天煤矿开采计划的优化方法及装置，实现了尽量提高露天煤矿生产的纯经济效益和尽量降低对生态环境的冲击的双重目标，进而实现了低碳生态化优化。
4.根据本技术的一个方面，提供了一种基于地质环境对露天煤矿开采计划的优化方法，包括：
5.根据露天煤矿的最终境界，利用地质-环境最优采场序列产生算法，生成地质-环境最优采场序列{c
*
}n，其中，地质-环境最优采场序列中每一个以最终境界结尾的子序列，即为一个可能的开采计划；
6.设任一开采计划的开采寿命为n年，第t年的采场对应于地质-环境最优采场序列{c
*
}n内的第i(t)个地质-环境最优采场c
*i(t)
，开采计划的地质-环境最优采场序列的子序列为{c
*
}n＝{c
*i(1)
，c
*i(2)
，......，c
*i(n)
}，其中，i(n)＝n，n为正整数，t≤n；
7.设t＝1，在地质-环境最优采场序列内，确定原煤量大于或等于预设产能域的下限且与预设产能域的下限最为接近的采场；
8.利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv1，其中，采剥量包括以下至少一种：原煤开采量、岩石剥离量和四纪层剥离量；
9.设t＝t+1，利用环境成本内生化的开采计划评价模型，计算第t年采场i(t) 的原煤开采量；
10.若采场i(t)的原煤开采量小于预设产能域的下限且未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场i(t)+1 的原煤开采量，直至满足预设可行计划条件；
11.若采场i(t)的原煤开采量小于预设产能域的下限且采场达到最终境界，确定满足预设可行计划条件，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命，
其中，开采寿命n＝t；
12.若采场i(t)的原煤开采量大于预设产能域的上限，确定未满足预设可行计划条件，优化终止；
13.若采场i(t)的原煤开采量位于预设产能域范围内，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，此时若采场未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场i(t)+1的原煤开采量，直至满足预设可行计划条件，若采场达到最终境界，确定满足预设可行计划条件，则生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命，其中，开采寿命 n＝t；
14.设t＝n-1，在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t年的原煤开采量；
15.若原煤开采量小于预设产能域的上限，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
；
16.判断采场是否达到最终境界，若采场未达到最终境界，确定未满足预设可行计划条件，则考虑t＝t+1的下一采场，直至满足预设可行计划条件；
17.若采场达到最终境界，确定满足预设可行计划条件，生成新的可行开采计划，并根据最大年产量计算基建投资和新的可行开采计划对应的第二净现值和开采寿命，其中，开采寿命n＝t；
18.根据可行开采计划对应的第一净现值和新的可行开采计划对应的第二净现值，确定最佳开采计划；
19.设t＝n-1，返回至利用环境成本内生化的开采计划评价模型，计算第t年采场i(t)的原煤开采量的步骤，计算第t年的原煤开采量，若原煤开采量大于预设产能域的上限，则设t＝t-1；
20.判断t是否大于0，若t大于0，返回至在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t 年的原煤开采量的步骤，直至满足预设可行计划条件；
21.若t小于或等于0，确定可行开采计划为最佳开采计划，输出输送最佳开采计划。
22.可选地，根据露天煤矿的最终境界，利用地质-环境最优采场序列产生算法，生成地质-环境最优采场序列{c
*
}n的步骤，具体包括：
23.根据露天煤矿的境界区域，确定多个模柱；
24.获取每个模柱的底面标高、顶面标高和每个模柱沿竖直方向的中心线；
25.根据底面标高、顶面标高和中心线，确定露天煤矿的最终境界；
26.根据最终境界和多个模柱，在最终境界内，确定多个锥体；
27.获取每个锥体的原煤量；
28.根据每个锥体的原煤量和预设原煤增量，在多个锥体中，确定多个目标锥体；
29.获取每个目标锥体的原煤量、剥离量和地表损毁面积；
30.根据每个目标锥体的原煤量、剥离量和地表损毁面积，计算每个目标锥体的指标值，其中，指标值指的是锥体的剥采比加上单位采煤量的地表挖损面积的综合值；
31.根据指标值，确定多个目标采场；
32.根据多个目标采场，生成地质-环境最优采场序列。
33.可选地，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
的步骤，具体包括：
34.获取第t年的原煤量、岩石量和四纪层量，以及第t-1年的原煤量、岩石量和四纪层量；
35.根据第t年的原煤量、岩石量和四纪层量和第t-1年的原煤量、岩石量和四纪层量，计算第t年的原煤开采量、岩石剥离量、四纪层剥离量；
36.利用环境冲击量化模型，计算第t年的土地挖损面积、排场压占土地面积、矸石堆压占土地面积、表土场压占土地面积和地面生产设施占用土地面积；
37.根据土地挖损面积、排场压占土地面积、矸石堆压占土地面积、表土场压占土地面积和地面生产设施占用土地面积，确定第t年的土地损毁面积；
38.根据土地损毁面积，利用环境成本量化模型，计算第t年的土地生态系统损坏环境成本；
39.计算第t年的温室气体排放量；
40.根据温室气体排放量，确定第t年的温室气体排放环境成本；
41.根据土地生态系统损毁环境成本和温室气体排放环境成本，确定第t年的综合环境成本；
42.计算第t年的净现值npv
t
。
43.可选地，根据土地损毁面积，利用环境成本量化模型，计算第t年的土地生态系统损坏环境成本的步骤，具体包括：
44.利用环境成本量化模型，计算第t年的土地生态系统的生态价值，其中，土地生态系统的生态价值包括以下至少一种：生物质生产价值、土壤固碳价值、植物固碳价值、释氧价值、空气净化价值、土壤保持价值、水源涵养价值和养分循环价值；
45.获取第t年的完成复垦作业的土地面积、第(t-1)年的完成复垦并恢复生态生产力的土地面积以及第(t-1)年的土地损毁面积；
46.第t年的土地生态系统损毁环境成本为：
47.c
land,t
＝(a
t-a
r,t-1
)(z
cplant
+z
o2
+z
air
+z
soil
+z
h2o
+z
neut
)+(a
t-a
t-1
)(z
yield
+z
csoil
)+a
r,tcrec
；
48.其中，a
t
为第t年末的土地损毁面积；a
r,t-1
为第(t-1)年末已累计完成复垦并恢复生态生产力的土地面积、a
r,t
为第t年内完成复垦作业的土地面积；a
t-1
为第(t-1)年末累计土地损毁面积；z
cplant
、z
o2
、z
air
、z
soil
、z
h2o
和z
neut
分别为单位面积土地一年的植物固碳价值、释氧价值、空气净化价值、土壤保持价值、水源涵养价值和养分循环价值；z
yield
和z
csoil
分别为单位面积土地的征地价格和土壤固碳价值；c
rec
为单位面积土地的复垦和养护成本。
49.可选地，根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命的步骤，具体包括：
50.根据最大年产量q
m-max
，确定基建投资i(q
m-max
)；
51.当t＝n时，根据第t年的净现值npv
t
和基建投资，计算可行开采计划的第一净现值为：npv＝npv
n-i(q
m-max
)；
52.其中，npvn为第t年的净现值；i(q
m-max
)为露天煤矿的基建投资。
53.可选地，根据可行开采计划对应的第一净现值和新的可行开采计划对应的第二净现值，确定最佳开采计划的步骤，具体包括：
54.将第一净现值与第二净现值进行比较；
55.若第二净现值大于第一净现值，将第二净现值对应的新的可行开采计划替换掉可行开采计划，即将新的开采计划作为最佳开采计划；
56.若第二净现值小于或等于第一净现值，可行开采计划保持不变，将第一净现值对应的可行开采计划作为最佳开采计划。
57.可选地，各项环境成本包括以下至少一项：生物质生产价值成本、土壤固碳价值成本、植物固碳价值成本、释氧价值成本、空气净化价值成本、土壤保持价值成本、水源涵养价值成本、养分循环价值成本和复垦养护成本。
58.根据本技术的另一方面，提供了一种基于地质环境对露天煤矿开采计划的优化装置，包括：
59.生成模块，用于根据露天煤矿的最终境界，利用地质-环境最优采场序列产生算法，生成地质-环境最优采场序列{c
*
}n，其中，地质-环境最优采场序列中每一个以最终境界结尾的子序列，即为一个可能的开采计划；
60.设置模块，用于设任一开采计划的开采寿命为n年，第t年的采场对应于地质-环境最优采场序列{c
*
}n内的第i(t)个地质-环境最优采场c
*i(t)
，开采计划的地质-环境最优采场序列的子序列为{c
*
}n＝{c
*i(1)
，c
*i(2)
，......，c
*i(n)
}，其中， i(n)＝n，n为正整数，t≤n；
61.设置模块，还用于设t＝1；
62.确定模块，用于在地质-环境最优采场序列内，确定原煤量大于或等于预设产能域的下限且与预设产能域的下限最为接近的采场；
63.计算模块，用于利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv1，其中，采剥量包括以下至少一种：原煤开采量、岩石剥离量和四纪层剥离量；
64.设置模块，还用于设t＝t+1；
65.计算模块，还用于利用环境成本内生化的开采计划评价模型，计算第t 年采场i(t)的原煤开采量；
66.优化模块，用于若采场i(t)的原煤开采量小于预设产能域的下限且未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场i(t)+1的原煤开采量，直至满足预设可行计划条件；
67.优化模块，还用于若采场i(t)的原煤开采量小于预设产能域的下限且达到最终境界，确定满足预设可行计划条件，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命，其中，开采寿命n＝t；
68.确定模块，还用于若采场i(t)的原煤开采量大于预设产能域的上限，确定未满足预设可行计划条件，停止优化；
69.优化模块，还用于若采场i(t)的原煤开采量位于预设产能域范围内，利用环境成
本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，此时若采场未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场 i(t)+1的原煤开采量，直至满足预设可行计划条件，若采场达到最终境界，确定满足预设可行计划条件，则生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命，其中，开采寿命n＝t；
70.设置模块，还用于设t＝n-1；
71.计算模块，还用于在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t年的原煤开采量；
72.计算模块，还用于若原煤开采量小于预设产能域的上限，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
；
73.判断模块，用于判断采场是否达到最终境界；
74.优化模块，还用于若采场未达到最终境界，确定未满足预设可行计划条件，则考虑t＝t+1的下一采场，直至满足预设可行计划条件；
75.生成模块，用于若采场达到最终境界，确定满足预设可行计划条件，生成新的可行开采计划，并根据最大年产量计算基建投资和新的可行开采计划对应的第二净现值和开采寿命，其中，开采寿命n＝t；
76.确定模块，还用于根据可行开采计划对应的第一净现值和新的可行开采计划对应的第二净现值，确定最佳开采计划；
77.设置模块，还用于设t＝n-1；
78.优化模块，还用于返回至利用环境成本内生化的开采计划评价模型，计算第t年采场i(t)的原煤开采量的步骤，计算第t年的原煤开采量，若原煤开采量大于预设产能域的上限，则设t＝t-1；
79.判断模块，还用于判断t是否大于0；
80.优化模块，还用于若t大于0，返回至在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t年的原煤开采量的步骤，直至满足预设可行计划条件；
81.输出模块，用于若t小于或等于0，输出最佳开采计划。
82.根据本技术又一个方面，提供了一种存储介质，其上存储有计算机程序，所述程序被处理器执行时实现上述基于地质环境对露天煤矿开采计划的优化方法。
83.根据本技术再一个方面，提供了一种电子设备，包括存储介质、处理器及存储在存储介质上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述基于地质环境对露天煤矿开采计划的优化方法。
84.借由上述技术方案，首先在最终境界内产生一个地质-环境最优采场序列，该序列的任何一个以最终境界结尾的子序列，构成一个可能的开采计划方案，而最优计划方案必然是序列的以最终境界结尾的子序列之一。其后，基于预设的产能域内求环境成本内生化的最佳开采计划的算法，求得经济与环境的综合净现值最大的最优计划方案，使得构建的最优开采计划方案，实现了达到了尽量提高露天煤矿生产的纯经济效益和尽量降低对生态
环境的冲击的双重目标，实现了低碳生态化优化。
85.上述说明仅是本技术技术方案的概述，为了能够更清楚了解本技术的技术手段，而可依照说明书的内容予以实施，并且为了让本技术的上述和其它目的、特征和优点能够更明显易懂，以下特举本技术的具体实施方式。
附图说明
86.此处所说明的附图用来提供对本技术的进一步理解，构成本技术的一部分，本技术的示意性实施例及其说明用于解释本技术，并不构成对本技术的不当限定。在附图中：
87.图1示出了本技术实施例提供的一种基于地质环境对露天煤矿开采计划的优化方法的流程示意图；
88.图2示出了本技术实施例提供的一种在最终境界内产生地质-环境最优采场序列的剖面示意图；
89.图3示出了本技术实施例提供的一种地质-环境最优采场序列的结构示意图；
90.图4示出了本技术实施例提供的一种基于地质环境对露天煤矿开采计划的优化装置的示意框图。
具体实施方式
91.下文中将参考附图并结合实施例来详细说明本技术。需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。
92.在本实施例中提供了一种基于地质环境对露天煤矿开采计划的优化方法，如图1所示，该方法包括：
93.步骤101，根据露天煤矿的最终境界，利用地质-环境最优采场序列产生算法，生成地质-环境最优采场序列{c
*
}n；
94.步骤102，设任一开采计划的开采寿命为n年，第t年的采场对应于地质
ꢀ‑
环境最优采场序列{c
*
}n内的第i(t)个地质-环境最优采场c
*i(t)
，开采计划的地质-环境最优采场序列的子序列为{c
*
}n＝{c
*i(1)
，c
*i(2)
，......，c
*i(n)
}；
95.步骤103，设t＝1，在地质-环境最优采场序列内，确定原煤量大于或等于预设产能域的下限且与预设产能域的下限最为接近的采场；
96.步骤104，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv1；
97.步骤105，设t＝t+1，利用环境成本内生化的开采计划评价模型，计算第 t年采场i(t)的原煤开采量；
98.步骤106，若采场i(t)的原煤开采量小于预设产能域的下限且未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场i(t)+1的原煤开采量，直至满足预设可行计划条件；
99.步骤107，若采场i(t)的原煤开采量小于预设产能域的下限且采场达到最终境界，确定满足预设可行计划条件，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和
开采寿命；
100.步骤108，若采场i(t)的原煤开采量大于预设产能域的上限，确定未满足预设可行计划条件，优化终止；
101.步骤109，若采场i(t)的原煤开采量位于预设产能域范围内，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，此时若采场未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场i(t)+1的原煤开采量，直至满足预设可行计划条件，若采场达到最终境界，确定满足预设可行计划条件，则生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命；
102.步骤110，设t＝n-1，在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t年的原煤开采量；
103.步骤111，若原煤开采量小于预设产能域的上限，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
；
104.步骤112，判断采场是否达到最终境界，若采场未达到最终境界，确定未满足预设可行计划条件，则考虑t＝t+1的下一采场，直至满足预设可行计划条件；
105.步骤113，若采场达到最终境界，确定满足预设可行计划条件，生成新的可行开采计划，并根据最大年产量计算基建投资和新的可行开采计划对应的第二净现值和开采寿命；
106.步骤114，根据可行开采计划对应的第一净现值和新的可行开采计划对应的第二净现值，确定最佳开采计划；
107.步骤115，设t＝n-1，返回至利用环境成本内生化的开采计划评价模型，计算第t年采场i(t)的原煤开采量的步骤，计算第t年的原煤开采量，若原煤开采量大于预设产能域的上限，则设t＝t-1；
108.步骤116，判断t是否大于0，若t大于0，返回至在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t年的原煤开采量的步骤，直至满足预设可行计划条件；
109.步骤117，若t小于或等于0，确定可行开采计划为最佳开采计划，输出输送最佳开采计划。
110.在该实施例中，根据新建或扩建的露天煤矿的最终境界，应用地质-环境最优采场序列产生算法，确定最终境界中的地质-环境最优采场序列{c
*
}n。需要说明的是，在最终境界p内，如果一个体积为v，地表挖损面积为a，且工作帮坡角不大于给定最大工作帮坡角θ的采场c
*
，是所有体积为v，地表面积为a且工作帮坡角不大于θ的采场中原煤量最大者，那么称采场c
*
为对于v、a和θ的“地质-环境最优采场”。
111.进一步地，令{c
*
}n为最终境界p内的地质-环境最优采场序列，序列中的n个采场按总体积升序排序，最后一个采场为p。如果{c
*
}n是完全嵌套序列，且其中任意两个相邻采场之间的增量足够小，那么，使综合净现值，即纯经济盈利现值减去环境成本现值后，最大的最优开采计划，即最优解序列，必然是{c
*
}n的一个以境界结尾的子序列。
112.进一步地，在最终境界中产生了地质-环境最优采场序列{c
*
}n后，{c
*
}n中的每一
而言，子序列{c
*
}1＝{c
6*
}所代表的计划方案，是一个开采寿命只有1年的方案，即1年就开采到最终境界，此方案显然不会是最佳方案，实际上是不可行的，不用考虑。对于一个现实矿山，依据境界中的原煤储量，可以估计一个原煤年生产能力的合理范围。例如，对于一个境界中原煤储量为2亿吨的矿山，原煤年生产能力在正常情况下不会小于400万吨、大于2000 万吨。因此，通过设置产能约束条件，把原煤年生产能力约束在一个相对合理范围，在子序列的评价中，一旦某年t的原煤开采量q
m,t
超越这一范围，该子序列就被看作不可行而不予考虑，直接进入下一个子序列的评价。这样又能大大减少所评价的子序列的数量。将原煤年生产能力的约束范围称为“产能域”。基于上述思路，设计出在预设的产能域内求环境成本内生化的最佳计划方案的算法。通过这一算法，将子序列的枚举和评价同时进行，并且在整个过程中，只需在内存中保留当前子序列和当前最佳子序列，不但可大幅降低了运算量，而且内存需求量很小。在整个枚举和评价过程结束后，所保留的当前最佳子序列即为最佳开采计划方案。
118.具体地，{c
*
}n为应用“地质-环境最优采场序列产生算法”产生的地质
‑ꢀ
环境最优采场序列；i(t)为计划方案第t年末的采场在{c
*
}n中的序号。设任一开采计划方案的开采寿命为n年，其第t年末的采场对应于地质-环境最优采场序列{c
*
}n内的第i(t)个地质-环境最优采场c
*i(t)
，即构成该开采计划的地质
‑ꢀ
环境最优采场序列{c
*
}n的子序列为{c
*
}n＝{c
*i(1)
，c
*i(2)
，......，c
*i(n)
}，其中， i(n)＝n，c
*i(n)
＝c
n*
＝境界p，n为正整数，t≤n。[q
l
，qu]为预设产能域。
[0119]
步骤1：设时间t＝1，即从第一年开始。在地质-环境最优采场序列{c
*
}n中找出这样一个采场，其原煤量不小于设定的产能域下限q
l
，且与q
l
最为接近。
[0120]
如果找到了这样一个采场，它在{c
*
}n中的序号为i(1)，即正在构建的计划方案在第1年末推进到采场。应用“环境成本内生化的开采计划评价模型”计算第1年的原煤开采量q
m,1
、岩石剥离量v
y,1
、四纪层剥离量v
4,1
、年末累计土地损毁面积a1、土地生态系统损毁环境成本c
land,1
、温室气体排放量q
g,1
、温室气体排放环境成本c
g,1
、年综合环境成本c
e,1
、npv1等，进入步骤2。
[0121]
进一步地，如果找不到这样一个采场，则在预设产能域[q
l
，qu]内无可行计划，优化算法终止。
[0122]
步骤2：置时间t＝t+1，即考虑下一年。
[0123]
步骤3：置年t的采场序号i(t)＝i(t-1)+1，i(t-1)是正在构建的开采计划方案在前一年末的采场在序列{c
*
}n中的序号。
[0124]
步骤4：应用“环境成本内生化的开采计划评价模型”计算年t的原煤开采量q
m,t
。
[0125]
步骤5：(a)如果q
m,t
《q
l
，此时，分两种情况：
[0126]
(1)若i(t)＝n，即采场是序列{c
*
}n中的最后一个采场，即最终境界，说明正在构建的开采计划方案已经抵达最终境界，得到了一个可行开采计划方案，其开采寿命n＝t年；应用“环境成本内生化的开采计划评价模型”，计算第t年的岩石剥离量v
y,t
、四纪层剥离量v
4,t
、年末累计土地损毁面积a
t
、土地生态系统损毁环境成本c
land,t
、温室气体排放量q
g,t
、温室气体排放环境成本 c
g,t
、年综合环境成本c
e,t
以及净现值npv
t
等，进入步骤6。
[0127]
(2)若i(t)＜n，说明年t的原煤开采量低于设定的产能域的下限，即该开采计划不可行，此时，置采场序号i(t)＝i(t)+1，即考虑一个更大的采场，返回到步骤4。
[0128]
(b)如果q
l
≤q
m,t
≤qu，说明该采场是正在构建的开采计划方案在第t 年末的可行采场，应用“环境成本内生化的开采计划评价模型”，计算第t 年的岩石剥离量v
y,t
、四纪层剥离量v
4,t
、年末累计土地损毁面积a
t
、土地生态系统损毁环境成本c
land,t
、温室气体排放量q
g,t
、温室气体排放环境成本c
g,t
、年综合环境成本c
e,t
以及净现值npv
t
等。此时，若i(t)＝n，那么正在构建的开采计划方案已经抵达最终境界，得到了一个可行开采计划方案，其开采寿命n＝t年，进入步骤6；否则，返回到步骤2。
[0129]
(c)如果q
m,t
》qu，说明年t的原煤开采量超出了设定的产能域的上限，无可行开采计划，算法终止。
[0130]
步骤6：至此，得到了第一个可行开采计划方案，该方案具有“最低原煤年产量”，即除最后一年外，每年的原煤开采量都刚刚满足设定的产能域下限q
l
。依据该开采计划的最大年产煤量q
m-max
和基建投资函数计算投资额 i(q
m-max
)，利用npv
n-i(q
m-max
)得到这一开采计划方案的总净现值npv。置该开采计划方案为当前计划方案，并把它保存为最佳方案。
[0131]
步骤7：置时间t＝n-1，即后退一年，其中，n为当前计划方案的开采寿命。
[0132]
步骤8：从t年开始构建新的可行计划方案，新方案在1～(t－1)年与当前计划方案相同。把当前计划方案第t年的采场序号增加1，即置i(t)＝i(t)+1。
[0133]
步骤9：应用“环境成本内生化的开采计划评价模型”计算年t的原煤开采量q
m,t
。
[0134]
步骤10：(a)如果q
m,t
≤qu，说明该采场是正在构建的开采计划方案第 t年末的可行采场，将它变为当前开采计划方案第t年末的采场，即将原来的采场替换为c
*i(t)
；应用“环境成本内生化的开采计划评价模型”，计算第t 年的岩石剥离量v
y,t
、四纪层剥离量v
4,t
、年末累计土地损毁面积a
t
、土地生态系统损毁环境成本c
land,t
、温室气体排放量q
g,t
、温室气体排放环境成本c
g,t
、年综合环境成本c
e,t
以及净现值npv
t
等。此时，若i(t)＝n，那么正在构建的开采计划方案已经抵达最终境界，得到了一个新的可行开采计划方案，其开采寿命n＝t年，进入步骤15；否则，进入步骤11。
[0135]
(b)如果q
m,t
＞qu，则说明年t的原煤开采量超出了设定的产能域上限，即该开采计划不可行。此时，置t＝t-1，即沿当前计划方案往回退一年。如果这时的t＞0，返回至步骤8；否则，所有可行开采计划方案的构建和评价完毕，进入步骤16。
[0136]
步骤11：置t＝t+1，即考虑下一年。
[0137]
步骤12：置年t的采场序号i(t)＝i(t-1)+1，其中，i(t-1)是正在构建的开采计划方案前一年的采场序号。
[0138]
步骤13：应用“环境成本内生化的开采计划评价模型”，计算年t的原煤开采量q
m,t
。
[0139]
步骤14：(a)如果q
m,t
＜q
l
，分两种情况：
[0140]
(1)i(t)＝n，说明该采场是序列{c
*
}n中的最后一个采场，即最终境界，那么，正在构建的开采计划方案已经抵达最终境界，得到了一个新的可行开采计划方案，其开采寿命n＝t年；应用“环境成本内生化的开采计划评价模型”，计算第t年的岩石剥离量v
y,t
、四纪层剥离量v
4,t
、年末累计土地损毁面积a
t
、土地生态系统损毁环境成本c
land,t
、温室气体排放量q
g,t
、温室气体排放环境成本c
g,t
、年综合环境成本c
e,t
以及净现值npv
t
等，进入步骤15。
[0141]
(2)i(t)＜n，即年t的原煤开采量低于设定的产能域下限，说明该开采计划不可行。此时，置采场序号i(t)＝i(t)+1，即考虑一个更大的采场，返回至步骤13。
[0142]
(b)如果q
l
≤q
m,t
≤qu，即采场是正在构建的开采计划方案第t年末的可行采场；应
用“环境成本内生化的开采计划评价模型”，计算第t年的岩石剥离量v
y,t
、四纪层剥离量v
4,t
、年末累计土地损毁面积a
t
、土地生态系统损毁环境成本c
land,t
、温室气体排放量q
g,t
、温室气体排放环境成本c
g,t
、年综合环境成本c
e,t
以及净现值npv
t
等。此时，若i(t)＝n，那么正在构建的开采计划方案已经抵达最终境界，得到了一个新的可行开采计划方案，其开采寿命n＝t 年，进入步骤15；否则，返回至步骤11。
[0143]
(c)如果q
m,t
＞qu，说明年t的原煤开采量超出了设定的产能域上限，即该开采计划不可行，算法半途而废，显示错误信息，进入步骤16输出到此为止的最佳计划方案。
[0144]
步骤15：得到了一个新的可行开采计划方案，依据该新的可行开采计划的最大年产煤量q
m-max
和基建投资函数计算投资额i(q
m-max
)，利用npv
n-i(q
m-max
) 得到这一新的可行开采计划方案的总净现值npv，即第二净现值。如果这一新的可行开采计划方案的总净现值npv大于保存的可行开采计划方案的总净现值npv，即第一净现值，将该新的可行开采计划方案保存为最佳计划方案，即把原可行开采计划方案替换掉；否则，原可行开采计划方案不变，作为最佳开采计划方案。进一步地，置该新的可行开采计划方案为当前最佳开采计划方案，返回至步骤7。
[0145]
步骤16：输出最佳开采计划方案，算法结束。
[0146]
通过上述方式，首先在最终境界内产生一个地质-环境最优采场序列，该序列的任何一个以最终境界结尾的子序列，构成一个可能的开采计划方案，而最优计划方案必然是序列的以最终境界结尾的子序列之一。其后，基于预设的产能域内求环境成本内生化的最佳开采计划的算法，求得经济与环境的综合净现值最大的最优计划方案，使得构建的最优开采计划方案，一方面能够同时给出最佳生产能力、最佳开采顺序和最佳开采寿命，实现了开采计划三要素的整体优化；另一方面，达到了尽量提高露天煤矿生产的纯经济效益和尽量降低对生态环境的冲击的双重目标，实现了低碳生态化优化。
[0147]
在本技术实施例中，进一步地，根据露天煤矿的最终境界，利用地质-环境最优采场序列产生算法，生成地质-环境最优采场序列{c
*
}n的步骤，具体包括：根据露天煤矿的境界区域，确定多个模柱；获取每个模柱的底面标高、顶面标高和每个模柱沿竖直方向的中心线；根据底面标高、顶面标高和中心线，确定露天煤矿的最终境界；根据最终境界和多个模柱，在最终境界内，确定多个锥体；获取每个锥体的原煤量；根据每个锥体的原煤量和预设原煤增量，在多个锥体中，确定多个目标锥体；获取每个目标锥体的原煤量、剥离量和地表损毁面积；根据每个目标锥体的原煤量、剥离量和地表损毁面积，计算每个目标锥体的指标值，其中，指标值指的是锥体的剥采比和单位采煤量的地表挖损面积的综合值；根据指标值，确定多个目标采场；根据多个目标采场，生成地质-环境最优采场序列。
[0148]
在该技术方案中，最终境界p确定后，无论开采顺序如何，最终都开采到p。因此，p是序列{c
*
}n中最大的，也是最后一个地质-环境最优采场，即 c
n*
＝p。令qn表示境界内的总煤量。一个地质-环境最优采场c
*
是所有原煤量为qn，且工作帮坡角不大于给定最大工作帮坡角θ的采场中，剥采比和单位采煤量的地表挖损面积均最小的那个采场。为叙述方便，以下把“单位采煤量的地表挖损面积”简称为“吨煤开采面积”。然而，求解剥采比和吨煤开采面积均最小的采场，是一个双目标优化问题。像所有多目标优化问题一样，不同目标之间往往是不一致的。求解此类问题的一般方法是把多目标转化为单目标。因此，本技术将采场的剥采比和吨煤开采面积综合为指标λ：
[0149][0150]
其中，w为采场的剥离量；q为采场的原煤量；a为采场的地表挖损面积；f为常量系数。
[0151]
这样，求原煤量为q的、剥采比和吨煤开采面积均最小的地质-环境最优采场，就可近似地转化为求指标λ最小的采场。其中，f可以理解为吨煤开采面积相对于剥采比而言，对采场“优劣性”的重要性系数，其取值可以用单位土地损毁面积的环境成本与单位剥离量的成本，即生产加上环境成本之比值估算。
[0152]
因此，本技术中产生地质-环境最优采场序列的锥体排除算法的基本思路为：从最终境界中，排除所有原煤量等于预设的相邻采场之间的原煤量增量δq的那些部分中，指标值λ最大的那一部分，剩余部分就是境界p内原煤量为q
n-1
＝(q
n-δq)的所有采场中，λ值最小的采场，即地质-环境最优采场 c
n-1*
；同理，从采场c
n-1*
中排除所有原煤量等于δq的那些部分中λ值最大的那一部分，得到的就是下一个原煤量比c
n-1*
小δq的地质-环境最优采场，即c
n-2*
；如此进行下去，直到剩余的部分的原煤量等于或小于δq，该剩余部分即为{c
*
}n中的最小采场c
1*
。由于每个采场是从一个更大的采场中排除一部分产生的，所以如此得到的采场序列一定是一个完全嵌套序列。
[0153]
具体地，根据露天煤矿的境界区域，确定多个模柱。获取多个模柱中，每个模柱的底面标高、顶面标高和沿竖直方向的中心线。将境界区域内的每一模柱的底面标高设置为最终境界p的边帮在同一模柱的中心线处的标高。将境界区域内的每一模柱的顶面标高设置为同一模柱中心线处的地表标高。把境界区域外的每一模柱的底面和顶面标高都设置为同一模柱中心线处的地表标高。如此设置后，模柱的顶面和底面之间的柱体的集合就构成了以柱状模型形式表示的最终境界p，亦即序列{c
*
}n中的最大采场。置当前采场为最终境界p。
[0154]
进一步地，根据最终境界和多个模柱的底面标高、顶面标高和中心线，确定最终境界中的多个锥体，需要说明的是，锥体可以近似看做是最终境界内可能开采的采场的形状。根据多个锥体的原煤量和预设原煤增量，将所有锥体都“扫描”一遍，得到一个锥体数组，对锥体数组按指标值从大到小排序。然后，逐次计算锥体数组中前nc个锥体的联合体的原煤量，需要注意的是，联合体中锥体之间的重叠部分只计一次，直到nc＝u时，其中，u为正整数，即锥体数组中前u个锥体的联合体的原煤量不大于1.1δq，且与δq最为接近。将锥体数组中前u个锥体从当前采场中排除，其中，排除某个锥体，就是把当前采场中模柱底面标高小于该锥体的锥壳标高的所有模柱的底面标高，提升到相应模柱处该锥体的锥壳标高，需要注意的是，如果锥壳标高大于地表标高，就提升到地表标高。这些锥体被排除后，得到一个新的、煤量缩小了约δq的地质-环境最优采场，保存新采场。如果新采场的原煤量小于或等于δq，说明已经完成了所有采场的产生，把产生的所有采场按总体积从小到大排序排列，就得到了地质-环境最优采场序列。
[0155]
在实际应用场景中，步骤1：如图2所示，为最终境界内产生地质-环境最优采场序列的剖面图。具体地，将境界p区域内的每一模柱的底面标高设置为境界p的边帮，如图2中“境界p”所指的粗虚线位置处，在同一模柱的中心线处的标高，或是设置为境界坑底在同一模柱的中心线处的标高；把境界p区域内的每一模柱的顶面标高设置为同一模柱中心线处
的地表标高；把境界p区域外的每一模柱的底面和顶面标高都设置为同一模柱中心线处的地表标高。如此设置后，模柱的顶面和底面之间的柱体的集合就构成了以柱状模型形式表示的最终境界p，亦即序列{c
*
}n中的最大采场。置当前采场为最终境界p。
[0156]
步骤2：置模柱序数k＝1。
[0157]
步骤3：如果模柱k的底面标高与顶面标高相等，说明该模柱在当前采场的范围之外，进入步骤7；否则，继续步骤4。
[0158]
步骤4：如图2中的锥a所示，将锥体顶点置于模柱k中心线上标高为 z＝z
l,k
+δh处，其中，z
l,k
是模柱k的底面标高。
[0159]
步骤5：计算锥体，即位于锥壳以下和当前采场以上的部分的剥离量、原煤量和地表面积，进而计算该锥体的指标λ值。如果锥体的原煤量不大于预设的相邻采场之间的原煤量增量δq，把锥体的这些参数值与锥体顶点坐标一起记录于一个“锥体数组”中，进入步骤6；否则，该锥体不予考虑，转至步骤7。
[0160]
步骤6：置z＝z+δh，即把锥体顶点沿模柱k的中心线上移一个步长δh，如图2中锥a之上的锥体。如果z≤z
u,k
+δh，其中，z
u,k
为地表在模柱k处的标高，返回至步骤5；否则，进入步骤7；
[0161]
步骤7：置模柱序数k＝k+1，即把锥体移动到下一个模柱。如果k≤nc，其中，nc为煤层柱状模型的模柱总数，返回至步骤3；否则，说明所有模柱已经被锥体“扫描”一遍，得到了一个锥体数组，进入步骤8。
[0162]
步骤8：对锥体数组按λ值从大到小排序。然后，对于nc＝1，2，......，u，其中，u为正整数，逐次计算锥体数组中前nc个锥体的联合体的原煤量，需要注意的是，联合体中锥体之间的重叠部分只计一次，直到nc＝u时，即锥体数组中前u个锥体的联合体的原煤量不大于1.1δq，且与δq最为接近。
[0163]
步骤9：将锥体数组中前u个锥体从当前采场中排除，其中，排除某个锥体，就是把当前采场中模柱底面标高小于该锥体的锥壳标高的所有模柱的底面标高，提升到相应模柱处该锥体的锥壳标高，需要注意的是，如果锥壳标高大于地表标高，就提升到地表标高。这些锥体被排除后，得到一个新的、煤量缩小了约δq的地质-环境最优采场，保存新采场。
[0164]
步骤10：如果新采场的原煤量大于δq，置当前采场为新采场，返回到步骤2；否则，已经完成了所有采场的产生，进入步骤11。
[0165]
步骤11：把产生的所有采场按总体积从小到大排序排列，就得到了地质
‑ꢀ
环境最优采场序列。输出结果，算法结束。
[0166]
需要说明的是，算法中锥体沿同一模柱的移动步长δh一般可取台阶高度 h。一般而言，δh越小，产生的采场越优，但运算时间越长。因此在软件设计中，可以由用户输入一个优化级别参数yj(yj为整数)，算法中的δh取h/yj。
[0167]
在本技术实施例中，进一步地，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
的步骤，具体包括：获取第t年的原煤量、岩石量和四纪层量，以及第t-1年的原煤量、岩石量和四纪层量；根据第t年的原煤量、岩石量和四纪层量和第t-1年的原煤量、岩石量和四纪层量，计算第t年的原煤开采量、岩石剥离量、四纪层剥离量；利用环境冲击量化模型，计算第t年的土地挖损面积、排场压占土地面积、矸石堆压占土地面积、表土场压占土地面积和地面生产设施占用
土地面积；根据土地挖损面积、排场压占土地面积、矸石堆压占土地面积、表土场压占土地面积和地面生产设施占用土地面积，确定第t 年的土地损毁面积；根据土地损毁面积，利用环境成本量化模型，计算第t 年的土地生态系统损坏环境成本；计算第t年的温室气体排放量；根据温室气体排放量，确定第t年的温室气体排放环境成本；根据土地生态系统损毁环境成本和温室气体排放环境成本，确定第t年的综合环境成本；计算第t年的净现值npv
t
。
[0168]
在该技术方案中，令q
*m,i
、v
*y,i
、v
*4,i
、a
*c,i
分别为序列{c
*
}n中的第i 个地质-环境最优采场c
i*
的原煤量、岩石量、四纪层量和地表挖损面积，其中， i＝1，2，......，n，n为正整数，q
*m,i
是考虑了原煤回采率后的量，即采出原煤量，损失的原煤计入v
*y,i
。为表述方便，以下把“地质-环境最优采场”简称为“采场”。
[0169]
不失一般性，假如从第1年开始，已经完成前(t-1)年的相关计算，正在进行第t年的计算。
[0170]
第t年末的采场为，是从第(t-1)年末的采场经过一年的开采扩展而来。因此，该开采计划方案在第t年的原煤开采量q
m,t
、岩石剥离量v
y,t
、四纪层剥离量v
4,t
为：
[0171][0172][0173][0174]
其中，指的是第t年的原煤量；指的是第(t-1)年的原煤量；指的是第t年的岩石量；指的是第(t-1)年的岩石量；指的是第t年的四纪层量；指的是第(t-1)年的四纪层量。
[0175]
进一步地，应用“环境冲击量化模型”，计算第t年末采场累计土地挖损面积a
c,t
、排土场累计压占土地面积a
p,t
、矸石堆累计压占土地面积a
g,t
、表土场累计压占土地面积a
b,t
以及地面生产设施累计占用土地面积a
s,t
。则第t年末累计土地损毁面积a
t
为：a
t
＝a
c,t
+a
p,t
+a
g,t
+a
b,t
+a
s,t
。
[0176]
进一步地，土地生态系统损毁的环境成本包括土地生态系统的生态价值损失和复垦与养护成本。土地生态系统的生态价值包括：生物质生产价值、土壤固碳价值、植物固碳价值、释氧价值、空气净化价值、土壤保持价值、水源涵养价值和养分循环价值。如果土地的生物质生产价值以征地费用计算并假设随用随征，那么该项价值损失按第t年新增的土地损毁面积计算。假设土壤碳库中的碳在土地被损毁的年内释放，土壤固碳价值的损失也按第t年新增的土地损毁面积计算。植物固碳价值、释氧价值、空气净化价值、土壤保持价值、水源涵养价值和养分循环价值的损失是持续性的，土地一旦被损毁，这些价值的损失从土地被毁开始一直持续到土地被复垦并恢复其生态生产力，所以第t年这六项生态价值损失需按这年的损毁土地保有量计算；第t年的损毁土地保有量等于第t年末累计土地损毁面积减去前一年末已完成复垦并恢复生态生产力的累计土地面积。因此，开采计划方案中，第t年的土地生态系统损毁环境成本c
land,t
为：
[0177]cland,t
＝(a
t-a
r,t-1
)(z
cplant
+z
o2
+z
air
+z
soil
+z
h2o
+z
neut
)+(a
t-a
t-1
)(z
yield
+z
csoil
)+a
r,tcrec
；
[0178]
其中，a
t
为第t年末的土地损毁面积；a
r,t-1
为第(t-1)年末已累计完成复垦并恢复生态生产力的土地面积、a
r,t
为第t年内完成复垦作业的土地面积；a
t-1
为第(t-1)年末累计土地损毁面积；z
cplant
、z
o2
、z
air
、z
soil
、z
h2o
和z
neut
分别为单位面积土地一年的植物固碳价值、释氧价值、空气净化价值、土壤保持价值、水源涵养价值和养分循环价值；z
yield
和z
csoil
分别为单位面积土地的征地价格和土壤固碳价值，需要说明的是，这些价值均应用“环境成本量化模型”计算；c
rec
为单位面积土地的复垦和养护成本。
[0179]
进一步地，开采计划方案第t年的温室气体排放量q
g,t
为：
[0180]qg,t
＝q
m,t
(em+e
x
)+v
y,tey
+v
4,t
e4；
[0181]
其中，em、e
x
、ey和e4分别为原煤开采、煤炭洗选、岩石剥离和四纪层剥离的单位作业量的温室气体排放量，以co2当量计，这些温室气体排放量均应用“环境冲击量化模型”计算。需要注意的是，温室气体排放量以co2当量计
[0182]
进一步地，开采计划方案中第t年的温室气体排放环境成本c
g，t
为：
[0183][0184]
其中，φ为c到co2的转换系数，转换系数的取值可以为3.6667；cc为碳成本，即捕捉和贮存大气中碳的单位成本。
[0185]
进一步地，第t年的综合环境成本c
e,t
为：c
e,t
＝c
land,t
+c
g,t
。
[0186]
进一步地，假设矿山销售的产品全部为经过洗选后的精煤，且销售收入和各项成本均发生在年末。该计划方案在第t年末的累计综合净现值npv
t
为：
[0187][0188]
其中，npv
t-1
为开采计划方案在第(t-1)年末的累计综合净现值；r
x
为洗选厂的回收率；p
t
为第t年的精煤售价；c
m,t
、c
x,t
、c
y,t
和c
4,t
分别为第t年原煤开采、洗选、岩石剥离和四纪层剥离的单位生产成本；cf为矿山企业的年不变成本；d为折现率。需要注意的是，累计综合净现值均为不计基建投资的净现值。
[0189]
需要说明的是，当t＝0时的初始条件为：i(0)＝0，q
*m,0
＝0，v
*y,0
＝0，v
*4,0
＝0， a0＝基建期损毁土地面积，a
r,0
＝0，npv0＝0。
[0190]
在本技术中，进一步地，根据土地损毁面积，利用环境成本量化模型，计算第t年的土地生态系统损坏环境成本的步骤，具体包括：利用环境成本量化模型，计算第t年的土地生态系统的生态价值，其中，土地生态系统的生态价值包括以下至少一种：生物质生产价值、土壤固碳价值、植物固碳价值、释氧价值、空气净化价值、土壤保持价值、水源涵养价值和养分循环价值；获取第t年的完成复垦作业的土地面积、第(t-1)年的完成复垦并恢复生态生产力的土地面积以及所述第(t-1)年的土地损毁面积；第t年的土地生态系统损毁环境成本为：
[0191]cland,t
＝(a
t-a
r,t-1
)(z
cplant
+z
o2
+z
air
+z
soil
+z
h2o
+z
neut
)+(a
t-a
t-1
)(z
yield
+z
csoil
)+a
r,tcrec
；
[0192]
在本技术实施例中，进一步地，根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命的步骤，具体包括：根据最大年产量q
m-max
，确定基建投资i(q
m-max
)；当t＝n时，根据第t年的净现值npv
t
和基建投资，计算可行开采计划的第一净现值
为：npv＝npv
n-i(q
m-max
)；其中，npvn为第t年的净现值；i(q
m-max
)为露天煤矿的基建投资。
[0193]
在该技术方案中，根据可行开采计划方案的最高年采煤量，确定基建投资，需要说明的是，该基建投资是开采计划方案的最高年采煤量q
m-max
的函数，一般为近似线性函数。进一步地，从第一年开始逐年计算，直到最后一年，即t＝n时，得到npvn。则当前可行开采计划方案的总综合净现值，即第一净现值为：npv＝npv
n-i(q
m-max
)；其中，i(q
m-max
)为矿山的基建投资。
[0194]
进一步地，作为上述基于地质环境对露天煤矿开采计划的优化方法的具体实现，本技术实施例提供了一种基于地质环境对露天煤矿开采计划的优化装置，如图4所示，该基于地质环境对露天煤矿开采计划的优化装置400包括：
[0195]
生成模块401，用于根据露天煤矿的最终境界，利用地质-环境最优采场序列产生算法，生成地质-环境最优采场序列{c
*
}n，其中，地质-环境最优采场序列中每一个以最终境界结尾的子序列，即为一个可能的开采计划；
[0196]
设置模块402，用于设任一开采计划的开采寿命为n年，第t年的采场对应于地质-环境最优采场序列{c
*
}n内的第i(t)个地质-环境最优采场c
*i(t)
，开采计划的地质-环境最优采场序列的子序列为{c
*
}n＝{c
*i(1)
，c
*i(2)
，......，c
*i(n)
}，其中，i(n)＝n，n为正整数，t≤n；
[0197]
设置模块402，还用于设t＝1；
[0198]
确定模块403，用于在地质-环境最优采场序列内，确定原煤量大于或等于预设产能域的下限且与预设产能域的下限最为接近的采场；
[0199]
计算模块404，用于利用环境成本内生化的开采计划评价模型，计算第t 年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv1，其中，采剥量包括以下至少一种：原煤开采量、岩石剥离量和四纪层剥离量；
[0200]
设置模块402，还用于设t＝t+1；
[0201]
计算模块404，还用于利用环境成本内生化的开采计划评价模型，计算第 t年采场i(t)的原煤开采量；
[0202]
优化模块405，用于若采场i(t)的原煤开采量小于预设产能域的下限且未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场 i(t)+1，计算采场i(t)+1的原煤开采量，直至满足预设可行计划条件；
[0203]
优化模块405，还用于若采场i(t)的原煤开采量小于预设产能域的下限且达到最终境界，确定满足预设可行计划条件，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命，其中，开采寿命n＝t；
[0204]
确定模块403，还用于若采场i(t)的原煤开采量大于预设产能域的上限，确定未满足预设可行计划条件，停止优化；
[0205]
优化模块405，还用于若采场i(t)的原煤开采量位于预设产能域范围内，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
，此时若采场未达到最终境界，确定未满足预设可行计划条件，则考虑下一个更大的采场i(t)+1，计算采场i(t)+1的原煤开采量，直至满足预设可行计划条件，
若采场达到最终境界，确定满足预设可行计划条件，则生成满足预设产能域的下限的可行开采计划，并根据最大年产能计算基建投资、可行开采计划对应的第一净现值和开采寿命，其中，开采寿命n＝t；
[0206]
设置模块402，还用于设t＝n-1；
[0207]
计算模块404，还用于在可行开采计划的基础上考虑下一个更大的采场 i(t)+1，利用环境成本内生化的开采计划评价模型，计算第t年的原煤开采量；
[0208]
计算模块404，还用于若原煤开采量小于预设产能域的上限，利用环境成本内生化的开采计划评价模型，计算第t年的采剥量、累计土地损毁面积、碳排放量、各项环境成本和净现值npv
t
；
[0209]
判断模块406，用于判断采场是否达到最终境界；
[0210]
优化模块405，还用于若采场未达到最终境界，确定未满足预设可行计划条件，则考虑t＝t+1的下一采场，直至满足预设可行计划条件；
[0211]
生成模块401，用于若采场达到最终境界，确定满足预设可行计划条件，生成新的可行开采计划，并根据最大年产量计算基建投资和新的可行开采计划对应的第二净现值和开采寿命，其中，开采寿命n＝t；
[0212]
确定模块403，还用于根据可行开采计划对应的第一净现值和新的可行开采计划对应的第二净现值，确定最佳开采计划；
[0213]
设置模块402，还用于设t＝n-1；
[0214]
优化模块405，还用于返回至利用环境成本内生化的开采计划评价模型，计算第t年采场i(t)的原煤开采量的步骤，计算第t年的原煤开采量，若原煤开采量大于预设产能域的上限，则设t＝t-1；
[0215]
判断模块406，还用于判断t是否大于0；
[0216]
优化模块405，还用于若t大于0，返回至在可行开采计划的基础上考虑下一个更大的采场i(t)+1，利用环境成本内生化的开采计划评价模型，计算第 t年的原煤开采量的步骤，直至满足预设可行计划条件；
[0217]
输出模块407，用于若t小于或等于0，输出最佳开采计划。
[0218]
需要说明的是，本技术实施例提供的一种基于地质环境对露天煤矿开采计划的优化装置所涉及各功能模块的其他相应描述，可以参考图1的对应描述，在此不再赘述。
[0219]
基于上述如图1所示方法，相应的，本技术实施例还提供了一种存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述如图1所示的基于地质环境对露天煤矿开采计划的优化方法。
[0220]
基于这样的理解，本技术的技术方案可以以软件产品或硬件产品或软硬件相结合的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是cd-rom，u盘，移动硬盘等)中，包括若干指令用以使得一台电子设备(可以是个人计算机，服务器，或者网络设备等)执行本技术各个实施场景的方法。
[0221]
基于上述如图1所示的方法，以及图4所示的基于地质环境对露天煤矿开采计划的优化装置实施例，为了实现上述目的，本技术实施例还提供了一种电子设备，具体可以为个人计算机、服务器、网络设备等，该电子设备包括存储介质和处理器；存储介质，用于存储计算机程序；处理器，用于执行计算机程序以实现上述如图1所示的基于地质环境对露天煤矿
开采计划的优化方法。
[0222]
可选地，该电子设备还可以包括用户接口、网络接口、摄像头、射频(radiofrequency，rf)电路，传感器、音频电路、wi-fi模块等等。用户接口可以包括显示屏(display)、输入单元比如键盘(keyboard)等，可选用户接口还可以包括usb接口、读卡器接口等。网络接口可选的可以包括标准的有线接口、无线接口(如蓝牙接口、wi-fi接口)等。
[0223]
本领域技术人员可以理解，本实施例提供的一种电子设备结构并不构成对该电子设备的限定，可以包括更多或更少的部件，或者组合某些部件，或者不同的部件布置。
[0224]
存储介质中还可以包括操作系统、网络通信模块。操作系统是管理和保存电子设备硬件和软件资源的程序，支持信息处理程序以及其它软件和/或程序的运行。网络通信模块用于实现存储介质内部各控件之间的通信，以及与该实体设备中其它硬件和软件之间通信。
[0225]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到本技术可以借助软件加必要的通用硬件平台的方式来实现，也可以通过硬件实现。
[0226]
本领域技术人员可以理解附图只是一个优选实施场景的示意图，附图中的单元或流程并不一定是实施本技术所必须的。本领域技术人员可以理解实施场景中的装置中的单元可以按照实施场景描述进行分布于实施场景的装置中，也可以进行相应变化位于不同于本实施场景的一个或多个装置中。上述实施场景的单元可以合并为一个单元，也可以进一步拆分成多个子单元。
[0227]
上述本技术序号仅仅为了描述，不代表实施场景的优劣。以上公开的仅为本技术的几个具体实施场景，但是，本技术并非局限于此，任何本领域的技术人员能思之的变化都应落入本技术的保护范围。