Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Эколого-экономические системы

..pdf
Скачиваний:
3
Добавлен:
05.02.2023
Размер:
1.71 Mб
Скачать

Занятие № 3. Задача оптимизации параметров процесса «Откачка нефти

с грунта» и ее автоматизированное решение

Задача оптимизации параметров технологического процесса ПОМ

«Откачка нефти с грунта» решена в несколько этапов:

Этап 1. Понимание системы. В процессе выполнения ПОМ «Откачка нефти с грунта» взаимодействуют следующие компоненты ЭЭС (рис. 3.1):

грунт и грунтовые воды, загрязненные нефтью;

вакуумный самосвал, выполняющий сбор нефти с грунта (самосвал);

шламовый амбар, предназначенный для временного хранения собранной загрязнённой нефтью почвы и жидкости;

фонд, расходуемый в процессе выполнения.

Грунт

Загружает

 

Разгружает

Шламовый

 

грунт

Самосвал

грунт

амбар

 

 

 

 

Vt

 

R1t

 

R2t

 

 

 

 

 

 

 

 

 

 

Затрачивает

Фонд Фt

Рис. 3.1. Взаимодействие компонентов ЭЭС в процессе выполнения ПОМ «Откачка нефти с грунта»

Этап 2. Формулировка цели моделирования системы. Необходимо

определить плановые показатели ПОМ «Откачка нефти с грунта» (затраты,

эффективность, предотвращенный ущерб и время выполнения) при

оптимальных параметрах технологического процесса;

Этап 3. Постановка задачи компьютерного моделирования

Логико-математическое описание процесса ПОМ «Откачка нефти с грунта». В качестве основных уровней модели процесса ПОМ «Откачка

11

нефти с грунта», характеризующих состояние ЭЭС, взяты следующие

уровни:

уровень заполнения вакуумного самосвала (далее самосвала) R1t;

уровень заполнения шламового амбара R2t;

уровень загрязнения Vt;

уровень фонда Фt.

Каждому из этих уровней соответствует определенная уровневая переменная (переменная состояния): объем вещества в самосвале, объем

загрязненного грунта нефтью, объем грунта в боне, объем фондов:

1) Объем вещества в самосвале Vсt являются переменной уровня. Он увеличивается в процессе загрузки, уменьшается в процессе разгрузки.

Изменение объема грунта t за момент времени t дается уравнением:

d/dt = PR1t TR1t ,

(3.1)

где t – объем грунта в самосвале, в м.куб.;

PR1t – объем грунта, загружаемого самосвалом за время t;

TR1t – объем грунта разгружаемого самосвалом за время t.

Причем, PR1t = t P1 , где P1 – скорость загрузки самосвала, в м.куб/ч,

3600

TR1t = t P2 , где P2 – скорость разгрузки самосвала, в м.куб/ч

3600

Формула (3.1) описывает различные состояния самосвала (загрузка,

разгрузка, ожидание), при условии:

если TR1t = 0, то модель описывает состояние загрузки самосвала;

если PR1t = 0, то модель описывает состояние разгрузки самосвала;

если PR1t = 0 и TR1t = 0, то модель описывает состояние ожидания;

2) Уровень заполнения шламового амбара Vаt являются

переменной

уровня. Он увеличивается при разгрузке самосвала.

 

Изменение объема грунта в амбаре t дается уравнением:

 

d/dt = TR1t ,

(3.2)

 

12

где t – объем грунта в амбаре, в м.куб.;

Nt – объем загружаемых самосвалов;

TR1t – объем грунта, разгружаемого самосвалом за время t.

3) Уровень загрязнения Vt являются переменной уровня. Он увеличивается за счет вредного воздействия производства, и уменьшается за счет выполнения ПОМ с помощью самосвала.

Изменение объема грунта загрязненного нефтью t, дается уравнением:

d/dt = ZR1t PR1t ,

(3.3)

где t – объем загрязненного грунта на лицензионном участке, в

м.куб.;

PR1t – объем грунта, загружаемого самосвалом за время t;

