Математические модели влияния особенностей призабойных зон на фильтрацию жидкости в многоскважинных системах тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат технических наук Харченко, Юрий Владимирович

• Харченко, Юрий Владимирович
• кандидат технических наук
• 2008

Оглавление диссертации кандидат технических наук Харченко, Юрий Владимирович

ГЛАВА 1. РАЗВИТИЕ МЕТОДА ХАРАКТЕРИСТИЧЕСКИХ ФУНКЦИЙ ДЛЯ ИССЛЕДОВАНИЯ ФИЛЬТРАЦИИ ЖИДКОСТИ В ПРИЗАБОЙНЫХ ЗОНАХ СКВАЖИН

1.1 Уравнения стационарной фильтрации жидкости в изотропных и анизотропных однородных и неоднородных средах

1.2 Применение метода характеристических функций к решению краевых задач фильтрации эюидкости в изотропных средах

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

1.4 Построение методом характеристических функций точного решения задачи фильтрации жидкости к круговой скважине в пласте с прямолинейной анизотропией при естественном контуре питания

1.4.1. Комплексные потенциалы течения

1.4.2 Характеристическая функция течения к круговой скважине в среде с прямолинейной анизотропией и точная формула расчета дебита при естественном контуре питания

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

1.5.1 с естественным контуром питания и гравийным фильтром

1.5.2 с центральным круговым контуром питания

1.5.3 с прямолинейным контуром питания 47 Основные результаты 1-ой главы

ГЛАВА 2. ИССЛЕДОВАНИЕ ОСОБЕННОСТЕЙ ФИЛЬТРАЦИИ ЖИДКОСТИ В ПРИЗАБОЙНОЙ ЗОНЕ СКВАЖИНЫ С КОНТЕЙНЕРНЫМ ГРАВИЙНЫМ ФИЛЬТРОМ 56 2.1 Простейшая плоскорадиальная модель фильтрационного потока к скважине и ее значение для практики разработки пластов

2.2 Применение плоскорадиальной модели фильтрационного потока к проектированию забоя скважины с контейнерным гравийным фильтром при линейном законе Дарси

Основные результаты 2-ой главы

ГЛАВА 3. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЛИЯНИЯ ОСОБЕННОСТЕЙ ПРИЗАБОЙНЫХ ЗОИ НА ФИЛЬТРАЦИЮ ЖИДКОСТИ В МНОГОСКВАЖИННЫХ СИСТЕМАХ

3.1 Интегральное представление Р -гармонических функций через аналитические

3.1.1 Переход к уравнению Гельмголы(а

3.1.2 Общее решение уравнения Гелъмголъца

3.1.3 Интегральное представление функций в функциональной последовательности

3.2 Особенности математического моделирования плоскопараллельной фильтрации к скважинам в пластах гармонической серии

3.3 Математическое моделирование фильтрации жидкости к скважинам в изотропных однородных и неоднородных средах методом наименьших квадратов

3.3.1 Общая постановка задачи

3.3.2 Нормированные фундаментальные решения

3.3.3 Применение метода функций Грина для моделирования течения к одиночным и групповым скважинам

3.3.4 Построение функции Грина методом распределённых источников

3.4 Числовые оценки зависимости дебитов скважин кольцевой батареи от скачка проницаемости в призабойных зонах 108 Основные результаты 3-ей главы

Введение диссертации (часть автореферата) на тему «Математические модели влияния особенностей призабойных зон на фильтрацию жидкости в многоскважинных системах»

Актуальность темы. При строительстве скважин призабойные зоны пласта по целому ряду причин приобретают новые фильтрационные свойства [1, 6, 8, 28, 30, 32, 54], отличающиеся от фильтрационных свойств пласта в местах, удаленных от скважин. К этому приводят технологические операции кислотной обработки призабойных зон, гидроразрыва пласта, гравийной обсыпки и н. др. Наиболее полно механизм изменения физических свойств горных пород в околоскважинпых зонах описан Н. Н. Михайловым в [54]. Главная из целого ряда причин, приводящая к изменению фильтрационных свойств призабойной зоны пласта, в ее кольматации [11, 16, 20, 34, 38, 42], определяемой процессом выноса с током флюида к фильтру скважины [10, 12, 18] продуктов разрушения пласта.

Наиболее сложная в конструктивном отношении часть водо-, нефте- и газодобывающих скважин относится к окончанию их ствола, располагающемуся в продуктивном пласте. Окончание скважины, принимающее продукт из коллектора и называемое активной частью ствола скважины, всегда снабжается фильтром той или иной конструкции [8, 10, 78, 79]. Главное назначение фильтра, во-первых, в том, чтобы предупредить так называемое «запескование» [3, 8, 11] ствола скважины, снижающее её эксплуатационные качества и приводящее в дальнейшем к необходимости проведения длительных и дорогостоящих ремонтных работ. Во-вторых, в том, чтобы обеспечить активную часть ствола скважины достаточно высокими «пропускными» для извлекаемого из пласта продукта свойствами.

Именно поэтому гидродинамический анализ влияния изменений околоскважинных зон пласта на производительность скважин является важным для практики разработки нефтегазовых месторождений [38]. В частности актуальной является задача развития таких методов расчета дебитов скважин, которые 1) учитывают отличие фильтрационных свойств призабойных зон скважин (ПЗС) от свойств остальной части пласта и,

2) позволяют определить максимально допустимые дебиты скважин, при которых риск кольматации ПЗС и запескования скважин минимален.

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

