Запишем уравнения, используя исходные данные и принятые значения для потоков вольфрама.

= 7050,75 кг

1) измельчение и обогащение: в твердосплавный концентрат переходит 7050,75 кг соответственно 0,2 и 0,8 количества молибдена, поступающего на операцию;

2) возгонка: в остатке от возгонки 0,3, в возгонах 0,7 количества молибдена;

3) выщелачивание: в растворе 0,98, в отвалах выщелачивания 0,02;

4) очистка от тяжелых металлов: в сульфидном кеке 0,05, в очищенном растворе 0,95;,.

5) нейтрализация: в кислом маточном растворе 0,01, в кристаллах тетрамолибдата аммония 0,99;

6) перекристаллизация: в маточном растворе 0,1, в кристаллах парамолибдата аммония (конечный продукт) 0,9. /

Зададимся производительностью по молибдену в конечном продукте 1000 т/год.

Обозначим потоки молибдена (т.е. количества его в различных материалах), т/год, неизвестными x-i; Хг, - ', Х^2, *13 в соответствии с рис.7.

Для определения неизвестных составим систему, состоящую их уравнений трех типов.

1. Уравнение, в котором задана производительность:

х-i з=1000.

2. Уравнения, в которых использованы соотношения между потоками. Это могут быть непосредственно содержащиеся в исходных данных сведения о доле ценного компонента от количества, поступающего на операцию, переходящей в каждый из продуктов этой операции, например:,

Х2 = 0,2 (xl + Х2), Х3 = 0,8 (Xi + X2), Х4 = 0,3 Х3,

Х5 = 0,7, Х3,

или вытекающие из этих сведений соотношения между потоками, выходящими с одной и той же операции, например:

Х3 = 4Х2,... .

Х5= 7/3 Х4.

Кроме того, возможно использование соотношения между потоками, относящимися к разным операциям схемы. Дело в том, что в реальных технологических схемах количества ценного компонента в тех или иных материалах (особенно в выводимых из схемы) часто задают как долю от количества ценного компонента в исходном материале или конечном продукте. Так, например, исходные данные. могли бы содержать,:; следующее указание: "Потери молибдена с, остатком выщелачивания составляют в среднем 0,6% от количества, поступающего с молибденитовым концентратом"; соответствующее, уравнение имело бы вид

Х7= 0,006 Xi.

Следует отметить, что в тех случаях, когда вместо данных, относящихся к потокам одной и той же операции, заданы соотношения между потоками разных операций, расчет методом "от операции к операции" невозможен.

3. Уравнения, отражающие равенство количества - ценного компонента, поступающего на операцию, количеству, выходящему с нее:;

Ха + Хз = Xi + *2 • или Хз = xl;

х4 + х5 = х3, Х6 + Х7 = Х4,

Х8 + Х9 = Хб + Х12 И Т.Д.

Легко убедиться, что общее число уравнений всех типов, которые можно составить, во много раз превышает число неизвестных, хотя система должна иметь единственное решение и соответственно число уравнений должно быть равно числу неизвестных. Причина, очевидно, состоит в том, что большая часть уравнений представляет собой линейные комбинации других; например, из

Х2 = 0,2 (Хт + Х2),

Х3 = 0,8 (Xi + Х2) следует

Ха + Хз = Xi + х2, или Хз = Xi;

Хз: Х2 = 0,8: 0,2, или Х3 = 4Х2.

В то же время решение возможно только при условии, что при числе уравнений, равном числу неизвестных, все уравнения линейно независимы. Поэтому после составления системы уравнений необходим тщательный контроль отсутствия в ней линейно зависимых (дублирующих) уравнений.

По-видимому, дублирования уравнений проще всего избежать, используя в системе, кроме уравнения, задающего производительность, только соотношения между количеством ценного компонента в потоке, выходящем с операции, и количеством, поступающим на эту операцию. Для схемы, показанной на рис.7, можно составить, например, следующую систему уравнений для определения потоков ценного компонента при производительности по конечному продукту, равной 1000:

1. Х2 = 0,2 (х-i + Х2), или Х2 = 0,25 Xi;

2. х3 = 0,8(х1 +х2);

3. х4 = о. з Хз;

4. Х5 = 0,7 Хз;

