б/ | = max J 6t | = max
KKn Ki<n
(1.26)
f n—2 »
(1.27)
где ta — табличное значение критерия Стьюдента при степени свободы п=2.
После отбрасывания дефектных точек вновь пересчитывают-ся коэффициенты уравнения притока и их погрешность, а также проводится проверка аномальных точек.
В случае недостаточного числа точек необходимо провести дополнительные опыты.
Пример 1.8. В табл. 1.20 приведены данные исследования скв. 2, уравнение притока которой рассчитывалось в примере 1.2. Проверить следует точку, соответствующую наименьшей депрессии.
Вычисляем: 6г=+22, +3, +1,6, ,+4,8 6~=7,25, 6~м= + 14,75. 5=10,1. При уровне надежности 0,9 с=\,АЪ. Следовательно, проверяемое значение отбрасывается, так как 14,75>1,45 • 10,1 = = 14,65.
^ После пересчета оставшихся данных получаем k = = 14,2 т-см2/кгс, Gkfk=4%. Таким образом, погрешность определения коэффициента продуктивности существенно снизилась.
Глава II
АДАПТАЦИОННЫЕ МЕТОДЫ РАСЧЕТА ПАРАМЕТРОВ ГАЗОПРОВОДА
Нестационарное изотермическое течение газа по горизонтальному газопроводу описывается системой нелинейных дифференциальных уравнений в частных производных вида [47]:
c
dt дх
(ИЛ)
dp
где R — универсальная газовая постоянная; Т — абсолютная! температура газа; g — ускорение свободного падения; р — давление в газопроводе; х—пространственная координата; t— временная координата; X — коэффициент гидравлического сопротивления; р — плотность газа; D — диаметр трубы; Q — расход газа; с — скорость звука в газе; z — коэффициент сжимаемости газа.
Задавшись соответствующими граничными и начальными условиями, можно оценить коэффициенты, входящие в уравнения.
При рассмотрении методов определения коэффициентов, входящих в дифференциальное уравнение в частных производных, предпочтение отдается методам рекуррентного оценивания, позволяющим разрабатывать алгоритмы, наилучшим образом приспособленные для реализации их на ЭВМ. Это прежде всего-метод стохастической аппроксимации — метод, наиболее близкий к существующим на газопроводе реальным условиям, когда измерения проводятся с аддитивными, центрированными и-независимыми случайными погрешностями.
В общем виде рассматривается объект, описываемый уравнением [15]
. ,v a»*(a, Q1 , x(s, t\ —^—J, |
—------- F[s U с x(s t\ —^
где х(s, t)—вектор, характеризующий состояние объекта в. момент t в точке 5; с — вектор неизвестных параметров.
30
Состояние объекта измеряется в дискретные моменты времени п—\, 2 ... и в конечном ряде точек г=1, 2 ..., расположенных на расстоянии Аг. Заменяя уравнение в частных производных уравнением в разностной форме, получим
x[r, n]=f(c, x[r—l, n — \], х[г, п— 1], х[г+ 1, Л Вводится показатель качества
Ф(Х[г, л], *[г, л])}
где М — оператор математического ожидания; х[г, п] — результат измерения.
Алгоритм оценки параметра с запишется в виде
F(c[n~ 1], дг[г— 1, п— 1],
Х[Г+1, Я— 1]),
где V— символ градиента; у(п) = \/п — переменный коэффициент.
Результаты применения метода стохастической аппроксимации для оценки неизвестных параметров позволяют взять этот метод за основу при разработке алгоритмов оценки фактического коэффициента гидравлического сопротивления линейной части магистрального газопровода [33].
Пример II.1. Алгоритм оперативной оценки коэффициента гидравлического сопротивления элементарного участка однони-точного газопровода.
Данный алгоритм позволяет получить оперативную оценку коэффициента гидравлического сопротивления элементарного участка в процессе эксплуатации магистрального газопровода. Предполагается, что общий расход газа измеряется в некоторой точке магистрального газопровода и в точках отбора газа потребителями.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.