Литературная справка. Несмотря на то, что, как уже сказано выше, ПЗС имеют фильтрациоиные свойства, отличающиеся от свойств пласта, в расчетах дебитов скважин это отличие, как правило, не учитывается. Это объясняется рядом причин. Во-первых, даже в простейшем случае плоскопараллельной фильтрации жидкости в однородном изотропном пласте, учет скачка проницаемости в ПЗС приводит к необходимости отдельного определения комплексных потенциалов течения за пределами ПЗС и внутри ПЗС со сложными граничными условиями сопряжения на границах ПЗС. Сложность и трудность решения задач такого рода можно видеть, например, по работе М. А. Лукомской [50], которая рассматривала задачи о притоке жидкости к скважинам в кусочно-однородных изотропных средах, и, в частности, случай скачка проницаемости на круговых и прямолинейных границах. Точные решения в [50] приведены лишь для случая, когда в области фильтрации есть одна, две границы раздела сред с разными проницаемостями в виде прямой или окружности. Задачам фильтрации в кусочно-однородных изотропных средах уделяется внимание в работах многих других авторов - О. В. Голубевой [27], В. П. Пилатовского [59], П. Я. Полубариновой-Кочнной [63], И. А. Прусова [64], В. А. Толпаева [82] и др. Но во всех перечисленных работах, за исключением [82], речь идет о задачах фильтрации в кусочно-однородных средах при наличии одной или двух границ раздела в виде или прямых или окружностей. Даже в этом случае, когда имеется только две границы раздела в виде окружностей, задача вызывает большие математические трудности, в чем можно убедиться по работам, посвященным фильтрационным теоремам о двух окружностях В. М. Радыгина [68] и А. Г. Ярмицкого [110]. В задаче же учета особых фильтрационных свойств призабойных зон группы скважин таких границ раздела со скачками проницаемости не одна, две, а значительно больше. Поэтому поиск точного аналитического решения такой задачи является практически безнадежным делом. Но для достаточно точного приближенного решения задачи об учете скачков проницаемости в ПЗС можно привлечь идею В. А. Толпаева [82] моделирования границ раздела кусочно-однородных сред эквипотенциалями течения. Именно по этому пути и пошел диссертант исследуя фильтрацию к кольцевой батарее скважин, в призабойных зонах каждой из которых есть скачок проницаемости.

Исследованиям фильтрации в анизотропных и кусочно-однородных анизотропных средах посвящено значительно меньшее количество работ. До сих пор не найдено точное аналитическое решение задачи Дюпюи для однородной среды с прямолинейной анизотропией. Приближенное решение задачи Дюпюи для прямолинейно-анизотропной однородной среды приведено в работах Е. С. Ромма [71], 10. Л. Соломко [76] и В. А. Толпаева и А. Д. Жернового [87]. Поэтому расчеты дебитов скважин в анизотропных пластах представляют актуальную проблему. Диссертант для расчета дебитов скважин в анизотропных пластах выделяет «локал» - призабойную зону, ограниченную круговым контуром скважины и некоторой близ расположенной эквипотенциалыо, называемой в диссертации естественным контуром питания, или иначе, границей естественной депрессионной воронки. Затем соискатель осуществляет «сращивание» решения в призабойной зоне в виде «локала» с решением за пределами «л окал а». На этом пути соискателю удалось решить ряд новых задач фильтрации жидкости в среде с прямолинейной анизотропией.

Важность учета при расчете течений особых фильтрационных свойств ПЗС определяется и тем, что ряд технологических операций, направленных на интенсификацию пефте- и газодобычи, весьма сильно изменяет проницаемость и пористость призабойных зон. В частности, сильно меняют проницаемость ПЗС операции кислотной обработки пласта [4], гидравлического разрыва [28], проникновение фильтрата при промывке пеной песчаной пробки скважины [4], оборудование скважинных фильтров гравийной обсыпкой [1, 6, 10, 30, 32, 36]. Исследованиям фильтрации к скважинам с трещинами гидроразрыва посвящены работы ряда авторов. В случае стационарной фильтрации однородной жидкости - В. М. Ентова и В. В. Мурзенко [31], в случае вертикальной трещины большой протяженности - А. Ф. Зазовского и Г. Т. Тодуа [37], в случае вертикальной эллиптической трещины гидроразрыва - В. В. Кадета и В. И. Селякова [41] и др. Исследованиям фильтрации к скважинам с гравийной обсыпкой фильтров специального внимания не уделялось, так как контейнерные гравийно-намывные фильтры [30] и гравийные фильтры [78] стали широко применяться в практике сравнительно недавно, примерно с 2000 г. Последнее и определило задачи второй главы диссертации.

Задачам фильтрации в неоднородных средах посвящено большое число работ. Наиболее важными в этом направлении являются работы Г. В. Голубева [23-25], Г. В. Голубева и Г. Г. Тумашева [26], В. Ф. Пивень [57, 58], И. А. Прусова [64] п н.др. Наибольшая сложность построения аналитических решений задач фильтрации в неоднородных изотропных средах в том, что для эллиптического уравнения фильтрации с переменным коэффициентом проницаемости бывает неизвестно первое фундаментальное решение [47, 80], описывающее течение к точечным особенностям типа сток (источник) и вихрь. Для построения первого фундаментального решения можно, в принципе, использовать известные [63] интегральные представления решений эллиптического уравнения с переменным коэффициентом проницаемости через аналитические функции комплексного переменного. Однако в случае переменного коэффициента проницаемости общего вида известными интегральными представлениями воспользоваться, ввиду сложности построения соответствующего ядра интегрального оператора, практически, не удается. Диссертанту удалось указать достаточно широкий класс переменных коэффициентов проницаемости, для которых названное интегральное представление удается построить значительно проще в духе методов С. Бергмана [13] и Г. И. Назарова [55, 56]. Авторской модификации методов С.Бергмана и Г.И.Назарова посвящен п. 3.1. диссертации. В частном случае, для так называемой «гармонической серии» законов изменения проницаемости, указанное автором интегральное представление сводится к известной [27] формуле перехода, позволяющей сводить задачи фильтрации к граничным задачам для уравнения Лапласа. На последнем пути автор исследовал возможности приближенных постановок краевых задач фильтрации в неоднородных средах гармонической серии.

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

Целью работы является развитие таких методов расчета дебитов скважин, которые 1) учитывают отличие фильтрационных свойств ПЗС от свойств остальной части пласта и, 2) позволяют определить максимально допустимые дебиты скважин, при которых риск кольматации призабойных зон и запескования скважин будет незначителен.

