Научная статья на тему 'УСТОЙЧИВОСТЬ АВТОМАТИЧЕСКОЙ СИСТЕМЫ РЕГУЛИРОВАНИЯ УРОВНЯ СТАБИЛЬНОГО КАТАЛИЗАТА В СТАБИЛИЗАЦИОННОЙ КОЛОННЕ УСТАНОВКИ КАТАЛИТИЧЕСКОГО РИФОРМИНГА'

УСТОЙЧИВОСТЬ АВТОМАТИЧЕСКОЙ СИСТЕМЫ РЕГУЛИРОВАНИЯ УРОВНЯ СТАБИЛЬНОГО КАТАЛИЗАТА В СТАБИЛИЗАЦИОННОЙ КОЛОННЕ УСТАНОВКИ КАТАЛИТИЧЕСКОГО РИФОРМИНГА Текст научной статьи по специальности «Математика»

CC BY
25
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АВТОМАТИЧЕСКАЯ СИСТЕМА РЕГУЛИРОВАНИЯ / КАТАЛИТИЧЕСКИЙ РИФОРМИНГ / МНОЖЕСТВО НЕОПРЕДЕЛЕННОСТЕЙ / НОМИНАЛЬНЫЙ ПОЛИНОМ / РОБАСТНАЯ УСТОЙЧИВОСТЬ / СЕМЕЙСТВО ПОЛИНОМОВ / СТАБИЛИЗАЦИЯ КАТАЛИЗАТА

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

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

STABILITY OF AN AUTOMATIC SYSTEM FOR CONTROLLING THE LEVEL OF STABLE CATALYZATE IN THE STABILIZATION COLUMN OF A CATALYTIC REFORMER

It is proposed to take into account the influence of uncontrolled disturbances on the level control of the catalyzate stabilization process by changing the transfer function coefficients of the automatic level control system in the stabilization column of the catalyzate stabilization unit of a catalytic reformer based on robust stability methods. To study the robust stability of an automatic system for controlling the level of the catalyzate stabilization process, a typical scheme of the control system is considered. To describe the problem in a general form, a general view of the transfer function of the automatic system for controlling the level of the catalyzate stabilization process is obtained. The problem of determining the robust stability of a closed automatic system for controlling the level of the catalyzate stabilization process is posed. A family of transfer functions of an automatic system for controlling the level of the catalyzate stabilization process with frequency uncertainty is considered. The problem of determining the robust stability of a closed-loop automatic system for controlling the level of the catalyzate stabilization process is reduced to setting the region of frequency uncertainty and constructing a hodograph of the robust Nyquist criterion. To determine the robust stability of the closed-loop automatic level control system of the catalyzate stabilization process, the nominal transfer function of the open-loop part of the automatic level control system and the function limiting the uncertainty area are given. Based on the construction of the hodograph of the robust Nyquist criterion, the closed automatic system for controlling the level of the catalyzate stabilization process has robust stability.

Текст научной работы на тему «УСТОЙЧИВОСТЬ АВТОМАТИЧЕСКОЙ СИСТЕМЫ РЕГУЛИРОВАНИЯ УРОВНЯ СТАБИЛЬНОГО КАТАЛИЗАТА В СТАБИЛИЗАЦИОННОЙ КОЛОННЕ УСТАНОВКИ КАТАЛИТИЧЕСКОГО РИФОРМИНГА»

УДК 681.5.037

Б01: 10.17277/уе81тк.2022.02.рр.235-241

УСТОЙЧИВОСТЬ АВТОМАТИЧЕСКОЙ СИСТЕМЫ РЕГУЛИРОВАНИЯ УРОВНЯ СТАБИЛЬНОГО КАТАЛИЗАТА В СТАБИЛИЗАЦИОННОЙ КОЛОННЕ УСТАНОВКИ КАТАЛИТИЧЕСКОГО РИФОРМИНГА

А. М. Джамбеков1, Б. С. Дмитриевский2, А. А. Терехова2

ГБПОУАО «Астраханский колледж вычислительной техники» (1), dzhambekovam@yandex.ru, Астрахань, Россия;

Кафедра «Информационные процессы и управление», terehova.aa@mail.tstu.ru; ФГБОУВО «ТГТУ» (2), Тамбов, Россия

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

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

Введение

