Фильтрация несмешивающихся жидкостей в призабойной зоне скважины тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

ВВЕДЕНИЕ.

ГЛАВА I. КОРРЕКТНОСТЬ ПОСТАНОВКИ И КАЧЕСТВЕННЫЕ СВОЙСТВА ОБОБЩЕННЫХ РЕШЕНИЙ НЕКОТОРЫХ ВЫРОЖДАЮЩИХСЯ

ЗАДАЧ ДВУХФАЗНОЙ ФИЛЬТРАЦИИ

§ I. Единственность и устойчивость обобщенных решений основных начально-краевых задач фильтрации . несмешиващихся жидкостей.

1. Исследование линейной сопряженной задачи

2. Теорема единственности.

3. Теорема о сходимости решений регуляризованных задач.

4. Теорема о непрерывной зависимости от коэффициентов и начально-краевых условий.

§ 2. О корректности постановки одномерной и осе-симметричной задач вытеснения при заданном перепаде давления.

§ 3. Качественные свойства обобщенных решений одномерных и осесимметричных задач вытеснения при заданных перепаде давления или суммарном потоке

1. Леммы о монотонной зависимости обобщенного решения от функции суммарного расхода

2. Лемма о влиянии функций капиллярного давления и суммарного расхода на характер поведения обобщенного решения.

3. Теорема о мажорантной и минорантной оценках обобщенного решения одномерной задачи вытеснения

§ 4. Асимптотическое поведение при бесконечном возрастании времени обобщенного решения одномерной задачи вытеснения с заданным суммарным расходом

ГЛАВА 2. ЧИСЛЕННОЕ И ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ОДНОМЕРНОЙ И ОСЕСИММЕТРИЧНОЙ ЗАДАЧ ВЫТЕСНЕНИЯ.

§ I. Численный расчет нестационарных задач вытеснения при заданных перепаде давлений или суммарном по

1. Описание алгоритма расчета одномерной задачи

2. Стационарные решения.

3. Результаты численных расчетов одномерной задачи

4. Осесимметричный случаи.

§ 2. Решение задачи вытеснения методом последовательной смены стационарных состояний (МПССС).

§ 3. Экспериментальное исследование вытеснения смачивающей жидкости несмачиващей.

1. Методика проведения эксперимента

2. Результаты эксперимента и?сравнение с расчетами

Глава 3. ВОЗМОЖНЫЕ СПОСОБЫ УСТРАНЕНИЯ ЭФФЕКТА КАПИЛЛЯРНОГО

ЗАПИРАНИЯ ВЫТЕСНЯЕМОЙ ФАШ.

§ I. Математическая модель воздействия виброисточника на призабойную зону скважины.

§ 2. Численное моделирование эффекта вибровоздействия конечно-разностным методом

§ 3. Стационарное решение задачи вытеснения в неоднородной пористой среде

§ 4. Упругий режим стационарной фильтрации с капиллярным запиранием смачивающей фазы

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

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

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

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

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

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

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

Целью данной работы является всестороннее исследование фильтрации двухфазной жидкости в призабойной зоне скважины в следующих аспектах:

- Построение и исследование корректности математической модели, описывающей этот процесс.

- Способы численно-аналитического решения получаемой системы дифференциальных уравнений.

- Экспериментальное моделирование вышеописанного явления вытеснения фильтрата бурового раствора цри вызове притока к скважине.

- Теоретическое обоснование возможных методов воздействия на призабойную зону для повышения дебита скважины.

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

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

Необходимо отметить, что к настоящему времени существует ряд работ, в которых приведены результаты изучения разнообразных задач, возникающих в теории фильтрации несмешивающихся жидкостей и ее приложениях. Можно упомянуть работы Л.С.Лейбензона /34/, Г.И.Баренблатта, В.М.Ентова, В.М.Рыжика /8/, И.А.Чарного /51/, В.Л.Данилова, Р.М.Каца /14/, Р.Коллинза /24/, А.Н.Коновалова /26/, А.Э.Шейдеггера /57/ и ряда других авторов.

Общепринятой моделью для исследования фильтрации несмешивающихся жидкостей является модель Маскета-Леверетта /3/, или, в одномерном случае, модель Рапопорта-Лиса /77/.

Уравнения, описывающие фильтрацию двухфазной жидкости в не-деформируемой изотропной однородной пористой среде выводятся из уравнений неразрывности для каждой фазы /24/

Jl + cLlv tf » О (I) d -t