Научная новизна и теоретическое значение работы определяются следующим:

1. Предложен новый метод математического моделирования фильтрации к скважине в изотропных и анизотропных средах на основе характеристических функций. С его помощью дана новая гидродинамическая интерпретация (р,д) — аналитических функций Г. Н. Положего, направленная на решения задач фильтрации в однородных средах с прямолинейной анизотропией.

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

3. Проведено подробное исследование простейшей плоскорадиальной модели фильтрационного потока. С помощью плоскорадиальной модели: 1) исследованы возможные относительная ошибка расчёта удельного дебита скважины, вызванные неточностью определения проницаемости пласта. 2) Предложен практический способ аппроксимации функции проницаемости к(г) радиально-неоднородной призабойной зоны скважины.

4. Предложена элементарная теория проектирования оптимального забоя скважины с контейнерным гравийным фильтром (КГФ). В частности:

1) выведены основные расчетные формулы для исследования линейной напорной фильтрации в призабойной зоне скважины с КГФ.

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

5. Выведена формула для интегрального представления потенциала скорости фильтрации в изотропных средах с коэффициентом проницаемости к(х,у) из серии решений уравнения Л{лЩx^)-g(x)■^¡Цx^) = 0 с произвольной наперёд задаваемой функцией (Заметим, что по аналогии с представленным в диссертации методом может быть сформулирован и другой результат - для случая, когда g = g{y) и л{^к[х,у))-g{y)-^k{x,y) = 0).

6. Предложена модификация известного в вычислительной электрофизике метода [39], который в диссертации назван методом распределенных источников. Метод распределенных источников в диссертации он доведён до практического алгоритма, сводящегося к решению СЛАУ с выведенными для неё коэффициентами.

7. Выведена формула для расчёта дебита отдельной скважины в круговой батарее из п скважин, учитывающая, что в призабойной зоне каждой скважины проницаемость к, не совпадает с проницаемостью к0 всего пласта. В частном случае к, = к0 из выведенной формулы вытекает известная формула В. Н. Щелкачёва [108]. С помощью выведенной формулы проведено исследование относительного прироста суммарного дебита скважин в кольцевой батарее, который можно получить в случае использования гравийной обсыпки фильтров скважин.

Подчеркнем, что развитый в диссертации аппарат можно применить не только к задачам теории фильтрации, но и для исследования процессов и явлений различной физической природы, описываемых уравнениями вида V = k(x,y,z)-grad (р; divv = 0.

Достоверность и обоснованность научных положений и результатов диссертации подтверждается корректным применением апробированного математического аппарата (векторный анализ, уравнения математической физики [80], теория аналитических и обобщенных аналитических функций [17, 47, 57, 58, 61, 63, 74]) и корректностью использования в вычислительных экспериментах апробированных специализированных программных сред (Maple 7, MathCAD 8, Visual С++ 2005, Matlab 5).

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

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

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

Результаты диссертационного исследования использовались в учебном процессе в СевКавГТУ при чтении студентам специальности 073000 «Прикладная математика» курса по теории функций комплексного переменного и ДПВ «Элементы гидродинамики и краевые задачи теории фильтрации».Полученные научные результаты могут быть использованы в разработках уточненных методик расчетов дебитов эксплуатационных скважин в таких организациях как: «Институт проблем нефти и газа РАН»; ОАО «ВНИИНефть им. акад. А.П. Крылова», г. Москва; ОАО «Газпром»; ОАО «СевКавНИПИГаз», г. Ставрополь и др.

Практическим внедрением результатов диссертационной работы в научно-исследовательские и проектные организации нефтегазодобывающей промышленности России является отчет по гранту Минобразования России АОЗ-2.8-235 (№ госрегистрации 01.20.04 12455).

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

СевКавГТУ (рук. д.ф.-м.н., профессор Толпаев В.А.).

Отдельные результаты диссертации докладывались и обсуждались на:

- 1-ой (29-30 марта 2001 г.) и 3-ей (17-19 апреля 2003 г.) региональных научных конференциях студентов и преподавателей «Проблемы компьютерных технологий и математического моделирования в естественных, технических и гуманитарных науках» (СевКавГТУ, гг. Ставрополь, Георгиевск);

- 2-ой научно-практической конференции «Совершенствование методов управления социально-экологическими процессами и их правовое регулирование» 17-18 мая 2001 г. (Ставропольский институт управления, г. Ставрополь);

- Междисциплинарном научном семинаре вузов Северо-Кавказского региона «Циклы» (СевКавГТУ, г. Ставрополь, 2002);

- 4-ой межрегиональной научной конференции «Студенческая наука — экономике России» 24-26 апреля 2003 г. (СевКавГТУ, г. Ставрополь);

- 9-ой Всероссийской научно-технической конференции «Информационные технологии в науке, проектировании и производстве» (Нижний Новгород, 2003г.);

- 4-ой (апрель 2004 г.) региональной научной конференции студентов и преподавателей «Математическое моделирование и информационные технологии» (СевКавГТУ, гг. Ставрополь, Георгиевск);

- 1-ой Международной научно-технической конференции «Инфотелекоммуникационные технологии в науке, производстве и образовании», апрель, 2004 г. (СевКавГТУ, гг. Ставрополь, Кисловодск);

- 4-ой Международной научно-практической конференции «Методы и алгоритмы прикладной математики в технике, медицине и экономике» (Новочеркасск, 2004г.);

- 5-ом Всероссийском симпозиуме по прикладной и промышленной математике (г. Кисловодск, 2004г.);

- 1-ой Международной Интернет-конференции «Проблемы повышения газонефтеотдачи месторождений и ПХГ на поздней стадии эксплуатации», 2004 г. (гг. Москва, Ставрополь);

- Воронежской зимней математической школе «Современные методы теории функций и смежные проблемы», декабрь 2007 г. (ВГУ, г. Воронеж);

Публикации. По теме диссертации опубликовано в соавторстве 17 работ [89-105], в том числе отчет (номер государственной регистрации 01.20.04 12455) по гранту АОЗ-2.8-235 Минобразования России и 5 работ [100

- 103, 105] в изданиях, рекомендованных ВАК РФ.

В совместных работах соискателю принадлежат выводы расчётных формул и разработка программ для выполнения вычислительных экспериментов. Руководителю - постановки проблемных задач и общее руководство.

Структура и объем работы. Материал диссертации изложен на 130 страницах машинописного текста. Диссертация состоит из введения, трех глав, заключения и списка литературы. Содержит 23 рисунка и 18 таблиц. Библиография включает 117 наименований, из которых 7 на иностранных языках. Главы начинаются с краткого вступления, в котором перечисляются ее основные цели и задачи, и заканчиваются формулировкой основных результатов. Формулы в текущем пункте имеют одинарную нумерацию. При ссылке на формулу из другого пункта применяется традиционная тройная нумерация: номер главы, номер пункта и номер формулы в пункте.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Харченко, Юрий Владимирович

8. Основные результаты диссертации внедрены в практику проектирования разработки нефтяных месторождений в виде научного отчета по гранту Минобразования России АОЗ-2.8-235 на тему «Математические модели влияния особенностей призабойных зон на фильтрацию жидкости в многоскважинных системах» (№ госрегистрации 01.20.04 12455).

ЗАКЛЮЧЕНИЕ

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

К основным результатам диссертации относятся следующие.

1. Предложен новый метод математического моделирования фильтрации однофазной несжимаемой жидкости к скважине в изотропных и анизотропных средах на основе метода характеристических функций. С его помощью указана новая гидродинамическая интерпретация (р,ц) — аналитических функций Г.Н. Положего, направленная на решения задач фильтрации в однородных средах с прямолинейной анизотропией.

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

3. Проведено подробное исследование простейшей плоскорадиальной модели фильтрационного потока. С помощью плоскорадиальной модели:

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

- предложен практический способ аппроксимации функции проницаемости к(г) радиалъно-неоднородной призабойной зоны скважины.

4. Предложена элементарная теория проектирования оптимального забоя скважины с контейнерным гравийным фильтром (КГФ). В частности:

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

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

4.3 Предлолсен алгоритм проектирования забоя нефтедобывающей скважины с КГФ, обеспечивающий минимальный риск кольматации призабойной зоны и запескования скважины.

5. Выведена формула, дающая интегральное представление для решений уравнения фильтрации жидкости в изотропных средах с коэффициентом проницаемости к{х,у), принадлежащем серии решений уравнения A{^Щx^)-g(x)■^k(x^y) = 0 с произвольной наперёд задаваемой ограниченной непрерывной функцией ^(д;). (По аналогии с представленным в диссертации методом можно вывести формулу и для другого случая, когда g = g{y) и л\!1к{х,у))-2{у)-^к{х,у) = 0).