Nt – объем загружаемых самосвалов за время t;

ZR1t – объем грунта, угнетенного загрязнителем за время t.

Формула (3.3) описывает различные состояния грунта, загрязненного нефтью (восстановление, загрязнение, ожидание), при условии:

если ZR1t=0, то модель описывает состояние удаления грунта;

если PR1t =0, то модель описывает состояние загрязнения грунта;

если PR1t =0 и ZR1t=0, то модель описывает фоновое состояние;

4) Фонды Фt являются переменой уровня.

Они уменьшаются за счет

затрат на использование технического средства самосвал.

Изменение размера Фt в момент времени t дается уравнением:

dФ/dt = Zt,

(3.4)

где Фt объем фондов в руб;

 

 

Zt затраты на использование самосвала за время t.

Величина Zt определяется по формуле:

 

t2

k

 

minZt

Ц i Nt dt ,

(3.5)

t1 i 1

13

Z t
Z t h

где Z(t) – затраты предприятия на выполнение мероприятия за время t;

Цi – цена использования 1 ед.ресурса в час. Если цена использования

ресурса зависит от его состояния, то Ц Ц (st) , где st - состояние ресурса;

Nt – объем используемых ресурсов в момент времени t; k – число типов используемых ресурсов.

Здесь переменные модели PR1t, TR1t и Nt определяют последовательность решений, изменяющих (t-1)-ю систему состояний ЭЭС в

(t)-ю.

Задача оптимизации параметров процесса ПОМ «Откачка нефти с грунта» состоит в определении оптимальных параметров технологического

 

 

 

 

 

 

процесса P ={P1,

P2,

N}

как параметров функции от t в интервале

t1 t t2

(при

заданных

начальных состояниях, ограничениях на

переменные

состояния

и

управления), минимизирующих заданный

«Мультипликативный обобщенный критерий-функционал»:

 

t 2

Z (t) T (t)

 

min МКФ(t)

 

dt .

 

ЭМ (t) ПЭУ (t) PR1(t) TR1(t)

 

t1

 

 

Он сформирован в результате обобщения следующих критериев-

функционалов:

1) Критерии, подлежащие минимизации:

«Затраты на выполнение ПОМ» Z(t) вычисляются по формуле:

Z t h Z t Ц

 

 

h N

,

c

 

 

3600

 

(3.6)

(3.7)

где затраты предприятия на выполнение природоохранного мероприятия в момент времени t+h;

затраты на использование самосвала в момент времени t;

Цс – цена использования 1 ед.самосвала в час. Если цена использования самосвала зависит от его состояния, то Ц с Ц с (st) , где st - состояние самосвала (загрузка, разгрузка, ожидание);

14

N – количество используемых самосвалов в час;

 

h – шаг времени моделирования;

 

«Время выполнения ПОМ» T(t) вычисляется по формуле:

 

T t h T t h ,

(3.8)

где T(t) – время выполнения мероприятия в момент времени t;

Время выполнение мероприятия до полного устранения загрязнения (

t 0 ) на участке равно:

t t2 t1

(0)

 

(0)

,

 

 

P1 N

P2 N

где t0 – начальный объем загрязненного грунта на лицензионном участке, в м.куб.;

P1 и P2 – скорость загрузки и разгрузки самосвала соответственно, в

м.куб/ч;

N – количество самосвалов, единиц в час;

2) Критерии, подлежащие максимизации:

«Эффективность мероприятия» ЭМ(t) вычисляется по формуле:

ЭМ (t h) ЭМ (t)

h P1 N

,

(3.9)

 

3600

 

 

где ЭМ(t) – эффективность мероприятия в момент времени t в м.куб/ч;

P1 – скорость загрузки самосвала, в м.куб/ч;

N – количество самосвалов, единиц в час;

– «Предотвращенный экологический ущерб» ПЭУ(t) вычисляется по формуле:

ПЭУ (t h) ПЭУ (t) Цпэу ЭМ (h) ,

(3.10)

где ПЭУ(t) – предотвращенный экологический ущерб, в руб.;

