392249ca-b931-434c-8215-e418712c00a3 168 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

Объединение матриц и построение поля риска

N_dist 370 :

Расстояние для расчета риска от центра оборудования, м

N_diap 8 :

На сколько интервалов разбиваем диапазон значений риска

N_delta 3 :

Шаг расчета (целые числа), м

N_dist N_dist N_delta / 0 round N_delta * :

N_dist 369

Корректировка расстояния расчета

field 13 matrix :

field 000014 mat :

Объявляем функцию для поиска в матрице значений

matr dist find_value matr dist N_delta / 0 round 1 + 2 el :

Добавление матрицы в поле

X_eq 100 :

Y_eq 100 :

Координаты оборудования (или начало трубопровода)

X_eq_end 200 :

Y_eq_end 200 :

Конец трубопровода

Построение матрицы координат нитки трубопровода

graf_pipe 12 matrix :

graf_pipe 0012 mat :

ОбЪявление двумерной матрицы

delt_X X_eq_end X_eq - abs :

delt_Y Y_eq_end Y_eq - abs :

long delt_X 2^delt_Y 2^+ sqrt :

delt_X N_delta / delt_Y N_delta / ≥ N_dots delt_X N_delta / 0 round : N_dots delt_Y N_delta / 0 round : if

N_dots 33

count_X 1 :

X X_eq : count_X N_dots ≤ X X X_eq_end X_eq - N_dots / 0 round + : graf_pipe count_X 1 el X : count_X count_X 1 + : 21 line for

graf_pipe 10001030106010901120115011801210124012701300133013601390142014501480151015401570160016301660232 mat

count_Y 1 :

Y Y_eq : count_Y N_dots ≤ Y Y Y_eq_end Y_eq - N_dots / 0 round + : graf_pipe count_Y 2 el X : count_Y count_Y 1 + : 21 line for

Переносим сюда матрицу риска от аварии на одном оборудовании

transfer 0.0049 0.0042 0.0032 0.0013 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.0049 0.0041 0.0021 0.0013 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.0049 0.004 0.002 0.0012 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.0049 0.0038 0.0019 0.0012 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0049 0.0037 0.0018 0.0012 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0049 0.0036 0.0017 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0049 0.0036 0.0016 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0047 0.0035 0.0016 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0046 0.0034 0.0015 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0045 0.0033 0.0014 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0044 0.0032 0.0014 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 01112 mat :

graf_all_all 12 matrix :

graf_all_all 0012 mat

ОбЪявление двумерной матрицы

Преобразуем ее в двумерную матрицу (исходную матрицу)

n 1 :

k 1 :