6. Предложена модификация известного в вычислительной электрофизике метода [52], который в диссертации назван методом распределенных источников. Метод распределенных источников предназначен для приближенного решения основной задачи теории фильтрации — задачи расчета дебитов одиночных и групповых скважин и доведён до практического алгоритма, требующего решения СЛАУ с выведенными для неё коэффициентами.

7. Выведена формула (3.4.9) для расчёта дебита отдельной скважины в круговой батарее из п скважин, учитывающая, что в призабойной зоне каждой скважины проницаемость кх не совпадает с проницаемостью к0 всего пласта. (В частном случае кх = к0 из выведенной формулы вытекает известная формула В. Н. Щелкачёва [129]). С помощью выведенной формулы проведено исследование относительного прироста суммарного дебита скважин в кольцевой батарее, если в призабойной зоне каждой скважины сделана гравийная обсыпка. Результаты исследования в диссертации представлены в табличном виде.

Список литературы диссертационного исследования кандидат технических наук Харченко, Юрий Владимирович, 2008 год

1. АКСЕНОВА, Н. А. Исследование и разработка техники, технологии заканчивания скважин с неустойчивыми коллекторами Текст. / Н. А. АКСЕНОВА // Автореферат дисс. на соиск. ученой степ. канд. техн. наук. - Тюмень, 2004. - 23 с.

2. АЛИЕВ, 3. С. Руководство по проектированию разработки газовых и газонефтяных месторождений Текст. / 3. С. АЛИЕВ, В. В. БОНДАРЕНКО // Республика Коми. г. Печора, издательство «Печорское время». 894 с.

3. АМИЯН, А. В. Промывка песчаной пробки пеной / А. В. АМИЯН, Н. П. ВАСИЛЬЕВА// М:, ВНИИОЭНГ, 1973. (Научно-технический обзор. Серия «Добыча»). — С. 97.

4. АМШЩ, В. А. Вскрытие и освоение нефтегазовых пластов Текст. /А. В. АМИЯН, Н. П. ВАСИЛЬЕВА // М.:, «Недра», 1980. - 380 с.

5. АРАВИН, В.И. Расчет плоской фильтрации в грунтах с криволинейной анизотропией Текст. / В.И. АРАВИН //Известия ВНИИГ им. Б.Е. Веденеева.- 1974. Т. 104.-С.З-9.

6. АШРАФЬЯН, М. О. Совершенствование конструкций забоев скважин Текст. / М. О. АШРАФЬЯН О. А. ЛЕБЕДЕВ, Н. М. САРКИСОВ // М, «Недра», 1987.-156 с.

7. БАРЕНБЛАТТ, Г.И. и др. Движение жидкостей и газов в природных пластах Текст. / Г.И. БАРЕНБЛАТТ, В.М. ЕНОТОВ, В.М. РЫЖИК // М.: Недра, 1984. 208 с.

8. БАСАРЫГИН, Ю.М. Закачивание скважин Текст. / Ю.М. БАСРЫГИН, А.И. БУЛАТОВ, Ю.М. ПРОСЕЛКОВ. М.: Недра, - 2000. - 670 с.

