ОПТИМИЗАЦИЯ ПРОЦЕССА ПОДОГРЕВА НЕФТЕГАЗОКОНДЕНСАТНОЙ СМЕСИ ПАРАМИ НАФТЫ В ТЕПЛООБМЕННИКЕ-КОНДЕНСАТОРЕ 10Е-03 УСТАНОВКИ ПЕРВИЧНОЙ ПЕРЕГОНКИ НЕФТИ

OPTIMIZATION OF THE PROCESS OF HEATING AN OIL AND GAS CONDENSATE MIXTURE BY NAPHTHA VAPOR IN THE 10E-03 HEAT EXCHANGER-CONDENSER OF THE INSTALLATION OF THE PRIMARY OIL DISTILLATION
Цитировать:
ОПТИМИЗАЦИЯ ПРОЦЕССА ПОДОГРЕВА НЕФТЕГАЗОКОНДЕНСАТНОЙ СМЕСИ ПАРАМИ НАФТЫ В ТЕПЛООБМЕННИКЕ-КОНДЕНСАТОРЕ 10Е-03 УСТАНОВКИ ПЕРВИЧНОЙ ПЕРЕГОНКИ НЕФТИ // Universum: технические науки : электрон. научн. журн. Худайбердиев А.А. [и др.]. 2022. 4(97). URL: https://7universum.com/ru/tech/archive/item/13532 (дата обращения: 15.08.2022).
Прочитать статью:
DOI - 10.32743/UniTech.2022.97.4.13532

 

АННОТАЦИЯ

Сформулирована целевая функция критерия оптимальности - удельной себестоимости подогреваемой нефтегазоконденсатной смеси 30%Н+70%ГК, включающая в себе изменения расхода электроэнергии и амортизационных отчислений в зависимости от производительности, условий работы теплообменника и расходуемой мощности для перекачки смеси. Построены кривые зависимостей удельной себестоимости подогретой смеси от поверхности нагрева аппарата и температуры смеси на выходе. Выявлены оптимальные условия эксплуатации теплообменника 10Е-03 Бухарского НПЗ: F = 445 м2, А = 245,8 сум/кг, Э = 250,5 сум/кг, Суд = 508,5 сум/кг и t = 77,0 оС при заданной производительности Gо = 105508,3 кг/ч.

ABSTRACT

The objective function of the optimality criterion is formulated - the unit cost of a heated oil and gas condensate mixture of 30%Oil+70%GC, which includes changes in electricity consumption and depreciation charges depending on the performance, operating conditions of the heat exchanger and the consumed power for pumping the mixture. Curves of the dependences of the unit cost of the heated mixture on the heating surface of the apparatus and the temperature of the mixture at the outlet are constructed. Optimal operating conditions of the heat exchanger 10E-03 of the Bukhara refinery have been identified: F = 445 m2, A = 245.8 sum/kg, E = 250.5 sum/kg, C = 508.5 sum/kg and t = 77.0 oС at a given productivity Go =105508.3 kg/h.

 

Ключевые слова: нефтегазоконденсатная смесь, перегонка, тяжелая нафта, нагревание, теплообменник, конденсатор, поверхность теплообмена, оптимизация, технологическая себестоимость.

Keywords: oil and gas condensate mixture, distillation, heavy naphtha, heating, heat exchanger, condenser, heat exchange surface, optimization, technological cost.

 

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

Как известно, критерий оптимальности является мерой количественной оценки оптимизируемого качества объекта - процесса или аппарата [3]. Выражение критерия оптимальности в виде экономической оценки, таких как производительность, себестоимость продукции, прибыль, рентабельность и пр., являются общей постановкой задачи оптимизации. В случаях, когда объектом оптимизации является часть технологического процесса, критерием оптималь-ности может служить технологическая характеристика, косвенно оценивающая экономичность работы аппарата - выход продукта, продолжительность процесса, расходы теплоносителей и их температуры на входе и выходе из теплообменников, разность температур между теплоносителями, тепловые коэффициенты и тепловая нагрузка аппарата. При этом в качестве варьируемых параметров рассматриваются расход сырья, температуры потоков теплоносителей на выходе и передаваемое аппаратом количество теплоты.

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

