Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov...

14
Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov * , Mariusz Kozlowski, Zlatogor Minchev IN COOPERATION WITH KCM AD, PLOVDIV*

Transcript of Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov...

Page 2: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

ZINC HYDROMETALLURGY DEMO VIDEO

PROBLEM DEFINITION

METHODOLOGY & SOFTWARE SELECTION

RESULTS

DISCUSSION

Page 4: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

THE MAIN GOAL IS TO OPTIMIZE THE CHARGING PROCESS OF THE ZINC METALLURGY FOLLOWING THE INPUT COMPONENTS MIXTURE

(“CONCENTRATE” – “CHARGE”) TECHNOLOGICAL PRICE - CT

L(x) =

CT = C1 - C2 - C3,

CTj - price of the j-th concentrate; xj – quantity of the j-th concentrate

→ CTj * xj → min,

Further on, for simplicity CTj and its three components are denoted

without the j-th index.

Page 5: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

C1 – depends on the Zn contents in the concentrate as follows:

IF [(1-p) * Zn < q THEN (Zn-q) * Cd /100) ELSE p * Zn * Cd /100 ]

Cd – daily price of Zn on London Metal Exchange (LME)

p – final agreed Zn content at official LME (85%)

q – minimum deduction of Fe at the official LME (8%)

Where:

Page 6: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

C2 – the correction that depends on the difference

between Cd and Cavg of the Zn:

Cd – daily price of Zn on London Metal Exchange (LME)

C2 = Eavg – (Cavg - Cd) / 10; Eavg – average expenses for 1 tone concentrate processing and transport, Eavg = 2000 $

C3 – penalty for Fe contents in the concentrate:

[IF Fe < q THEN 0 ELSE (Fe- q)* s]

q – minimum deduction of Fe at the official LME (8%)

s – the penalty for 1 % Fe contents in the concentrate, s = 2$

Page 7: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

LINEAR OPTIMIZATION - assuming that the components mixture is a linear combination of the concentrate components

MS EXCEL SOLVER – well-known and commonly available

Page 8: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

Input concentrate components n, n=21

Page 9: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

LINEAR OPTIMIZATION INTEGER LINEAR OPTIMIZATION

49

3.9

2

Tota

l Pri

ce (

CT)

49

3.1

8

Tota

l Pri

ce (

CT)

The Integer optimization provides results that require discretization change in accordance with the real production system mechanical capabilities improvements.

Page 10: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

No Boundary Conditions of

the charge components

min Price = 543.55 $/t

Linea

r m

od

el

No Boundary Conditions of

the charge components

min Price = 549.61 $/t Lin

ear In

t mo

del

divisib

le on

1%

No Boundary Conditions of

the charge components

min Price = 558.03 $/t

Linea

r Int m

od

el d

ivisible o

n 5

%

No Boundary Conditions of

the charge components

min Pb+Cu+SiO2= 4.40 %

Linea

r m

od

el

No Boundary Conditions of

the charge components

Linea

r Int m

od

el d

ivisible o

n 1

%

No Boundary Conditions of

the charge components

Linea

r Int m

od

el d

ivisible o

n 5

%

Price < 560 $/t

min Pb+Cu+SiO2= 4.88 %

Price < 560 $/t

Price < 560 $/t

min Pb+Cu+SiO2= 5.15 %

I

Page 11: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

No Boundary Conditions of

the charge components

min Pb+Cu+SiO2= 4.40 %

Linea

r m

od

el

No Boundary Conditions of

the charge components

Linea

r Int m

od

el d

ivisible o

n 1

%

No Boundary Conditions of

the charge components

Linea

r Int m

od

el d

ivisible o

n 5

%

Price < 560 $/t

min Pb+Cu+SiO2= 4.88 %

Price < 560 $/t

Price < 560 $/t

min Pb+Cu+SiO2= 5.15 %

No Boundary Conditions of

the charge components

min Fe = 7.64 %

Linea

r m

od

el

No Boundary Conditions of

the charge components

Linea

r Int m

od

el d

ivisible o

n 1

%

No Boundary Conditions of

the charge components

Linea

r Int m

od

el d

ivisible o

n 5

%

Price < 560 $/t

Price < 560 $/t

Price < 560 $/t

Pb+Cu+SiO2 < 5.20 %

Pb+Cu+SiO2 < 5.20 %

Pb+Cu+SiO2 < 5.20 %

min Fe = 7.82 %

min Fe = 7.98 %

II

Page 12: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

No Boundary Conditions of

the charge components

min Price = 543.55 $/t

Linea

r m

od

el

No Boundary Conditions of

the charge components

min Price = 549.61 $/t Lin

ear In

t mo

del

divisib

le on

1%

No Boundary Conditions of

the charge components

min Price = 558.03 $/t

Linea

r Int m

od

el d

ivisible o

n 5

%

Lower & Upper Boundary Conditions

Linea

r m

od

el

Lower & Upper Boundary Conditions

Linea

r Int m

od

el d

ivisible o

n 1

%

Lower & Upper Boundary Conditions

Linea

r Int m

od

el d

ivisible o

n 5

%

min Price = 548.83 $/t

min Price = 553.41 $/t

No feasible solution

III

Page 13: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc

The developed mathematical model based on liner programming techniques for industrial support of zinc metallurgy charging process successfully demonstrated the practical problems between theoretical and practical industrial requirements meeting. Evidently the integer liner programming is a good tool for industrial application together with MS Excel software. Unfortunately, not always a feasible solution is possible to be found

Page 14: Doychin Boyadzhiev, Zlatko Varbanov, Alexander Peltekov ...esgi95.fmi.uni-sofia.bg/presentations/03.pdfmariusz kozlowski, zlatogor minchev in cooperation with kcm ad, plovdiv* zinc