9. БАСНИЕВ, К. С. Нефтегазовая гидромеханика Текст. / К. С. БАСНИЕВ, Н. М. ДМИТРИЕВ, Г. Д. РОЗЕНБЕРГ Москва - Ижевск: Институт компьютерных исследований, — 2003. — 480 с.

10. БАШКАТОВ, А. Д. Сооружение гравийных фильтров за рубежом Текст.119

11. А. Д. БАШКАТОВ, М. И. ФАЗУЛЛИН, Е. Н. ДРЯГАЛИН // Техника и технология геолого-разведочных работ; организация производства. Обзор/ ВНИИ экономики минерального сырья и геологоразведочных работ (ВИЭМС) М.: - 1985. - 61 с.

12. БАТТТКАТОВ, А. Д. Предупреждение пескования скважин Текст. / А. Д. БАШКАТОВ // М.:, «Недра» 1991. - С. 176.

13. БАШКАТОВ, А. Д. Технология оборудования скважин многосекционными гравийными фильтрами Текст. / А. Д. БАШКАТОВ, Б. Л. СУШАНСКИЙ, В. И. ТАРАБУКИН //Спец. строит, работы. М. -1988.-Вып. 6.

14. БЕРГМАН, С. Интегральные операторы в теории линейных уравнений с частными производными Текст. / С. БЕРГМАН // М.:, «Мир». - 1964.

15. БОРИСОВ, Ю. П. Особенности проектирования разработки нефтяных месторождений с учетом их неоднородности Текст. / Ю. П. БОРИСОВ, 3. К. РЯБИНИНА, В. В. ВОИНОВ М.:, «Недра», - 1976. - 286 с.

16. БУЙКИС, A.A. Моделирование процессов фильтрации в слоистых средах методом консервативного осреднения / A.A. БУЙКИС // Диссертация д-ра физ.-мат. наук. Рига, 1987.-358 с.

17. ВЕКУА, И.Н. Обобщенные аналитические функции Текст. / И.Н. ВЕКУА// М.: «Физматиз», 1959 628 с.

18. ГАВРИЛКО, В.М. Фильтры буровых скважин Текст. / В.М. ГАВРИЛКО, B.C. АЛЕКСЕЕВ// М.: Недра, - 1985. - 334 с.

19. ГАРНАЕВ, А. Ю. MS Excel 2002: разработка приложений Текст. / А. Ю. ГАРНАЕВ// СПб.: БХВ-Петербург, 2004. - 768 с.

20. ГЛАДЫШЕВ, Ю.А. О методе перехода при решении задач фильтрации в пластах с переменными по простиранию мощностью и проницаемостью. Текст. / Ю.А. ГЛАДЫШЕВ // Гидромеханика. М.: МОПИ им. Н.К.Крупской, - 1974. - вып. 3. - С.217-221.

21. ГЛАДЫШЕВ, Ю.А. Построение потенциальных стационарных течений идеальной жидкости в искривлённом слое переменной толщины методом перехода Текст. / Ю.А. ГЛАДЫШЕВ // Труды МОПИ им. Н.К. Крупской. 1964. - Т. 142., вып. 5. - С.39-48.

22. ГОЛУБЕВ, Г.В. О некоторых точных решениях задачи об определении поля давлений в неоднородном нефтяном пласте Текст. / Г.В. ГОЛУБЕВ // Известия вузов. Нефть и газ. 1966. — №2. — С. 86-87.

23. ГОЛУБЕВ, Г.В. Об одном методе определения поля давлений в неоднородной пористой среде Текст. / Г.В. ГОЛУБЕВ // Известия АН СССР.МЖГ.- 1967.-№3.- С. 180-182.

24. ГОЛУБЕВ, Г.В. Определение поля давлений в кусочно-однородных пластах различной формы Текст. / Г.В. ГОЛУБЕВ // Труды унимверситета / Казанский ун-т. 1958. — Т.118. -Кн.2. - С.166-192.

25. ГОЛУБЕВ, Г.В. Фильтрация несжимаемой жидкости в неоднородной пористой среде Текст. / Г.В. ГОЛУБЕВ, Г.Г. ТУМАШЕВ // Казань: Изд-во Казанского ун-та, 1972. — 195 с.

26. ГОЛУБЕВА, О.В. Курс механики сплошных сред Текст. / О.В. ГОЛУБЕВА // М.: Высшая школа 1972. - С.364

27. ГУКАСОВ, Н. А. Теория и практика добычи газожидкостных смесей Текст. / Н. А. ГУКАСОВ, Г. Г. КУЧЕРОВ // М.: ООО «ИРЦ Газпром», -2005.-306 с.

28. ДОМБРОВСКИЙ, Г.А. Метод аппроксимации адиабаты в теории плоских121течений газа Текст. / Г.А. ДОМБРОВСКИЙ // М.: «Наука» 1964.

29. ДУБЕНКО, Д. В. Повышение эффективности работы скважин ПХГ путем совершенствования технологии сооружения гравийно-намывных фильтров / Д. В. ДУБЕНКО // Автореферат дисс. канд. техн. наук. -Краснодар, 2003. 23 с.

30. ЕНТОВ, В.М. Стационарная фильтрация однородной жидкости в элементе разработки нефтяного пласта с трещиной гидроразрыва Текст. / В.М. ЕНТОВ, В.В. МУРЗЕНКО// Изв. РАН. МЖГ. 1994.- №1. - С.104-112.

31. ЕРМИЛОВ, О. М. Геотехнологические основы разработки газовых месторождений крайнего севера Текст. / О. М. ЕРМИЛОВ // Автореферат диссертации на соискание ученой степени доктора технических наук. Уфа, 1991. 48 с.

32. ЕРМИЛОВ, О. М. Эксплуатация газовых скважин Текст. / О. М. ЕРМИЛОВ, 3. С. АЛИЕВ, В.В.РЕМИЗОВ, Л. С. ЧУГУНОВ // М.:, «Наука», 1995.-359 с.

