一种确定选煤厂品种煤装车计划的动态价值寻优方法与流程

1.本发明涉及选煤厂配煤技术领域，具体涉及一种确定选煤厂品种煤装车计划的动态价值寻优方法，适用于大中型选煤厂品煤种装车计划的制定。


背景技术：

2.选煤厂通过洗选将出矿原煤分成不同品质的煤，装入筒仓中，然后根据市场需求将筒仓中的煤按照要求的比例掺配成品煤种和电煤销售，其中品种煤的价值更高，煤质要求也更高，提高品煤种的销量可显著提高选煤厂的经济效益。由于品煤种的价格随市场波动，品种煤的定价一般是按照达标和不达标考虑，这就导致在掺配不同品种煤时，品种煤的有效价值也是波动的，给产品煤的装车计划制定带来了困难，大多选煤厂在制定品种煤装车计划时，仅凭经验，没有实现品种煤的价值效益最大化。


技术实现要素：

3.为了克服上述现有技术存在的缺点，本发明提供了一种确定选煤厂品种煤装车计划的动态价值寻优方法，适用于大中型选煤厂的当日、当月、季度或者年度装车计划的制定。
4.为达到上述目的，本发明采用如下技术方案：
5.一种确定选煤厂品种煤装车计划的动态价值寻优方法，包括以下步骤：
6.第一步：取得需要装车品种煤的基本信息以及各筒仓的煤质和煤量信息；
7.第二步：判别不同混煤的品种煤价值；
8.第三步：求解所有穷举中判定为品种煤的混煤的最大装车趟数；
9.第四步：确定最优装车品种煤及对应装车趟数、各筒仓需要的煤量；
10.第五步：装完第四步推荐的煤种煤后，在剩余筒仓煤量下，重复第三步到第四步，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计不同品种煤的装车趟数，得到总的品种煤装车趟数。
11.本发明进一步的改进在于，第一步中，具体如下：
12.(1)品种煤相关信息
13.品种煤k，k表示品种煤的序号，选煤厂通常有两种品种煤，即k最大为2；
14.mtmbk，表示第k种品种煤的全水分控制范围，％；mtmink表示全水分控制下限，mtmaxk表示全水分控制上限；
15.admbk，表示第k种品种煤的干基灰分控制范围，％；admink表示干基灰分控制下限，admaxk表示干基灰分控制上限；
16.qmbk，表示第k种品种煤的收到基水分控制范围，kcal/kg；qmink表示收到基低位发热量控制下限，qmaxk表示收到基低位发热量控制控制上限；
17.根据用户需求对mtmbk、admbk和qmbk指标进行单一指标或者组合指标进行限定；
18.pk，表示第k品种煤煤质达标后的原煤价格，元/吨；
19.dwk，表示第k种品种煤每趟车的装车煤量，吨；
20.(2)各筒仓信息：
21.筒仓i：表示第i个筒仓；
22.mtci，表示第i个筒仓煤样的全水分，％；
23.adci，表示第i个筒仓煤样的干基灰分，％；
24.qci，表示第i个筒仓煤样的收到基低位发热量，kcal/kg；
25.煤量mli，表示第i个筒仓的初始煤样重量，t；
26.(3)装车趟数：
27.zzcts，计划的总装车趟数，趟；
28.n，根据筒仓煤量推荐的品种煤装车趟数，趟。
29.本发明进一步的改进在于，第二步中，具体方法如下
30.201)穷举筒仓在不同比例掺配时的混煤煤质，具体见式(1)到式(5)，r1、r2、r3分别表示三个仓的配比，％，mth、aarh、qh表示掺配后混煤的全水分、收到基灰分、收到基低位发热量，adh表示掺配后混煤的干基灰分，j表示三个筒仓的掺配比例按照1％步长变化计算的次数：
31.mthj＝(r1j
×
mtc1+r2j
×
mtc2+r3j
×
mtc3)/100
ꢀꢀꢀ
(1)
32.aarhj＝(r1j
×
aarc1+r2j
×
aarc2+r3j
×
aarc3)/100
ꢀꢀ
(2)
33.qhj＝(r1j
×
qc1+r2j
×
qc2+r3j
×
qc3)/100
ꢀꢀꢀ
(3)
34.r1j+r2j+r3j＝100 (4)
35.adhj＝aarhj
×
100/(100-mthj)
ꢀꢀꢀ
(5)
36.202)判断每一个掺配比例下的混煤煤是“品种煤1”、“品种煤2”，还是“非品种煤”。
37.本发明进一步的改进在于，第二步的202)中，a)如果mthj、adhj、qhj的指标满足品种煤1，则标记为“品种煤1”，该混煤的单价pd1j＝p1/qhj；
38.b)如果mthj、adhj、qhj的指标满足品种煤2，则标记为“品种煤2”，该混煤的单价pd2j＝p2/qhj；
39.c)如果混煤的煤质指标既满足品种煤1又满足品种煤2的要求，则对比每大卡热值的价值判断该煤种属于哪种煤，具体措施如下：
40.如果该煤种的pd1j≥pd2j，则最终标记该煤种为“品种煤1”，否则标记为“品种煤2”，
41.d)煤质既不满足品种煤1也不满足品种煤2的要求，则该混煤标记为：“非品种煤”。
42.本发明进一步的改进在于，第三步中，求解方法如下：
43.301)求解第j种掺配比例下的最大装车趟数njmax，j从1开始；
44.a)假定在该种掺配方案下的最大装车趟数njmax＝1，
45.b)在njmax装车趟数时需要的每个筒仓的煤量mlnij，t，具体见式(6)
46.mlnij＝rij
×
dwk
×
njmax
ꢀꢀꢀꢀ
(6)
47.mlnij：表示第j种掺配比例下，配置第j种混煤时需要的第i个筒仓的煤量，t；
48.dwk：表示第j种掺配比例下，当混煤判定为品种煤1时，dwk对应品种煤1的装车吨位，当判定为品种煤2时，dwk对应品种煤2的装车吨位；
49.c)如果mlnij≤vgi，最大装车趟数增加1趟，则njmax＝njmax+1，再次从a)到b)循
环找出最大njmax，即各个筒仓在满足煤量的条件下，最多能装多少趟该品种煤；
50.d)确定在各筒仓装完njmax趟品种煤时剩余煤量mlsij，和品种煤的总价值pjmax，具体计算分别见式(7)和式(8)
51.mlsij＝mli-mlnij
ꢀꢀꢀꢀ
(7)
52.pjmax＝njmax
×
dwk
×
pk/10000
ꢀꢀꢀ
(8)
53.302)重复步骤301)，直到求解得出所有穷举中判定为“品种煤1”、“品种煤2”的最大装车趟数njmax和总价值pjmax。
54.本发明进一步的改进在于，第四步中，求解方法如下：
55.情况1：最大装车趟数只对应一个掺配比例，则第三步中装车趟数最多的确定为最优装车品种煤；
56.情况2：当装车趟数最多的对应多个配比时，确定总价值最高的为最优装车品种煤；
57.情况3：当装车趟数最多的对应多个配比时，且总价值最高的也有多种配比时，确定混煤发热量qhj最低的为最优装车品种煤；
58.情况4：当前最大装车趟数均为0，则表示本次无法配出品种煤；
59.满足情况1或者情况2或者情况3时，标记该计算工况为最优装车工况，则装完本次推荐装车计划后各筒仓需要的煤量mlnif，t；计算参见式(9)：
60.mlnif＝rif
×
dwf
×
nmaxf
ꢀꢀꢀꢀ
(9)
61.rif：第f次求解当前煤量下的最优装车品种煤及对应最大装车趟数时，最优装车品种煤对应的各筒仓的掺配比例，％；
62.dwf：第f次求解当前煤量下的最优装车品种煤及对应最大装车趟数时，最优装车品种煤对应的每趟车的装车煤量，t；
63.nmaxf：第f次求解当前煤量下的最优装车品种煤及对应最大装车趟数时，最优装车品种煤对应的最大装车趟数，t；
64.装完本次推荐最优装车品种煤，此时各筒仓的剩余煤量mlsif，t；计算参见式(10)：
65.mlsif＝mli
‑ꢀ
mlnif
ꢀꢀꢀꢀ
(10)。
66.本发明进一步的改进在于，第五步中，求解下一个推荐装车品种煤及趟数如下：
67.当mls1f+mls2f+mls3f》0，nmaxf《zzcts，在当前筒仓煤量下即mlif＝mlsif，重复第三步到第四步的步骤401)，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计品种煤1和品种煤2的装车趟数，得到总的品种煤装车趟数nf。
68.本发明进一步的改进在于，第五步中，求解最终推荐的品种煤装车趟数n如下：
69.如果推荐的品种煤的总装车趟数nf≤zzcts，则品种煤的推荐装车趟数为n＝nf，否则品种煤的推荐装车趟数为n＝zzcts。
70.本发明至少具有如下有益的技术效果：
71.1、通过动态配煤过程中不同配比混煤的价值最优法确定品种煤的具体分类，实现了品种煤系列和价值的精准定位。
72.2、根据品种煤的煤质要求和销售价格以及筒仓存煤情况，构建了装车计划品种煤
的精准配煤模型，实现装车计划品种煤装车的精细化。
73.3、实现品种煤按梯度价格销售、在配煤过程中价值动态变化下品种煤产量和销售效益最大化，提高选煤厂经济效益。
附图说明
74.图1为不同混煤的品种煤类别及价值的判定方法流程图。
75.图2为不同品种煤的最大装车趟数确定方法流程图。
76.图3为装车计划品种煤最大装车趟数确定方法流程图。
具体实施方式
77.下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
78.本发明提供的一种确定选煤厂品种煤装车计划的动态价值寻优方法，包括以下步骤：
79.第一步：取得需要装车品种煤的基本信息以及各筒仓的煤质和煤量信息，具体如下：
80.(1)品种煤相关信息
81.品种煤k，k表示品种煤的序号，如选煤厂有品种煤1和品种煤2两种品种煤，即k的最大值为2。
82.mtmbk，表示第k种品种煤的全水分控制范围，％；mtmin表示全水分控制下限，mtmax表示全水分控制上限；
83.admbk，表示第k种品种煤的干基灰分控制范围，％；admin表示干基灰分控制下限，admax表示干基灰分控制上限；
84.qmbk，表示第k种品种煤的收到基水分控制范围，kcal/kg；qmin表示收到基低位发热量控制下限，qmax表示收到基低位发热量控制控制上限；
85.可以根据用户需求对mtmbk、admbk、qmbk等指标进行单一指标或者组合指标进行限定；pk，表示第k种品种煤煤质达标后的原煤价格，元/吨；
86.dwk，表示第k种品种煤每趟车的装车煤量，吨；
87.(2)各筒仓信息：
88.筒仓i：表示第i个筒仓；
89.mtci，表示第i个筒仓煤样的全水分，％；
90.adci，表示第i个筒仓煤样的干基灰分，％；
91.qci，表示第i个筒仓煤样的收到基低位发热量，kcal/kg；
92.煤量mli，表示第i个筒仓的初始煤样重量，t。
93.(3)装车趟数：
94.zzcts，计划的总装车趟数，趟。
95.n，根据筒仓煤量推荐的品种煤装车趟数，趟。
96.第二步：判别不同混煤的品种煤类别及价值，具体方法如下并见图1；
97.201)穷举筒仓在不同比例掺配时的混煤煤质，具体见式(1)到式(5)，r1、r2、r3分别表示三个仓的配比，％，mth、aarh、qh表示掺配后混煤的全水分、收到基灰分、收到基低位发热量，adh表示掺配后混煤的干基灰分，j表示三个筒仓的掺配比例按照1％步长变化计算的次数：
98.mthj＝(r1j
×
mtc1+r2j
×
mtc2+r3j
×
mtc3)/100
ꢀꢀꢀ
(1)
99.aarhj＝(r1j
×
aarc1+r2j
×
aarc2+ r3j
×
aarc3)/100
ꢀꢀ
(2)
100.qhj＝(r1j
×
qc1+r2 j
×
qc2+ r3j
×
qc3)/100
ꢀꢀ
(3)
101.r1j +r2j +r3j ＝100
ꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
102.adhj＝ aarhj
×
100/(100
‑ꢀ
mthj)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
103.202)判断每一个掺配比例下的煤种是“品种煤1”、“品种煤2”，还是“非品种煤”：
104.a)如果mthj、adhj、qhj的指标满足品种煤1，则标记为“品种煤1”，该混煤的单价pd1j＝p1/qhj；
105.b)如果mthj、adhj、qhj的指标满足品种煤2，则标记为“品种煤2”，该混煤的单价pd2j＝p2/qhj；
106.c)如果混煤的煤质指标既满足品种煤1又满足品种煤2的要求，则对比每大卡热值的价值判断该煤种属于哪种煤，具体措施如下：
107.如果该煤种的pd1j≥pd2j，则最终标记该煤种为“品种煤1”，否则标记为“品种煤2”，
108.d)煤质既不满足品种煤1也不满足品种煤2的要求，则该混煤标记为：“非品种煤”；
109.第三步：求解所有穷举中判定为“品种煤1”、“品种煤2”的混煤的最大装车趟数njmax，求解方法如下及图2：
110.301)求解第j种掺配比例下的最大装车趟数njmax，j从1开始；
111.a)假定在该种掺配方案下的最大装车趟数njmax＝1，
112.b)在njmax装车趟数时需要的每个筒仓的煤量mlnji，t，具体见式(6)
113.mlnij＝rij
×
dwk
×
njmax
ꢀꢀꢀꢀ
(6)
114.mlnij：表示第j种掺配比例下配置第j种混煤时需要的第i个筒仓的煤量，t；
115.dwk：表示第j种掺配比例下，当混煤判定为品种煤1时，dwk对应品种煤1的装车吨位，当判定为品种煤2时，dwk对应品种煤2的装车吨位；
116.c)如果mlnij≤vgi，最大装车趟数增加1趟，则njmax＝njmax+1，再次从a)到b)循环找出最大njmax，即各个筒仓在满足煤量的条件下，最多能装多少趟车该品种煤；
117.d)确定在各筒仓装完njmax趟品种煤时剩余煤量mlsij(吨)，和品种煤的总价值pjmax(万元)，具体计算分别见式(7)和式(8)
118.mlsij＝mli-mlnij
ꢀꢀꢀꢀ
(7)
119.pjmax＝njmax
×
dwk
×
pk/10000
ꢀꢀꢀ
(8)
120.302)重复步骤301)，直到求解得出所有穷举中判定为“品种煤1”、“品种煤2”在最大装车趟数njmax和总价值pjmax；
121.第四步：第f次确定最优装车品种煤及对应装车趟数nmaxf、各筒仓需要的煤量
mlnif，t，并最终得出品种煤的最终推荐装车趟数，其中f从1开始，f表示第f次求解最优装车趟数及煤种，求解方法如下及图3：
122.401)确定最优装车品种煤及对应装车趟数
123.情况1：最大装车趟数只对应一个掺配比例，则第三步中装车趟数最多的确定品种煤确定为最优装车品种煤，；
124.情况2：当装车趟数最多的对应多个配比时，确定总价值最高的为最优装车品种煤；
125.情况3：当装车趟数最多的对应多个配比时，且总价值最高的也有多种配比时，确定混煤发热量qhj最低的为最优装车品种煤；
126.情况4：当前最大装车趟数均为0，则表示本次无法配出品种煤。
127.满足情况1或者情况2或者情况3时，标记该计算工况为最优装车工况，则装完本次推荐装车计划后各筒仓需要的煤量mlnif，t；计算参见式(9)：
128.mlnif＝rif
×
dwf
×
nmaxf
ꢀꢀꢀꢀ
(9)
129.rif：第f次求解时，最优装车品种煤对应的各筒仓的掺配比例，％；
130.dwf：第f次求解时，最优装车品种煤对应的每趟车的装车煤量，t；
131.装完本次推荐最优装车品种煤，此时各筒仓的剩余煤量mlsif，t；计算参见式(10)：
132.mlsif＝mli
‑ꢀ
mlnif
ꢀꢀꢀꢀ
(10)
133.402)求解下一个推荐装车品种煤及趟数。
134.当mls1f+mls2f+mls3f》0，nf《zzcts，在当前筒仓煤量下即mlif＝mlsif，重复第三步到第四步的步骤401)，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计品种煤1和品种煤2的装车趟数，得到总的品种煤装车趟数nf。
135.403)求解最终推荐的品种煤装车趟数n：
136.如果推荐的品种煤的总装车趟数nf≤zzcts，则品种煤的推荐装车趟数为n＝nf，否则品种煤的推荐装车趟数为n＝zzcts。
137.实施例
138.一种确定选煤厂品种煤装车计划的动态价值寻优方法，包括以下步骤：
139.第一步：取得需要装车品种煤的基本信息以及各筒仓的煤质和煤量信息，具体如下：
140.(1)品种煤相关信息
141.品种煤k，k表示品种煤的序号，选煤厂通常有两种品种煤，即k＝2。
142.mtmbk，表示第k种品种煤的全水分控制范围，％；
143.admbk，表示第k种品种煤的干基灰分控制范围，％；
144.qmbk，表示第k种品种煤的收到基水分控制范围，kcal/kg；
145.pk，表示第k品种煤煤质达标后的原煤价格，元/吨；
146.dwk，表示第k种品种煤每趟车的装车煤量，吨；
147.选煤厂通常有两种品种煤，分别为品种煤1和品种煤2。
148.(2)各筒仓信息：
149.筒仓i：表示第i个筒仓；
150.mtci，表示第i个筒仓煤样的全水分，％；
151.adci，表示第i个筒仓煤样的干基灰分，％；
152.qci，表示第i个筒仓煤样的收到基低位发热量，kcal/kg；
153.煤量mli，表示第i个筒仓的初始煤样重量，t。
154.(3)装车趟数：
155.zzcts，计划的总装车趟数，趟。
156.n，根据筒仓煤量推荐的品种煤装车趟数，趟。
157.本次品种煤、各筒仓以及计划总装车趟数信息如下：
158.(1)品种煤相关信息
[0159][0160]
(2)各筒仓信息：
[0161] 筒仓1筒仓2筒仓3 块精煤末精煤混煤筒仓煤量vgi/吨7000800015000全水分mtci/％14.815.718.3干基灰分adci/％4.206.3010.68收到基低位发热量qci/kcal/kg604058275338
[0162]
(3)装车趟数：
[0163]
本次计划的总装车趟数zzcts＝14趟
[0164]
第二步：判别不同混煤的品种煤价值，具体方法如下并见图1；
[0165]
201)穷举筒仓在不同比例掺配时的混煤煤质，具体见式(1)到式(5)，r1、r2、r3分别表示三个仓的配比，％，mth、aarh、qh表示掺配后混煤的全水分、收到基灰分、收到基低位发热量，adh表示掺配后混煤的干基灰分，j表示三个筒仓的掺配比例按照1％步长变化计算的次数：
[0166]
mthj＝(r1j
×
mtc1+r2j
×
mtc2+r3j
×
mtc3)/100
ꢀꢀꢀ
(1)
[0167]
aarhj＝(r1j
×
aarc1+r2j
×
aarc2+r3j
×
aarc3)/100
ꢀꢀ
(2)
[0168]
qhj＝(r1j
×
qc1+r1 j
×
qc2+ r1 j
×
qc3)/100
ꢀꢀ
(3)
[0169]
r1j +r2j +r3j ＝100
ꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0170]
adhj＝ aarhj
×
100/(100
‑ꢀ
mthj)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0171]
202)判断每一个掺配比例下的煤种是“品种煤1”、“品种煤2”，还是“非品种煤”：
[0172]
a)如果mthj、adhj、qhj的指标满足品种煤1，则标记为“品种煤1”，该混煤的单价pd1j＝p1/qhj；
[0173]
b)如果mthj、adhj、qhj的指标满足品种煤2，则标记为“品种煤2”，该混煤的单价pd2j＝p2/qhj；
[0174]
c)如果混煤的煤质指标既满足品种煤1又满足品种煤2的要求，则对比每大卡热值的价值判断该煤种属于哪种煤，具体措施如下：
[0175]
如果该煤种的pd1j≥pd2j，则最终标记该煤种为“品种煤1”，否则标记为“品种煤2”，
[0176]
d)煤质既不满足品种煤1也不满足品种煤2的要求，则该混煤标记为：“非品种煤”；
[0177]
本次按照20％的掺配比例进行穷举，则最终确定不同掺配比例下混煤的品种煤类别如下表：
[0178]
[0179][0180]
第三步：求解所有穷举中判定为“品种煤1”、“品种煤2”的混煤的最大装车趟数njmax，求解方法如下及图2：
[0181]
301)求解第j种掺配比例下的最大装车趟数njmax，j从1开始；
[0182]
a)假定在该种掺配方案下的最大装车趟数njmax＝1，
[0183]
b)在njmax装车趟数时需要的每个筒仓的煤量mlnji，t，具体见式(6)
[0184]
mlnij＝rij
×
dwk
×
njmax
ꢀꢀꢀꢀ
(6)
[0185]
mlnij：表示第j种掺配比例下，配置第j种混煤时需要的第i个筒仓的煤量，t；
[0186]
dwk：表示第j种掺配比例下，当混煤判定为品种煤1时，dwk对应品种煤1的装车吨位，当判定为品种煤2时，dwk对应品种煤2的装车吨位；
[0187]
c)如果mlnij≤vgi，最大装车趟数增加1趟，则njmax＝njmax+1，再次从a)到b)循环找出最大njmax，即各个筒仓在满足煤量的条件下，最多能装多少趟车该品种煤；
[0188]
d)确定在各筒仓装完njmax趟品种煤时剩余煤量mlsij，和品种煤的总价值pjmax，具体计算分别见式(7)和式(8)
[0189]
mlsij＝mli-mlnij
ꢀꢀꢀꢀ
(7)
[0190]
pjmax＝njmax
×
dwk
×
pk/10000
ꢀꢀꢀ
(8)
[0191]
302)重复步骤301)，直到求解得出所有穷举中判定为“品种煤1”、“品种煤2”在最大装车趟数njmax和总价值pjmax；
[0192]
本次按照20％的掺配比例进行穷举，计算第三步表中判定为品种煤的最大装车趟数njmax及装车总价值pjmax的计算结果如下：
[0193][0194][0195]
第四步：第f次确定最优装车品种煤及对应装车趟数nmaxf、各筒仓需要的煤量mlnif，t，并最终得出品种煤的最终推荐装车趟数，其中f从1开始，f表示第f次求解最优装
车趟数及煤种，求解方法如下及图3：
[0196]
401)确定最优装车品种煤及对应装车趟数
[0197]
情况1：最大装车趟数只对应一个掺配比例，则第三步中装车趟数最多的确定品种煤确定为最优装车品种煤，；
[0198]
情况2：当装车趟数最多的对应多个配比时，确定总价值最高的为最优装车品种煤；
[0199]
情况3：当装车趟数最多的对应多个配比时，且总价值最高的也有多种配比时，确定混煤发热量qhj最低的为最优装车品种煤；
[0200]
情况4：当前最大装车趟数均为0，则表示本次无法配出品种煤。
[0201]
满足情况1或者情况2或者情况3时，标记该计算工况为最优装车工况，则装完本次推荐装车计划后各筒仓需要的煤量mlnif，t；计算参见式(9)：
[0202]
mlnif＝rif
×
dwf
×
nmaxf
ꢀꢀꢀꢀ
(9)
[0203]
rif：第f次求解时，最优装车品种煤对应的各筒仓的掺配比例，％；
[0204]
dwf：第f次求解时，最优装车品种煤对应的每趟车的装车煤量，t；
[0205]
装完本次推荐最优装车品种煤，此时各筒仓的剩余煤量mlsif，t；计算参见式(10)：
[0206]
mlsif＝mli
‑ꢀ
mlnif
ꢀꢀꢀꢀ
(10)
[0207]
402)求解下一个推荐装车品种煤及趟数。
[0208]
当mls1f+mls2f+mls3f》0，nf《zzcts，在当前筒仓煤量下即mlif＝mlsif，重复第三步到第四步的步骤401)，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计品种煤1和品种煤2的装车趟数n1和n2，得到总的品种煤装车趟数nf。
[0209]
403)求解最终推荐的品种煤装车趟数n：
[0210]
如果推荐的品种煤的总装车趟数nf≤zzcts，则品种煤的推荐装车趟数为n＝nf，否则品种煤的推荐装车趟数为n＝zzcts。
[0211]
首先根据第三步表格中的装车总价值pjmax，混煤8、混煤9、混煤10以及混煤13的总价值均为1296.00万元，其中混煤10的发热量最低，推荐为第一次的优先装车方案，即推荐按照混煤10的配比装4趟品种煤2，
[0212]
即：n2＝4，
[0213]
装完4趟混煤10以后，各筒仓剩余煤量mls11＝88吨、mls21＝8000吨、mls31＝4632吨，
[0214]
第二次计算各品种煤的最大装车趟数njmax2，具体结果如下：
[0215][0216]
根据上表，本次混煤17和混煤18的最大装车趟数均为2趟，装车总价值也均为648万元，则推荐发热量较低的混煤18为第二次推荐品种煤装车，即本次推荐再装2趟品种煤2，配煤比例为混煤18。
[0217]
装完2趟混煤18以后，各筒仓剩余煤量mls12＝88吨、mls22＝2816吨、mls32＝1176吨，
[0218]
第三次计算各品种煤的最大装车趟数njmax3，根据下表可见已不能装任何品种煤：
[0219][0220]
根据第一次计算结果和第二次结算结果，最终推荐的品种煤装车计划如下：
[0221][0222]
[0223]
本次可装品种煤的趟数nf＝n2＝4+2＝6，
[0224]
本次nf《zzcts，所以最终推荐品种煤装车趟数n＝nf＝6
[0225]
即推荐装6趟品种煤2，总品种煤价值为1944万元。
[0226]
虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