-mlf + oLvT/ = 0 (2) дЬ и обобщенного закона Дарси (без учета сил гравитации) = -lHvp <3)

--Kvp, (4) ь —>

Здесь р^, р , , - давления и скорости фильтрации для смачивающей и несмачивагацей фаз соответственно, S - насыщенность несмачивающей фазой, \УЪ - постоянная пористость, , у - относительные фазовые проницаемости,

К » К^ - отнесенные к вязкостям абсолютные проницаемости фаз (т.е. проницаемости при полном насыщении пор данной жидкостью, К = Jt*j - вязкости, тЬ - абсолютная проницаемость). Функции , ^ обладают следующими свойствами /34/, с.104: l(o)=to=o, <»

Заметим, что, без ограничения общности, пористость можно считать равной единице (если не отмечено особо).

Система уравнений (I) - (2) замыкается форь/улой Лапласа для капиллярного давления, которую согласно Леверетту /71/ можно представить в виде р-р* = рк(5) (6) где рк - функция капиллярного давления, обладающая свойствами: dLpK/oljS >0 ; рк(0) =0

При качественном исследовании системы (I) - (4), (6) в плоском и пространственном случаях эффективным оказалось использование в качестве искомых функций насыщенности S и "приведенного" давления Р , определяемого формулой

P = (7) о

В терминах функций Я , Р система (I) - (4), (6) редуцируется к системе двух уравнений /3/ lw(a(s)VS +i„vP) (8) Н

MKCs)vP)=o (9) где ,

В силу свойств функций относительных фазовых проницаемостей (5) справедливо

К>о , 5е[од] (ю)

CL(0) =01(1)-О, CL>0 , S^CO.l) (П)

Таким образом, система уравнений (8), (9) представляет собой квазилинейную систему составного типа, состоящую из равномерно эллиптического уравнения для Р и вырождающегося параболического уравнения для S . Фильтрационное течение рассматривается в заданной ограниченной, вообще говоря многосвязной, области i = с кусочно-гладкой границей

Пусть й^й x(0/T)fQ - замыкание , И/ - внешняя нормаль к Г . Граница Г включает в себя три компоненты: Гр ,

Г г

1 ± » 1 I

Компонента Г^ состоит из внешней границы зоны тампонажа, существовавшей к начальному моменту времени (контур питания),

Г^ соответствует эксплуатационным скважинам, a FJ, состоит из непроницаемой поверхности.

Перейдем к постановке краевых условий на Г . На Г0 задается условие непротекания вытесняемой смачивающей фазы и смеси, имеющие вид

0, Гот = Г>[0,Т] (к)

I от

На П, , согласно /77/, задается условие отсутствия потока вытесняемой фазы tfjtln =0, rCT-rt"Co,n (и

I ДГЦ и расход или давление вытесняющей фазы хД)еГ,

ZT = coit^-k >0

15)

Заметим, что в случае, когда на задано условие (14), г( можно включить в Г.* • ат • от

Наиболее сложной является формулировка краевого условия на йце Леверетт /71/ отмечал наличие так называемого концевого эффекта при вытеснении одной жидкости другой. Сущность его кратко состоит в следующем.

Пусть смачивающая жидкость вытесняется несмачивающей. При подходе несмачивающей вытесняющей фазы к выходному сечению в силу наличия капиллярного скачка давление больше на стороне несмачивающей жидкости. Однако после выхода последней в свободное пространство кривизна границы раздела между жидкостями невелика и давления в фазах выравниваются. Следовательно из условия непрерывности давлении слева и справа от выходного сечения следует равенство давлений в фазах на выходе, фтоге это приводит к накоплению вытесняемой смачивающей фазы в области прилегающей к выходу. Аналогичная ситуация имеет место и при вытеснении несмачивающей жидкости смачивающей. Это явление экспериментально исследовалось Перкинсом /76/, Кайтом и Рапопортом /69/, Батики и др. /62/. Численные расчеты задач фильтрации двухфазной жидкости с учетом "концевого" эффекта приведены в работах /13, 58, 65, 66, 72/.

Однако следует отметить, что во всех этих работах главным образом рассматривалось вытеснение несмачивающей фазы смачиваю

- ю щей.

Теоретическое исследование стационарной задачи о вытеснении смачивающей жидкости несмачивающей с учетом концевого эффекта дано В.И.Пеньковским /37/ . Это решение, получаемое при изучении соответствующей краевой задачи, обнаружено ранее С.Н.Бузиновым /II/ как класс частных решений стационарных уравнений двухфазной фильтрации. Некоторое обобщение решения Бузинова на случай фильтрации жидкостей с вязко-пластическими свойствами дано в работе /32/.

Таким образом, при совместном истечении фаз краевое условие на выходе имеет вид

Р^т'Р^т"*"60"-*^0 (16)

Дяя насыщенности ставится начальное условие вида s(x,o) = 5°(x), хб «0, хеГ, (17)

В терминах функций & , Р краевые условия (12) - (16) примут вид

KvP-kL =(clvs lHvP)-Tt|p = о (I8)

1 от I от

-Plrv-pc CI9)