ЭМ(h) – эффективность мероприятия, объем восстановленного грунта в момент времени h, м.куб.;

Цпэу – стоимость 1 ед. м.куб. восстановленного грунта;

Объем загружаемого грунта вычисляется по формуле:

15

PR1(t h) = PR1(t) h P1 N , 3600

где P1 – скорость загрузки;

Объем разгружаемого грунта вычисляется по формуле:

TR1(t h)=TR1(t)

h P2 N

,

3600

 

 

где P2 – скорость разгрузки;

Количество самосвалов/ часов вычисляется по формуле:

N(t h) = N(t) h N . 3600

Начальные состояния, ограничения на переменные модели:

(3.11)

(3.12)

(3.13)

начальные состояния подсистем ЭЭС: объем загрязненного грунта на лицензионном участке t0; объем грунта в самосвале Vct0; объем грунта в шламовом амбаре t0, объем фонда Фt0, цена использования самосвала за час работы Цс; стоимость 1 ед. восстановленного грунта Цпэу;

целевые состояния подсистем ЭЭС по завершению ТП ПОМ: tH=0,

tH =0 и tH t0;

ограничения на переменные состояния ЭЭС:

Vci t Vci t max ,

(3.14)

где Vct – объем грунта в самосвале;

 

Vci t max – предельно-допустимый объем грунта в самосвале;

 

i – тип самосвала, для каждого типа самосвала свой объем кузова;

 

t t max ,

(3.15)

где t – объем грунта в амбаре;

 

t max – предельно-допустимый объем бона;

 

ограничения на режимы функционирования самосвала:

 

P1i min P1i P1i max ,

(3.16)

где P1i – скорость загрузки самосвала типа i;

 

 

16

P1min и P1max – минимально и максимально допустимые скорости загрузки самосвала типа i соответственно;

P2i min P2i P2i max ,

(3.17)

где P2i – скорость разгрузки самосвала;

P2i min и P2i max – минимально и максимально допустимые скорости

разгрузки самосвала соответственно;

 

ограничения на число доступных самосвалов в ТП ПОМ:

 

N min N N max ,

(3.18)

где N – число самосвалов в час,

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

В целом, состояние эколого-экономической системы в процессе выполнения природоохранного мероприятия «Откачка нефти с грунта» характеризуется приведенными выше дифференциальными уравнениями

(3.2)-(3.5). Динамика вычисляемых эколого–экономических показателей данного природоохранного мероприятия определяется дифференциальными уравнениями (3.2)-(3.13).

Предложенное выше логико-математическое описание природоохранного мероприятия «Откачка нефти с грунта» требует учета при синтезе моделей компонентов и моделирующего алгоритма с целью автоматизированного решения поставленной задачи оптимального управления ЭЭС (3.6), (3.14)-(3.18);

Этап 4. Разработка моделей компонентов ЭЭС в формате МКЦ

В процессе выполнения ПОМ «Откачка нефти с грунта» задействованы следующие компоненты ЭЭС:

компонент природной среды Грунт;

техническое средство Самосвал;

компонент природной среды Шламовый амбар;

компонент Фонд;

17

конвертеры, в т.ч.:

затраты на выполнение ПОМ;

эффективность мероприятия;

предотвращенный экологический ущерб;

время выполнения ПОМ;

компенсирующее воздействие на КПС Грунт;

управляемый ключ и условие, подающее ему сигнал.

Разработаны условные обозначения и математические модели компонентов ЭЭС (табл. 3.1) и конвертеров ПОМ «Откачка нефти с грунта»

(табл. 3.2), в их 3-х аспектах: математическом, топологическом и физическом. Моделей компонентов и конвертеров выполнены с учетом выше предложенного логико-математического описания ПОМ «Откачка нефти с грунта» (3.1)-(3.16).

Табл. 3.1.

Модели компонентов ЭЭС в процессе ПОМ «Откачка нефти с грунта»

 

Название и

 

Математическая

Параметры

изображение

 

модель

 

 

 

 

 

 

 

 

 

 

1.

 

 

 

 

 