技术特征：
1.一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，包括以下步骤：第一步：取得需要装车品种煤的基本信息以及各筒仓的煤质和煤量信息；第二步：判别不同混煤的品种煤价值；第三步：求解所有穷举中判定为品种煤的混煤的最大装车趟数；第四步：确定最优装车品种煤及对应装车趟数、各筒仓需要的煤量；第五步：装完第四步推荐的煤种煤后，在剩余筒仓煤量下，重复第三步到第四步，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计不同品种煤的装车趟数，得到总的品种煤装车趟数。2.根据权利要求1所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第一步中，具体如下：(1)品种煤相关信息品种煤k，k表示品种煤的序号，选煤厂通常有两种品种煤，即k最大为2；mtmbk，表示第k种品种煤的全水分控制范围，％；mtmink表示全水分控制下限，mtmaxk表示全水分控制上限；admbk，表示第k种品种煤的干基灰分控制范围，％；admink表示干基灰分控制下限，admaxk表示干基灰分控制上限；qmbk，表示第k种品种煤的收到基水分控制范围，kcal/kg；qmink表示收到基低位发热量控制下限，qmaxk表示收到基低位发热量控制控制上限；根据用户需求对mtmbk、admbk和qmbk指标进行单一指标或者组合指标进行限定；pk，表示第k品种煤煤质达标后的原煤价格，元/吨；dwk，表示第k种品种煤每趟车的装车煤量，吨；(2)各筒仓信息：筒仓i：表示第i个筒仓；mtci，表示第i个筒仓煤样的全水分，％；adci，表示第i个筒仓煤样的干基灰分，％；qci，表示第i个筒仓煤样的收到基低位发热量，kcal/kg；煤量mli，表示第i个筒仓的初始煤样重量，t；(3)装车趟数：zzcts，计划的总装车趟数，趟；n，根据筒仓煤量推荐的品种煤装车趟数，趟。3.根据权利要求2所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第二步中，具体方法如下201)穷举筒仓在不同比例掺配时的混煤煤质，具体见式(1)到式(5)，r1、r2、r3分别表示三个仓的配比，％，mth、aarh、qh表示掺配后混煤的全水分、收到基灰分、收到基低位发热量，adh表示掺配后混煤的干基灰分，j表示三个筒仓的掺配比例按照1％步长变化计算的次数：mthj＝(r1j
×
mtc1+r2j
×
mtc2+r3j
×
mtc3)/100
ꢀꢀꢀ
(1)aarhj＝(r1j
×
aarc1+r2j
×
aarc2+ r3j
×
aarc3)/100
ꢀꢀ
(2)qhj＝(r1j
×
qc1+r2 j
×
qc2+ r3j
×
qc3)/100
ꢀꢀ
(3)
r1j +r2j +r3j ＝100
ꢀꢀꢀꢀꢀꢀꢀꢀ
(4)adhj＝ aarhj
×
100/(100
‑ꢀ
mthj)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)202)判断每一个掺配比例下的混煤煤是“品种煤1”、“品种煤2”，还是“非品种煤”。4.根据权利要求3所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第二步的202)中，a)如果mthj、adhj、qhj的指标满足品种煤1，则标记为“品种煤1”，该混煤的单价pd1j＝p1/qhj；b)如果mthj、adhj、qhj的指标满足品种煤2，则标记为“品种煤2”，该混煤的单价pd2j＝p2/qhj；c)如果混煤的煤质指标既满足品种煤1又满足品种煤2的要求，则对比每大卡热值的价值判断该煤种属于哪种煤，具体措施如下：如果该煤种的pd1j≥pd2j，则最终标记该煤种为“品种煤1”，否则标记为“品种煤2”，d)煤质既不满足品种煤1也不满足品种煤2的要求，则该混煤标记为：“非品种煤”。5.根据权利要求3所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第三步中，求解方法如下：301)求解第j种掺配比例下的最大装车趟数njmax，j从1开始；a)假定在该种掺配方案下的最大装车趟数njmax＝1，b)在njmax装车趟数时需要的每个筒仓的煤量mlnij，t，具体见式(6)mlnij＝rij
×
dwk
×
njmax
ꢀꢀꢀꢀ
(6)mlnij：表示第j种掺配比例下，配置第j种混煤时需要的第i个筒仓的煤量，t；dwk：表示第j种掺配比例下，当混煤判定为品种煤1时，dwk对应品种煤1的装车吨位，当判定为品种煤2时，dwk对应品种煤2的装车吨位；c)如果mlnij≤vgi，最大装车趟数增加1趟，则njmax＝njmax+1，再次从a)到b)循环找出最大njmax，即各个筒仓在满足煤量的条件下，最多能装多少趟该品种煤；d)确定在各筒仓装完njmax趟品种煤时剩余煤量mlsij，和品种煤的总价值pjmax，具体计算分别见式(7)和式(8)mlsij＝mli-mlnij
ꢀꢀꢀꢀ
(7)pjmax＝njmax
×
dwk
×
pk/10000
ꢀꢀꢀ
(8)302)重复步骤301)，直到求解得出所有穷举中判定为“品种煤1”、“品种煤2”的最大装车趟数njmax和总价值pjmax。6.根据权利要求5所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第四步中，求解方法如下：情况1：最大装车趟数只对应一个掺配比例，则第三步中装车趟数最多的确定为最优装车品种煤；情况2：当装车趟数最多的对应多个配比时，确定总价值最高的为最优装车品种煤；情况3：当装车趟数最多的对应多个配比时，且总价值最高的也有多种配比时，确定混煤发热量qhj最低的为最优装车品种煤；情况4：当前最大装车趟数均为0，则表示本次无法配出品种煤；满足情况1或者情况2或者情况3时，标记该计算工况为最优装车工况，则装完本次推荐装车计划后各筒仓需要的煤量mlnif，t；计算参见式(9)：
mlnif＝rif
×
dwf
×
nmaxf
ꢀꢀꢀꢀ
(9)rif：第f次求解当前煤量下的最优装车品种煤及对应最大装车趟数时，最优装车品种煤对应的各筒仓的掺配比例，％；dwf：第f次求解当前煤量下的最优装车品种煤及对应最大装车趟数时，最优装车品种煤对应的每趟车的装车煤量，t；nmaxf：第f次求解当前煤量下的最优装车品种煤及对应最大装车趟数时，最优装车品种煤对应的最大装车趟数，t；装完本次推荐最优装车品种煤，此时各筒仓的剩余煤量mlsif，t；计算参见式(10)：mlsif＝mli
‑ꢀ
mlnif
ꢀꢀꢀꢀ
(10)。7.根据权利要求6所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第五步中，求解下一个推荐装车品种煤及趟数如下：当mls1f+mls2f+mls3f>0，nmaxf<zzcts，在当前筒仓煤量下即mlif＝mlsif，重复第三步到第四步的步骤401)，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计品种煤1和品种煤2的装车趟数，得到总的品种煤装车趟数nf。8.根据权利要求7所述的一种确定选煤厂品种煤装车计划的动态价值寻优方法，其特征在于，第五步中，求解最终推荐的品种煤装车趟数n如下：如果推荐的品种煤的总装车趟数nf≤zzcts，则品种煤的推荐装车趟数为n＝nf，否则品种煤的推荐装车趟数为n＝zzcts。

技术总结
本发明公开了一种确定选煤厂品种煤装车计划的动态价值寻优方法，包括：第一步：取得需要装车品种煤的基本信息以及各筒仓的煤质和煤量信息；第二步：判别不同混煤的品种煤价值；第三步：求解所有穷举中判定为“品种煤1”、“品种煤2”的混煤的最大装车趟数；第四步：确定最优装车品种煤及对应装车趟数、各筒仓需要的煤量；第五步：装完第四步推荐的煤种煤后，在剩余筒仓煤量下，重复第二步到第四步，求解下一个最优装车品种煤及对应装车趟数，直到新增品种煤装车趟数为0或者总品种煤装车趟数大于总装车趟数；分别统计品种煤1和品种煤2的装车趟数，得到总的品种煤装车趟数。本发明能够实现品种煤的价值最大化，实现选煤厂的品种煤效益最大化。最大化。最大化。


技术研发人员：刘家利 张喜来 晋中华 杨忠灿 屠竞毅 方顺利 杜智华 李炎 李宇航 周广钦 李仁义 孙军 郭洋洲 张森 张向宇 贾子秀
受保护的技术使用者：西安热工研究院有限公司
技术研发日：2023.07.18
技术公布日：2023/10/20