В работе [1] выполнено моделирование автоматической системы регулирования (АСР) температуры низа стабилизационной колонны (СК) блока стабилизации катализата (БСК) установки каталитического риформинга (КР) [2]. Октановое число продуктовой смеси БСК (стабильного катализата) определяет качество всей цепочки процесса КР [3]. Для эффективного управления БСК помимо разработки эффективных алгоритмов управления, обеспечивающих достижение экономического или иного эффекта, необходимо обеспечение устойчивого режима работы АСР технологических параметров (температуры, давления, расхода, уровня и пр.) [4]. Одним из основных технологических параметров БСК является уро-

вень катализата в СК [5]. На регулирование уровня в СК оказывают влияние неконтролируемые возмущения: расход сбрасываемого газа стабилизации, перепад температуры в СК [6].

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

Постановка задачи

Для исследования робастной устойчивости АСР уровня в СК рассмотрим типовую схему, представленную на рис. 1, где обозначены: задающее воздействие г(0; ошибка управления е(/); управление и(/); выход у((); регулятор (ПИД-регулятор) Я; объект управления (процесс в СК) Ои [1]. В качестве регулируемого параметра рассмотрен уровень в СК. Задающее воздействие г(0 - ступенчатое изменение положения (хода) регулирующего органа (задвижки) на линии (трубопроводе) стабильного катализата [1].

Передаточная функция разомкнутой части АСР уровня в СК, определенная на основе экспериментальной переходной кривой для установки КР Л-35-11/1000, имеет вид [8]

25

Ж (*) =-2---(1)

1 + 2,285 +1,1952 + 4,1653

Для описания задачи в общем виде получим передаточную функцию (1)

к

Ж (5) =

2 3 ао + а^ + +

(2)

где к - коэффициент усиления, к = 25; а, - постоянные коэффициенты; 5 - оператор Лапласа.

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

Пусть разомкнутая часть АСР уровня в СК описывается семейством скалярных передаточных функций

Ж (5) = ВД +Д(5), (3)

где Жо(5) - номинальная передаточная функция (2); Д(5) - частотная неопределенность, удовлетворяющая условию

Д(/ю)<у| Н (уев) (4)

г(0

е(/)

Я

и(0

Ои

у(0

Рис. 1. Структурная схема АСР уровня в СК

при всех ю для некоторой функции H(s), такой, что H ^(s) - устойчивая реализуемая дробно-рациональная функция.

Необходимо знать, будет ли робастно устойчива АСР уровня в СК, замкнутая единичной обратной связью (см. рис. 1). Как известно, при отсутствии неопределенностей вопрос решается с помощью годографа Найквиста. Рассмотрим его робастную модификацию применительно к настоящей работе.

При использовании робастного критерия Найквиста предполагается, что все передаточные функции из семейства (3) имеют одинаковое число N неустойчивых полюсов с Wq (s) при всех допустимых A(s).

Построим годограф

Wo) = W^ -1, 0

H (»1

< ю < ж.

(5)

Исходя из робастного критерия Найквиста замкнутая АСР уровня в СК робастно устойчива тогда и только тогда, когда годограф Ж (_/ю) охватывает круг с центром в точке (-1, ] 0) и радиусом V N/2 раз против часовой стрелки, при этом не пересекая его.

Задача заключается в определении робастной устойчивости замкнутой АСР уровня в СК (см. рис. 1) при заданной номинальной передаточной функции разомкнутой части АСР Ж (я), частотной неопределенности Д(5) и подтверждении (не подтверждении) выполнения условия робастного критерия Найквиста, на основании которого делаем вывод о робастной устойчивости (неустойчивости) замкнутой АСР уровня в СК.

Определим робастную устойчивость замкнутой АСР уровня в СК.

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

С учетом (1) имеем номинальную передаточную функцию разомкнутой АСР уровня в СК

ад=-25---.

1 + 2,285 +1,1952 + 4,1653

Передаточная функция имеет два неустойчивых полюса: 51,2 ~ 0,059 ±] 0,769; один устойчивый 53 ~ -0,404.

Зададим функцию Н в виде

„2

H-\s) =

1 + 2s + s

1 + 3s + 4s2 + 2s3

С учетом (4) зададим частотную неопределенность

,1 + 3/ю + 4(/ю)2 + 2(/ю)3

A(jra)<v| H (jra) = v

1 + 2/ю+(/ю)2

1 - 4ю2 + jfcrn - 2ю3

1 - ю + j 2ю

(1 - 4ю2 ) +(3ю- 2ю3 )

(-«2 )2

+ 4ю

2

Все передаточные функции из семейства (3) имеют одинаковое число N = 2 неустойчивых полюсов с (5) при всех допустимых Д(5).

= v

= V

Im

40 30 20 10 0 -10 -20

/

! [ (0=оо <й=С

у

< /

-30

-10 0 10 20 30 40 50 60 Рис. 2. Годограф робастного критерия Найквиста Построим годограф робастного критерия Найквиста

_25_+1

= Wo0ю) +1 -1 = 1 + 2,28]ю + 1,19(/ю)2 + 4,16(/ю)3 -j, 0

Re

H (И

II

( - 4ю2 )+(зю-2ю3 )

< ю < да.

(-«2)

+ 4ю2

Результаты моделирования приведены на рис. 2.

Поскольку выполняется условие робастного критерия Найквиста, делаем вывод о робастной устойчивости замкнутой АСР уровня в СК.

Заключение

Таким образом, на основе построения годографа робастного критерия Найквиста замкнутая АСР уровня в СК обладает робастной устойчивостью. Данный результат может быть использован при разработке АСР расхода орошения в СК и выборе настроек регуляторов.

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

1. Джамбеков, А. М. Локальный ПИД-регулятор стабилизации катализата / А. М. Джамбеков, И. А. Щербатов // Современные технологии. Системный анализ. Моделирование. - 2016. - № 3 (51). - С. 98 - 104.

2. Ancheyta, J. Modeling and Simulation of Catalytic Reactors for Petroleum Refining / J. Ancheyta. - Hoboken : Wiley, 2011. - 528 p.

3. Smith, J. M. Chemical Engineering Kinetics / J. M. Smith. - St. Louis : McGraw-Hill, 1981. - 676 p.

4. Gumen, M. I. Increasing of the Efficiency of the Reforming LG-35-11/300 / M. I. Gumen // Petroleum Processing and Petrochemistry. - 2001. - No. 11. - P. 54 - 57.

5. Fluid Catalytic Cracking (FCC) Process Modeling, Simulation and Control / C. I. C. Pinheiro, J. L. Fernandes, L. Domingues [et al.] // Industrial Engineering Chemistry Research. - 2012. - Vol. 51 (1). - P. 1 - 29. doi: 10.1021/IE200743C

6. Weekman, V. A. Model of Catalytic Cracking Conversion in Fixed, Moving and Fluid-Bed Reactors / V. A. Weekman // Industrial and Engineering Chemistry Process Design and Development. - 1968. - Vol. 7 (1). - P. 90 - 95. doi: 10.1021/I260025A018

7. Ефремов, В. В. Расчет робастно-устойчивости автоматической системы регулирования температуры водообогреваемой спецодежды водолаза / В. В. Ефремов, А. Г. Шелудько, А. А. Анисимов // Дизайн и технологии. - 2013. - № 34 (76). -С. 75 - 79.

8. Cristea, M. V. Simulation and Model Predictive Control of a UOP Fluid Catalytic Cracking Unit / M. V. Cristea, S. P. Agachi, V. Marinoiu // Chemical Engineering and Processing: Process Intensification. - 2003. - Vol. 42, Issue 2. -P. 67 - 91. doi: 10.1016/S0255-2701(02)00017-X

9. Поляк, Б. Т. Робастная устойчивость и управление / Б. Т. Поляк, П. С. Щербаков. - М. : Наука, 2002. - 303 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

10. Nguang, S. K. Robust Stabilization of a Class of Time-Delay Nonlinear Systems / S. K. Nguang // IEEE Transactions on Automatic Control. - 2000. - Vol. 45, No. 4. - P. 756 - 762. doi: 10.1109/9.847117

11. Wu, L. Delay-Dependent Robust Stabilization for Uncertain Singular Systems with Discrete and Distributed Delays / L. Wu, W. Zhou // Journal of Control Theory and Applications. - 2008. - Vol. 6. - P. 171 - 176. doi: 10.1007/S11768-008-6166-3

12. Мартьянов, Е. И. Проблемно-ориентированная система управления и оптимизации основных параметров технически сложных систем / Е. И. Мартьянов, С. В. Карпушкин, В. В. Алексеев // Вестн. Тамб. гос. техн. ун-та. - 2021. - Т. 27, № 3. - С. 336 - 344. doi: 10.17277/vestnik.2021.03.pp.336-344

13. Математические модели многосвязных объектов управления / М. Н. М. Саиф, В. Г. Матвейкин, Б. С. Дмитриевский [и др.] // Вестн. Тамб. гос. техн. ун-та. - 2019. - Т. 25, № 1. - С. 53 - 62. doi: 10.17277/vestnik.2019.01.pp.053-062

Stability of an Automatic System for Controlling the Level of Stable Catalyzate in the Stabilization Column of a Catalytic Reformer

1 2 2 A. M. Dzhambekov , B. S. Dmitrievsky , A. A. Terekhova

Astrakhan College of Computer Engineering (1), dzhambekovam@yandex.ru, Astrakhan, Russia; Department of Information Processes and Management, terehova.aa@mail.tstu.ru; TSTU (2), Tambov, Russia

Keywords: automatic control system; catalytic reforming; a lot of uncertainty; nominal polynomial; robust stability; family of polynomials; catalyst stabilization.

Abstract: It is proposed to take into account the influence of uncontrolled disturbances on the level control of the catalyzate stabilization process by changing the transfer function coefficients of the automatic level control system in the stabilization column of the catalyzate stabilization unit of a catalytic reformer based on robust stability methods. To study the robust stability of an automatic system for controlling the level of the catalyzate stabilization process, a typical scheme of the control system is considered. To describe the problem in a general form, a general view of the transfer function of the automatic system for controlling the level of the catalyzate stabilization process is obtained. The problem of determining the robust stability of a closed automatic system for controlling the level of the catalyzate stabilization process is

posed. A family of transfer functions of an automatic system for controlling the level of the catalyzate stabilization process with frequency uncertainty is considered. The problem of determining the robust stability of a closed-loop automatic system for controlling the level of the catalyzate stabilization process is reduced to setting the region of frequency uncertainty and constructing a hodograph of the robust Nyquist criterion. To determine the robust stability of the closed-loop automatic level control system of the catalyzate stabilization process, the nominal transfer function of the open-loop part of the automatic level control system and the function limiting the uncertainty area are given. Based on the construction of the hodograph of the robust Nyquist criterion, the closed automatic system for controlling the level of the catalyzate stabilization process has robust stability.

References

1. Dzhambekov A.M., Shcherbatov I.A. [Local PID-controller for catalyzate stabilization], Sovremennyye tekhnologii. Sistemnyy analiz. Modelirovaniye [Modern technologies. System analysis. Modeling], 2016, no. 3 (51), pp. 98-104. (In Russ., abstract in Eng.)

2. Ancheyta J. Modeling and Simulation of Catalytic Reactors for Petroleum Refining, Hoboken: Wiley, 2011, 528 p.

3. Smith J.M. Chemical Engineering Kinetics, St. Louis: McGraw-Hill, 1981, 676 p.

4. Gumen M.I. Increasing of the Efficiency of the Reforming LG-35-11/300, Petroleum Processing and Petrochemistry, 2001, no. 11, pp. 54-57.

5. Pinheiro C.I.C., Fernandes J.L., Domingues L. [et al.] Fluid Catalytic Cracking (FCC) Process Modeling, Simulation and Control, Industrial Engineering Chemistry Research, 2012, vol. 51 (1), pp. 1-29, doi: 10.1021/IE200743C

6. Weekman V.A. Model of Catalytic Cracking Conversion in Fixed, Moving and Fluid-Bed Reactors, Industrial and Engineering Chemistry Process Design and Development, 1968, vol. 7 (1), pp. 90-95, doi: 10.1021/I260025A018

7. Yefremov V.V., Shelud'ko A.G., Anisimov A.A. [Calculation of robust stability of an automatic temperature control system for a diver's water-heated overalls], Dizayn i tekhnologii [Design and technologies], 2013, no. 34 (76), pp. 75-79. (In Russ.)

8. Cristea M.V., Agachi S.P., Marinoiu V. Simulation and Model Predictive Control of a UOP Fluid Catalytic Cracking Unit, Chemical Engineering and Processing: Process Intensification, 2003, vol. 42, issue 2, pp. 67-91, doi: 10.1016/S0255-2701(02)00017-X

9. Polyak B.T., Shcherbakov P.S. Robastnaya ustoychivost' i upravleniye [Robust stability and control], Moscow: Nauka, 2002, 303 p. (In Russ.)

10. Nguang S.K. Robust Stabilization of a Class of Time-Delay Nonlinear Systems, IEEE Transactions on Automatic Control, 2000, vol. 45, no. 4, pp. 756-762, doi: 10.1109/9.847117

11. Wu L., Zhou W. Delay-Dependent Robust Stabilization for Uncertain Singular Systems with Discrete and Distributed Delays, Journal of Control Theory and Applications, 2008, vol. 6, pp. 171-176, doi: 10.1007/S11768-008-6166-3

12. Mart'yanov Ye.I., Karpushkin S.V., Alekseyev V.V. [Problem-oriented control system and optimization of the main parameters of technically complex systems], Transactions of the Tambov State Technical University, 2021, vol. 27, no. 3, pp. 336-344, doi: 10.17277/vestnik.2021.03.pp.336-344 (In Russ., abstract in Eng.)

13. Saif M.N.M., Matveykin V.G., Dmitriyevskiy B.S., Bashkatova A.V., Mamontov A.A. [Mathematical models of multiply connected control objects], Transactions of the Tambov State Technical University, 2019, vol. 25, no. 1, pp. 53-62, doi: 10.17277/vestnik.2019.01.pp.053-062 (In Russ., abstract in Eng.)

Stabilität der automatischen Niveauregulierung des stabilen Katalysators in der Stabilisierungskolonne der katalytischen Reformieranlage

Zusammenfassung: Es ist vorgeschlagen, den Einfluss unkontrollierter Störungen auf die Niveauregelung des Katalysatorstabilisierungsprozesses durch Veränderung der Übertragungsfunktionskoeffizienten der automatischen Niveauregelung in der Stabilisierungskolonne der Katalysatorstabilisierungseinheit der katalytischen Reformierungseinheit basierend auf robusten Stabilitätsmethoden zu berücksichtigen. Um die robuste Stabilität eines automatischen Systems zur Steuerung des Niveaus des Katalysatorstabilisierungsprozesses zu untersuchen, ist ein typisches Schema des Steuerungssystems betrachtet. Um das Problem in allgemeiner Form zu beschreiben, ist ein allgemeiner Überblick über die Übertragungsfunktion des automatischen Systems zum Steuern des Niveaus des Katalysatorstabilisierungsprozesses erhalten. Das Problem besteht darin, die robuste Stabilität eines geschlossenen automatischen Systems zur Steuerung des Niveaus des Katalysatorstabilisierungsprozesses zu bestimmen. Es ist eine Familie von Übertragungsfunktionen eines automatischen Systems zum Steuern des Niveaus des Katalysatorstabilisierungsprozesses mit Frequenzunsicherheit betrachtet. Um die robuste Stabilität des automatischen Niveauregelsystems mit geschlossenem Regelkreis des Katalysatorstabilisierungsprozesses zu bestimmen, sind die nominelle Übertragungsfunktion des offenen Regelkreisteils des automatischen Niveauregelungssystems und die den Unsicherheitsbereich begrenzende Funktion angegeben.

Stabilité du système de contrôle automatique du niveau du catalyseur stable dans la colonne de stabilisation de l'installation de reformage catalytique

Résumé: Est proposée la prise en compte de l'effet des perturbations incontrôlées sur la régulation du niveau du processus de stabilisation du catalyseur en modifiant les coefficients de la fonction de transfert de l'installation de reformage catalytique basée sur des méthodes de la stabilité robuste. Est étudié le schéma type du système de la régulation. Est obtenu une vue d'ensemble de la fonction de transfert du système automatique de régulation du niveau du processus de la stabilisation du catalyseur. Est examinée la famille des fonctions de transfert du système avec une incertitude de fréquence. Pour la détermination de la stabilité du système automatique fermé est définie la fonction de transfert nominale de la partie ouverte du système automatique de régulation du niveau et la fonction limitant la zone d'incertitude. A la base de la construction du godographe du test robuste de Nyquist, le système automatique fermé du contrôle du niveau du processus de stabilisation de catalyseur possède une stabilité robuste.

Авторы: Джамбеков Азамат Матифулаевич - кандидат технических наук, преподаватель, ГБПОУ АО «Астраханский колледж вычислительной техники», Астрахань, Россия; Дмитриевский Борис Сергеевич - доктор технических наук, профессор кафедры «Информационные процессы и управление»; Терехова Анастасия Андреевна - аспирант кафедры «Информационные процессы и управление», ФГБОУ ВО «ТГТУ», Тамбов, Россия.

i Надоели баннеры? Вы всегда можете отключить рекламу.