VN1 – объем загрузки

 

 

 

 

 

VN2 – объем разгрузки

 

 

 

 

 

 

 

 

 

 

 

 

VN4 – объем заполнения самосвала

 

 

 

 

 

Модель загрузка:

N3 – связь для параметризации

 

Самосвал

 

VN4 t h VN4(t) VN1(h)

N5 – связь для обмена информации и

 

 

 

 

 

Модель разгрузка:

передачи параметров

 

 

N3

 

 

 

 

 

 

 

 

N1

 

 

N4

VN4 t h VN4(t) VN1(h)

Параметры модели самосвал,

 

 

 

 

 

 

 

 

 

 

Самосвал

 

 

 

 

N2

N5

 

загружаются из БД:

 

 

 

 

 

Модель ожидание:

P1 – скорость загрузки

 

 

 

 

 

VN4 t h VN4(t)

 

 

 

 

 

P2 – скорость разгрузки

 

 

 

 

 

 

 

 

 

 

 

 

Vmax – вместимость кузова

 

 

 

 

 

 

N – количество используемых

 

 

 

 

 

 

самосвалов данного типа в час

 

 

 

 

 

 

 

 

 

 

 

 

 

18

2.

Грунт

 

 

 

 

VN1(h) – скорость загрузки,

 

 

 

 

 

N2

 

 

 

 

компенсирующее воздействие

 

 

 

 

 

 

 

 

 

 

 

N1

 

 

N3

 

VN2 – начальный объем загрязнения,

 

Грунт

 

VN3 t h VN3(t) VN1(h)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

параметр загружается из БД

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VN3(t+h) – объем загрязненного грунта

 

 

 

 

 

 

 

 

 

 

в момент времени t+h

 

 

 

 

 

 

 

 

 

 

 

3.

 

 

 

 

 

 

 

 

 

VN1(h) – объем разгрузки,

 

 

 

 

 

 

 

 

 

 

 

Шламовый амбар

 

компенсирующее воздействие

 

 

 

 

 

 

 

 

 

 

 

 

VN3 t h VN3(t) VN1(h)

VN2 – начальный объем заполнения,

 

 

 

 

 

N2

 

 

 

параметр из БД

 

 

 

 

 

 

 

 

 

 

N1

 

 

 

 

 

 

N3

 

 

 

 

 

 

Амбар

 

 

 

 

VN3(t+h) – объем загрязненного грунта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в амбаре в момент времени t+h

 

 

 

 

 

 

 

 

 

 

 

4.

Фонд

 

 

 

 

VN1(h) – затраты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N2

 

 

 

VN3 t h VN3(t) VN1(h)

VN2 – начальный объем фонда,

 

N1

 

 

 

N3

 

параметр загружается из БД

 

 

 

Фонд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VN3(t+h) – объем фонда в момент t+h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.

Ресурс

 

 

 

 

VN1 – потребление ресурса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VN3 VN2 VN1

VN2 – начальный объем ресурса,

 

 

 

 

 

N2

 

 

 

параметр загружается из БД

 

 

 

 

 

 

 

 

 

 

N1

 

 

 

N3

 

 

 

 

 

 

 

 

Ресурс

 

 

VN3 – объем свободного ресурса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19

6.Переключатель режима модели Самосвал

 

UK

 

N1

U

 

 

 

 

N2

 

N3

К1

N5

 

N4

К2

N6

Если VN2t =[1,0],

то режим загрузки:

VN5t = VN3t,

VN6t = 0.

Если VN2t =[0,1],

то режим разгрузки:

VN5t = 0,

VN6t = VN4t.

Если VN2t =[0,0],

то режим ожидания:

VN5t = 0,

VN6t = 0.

U – условия переключения режима работы модели. Если для этого требуется время, то используется переключение с задержкой

VN1 – уровень наполнения самосвала

N2 – сигнальный выход

VN2 – вектор-сигнал для ключей К1,

К2, которые зависят от величины VN1 N4 – связь с КПС Грунт

N5 – связь со шламовым амбаром для локализации загрязнения

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]