20) civs - i,vP)-£ I p = V(x Д € Г2т (21)

ГГ или

I ат (22)

1 лт 0

Здесь VifcyO " заданная функция.

Начально-краевую задачу для системы (8), (9) с условиями (17) - (21) назовем задачей с заданным суммарным расходом, а с условиями (17) - (20), (22) - задачей с заданным перепадом давлений. Отметим, что, по терминологии авторов / 3/, задача (8), (9), (17) - (21) относится к основным начально-краевым задачам двухфазной фильтрации несмешивающихся жидкостей, следовательно к ней применимы результаты полученные в /3/ для этого класса задач.

Задача с заданным перепадом в данной постановке ранее не изучалась.

Для случая одномерного течения система (8), (9) сводится к одному уравнению для насыщенности имеющего вид где

V(t) = - К v Р - суммарная скорость фильтрации, которая в зависимости от граничных условий является заданной величиной или функционалом от S . Начально-краевые условия (17) ~ (22) при одномерном движении специализируются к виду

H'VSJUt-o <*> ж,0) = S'Cx), S ° (о) = 0. <25>

Если выполнено условие (21), то V= ~Q("t) » гДе Q ^ 0 - известная функция. При условии (22) суммарный расход определяется следующим образом. Т.к. 1/+ l^ = V("t) , то из формул (3), (4), (6) следует c ol£ Dx

Интегрируя последнее равенство по ос , с учетом краевых условий (13) и (16), получим vw=-(AP-F(S(m)))/3(0 (26) где Р = ;]« = [р(5)^х^нК"!дрнро-р,>0.

В работе также будет рассматриваться случай осесимметричного притока к скважине, поэтому проведем необходимые преобразования системы уравнений осесимметричной фильтрации, которая включает в себя (6) и i^u L& a ^-t dx\ UtrW

27)

Здесь •x-'fc2', *Ъ - радиус-вектор, отсчитываемый от оси скважины до границы зоны тампонажа, - радиус скважины. Складывая уравнения (27) и аналогично одномерному случаю проводя соответствующие преобразования, окончательно нетрудно получить уравнение дяя насыщенности в виде

4 1)t ^xV ^х v v где в случае заданного перепада давлений суммарная скорость определяется из соотношения

Здесь 3*.= Jj(5)/xix.

Легко видеть, что уравнения (23), (28), функционалы (26), (29) и начально-краевые условия (24), (25) для случаев одномерной и осесимметричной фильтрации можно записать в следующей общей форме s(*,o) = S°(x) , (33) где ^ = L = 0 в случае одномерной фильтрации, t - i , ^ > О при осесимметричном притоке, 3- = ] p(£)/xLaX. Подразумевается, что осуществлено преобразование -L.

Здесь уместно привести некоторые результаты, относящиеся к исследованию начально-краевых задач фильтрации несмешивакхцихся жидкостей.

Г.В.Алексеев и Н.В.Хуснутдинова /2/ рассмотрели первую краевую задачу при заданных насыщенностях и задачу Коши для уравнения (23) при известном суммарном расходе. Ими показано, что обобщенные решения этих задач являются непрерывными функциями. Аналогичный вопрос для случая заданных перепада давлений и насыщенностей исследовался в работе /50/.

В плоском и пространственном случаях начально-краевые задачи двухфазной фильтрации исследовались в работах С.Н.Антонцева и

B.Н.Монахова (см. в /3/ подробную библиографию), С.Н.Кружкова и

C.М.Сукорянского /28/, /29/, в /30/, /35/.

