![]() Главная страница Случайная страница КАТЕГОРИИ: АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника |
Использование принципа суперпозиции в расчетах внедрения краевой воды в газовую залежь круговой формы.
Имеется залежь радиусом Rз, заданы Qдобст(t), Требуется рассчитать qв(t), Qв(t), ¶2Р/¶r2 +(¶Р/¶r)/r =1/c× ¶Р/¶t - линейное уравнение. Воспользуемся принципом суперпозиции для решения линейного уравнения Qв(tn)= Предположим, что на укрупненной скважине наложены n одинаковых скважин, работающие с своим дебитом Dqвj, но разные промежутки времени. Т.о. мы можем использовать принцип наложения скважин, у каждой которой постоянный дебит, но разные по значению. Т.е. в пространстве мы имеем несколько скважин, но совмещенных вместе Pн-Р(Rз, t)=SDPj; j=1, n (2) Pн-Р(Rз, t)=Pн-(mв/(2p× k× h))× S[Dqвj× fon-foj-1=χ (tn-tj-1)/R2З tj-1=S Dtj=(j-1) Dt (4) Принимаем при t0=0 => Dqв0=0 (5) Выделим из (3) последние слагаемое: Pн-Р(Rз, tn)=Pн-(mв/(2pkh))× S[Dqвj× qвj=qвj-1+Dqвj Теперь предположим, что все показатели процесса поступления воды ч/з стенку укрупненной скважины нам известны на момент времени tn-1. Требуется определить показатели на момент времени tn-? Qв(tn)=Qв(tn-1)+[qв(n-1)+Dqвn]× Dt (7) Запишем уравнение материального баланса для ВНР с учетом (7)
Ур. (6) описывает потери давления в области м/у контуром водонапорной системы и контуром начальной газоносности.
Р(Rз, t)-[ ln[Rз/R(t)]× [qвn-1+Dqвn] (10) Контур укруп.скв. не изменен во времени, контур газоносности изменяется. Приравниваем (6) и (10) уравнения и подставляем ур.(8), получим Исключая Р из (12) с учетом (9) и (5) получим: Рн-(mв/(2pkh))(× S[Dqвj×
+rв× g× y(t)+(mв/(2pkh))× ln(Rз/R(tn))× [qвn-1+Dqвn] (11) (11) квадратное отн-но Dqвn Dqв(tn)=b/(2× a)-[(b2/(4× a2)-c/a]0.5 (12) а=mвDt× (2× p× kВ× h))ln[Rз/R(tn)]) b=РнDt-(Dt× qв(tn-1)mв/(2× p× kB× h))× ln[Rз/R(tn)]+L× mв× /(2× p× k× h)- -(Dtmв/(2× p× k× h))× S[Dqвj× +(L× mв/(2× p× kB× h))× ln[Rз/R(tn)]- rв× g× y(t) Dt c=Рн× L-(L× mв/(2× p× k× h))× S[Dqвj× L= d=(Рн В (12) входят параметры на момент времени tn: R(tn), y(tn), z(tn). Поэтому решение производят методом последовательных итераций. В 1-м приближении: R(1)(tn)=R(tn-1); y(1)(tn)=y(tn-1); z(1)(tn)=z(tn-1) Далее вычисляют коэф-ты а, b, с, но предварительно d и L и находят Чтобы уточнить R(tn), используют соотношение QB(tn)=π [R2з-R2(tn)]mh( Из (13)определяют R(1)(tn). Чтобы уточнить у(1)(tn) надо с помощью структурной карты и карты емкостного параметра (α mh) устанавливают зависимость у=f( ymax=H-этаж газонос-и Итерации ведутся до сходимости Р. Рез-ты расчетов сравнивают по давлениям: ½ Р(2)(tn)-Р(1)(tn)½ £ e
|