33. ЕРМИЛОВ, О. М. Разработка крупных газовых месторождений в неоднородных коллекторах Текст. / О. М. ЕРМИЛОВ, В. Н. МАСЛОВ, Е. М. НАНИВСКИЙ// М.:, «Недра», 1987. - 207 с.

34. ЖЕЛТОВ, Ю. П. Разработка нефтяных месторождений Текст. / Ю. П. ЖЕЛТОВ // М.: ОАО «Издательство «Недра», 1998. 365 с.

35. ЗАЗОВСКИЙ, А.Ф. О стационарном притоке жидкости к скважине с вертикальной трещиной гидроразрыва большой протяженности Текст. /

36. A.Ф. ЗАЗОВСКИЙ, Г.Т. ТОДУА // Изв. АН СССР. МЖГ. 1990. - № 4. - С. 107-116.

37. ИЛЬИН, В.А. Численные методы решения задач электрофизики Текст. /

38. B.А. ИЛЬИН// М.: «Наука», 1985-334с.

39. ИЛЬИН, В.А. Основы математического анализа Текст. / В.А. ИЛЬИН, Э.Г. ПОЗНЯК// М.: Наука-Физматлит, 2000. - 4.1,2. - 616-447 с.

40. КАДЕТ, В.В. Фильтрация флюида в среде, содержащей эллиптическую трещину гидроразрыва Текст. /В.В. КАДЕТ, В.И. СЕЛЯКОВ// Известия вузов. Нефть и газ. 1988. - №5.-С.54-60.

41. КАЗАРЯН, В. А. Разработка технологии повышения продуктивности скважин с использованием реагентных методов разглинизации Текст. / В. А. КАЗАРЯН // Автореферат дисс. на соиск. учёной степени канд. техн. наук. М.: ВНИИГАЗ, 1995. 19 с.

42. КОЛЛИНЗ, Р. Течение жидкостей через пористые материалы Текст. / Р. КОЛЛИНЗ //М.: Мир, 1964,-350 с.

43. КОРОТЕНКО, В.А. Повышение добывных возможностей нефтяных скважин, эксплуатирующих слоистый коллектор Текст. / В.А. КОРОТЕНКО // Автореферат диссертации на соискание учёной степени кандидат технических наук — ТюмГНГУ, Тюмень, 2007 24 с.

44. КОЧИНА, И.Н. Фильтрация через глинистые корки Текст. / И.Н. КОЧЕНА, H.H. МИХАЙЛОВ // Изв. АН СССР. МЖГ. 1976. - №6. -С.70-75.

45. ЛАВРЕНТЬЕВ, М.А. Методы теории функций комплексного переменного Текст. / М.А. ЛАВРЕНТЬЕВ, Б.В. ШАБАТ// М.: Наука, -1973.-736 с.

46. ЛЕЙБЕНЗОН, Л.С. Движение природных жидкостей и газов в пористой среде Текст. / Л.С. ЛЕЙБЕНЗОН // М.-Л., 1947. - 244 с.

47. ЛУКОМСКАЯ, М.А. Решение некоторых задач о притоке жидкости к скважинам Текст. / М.А. ЛУКОМСКАЯ // ПММ. 1947. - Т. XI. - вып. 6.

48. МАСКЕТ, М. Течение однородных жидкостей в пористых средах Текст. /М. МАСКЕТ// М.-Л.: ГНТИ, 1949.-628 с.

49. МИЗЕС, Р. Математическая теория течений сжимаемой жидкости Текст. / Р. МИЗЕС //М.:, Иностранная литература, 1961 г.

50. МИРЗАДЖАНЗАДЕ, А. X. Гидродинамика в бурении Текст. / А. X. МИРЗАДЖАНЗАДЕ, В. М. ЕНТОВ // М.:, «Недра», 1985. 196 с.

51. МИХАЙЛОВ, H.H. Изменение физических свойств горных пород в околоскважинных зонах Текст. / H.H. МИХАЙЛОВ// М.: Недра, 1987. - 152 с.

52. НАЗАРОВ, Г. И. О двух классах точных общих решений уравнений типа Чаплыгина в газовой динамике Текст. / Г. И. НАЗАРОВ //Ученые записки № 73 «Проблемы гидромеханики и теории упругости». Вып. 2. Изд-во Томского университета. Томск, 1973.

53. НАЗАРОВ, Г. И. Точное решение уравнений газовой динамики Текст. / Г. И. НАЗАРОВ ////Известия АН СССР. МЖГ, № 3, 1968.

54. ПИВЕНЬ, В.Ф. К теории осесимметричных обобщенных аналитических функций в динамических процессах Текст. / В.Ф. ПИВЕНЬ // Докл. АН СССР. 1990. -Т.313. -№6.-С. 1424-1426.

55. ПИВЕНЬ, В.Ф. О теории двумерных процессов в слоях переменной проводимости, характеризуемых степенью гармонической функции Текст. / В.Ф. ПИВЕНЬ // ДАН. 1995. - Т.344. - №5. - С. 327-629.

56. ПИЛАТОВСКИИ, В.П. Основы гидромеханики тонкого пласта Текст. / В.П. ПИЛАТОВСКИЙ// М.: Недра, 1966.-317 с.

57. ПИРВЕРДЯН, A.M. Нефтяная подземная гидравлика Текст. / A.M. ПИРВЕРДЯН// Баку: Азнефтеиздат, 1956.-332 с.

58. ПОЛОЖИЙ, Г.Н. Обобщение теории аналитических функций комплексного переменного Текст. / Г.Н. ПОЛОЖИИ // Киев: Изд-во Киевск. ун-та, 1965. 442 с.

59. ПОЛОЖИЙ, Г.Н. Теория и применение р — аналитических и (р, q)— аналитических функций Текст. / Г.Н. ПОЛОЖИЙ // Киев, «Наукова думка», 1973 423 с.

