Файл: Герман, В. Т. Построение информационной системы управления технологическими процессами добычи и подготовки газа научно-экономический обзор.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 03.11.2024
Просмотров: 39
Скачиваний: 0
можно определить дебит скважин с |
вероятностью |
J3 = 0,9 5 , откло |
|
нением |
6 = 0,15 и точностью |
6уз= 0,05 . |
|
N = 36.
Это значит, что для оценки с точностью 5% среднего дебита скважин достаточно из всей совокупности выделить для обследования 36 сква жин.
ОПРЕДЕЛЕНИЕ ПОСЛЕДОВАТЕЛЬНОСТИ СБОРА ТЕХНОЛОГИЧЕСКОЙ ИНФОРМАЦИИ ПРИ ЦЕНТРАЛИЗОВАННОМ КОНТРОЛЕ
И УПРАВЛЕНИИ
Реальная оценка состояния процессов добычи и подготовки газа
возможна только по результатам анализа большого набора различной промысловой технологической информации. Поток информации при цен
трализованном контроле и управлении представляет собой один из видов дискретного потока информации, который называется опросным.
При этом |
необходимо |
определить |
периодичность |
сбора информации о |
каждом из |
параметров процессов |
добычи и подготовки газа, а также |
||
последовательность |
опроса каждого параметра в |
отдельности. |
Рассмотрим построение программ сбора информации или рассредо
точенную во времени |
последовательность опроса параметров техноло |
|
гических процессов [ |
5 ]. Интервал времени, через который повторяе |
|
тся вся программа, назовем циклом программы - Т. Цикличность |
вы |
|
бирается, исходя из |
числа контролируемых на пульте параметров |
и |
их измерений. Так, если интервал времени между замерами парамет - ров составляет несколько минут (часов), то временем цикла можно
выбрать час (сутки). Необходимо построить такую последовательность опроса, при которой на протяжении цикла число опросов в каждый ин тервал времени было примерно одинаковым, т .е . чтобы информация со
биралась равномерно. |
Эту |
задачу можно формализовать |
следующим |
|
образом. |
|
|
|
|
Обозначим наличие опроса через "I", а отсутствие |
его - через |
|||
"0П. Тогда, разбив цикл Т |
на интервалы времени (для часового цик |
|||
ла - минуты, суточного - |
часы), можно еоотавить по каждому пара - |
|||
метру набор из единиц и нулей. |
|
|
||
Имея такие наборы по каждому параметру |
и сложив их по минутам |
|||
(часам), получим определенное число опросов |
по минутам (часам) . |
|||
Обозначим эти суммы |
Sj |
( j = I , 2 , . . . , 6 0 ) . |
Циклический сдвиг впра |
во или влево на произвольное число минут (часов) будем называть
преобразованием набора. При этом число единиц и интервал между ни ми в наборе для каждого параметра остаются, меняются лишь моменты
начала опроса. Задачу можно сформулировать следующим образом: пре образовать все наборы таким образом, чтобы получить равномерное распределение сумм опросов по минутам (часам). Для решения постав
ленной задачи можно применить следующий алгоритм. |
|
|
|||||||||||
|
Представим |
набор, |
соответствующий |
|
-му параметру |
в виде |
|||||||
вектора |
Xj |
с |
компонентами |
X - |
( I |
= I , |
. . . , N ) часового ци |
||||||
кла |
j |
= 1 , 2 , . . . , |
60, |
тогда Xi=X£(Xf, ,Х£е, - / . « ) • |
Величины Ху |
||||||||
принимают значения 0 или I . |
Подразумевая под |
преобразованием век |
|||||||||||
тора циклический сдвиг его компонент, обозначим Х[ |
1 (X;, |
) |
|||||||||||
вектор, полученный из вектора Х;(Х;, ,Х^г ; ..., |
Xt-fi0) |
путем цикли |
|||||||||||
ческого сдвига компонент, проделанного |
К; |
раз, |
К 6 1 6 |
60. |
|||||||||
|
Расстояние |
между векторами |
X; (Х;,(...)и Ху (X;,,...^Определяется |
||||||||||
по известной |
формуле______________________________________ |
|
|||||||||||
|
|
j3 (Xj ?Xj)=^Xji- X;;^ + (Xj2-X^z) + ■• • + (Xjeo ~ Хсео)г , |
|
||||||||||
|
Через p (Xjt Xj |
J) |
обозначают расстояние |
между векторами X; |
|||||||||
и |
Xj |
после |
каждого |
Xj |
-го |
преобразования вектора Xj . |
Далее |
||||||
фиксируют вектор |
X, |
, |
а вектор |
Хг |
подвергают |
преобразованию. |
После каждого преобразования необходимо определить расстояние меж ду фиксированным и преобразованным векторами. Получив 59 значений
у э (х ,, x |
/ Kf)), |
выбирают максимальное (если |
их несколько, то |
любое |
||||||||
из них). Определение вектора |
Х2 К;1 |
по принципу максимума рассто |
||||||||||
яния между векторами |
X, и |
Х2 |
дает |
рассредоточение |
отличных от |
|||||||
нуля компонент вектора |
Хг |
по отношению |
X, . |
После |
этого |
нахо |
||||||
дят |
сумму векторов |
X, |
и х 2*г5 |
и аналогично поступают |
с векторами |
|||||||
Х1 + Х2 К^ |
, подвергая преобразованию вектор Х3 . |
Под суммой пони |
||||||||||
мают обычную алгебраическую сумму векторов. |
|
|
|
|||||||||
|
По^ле^ний |
шаг - |
нахождение максимума расстояния между векто - |
|||||||||
рами |
£ |
X; |
и XN |
|
и определение вектора X(nKn |
. Б |
результате |
|||||
получают равномерную загрузку каждого интервала времени |
внутри |
|||||||||||
цикла. На этом заканчивается решение поставленной задачи. |
|
|||||||||||
|
Приведенный алгоритм реализуется на ЭЦВМ и дает хорошие ре - |
|||||||||||
зультаты. Блок-схема приведена на рис.7. |
|
|
|
|
||||||||
|
Пример. Рассмотрим группу скважин, |
подключенных |
к |
пульту |
||||||||
телеконтроля. |
Для каждой скважины известен интервал времени |
между |
22
Рис.7. |
Блок-схема алгоритма оптима |
/ |
Ввод исходных данных |
|||||
льного |
рассредоточения |
сбора |
техно |
|||||
логической информации |
|
|
|
\ |
|
|||
|
|
|
|
|
|
|
Циклический сдвиг компо |
|
опросами (табл.9 ). В таблице |
нали - |
2нент исходных данных |
||||||
3 |
t |
|
||||||
чие опроса обозначено "I", |
отсутст |
|
||||||
вие "О". Сумма опросов |
по |
часам |
Sj |
Определение расстояния |
||||
|
между векторами р |
|||||||
|
|
|
|
|
|
к. |
* |
|
|
|
|
|
|
|
Запоминание величин |
||
/>[(xf + x 2+x3) tx * ]= V 0 ; |
|
|
|
|
расстояний |
р |
||
|
|
|
|
f |
|
|||
f[(K1+x2t x 3) , x ^3 i w |
; |
|
|
|
Li. |
|
||
|
|
|
ного р |
|
||||
|
|
|
|
|
|
|
Определение максималь |
|
|
|
|
|
|
|
Li. |
+ |
|
4 ( V |
x 2+x3)4x ?]=V 02 . |
|
|
|
Сложение преобразован |
|||
|
|
|
ных векторов |
|||||
|
|
|
|
|
|
L |
♦ |
|
|
|
|
|
|
|
Печать полученной после |
||
|
|
|
|
|
|
довательности |
опроса |
|
|
|
|
|
|
|
|
I |
|
Дальнейшие преобразования |
век |
|
( конец) |
|
||||
|
|
|
тора X будут приводить к ранее рас смотренным сочетаниям. Сопоставляя полученные величины расстояний, мо
жно заключить, что максимальное рассредоточение векторов получает
ся при К = *■ (табл.1 0 ). |
|
Т а б л и ц а |
9 |
Последовательность опроса газовых скважин |
|
№________________ Опросы через время, час______________ ____
ишы" I 2 3 И 5 б 7 8 9 10 |
II 12 13 14 |
15 |
16 |
17 |
18 |
1920 21 |
22 25 24 |
|||
1 |
I 0 0 I 0 0 I 0 0 I . 0 |
О I О О |
I |
О О |
I |
О О I 3 О |
||||
2 |
I I I I I I I I I I I I I I I I I I I I I I I I |
|||||||||
3 |
0 I 0 I 0 I 0 I 0 1 |
О I О I О |
I |
О I |
О |
I 0 1 О I |
||||
4 0 0 0 1 0 0 0 1 0 0 |
О I О О О |
I О О |
О |
I О О О I |
||||||
Sj, |
2 2 1 ^ 1 2 2 3 1 3 |
I |
3 2 |
2 I |
4 |
I |
2 |
2 |
3 1 |
3 I 3 |
23
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
|
10 |
|
||||
|
Преобразованная последовательность |
опросов |
газовых |
|
|
|
||||||||||
|
|
|
|
|
скважин |
|
|
|
|
|
|
|
|
|
|
|
№ |
|
|
Опросы через |
время, |
час |
|
|
|
|
|
|
|
2h |
|||
сква |
г 2 |
3 L 5 |
б ? 8 9 10 II |
I2 |I3 |W |
15 |
1б|17 |
18 |
19 |
20 |
21 |
22 |
2 3 |
||||
жины |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
I 0 0 I 0 0 I 0 0 I 0 |
0 I 0 0 I 0 0 |
I 0 0 |
I 0 0 |
||||||||||||
2 |
I I |
1 I I |
I X I I I I |
I |
I |
I |
I |
I |
1 |
I |
I |
I |
I |
I |
I |
1 |
|
|
т |
т |
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
0 I 0 I 0 I 0 I 0 I 0 |
I 0 I 0 I 0 I 0 Т 0 1 |
0 I |
|||||||||||||
k |
0 0 I 0 0 0 I 0 0 0 I |
0 0 0 Т |
0 0 0 I 0 0 0 I |
0 |
||||||||||||
|
2 2 2 3 I 2 3 2 I 3 2 |
2 2 2 2 3 I 2 |
3 |
2 |
I |
3 |
2 |
2 |
||||||||
|
|
|
АДАПТАЦИЯ ПРОЦЕССА СБОРА |
|
|
|
|
|
|
|||||||
|
|
ТЕХНОЛОГИЧЕСКОЙ ИНФОРМАЦИИ |
|
|
|
|
|
|||||||||
|
Увеличение объема данных, характеризующих состояние и режимы |
|||||||||||||||
технологического оборудования, требует разработки алгоритмов, ко |
||||||||||||||||
торые повысили бы эффективность процесса управления сбором инфор |
||||||||||||||||
мации, когда система подвержена воздействию множества дестабили |
- |
|||||||||||||||
зирующих факторов, причем степень |
и время |
воздействия |
или |
|
сами |
|||||||||||
факторы |
заранее неизвестны. |
Одним из методов повышения эффектив |
||||||||||||||
ности является |
применение принципа адаптации. |
|
|
|
|
|
|
|
||||||||
|
Процесс адаптации системы проявляется в изменении состояния, |
|||||||||||||||
структуры, параметров, алгоритма работы и т .д . |
[3,4-]. |
|
|
|
|
|||||||||||
|
Алгоритмы адаптивной дискретизации |
позволяют |
производить вы |
борку данных в моменты поступления существенной для диспетчера ин формации об исследуемом управляемом процессе, когда изменяются па раметры потоков информации. Для выбора адаптивного алгоритма уп равления сбором информации необходимо знать статистику источников информации. Измерительные сигналы процессов добычи и подготовки газа в силу наличия инерционных элементов можно рассматривать как совокупность непрерывных функций времени с ограниченной вариацией,
имеющих на интервале определенные ограниченные первые производные.
Для такого класса функций оптимальны функции - f i t ) , описывающие сигнал в виде ряда
(1)
L- U
2k
Устройство адаптивной дискретизации должно контролировать со ответствие сигнала f ( f ) его модели вида (I ):
SuP l f ( t ) - P ( i ) j = Sup £ ( t ) - £ gon |
( 2 ) |
|||
[0,T ] |
t e [ 0 , T ] |
, |
||
|
||||
где О,Т - интервал адаптивной дискретизации; |
В(i),C^0R - те- |
|||
кущая и допустимая ошибки приближения. |
|
|
||
В момент времени Т, |
когда условие |
(2) нарушается, выбирают |
информацию об измеряемом сигнале. Для создания устройства необхо
димо выбрать систему базисных функций |
( f ) , найти способ |
вы |
|||||
числения коэффициентов |
С; |
и алгоритм вычисления и |
контроля |
||||
ошибки приближения |
согласно |
выражению (2 ) . Этот алгоритм |
будет |
||||
одновременно определять интервал дискретизации, т .е . управлять |
|||||||
сбором |
информации. |
Поскольку |
на сигнал |
f ( t ) наложена |
случайная |
||
помеха, |
причем с |
более |
высокочастотным |
спектром, то |
желательно, |
чтобы алгоритм обладал фильтрующими свойствами и не включал one -
рации дифференцирования |
сигнала. С |
этой |
точки зрения |
в качест |
ве базисных функций ряда |
(I ) могут |
быть |
выбраны системы |
ортогона |
льных полиномов. Для упрощения технической реализации |
алгоритмов |
|||
адаптивной дискретизации, повышения точности работы и расширения |
||||
класса функций, допускающих обработку с |
помощью такого |
алгоритма, |
в качестве базиса целесообразно выбрать ортогональную систему по
линомов Хаара. |
Она обладает следующим свойством: любая непрерывная |
|||||||||
функция f |
(х ) |
на интервале [ о ;l] |
разлагается |
в равномерно |
сходя |
|||||
щийся ряд |
вида |
(I) |
по функциям системы; |
замена переменных |
~ |
|||||
преобразует интервал |
O iX ^ I |
в интервале |
. |
|
||||||
Система Хаара состоит |
из |
кусочно-постоянных функций и |
для |
|||||||
всех натуральных |
|
m > 0 |
и |
|
„ m - 1 |
|
|
|||
|
|
2 |
определена следующим |
|||||||
образом: |
|
|
m~l |
|
|
|
|
|
|
|
|
|
|
г |
При |
X & В mj |
|
|
|
||
Xmj( x H |
- 2 |
г |
при |
X& 6 +mj |
|
|
( 3) |
|||
|
|
|
|
при |
|
е mj |
|
|
|
|
где |
|
|
|
j |
|
, причем |
Х00 W |
= 1- |
|
|
|
m- 1 |
) |
2 ш-,- ] |
|
|
|||||
Преимущество аппроксимации по |
системе (3) |
заключается в прос |
тоте технической реализации вычисления коэффициентов Cl ряда ( I ) .
Недостатком системы (3) является наличие внутренних точек разрыва функций, образующие систему. Для устранения его предлагается ие-
25