Для уравнений вида (23) корректность постановки ряда начально-краевых задач, дифференциальные свойства обобщенных решений этих задач изучались Аронсоном /60/, /61/, Кружковым /30/, Олей-ник О.А. и др. /36/, Калашниковым А.С. /22/. Частные решения ряда вырождающихся задач двухфазной фильтрации строились в /7, 39, 41, 42, 43/. Теорема единственности обобщенного решения в случае первой краевой задачи для "приведенного" давления и насыщенности анонсирована в /29/. Единственность и устойчивость обобщенных решений задач двухфазной фильтрации в данной выше постановке (заданных на границе области фильтрации давлений и/или потоков фаз) была доказана лишь в некоторых частных случаях: либо наперед заданной во всей области суммарной скорости фильтрации, либо в случае, когда класс решений таков, что функцияcC^SCx^t)) суммируема с некоторой положительной степенью, т.е. при наличии вырождения на множестве меры нуль /3/.

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

В последние три десятилетия, в связи с большой практической важностью, интенсивно развиваются методы численного решения задач двухфазной фильтрации. По численным расчетам на ЭВМ разнообразных задач фильтрации многофазной жидкости имеется весьма обширная библиография. В качестве примера можно привести сборники /54/ - /56/. Наряду с многочисленными работами прикладного характера' существует ряд работ, посвященных теоретическому обоснованию решения задач двухфазной фильтрации конечно-разностным методом, например, /3/, /28/, /25/, /46/.

Различные трудности, возникающие при численной реализации этих задач, достаточно полно освещены в работах А.Н.Коновалова /26/, /27/.

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

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

Таким методом, применяемым в настоящей работе, является метод последовательной смены стационарных состояний, развитый И.АЛар-ным /52/ и затем усовершенствованный Г.И.Баренблаттом /9/ в более общий метод интегральных соотношений. Расчеты по этому методу автомодельных задач двухфазной фильтрации проводились в /31, 51 /.

В данной работе метод последовательной смены стационарных состояний впервые применен к расчету одноразмерных неавтомодельных нестационарных задач двухфазной фильтрации.

Довольно длительную историю имеет использование метода вибровоздействия в технологии добычи нефти и газа. Одно из самых ранних упоминаний сделано в /73/. Подробное описание различных аспектов использования этого метода дано в /12/. Исследование влияния акустических или ультразвуковых полей, получаемых цри вибровоздействии, на проницаемость пористой среды проводилось в работах /49/, /53/, /67/.

Основы математического моделирования с целью объяснения влияния способа вибровоздействия на фильтрацию не см вшивающихся жидкостей в призабойной зоне скважины заложены в /38/.

В настоящей работе излагаются результаты дальнейших исследований по построенной в /38/ математической модели.

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

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

Второй параграф главы I посвящен исследованию корректности постановки одномерной и осесимметричной задач вытеснения цри заданном перепаде давлений. Доказаны теоремы существования обобщенного решения, а для частного вида начальной функции и единственности, обнаружен ряд свойств обобщенного решения. Например, следует отметить, что принцип максимума имеет вид, отличный от доказанного ранее /3/, и это отличие состоит в том, что предельное значение насыщенности сверху определяется величиной перепада давлений.

В § 3 этой же главы изучены различные качественные свойства обобщенных решений задач при заданных перепаде давлений или суммарном расходе в случае одноразмерных течений, обнаруженные в ходе дальнейших исследований.

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

Если два расхода Q1(t) и связаны неравенством

G^ ^ Q^ , то аналогичное неравенство имеет место и дяя соответствующих обобщенных решений. Далее, если функция Q(-fc.) является монотонно неубывающей функцией времени, то и обобщенное решение есть монотонно неубывающая функция времени. В качестве следствия приведена оценка решения нестационарной задачи через решение стационарной. Выяснено влияние функции капиллярного давления на структуру обобщенного решения. Получены оценки сверху и снизу для обобщенных решений через некоторые известные функции.

Следует отметить, что для уравнений нелинейной нестационарной однофазной фильтрации аналогичные теоремы сравнения получены В.М.Ентовым в работе /19/.

В § 4 для решения вырождающейся одномерной задачи вытеснения при известном суммарном расходе доказано свойство стабилизации при бесконечном возрастании времени к стационарному решению с оценкой скорости сходимости.

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

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

В § 3 дана экспериментальная проверка полученных теоретических результатов. В нем изложена методика проведения эксперимента, приведены его результаты и дано сравнение экспериментальных данных с результатами численных расчетов.

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

В § 3 посредством анализа стационарных решений для случая кусочно-однородной пористой среды изучена еще одна возможность устранения капиллярного запирания фильтрата бурового раствора.