60. ПОЛУБАРИНОВА-КОЧИНА, П.Я. Теория движения грунтовых вод Текст. / П.Я. ПОЛУБАРИНОВА-КОЧИНА // М.: Недра, 1977.-664 с.

61. ПРУСОВ, И.А. Двумерные краевые задачи фильтрации Текст. / И.А. ПРУСОВ // Минск: Изд. университетское, 1987. - 182 с.

62. ПРУСОВ, И.А. Вывод основных уравнений фильтрации жидкости в анизотропной среде Текст. / И.А. ПРУСОВ, И.А. ВЕРЕМЧУК // Изв. АН БССР. Серия физ.-мат. наук. Минск, 1974.-№1.-С.109-112.

63. ПЫХАЧЕВ, Г.Б. Подземная гидравлика Текст. / Г.Б. ПЫХАЧЕВ, ИСАЕВ Р.Г. // М.: Недра, 1973.-360с.

64. РАБИНОВИЧ, Н.Р Инженерные задачи механики сплошной среды в бурении Текст. / Н.Р. РАБИНОВИЧ II- М.: Недра, 1989. 270 с.

65. РАДЫГИН, В.М. Фильтрационная теорема о двух окружностях Текст. /125

66. В.M. РАДЫГИН // Задачи гидродинамики при усложненных моделях среды МОИП. М.: Наука, 1985. -С. 18-23.

67. РАДЫГИН, В.М. Применение функций комплексного переменного в задачах физики и техники Текст. / В.М. РАДЫГИН, ГОЛУБЕВА О.В. //М.: Высшая школа, 1983.-160с.

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

69. РОММ, Е.С. Фильтрационные свойства трещиноватых горных пород Текст. / Е.С. РОММ //- М.: Недра, 1966 238 С.

70. САМАРСКИЙ, A.A. Введение в численные методы Текст. / A.A. САМАРСКИЙ //- М.: «Наука», 1982 271с.

71. САМАРСКИЙ, A.A. Численные методы Текст. / A.A. САМАРСКИЙ, A.B. ГУЛИН// М.: «Наука», 1989 430с.

72. СВЕШНИКОВ, А.Г. Теория функций комплексной переменной Текст. / А.Г. СВЕШНИКОВ, А.Н. ТИХОНОВ// М.: Наука, 1970. - 304 с.

73. СМИРНОВ, С. В. Технология повышения эффективности работы скважин с помощью устройства гидропескоструйной щелевой винтовой перфорации Текст. / C.B. СМИРНОВ, И. РЖАНОВ, А. ГИЛЬДЕЕВ //НефтьГаз промышленность. -2007.-№3 (31), май. С. 54-57.

74. СОЛОМКО, Ю.Л. О моделировании работы скважины в анизотропном грунте Текст. / Ю.Л. СОЛОМКО // МОИП. Избранные задачи гидродинамики. М.: Наука. - 1977. - С.92-96.

75. Справочник по специальным функциям Текст. / Под ред. Абрамовича М. и Стигана И. М.: Наука, 1979. - 830 с.

76. СТО Газпром РД 2.1 154 - 2005. Ввод в действие с 01.03.2005. Утв. зам. пред. правления ОАО «Газпром» Ананенковым А. Г. / Фильтры гравийные газовых скважин. Оборудование, материалы, технология сооружения. -М.: ООО «ИРЦ Газпром», 2005. - 13 с.

77. ТИМАШЕВ, Г.В. Скважинные фильтры (по патентным материалам зарубежных стран) Текст. / Г.В. ТИМАШЕВ, Т.А. АТАКУЛОВ, О.Ж.126

78. КАЛНИН, A.A. ГОРОШКО// ВНИИЭГАЗПРОМ. НТО. Серия "Разработка и эксплуатация газовых и газоконденсатных месторождений". -М., 1977. вып. 13. - С.45.

79. ТИХОНОВ, А.Н. Уравнения математической физики Текст. / А.Н. ТИХОНОВ, A.A. САМАРСКИЙ//М.: Наука, 1972.-735 с.

80. ТОЛПАЕВ, В.А. О построении точных решений задач напорной фильтрации в некоторой серии анизотропных сред Текст. / В.А. ТОЛПАЕВ // Известия СКНЦ ВШ. Естественные науки. 1979.-№4.-С. 33-36.

81. ТОЛПАЕВ, В.А. Решение задач фильтрации в кусочно-неоднородных средах методом моделирования границ раздела эквипотенциалями течения Текст. / В.А. ТОЛПАЕВ // Известия вузов. Северо-Кавказский регион. Естественные науки. 1999. - № 4. -С.39-43.

82. ТОЛПАЕВ, В.А. Связь плоскопараллельных течений в изотропных и анизотропных грунтах Текст. / В.А. ТОЛПАЕВ // Гидромеханика. М.: МОПИ им. Н.К. Крупской, 1975.-вып. 4.-С. 11-16.

83. ТОЛПАЕВ, В.А. Уравнения линейной напорной плоскопараллельной фильтрации в анизотропных средах Текст. / В.А. ТОЛПАЕВ // Известия СКНЦВШ. Естественные науки. 1984.-№2.-С.45-49.

84. ТОЛПАЕВ, В.А. Фильтрация жидкости в анизотропных и неоднородных грунтах Текст. / В.А. ТОЛПАЕВ // Учебное пособие. Ставрополь: СевКавГТУ, 2000. - 196 с.

85. ТОЛПАЕВ, В.А. Численно-аналитические методы расчета дебитов одиночных и групповых скважин в неоднородных средах Текст. / В.А. ТОЛПАЕВ // Известия вузов. Северо-Кавказский регион. Естественные науки. 2000. - №1. - С.53-57.

86. ТОЛПАЕВ, В.А. Решение краевой задачи Дюпюи для среды с прямолинейной анизотропией Текст. / В.А. ТОЛПАЕВ, А. Д. ЖЕРНОВОЙ// Известия СКНЦ ВШ. Естественные науки. 1988. -№4.-С.80-87.

