Научная статья на тему 'Автоматизация технологических процессов добычи гидроминерального сырья'

Автоматизация технологических процессов добычи гидроминерального сырья Текст научной статьи по специальности «Энергетика и рациональное природопользование»

CC BY
278
29
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГИДРОЛИТОСФЕРНЫЕ ПРОЦЕССЫ / ДОБЫВАЮЩИЕ СКВАЖИНЫ / ОПТИМИЗАЦИЯ / ДЕПРЕССИОННАЯ ВОРОНКА / HYDROLITHOSPHERE PROCESSES / PRODUCTION WELLS / OPTIMIZATION / DEPRESSION FUNNEL

Аннотация научной статьи по энергетике и рациональному природопользованию, автор научной работы — Першин Иван Митрофанович, Малков Анатолий Валентинович, Веселов Геннадий Евгеньевич, Першин Максим Иванович

Проблема рационального природопользования в области добычи гидроминерального сырья требует применения новых технологических решений, обеспечивающих экологическую безопасность эксплуатации гидроминеральной базы при возрастающих потребностях в гидроминеральном сырье. Используя принципы и методы систем с распределенными параметрами, удается (при заданных полях депрессионной воронки отражающих желаемый технологический режим эксплуатации месторождения минеральных вод) обеспечить увеличение дебита. Желаемый технологический режим обеспечивается системой управления параметрами гидролитосферного процесса. Входным воздействием на объект управления (гидролитосферный процесс) служат дебиты добывающих скважин. Функцией выхода понижение уровня в зонах расположения скважины. Проблема выбора исполнительного устройства (определение числа добывающих скважин и их расположения) является одной из центральных при решении задачи проектирования распределенных систем управления гидролитосферными процессами. Добывающие скважины оказывают существенное воздействие на понижение уровня в зонах их расположения. Математически описать это воздействие можно с использованием соотношения Джейкоба, использующего значения гидравлических параметров водоносных горизонтов. Определение этих параметров сложная и дорогостоящая задача. В статье рассматривается технология проектирования системы управления одной добывающей скважиной и переход к распределенным системам управления гидролитосферными процессами. Приводится структура и методика вычисление параметров аппроксимирующего звена описывающего изменение уровня по пространственной координате. Для определения параметров рассматриваемого звена проводятся либо экспериментальные исследования, либо численное моделирование объекта управления. Используя аппроксимирующее звено, в статье рассматривается решение задачи выбора числа добывающих скважин исходя из обеспечения максимального дохода за выбранный срок эксплуатации месторождения. При этом было задано одинаковым понижение уровня в зонах расположения добывающих скважин, обеспечивающее заданный технологический режим эксплуатации месторождения, а так же геометрические границы расположения скважин. Безусловно, выбор числа добывающих скважин может осуществляться исходя и из других критериев. Проектирование распределенных систем управления совокупностью добывающих скважин рассмотрено в литературе по синтезу систем с распределенными параметрами.

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

Похожие темы научных работ по энергетике и рациональному природопользованию , автор научной работы — Першин Иван Митрофанович, Малков Анатолий Валентинович, Веселов Геннадий Евгеньевич, Першин Максим Иванович

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

AUTOMATION OF HYDROMINERAL RAW MATERIALS PRODUCTION TECHNOLOGICAL PROCESSES

The problem of rational nature management in the field of mining of hydromineral raw materials requires the use of new technological solutions that ensure the environmental safety of operating the hydromineral base with increasing demands for hydromineral raw materials. Using the principles and methods of systems with distributed parameters, it is possible (with the given fields of the depression funnel reflecting the desired technological mode of operation of the mineral water field) to increase the flow rate. The desired process mode is provided by the system for controlling the parameters of the hydro-lithospheric process. The input impact on the control object (hydro-lithospheric process) is the flow rates of producing wells. The exit function is to lower the level in the zones where the well is located. The problem of choosing an actuator (determining the number of production wells and their location) is one of the central ones in solving the design problem of distributed control systems for hydro-lithospheric processes. Production wells have a significant impact on lowering the level in the areas of their location. Mathematically, this effect can be described using the Jacob ratio, using the values of the hydraulic parameters of the aquifers. Determining these parameters is a complex and expensive task. The article discusses the technology of designing a single production well control system and the transition to distributed control systems for hydro-lithospheric processes. The structure and methodology for calculating the parameters of the approximating link describing the change in the level of the spatial coordinate is given. To determine the parameters of the considered link, either experimental studies or numerical modeling of the control object are carried out. Using an approximation link, the article discusses the solution to the problem of selecting the number of producing wells based on ensuring maximum income for a selected field life. At the same time, the same level of reduction was set in the zones of production wells, ensuring the specified technological mode of field operation, as well as the geometric boundaries of the wells. Of course, the choice of the number of producing wells can be made on the basis of other criteria. Designing distributed control systems for a set of production wells is discussed in the literature on the synthesis of systems with distributed parameters.

Текст научной работы на тему «Автоматизация технологических процессов добычи гидроминерального сырья»

Ильюшин Юрий Валерьевич - Санкт-Петербургский горный университет; e-mail: Ilyushin_YuV@pers.spmi.ru; 199106, Санкт-Петербург, Васильевский остров, 21 линия, д. 2; тел.: 89111783523; кафедра системного анализа и управления; доцент.

Першин Иван Митрофанович - e-mail: Pershin_IM@pers.spmi.ru; тел.: 89187903619; кафедра системного анализа и управления; профессор.