В четвертом параграфе главы 3 на основе полученных стационарных решений оценено влияние сжимаемости вытесняющей фазы на расход этой фазы и распределение насыщенности в режиме капиллярного запирания вытесняемой фазы.

Материалы диссертации изложены в работах / 5, 16-18/, а также ряде отчетов по НИР, выполненных в лаборатории фильтрации Института гидродинамики СО АН СССР им. М.А.Лаврентьева.

Результаты диссертации докладывались на республиканской конференции по нелинейным задачам математической физики (г.Донецк,

1983), на семинаре теоретического отдела Института гидродинамики СО АН СССР под руководством члена-корреспондента АН СССР Л.В.Овсянникова, на семинаре "Краевые задачи механики сплошной среды" в Институте гидродинамики СО АН СССР под руководством профессора В.Н.Монахова; на семинаре "Большие задачи математической физики" в Новосибирском госуниверситете под руководством профессора А.Н.Коновалова, на семинарах лаборатории фильтрации Института гидродинамики СО АН СССР, на семинаре по механике жидкости и газа под руководством академика АН АзССР А.Х.Мирзаджанзаде и академика С.А.Христиановича в Институте проблем механики АН СССР (г.Москва),

Автор выражает искреннюю благодарность д.ф.-м.н. С.Н.Антонце-ву и к.ф.-м.н. В.И.Пеньковскому за научное руководство и поддержку, а также помощь в работе.

- 103-ЗАК1ЮЧЕНИЕ

Результаты проведенного исследования:

1) Сформулирована математическая модель для описания фильтрации несмешивающихся жидкостей в призабойной зоне скважины и исследована корректность постановок соответствующих начально-краевых задач.

2) В одномерном и плоском осесимметричном случаях обнаружен ряд качественных свойств обобщенных решений задач вытеснения смачивающей фазы несмачивагацей в однородной пористой среде с учетом выходного концевого эффекта.

3) Предложен способ численного расчета и проведено численное исследование одномерной и плоской осесимметричной задач вытеснения смачивающей фазы несмачивающей при заданных перепаде давлений и суммарном потоке фаз. Изучено влияние различных физических параметров на характер процесса вытеснения.

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

5) На основе известной математической модели проведено численное исследование воздействия виброисточника на призабойную зону скважины.

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

7) Изучено влияние сжимаемости вытесняющей фазы на характеристики фильтрационного течения в режиме капиллярного запирания вытесняемой фазы.

8) Проведенное исследование позволяет дать некоторые рекомендации по оптимальному выбору параметров промывочных жидкостей: а) Необходимо применять промывочные жидкости с наименее возможным значением коэффициента поверхностного натяжения их фильтрата. б) Повышения дебита скважины можно достичь использованием таких промывочных жидкостей, фильтрат которых изменяет смачивающие свойства пород в окрестности скважины в сторону их гидрофобиза-ции. в) Вязкость фильтрата бурового раствора должна быть минимально; возможной.

1. Абрашин В.М. Разностные схемы для параболических уравнений с нелинейным вырождением. - Дифференц.уравнения,1976, т.12, гё 8, с.1470-1484.

2. Алексеев Г.В., Хуснутдинова Н.В. О разрешимости первой краевой задачи для уравнения одномерной фильтрации двухфазной жидкости. ДАН СССР, 1972, т.202, гё 2, с.310-312.

3. Антонцев С.Н., Кажихов А.В., Монахов В.Н. Краевые задачи механики неоднородных жидкостей. Новосибирск: Наука,1983, - 315 с.

4. Антонцев С.Н. О локализации решений нелинейных вырождающихся эллиптических и параболических уравнений. ДАН СССР, 1981, т.260, 1Гз 6, с.1289-1293.

5. Антонцев С.Н., Доманский А.В. Единственность обобщенных решений вырождающихся задач двухфазной фильтрации. Новосибирск: Численные методы механики сплошной среды, 1984, т.15, № 6, с.15-27.

6. Бан А. и др. Влияние свойств горных пород на движение в них жидкости. М.: Гостоптехиздат, 1962, - 257 с.

7. Баренблатт Г.И. Об одном классе точных решений плоской одномерной задачи нестационарной фильтрации газа в пористой среде. ПММ, 1953, т.17, № 6, с.739-742.

8. Баренблатт Г.И., Ентов В.М., Рыжик В.М. Теория нестационарной фильтрации жидкости и газа. М.: Недра, 1972, - 287 с.