С учетом производительности теплообменника G и ограничения на технологические показатели сырья и греющего теплоносителя, критерий оптимальности в общем виде можно выразить в виде [4]:

R = f (G, Kс, Kз, Эз),                                                                (1)

где Кс - технологические показатели сырья; Кз и Эз - удельные капитальные и эксплуатационные затраты для осуществления процесса, отнесенные к одному году нормативного срока окупаемости Тн.

Входящие в (1) экономические величины Кз и Эз придают универсаль-ность данному критерию и позволяют применять для оптимизации конструкций аппаратов, независимо от их назначения и особенностей.

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

Стех = СоGо + СтGт + Сэ(Nн+Nд) + Аа Fа+ Акн Fкн + Ан(Nн+Nд),                  (2)

где Со, Ст и Сэ - соответственно, себестоимость сырья, греющего теплоносителя и электроэнергии; Gо и Gт - расходы сырья и теплоносителя; Nн и Nд - мощность насоса для перекачки нефти и дистиллятов фракций; Fт и Fкн - поверхность теплопередачи теплообменников и конденсаторов; Аа и Ан - амортизационые отчисления для технологических аппаратов и насосов.

Отметим, что нефть в ходе подогрева не подвергается технологической обработке. Поэтому стоимость сырой нефти Со зависит от ее качества и не зависит от режима работы теплообменного оборудования. Необходимо также учесть, что отходящие из ректификационной колонны потоки дистиллятов топливных фракций и кубовой остаток подлежат охлаждению до температуры их хранения в парках НПЗ [1,6]. Исходя из этого, для повышения тепловой эффективности нефтеперегонной установки, эти горячие потоки используются для последовательного многоступенчатого предварительного подогрева нефти, подаваемую в змеевиковую печь. Поэтому затраты, связанные с использова-нием горячих потоков, не оказывают влияния на технологическую себестои-мость подогретой нефти в теплообменных аппаратах. Помимо этого, заработная плата персонала по техническому обслуживанию аппаратов также не зависит от интенсивности эксплуатации оборудования. В силу вышеописанных обстоя-тельств, затраты, связанные с покупкой нефти, теплоносителей и зарплата технического персонала в состав выражения критерия оптимальности исследуемого процесса (2) не включаются [4]:

Стех = Сэ(Nн+Nд) + Аа Fа+ Акн Fкн + Ан(Nн+Nд).                             (3)

Как известно, трубчатые теплообменные аппараты трех блоков предва-рительного подогрева сырья нефтеперегонной установки имеют различные конструкции и производительность [1,2,4,6,7]. По этой причине, для выявления оптимального состава блоков теплообменников нефтеперегонной установки и разработка ее энергосберегающей технологической схемы в качестве критерия оптимальности целесообразно принять удельную технологическую себестои-мость нагретого сырья Суд = Стех/Gо. В этом случае (3) можно выразить в виде:

Суд = 1/Gо[Сэ(Nн+Nд) + АаFа + АкнFкн + Ан (Nн + Nд)].                   (4)

Сопоставительная оценка влияния статьи затрат на технологическую себестоимость подогретой нефти проводится путем анализа уравнений для расчета параметров, входящих в выражение целевой функции критерия оптимальности (4).

Мощность насоса N (кВт) для перекачки технологических потоков (нефти и дистиллятов фракций) по трубкам теплообменных аппаратов можно определить по известному выражению [8]:

N = (Gо .DP)/(1000 r hн),                                                         (5)

где Gо - массовый расход потока, кг/с; DP - гидравлическое сопротивление тракта перекачки потока, Па; r - плотность потока, кг/м3; hн - к.п.д. насоса.

Величина потери давления DP для преодоления сил внутреннего трения в теплообменных трубках аппарата определяется по известной формуле [8]:

DP = 0,5u2r (lLобщ/dэквji),                                                 (6)