Ilyushin Yury Valerievich - St. Petersburg Mining University; e-mail: Ilyushin_YuV@pers.spmi.ru; 199106, St. Petersburg, Vasilievsky Island, 21 line 2; phone: +79111783523; the department of system analysis and control; associate professor.

Pershin Ivan Mitrofanovich - e-mail: Pershin_IM@pers.spmi.ru; phone: +79187903619; the department of system analysis and control; professor.

УДК 28;50 DOI 10.23683/2311-3103-2018-5-198-210

И.М. Першин, А.В. Малков, Г.Е. Веселов, М.И. Першин

АВТОМАТИЗАЦИЯ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ ДОБЫЧИ ГИДРОМИНЕРАЛЬНОГО СЫРЬЯ

Проблема рационального природопользования в области добычи гидроминерального сырья требует применения новых технологических решений, обеспечивающих экологическую безопасность эксплуатации гидроминеральной базы при возрастающих потребностях в гидроминеральном сырье. Используя принципы и методы систем с распределенными параметрами, удается (при заданных полях депрессионной воронки отражающих желаемый технологический режим эксплуатации месторождения минеральных вод) обеспечить увеличение дебита. Желаемый технологический режим обеспечивается системой управления параметрами гидролитосферного процесса. Входным воздействием на объект управления (гидролитосферный процесс) служат дебиты добывающих скважин. Функцией выхода - понижение уровня в зонах расположения скважины. Проблема выбора исполнительного устройства (определение числа добывающих скважин и их расположения) является одной из центральных при решении задачи проектирования распределенных систем управления гидролитосферными процессами. Добывающие скважины оказывают существенное воздействие на понижение уровня в зонах их расположения. Математически описать это воздействие можно с использованием соотношения Джейкоба, использующего значения гидравлических параметров водоносных горизонтов. Определение этих параметров - сложная и дорогостоящая задача. В статье рассматривается технология проектирования системы управления одной добывающей скважиной и переход к распределенным системам управления гидролитосферными процессами. Приводится структура и методика вычисление параметров аппроксимирующего звена описывающего изменение уровня по пространственной координате. Для определения параметров рассматриваемого звена проводятся либо экспериментальные исследования , либо численное моделирование объекта управления. Используя аппроксимирующее звено, в статье рассматривается решение задачи выбора числа добывающих скважин исходя из обеспечения максимального дохода за выбранный срок эксплуатации месторождения. При этом было задано одинаковым понижение уровня в зонах расположения добывающих скважин, обеспечивающее заданный технологический режим эксплуатации месторождения, а так же геометрические границы расположения скважин. Безусловно, выбор числа добывающих скважин может осуществляться исходя и из других критериев. Проектирование распределенных систем управления совокупностью добывающих скважин рассмотрено в литературе по синтезу систем с распределенными параметрами.

Гидролитосферные процессы; добывающие скважины; оптимизация; депрессионная воронка.

I.M. Pershin, A.V. Malkov, &Е. Veselov, M.I. Pershin

AUTOMATION OF HYDROMINERAL RAW MATERIALS PRODUCTION TECHNOLOGICAL PROCESSES

The problem of rational nature management in the field of mining of hydromineral raw materials requires the use of new technological solutions that ensure the environmental safety of operating the hydromineral base with increasing demands for hydromineral raw materials. Using the principles and methods of systems with distributed parameters, it is possible (with the given fields of the depression funnel reflecting the desired technological mode of operation of the mineral water field) to increase the flow rate. The desired process mode is provided by the system for controlling the parameters of the hydro-lithospheric process. The input impact on the control object (hydro-lithospheric process) is the flow rates ofproducing wells. The exit function is to lower the level in the zones where the well is located. The problem of choosing an actuator (determining the number ofproduction wells and their location) is one of the central ones in solving the design problem of distributed control systems for hydro-lithospheric processes. Production wells have a significant impact on lowering the level in the areas of their location. Mathematically, this effect can be described using the Jacob ratio, using the values of the hydraulic parameters of the aquifers. Determining these parameters is a complex and expensive task. The article discusses the technology of designing a single production well control system and the transition to distributed control systems for hydro-lithospheric processes. The structure and methodology for calculating the parameters of the approximating link describing the change in the level of the spatial coordinate is given. To determine the parameters of the considered link, either experimental studies or numerical modeling of the control object are carried out. Using an approximation link, the article discusses the solution to the problem of selecting the number ofproducing wells based on ensuring maximum income for a selected field life. At the same time, the same level of reduction was set in the zones of production wells, ensuring the specified technological mode of field operation, as well as the geometric boundaries of the wells. Of course, the choice of the number of producing wells can be made on the basis of other criteria. Designing distributed control systems for a set ofpro-duction wells is discussed in the literature on the synthesis of systems with distributed parameters.

Hydrolithosphere processes; production wells; optimization; depression funnel.

Состояние проблемы. Гидролитосфера - это одна из основных оболочек биосферы, от сохранения целостности которой во многом зависят все живые существа. Неблагоприятная ситуация с поверхностными водами, утилизация отходов, интенсивная добыча минеральных ресурсов - основные причины антропогенного воздействия на гидролитосферу, следствием этого воздействия является истощение и загрязнение водных ресурсов [1-4]. Управление технологическими процессов добычи гидроминерального сырья предполагает решение следующих задач: совершенствование методов определения гидравлических параметров водоносных горизонтов; совершенствование методов построения математических моделей рассматриваемых процессов (принципы и методы описания математических моделей рассматриваемых объектов приведены в [5-10]); совершенствование принципов и методов управления рассматриваемыми процессами. Входным воздействием на объект управления служат дебиты добывающих скважин. Функцией выхода - понижение уровня в зонах расположения скважины. Поскольку определение гидравлических параметров водоносных горизонтов - сложная и дорогостоящая задача, то рассматривается задача аппроксимации динамических и статических характеристик объекта с использованием экспериментальных данных.

