Методы повышения эффективности процессов добычи и транспорта газа, страница 43


]tNk{t)dt = bN.

.0

;,   Для случая фильтрации

j    При дифференцировании bo = O,   bx = 1,   Ь2 = Ь3 =

Для того чтобы   получить   прогнозированное   значение на интервал, надо положить Ьо=1, Ьх = {—Д*),. Ь2=(—At)2, ..., bN

 () Дисперсия ошибки

G2Q=2nGxx~-,                                                                                       (Ш.34)

где Gxx — спектральная плотность случайных ошибок.

В связи с тем, что при решении промысловых задач прихо­дится сталкиваться с обработкой экспериментального или стати­стического материала, снятого в дискретные моменты или в дискретных точках, приходится интегральный оператор заменить дискретным при t=nAt:

 (111.35)

t=/z

где т — число интервалов памяти; п — номер точки, соответст­вующей моменту, для которой вычисляется сумма; i — текущий номер точки; Q[iAt] — значение реализации в iточке; А^ — шаг дискретности; Q'[nAt] — дифференцированное значение случай­ной функции.

K(t) определяется из выражения K(t) = ii,0-hni(t), где цо = — 6/Т2, 1^! =—123. Весовые функции представляют собой поли­ном 1-го порядка


п

S

t=nm


 m2


QlnAt],


где а и b — коэффициенты весовой функции. jxo=— , j^i==------------------------------------------------------------------------- .

Коэффициенты весовой функции для полиномов до 4-го по­рядка приведены в работах [27, 40].

6   Зак. 2194                                                                                                                             81


С целью проверки возможности применения этого метода в данном случае проводилось его экспериментальное исследование на аналоге-модели.

Предварительная расчетная схема процесса строилась с уче­том следующих предпосылок.

1.  Существует непроницаемая перемычка и переток между
пластами отсутствует.

2.  В начальный период разработки забойное давление равно»
пластовому давлению пласта с меньшим давлением.

Эти предположения с достаточной для практики точностью характеризуют процесс.

ZQ, см3/с 2,0 2,5 3,0

20      W     60     Ь,мин Йис. Ш.5. Зависимость 2Q от /

3,5

10 20 дО W 50 60 70 4 мин Рис. III.6. Зависимость dQ/dt от t.

В данном случае приток жидкости или газа в скважину про­исходит из более мощного пласта. В процессе разработки место­рождений пластовое давление работающего пласта падаетг следовательно, уменьшается и дебит скважины. Для поддержа­нии заданных отборов необходимо снижать забойное давление. Снижение забойного давления вызовет в определенный момент-подключение пласта с меньшим давлением. Дебит этого пласта будет возрастать по мере проводимых мероприятий по снижению» забойного давления.

В лабораторных условиях аналогами пластов служили ко­лонки с пористой средой, наполненные газом. Давление газа в колонках контролировалось манометрами. Дебит в общем кол* лекторе регистрировался с помощью газовых часов. Первона­чально к коллектору подключалась колонка с большим давле* нием. Замеры снимались через каждые 30 с. В момент снижения давления в первой колонке до величины давления во второй-к коллектору подсоединялась вторая колонка. Регистрировался момент ее включения и продолжалось снятие замеров по газо­вым часам, которые в данном случае регистрировали суммар­ный дебит двух колонок.

По результатам эксперимента был построен график зависи­мости суммарного дебита от времени (рис. Ш.5). Статистиче­ское дифференцирование полученной кривой осуществлялось на-ЭВМ «Наири-2».

На основании полученных данных была построена зависи­мость, показанная на рис. III.6. Эта зависимость характеризует

82


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

ОПРЕДЕЛЕНИЕ СОРБЦИОННОЙ  СПОСОБНОСТИ ПОРОДЫ