9. Баренблатт Г.И. О некоторых приближенных методах в теории одномерной неустановившейся фильтрации жидкости при упругом режиме. Изв.АН СССР, ОТН, 1954, № 9, с.35-50.

10. Бэр Я., Заславски Д., Ирмей С. Физико-математическое основы фильтрации воды. М.: Мир, 1979, - 451 с.- 106

11. Бузинов С.Н. К вопросу об определении остаточной перенасыщенности. ДАН СССР, 1957, т.116, Ш I, с.28-31.

12. Гадиев С.М. Использование вибрации в добыче нефти. М.: Недра, 1977, - 200 с.

13. Гусейнов П.М., Яницкий П.А. Концевой эффект при двухфазной фильтрации в условиях малой начальной обводненности. Изв. АН СССР, МЖГ, 1979, гё I, с.61-66.

14. Данилов В.Л., Кац Р.С. Гидродинамические расчеты взаимного вытеснения жидкостей в пористой среде. М.: Недра,1980, -264 с.

15. Данилов В.Л., Коновалов А.Н., Якуба С.И. Об уравнениях и краевых задачах теории двухфазннх фильтрационных течений в пористой среде. ДАН СССР, 1968, т.183, № 2, с.307-310.

16. Ентов В.М. Теоремы сравнения для уравнений нестационарной фильтрации. М.: ПММ, 1965, т.29, № I, с.200-205.- 107

17. Зубов Н.В., Цыбульский Г.П. Новый метод решения одномерной задачи о вытеснении нефти водой с учетом капиллярных сил. -Изв. АН СССР, ЖГ, 1975, 3, с.83-91.

18. Индельман П.В., Кац P.M. Численное моделирование двухфазной фильтрации с учетом концевого эффекта. Изв.АН СССР, MKT, 1979, № I, с.175-179.

19. Калашников А.С. 0 непрерывной зависимости обобщенных решений уравнения нестационарной фильтрации от функции, задающей режим потока. ПММ, 1978, т.42, № I, с.183-185.

20. Калашников А.С. Об условиях положительности и ограниченности решений второй краевой задачи для параболических уравнений. -Сиб. мат.журнал, 1981, т.22, II 6, с. 72-80.

21. Коллинз Р. Течение жидкостей через пористые материалы. М.: Мир, 1964, - 349 с.

22. Коновалов А.Н. О некоторых вопросах, возникающих при численном решении задач фильтрации двухфазной несжимаемой жидкости. В кн.: Тр.Мат.ин-та АН СССР,1973, т.122, с.2-23.

23. Коновалов А.И. Задачи фильтрации многофазной несжимаемой жидкости. Новосибирск; изд. Новосиб.ун-та, 1972, - 128 с.

24. Коновалов А.Н., Коробицына Э.Л. К обоснованию метода расщепления по физическим процессам для задач фильтрации. Новосибирск; ИТПМ СО АН СССР, Препринт, 1978, - 21 с.

25. Кружков С.Н., Сукорянский С.М. Краевые задачи для систем уравнений двухфазной фильтрации: постановка задачи, вопросы разрешимости, обоснование приближенных методов. Мат. сб., 1977, т.104, Ш I, с.69-88.

26. Кружков. С.Н. О некоторых вопросах теории нелинейных дифференциальных уравнений связанных с задачами фильтрации. Успехи матем.наук, 1978, т.33, № 3, с.148.

27. Кружков С.Н. Результаты о непрерывности решений параболических уравнений и некоторые их применения. .- Матем.заметки, 1969, 6, № I, с.97-107.

28. Куванышев У.П. 0 фильтрации двухфазной жидкости в пористой среде с учетом фазовой проницаемости и капиллярных сил. -Тр.Тат.нефтьНИИ, 1964, вып.6, с.321-332.

29. Кудрявцев Г.В., Фаткуллин А.Х. Расчет остаточной насыщенности в неоднородном пласте при взаимном вытеснении жидкостей с вязко-пластическими свойствами. -ПМТФ, 1970, J£ 4, с.189-192.

30. Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического типа. М.: Наука, 1967, - 736 с.

31. Лейбензон Л.С. Движение.природных жидкостей и газов в пористой среде. М.-Л.: Гостехиздат, 1947, - 244 с.

32. Мороз П.А., Сукорянский С.М. Некоторые проблемы решения задач двухкомпонентной фильтрации жидкости в пласте. М.: Науч.тр.Ин-та автоматики, М-во приборостсредств автомат, и систем упр.СССР, 1975, J& 7, с.75-86.