X 1 : X transfer rows transfer cols * (≤ X X 1 + : graf_all_all X 1 el X N_delta * N_delta - : graf_all_all X 2 el transfer n k el : n transfer rows / 1 - 0 ≡ n 0 : k k 1 + : 21 line n n : if n n 1 + : 41 line for

Получаем переменную Val с значением риска на заданном расстоянии dist

Строим матрицу в формате (координата Х, координата У, значение риска) для дальнейших выборок и суммирования с другими матрицами. Исходим, что поле риска вокруг оборудования круговое.

Здесь будут перечислены все точки поля, которое имеет размер N_dist по вертикали и по горизонтали

prom X_eq 3 nthroot Y_eq 4 nthroot * :

Идентификатор точки

count 1 :

X 0 : X N_dist ≤ (X X N_delta + : Y 0 : Y N_dist ≤ (Y Y N_delta + : field count 1 el X : field count 2 el Y : 11 line dist_eq_to_dot X X_eq - (2^Y Y_eq - (2^+ sqrt : dist_eq_to_dot N_dist > dist_eq_to_dot N_dist : dist_eq_to_dot dist_eq_to_dot : if field count 4 el prom ≠ field count 3 el graf_all_all dist_eq_to_dot find_value field count 3 el + : field count 4 el prom : 21 line prom prom : if count count 1 + : 21 line 11 line 11 line 11 line 11 line 51 line for for

Здесь закончили добавление риска от одного оборудования к общей матрице

Добавление матрицы в поле

X_eq 250 :

Y_eq 250 :

Координаты оборудования

Переносим сюда матрицу риска от аварии на одном оборудовании

transfer 0.0049 0.0042 0.0032 0.0013 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.0049 0.0041 0.0021 0.0013 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.0049 0.004 0.002 0.0012 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.0049 0.0038 0.0019 0.0012 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0049 0.0037 0.0018 0.0012 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0049 0.0036 0.0017 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0049 0.0036 0.0016 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0047 0.0035 0.0016 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0046 0.0034 0.0015 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0045 0.0033 0.0014 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0 0.0044 0.0032 0.0014 0.0011 0.001 0.001 0.001 0.001 0.001 0.001 0.001 01112 mat :

graf_all_all 12 matrix :

graf_all_all 0012 mat

ОбЪявление двумерной матрицы

Преобразуем ее в двумерную матрицу (исходную матрицу)

n 1 :

k 1 :

X 1 : X transfer rows transfer cols * (≤ X X 1 + : graf_all_all X 1 el X N_delta * N_delta - : graf_all_all X 2 el transfer n k el : n transfer rows / 1 - 0 ≡ n 0 : k k 1 + : 21 line n n : if n n 1 + : 41 line for

Получаем переменную Val с значением риска на заданном расстоянии dist

Здесь будут перечислены все точки поля, которое имеет размер N_dist по вертикали и по горизонтали

prom X_eq 3 nthroot Y_eq 4 nthroot * :

Идентификатор точки

count 1 :

X 0 : X N_dist ≤ (X X N_delta + : Y 0 : Y N_dist ≤ (Y Y N_delta + : field count 1 el X : field count 2 el Y : 11 line dist_eq_to_dot X X_eq - (2^Y Y_eq - (2^+ sqrt : dist_eq_to_dot N_dist > dist_eq_to_dot N_dist : dist_eq_to_dot dist_eq_to_dot : if field count 4 el prom ≠ field count 3 el graf_all_all dist_eq_to_dot find_value field count 3 el + : field count 4 el prom : 21 line prom prom : if count count 1 + : 21 line 11 line 11 line 11 line 11 line 51 line for for

Здесь закончили добавление риска от одного оборудования к общей матрице

Если объект является линейным, то этот цикл нужно повторить для точек вдоль линии

Построение поля риска

Количество линий (областей)

Размер точки

N_line 10 :

S_dot 30 :

field_sort field 3 csort :

k field_sort rows :

max_risk field_sort k 3 el :

min_risk 0 :

n 1 : n k ≤ (min_risk 0 ≡ (& n n 1 + : min_risk field_sort n 3 el : for

min_risk field 3 col min :

max_risk field 3 col max :

min_risk 0.002

max_risk 0.0059

wid max_risk min_risk - N_line / :

Ширина одного диапазона значений для построения линии

Минимальные и максимальные границы диапазона записываем в матрицу. Вверху верхняя граница, внизу - нижняя граница диапазона

diap 210 matrix :

diap 00000000000000000000210 mat

x 1 : x N_line ≤ x x 1 + : diap 2 x el max_risk x 1 - (wid * - : diap 1 x el max_risk x (wid * - : 21 line for

diap 0.0055 0.0051 0.0047 0.0043 0.004 0.0036 0.0032 0.0028 0.0024 0.002 0.0059 0.0055 0.0051 0.0047 0.0043 0.004 0.0036 0.0032 0.0028 0.0024 210 mat

sort 12 matrix :

sort 0012 mat

sort_1 12 matrix :

sort_1 0012 mat

sort_2 12 matrix :

sort_2 0012 mat

sort_3 12 matrix :

sort_3 0012 mat

sort_4 12 matrix :

sort_4 0012 mat

sort_5 12 matrix :

sort_5 0012 mat

sort_6 12 matrix :

sort_6 0012 mat

sort_7 12 matrix :

sort_7 0012 mat

sort_8 12 matrix :

sort_8 0012 mat

sort_9 12 matrix :

sort_9 0012 mat

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

Получаем матрицу в формате (координата Х, координата У) для построения графика. Это координаты соответствуют выбранному диапазону значений риска.

n 0 :

mb 0 :

m 0 :

mc 0 :

ma 0 :

md 0 :

me 0 :

mh 0 :

mf 0 :

mg 0 :

count 1 :

X 0 : X N_dist ≤ (X X N_delta + : Y 0 : Y N_dist ≤ (Y Y N_delta + : field count 3 el diap 11 el ≥ (field count 3 el diap 21 el ≤ (& (sort n 1 + 1 el X : sort n 1 + 2 el Y : sort n 1 + 3 el . : 31 line sort n 1 + 4 el S_dot : sort n 1 + 5 el red : n n 1 + : 41 line n n : if field count 3 el diap 12 el ≥ (field count 3 el diap 22 el ≤ (& (sort_1 m 1 + 1 el X : sort_1 m 1 + 2 el Y : sort_1 m 1 + 3 el . : 31 line sort_1 m 1 + 4 el S_dot : sort_1 m 1 + 5 el purple : m m 1 + : 41 line n n : if field count 3 el diap 13 el ≥ (field count 3 el diap 23 el ≤ (& (sort_2 ma 1 + 1 el X : sort_2 ma 1 + 2 el Y : sort_2 ma 1 + 3 el . : 31 line sort_2 ma 1 + 4 el S_dot : sort_2 ma 1 + 5 el fuchsia : ma ma 1 + : 41 line n n : if 21 line field count 3 el diap 14 el ≥ (field count 3 el diap 24 el ≤ (& (sort_3 mb 1 + 1 el X : sort_3 mb 1 + 2 el Y : mb mb 1 + : 31 line n n : if field count 3 el diap 15 el ≥ (field count 3 el diap 25 el ≤ (& (sort_4 mc 1 + 1 el X : sort_4 mc 1 + 2 el Y : mc mc 1 + : 31 line n n : if field count 3 el diap 16 el ≥ (field count 3 el diap 26 el ≤ (& (sort_5 md 1 + 1 el X : sort_5 md 1 + 2 el Y : md md 1 + : 31 line n n : if field count 3 el diap 17 el ≥ (field count 3 el diap 27 el ≤ (& (sort_6 me 1 + 1 el X : sort_6 me 1 + 2 el Y : me me 1 + : 31 line n n : if field count 3 el diap 18 el ≥ (field count 3 el diap 28 el ≤ (& (sort_7 mf 1 + 1 el X : sort_7 mf 1 + 2 el Y : mf mf 1 + : 31 line n n : if field count 3 el diap 19 el ≥ (field count 3 el diap 29 el ≤ (& (sort_8 mg 1 + 1 el X : sort_8 mg 1 + 2 el Y : mg mg 1 + : 31 line n n : if field count 3 el diap 110 el ≥ (field count 3 el diap 210 el ≤ (& (sort_9 mh 1 + 1 el X : sort_9 mh 1 + 2 el Y : mh mh 1 + : 31 line n n : if count count 1 + : 91 line 21 line for for

count 1 : count field rows ≤ count count 1 + : field count 3 el diap 11 el ≥ (field count 3 el diap 21 el ≤ (& (sort n 1 + 1 el field count 1 el : sort n 1 + 2 el field count 2 el : sort n 1 + 3 el . : 31 line sort n 1 + 4 el S_dot : sort n 1 + 5 el red : n n 1 + : 41 line n n : if field count 3 el diap 12 el ≥ (field count 3 el diap 22 el ≤ (& (sort_1 m 1 + 1 el field count 1 el : sort_1 m 1 + 2 el field count 2 el : sort_1 m 1 + 3 el . : 31 line sort_1 m 1 + 4 el S_dot : sort_1 m 1 + 5 el purple : m m 1 + : 41 line n n : if field count 3 el diap 13 el ≥ (field count 3 el diap 23 el ≤ (& (sort_2 ma 1 + 1 el field count 1 el : sort_2 ma 1 + 2 el field count 2 el : sort_2 ma 1 + 3 el . : 31 line sort_2 ma 1 + 4 el S_dot : sort_2 ma 1 + 5 el fuchsia : ma ma 1 + : 41 line n n : if field count 3 el diap 14 el ≥ (field count 3 el diap 24 el ≤ (& (sort_3 mb 1 + 1 el field count 1 el : sort_3 mb 1 + 2 el field count 2 el : mb mb 1 + : 31 line n n : if field count 3 el diap 15 el ≥ (field count 3 el diap 25 el ≤ (& (sort_4 mc 1 + 1 el field count 1 el : sort_4 mc 1 + 2 el field count 2 el : mc mc 1 + : 31 line n n : if field count 3 el diap 16 el ≥ (field count 3 el diap 26 el ≤ (& (sort_5 md 1 + 1 el field count 1 el : sort_5 md 1 + 2 el field count 2 el : md md 1 + : 31 line n n : if field count 3 el diap 17 el ≥ (field count 3 el diap 27 el ≤ (& (sort_6 me 1 + 1 el field count 1 el : sort_6 me 1 + 2 el field count 2 el : me me 1 + : 31 line n n : if field count 3 el diap 18 el ≥ (field count 3 el diap 28 el ≤ (& (sort_7 mf 1 + 1 el field count 1 el : sort_7 mf 1 + 2 el field count 2 el : mf mf 1 + : 31 line n n : if field count 3 el diap 19 el ≥ (field count 3 el diap 29 el ≤ (& (sort_8 mg 1 + 1 el field count 1 el : sort_8 mg 1 + 2 el field count 2 el : mg mg 1 + : 31 line n n : if field count 3 el diap 110 el ≥ (field count 3 el diap 210 el ≤ (& (sort_9 mh 1 + 1 el field count 1 el : sort_9 mh 1 + 2 el field count 2 el : mh mh 1 + : 31 line n n : if count count 1 + : 111 line for

rr field rows :

count 1 : count rr 1 - ≤ (count count 1 + : field count 3 el diap 11 el ≥ (field count 3 el diap 21 el ≤ (& (sort n 1 + 1 el field count 1 el : sort n 1 + 2 el field count 2 el : sort n 1 + 3 el . : 31 line sort n 1 + 4 el S_dot : sort n 1 + 5 el red : n n 1 + : 41 line n n : if for

count 1 : count rr 1 - ≤ (count count 1 + : field count 3 el diap 12 el ≥ (field count 3 el diap 22 el ≤ (& (sort_1 m 1 + 1 el field count 1 el : sort_1 m 1 + 2 el field count 2 el : sort_1 m 1 + 3 el . : 31 line sort_1 m 1 + 4 el S_dot : sort_1 m 1 + 5 el purple : m m 1 + : 41 line n n : if for

sort sort_1 sort_2 sort_3 sort_4 sort_5 sort_6 sort_7 sort_8 sort_9 10 1 sys