87. ТОЛПАЕВ, В.А. Гидродинамические особенности течения жидкости в призабойной зоне скважины Текст. / В.А. ТОЛПАЕВ, В.В. ЗАХАРОВ // Вестник СевКавГТУ. Серия «Физико-химическая». — Ставрополь, 2003. -№1(7).-С. 120-127.

88. ХАРЧЕНКО, Ю.В. Построение точных решений задач фильтрации жидкости к скважине методом характеристических функций /В.А. ТОЛПАЕВ, Ю.В. ХАРЧЕНКО // Сборник научных трудов СевКавГТУ.

89. Серия «Естественнонаучная». Ставрополь, 2002.-вып.5.-С.98-101.

90. ХАРЧЕНКО, Ю.В. Потенциалы фильтрационных течений от источников и стоков в областях круговой формы, полуплоскости, квадранта /В.А. ТОЛПАЕВ, Ю.В. ХАРЧЕНКО // Вестник СевКавГТУ. Серия «Физико-химическая» Ставрополь, 2003. - №1(7). - С. 150-158

91. ХАРЧЕНКО, Ю.В. Вычислительные эксперименты по оптимизации размещения нефтедобывающих скважин в круговом пласте /В.А. ТОЛПАЕВ, B.C. КИРИЛЛОВ, Ю.В. ХАРЧЕНКО //Сборник научных трудов СевКавГТУ. Серия «Нефть и газ». Ставрополь, 2000. - вып. 3. -С. 131-136.

92. ХАРЧЕНКО, Ю.В. Математическая модель циклического взаимодействия скважин с индивидуальными фильтрационными свойствами призабойных зон /В.А. ТОЛПАЕВ, B.C. КИРИЛЛОВ, Ю.В. ХАРЧЕНКО // Периодическое издание «Циклы». Ставрополь,129

93. СевКавГТУ, 2000. вып. 2. - С.57-62.

94. ХАРЧЕЬЖО, Ю. В. Значение простейшей модели плоско-радиального фильтрационного потока для практики разработки пластов. /В.А. ТОЛПАЕВ, Ю.В. ХАРЧЕНКО // Научно-технический журнал «Нефтепромысловое дело»- 2005 №10. - С. 9 - 13.

95. ХАРЧЕЬЖО, Ю. В. Математическое моделирование течений к одиночным и групповым скважинам. /В.А. ТОЛПАЕВ, Ю.В. ХАРЧЕНКО // Научно-технический журнал «Автоматизация, телемеханизация и связь1 в нефтяной промышленности», М.:, —2004 — №11.

96. ХАРЧЕЬЖО, Ю. В. Числовые оценки зависимости дебитов скважин кольцевой батареи от скачка проницаемости в призабойных зонах /В.А. ТОЛПАЕВ, Ю.В. ХАРЧЕЬЖО // Ежемесячный научно-технический журнал «Нефтепромысловое дело», Москва, 2005. - № 11.— С. 42-44

97. ХАРЧЕЬЖО, Ю. В. Влияние проницаемости гравийного фильтра на дебит буровой скважины при линейном законе Дарси Текст. / В.А. ТОЛПАЕВ, Ю.В. ХАРЧЕЬЖО, В.В. ЗАХАРОВ // Известия вузов. Северо-Кавказский регион. Естественные науки, — 2003. — № 3 — С. 36-41.

98. ХАРЧЕЬЖО, Ю. В. Закон Дарси для высокопроницаемых сред /В. А. ТОЛПАЕВ, А. Е. ЬДИЦКУН, Ю.В. ХАРЧЕЬЖО // НТЖ «Обозрение прикладной и промышленной математики» М.:, — 2004. — Том 11, вып. 3. -С. 677-680

99. ШАХНАЗАРОВ, А. А. Крепление призабойной зоны скважин Текст. / А. А. ШАХНАЗАРОВ // М., Гостоптехиздат, 1959. 83 с.

100. ШЕЙДЕГГЕР, А.Э. Физика течения жидкостей через пористые среды Текст. / А.Э. ШЕЙДЕГГЕР //М.: Гостоптехиздат, 1960. 250 с.

101. ЩЕЖАЧЕВ, В.Н Подземная гидравлика Текст. / В.Н. ЩЕЛКАЧЕВ, Б.Б. ЛАПУК// Москва-Ижевск: НИЦ «Регулярная и хаотическая динамика», -2001- с.736

102. ЯНКЕ, Е. Специальные функции Текст. / Е. ЯНКЕ, Ф.ЭМДЕ, Ф.ЛЁШ// -М.: Наука, 1977.-342с.

103. ЯРМИЦКИЙ, А.Г. Фильтрационная теорема о двух окружностях Текст. / А.Г. ЯРМИЦКИЙ // Изв. АН СССР. МЖГ. 1986. - № 4. - С.76-82.

104. MARCUS H. The permeability of sample of an anisotropy medium // J. Geophys. Res. 1962. V. 67. №13, p. 5215-5225.

105. MEEGODAN.J., KINGI.P., ARULANDAN K. An expression for the permeability of anisotropy granular media // Int. J. number, anal, methods in geomechanics. 1989. V. 13. p. 575-598.

106. MUSKAT M. Physical principles of oil production. New York. McGraw-Hill. 1949.922 р.

107. MUSKAT M. The flow of homogenous fluids through porous media. Ann. Arbour. Mich. Edwards. 1946. 736 p.

108. NIKOLAEVSKIJ V.N. Mechanics of Porous and Fractured Media. Singapore: World Scientific. 1990.472р.

109. PIVEN' V.R The theory of two-dimensional processes in inhomogeneous layers with power law of their conductivity variation // J. Appl. Maths. Mechs. 1997, Vol. 61, №4, P. 577-586.

110. SCHEIDEGGER A.E. The physics of flow through porous media. Univ. of Toronto Press. 1974, 3d edition. 353 p.