33. Олешик О.А., Калашников А.С., Чжоу Юй-Линь. Задача Коши и краевые задачи для уравнений типа нестационарной фильтрации. Изв. АН СССР, 1958, серия матем.,т.22,№ 5,с.667-704.

34. Пеньковский В.И. Концевой эффект капиллярного запирания вытесняемой фазы при фильтрации несмешивавдихся жидкостей. -Изв.АН СССР, ЖГ, 1983, № 6, с.184-187.

35. Пеньковский В.И. Модель процесса очистки газоносного пласта вибровоздействием. В кн.: Механика быстропротекащих процессов, Новосибирск, изд.Ин-та гидродинамики СО АН СССР,1981, с.84-93 (Динамика сплошной среды, вып.51).

36. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. М.: Наука, 1977, - 664 с.- 109

37. Развитие исследований по теории фильтрации в СССР/Под редакцией П.Я.Полубариновой-Кочиной. М.: Наука, 1969, - 545 с.

38. Рахимкулов Й.Х., Швидлер М.И. Автомодельная задача о совместном движении нефти и воды. Изв. АН СССР, ОТН, 1962, № 2,с.136-137.

39. Рыжик В.М., Чарный И.А., Чэнь-Чжун-Сян. О некоторых точных решениях уравнений нестационарной фильтрации двухфазной жидкости. Изв.АН СССР, ОТН, Мех. и машиностр.,1961, Ш I, с.121-126.

40. Рыжик В.М. О капиллярной пропитке водой нефтенасыщенного гид-фильного пласта. Изв.АН СССР, ОТН, Мех. и машиностр., I960, № 2, с.I49-I5I.

41. Самарский А.А. Теория разностных схем. М.: Наука, 1977, - 636 с.

42. Соболев СЛ. Некоторые применения функционального анализа в математической физике. Новосибирск: изд.СО АН СССР,1962, -265 с.

43. Сукорянский С.М. О сходимости одного численного метода решения уравнений одномерной двухкомпонентной фильтрации жидкости в нефтяном пласте. Научн.тр.Ин-та автоматики. М-во при-боростр.средств автомат, и систем управ.СССР,1975,вып.7,с.87-92.

44. Фридман А. Уравнения с частными производными параболического типа. М.: Мир, 1968, - 427 с.

45. Царев В.П. Восстановление проницаемости коллекторов, затам-понированных пресными промывочными жидкостями. Газовая промышленность, 1976, JU? 12, с.21-24.

46. Цхай А.А. 0 разрешимости одной задачи одномерной фильтрации двухфазной жидкости. В кн.: Математические проблемы гидродинамики, Новосибирск, изд. Ин-та гидродинамики СО АН СССР, 1983, с.173-179 (Динамика сплошной среды, вып.59).

47. Чарный И.А. Подземная гидрогазодинамика. М.:1963, - 396 с.

48. Чарный И.А. Метод последовательной смены стационарных состояний и его приложения в задачам нестационарной фильтрации жидкостей и газов. Изв. АН СССР, ОТН, № 3, 1949, с.323-343.

49. Черский II.В. и др. Влияние ультразвуковых полей на проницаемость горных пород при фильтрации воды. ДАН СССР, 1977,т.232, Jfe I, с.201-204.

50. Численные методы решения задач фильтрации несжимаемой жидкости/Под ред.А.Н.Коновалова.-Новосибирск: изд.ВЦ СО АН СССР, 1972, 227 с.

51. Численные методы решения задач фильтрации несжимаемой жидкости / Под ред.А.Н.Коновалова. Новосибирск: изд. ВЦ СО АН СССР, 1975, - 318 с.

52. Численное решение задач фильтрации многофазной несжимаемой жидкости /Под ред.А.Н.Коновалова. -Новосибирск: изд.ВЦ СО АН СССР, 1977, 239 с.

53. Шейдеггер А.Э. Физика течений жидкостей через пористые среды.-М.: Гостоптехиздат, I960, 249 с.

54. Швидлер М.И., Леви Б.И. Одномерная фильтрация несмешивающихся жидкостей. М.: Недра, 1970, - 155 с.59. tackle. G.E. Ehtizical ResU"tiVLtij Lotj otfi Aid in, 0&^егтл)гии| Soiue Яе-бегУо1г -TW>. AIME, 1942, v. 4A6, p.54-61.