Методика выбора числа добывающих скважин, использующая известное решение Джейкоба, рассмотрена в [1, 11] . Эта методика ориентирована на вычисление достаточно большого число гидродинамических параметров объекта с помощью обработки большого объема опытно-фильтрационных работ. Рассматриваемая в статье методика использует данные реакции объекта управления на статическое и динамическое входные воздействия, что существенно упрощает решение поставленной задачи.

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

У

ь „х

Скважина

h jfx.y.z.t)

H2(x,y,z,T)

х

-X

"Zn

Рис. 1. Схема месторождения Математическая модель рассматриваемого объекта записывается в виде [7-9].

Грунтовые воды

д/1( х, у, г,т) д2к1( х, у, г,т) д2\( х, у, г,т) д2\( х, у, г,т)

1,х 1,у 1,г 2

дт

дН2 (х, y, z, т) дт

дх2 1,у ду

0 < х < Lx;0 < у < Ly ;0 < г < Ь^. Эксплуатационный пласт

д2 H2( х, y, г,т) + k д2 H2( х, y, z,t)

+ k

дг'

д2 H2 (х, y, z, т)

дгУ

дх2 2'у ду2 2

+¥ (т)-д( х0, у0, г0); 0 < х < Ь ;0 < у < Ь ;0 < г < Ь . Граничные условия (условия Дарси) между пластами задаются в виде Грунтовые воды - верхний пласт

/ (х, ^ ^ т) = ^ ( х y, Ьч,т) + Ь1 ■ (Н2 (x, У,0, т) - ^ х У, Ьг1,т)) ,

Н2 (х,у,0, т) = Н2 (х,у, 0, т) -Ь ■ (Н2 (х,у, 0,т)-/ (х,у,Ь^, т)).

где Ь1 = 0,00003 сут.-1 - параметр перетекания. Нижняя граница пласта:

дН2 (х, у, Ь , т)/ дг = 0 Боковые грани: / (0, у, г, т) = / 0; Н 2 (0, у, г, т) = Н 2 01,

д/\(Ьх, y, Z, т) / дх = 0; Н 2 (Ьх, У, г т) = Н 2,02.

/(х,0,г,т) = /(х,Ьу,г,т) = / 0,

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

Н2 (х,0, г, т) = Н2 (х, Ly, г, т) = Н2 01 - (Н2 01 - Н2 02) • х /Lx,

где ^,0=г (0<z<Lz1), Н201=85 м., Н202=65 м. - начальные состояния невозмущенных грунтовых вод и пласта.

Физические параметры объекта [1]: ^ - напор в горизонте грунтовых вод; Н2 - напор в изучаемом водоносном горизонте; коэффициенты фильтрации по соответствующим координатам: к1х= 0.198 м/сут., ^,у= 0.196 м/сут.,k1,!:= 0.148 м/сут., 0.146 м /сут.ЛЯу= 0,148 м /сут.Л2,г= 0.024 м/сут.; п2= 0.000581/м. - уп-ругоемкость пласта; V(т) - понижение напора, вызванное воздействием (дебитом) добывающей скважины; 5(х0, у0, г0) - функция, равная единице, если х=х0, у=у0, z=z0, и равная нулю в других случаях; х,у,г - пространственные координаты; т - время.

Геометрические данные месторождения приведены в табл. 1.

Таблица 1

Геометрические данные

Lx Ly Lz1 Lz2

350 м. 300 м. 45 м. 150 м.

При моделировании объекта управления будем полагать, что число точек дискретизации по координатам: х^ (N^13), у^ (Ыу=13),г^-( N^=9, N^=9). Соответственно шаги дискретизации были заданы в виде: Лх^х/(^-1); Лу=Ьу/(^-1); Лг1=Ьг1/(^1-1); Az2=Lz2/(Nz2-1); Добывающая скважина осуществляет забор гидроминерального сырья в «точке» Nx0=7, N^=7, Nz20=5. Входным воздействием на объект управления служит дебит добывающей скважины Q(т), который связана с функцией V(т) следующим соотношением

V(т)=Кv•Q(т). (1)

Поскольку имеется рабочая скважина, с использованием которой осуществляют добычу гидроминерального сырья. При этом дебит рассматриваемой скважины составляет 1000м3/сут. Проведем эксперимент на реальном объекте - скачком увеличим дебит на 100м3/сут. (1000+100) .

(100м3/сут.=100/(3600^24)=0.0011574м3/сек.).

Понижение уровня в установившемся режиме, в области расположения заборного устройства скважины, составило 3.5м. Коэффициент передачи объекта управления^, равный отношению функции выхода (понижения уровня) к входному воздействию (0.001157)

К:=3.5/0.0011574=3024.019.

Верификация математической модели объекта управления. Процедура верификации распадается на следующие этапы:

1. Выберем в формуле (1) значение К^0.0308. Используя численную модель объекта управления, промоделируем гидролитосферный процесс при входном воздействии Q(т)=-I00 м3/сут. (знак - в математической модели показывает, что мы извлекаем из пласта гидроминеральное сырье). Подавая входное воздействие Q Ю(т) = -100 м3/сут.=-100/(3600-24)=-0.001157 м3/сек.), определим изменение уровня в области расположения скважины (х=х0, у=у0, г=г0). Это изменение составило -2.408456 м. (график рис. 2).

15 30 45 60 Время в сутках

Рис. 2. График изменения уровня

2. Вычислим поправку к выбранному коэффициенту К . Эта поправка будет равна Кр = 3.5 м./2.408456 м. = 1.45213 .

Приведенная выше формула (1) преобразуется к виду:

У(т)=КуКр ^(т) =0.0308-1.45213^(т). (2)

Используя численную модель объекта управления, с учетом (2 ), промоделируем гидролитосферный процесс при входном воздействии Р(т)=-100 м3/сут. По результатам моделирования построен график, приведенный на рис. 3.

Рис. 3. График изменения уровня в скорректированной модели пласта в точке

х=хо у=уа г=го

Используя результаты моделирования, определим статический коэффициент передачи математической модели объекта управления К1= (изменение уровня в установившемся режиме)/(входное воздействие).

К1=-3.499097/(-0.001157)=3023.239.

Статические коэффициенты передачи реального объекта и математической модели рассматриваемого объекта близки.

На рис. 4 приведены значения рассогласования уровней в точках дискретизации (разница между полем уровней в установившемся режиме при воздействии Р=100 м3/сут. и полем уровней невозмущенного состояния пласта).

Дх

Дх

z =75м.

х

ду-i л у-:

■ (.0 - 0.0 - 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 од

0.000000 -0.011467 -0.022900 -0.034013 -0.043989 -0.051330 -0.054149 -0.051330 0.043939 -0.034013 -0.022900 -0.011467 0.000000

0.000000 -0.013859 -0.043151 -0.072734 -0.095104 -0.114593 -0.122332 -0.114598 -0.096104 -0.072734 -0.048161 -0.023869 0.000000

0.000000 -0.0377« -0.077663 -0.121000 -0.166Ш -02 03344 ■ 0.230288 -0.208844 -0.166880 -0.121000 -0.077668 -0.037792 0.000000

0.000000 -0.052753 -0.111403 -0.182301 -0 J 70948 -0J74037 -0.450920 -0.374037 -02 70943 -0.182301 -0.111403 -0.052753 0.000000

0.000000 0.000000 -0.066064 -0.072 1 31 -0.14422 9 -0.16146 5 -0.2=2122 -0J 99357 -0.420=88 -0.57070 7 -0.697409 -1.261722 ■1.068271 1-3.4990971 ■0.697409 -1.261722 0.42 0533 0.57(17(17 -0.2=2122 -0J 99357 0.14422 9 -0.16146 5 .0.0661)64 -0.072181 0.000000 0.000000

о.оооою -0.066022 -0.144150 -0J 52017 -0.420465 -0.697277 -1.068135 -0.597277 -0.420155 -0J 52017 -0.144150 -0.056Ш2 0.000000

0.000000 -0.052 675 0.111264 -0.182109 -0.270724 -0.373797 -0.450674 -0J73797 -ОЛ70724 -0.182109 -0.111264 -0.052 675 0.000000

0.000000 -0.037693 -0.077433 -0.120750 -0.155503 -02 03545 -0329984 -0 J 08545 -0.166603 -0.120750 -0.077488 -0.037693 0.000000

0.000000 -0.023775 -0.047990 -0.072 507 -0.095842 -0.114317 -0.12204= -0.114317 -0.09=342 -0.072 507 -0.047990 -0.023775 0.Ю0000

0.000000 0.0

-0.011407 0.0

-0.022 7Я 0.0

-0.033871 0.0

-0.043826 0.0

-0.051155 0.0

-0.053970 0.0

-0.051155 0.0

-0.043326 0.0

-0.033871 0.0

-0.022792 0.0

-0.011407 0.0

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

0.000000 0.0

Рис. 4. Рассогласования уровней

Определение параметров аппроксимирующего звена. Определим реакцию объекта на динамическое входное воздействие

Q(т)=0.00П574•sm(Wl •х),ю1=0.00001.

По результатам моделирования построен график, приведенный на рис. 5

(3)

Рис. 5. Реакция объекта на динамическое входное воздействие в точке

Х=Х0, у=у0, 2=20

С использованием динамических характеристик (см. рис. 5), определим сдвиг по фазе сигнала выхода относительно входного сигнала (период колебаний входного сигнала составляет 174.53 час.) Дф1 = -2^14/174.53 = -0.54 рад.

Методика определения параметров аппроксимирующего звена. Исследования, приведенные в [12, 13] показывают, что для описания динамических характеристик рассматриваемых объектов может быть рекомендовано следующее аппроксимирующее звено:

Жа (^ =7К ' еХР М) , "

р +1 ^ а

где ^-параметр, учитывающий радиус влияния скважины, 8-оператор Лапласа. В рассматриваемом случае О=(то^у)2, К и а параметры, значения которых определяются в соответствии со следующей методикой:

1. Полагая в (1) ю=0 и приравнивая статический коэффициент усиления аппроксимирующего звена значению Кь получим:

3023.239 = — • ехр(— в), Р = ((*/300)2)12 . в + 1 ^ н»

Решая уравнение, получим К=3087.05719458.

2. Полагая в (1) 8=]ю (ю= , запишем соотношение для определения фазы аппроксимирующего звена

Аф =— 1т(вш1) — arctan(Im(вШ1)/(ве(ви1) +1)), Рю1 = ^ + О ^. (4)

3. Подставляя в исходное уравнение (4) «¡=0.00001, Дф1= - 0.54 и решая численным методом, получим: а=0.0000525.

Аппроксимирующее звено, для рассматриваемого выше объекта, записывается в виде

3087.05719458 / ч ( з , . V2 (5)

(з) =--ехр(—Р), Р= -+ (^/300)2 . (5)

а Р +1 ^ ' Р \ 0.0000525 )

Определение влияния дебита скважины на различные участки пласта.

Рассмотренное выше звено описывает взаимосвязь понижения уровня (в окрестности г0=1) и дебита скважины (рис. 6).

¡Л

Скважина ^^

Колодец /

г

-ь-

г0

W (s) exp {-ß-(Го + k • r)), ß = \s + G I . (6)

ß +1 V a )

W (r) = -^1 • exP {-ß • (го + k •r)), ß = {GГ •

Рис. 6. Схема для расчета понижения уровня в произвольной точке

Введем в аппроксимирующее звено (3) составляющую Ьг, учитывающую изменение уровня от г в статическом режиме (см. рис. 6).

К ■■ ехр(—Р-(гп + k• г)), р = {3

Запишем соотношение (6) в статическом состоянии К

Р+1

где k - параметр, значение которого определяется с использованием экспериментальных данных, или результатов моделирования объекта управления.

Для рассмотренного выше процесса, полагая г0=1, статическая составляющая может быть записана в виде

3087.05719458 ^ Л п /^/300)П>2 (7)

(г) =-——---ехр (—Р-(1 + k • г)), Р = ((я7 300) ) . (7)

Для определения значения параметра к воспользуемся вычисленными значениями рассогласования уровней ЛН2(х=х0,уу^0) , (у=1., ^у) (см. рис. 4 и рис. 7).

0.000000 -0.072181 -0.151465

-0.570707 -1.251723

|ЗТ&М>97|

-1.261722 -0.570707 ■ -0:25ВД=7 -0,161465 -0.072 1 81 0.000000

Рис. 7. Рассогласования уровней (^Н2(х=Хо,Уу,2о))

Вычисление параметра к будем осуществлять исходя из минимума суммы квадратов отклонений функции (7) в точках дискретизации (yY = 1...Nyo), от вычисленных значений, приведенных на рис. 7.

• exp(— ß • (r • k + ro )• Q-AH2(Xo, yy, Zo)2 ^ min .

7=1 ß + 1

Вычисленное значение параметра k = 3.724976.

Взаимосвязь параметров депрессионной воронки в статическом режиме, по координате у (r+r0 ),(см. рис. 6, рис. 1) для рассматриваемого месторождения записываются в виде

AH 2(r + ro) = 308705719458 •G • exp (—р. (r • 3.724976 + -)).

ß +1

(8)

Рис. 8. График рассогласования уровней

Полученная модель депрессионной воронки (8) достаточно точно описывает процесс ЛН2(х=х0,уу^0) (см. рис. 8 и рис. 7). При проектировании распределенной системы управления одной из центральных задач является задача выбора оптимального число добывающих скважин.

Выбор оптимального числа добывающих скважин. При выходе на стационарный режим понижение уровня в точке установки /-ой скважины определяется суммой влияний /-й добывающей скважины, всех остальных добывающих скважин рассматриваемого пласта, имеющих дебит Q1. При этом, параметры К и г0 , для каждой /-й скважины могут быть разные, соответственно К1 и г01. Рассматриваемое влияние, в соответствии с (8) описывается следующим соотношением:

K • Q t v n K, • Q,

AHi = -ß-^ • exp (— ß • (Год > + £ jf-ß + 1 j=1. j *• ß + 1

• exp

(— ß • (-j + j ß = (G)12, (9)

(1=1...п),

где АН1 - понижение уровня в рассматриваемой 1-й скважине (1=1...п); г , -расстояние от 1-й до _|-й взаимодействующей скважины (рис. 9); п - число скважин; Q1 - дебит 1 - ой скважины; QJ - дебит ]-ой скважины.

Рис. 9. Схема определения riv

Полагая

С = — • exp (- ß • (Год ))' CUj = ^

в +1 ' Ь] в +1 соотношение (9) преобразуется к виду

• exP (

(- ß • (1, + r0,i)), ß = (G)12,

ah = ci • q + x с,, • qj, c.

(10)

j=1, j * i

Запишем уравнение (10) в матричном виде

AHj "

ah2 =

_AH_

с с -с с,, с, с,_

& &

Преобразуя, получим матричное уравнение для определения дебитов добывающих скважин

q1

q2 =

q

с с ••• с с с ••• с

2,1 ' 2 ' 2,n

с с •••с

n,1 ' n,2> n

-1 ahj"

ah2

_ah_

(11)

Матричное уравнение (11) позволяет определить дебиты добывающих скважин, исходя из заданных значений АН1 - понижений уровня (1=1.. .п). Суммарный дебит может быть определен из соотношения

Q = Qj

j=1

(12)

Положим, что имеется предприятие, осуществляющее добычу гидроминерального сырья, которое продает потребителям. Для рассматриваемого месторождения минеральных вод: заданное понижение уровня составляет 4м.; г01=1м.; цена 1 м3 добытого сырья, за вычетом НДС, составляет 250 руб.; налог на недропользование 7.5 % ; затраты на бурение, обустройство и содержание одной скважины за 10 лет эксплуатации составляют 25 млн. руб. Прибыль за десять лет эксплуатации рассматриваемого месторождения составит:

П=д-250-(365-10)-0.925 - 25000000-п. (13)

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

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

Методика решения задачи оптимизации. Рассматриваемая методика распадается на следующие этапы:

1. Для заданного значения Ly и выбранного числа скважин п (полагая, что крайние скважины (1 и п) расположены на расстоянии 2-Ду от границ (см. рис. 9)), определим Ду= Ьу/( п+3).

2. По соотношениям (11)-(12) вычисляем дебиты добывающих скважин и суммарный дебит.

3. Используя соотношение (13) определим прибыль за десять лет эксплуатации месторождения. По результатам расчетов построены графики, которые приведены на рис. 10.

1 3 5 7 & 11 и 15 17 ЧИСЛО скважин

Рис. 10. Результаты расчетов

Как следует из расчетов (см. рис. 10) при заданном значении Ly=300 м, и заданных условиях (AHj = 4м., г0>1=1м., цена 1 м3 добытого сырья, за вычетом НДС, составляет 250 руб; налог на недропользование 7.5 %; затраты на бурение, обустройство и содержание одной скважины за 10 лет эксплуатации составляют 25 млн. руб.) оптимально выбрать 7 добывающих скважин. При этом суточный дебит составит 457.17 м3, а прибыль за 10 лет эксплуатации месторождения 210.87 млн. руб. (см. рис.10). Вычисленное значение шага между добывающими скважинами Ду=30м. (см. рис. 9).

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

Методика синтеза распределенных систем управления достаточно подробно изложена в [13-20]. Эта методика предполагает, что определено число добывающих скважин и их расположение. Если имеется возможность экспериментального исследования, то для оценки динамических характеристик, выбираются две пространственные моды входного воздействия и исследуется поведение объекта [1, 16]. Если не имеется возможности экспериментального исследования, то используют численную модель объекта управления. В результате реализации методики синтезируется распределенный регулятор, который управляет дебитами всех добывающих скважин.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Малков А.В., Першин И.М., Помеляйко И.С., Першин М.И., Уткин В.А., и др. Кисловод-ское месторождение углекислых минеральных вод. Системный анализ, диагностика, прогноз и управление. - М.: Наука, 2015. - 283 с.

2. Малков А.В., Помеляйко И.С., Першин И.М. Особенности эколого-гидрогеологического мониторинга курортных территорий на примере Кисловодска // Вопросы курортологии, физиотерапии и лечебной физической культуры. - 2018. - Т. 95. - Вып. 2. 2'. - С. 80-82.

3. Першин И.М.,Малков А.В., Першин М.И. Оперативное и стратегическое управление режимами эксплуатации гидролитосферных объектов // Недропользование XXI век. - Январь 2014. - № 6а (44). - С. 40-47.

4. Помеляйко И.С., Першин И.М. Системный анализ экологического состояния зоны ги-пергенеза курорта Кисловодска // Вестник Северо-Кавказского федерального университета. - 2013. - № 3 (36). - С. 74-80.

5. Першин И.М, Малков А.В., Першин М.И. Математические модели гидролитосферных процессов // Недропользование XXI век. - 2018. - № 3. - С. 56-65.

6. Martirosyan A.V., Martirosyan K.V., Pershin I.M. Analysis of the Caucasus Mineral Waters' Field's Modeling // Modern Applied Science. - 2015. - Vol. 9. - P. 204.

7. Pershin I.M., Drovosekova T.I. Peculiarities of modelling hydro-lithospheric processes in the region of KavkazskiyeMineralnyeVody (caucasus mineral springs) // Proceedings of the 19th International Conference on Soft Computing and Measurements, SCM 2016. 7519732.

- P 215-217.

8. Pershin I.M., Pervukhin D.A., Ilyushin Y. V., Afanaseva O.V. Design of distributed systems of hydrolithosphere processes management // A synthesis of distributed management systems. IOP Conf. Series: Earth and Environmental Science 87. - 2017. 032029. - P. 1-7.

9. Drovosekova T.I., Pershin I.M. Peculiarities of modelling hydro-lithospheric processes in the region of Kavkazskiye Mineralnye Vody (caucasus mineral springs) // Proceedings of the 19th International Conference on Soft Computing and Measurements, SCM 2016. 7519732.

- P. 215-217.

10. Малков А.В., Першин И.М., Першин М.И. Оптимизация количества водозаборных скважин при эксплуатации месторождений минеральных вод // Недропользование XXI век.

- Февраль 2015. - № 1 (51). - С. 100-104.

11. Першин И.М., Веселое Г.Е., Першин М.И. Методы аппроксимации передаточных функций распределенных объектов // Системный синтез и прикладная синергетика: Сб. научных трудов VII Всероссийской научной конференции. - 2015. - С. 106-117.

12. Першин М.И. Исследование погрешностей динамических характеристик распределенных объектов при аппроксимации // Современная наука и инновации. - 2014. - Вып. № 4 (8). - С. 46-50.

13. Малков А.В., Першин И.М. Системы с распределенными параметрами. Анализ и синтез.

- М.: Научный мир, 2012. - 476 с.

14. Малков А.В., Першин И.М. Синтез распределенных регуляторов для систем управления гидролитосферными процессами. - М.: Научный мир, 2007. - 364 с.

15. Першин И.М., Дубогрей В.Ф, Малков А.В. Методика синтеза распределенных систем управления режимами эксплуатации месторождений минеральных вод // Известия вузов. Геология и разведка. - 2012. - № 2. - С. 74-78.

16. Веселов Г.Е., Першин М.И. Проектирование распределенных систем управления гидро-литосферными процессами // Известия вузов. Геология и разведка. - 2016. - № 1.

- С. 99-105.

17. Першин И.М., Кузьмин Н.Н., Малков А.В. Формирование целевых функций в задачах управления гидролитосферными процессами // 5-я Российская мультиконференция по проблемам управления: Материалы конференции «Информационные технологии в управлении» (ИТУ-2012). - 2012. - С. 622-632.

18. Малков А.В., Першин И.М. Синтез распределенных регуляторов для систем управления гид-ролитосферными процессами // Мир лингвистики и коммуникации. - 2007. - № 2. - С. 235.

19. Першин И.М., Веселов Г.Е., Першин М.И. Синтез распределенных систем управления гидролитосферными процессами месторождений минеральных вод // Известия ЮФУ. Технические науки. - 2014. - № 8 (157). - С. 123-137.

20. Pershin I.M., Pervukhin D.A., Ilyushin Y.V., Afanaseva O.V. Design of distributed systems of hydrolithospere processes management // Selection of optimal number of extracting wells IOP Conf. Series: Earth and Environmental Science 87. - 2017. 032030. - P. 1-6.

REFERENCES

1. Malkov A.V., Pershin I.M., Pomelyayko I.S., Pershin M.I., Utkin V.A., i dr. Kislovodskoe mestorozhdenie uglekislykh mineral'nykh vod. Sistemnyy analiz, diagnostika, prognoz i upravlenie [Kislovodsk deposits of acidulous waters. System analysis, diagnostics, forecasting and management]. Moscow: Nauka, 2015, 283 p.

2. Malkov A.V., Pomelyayko I.S., Pershin I.M. Osobennosti ekologo-gidrogeologicheskogo monitoringa kurortnykh territoriy na primere Kislovodska [Features of ecological and hydro-geological monitoring of resort areas on the example of Kislovodsk], Voprosy kurortologii, fizioterapii i lechebnoy fizicheskoy kul'tury [Questions of balneology, physiotherapy and physical therapy], 2018, Vol. 95, Issie 2. 2', pp. 80-82.

3. Pershin I.M., Malkov A.V., Pershin M.I. Operativnoe i strategicheskoe upravlenie rezhimami ekspluatatsii gidrolitosfernykh ob"ektov [Operational and strategic management modes of operation hydrolithosphere objects], Nedropol'zovanie XXI vek [Subsoil use XXI century], January 2014, No. 6a (44), pp. 40-47.

4. Pomelyayko I.S., Pershin I.M. Sistemnyy analiz ekologicheskogo sostoyaniya zony gipergeneza kurorta Kislovodska [System analysis of the ecological state of the hypergenesis zone of Kislovodsk], Vestnik Severo-Kavkazskogo federal'nogo universiteta [Bulletin of the North Caucasus Federal University], 2013, No. 3 (36), pp 74-80.

5. Pershin I.M, Malkov A. V., Pershin M.I. Matematicheskie modeli gidrolitosfernykh protsessov [Mathematical models hydrolithosphere processes], Nedropol'zovanie XXI vek [Subsoil use XXI century], 2018, No. 3, pp. 56-65.

6. Martirosyan A.V., Martirosyan K.V., Pershin I.M. Analysis of the Caucasus Mineral Waters' Field's Modeling, Modern Applied Science, 2015, Vol. 9, pp. 204.

7. Pershin I.M., Drovosekova T.I. Peculiarities of modelling hydro-lithospheric processes in the region of KavkazskiyeMineralnyeVody (caucasus mineral springs), Proceedings of the 19th International Conference on Soft Computing and Measurements, SCM 2016. 7519732, pp. 215-217.

8. Pershin I.M., Pervukhin D.A., Ilyushin Y. V., Afanaseva O.V. Design of distributed systems of hydrolithosphere processes management, A synthesis of distributed management systems. IOP Conf. Series: Earth and Environmental Science 87, 2017. 032029, pp. 1-7.

9. Drovosekova T.I., Pershin I.M. Peculiarities of modelling hydro-lithospheric processes in the region of Kavkazskiye Mineralnye Vody (caucasus mineral springs), Proceedings of the 19th International Conference on Soft Computing and Measurements, SCM 2016. 7519732, pp. 215-217.

10. Malkov A.V., Pershin I.M., Pershin M.I. Optimizatsiya kolichestva vodozabornykh skvazhin pri ekspluatatsii mestorozhdeniy mineral'nykh vod [Optimization of the number of water wells in the operation of mineral water fields], Nedropol'zovanieXXI vek [Subsoil use XXI century], February 2015, No. 1 (51), pp. 100-104.

11. Pershin I.M., Veselov G.E., Pershin M.I. Metody approksimatsii peredatochnykh funktsiy raspredelennykh ob"ektov [Methods of approximation of transfer functions of distributed objects], Sistemnyy sintez i prikladnaya sinergetika: Sb. nauchnykh trudov VII Vserossiyskoy nauchnoy konferentsii [System synthesis and applied synergetics: Proceedings of the VII all-Russian scientific conference], 2015, pp. 106-117.

12. Pershin M.I. Issledovanie pogreshnostey dinamicheskikh kharakteristik raspredelennykh ob"ektov pri approksimatsii [Investigation of errors of dynamic characteristics of distributed objects in the approximation], Sovremennaya nauka i innovatsii [Modern science and innovation], 2014, Issue No. 4 (8), pp. 46-50.

13. Malkov A.V., Pershin I.M. Sistemy s raspredelennymi parametrami. Analiz i sintez [Systems with distributed parameters. Analysis and synthesis]. Moscow: Nauchnyy mir, 2012, 476 p.

14. Malkov A.V., Pershin I.M. Sintez raspredelennykh regulyatorov dlya sistem upravleniya gidrolitosfernymi protsessami [The synthesis of distributed controllers for control systems hydrolithosphere processes]. Moscow: Nauchnyy mir, 2007, 364 p.

15. Pershin I.M., Dubogrey V.F, Malkov A.V. Metodika sinteza raspredelennykh sistem upravleniya rezhimami ekspluatatsii mestorozhdeniy mineral'nykh vod [Methods of synthesis of distributed control systems modes of operation of mineral water deposits], Izvestiya vuzov. Geologiya i razvedka [Proceedings of Higher Schools. Geology and Exploration], 2012, No. 2, pp. 74-78.

16. Veselov G.E., Pershin M.I. Proektirovanie raspredelennykh sistem upravleniya gidrolitosfernymi protsessami [Design of distributed control systems hydrolithosphere processes], Izvestiya vuzov. Geologiya i razvedka [Proceedings of Higher Schools. Geology and Exploration], 2016, No. 1, pp. 99-105.

17. Pershin I.M., Kuz'min N.N., Malkov A.V. Formirovanie tselevykh funktsiy v zadachakh upravleniya gidrolitosfernymi protsessami [The formation of the objective functions in problems of control hydrolithosphere processes], 5-ya Rossiyskaya mul'tikonferentsiya po problemam upravleniya: Materialy konferentsii «Informatsionnye tekhnologii v upravlenii» (ITU-2012) [5th Russian multi-conference on management: Proceedings of the conference "Information technologies in management" (ITU-2012)], 2012, pp. 622-632.

18. Malkov A.V., Pershin I.M. Sintez raspredelennykh regulyatorov dlya sistem upravleniya gidrolitosfernymi protsessami [The synthesis of distributed controllers for control systems hydrolithosphere processes], Mir lingvistiki i kommunikatsii [The world of linguistics and communication], 2007, No. 2, pp. 235.

19. Pershin I.M., Veselov G.E., Pershin M.I. Sintez raspredelennykh sistem upravleniya gidrolitosfernymi protsessami mestorozhdeniy mineral'nykh vod [Synthesis of distributed control systems hydrolithosphere processes of mineral waters], Izvestiya YUFU. Tekhnicheskie nauki [Izvestiya SFedU. Engineering Sciences], 2014, No. 8 (157), pp. 123-137.

20. Pershin I.M., Pervukhin D.A., Ilyushin Y.V., Afanaseva O.V. Design of distributed systems of hydrolithospere processes management, Selection of optimal number of extracting wells IOP Conf. Series: Earth and Environmental Science 87, 2017. 032030, pp. 1-6.

Статью рекомендовал к опубликованию д.т.н. А.Б. Чернышев.

Першин Иван Митрофанович - Северо-Кавказский федеральный университет, филиал в г. Пятигорске; e-mail: ivmp@yandex.ru; 357500, г. Пятигорск, просп. 40 лет октября, 56; тел.: 88793973927; кафедра управления в технический системах; зав. кафедрой; д.т.н.; профессор.

Малков Анатолий Валентинович - ООО «Нарзан-гидроресурсы»; e-mail: Anatol.Malkov@yandex.ru; 357700, г. Кисловодск, ул. Кирова, 43; тел.: +78793729139; директор; к.г-м.н.; д.т.н.

Веселов Геннадий Евгеньевич - Южный федеральный университет; e-mail: gev@sfedu.ru; 347900, г. Таганрог, ул. Чехова, 2; тел.: +78634360450; Институт компьютерных технологий и информационной безопасности; директор; д.т.н.; профессор.

Першин Максим Иванович - e-mail: maksimpershin@bkmail.ru; тел.: 89280093030; Институт компьютерных технологий и информационной безопасности; аспирант.

Pershin Ivan Mitrofanovich - North-Caucasian Federal University, a branch in the town of Pyatigorsk; e-mail: ivmp@yandex.ru; 357500, Pyatigorsk, ave. 40 years on October 56; phone: +78793973927; the department of management of technical systems; head of department; dr. of eng. sc.; professor.

Malkov Anatoly Valentinovitch - LLC "Narzan-hydro"; e-mail: Anatol.Malkov@yandex.ru; Kislovodsk, Kirov Street, 43; phone: +78793729139; Director; cand. of geological-mineralogical sc.; dr. of eng. sc.; professor.

Veselov Gennady Evgen'evich - Southern Federal University; e-mail: gev@sfedu.ru; 2, Chekhov street, Taganrog, 347900, Russia; phone: +78634360450; Institute of Computer Technology and Information Security; director; dr. of eng. sc.; professor.

Pershin Maksim Ivanovich - e-mail: maksimpershin@bkmail.ru; phone: +79280093030; Institute of Computer Technology and Information Security; postgraduate student.

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