где u = G/(0,785dвн2r) - скорость потока в трубках аппаратов, м/с; dвн - внутренний диаметр трубок, м; l = f(Re) - коэффициент трения, определяемый в зависимости от режима движения сырья в трубках по числу Re; Re = (udвнr)/µ - число Рейнольдса; µ - динамический коэффициент вязкости сырья, Па.с; Lобщ = n.l - общая длина трубок, м; n - число труб в аппарате, шт.; l - рабочая длина одной трубы, м; åji - суммарный коэффициент местных сопротивлений.

С учетом производительности теплообменных аппаратов по сырью Gо их теплопередающая поверхность Fа определяется по выражению

Fа = Q/(К Dtср) = Gо (cвых tвых - cвх tвх)/(К Dtср),                          (7)

где Q = Gо(cвых tвых - cвх tвх) - тепловая нагрузка аппарата, Вт; cвх и cвых - теплоемкость сырья при температурах его входа в аппарат tвх и на выходе из него tвых, Дж/(кг.оС); К - коэффициент теплопередачи в аппарате, Вт/(м2 оС); Dtср - полезная разность температур, оС.

Теплоемкость нефтяного сырья c (кДж/кг.оС), с учетом его температуры Т и относительной плотности , определяется по формуле [1,2,4]:

.                                              (8)

Коэффициенты теплоотдачи от греющего теплоносителя к стенке труб a1 и от стенки к подогреваемой жидкости a2, а также коэффициент теплопередачи К в теплообменных аппаратах рассчитывается по уточненной методике [9], с учетом условий работы, используя температурные изменения показателей свойств сырья - плотность ρ, вязкость ν, и μ, теплоемкость с, теплопроводность λ и др.

Отметим, что в настоящее время на нефтеперерабатывающих заводах (НПЗ) амортизационные отчисления Аа принимаются как условно-постоянная величина от стоимости аппаратов Ца. В действительности же показатель Аа является переменной величиной и зависит от интенсивности работы теплообменных аппаратов Т [4]:

Аа = (ЕнЦа)/24Т Fа = (EнЦа)/24Тн [Gо (cвыхtвых - cвхtвх)/К.Dtср],                         (9)

где Ен = 0,15 - нормативный коэффициент эффективности капитальных вложений в отрасль промышленности; Цт - стоимость аппарата, сум.

Аналогичным образом, амортизационные отчисления для насосов Ан [4]

Ан = (ЕнЦн)/24 Тн  N = (EнЦт)/24 Тн [(G .DP)/(1000 r hн)],                          (10)

где Цн - стоимость насоса, сум.

Ограничения в области исследования целевой функции критерия оптимальности устанавливаются по температуре подогретого сырья на выходе блока теплообменников tогр (tвых£ 220÷240 °C) и специфики работы змеевиковой печи (минимальная температура сырья на входе в печь tмин £ 120÷150 °C) [1,2,6].

Таким образом, целевая функция критерия оптимальности подогрева нефтяного сырья углеводородными теплоносителями в кожухотрубчатом теплообменном аппарате можно сформулировать как систему уравнений:

Решение системы уравнений (11) сводится к выявлению оптимальных условий эксплуатации теплообменников блока подогрева нефти, обеспечива-ющие минимальную технологическую себестоимость подогрева сырья.

Результаты. Целевая функция критерия оптимальности подогрева нефти (11) исследовалась применительно к условиям эксплуатации промышленного горизонтального кожухотрубчатого теплообменника-конденсатора 10Е-03 пер-вого этапа предварительного подогрева сырья Бухарского НПЗ. В данном аппарате осуществляется подогрев технологической смеси, состоящей из 30 % нефти и 70 % газового конденсата (30%Н+70%ГК) теплом конденсирующихся паров тяжелой нафты, отводимых с верхней части атмосферной ректифика-ционной колонны нефтеперегонной установки [10].

Теплообменник 10Е-03 имеет следующие конструктивные параметры: d = 20/25 мм, l = 6000 мм, n = 1106 шт. Поверхность теплопередачи аппарата по величине dср равняется F = 468,3 м2. Стоимость аппарата по данным завода составляет Цн = 462842 тыс. сум. Установленная мощность насоса для перекачки нефтегазоконденсатной смеси через восемь последовательно соединенных теплообменников блока равна Nн = 250 кВт, стоимость насоса составляет Цн = 224432 тыс. сум. При этом величина амортизационных отчислений составляют, соответственно:

- для теплообменника 10Е-03 по (9):

Ат = (Ен.Цт)/(24.Тн.Fт) = (0,14.462842.103)/(24.340.468,3) = 16,97 сум/м2.

- для насоса по (10):

Ан = (EнЦн)/24ТнNн = (0,14. 224432.103)/(24.340.250) = 15,4 сум/кВт,

Нефтегазоконденсатная смесь поступает в аппарат при температуре tвх = 49 оС, где она подогревается до tвых = 96 оС конденсирующимися парами тяжелой нафты при Т1 = 165 оС. Коэффициент теплопередачи К (Вт/м2 К) в аппарате определялся по методике [9], с использованием уточненных значений температурных изменений показателей свойств сырья и теплоносителя [1,2,4,11-13].

Исследование проведены при следующих технологических параметрах аппарата: производительность смеси Gо = 105508,3 кг/ч и расход паров нафты Gд = 61898 кг/ч, варьируя температурой подогрева смеси в диапазоне 49÷96 оС.

Для данного варианта исследований вид целевой функции (4) упрощается и принимает вид

Суд = 1/Gо [СэNн + АаFа + АнNн].                                               (12)

Результаты расчета целевой функция критерия оптимальности процесса в аппарате 10Е-03 (11) при значениях Gо , Gд, tвх, tвых и Т1 приведены в таблице 1. 

Таблица 1.

Результаты расчета целевой функция критерия оптимальности подогрева смеси 30%Н+70%ГК парами тяжелой нафты в теплообменнике 10Е-03 при Gо = 105508 кг/ч, Gд= 61898 кг/ч и tкн = 165 о

п/п

tтн ,

оС

Δtср,

оС

K , Вт/(мС)

F,

м2

Ат F., сум/кг

СэNн+AнNн, сум/кг

Cуд , сум/кг

1

49,0

82,48

131,94

390,68

214,619

250,456

476,67

2

52,5

81,08

132,02

397,20

218,201

250,456

480,44

3

56,0

82,50

132,09

390,16

214,332

250,456

476,36

4

59,5

80,75

132,17

398,39

218,855

250,456

481,13

5

63,0

79,00

132,24

406,99

223,580

250,456

486,1

6

66,5

77,25

132,31

415,99

228,522

250,456

491,32

7

70,0

75,50

132,38

425,40

233,694

250,456

496,77

8

73,5

73,75

132,45

435,27

239,114

250,456

502,49

9

77,0

72,00

132,52

445,62

244,799

250,456

508,48

10

80,5

70,25

132,59

456,48

250,769

250,456

514,77

11

84,0

68,50

132,65

467,91

257,046

250,456

521,39

12

87,5

66,75

132,72

479,94

263,655

250,456

528,35

13

91,0

65,00

132,78

492,62

270,622

250,456

535,70

14

94,5

63,25

132,85

506,01

277,976

250,456

543,45

 

По данным табл. 1 построены кривые изменения удельной технологи-ческой себестоимости Суд подогрева нефтегазоконденсатной смеси в диапазоне повышения ее температуры 49÷96 оC при постоянной температуре конденса-ции паров тяжелой нафты tкн = 165 оС (рис. 1).

 

Рисунок 1. Изменение удельной себестоимости Суд подогрева нефтегазоконденсатной смеси парами тяжелой нафты в теплообменнике 10Е-03 от повышения температуры смеси t при Gо= 105508 кг/ч, Gд= 61898 кг/ч и tкн = 165 оС

 

Как видно из рисунка 1, с повышением температуры смеси удельная технологическая себестоимость ее подогрева Суд в аппарате интенсивно растет от 476,67 до 543,45 сум/кг с криволинейной закономерностью. Скачок вели-чины Суд в начальной участке кривой поясняется изменением гидродинами-ческого режима потока в аппарате от значения перепада температуры Δtср между теплоносителями.