55. A^on-m D.&. Re<ju£axltij ptopetiles of tflovJs -ttaoujh, pozoiAb №aAlol.-SIAM З.АррВ. Ми/IL,

56. АгоибОYi D.&. RegiJ!aK,L"hj pzopeciu-b ^loW-Ь -th/towcjkрогомл vneotbof. -tta Lvule*,|cice,.-A^ck. RaAlon. Meek,. cHrwl Amu?., 1370, v. 3?, A/1 ? p. 1 -10.

57. BctAye^ R C.^Hoolgut* R Reuben, D. B.1.rLenptA-ktvijj taPaAiVe petmeafcfc-tg avidL We-t-toL^L^ <jW uvi-yiecxdij --ytcdre di-bp^OLceVYien/t h^O/bUrtemm-t-b. Soc. Pe-ttoC. Engwj.3., 1381, v.21,^ p.296-30<?.

58. Bi>ucfc8ej 3., LeveJce-tt M. С. MeckcuU-sm, F-CwLfll Du-p&KuewieH/f: Lvb

59. Donj£eu> DJt., В£а1г P. M., VonjwJt, R.l Coi-tcM^a-Lon, of. Luweoac Vlaienflood Behave I^e&AcLlvig-Lk E ffect of

60. СирМаъц Рге^иле-Тгап^. А1МБ,19«,у.&15,р.96-10г.

61. Расгои?а»ьЯб H.V., Скеи, W.I. Ivtf<£иелле o|upon, It^uuL |£oW tWwcjk potow-ь vvtedU.-WJ-li'tci-bovtLC-^ 1969, v.?, p. 135 -199.

62. Pe-Uo£. TecUof., N/10 , p. И -50.

63. Leve^ett M.t. P«£ow/ of -WoiteJt M ix Web -tko-ugk, Цисоп-ьо fioloited Soividl^.-^n^ А1МЕ,1339лу.132, p1*1.

64. Le-Vett-eA-t M.C. Capl££oouj BekotVuOt ivu Рог01Л-Ъ SofcuJ*.-TW*. А1МЕ,1Н1,У.Ш,р.152-169.

65. Rc. Hwevb C.R. A tfuYYiettcoi^ 9>o£wti,oVb of ike LimcUl Dl-b-рВмсгмлиЛ: E^unKovt'Witk р^г^ияс. -Ttom^. AIME , 1953, W.73.(1ei:2getL H.A. Р^<£и\гоыхяу Repoit on, Ftoictwti.ttcj ^

66. V^iDLitovb AN'ew Me-tkocL of SttrniA-Eatim-3. РеЛгоЬ. TeckW., 195£; Д/ll, p. 13-16.

67. Моо1& Т.£<to£odL R.L. Tk Effect

68. CaptWcotlttj ok tkfi. Di-6p^0Lceyv?evut of Oil ify Y/ni&i,.-Ptod. Mountb£j; 195g, v.jeo, л/±0^ p.20-30.

69. Muvigowt V. In+e^foicUi Effects iv\, Lvmu-bctUfe. Liquid Lic^wlcl DL-*>p£ace,wuLyi/t Ы Рогои-ъ M^ckoi .О. Soc. Ре1го£. Ewj^. P.M. Avo, InventL^a-tLou, of tk.e Ro£e of CapLHcu^рогса^ LK LaUotato'Uj

70. Wat en. Р€оо<А*.- 3. Petfcot. TeeUo?.,1./11, p. 49 -51.

71. Papopot4: L. A v Lea* W.l Ptop^Jttte-b Linecat Waten-ffcoocl*.-Тгакб. AIM E , 1953 , v. 19*.78$lckajicUokt, 3.G-. Л a£. Lai^cx-^ij dete'tMxn.oi-kLow, of t^Pitlve-Тъяи*. AIM E,195fc, V.19S, p.l8*-I9G.

72. Wotgvi£R. 0.R. EffecA of Intend с Tev\>> itm, оуь DLvptauL

73. YyuzmA. Eff:clevui^.- Soc. Ре1го£. Eyijmj.I, 19CG, V. С |\M.

74. V/gcfcoff R.D. fBot*e.t H.R TV -f£ow of Goi* Luj,uU Mux-tirn-ь tkjiougVb ITncon-bO'EudlQttecL Somd 6. Pk^-bi-c^ 1936^ v.? ,hld 325 - 3^5*.о- 113 -0.2 0.4 х. 0.6 0.80.208 1.00.21. Рис.104 ос. 0.60.204 В 0.60.8