5. Хб = 0,98 Х4;

б ху = 0,02 х4;

7. х8 = 0,05 (х6 + х12);

8. Х9 = 0,95(Х6 + Х12);

э. хю = 0,01 х9;

10. Хц =0,99 Xg;

11. х12 = о,1 (хб + хц);

12. Х13 = 0,9(Х5 + Хц);

13. Х13 = 1000.

Подобную систему уравнений можно довольно легко решить вручную, путем последовательной подстановки, сокращения переменных при вычитании одних уравнений из других и т.д. Однако если число неизвестных очень велико (а в реальных схемах число потоков может достигать многих десятков и даже сотен), возрастает трудоемкость расчетов и вероятность ошибок. В подобных случаях для сокращения числа уравнений можно рекомендовать обозначать неизвестными не количества ценного компонента в каждом из потоков, а суммарные количества, поступающие на операции схемы; при этом, очевидно, число уравнений на 1 больше числа операций. Например, для той же схемы (рис.7) получаем:

1. У1 = G ucx + 0,2 /!, или 0,8 у! = G исх;

2. у2 = 0,8 yi;

3. уз = о, з у2;

4. у4 = 0,98 Уз + 0,1 Уб',

5. у5 = 0,95 у4;

6. у6 = 0,7 у2 + 0,99 у5;

7. G исх = 0,9 у6.

После решения подобной системы уравнений расчет количества ценного компонента в каждом из потоков не вызывает затруднений.

Однако наиболее эффективным способом преодоления трудностей расчета сложных технологических схем является использование компьютеров.

В отличие от человека, выбирающего для каждой конкретной системы уравнений наиболее рациональный путь решения, в программах для цифровых вычислительных машин можно использовать только универсальные, единые для всех систем линейных уравнений способы вычисления. Среди таких способов наиболее распространены метод Гаусса с выбором главного элемента столбца или строки и метод обращения матрицы [2].

До начала вычислений необходимо ввести исходные данные: при использовании метода Гаусса - расширенную матрицу коэффициентов системы линейных уравнений (значения коэффициентов при неизвестных и свободные члены каждого из уравнений):

а при использовании метода обращения матрицы - отдельно квадратную матрицу коэффициентов и вектор-столбец свободных членов:

Для составления матрицы коэффициентов члены уравнений, содержащие неизвестные, необходимо расположить слева от знака равенства в порядке возрастания индекса неизвестного, оставив справа только свободные члены; отсутствующие неизвестные вносят в уравнения с коэффициентами, равными нулю. После этого можно приступить к вводу матрицы коэффициентов.

Однако применительно к системе уравнений, получаемой при описании распределения ценного компонента по потокам технологической схемы, этот способ нерационален, так как подавляющая часть коэффициентов равна нулю. Например, при вводе коэффициентов составленной ранее системы из 13 уравнений первые две строки должны быть записаны в следующем виде: - 0,25 Xi + 1 Х2 + О Х3 + О Х4 + О Х5 + О Х6 + О Х7 + О Х8 + О Х9 +

+ о х10 + о хп + ох12+ о х13 = о;

- 0,8 х1 - 0,8 х2 + 1 хЗ + 0 х4 + 0 х5 + 0 хб + 0 х7 + 0 х8 + О Х9 +

+ о Хю+ о х-и +о X-I2 + о х13 = о,

а первые две строки расширенной матрицы коэффициентов соответственно

-0,25 1 000000000000 - 0,8 - 0,8 000000000000

При большом числе неизвестных количество вводимых нулей становится громадным: например, при 40 неизвестных расширенная матрица состоит из 40 х 41 = 1640 коэффициентов, из которых более 1500 будут равны нулю. Очевидно, что ввод подобной матрицы настолько трудоемок и неизбежно сопровождается таким количеством ошибок, что превращается в сложную задачу.

Это затруднение устраняется, если ввод матрицы осуществлять в два этапа: сначала заполнить всю матрицу нулями (эта операция выполняется очень легко), а затем ввести ненулевые коэффициенты, заменяя ими нули.

Ниже описан расчет балансов по ценному компоненту методом Гаусса с помощью программы на языке BASIC и методом обращения матрицы с помощью табличного процессора EXCEL [3, 4].