По данным таблицы 1 также построена зависимость изменения удельной технологической себестоимости Суд подогрева смеси от величины поверхности нагрева теплообменника Fа (рис. 2). Как видно из рисунка 2, с увеличением потребной поверхности нагрева аппарата Fа в пределах от 391 до 506 м2 удельная себестоимость подогрева смеси Суд в нем интенсивно растет по наклонной кривой от 476,7 до 543,5 сум/кг.

 

Рисунок 2. Изменение удельной себестоимости Суд подогрева нефтегазоконденсатной смеси парами тяжелой нафты от величины поверхности нагрева Fа теплообменника 10Е-03 при Gо= 105508 кг/ч, Gд= 61898 кг/ч, tвх = 49 оС и tвых = 96 оС

 

Характер изменений составляющих удельную себестоимость подогрева рабочей смеси Суд в теплообменнике - энергетических затрат для осуществле-ния процесса Э = СэNн+AнNн и амортизационных отчислений на оборудования А = АаFа отражено на рисунке 3. Анализ значений этих составляющих себестои- мости подогрева смеси сводится к следующему. Для перекачки смеси по трубкам аппарата в количестве Gо = 105508,3 кг/ч потребуется Nн = 16,1 кВт мощности. При себестоимости электроэнергии Сэ = 440,52 сум/кВт для промышленных предприятий, энергетические затраты для осуществления процесса подогрева смеси в аппарате составляют постоянную величину Э = 250,456 сум/кг (кривая 2 рис. 3).

 

Рисунок 3. Изменение энергетических затрат Э и амортизационных отчислений А при подогрева нефтегазоконденсатной смеси парами тяжелой нафты в теплообменника 10Е-03 при Gо= 105508 кг/ч и Gд= 61898 кг/ч

 

В данном случае величина амортизационных отчислений на теплообменник А, которая зависит от режима его работы, интенсивно возрастает по наклонной кривой от 214,62 до 277,98 сум/кг (кривая 1 рис. 3). По нашему мнению точка пересечений кривых 1 и 2, где значение F = 445 м2, Э = 250,5 сум/кг, А = 244,8 сум/кг, Суд = 508,5 сум/кг и температура нагретой смеси на выходе из аппарата tвых =77,0 оС, характеризуют оптимальные условия эксплуатации теплообменника 10Е-03 при заданной его производительности Gо = 105508,3 кг/ч.

Заключение. Таким образом, благодаря исследованию целевой функции оптимальности подогрева нефтегазоконденсатной смеси в кожухотрубчатом аппарате выявлены оптимальные величины технологических параметров про-цесса в промышленном теплообменнике 10Е-03 Бухарского НПЗ: F = 445 м2, А = 245,8 сум/кг, Э = 250,5 сум/кг, Суд = 508,5 сум/кг и t = 77,0 оС. Предложенная методика оптимизации подогрева нефтегазоконденсатного сырья позволяет проводить анализ энергетической (тепловой) эффективности применяемых теплообменников в трех блоках тепловой подготовки углеводородного сырья к первичной перегонке, а также подобрать наиболее эффективные режимы их эксплуатации.

 

Список литературы:

  1. Глаголева О.Ф., Капустин В.М., Гюльмисарян Т.Г. и др. Технология переработки нефти. В 2-х частях. Часть I. Первичная переработка нефти/ Под ред. О.Ф. Глаголевой и В.М. Капустина. - М.: Химия, КолосС, 2006. - 400 с.
  2. Мановян А.К. Технология первичной переработки нефти и природного газа. Учебное пособие для вузов. 2-е изд. - М.: Химия, 2001. - 568 с.
  3. Бояринов А.И., Кафаров В.В. Методы оптимизации в химической технологии. Изд. второе, перераб. и доп. - М.: Химия, 1975. - 576 с.
  4. Худайбердиев А.А. Интенсификация подогрева нефтяного сырья. Монография. - Ташкент: Navroz, 2019. - 213 с.
  5. Худайбердиев А.А. Оптимизация процесса нагревания углеводородного сырья в трубчатом теплообменном аппарате// Химическая технология и экология в нефтяной и газовой промышленности: Материалы III Между-народной НПК: Булатовские Чтения: Сборник статьей/ Под общ. ред. проф. О.В. Савенок. В 5 т. - Краснодар: Юг, 31 марта 2019. - Т.4. - С.144-146.
  6. Ахметов С.А., Сериков Т.П., Кузеев И.Р., Баязитов И.М. Технология и оборудование процессов переработки нефти и газа: учебное пособие/ Под ред. С.А. Ахметова. - СПб.: Недра, 2006. - 868 с.
  7. Скобло А.И., Молоканов Ю.К., Владимиров А.И., Щелкунов В.А. Процессы и аппараты нефтегазопереработки и нефтехимии: Учебник для вузов. 3-е изд., перераб. и доп. - М.: ООО "Недра-Бизнесцентр", 2000. - 677 с.
  8. Павлов К.Ф., Романков П.Г., Носков А.А. Примеры и задачи по курсу процессов и аппаратов химической технологии. Учебное пособие для вузов/ Под ред. П.Г. Романкова. - 10-е изд., перераб. и доп. - Л.: Химия, 1987. - 576 с.
  9. Khudayberdiev A.A., Rakhimdjanova Sh.S. Technique for refined calculation of heat exchange during condensation of vapors of fractions of oil and gas condensate mixture in the shell and pipe apparatus. Scientific and technical journal of Namangan institute of engineering and technology. – Namangan, 2021. - Vol 6 Issue (4) 2021. - р. 126-133.
  10. Технологический регламент установки атмоферной перегонки смеси газоконденсата и нефти и фракционирования гидроочищенной нафты Бухарского НПЗ. - TR 16472899-001: 2009.
  11. Ахметов С.А. Технология глубокой переработки нефти и газ: учебное пособие для вузов. Изд. 2-е, перераб. и доп. - СПб.: Недра, 2013. - 544 с.
  12. Рабинович Г.Г., Рябых П.М., Хохряков П.А. и др. Расчеты основных процессов и аппаратов нефтепереработки: Справочник/ Под ред. Е.Н. Судакова. - 3-е изд., перераб. и доп. - М.: Химия, 1979. - 568 с.
  13. Григорьев Б.А., Богатов Г.Ф., Герасимов А.А. Теплофизические свойства нефти, нефтепродуктов, газовых конденсатов и их фракций/ Под общ. ред. проф. Б.А. Григорьева. - М.: Издательство МЭИ, 1999. - 372 с.
Информация об авторах

д-р техн. наук, главный научный сотрудник лаборатории «Процессы и аппараты химической технологии» Института общей и неорганической химии АН Республики Узбекистан, Узбекистан, г. Ташкент

Doctor of Technical Sciences, Chief Researcher of the Laboratory "Processes and Apparatuses of Chemical Technology" of the Institute of General and Inorganic Chemistry of the Academy of Sciences of the Republic of Uzbekistan, Uzbekistan, Tashkent

д-р техн. наук, профессор кафедры “Информатика, автоматизация и управления” Ташкентского химико-технологического института, Узбекистан, г. Ташкент

Doctor of Engineering Science, prof. “Informatics, Automation and Control” department of the Tashkent chemical-technological institute, Uzbekistan, Tashkent

базовый докторант кафедры Процессы и аппараты химической технологии, Ташкентского химико-технологического института МВиССО Республики Узбекистан, Узбекистан, г. Ташкент

Doctoral student of the Department Processes and Apparatuses of Chemical Technology of the Tashkent Chemical-Technological Institute, Republic of Uzbekistan, Uzbekistan, Tashkent

докторант Института общей и неорганической химии АН РУз, Узбекистан, г. Ташкент

Doctoral student of the Institute of General and Inorganic Chemistry of the Academy of Sciences of the Republic of Uzbekistan, Uzbekistan, Tashkent

Журнал зарегистрирован Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор), регистрационный номер ЭЛ №ФС77-54434 от 17.06.2013
Учредитель журнала - ООО «МЦНО»
Главный редактор - Ахметов Сайранбек Махсутович.
Top