Бесплатный автореферат и диссертация по наукам о земле на тему
Комплексное решение проблемы оценки коэффициента вытеснения нефти в различных геолого-технологических условиях
ВАК РФ 25.00.12, Геология, поиски и разведка горючих ископаемых

Автореферат диссертации по теме "Комплексное решение проблемы оценки коэффициента вытеснения нефти в различных геолого-технологических условиях"

На правах рукописи

Хижняк Григорий Петрович

КОМПЛЕКСНОЕ РЕШЕНИЕ ПРОБЛЕМЫ ОЦЕНКИ КОЭФФИЦИЕНТА ВЫТЕСНЕНИЯ НЕФТИ В РАЗЛИЧНЫХ ГЕОЛОГО-ТЕХНОЛОГИЧЕСКИХ УСЛОВИЯХ

25.00.12 - Геология, поиски и разведка нефтяных и газовых месторождений

Автореферат диссертации на соискание ученой степени доктора технических наук

1 2 [Г.? ¿0:2

Пермь 2012

005012332

Работа выполнена в Федеральном государственном бюджетном учреждении высшего профессионального образования «Пермский национальный исследовательский политехнический университет» и филиале ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть» в г. Перми

Научный консультант:

доктор геолого-минералогических наук Галкин Сергей Владиславович

Официальные оппоненты:

доктор технических наук, профессор Силаев Валерий Аркадьевич

доктор технических наук, профессор Крысин Николай Иванович

доктор физико-математических наук, профессор Долгаль Александр

Сергеевич

Ведущее предприятие - Научно-исследовательское, проектное и производственное предприятие «Недра» (г. Пермь).

Защита состоится 29 марта 2012 года в 15 часов на заседании совета по защите докторских и кандидатских диссертаций Д 212.188.03 при Пермском национальном исследовательском политехническом университете по адресу: 614990, г. Пермь, Комсомольский проспект, 29, ауд. 423 б.

С диссертацией можно ознакомиться в научной библиотеке ПНИПУ.

Автореферат разослан «_» февраля 2012 г.

Ученый секретарь

совета, доктор геолого-минералогических наук, профессор

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Актуальность работы. В настоящее время в соответствии с требованиями ОСТ 39-195-86 основным методом определения коэффициента вытеснения (Кет) является метод лабораторного моделирования. В качестве моделей пласта используются составные образцы керна изучаемого объекта, насыщенные нефтью и содержащие остаточную воду. Одним из важнейших параметров, влияющих на коэффициент вытеснения, является проницаемость коллектора. Однако при ее обосновании встречается ряд проблем.

В настоящее время проницаемость определяется по керну, по геофизическим (ГИС) и гидродинамическим (ГДИ) исследованиям скважин. Наиболее надежно отражающей пластовые условия считается проницаемость, определенная по гидродинамическим исследованиям. В случае достаточного количества исследований фазовая проницаемость по нефти при ГДИ берется за основу. Однако процент исследованных скважин методом ГДИ в среднем составляет порядка 30 % от общего числа скважин, часто не превышая 10 %.

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

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

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

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

Цели исследований заключаются:

-в использовании материалов ГИС для получения достоверных данных о значениях коэффициентов проницаемости и вытеснения;

- в научном обосновании количественных характеристик микроструктуры порового пространства пород-коллекторов с целью оценки коэффициента вытеснения без непосредственного физического моделирования процесса нефтевытеснения на лабораторных установках;

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

Основные задачи исследований определяются поставленными целями и формулируются следующим образом:

-разработка моделей прогноза проницаемости пород-коллекторов нефтяных месторождений Пермского Прикамья по данным ГИС;

- научное обоснование использования коэффициента подвижности для оценки коэффициента вытеснения в различных геолого-технологических условиях;

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

Научная новизна

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

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

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

Защищаемые положения:

-метод оценки коэффициента проницаемости пород-коллекторов по данным ГИС с помощью многомерных моделей;

-способ определения застойных зон порового пространства пород-коллекторов методом «метки»;

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

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

Практическая ценность работы заключается:

- в повышении точности оценок извлекаемых запасов нефти за счет использования результатов геолого-петрофизической интерпретации данных ГИС для моделирования пространственного изменения коэффициента вытеснения;

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

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

Реализация результатов исследования. Исследования по теме диссертации выполнялись с 1974 по 2011 г. Результаты исследований реализованы в двух авторских свидетельствах [2, 13], патентах на изобретение [22] и полезную модель [37], ряде методических рекомендаций и нашли применение практически во всех отчетах по подсчету запасов, проектированию и анализу разработки нефтяных месторождений Пермского Прикамья в 1980-2011 гг.

Исходные данные. Для решения поставленных задач выполнены лабораторные исследования более 2500 образцов керна пород-коллекторов, отобранных более чем из 400 нефтяных залежей северо-востока Волго-Уральской нефтегазоносной провинции, с привлечением данных по месторождениям Западной Сибири, Свердловской, Оренбургской областей и Красноярского края методами физического моделирования процесса нефте-вытеснения, «метки», фильтрации взаиморастворимых жидкостей.

Апробация работы и публикации. Материалы диссертации докладывались и обсуждались на всесоюзном семинаре «Современные проблемы и математические методы теории фильтрации» /М., 1984 г./; на заседании Межведомственного совета по лабораторным методам, применяющимся в нефтегазовой геологии (группа по коллекторам) /Львов, 1984 г./; на всесоюзном семинаре «Современные проблемы и математические методы теории фильтрации» ПЛ., 1984 гУ; на межсекционном семинаре «Математическое

описание микронеоднородных сред и расчет их физических свойств» /М., 1985 г./; на научно-технической конференции «Совершенствование методов поисков, разведки и разработки нефтяных и газовых месторождений» /Пермь, 1989 г./; на конференции «Основные направления научно-технического прогресса в развитии нефтяной промышленности Пермского Прикамья» /Пермь, 1989 г./; на краевой научно-технической конференции студентов, аспирантов и молодых ученых «Проблемы комплексного освоения месторождений полезных ископаемых в Пермском крае» /Пермь, 2007 г./; на международной научно-технической конференция «Нефтегазовое и горное дело» /Пермь, 2009, 2010, 2011 гг./; на научно-технической конференции «Актуальные проблемы геологических исследований и разработки месторождений Пермского края» /Пермь, 2010 г./.

Результаты работ использованы в 5 внедренных в ООО «ЛУКОЙЛ-ПЕРМЬ» методических руководствах [4,9,10,17,18].

Основные положения диссертации опубликованы более чем в 50 научных работах, изданных в Москве, Минске, Уфе, Перми. За выполненные научные разработки автору присужден диплом Пермского обкома комсомола, Обкома профсоюза рабочих нефтяной и газовой промышленности.

Автор приносит глубокую благодарность научному консультанту доктору геолого-минералогических наук, профессору С.В. Галкину.

Большой вклад в постановку и решение задачи обоснования коэффициента вытеснения нефти внесли доктор технических наук, профессор Б.И. Тульбович и доктор геолого-минералогических наук, профессор Ю.В. Шурубор|, которым автор выражает искреннюю благодарность.

Искреннюю признательность за консультации и поддержку в работе над диссертацией автор выражает H.A. Лядовой, A.B. Распопову, В.И. Галкину, A.B. Растегаеву, A.C. Некрасову, В.Б. Бейзману, М.Ф. Серкину, В.К. Червяковой. Научным исследованиям способствовала действенная помощь сотрудников Центра исследования керна и пластовых флюидов филиала ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть».

Объем и структура работы. Диссертация состоит из введения, шести глав и заключения, изложена на 259 страницах машинописи, включая 63 рисунка и 31 таблицу. Список литературы состоит из 218 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе «Анализ современного опыта оценки коэффициента вытеснения при проектировании разработки нефтяных месторождений» рассмотрена история определения коэффициента вытеснения нефти (Квт), ведущаяся в Пермском крае с 1970 г.

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

В.И. Азаматовым и Н.М. Свихнушиным (1976 г.) рассмотрены возможности оценки коэффициента вытеснения нефти водой по результатам исследований скважин микрозондами и микроэкранированными зондами на основе использования параметров промытой части пласта. В.М. Гусевым и Б.И. Тульбовичем (1984 г.) для терригенных коллекторов Пермского Прикамья, представленных песчаниками с глинистым цементом, предложена зависимость Квт от двойного разностного параметра Ыпг, определяемого по

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

В 1978 г. вышел первый отраслевой стандарт ОСТ 39-070-78 на определение Квт нефти водой в лабораторных условиях, в разработке которого приняли участие специалисты Всесоюзного нефтегазового научно-исследовательского института (ВНИИ), Башкирского научно-исследовательского проектного института (БашНИПИнефть), Сибирского научно-исследовательского института нефтяной промышленности (СибНИИНП) и Украинского государственного научно-исследовательского и проектного института нефтяной промышленности (УкргипроНИИнефть).

С 1986 г. и по настоящее время руководящим документом, регламентирующим порядок определения коэффициента вытеснения, является ОСТ 39-195-86, разработанный при участии одного из ведущих петрофизи-ков Пермского Прикамья, доктора технических наук, профессора Бориса Из-раилевича Тульбовича. В 1988-1989 гг. автором диссертации совместно с В.Г. Михневичем и Б.И. Тульбовичем [9,11] опубликованы результаты обработки данных по определению коэффициента вытеснения для 39 залежей в терригенных и 22 в карбонатных палеозойских отложениях Пермского

Прикамья. В расчетных формулах в качестве параметров-аргументов были использованы: коэффициент подвижности (Кпр/рн), представляющий отношение проницаемости пористой среды к вязкости нефти; и отражающий реологические свойства нефти и прочность межфазных пленок комплексный параметр (К</Ка) в виде отношения концентраций смол силикагелевых и асфальтенов в нефти. В 1994 г. теми же авторами было показано [17], что точность расчетных оценок коэффициента вытеснения возрастает, если использовать свои уравнения для нижнекаменноугольных терригенных, девонских терригенных, среднекаменноугольных карбонатных и турнейских карбонатных коллекторов. К 1996 г. экспериментальный материал по определению коэффициента вытеснения значительно увеличился, что позволило оценивать его расчетным способом более дифференцированно, с учетом возраста продуктивных отложений и их структурно-тектонической приуроченности [18]. Конкретные зависимости приведены в табл. 1 и на рис. 1.

Таблица 1

Расчетные формулы для определения коэффициента вытеснения нефти водой (терригенные и карбонатные коллекторы Пермского Прикамья)

Структурно-тектонические элементы Количество Расчетная формула Погрешность, %

месторождений опытов абсолютная отно-ситель ная

Каширо-ве рейские карбонатные пласты

ПП 7 10 Кет = 0,02161п (Кпо/ц,,) + 0,6884 4,1 7,1

Башкирские карбонатные пласты

БКВ 6 24 Квт = 0,01361п (Кщ,/ |хн) + 0,6231 2,8 4,9

БС 6 11 Квт = 0,04641п (Кщ, / цн ) + 0,7513 3,4 6,4

ПС б 27 КЕХ = 0,02421п (Кго / ц„ ) + 0,6554 4,0 7,4

ВКВ 4 12 Кет = 0,04801п(Кго/ц„) + 0,7309 2,2 4,5

сд 7 13 Квт = 0,04341п(Кго/дн) +0,7483 2,6 4,7

ПП 29 87 Квт = 0,02541п (Кщ, / цн ) + 0,6552 3,0 5,6

Тульский и бобриковский терригенн ые пласты

БКВ 11 29 К„ = 0,03001п (Кго / цн ) + 0,7227 4,0 6,3

Окончание табл. 1

Структурно-тектонические элементы Количество Расчетная формула Погрешность, %

месторождений опытов абсолютная отно-ситель ная

БС 12 30 Квт = 0,0320ln(Km/u„) + 0,7057 3,0 4,9

ПС 5 19 Квт = 0,0593 1п(Кпп/ц„) + 0,7873 3,7 6,2

ВКВ 8 21 Квт = 0,0308 1п(Кт/ц„) + 0,7193 4,2 7,1

сд 5 19 Квх = 0,02841п (К™ / ц„) + 0,7487 3,8 5,8

ПП 40 112 Квт = 0,03801п (Кго/ д„) + 0,7343 3,8 6,4

Малиновский терригенный пласт

БКВ б 13 Квт = 0,0425 ln (Кт / ц„) + 0,7536 2,5 3,8

ПП 11 16 Квт = 0,05871п(Кт/ц„) + 0,7837 3,1 4,8

Турнейский и фаменский карбонатные пласты

БКВ 4 9 Квх = 0,0187 1п (Кпо / ц„ ) + 0,7154 2,9 4,9

сд в целом 7 14 Квт = 0,06191п (К„р/ цн) + 0,8432 3,5 6,2

север 4 5 Квт = 0,0366 1п(Кт/ц„) + 0,7814 1,3 2,1

юг 3 9 КВТ = 0,0850 1п (Кго/ ц„) + 0,9146 2,5 4,5

ПП 20 47 Квт = 0,03021п (Кго / ц„ ) + 0,7420 4,0 6,7

Девонские терригенные пласты

ВКВ 3 8 Квт = 0,09561п (К™ / ц„ ) + 0,9530 3,1 4,9

ПП 5 11 Кет = 0,0891 1п (Кго/ ц„) + 0,9458 3,3 5,2

Примечание: БКВ - Бымско-Кунгурская впадина; БС - Башкирский свод; ПС — Пермский свод; ВКВ - Верхнекамская впадина; СД - Соликамская депрессия; ПП -Пермское Прикамье в целом. Проницаемость - в мкм2, вязкость - в мПа-с.

Коэффициент подвижности, мкм2/мПа-с

Коэффициент подвижности, мкм2/мПа-с

Коэффициент подвижности, мкм2/мПас

Рис. 1. Зависимости коэффициента вытеснения нефти водой от коэффициента подвижности по пластам для различных тектонических структур

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

Из табл. 1 видно, что в большинстве случаев уравнения регрессии, относящиеся к отдельным структурно-тектоническим элементам, дают оценки коэффициента вытеснения более точные, чем уравнения, относящиеся к Пермскому Прикамью в целом.

Во второй главе «Использование материалов ГИС для количественной оценки свойств пород-коллекторов при подсчете запасов» рассмотрены возможности применения многомерных зависимостей для оценки проницаемости.

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

Для определения проницаемости Кцр по данным ГИС чаще всего используют парные связи проницаемости с пористостью (Кп), параметром насыщения (Рн), глинистостью, параметром самопроизвольной поляризации (апс), относительным удельным сопротивлением (рп/рв). Б.Ю. Вендель-штейн, Л.В. Сиротенко и др. показали, что для повышения точности получаемых оценок проницаемости нужно пользоваться многомерными зависимостями. В идеале каждая из таких зависимостей должна относиться к конкретному пласту конкретного месторождения или даже к отдельному участку такого пласта.

Задача оценки проницаемости по данным ГИС решалась для терриген-ных отложений визейского яруса Пермского Прикамья [36] с применением уравнений множественной регрессии, учитывающих информацию о пористости Кп, относительной глинистости (г|гл), относительном удельном сопротивлении рп/рв, т.е. о тех параметрах, которые определяются не только при лабораторном исследовании керна, но и по общепринятым отраслевым методикам обработки материалов ГИС. Установлено, что уравнения приемлемой точности действительно можно найти, если с многомерными уравнениями регрессии работать под определенные пласты конкретных месторождений. Например, получены следующие зависимости:

1) для Уньвинского месторождения

Ьо Кпр = 3,0650 I* Кп -1,9514 Т1гл+ 0,5448(р„/рв) -4,8725; (2.1)

2) для Кокуйского месторождения

^ Кпр= 0,0061 Кп - 2,2390 т]гл + 0,7562 (рп/рв) -2,1399; (2.2)

3) для Ножовского месторождения

Ьё Кпр = 0,0684 Кп - 0,7901 т]гл + 0,7992 (рп/рв) - 1,2060. (2.3) Эти зависимости, характеризуются высокими значениями коэффициента детерминации (Л2), для К„р соответственно равными 0,87, 0,79 и 0,79. Наглядное представление о точности одной из них дает рис. 2.

1

3 10 000

7

о

£ 1000

5

1 100

1

л н 10

í

I 1

1 с 0.1

« I • •

* •

• • *нт корреляции —

г-0>» ........

0.1

1 10 100 1000 Проницаемость, измеренная по керну, 10Г1 мкм1

10 000

Рис. 2. Соотношение измеренной и расчетной проницаемостей терригенных отложений визейского яруса Уньвинского месторождения (проницаемость рассчитана по уравнению (2.1) с использованием значений Кп, Т1гл, рп/рв, определенных по керну

в лаборатории)

т • •- Уйминекое м-е * - Чаимимемм мл

* • Д-Юрчукеиоеи-в Коэффициент корреляции г >0,99

Проницаемость, измеренная по керну, 1СГ3 мкм1

Рис. 3. Соотношение измеренной и расчетной проницаемостей терригенных отложений визейского яруса нефтяных месторождений Соликамской депрессии (проницаемость рассчитана по уравнению (2.1) с использованием значений Кп, ЛгтРг/Рв.

оцененных по данным ГИС)

Если для рис. 2 расчетные значения проницаемости получены с использованием К„, т)гл> Рп/Рв) определенных по керну в лаборатории, то рис. 3 отображает ситуацию, когда КП1 Т1ГЛ) рп/рв оценивались по данным ГИС.

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

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

Поскольку коэффициент вытеснения зависит в основном от динамических свойств порового пространства, определяемых микроструктурой этого пространства, требуются специальные петрофизические исследования, позволяющие детально изучать микроструктуру породы-коллектора. Ю.А. Чизмаджев и Л.М. Марморштейн рекомендовали для изучения микроструктурных особенностей пористых сред в динамическом режиме использовать методы «метки» и фильтрации взаиморастворимых жидкостей (сме-симое вытеснение).

Вполне очевидно, что одним из факторов, оказывающих определяющее влияние на коэффициенты уравнений, которыми описывается зависимость коэффициента вытеснения от коэффициента подвижности, является наличие в поровом пространстве так называемых застойных зон. Для определения этого параметра используется метод «метки» [1-4, 7, 8, 15]. Толчком к освоению и развитию метода послужило ознакомление с моделью пористой среды в виде совокупности ячеек идеального перемешивания с застойными зонами, предложенной Ю.А. Чизмаджевым. Авторское свидетельство [2] на способ определения застойных зон в породах-коллекторах методом «метки» получено.

Пористая среда представляется набором одинаковых ячеек, каждая из которых имеет объем V и протяженность /. Доля этого объема, равная а„, занята застойными зонами, а остальная часть (1 - а„) - проточная среда, в которой происходит идеальное перемешивание. Объемная скорость переноса вещества в проточной среде равна дм, объемная скорость обмена между этой средой и застойными зонами - ры, линейная скорость потока - им. Механизм обмена между проточной и застойной зонами может быть как конвективным, так и диффузионным. При выполнении исследований породы-коллектора методом «метки» в образец вводят метку - небольшое количество динамически нейтральной примеси, которая не меняет свойств потока, и далее следят (в случае газовой метки с помощью хроматографа) за измене-

нием рассматриваемой в качестве функции времени / концентрации ф) примеси (размыванием «метки») на выходе из пористой среды. Форма выходной кривой связана со структурными параметрами образца. Согласно Ю.А. Чизмаджеву распределение концентрации метки на выходе образца, длина которого равна Ь, описывается уравнениями:

ф) = а ехр[—¿>(/ - с)2] при / < с, (3.1)

ф) = а ехр[-г»(/ - с)2] + (1 ехр(-уО при г > с, (3.2)

где

\-k2x-

1

ъ=

1

И..

Рм_

и«

(3.3)

у^ Р» "м <7« «,/

На основе совместного решения уравнений (3.3) по значениям коэффициентов а, ¿>, с, у, получаемым при аналитической аппроксимации экспериментальной кривой размывания «метки», можно найти оценки доли застойных зон ац, коэффициента обмена скорости обмена р„, длины ячейки идеального перемешивания / [2, 3, 4]. Получены следующие формулы:

1

а., =-

Ь.

<7м

, а 2 л 1

2Ьс:

Ч

Л'

(3.4)

(3.5)

Ры ='

2 Ьс2с1

1 а №

(3.6)

/=-

2 Ьс

(3.7)

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

Для пород-коллекторов обычно наблюдается достаточно тесная корреляционная связь между Кет и КпрГ, которая в данном случае имеет вид:

Квт = 0,1370 1§КпрГ +0,3107, Л2 = 0,75.

Однако анализ этой связи показал, что данное корреляционное поле состоит из двух подполей. Для терригенных пород значения Квт в зависимости от проницаемости обычно выше, чем для карбонатных (рис. 4), что хорошо подтверждается следующими уравнениями регрессии:

- для терригенных пород

К„т = 0,1211 \ё КпрГ + 0,3822, Я2 = 0,77;

- для карбонатных

Квт = 0,1075 ^КпрГ + 0,3008, Л2 = 0,76.

Рис. 4. Зависимость коэффициента вытеснения от проницаемости

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

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

Квт =-1,1906 ам +0,9140, г = -0,82.

При сохранении общей тенденции к уменьшению Кет при увеличении ам наблюдается дифференциация их соотношений для терригенных и карбонатных пород:'

- для терригенных пород: Квт =-1,4931 а„+ 1,0201, г = -0,83,

-для карбонатных пород: Кст = -0,7410 ам + 0,7135, г = -0,84.

Угловые коэффициенты приведенных уравнений отличаются почти в два раза.

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

Клит = 1,283 - 0,03718 Кпр - 2,114 аи.

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

Кет = 0,8264 КЛ1ГГ + 0,0905, г = 0,79 - терригенные породы;

Квт = 0,3920 Клит + 0,2616, г = 0,77 - карбонатные породы.

Построенная модель позволяет оценить влияние коэффициента литологии Клиг на коэффициент вытеснения Квт для терригенных и карбонатных пород.

Коэффиииент литологии, д.еа.

Рис. 5. Зависимость коэффициента вытеснения от коэффициента литологии

Полученный с помощью метода «метки» набор характеристик микроструктуры коллектора (ам, рм, рм/<7м) не является окончательным для исследования интересующей нас проблемы обусловленности коэффициента вытеснения микроструктурными особенностями коллектора. Возможно расширение этого набора за счет включения в него дополнительных количественных характеристик микроструктуры, которые, в частности, можно получить на основе исследований образцов методом смесимого вытеснения (СВ).

Интерпретация результатов исследования микроструктуры породы-коллектора методом «метки» основана на достаточно упрощенном представ-

лении о поровом пространстве как совокупности ячеек с проточными и застойными зонами. Совершенно очевидно, что различия между проточными и застойными зонами имеют не абсолютный, а относительный характер, а степень «размытости» интересующих нас различий зависит от распределения поровых каналов по размерам. Сведения об этом распределении можно получить, исследуя образцы пород-коллекторов различными методами: растровой электронной микроскопии, ртутной порометрии, капилляриметрии, центрифугирования [7, 8, 15, 20], смесимого вытеснения. Как наиболее полно отражающий динамические свойства порового пространства использовался метод смесимого вытеснения.

Метод смесимого вытеснения предложен Т. Клинкенбергом и развит в более поздних работах Л.М. Марморштейна и Ю.Б. Меклера. Сущность метода состоит в предварительном насыщении изучаемого образца пористой среды жидкостью, химически не взаимодействующей с веществом образца, последующем вытеснении ее другой жидкостью, полностью растворяющейся в первой и практически не отличающейся от нее по плотности и вязкости. Измерив концентрации С(/) вытесняющей жидкости в струе фильтрата как функции времени / и соответствующим образом обработав результаты исследования, нетрудно рассчитать оценки среднего радиуса поровых каналов г и коэффициента Траска ТЛ, равного отношению 75- и 25-процентных квартилей интегральной функции распределения радиусов г поровых каналов и являющегося одной из наиболее употребительных характеристик неоднородности пористых сред.

Отличительной особенностью применяемого варианта метода [5, 6, 10] является использование вытеснения гексана гептаном (бесцветные жидкости с различными показателями преломления).

Возможности метода СВ не исчерпываются нахождением оценок г и 77?. Обратившись к «двухфазной диффузионной модели» пористой среды, предложенной Ю.А. Чизмаджевым, метод позволяет определить долю а застойной фазы порового пространства, скорость р обмена и эффективный коэффициент диффузии £> [5, 6, 10, 19]. Пусть сс - относительная доля застойной фазы в объеме V порового пространства, (I - а)- доля проточной фазы, и - средняя скорость течения жидкости по обеим фазам, с1и с2 - концентрации метки в проточной и застойной фазах соответственно; обозначим V = р/[( 1 - а)у], у =р/ах. Тогда уравнения материального баланса при фильтрации вдоль оси х имеют вид:

л

дс, „9 с, дс, р

—¡- = й——+---

д1 дх2 дх (1 - а)у

(3.8)

при граничных условиях:

с,(х = 0,0 = 1, с1(х = со,0 = 0. (3.9)

Применяя к уравнениям системы (3.8) и граничным условиям (3.9) преобразование Лапласа

Ш = ¡е^т,

сводим систему (3.8) к обыкновенному дифференциальному уравнению 2-го порядка

(3.10)

с граничными условиями:

с,(д: = 0,о) = 1/а, С!(* = оо,о) = 0. (3.11)

В дальнейшем, используя обратное преобразование Лапласа и теорему Коши о вычетах, мы получили [19] аналитическое решение задачи фильтрации двух взаиморастворимых жидкостей в пористой среде с застойной фазой, имеющее вид:

ехр{хи/2Р}

К"

• г

БШ

4Ъ\ г-к/й

+2

м,

вт

XV sг-s¡+v2(s2-*:/£>)

у/ЗУ 1+У2

¿V

(3.12)

где г и V - переменные интегрирования.

Отметим, что формула (3.12) в качестве частных случаев охватывает задачи конвективной диффузии (р = 0) и молекулярной диффузии (р = 0, и = 0). Путем подбора синтетического профиля концентрации, минимально (по величине среднеквадратического отклонения) отличающегося от полученного при исследовании конкретного образца методом смесимого вытеснения, оцениваются следующие параметры этого образца: доля а застойной фазы, скорость р ее обмена с проточными каналами, эффективный коэффициент диффузии £>.

Возможные практические применения формулы (3.12) связаны с тем, что в настоящее время разрабатываются и испытываются методы повышения нефтеизвлечения путем закачки в пласт различных химреагентов [12, 24], растворов полимеров [25,26], реализации вытеснения нефти в условиях, близких к смесимому [32]. Ввиду высокой стоимости реагентов на практике

обычно применяют их оторочки. Оптимальный размер оторочки зависит не только от геометрии межскважинного пространства, гравитационных эффектов, коэффициентов подвижности и т.д., но и от наличия микронеодно-родностей порового пространства, в частности, тупиковых пор или зон с относительно небольшим массообменом, которые называют застойными. Доля таких зон может достигать 60 % объема порового пространства [3]. Поэтому обеспечиваемая формулой (3.12) возможность получения оценки а важна как в теоретическом, так и в прикладном отношении. Характеристики микроструктуры порового пространства могут найти эффективное применение при проектировании, анализе и регулировании разработки коллекторов со сложными условиями нефтеизвлечения [14, 16] и так называемых сложнопо-строенных коллекторов, связанных с продуктивными пластами, в объемах которых имеются «высокопроницаемые канально-дренажные системы» (С.О. Денк, 1999 г.), расчленяющие пласты на занимающие большие доли их объемов блоки пород с более низкой (в пределе, почти нулевой) проницаемостью, которые в процессе разработки способны «подпитывать» канально-дренажные зоны дополнительными порциями углеводородов.

В данной работе среди результатов метода смесимого вытеснения наибольший интерес представляет коэффициент Траска. Комплексирование коэффициента Траска с параметром ам, оцененным методом «метки», позволяет оценивать коэффициенты уравнения Кет = А1-Ьп(Кпода) + А2 без проведения большого количества прямых измерений коэффициента вытеснения на образцах керна.

Перспективы использования результатов исследования структуры порового пространства коллектора методами «метки» и смесимого вытеснения при расчетах коэффициента вытеснения оценены на основе обработки данных по совокупности из 98 образцов керна терригенных и карбонатных коллекторов. Для каждого из образцов были определены проницаемость и доля застойных зон методом «метки». Методом смесимого вытеснения изучены 68 образцов и для них определен коэффициент Траска, для 50 имелись значения коэффициента подвижности и коэффициента вытеснения нефти водой.

Для всех групп исследованных образцов известны уравнения связи коэффициента вытеснения Квт с коэффициентом подвижности Кподв, приведенные в табл. 1, в которой видно, что уравнения всегда имеют вид Кет = А1(Ьп Кподв) + А2, но набор коэффициентов (А1, А2) для каждой группы свой. Сопоставление значений коэффициента А1 с характеризующими соответствующие группы значениями ам показало, что имеется достаточно тесная (г = -0,95) обратная линейная зависимость между А1 и ам (рис. 6, а, значения аббревиатур - см. табл. 1).

Из сопоставления значений коэффициента А2 с соответствующими значениями отношения ТЛ/а„ следует, что А2 связано тесной (г = 0,97) прямой линейной зависимостью с соотношением ТН/ам (рис. 6, 6).

од

0Д5

03 03$

Доля застойных зон а„, д.ея.

0,4

О,<5

ЗЛ 3,6 3,7 3,8 3,9 4 4,1 4,2 43 4,4 4,5 4,6 4,7 4,8 4,9 5

ТК/а„,д.ед.

б

Рис. 6. Зависимость коэффициента А1 от доли застойных зон а„ (а) и коэффициента А2 от (7У?/ам) (б)

Следовательно, формула для определения коэффициента вытеснения в качестве функции а„, 77?, КП0Дв должна иметь вид:

Квх = В1(ам ■ Ьп Кподв) + В2(Ьп Кподв) + ВЗ (7К/схн) + В4. (3.13)

Обработка методом наименьших квадратов данных указанных выше 50 образцов терригенных и карбонатных коллекторов показала, что при размерности коэффициента подвижности мкм2/(мПас) для В1, В2, ВЗ, В4 должны быть приняты такие оценки:

В1 = 0,0050, В2=-0,0590, ВЗ = 0,0123, В4 = 0,1675. (3.14)

Результаты использования уравнения (3.13) и соотношений (3.14) для расчета средних значений коэффициента вытеснения для каждой из рассмотренных групп образцов отражены в табл. 2, где расчетные значения Кв-г(расч) сопоставлены со значениями К8Т, определенными для тех же групп путем моделирования процесса нефтевытеснения на образцах керна. Соответствующая обработка данных табл. 2 показывает, что теснота линейной связи Квт и КВТ(расч) характеризуется коэффициентом корреляции г - 0,92.

Таблица 2

Сравнение экспериментальных и расчетных значений коэффициента

вытеснения

Тектоническая Пласт Кол-во образцов, шт. Коэффициент вытеснения Коэффициент вытеснения

структура экспериментальный расчетный

Квт, д. ед. Квт(расч1> Д- ед.

БКВ В3В4 8 0,641 0,588

БКВ Бш 9 0,556 0,574

БКВ Тл+Бб 14 0,506 0,557

БКВ Т-Фм 10 0,657 0,595

БС Тл+Бб 3 0,505 0,527

сд Тл+Бб б 0,493 0,531

Высокое значение г свидетельствует о том, что коэффициент вытеснения может рассчитываться по определенному в результате интерпретации данных ГИС коэффициенту подвижности с использованием оценок аы и 77?, для получения которых достаточно изучить сравнительно небольшое количество образцов керна методами «метки» и смесимого вытеснения.

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

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

Эксперименты по определению коэффициента вытеснения К^ проводились в филиале ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть» в г. Перми в соответствии с требованиями ОСТа 39-195-86 на моделях пласта - составных образцах из керна изучаемого объекта. Порядок компоновки модели принимался таким, что по направлению вытеснения каждый последующий образец имел меньшую проницаемость. Авторское свидетельство [13] на конструкцию кернодержателя получено Г.П. Хижняком и др. в 1991 г., патент на полезную модель [37] - в 2011 г. Опыты проводились на установках УИПК-1М и УИК с нефтями конкретных пластов. Нагнетание вытесняющей воды проводилось при скоростях, близких к реальным, непрерывно до полного обводнения выходящей жидкости. Коэффициент вытеснения рассчитывался по формуле

у

■у _ ' выт.нефти ' нач

где ^выт.нефти - объем вытесненной из керна нефти, Утч - объем нефти, первоначально содержащейся в керне.

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

Для месторождений Урало-Поволжья значения проницаемости и вязкости нефтей различных одновозрастных объектов разработки меняются в широких пределах. В этих условиях более корректно использование в качестве аргумента параметра «коэффициент подвижности» (КПОДв), равного отношению проницаемости коллектора по нефти КпрН к вязкости нефти (|1Н):

Кподв (4-1)

Ин

••

К_. - 0.0953 Ьп К' _ + 0,6721

Л

м

аа оз о,4

Коэффициент газопроницаемости, мкмг

03

0,7

* оГ

5 0,5 ш

1 СМ

ь

3

а о^ 5

£ .

К.,» 0,0468 Ьп Кпр11 +0,6419 Яг= 0,98

0.05 0,1 0,13

Коэффициент проницаемости по нефти, мкм2

»-

г 1^ = 0,03 Я2= 0,9( 'Ц, +0,756 Я

*

4 0.6

и"

¡«3 8 з

а

Й оз 3

?0.1

0,002 0,«Н ОДИ 0,008 0,01 0,012

Коэффициент подвижности, мкм/мПа-с

0,016

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

Зависимости для продуктивных отложений месторождений имеют общий вид (рис. 7, в) с высокими коэффициентами корреляции [28,30]:

КВТ = А1 ЬпКП0дв + А2. (4.2)

Применение параметра КПОДв для визейских терригенных отложений одного из месторождений позволило перейти от двух разных зависимостей Квт для нефтенасыщенных и битуминозных образцов (рис. 8, а) к единой зависимости (рис. 8,6).

0,7

5

1 0^5

04

I 0,55

К„64,949 К^-23418 К^ + 3,6889К„рг + 0,427«

Я-0,98

К

I

0,5 ■

0,45 ■

О - нефтенасышенные

• - битуминозные

А * битуминозные (соседнее ы-«)

0,4

0,05 0,1 0,15 0,2 0,25 03

Абсолютная газопроницаемость, мкм2

а

0,35

0,4

0,7

$ 0.65

ВС

I Оуб

| 0^5

е оз н

| 0,45

М. 0,4

й 035

03

• _

: У« к„-о,о ПЬпК^! Л2= 0,89 ^и +0,71 31

(

,»— -

О - нефтендсышенныс образцы • - битуминозные образцы . 1 . , , , 1 . ; , ,

ода

0,12

044

0,04 ОД« ОЛв 0,1

Коэффициент подвижности, мкм:/мПаас

б

Рис. 8. Зависимости коэффициента вытеснения от абсолютной газопроницаемости (а), коэффициента подвижности (б) для визейских терригенных отложений месторождения Пермского свода

Для месторождений Пермского края обоснованы и в течение длительного времени опробованы [31, 33] надежные методики оценки К„ как опытным путем в лабораторных условиях по керну, так и на основе установ-

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

Кз^АМпСКпрн/ц^ + Аг, (4.3)

где Кпр„/|д.0- коэффициент подвижности; цо - относительная вязкость, равная отношению вязкости нефти цн к вязкости вытесняющей воды цв. Нетрудно заметить, что параметр цо допускает использование и пресной воды, что является частным случаем, при котором Цо = Цн-

При недостаточном для установления величины коэффициента вытеснения объеме керна изучаемого месторождения привлекаются данные по месторождениям-аналогам [30]. Разработка нефтяных месторождений в Пермском крае осуществляется с 1929 г. Месторождения имеют широкий диапазон геолого-физических характеристик пластов-коллекторов и физико-химических свойств насыщающих их флюидов. За это время накоплен большой объем фактического материала, полученный при лабораторных исследованиях керна и нефтей, что создало основу для применения метода аналогии при обосновании Квт.

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

а) принадлежность к одновозрастным отложениям;

б) зависимости фильтрационно-емкостных характеристик, построенные для этих объектов, обладают высокими значениями коэффициента детерминации.

Близость фильтрационно-емкостных свойств (ФЕС) на примере продуктивных визейских терригенных отложений одного из месторождений иллюстрируется зависимостями пористости (рис. 9, а) и остаточной водонасы-щенности (рис. 9, б) от абсолютной газопроницаемости. Параметры ФЕС находятся в общем поле корреляции, описываются едиными зависимостями с высокими коэффициентами детерминации 0,77 и 0,85. Это позволяет отнести выбранные месторождения к одной группе объектов-аналогов и использовать данные по этим объектам для обоснования коэффициента вытеснения любого из них. На следующем этапе строится зависимость коэффициента вытеснения от коэффициента подвижности для визейских терригенных отложений объектов-аналогов как функция подвижности (рис. 9, в). Распределение имеет вид (4.2) с высоким коэффициентом детерминации Я2 = 0,96. По значениям проницаемости и вязкости нефти рассчитывается коэффициент подвижности, и для него с зависимости снимается значение коэффициента вытеснения пласта в целом или его отдельных объектов разработки.

т ■к» = M865 Knprw< Я2 = 0,77 <6

т

т

1 * *- ' ( < 3 "V > ,'ш А

Ц-

1 1 Пласт Тл2, Юрчукское м-е Пласт Бб, Юрчукское м-е Пласт Мл, Юрчукское м-е Пласт Бб, Чашкинское м-е. Пласт Тл2. Чашкинское м-с

0,(11 ОД 1 10 100 1000 10000

Абсолютная газопроницаемость, мкм:-Ш '

а

Lg KQB - 0,0091(lgKnpr)< - 0,0068(lgKnpr)3- 0,1336(igKnpr) 2 --0,16(lgK|ipr) +1,8875, J?2» 0,85 ]

А < <

O

, 0

- • пласт Бб, Логовское м-е о пласт 1л2, Юрчукское м-е ■ пласт Бб, Юрчукское м-е о пласт Мл, Юрчукское м-е * пласт Бб, Чашкинское м-е д пласт Тл2, Чашкинское м-е | A

» 2 a 1 «

Jo i У >

0,01 0,1 1 10 100 1000 10000 Абсолютная газопроницаемость, мкм *10 "3

б

Коэффициент подвижности, мкм2/мПа-с

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

и месторождений-аналогов

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

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

Разработанные методические подходы к обоснованию коэффициента вытеснения позволили получить сводку уравнений Квт - КпрН/ц0 для 170 месторождений, находящихся в распределенном фонде месторождений группы предприятий «ЛУКОЙЛ-ПЕРМЬ» и принадлежащих 10 крупным структурам Пермского края. Полученные уравнения, некоторые из которых приведены на рис. 10, использовались в отчетах по подсчету запасов и проектных документах на разработку, созданных в 2001-2011 гг. В качестве примера в табл. 3 представлены данные по одной из крупнейших структур - Пермскому своДУВ пятой главе «Обоснование коэффициента вытеснения через проницаемость по ГИС» рассматривается метод оценки коэффициента вытеснения по уравнениям, часть из которых приведена в табл. 3, при отсутствии достоверных сведений о проницаемости по нефти, когда за основу принимается проницаемость по ГИС.

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

Такой подход наиболее актуален для порово-трещинных коллекторов. Определяемая для них по результатам гидродинамических исследований проницаемость по нефти Кпрш как правило, больше, чем проницаемость их поровой части. Если за основу берется проницаемость по ГДИ, то коэффициент вытеснения, снимаемый с зависимости от подвижности и относящийся к поровому типу коллектора, оказывается завышенным из-за более высокого значения проницаемости, характерного для другого, более сложного типа коллектора. К недостаткам использования проницаемости по гидродинамическим исследованиям относится и ограниченность количества гидродинамических исследований скважин.

о,о< о ди. одм о,1 Кюффвцвевт аодвкжяости, мкм2*мПя*с

а

0,14

1^,-0,0198 и Кврд/|1, + 0Д»23

• -пркшаводя ■ — пластовая вода

от одг о/а мм од» еле одп |й»ффнцаевт оадваасяости, нкм'*мПа<с/мПа'С

б

оде

од»

ОА

у ш 0Д)3951л Кпр,/йн + 0,7386

я2* ом

у- 0,03131ДКпр1/Ци +0,7608

о - несторождевкя Берояшкхоп выступа Соликамской лярсссп несторождения ЧервушвиавоИ груиоы Башкирского сват _

0,1

03 03 0,1 Коэффициент палвижмосгы, нкм2'мПа*с

03

од;

Рис. 10. Зависимости коэффициента вытеснения от коэффициента подвижности: а -для башкирских карбонатных и бобриковских терригенных отложений Сибирского месторождения; б - для бобриковских терригенных отложений Шатовского месторождения; в-для визейских терригенных отложений месторождений Березниковско-го выступа Соликамской депрессии и Чернушинской группы Башкирского свода

Таблица 3

Коэффициенты А1 и А2 и диапазоны изменения коэффициентов подвижности и вытеснения для месторождений Пермского свода

Месторождение Пласты ВЗВ4 Коэффициенты уравнения Квт=А1ьп(Кпр„/д0)+А2 Диапазоны изменения параметров

Коэфф подви ициент жности 2 1, мкм Коэффициент вытеснения

А1 0,0412 А2 (Кподв (Квт), лед.

Баклзновское 0,7303 0,01349 0,14022 0,538 0,640

Бш 0,0380 0,7485 0,00027 0.01566 0,430 0.579

Тл2. Тл-Бб 0,0262 0.7242 0,00001 0,28538 0,400 0,701

Горское Тл-Бб, Бб 0,0262 0.7242 0,00001 0,28538 0.400 0.701

Т 0,0399 0,7167 0?00357 0.08668 0.450 0,604

Маячное Бш 0,0380 0,7485 0,00027 0,01566 0,430 0,579

Тл2 0,0262 0,7242 0.00001 0.28538 0.400 0,701

Т 0,0399 0,7167 0.00357 0.08668 0,450 0.604

Межевское Бш 0,0270 0.6755 0,00299 0,03005 0,520 0,602

Тл2 0,0210 0,7261 0,01608 0,14323 0,628 0,681

Осинское Бш, Срп 0.0229 0.6455 0,00061 0,01167 0,464 0,562

ДО 0,0860 0,9174 0?00616 0,06150 0,439 0,672

Полазненское ВЗВ4+Бш 0,0270 0,6755 0,00299 0.03005 0.520 0.602

Тл-ББ 0.0210 0,7261 0,01608 0,14323 0,628 0.681

Т 0,0198 0,6556 0,00400 0,08668 0,528 0,604

Рассветное Бш 0,0380 0,7485 0,00027 0.01566 0,430 0.579

Тл-Бб 0,0262 0,7242 0.00001 0.28538 0.400 0,701

т 0,0860 0.9174 0,00616 0,06150 0,439 0.672

Ярино- Каменноложское Бш 0,0257 0.6739 0,00374 0.03005 0.525 0,602

Тл-Бб 0.0270 0,7131 0,00026 0,13667 0,498 0,660

Т 0,0249 0,6693 0,00400 0,08668 0,528 0.604

Зоринское Тл2а 0,0210 0,7261 0,01430 0,14323 0,628 0.681

ДО 0,0860 0,9174 0,00616 0,06150 0.439 0,672

Козубаевское Тл-Бб 0,0262 0,7242 0,00001 0,28538 0,400 0.701

Кузнецовское Тл2 0,0210 0,7261 0,01608 0,14323 0,628 0.681

т 0,0198 0,6556 0,00400 0,08668 0.528 0,604

ДО 0,0860 0,9174 0,00616 0.06150 0,439 0.672

Куэтымсхое Тл-Бб 0,0270 0,7131 0,00026 0,13667 0,498 0.660

Т 0,0249 0,6693 0,00400 0,08668 0.528 0.604

Лобановское Бш 0,0380 0,7485 0.00027 0,01566 0,430 0,579

Тл2-б, Тл-Бб 0,0262 0,7242 0,00001 0.28538 0,400 0.701

Северокамское ВЗВ4+Бш 0,0131 0.6335 0,00070 0.01653 0.515 0,609

Тл2 0.0210 0,7261 0.01430 0.14323 0.628 0,681

ДО 0,0860 0,9174 0,00616 0,06150 0,439 0,672

Шеметинское Бш 0,0270 0,6755 0,00299 0.03005 0.520 0.602

Тл2, Тл2-а, Тл2-б, Мл 0,0210 0,7261 0,01608 0,14323 0,628 0,681

Т 0,0198 0,6556 0,00400 0.08668 0,528 0.604

Беляевское Тл2-б 0,0262 0,7242 0,00001 0,28538 0.400 0,701

Указанных недостатков можно избежать при определении проницаемости по геофизическим исследованиям Кпр(ГИС). В отличие от гидродинамических геофизическими методами исследования охвачена практически каждая скважина на месторождении. При обосновании К„р(ГИС) [33] используется зависимость «пористость Кп - проницаемость КпрГ», описываемая уравнением вида

Кп = А(Кпрг)л, (5.1)

где КпрГ- абсолютная газопроницаемость образцов керна порового типа.

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

Кпрн = С (КпрГ) • (5.2)

Примеры переходных зависимостей для визейских терригенных и карбонатных отложений месторождений Башкирского свода представлены на рис. 11. Зависимости, построенные по выборкам для одиночных образцов и моделей пласта, обладают очень высокими коэффициентами детерминации: -Л2 = 0,98 - для терригенных пластов Тл, Бб, Мл [33], В? = 0,96 - для карбонатных пластов КВЬ В3В4, Бш, Т-Фм.

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

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

10000

—3 2

Абсолютная газопроницаемость, 10 мкм

10 000

100

Абсолютная газопроницаемость, 10"3 мкм2

Рис. П. Сопоставление фазовой проницаемости по нефти с абсолютной газопроницаемостью для визейских терригенных(д) и карбонатных (б) отложений месторождений Башкирского свода

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

Абсолютам газопроницаемость, 10~3мкм2

Рис. 12. Сопоставление переходных зависимостей для терригенных и карбонатных пород Башкирского свода

Для определения коэффициента вытеснения нефти пласта Т-Фм Сибирского месторождения проведены опыты на моделях коллектора порового типа. Результатом является зависимость (рис. 13), описываемая уравнением:

Кет= 0,0276 Ьп КПОдв+ 0,6768. (5.3)

Л

4

К„ - 0,027« Ьп 1С„ри /щ, + 0,6768 К1 =0,88

О Ш 0,02 0,03 ОМ 0,05 0,0« 0,07 0/Ш

Коэффициент подвижности, мкм^/мПа-с

Рис. 13. Зааисимость коэффициента вытеснения от коэффициента подвижности для турнейско-фаменских отложений Сибирского месторождения

В свою очередь переходная зависимость от проницаемости по ГИС к фазовой проницаемости по нефти КпрИ для турнейско-фаменских отложений месторождений Соликамской депрессии, полученная на основе экспериментальных данных для одиночных образцов и моделей пласта (рис. 14), описывается уравнением [34]:

Кпр„ = 0,158 (Кпр(ГИС))1'266. (5.4)

Расчеты показали, что использование проницаемости по ГДИ ведет к завышению Квт для Сибирского поднятия на 12,5 % относительных, а для Родыгинского - на 7,5 %.

1000

% V

а В

1 м §•

а

I м»

I

0.001

1 1

- Я Архангельское местор. : Гагаринсжое мсстор. А Гежсхос местор. * Логове юзе местор, О Маговское местор. * Мыськнсюе мсстор. О Озерное мсстор. • Ли\товсхое х-е ■ Сибирское местор. А Чашкинсте местор. □ Шсршневсюе мсстор, • Юрчуксное местор. А Унъаянсхое м-с —- линия равных значений пл. Т-Фм 1

Я2- »ЦрГ 0Я ^ ааГ

1 о ■ •

У ' п <И НГ

** •

у' X ........ в1

0,001

0,01

0,1 I 10

Абсолютная газопроницаемость, Ю^мкм2

Рис. 14. Переходная зависимость для турнейско-фаменских отложений месторождений Соликамской депрессии

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

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

Очевидно, что в процессе разработки нефтяного месторождения на КИН влияет большое количество факторов как геологических, обусловлен-

ных особенностями эксплуатационных объектов, так и технологических, определенных условиями разработки [27].

Оценка тесноты связи между коэффициентом вытеснения и КИН проводилась с учетом стадийного характера разработки нефтяных месторождений. Отнесение состояния разработки объекта к какой-либо стадии принимается согласно [29] следующим образом.

Первая стадия разработки - стадия разбуривания скважин основного фонда и освоения системы заводнения для поддержания пластового давления. Основные характеристики: обводнённость 0-13 %, выработка 20 % извлекаемых запасов нефти. Средний темп отбора жидкости составляет от 0 до 2 %.

Вторая стадия разработки — стадия относительно высокого уровня добычи нефти. Характеризуется ростом обводнённости от 13 до 40 %, средний темп отбора жидкости составляет б %, выработка равна 20-45 %.

Третья стадия разработки - стадия значительного падения добычи нефти при разработке нефтяного объекта. Характерна значительная обводненность - 35-55 %, выработка равна 38-65 %. Средний темп отбора жидкости изменяется от 9 до 11 %.

Четвёртая стадия разработки - стадия разработки нефтяного эксплуатационного объекта, характеризующаяся низкими уровнями добычи нефти, постепенным сокращением действующего фонда скважин, при обводнённости больше 50 % и выработке больше 60 %, со значительным падением среднего темпа отбора жидкости.

К поздним стадиям разработки отнесены третья и четвертая стадии.

Сопоставление коэффициентов вытеснения и извлечения нефти проведено для выборки, включающей терригенные и карбонатные продуктивные отложения 135 месторождений Пермского края, разрабатываемых как с поддержанием пластового давления (ППД), так и на естественном режиме (Е), для всех стадий разработки, а также объекты, не разрабатываемые в настоящее время [31].

Рассматривались следующие объекты разработки: верхнекаменноугольные карбонатные верейские и башкирские залежи; визейские терригенные залежи; карбонатные турнейско-фаменские залежи; терригенные девонские залежи.

Для всех залежей генеральной совокупности (п = 666) связь между коэффициентами КИН и Кет существует, но недостаточно тесная - средняя относительная погрешность (Дом) аппроксимации данных линейной зависимостью

КИН = А1-Квх + А2 (6.1)

равна 31,1 %.

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

блюдаются для залежей, разрабатываемых с ППД (144 объекта). При этом положительная корреляция для величины КИН с ростом Кет значительно усиливается на 3-4-й стадиях разработки - средняя относительная погрешность линейной аппроксимации До™ = 12,3 %, на первых двух стадиях разброс данных существенно выше - Дощ =17,1 %. Объясняется это тем, что на 1-2-й стадиях большие, по сравнению с 3-4-й стадиями, отклонения связаны с неточной оценкой КИН. На 3-4-й стадиях КИН, как правило, уточняется, что ведет к усилению связи.

Для 70 эксплуатационных объектов, разрабатываемых на естественном режиме, теснота корреляционной связи КИН-КВТ существенно ниже, что обусловлено кроме прочего большим диапазоном для них изменения коэффициента охвата пластов (Кои,) процессом вытеснения. При этом режиме наибольшее влияние на КИН оказывают продуктивность пластов и плотность сетки скважин.

Аналогичный результат получен при анализе разработки турнейских карбонатных продуктивных отложений месторождений Пермского края. В анализе участвовало 66 залежей 45 месторождений, из них 29 залежей разрабатываются с применением ППД, и 37 - на естественном режиме.

При выделении групп в турнейских залежах, как и в визейских терри-генных, зависимости вида (6.1) имеют различные коэффициенты А1 и А2, характеристики аппроксимации для различных режимов и стадий разработки приведены в табл. 4.

Таблица 4

Характеристики аппроксимации зависимости КИН от Квт уравнением КИН = А1-Квт + А2 для турнейских залежей Пермского края

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

А) А2 Квтшт КггПЖ КИН™, КИНетау

ППД+Е 1-4 66 0,7949 -0,1177 0,42 0,421 0,708 0,100 0.517 18.6

Е 1-4 37 1,1478 -0,3396 0,48 0,494 0,708 0,100 0,517 25,0

ППД 1-4 29 0.5313 0,0526 0.48 0,421 0,647 0,200 0,440 10,7

1-2 12 0.5989 -0.0017 0,28 0.508 0.613 0.200 0,432 17,7

3-4 17 0.6017 0.0198 0,76 0,421 0,647 0283 0,440 6.5

Примечание: ППД - поддержание пластового давления; Е - естественный режим.

Для залежей, разрабатываемых с поддержанием пластового давления и находящихся на 3-й и 4-й стадиях разработки, зависимость КИН от Квт обладает самым низким значением относительной погрешности линейной аппроксимации, равной Дота = 6,5 % (рис. 15).

04 <М5 М

04$

3 ОД

В* 0«

0.15.0,1

оде о

«,4 0,45 0,5 0Д5 ОД 0,65 0,7

С^яла

Рис. 15. Зависимость КИН от коэффициента вытеснения Квт для турнеиских

карбонатных залежей, разрабатываемых с поддержанием пластового давления на 3-4-й стадии

Таким образом, в результате проведенных исследований могут быть сделаны следующие выводы.

1. Важной характеристикой при оценке КИН является коэффициент вытеснения Кет. На примере визейских терригенных и турнейских карбонатных залежей установлены значимые различия для залежей, разрабатываемых с применением ППД и на естественном режиме. Теснота связи КИН и коэффициента вытеснения Квт значительно усиливается при применении систем с ППД.

2. При малом значении свободного члена А2 в линейной зависимости КИН-КВТ коэффициент А1, равный тангенсу угла наклона прямой, характеризует величину коэффициента охвата Кохв. По результатам выполненного анализа КоХВ для визейских терригенных и турнейских карбонатных залежей находится в диапазоне (0,6-0,7), что соответствует фактически достигаемым значениям при существующих технологиях разработки в Пермском крае и свидетельствует о надежности разработанных и реализованных методических подходов к обоснованию коэффициентов вытеснения.

ЗАКЛЮЧЕНИЕ

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

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

~"Даг11 -«V у»<ШПх+Ш9»

г-1 1,76 •

о •

t •

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

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

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

5. Разработан метод оценки коэффициента вытеснения через проницаемость по ГОС. что расширяет сферу использования материалов ГИС, распространяя ее не только на подсчет геологических, но и извлекаемых запасов.

6. Разработаны модели оценки коэффициента вытеснения нефти как функции коэффициента подвижности, определенного с использованием фазовой проницаемости и относительной вязкости, для различных типов коллекторов и различной их тектонической принадлежности, позволяющие рекомендовать значения коэффициента вытеснения для всех 414 пластов 170 месторождений Пермского края, находящихся на 01.01.2012 г. в распределенном фонде группы предприятий «ЛУКОЙЛ-ПЕРМЬ».

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

Разработанные методические подходы успешно реализованы при оценке коэффициента вытеснения для месторождений в других регионах России: Восточно-Придорожное (ООО «ЛУКОЙЛ-Западная Сибирь» ТПП «Кога-лымнефтегаз») [23], Харьяжское (ОАО «ПечорНИПИнефть»), Дачно-Репинское (ООО «Оренбургнефтегаз»), Юрюбчено-Тохомская зона (Красно-ярскКНИИГИМС), Сабардинское и Сухореченское (ОАО Уралнефть г. Екатеринбург), Северо-Алексеевское (ООО «ЛУКОИЛ-Нижневолжскнефть», г. Волгоград) [32]. Помимо названных организаций обоснование коэффициента вытеснения нашло применение в проектных документах на разработку, выполненных в ЗАО «ИНКОНКО» (г. Москва), ОАО «Центральная геофизическая экспедиция» (г. Москва), ООО «Компания нефтяных инженеров» (г. Москва), ОАО «Гипровостокнефтъ» (г. Самара).

Список работ, опубликованных по теме диссертации

Положения диссертации опубликованы примерно в 60 научных работах. Основными являются следующие научные публикации:

1. Тульбович Б.И., Хижняк Г.П. Динамический метод изучения застойных зон порового пространства коллекторов // Геология, разработка, бурение и эксплуатация нефтяных месторождений Пермского Приуралья: сб. науч. тр. / ПермНИПИнефть. -Пермь, 1976.-С. 34-39.

2. A.C. СССР № 706750, кл G01N15/08. Способ определения параметров тупиковых зон порового пространства пород-коллекторов / Б.И. Тульбович, Г.П. Хижняк. № 2646947; заявл. 18.07.78, опубл. 30.12.79. БИ№ 48.

3. Рыжик В.М., Тульбович Б.И., Хижняк Г.П. Метод «метки» определения застойных зон порового пространства пород-коллекторов // Инженерно-физический журнал. -1980. - Т. 39. - № 5. - С. 882-886.

4. Методические рекомендации по определению параметров застойных зон порового пространства пород-коллекторов/сост. В.М. Рыжик, Б.И. Тульбович, Г.П. Хижняк I ПермНИПИнефть. - Пермь, 1980. -17 с.

5. Шурубор Ю.В., Тульбович Б.И., Хижняк Г.П. Определение доли застойных зон порового пространства и характеристик их массообмена при фильтрации взаиморастворимых жидкостей // Современные проблемы и математические методы теории фильтрации: тез. докл. всесоюз. сем. / Ин-т проблем механики АН СССР. - М., 1984. -С. 107-108.

6. Хижняк Г.П. Изучение структуры порового пространства пород-коллекторов методом фильтрации взаиморастворимых жидкостей И Математические методы описания горных пород и расчета их эффективных свойств: сб. науч. тр. / МОИП. - М.: Наука, 1986.-С. 50-58.

7. Комплексное исследование низкопроницаемых терригенных пород в связи с задачами повышения нефтеотдачи / Б.И. Тульбович, О.В. Кочкин, В.П. Митрофанов, Г.П. Хижняк//Ускорение научно-технического прогресса при поиске, разведке и разработке нефтяных и газовых месторождений: тез. докл. / Перм. политехи, ин-т. -Пермь, 1987.-С. 81-82.

8. Анализ структурных групп низкопроницаемых терригенных пород в связи с применением химических реагентов / О.В. Кочкин, В.П. Митрофанов, Б.И. Тульбович, Г.П. Хижняк // Проблемы освоения трудноизвлекаемых запасов нефти Пермского Приуралья: сб. науч. тр./ИГиРШ. -М., 1988. - С. 38-41.

9. Методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом/сост. В.Г. Михневич, Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. - Пермь, 1988. -12 с.

10. Методические рекомендации по изучению структуры порового пространства пород-коллекторов при фильтрации взаиморастворимых жидкостей / сост. Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. - Пермь, 1988. - 22 с.

11. Применение обобщенных петрофизических зависимостей при подсчете балансовых и извлекаемых запасов/Б.И. Тульбович, В.Г. Михневич, В.П. Митрофанов, В.Б. Бейзман, Г.П. Хижняк II Проблемы геологии и разработки нефтяных месторождений в районах с истощающимися ресурсами: сб. науч. тр. / БашНИПИнефть. - Уфа, 1989. - Вып. 79. - С. 117-123.

12. Лабораторные исследования повышения эффективности нефтеизвлечения I З.Р. Борсуцкий, Л.В. Казакова, Б.И. Тульбович, Г.П. Хижняк // Основные направле-

ния научно-технического прогресса в развитии нефтяной промышленности Пермского Прикамья: тез. докл. регион, хонф. / ПермНИПИнефть. - Пермь, 1989. - С. 78-81.

13. А. с. СССР № 1656410, кл <301М5/08. Кернодержатеяь / З.Р. Борсуцкий, О.В. Кочкин, Г.П. Хижняк. № 4483074; заявл. 12.09.88, опубл. 15.06.91. БИ № 22.

14. Особенности петрофизических свойств гидрофобных пород / В.П. Митрофанов, Б.И. Тульбович, О.В. Кочкин, Г.П. Хижняк И Геология нефти и газа. -1992. -№7.-С. 25-28.

15. Митрофанов В.П., Хижняк Г.П., Кочкин О.В. Дифференцированная оценка коэффициента вытеснения низкопроницаемых карбонатных коллекторов // Дифференциация запасов и ресурсов нефти: сб. науч. тр. / ИГиРГИ. - М., 1992.-С. 179-185.

16. Влияние гвдрофобности пород на их петрофизические свойства / В.П. Митрофанов, Б.И. Тульбович, О.В. Кочкин, Г.П. Хижняк // Проблемы повышения эффективности разработки, подсчета запасов и извлечения углеводородов в Пермском Прикамье: сб. науч. тр. / ПермНИПИнефть. - Пермь, 1992. - С. 193-205.

17. Методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом для продуктивных отложений Пермского Приуралья / сост. В.Г. Михневич, Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. - Пермь, 1994.-12 с.

18. Временные методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом для разновозрастных отложений структурно-тектонических элементов Пермского Прикамья / сост. Г.П. Хижняк, В.П. Митрофанов / ПермНИПИнефть. - Пермь, 1996. -15 с.

19. Хижняк Г.П. Определение параметров застойных зон порового пространства при фильтрации взаиморастворимых жидкостей II Геология, разработка, бурение н эксплуатация нефтяных месторождений Пермского Прикамья: сб. науч. тр. / ОАО ПермНИПИнефть.-Пермь, 1999.-Вып. 2.-С. 102-110.

20. О кондициях фильтрационно-емкостных свойств нефтенасыщенных пород / В.П. Мигрофанов, Л.И. Чернова, В.Б. Бейзман, Г.П. Хижняк // Геология, геофизика и разработка нефтяных месторождений. -1999. -№ 10. -С. 14-18.

21. Хижняк Г.П. К вопросу определения смачиваемости пород-коллекторов // Геология, геофизика и разработка нефтяных месторождений. - 2002. - N 8. - С. 44-47.

22. Пат 2216723 РФ, кл. СЮШ13/00, С0Ш13/02, Способ определения смачиваемости пористых материалов / Хижняк Г.П.; заявитель Общество с ограниченной ответственностью «ПермНИПИнефть», патентообладатель Общество с ограниченной ответственностью "ЛУКОЙЛ-Инжиниринг". - № 2002114555/28; заявл. 03.06.2002; опубл. 20.11.2003, Бюл. № 32.

23. К оценке коэффициента вытеснения нефти водой в продуктивных отложениях Западно-Сибирского региона / Митрофанов В.П., Хижняк Г.П. // Геология, геофизика и разработка нефтяных и газовых месторождений. -2004. -Кг 10.-С. 69-75.

24. Изучение изменения проницаемости пород-коллекторов при воздействии культуральной и биоцидсодержащей жидкостей на моделях пласта / Л.В. Закшевская, Г.А. Шарова, Г.П. Хижняк, К.К. Стяжкин, И.Ю. Артемкина И Нефтепромысловое дело. -2006. — № 10. - С. 33-38.

25. Хижняк Г.П, Татаринов И.А., Спасибко А.В. Применение биополимера БП-92 при лабораторном определении коэффициента нефгевытеснения турнейских отложений Ашугайского месторождения // Геология, геофизика и разработка нефтяных и газовых месторождений. -2007. -№ 1. -С. 50-54.

26. Результаты лабораторных исследований и опытно-промысловых работ по закачке полимеров в пласт / A.B. Распопов, Г.П. Хижняк, ВА. Мордвинов, И.Р. Юшков// Геология, геофизика и разработка нефтяных и газовых месторождений. -2007.-№10.-С. 60-63.

27. Галкин C.B., Хижняк Г.П., Ефимов A.A. Влияние геолого-технологических показателей на коэффициенты извлечения нефти (на примере залежей Пермского края) // Научные исследования и инновации / Перм. гос. техн. ун-т. - Пермь, 2008. -Т. 1.- №4. -С. 78-83.

28. Хижняк Г.П., Лядова НА Определение коэффициента нефтевытеснения различными агентами для пород-коллекторов Пермского края // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2008. - № 9. - С. 49-54.

29. Оценка коэффициентов извлечения нефти для месторождений Пермского края на основе статистических моделей / C.B. Галкин, Т.Б. Поплаухина, A.B. Распопов, Г.П. Хижняк //Нефтяное хозяйство. -2009. -№ 4. - С. 38-39.

30. Распопов AB., Хижняк Г.П. Обоснование коэффициента вытеснения нефти водой с привлечением результатов исследований объектов-аналогов // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2009. - № 6. - С. 39-43.

31. Опыт применения методики оценки коэффициента нефтевытеснения при проектировании разработки нефтяных месторождений Пермского края / Г.П. Хижняк, Т.Б. Поплаухина, C.B. Галкин, A.A. Ефимов // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2009. - № 8. -С. 42-45.

32. Хижняк ГЛ., Юркив Н.И. Применение водо- и нефтерастворимых оторочек для увеличения коэффициента нефтеизвлечения на поздней стадии разработки // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2009. - № 9. -С. 48-51.

33. Хижняк Г.П., Распопов AB., Ефимов АА. Методические подходы при обосновании коэффициента вытеснения нефти в различных геолого-физических условиях // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2009. -№ 10.-С. 32-35.

34. Ефимов A.A., Хижняк Г.П. Методические подходы при обосновании коэффициента вытеснения нефти турне-фаменских отложений Сибирского месторождения // Научные исследования и инновации / Перм. гос. техн. ун-т. - Пермь, 2010. -Т. 4.-№1.-С. 21-26.

35. Хижняк Г.П., Распопов AB., Ефимов A.A. Эффективность вытеснения нефти пластовыми водами по данным лабораторных исследований на кернах // Нефтяное хозяйство. - 2011. - № 10.-С. 60-61.

36. Хижняк Г.П. Использование материалов геофизических исследований скважин для оценки проницаемости терригенных пород // Нефтяное хозяйство. - 2011. -№ 6. - С. 86-68.

37. Пат № 111664 РФ, кл. G01N15/08, Кернодержатель / Хижняк Г.П., Юшков В.В.; заявитель и патентообладатель Общество с ограниченной ответственностью "ЛУКОЙЛ-Инжиниринг" - № 2011127983; заявл. 07.07.2011; опубл. 20.12.2011, Бюл. № 35.

Подписано в печать 21.12.2011. Формат 60x90/16. Усл. печ. л. 2,5. Тираж 120. Заказ 246/2011.

Содержание диссертации, доктора технических наук, Хижняк, Григорий Петрович

Список таблиц.

Список рисунков.

Введение

1 Анализ современного опыта оценки коэффициента вытеснения при проектировании разработки нефтяных месторождений.

2 Использование материалов ГИС для количественной оценки свойств пород-коллекторов при подсчете запасов

2.1 О статических и динамических свойствах пород-коллекторов

2.2 Определение проницаемости пород-коллекторов по данным геофизических исследований скважин.

3 Использование результатов исследования структуры порового пространства коллектора в расчетах коэффициента вытеснения.

3.1 Метод «метки» и его использование для определения параметра «доля застойных зон в поровом пространстве»

3.1.1 Развитие модели Ю.А. Чизмаджева.

3.1.2 Физическое моделирование пористой среды с застойными зонами.

3.1.3 Выделение вклада образца в суммарную кривую размывания «метки».

3.1.4 Применение способа «метки» для исследования образцов керна

3.1.5 Методика оценки влияния коэффициента литологии на коэффициент вытеснения терригенных и карбонатных пород

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

3.2.1 Обзор методов исследования структуры порового пространства пород-коллекторов.

3.2.2 Возможности метода смесимого вытеснения

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

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

3.2.5 Прикладное значение результатов исследования образцов пород-коллекторов методом смесимого вытеснения.

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

4 Методические подходы обоснования коэффициента вытеснения нефти.

4.1 Параметр «коэффициент подвижности».

4.2 Эффективность вытеснения нефти пластовыми водами по данным лабораторных исследований на кернах.

4.3 Метод аналогии.

5 Обоснование коэффициента вытеснения через проницаемость по ГИС.

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

Введение Диссертация по наукам о земле, на тему "Комплексное решение проблемы оценки коэффициента вытеснения нефти в различных геолого-технологических условиях"

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

Актуальность работы. В настоящее время в соответствии с требованиями ОСТ 39-195-86 основным методом определения коэффициента вытеснения (Квт) является метод лабораторного моделирования. В качестве моделей пласта используются составные образцы керна изучаемого объекта, насыщенные нефтью и содержащие остаточную воду. Одним из важнейших параметров, влияющих на коэффициент вытеснения, является проницаемость коллектора. Однако при ее обосновании встречается ряд проблем.

В настоящее время проницаемость определяется по керну, по геофизическим (ГИС) и гидродинамическим (ГДИ) исследованиям скважин. Наиболее надежно отражающей пластовые условия считается проницаемость, определенная по гидродинамическим исследованиям. В случае достаточного количества исследований фазовая проницаемость по нефти при ГДИ берется за основу. Однако процент исследованных скважин методом ГДИ в среднем составляет порядка 30 % от общего числа скважин, часто не превышая 10 %.

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

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

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

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

Цели исследований заключаются:

- в использовании материалов ГИС для получения достоверных данных о значениях коэффициентов проницаемости и вытеснения;

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

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

Основные задачи исследований определяются поставленными целями и формулируются следующим образом:

- разработка моделей прогноза проницаемости пород-коллекторов нефтяных месторождений Пермского Прикамья по данным ГИС;

- научное обоснование использования коэффициента подвижности для оценки коэффициента вытеснения в различных геолого-технологических условиях;

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

Научная новизна.

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

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

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

Защищаемые положения:

- метод оценки коэффициента проницаемости пород-коллекторов по данным ГИС с помощью многомерных моделей;

- способ определения застойных зон порового пространства пород-коллекторов методом «метки»;

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

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

Практическая ценность работы заключается:

- в повышении точности оценок извлекаемых запасов нефти за счет использования результатов геолого-петрофизической интерпретации данных ГИС для моделирования пространственного изменения коэффициента вытеснения;

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

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

Реализация результатов исследования. Исследования по теме диссертации выполнялись с 1974 по 2011 г. Результаты исследований реализованы в двух авторских свидетельствах [4, 6], патентах на изобретение [138] и полезную модель [139], ряде методических рекомендаций и нашли применение практически во всех отчетах по подсчету запасов, проектированию и анализу разработки нефтяных месторождений Пермского Прикамья в 1980-2011 гг.

Исходные данные. Для решения поставленных задач выполнены лабораторные исследования более 2500 образцов керна пород-коллекторов, отобранных из более 400 нефтяных залежей северо-востока

Волго-Уральской нефтегазоносной провинции, с привлечением данных по залежам Западной Сибири, Свердловской, Оренбургской областей и Красноярского края методами физического моделирования процесса нефтевытеснения, «метки», фильтрации взаиморастворимых жидкостей.

Апробация работы и публикации. Материалы диссертации докладывались и обсуждались на всесоюзном семинаре «Современные проблемы и математические методы теории фильтрации» /М., 1984 г./; на заседании Межведомственного совета по лабораторным методам, применяющимся в нефтегазовой геологии (группа по коллекторам) /Львов, 1984 г./; на всесоюзном семинаре «Современные проблемы и математические методы теории фильтрации» /М., 1984 г./; на межсекционном семинаре «Математическое описание микронеоднородных сред и расчет их физических свойств» /М., 1985 г./; на научно-технической конференции «Совершенствование методов поисков, разведки и разработки нефтяных и газовых месторождений» /Пермь, 1989 г./; на конференции «Основные направления научно-технического прогресса в развитии нефтяной промышленности Пермского Прикамья» /Пермь, 1989 г./; на краевой научно-технической конференции студентов, аспирантов и молодых ученых «Проблемы комплексного освоения месторождений полезных ископаемых в Пермском крае» /Пермь, 2007 г./; на международной научно-технической конференция «Нефтегазовое и горное дело» /Пермь, 2009, 2010, 2011 гг./; на научно-технической конференции «Актуальные проблемы геологических исследований и разработки месторождений Пермского края» /Пермь, 2010 г./. о

Результаты работ использованы в 5 внедренных в ООО «ЛУКОЙЛ-ПЕРМЬ» методических руководствах [35, 116, 118-120].

Основные положения диссертации опубликованы более чем в 50 научных работах, изданных в Москве, Минске, Уфе, Перми. За выполненные научные разработки автору присужден диплом Пермского обкома комсомола, Обкома профсоюза рабочих нефтяной и газовой промышленности.

Автор приносит глубокую благодарность научному консультанту доктору геолого-минералогических наук, профессору C.B. Галкину.

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

Ю.В. Шурубор|, которым автор выражает искреннюю благодарность.

Искреннюю признательность за консультации и поддержку в работе над диссертацией автор выражает H.A. Лядовой, A.B. Распопову, В.И. Галкину, A.B. Растегаеву, A.C. Некрасову, В.Б. Бейзману, М.Ф. Серкину, В.К. Червяковой. Научным исследованиям способствовала действенная помощь сотрудников Центра исследования керна и пластовых флюидов филиала ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть».

Объем и структура работы. Диссертация состоит из введения, шести глав и заключения, изложена на 259 страницах машинописи, включая 63 рисунка и 31 таблицу. Список литературы состоит из 218 наименований.

Заключение Диссертация по теме "Геология, поиски и разведка горючих ископаемых", Хижняк, Григорий Петрович

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

5. Разработан метод оценки коэффициента вытеснения через проницаемость по ГИС, что расширяет сферу использования материалов ГИС, распространяя ее не только на подсчет геологических, но и извлекаемых запасов.

6. Разработаны модели оценки коэффициента вытеснения нефти как функции коэффициента подвижности, определенного с использованием фазовой проницаемости и относительной вязкости, для различных типов коллекторов и их тектонической принадлежности, позволяющие рекомендовать значения коэффициента вытеснения для всех 414 пластов 170 месторождений Пермского края, стоящих на балансе «ЛУКОЙЛ-ПЕРМЬ» на 01.01.2012 г.

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

Разработанные методические подходы успешно реализованы при оценке коэффициента вытеснения месторождений в других регионах России: Восточно-Придорожное (ООО «ЛУКОЙЛ-Западная Сибирь» ТПП «Когалымнефтегаз»), Харьяжское (ОАО «ПечорНИПИнефть»), Дачно-Репинское (ООО «Оренбургнефтегаз»), Юрюбчено-Тохомская зона (КрасноярскКНИИГИМС), Сабардинское и Сухореченское (ОАО Уралнефть г. Екатеринбург), Северо-Алексеевское (ООО «ЛУКОЙЛ-Нижневолжскнефть», г. Волгоград). Помимо названных организаций обоснование коэффициента вытеснения нашло применение в проектных документах на разработку, выполненных ЗАО «ИНКОНКО» (г. Москва), ОАО «Центральная Геофизическая Экспедиция» (г. Москва), ООО «Компания нефтяных инженеров» (г. Москва), ОАО «Гипровостокнефть» (г. Самара).

Библиография Диссертация по наукам о земле, доктора технических наук, Хижняк, Григорий Петрович, Пермь

1. Авдусин П.П., Цветкова М.А. О структуре поровых пространств песчаных коллекторов. - Нефт. хоз-во. - 1938. - № 6.

2. Адсорбционный метод изучения пористой структуры и установка для изучения адсорбции углеводородных газов // Проблемы нефти и газа Тюмени / Г.М. Ярышев, В.А. Питателев, Г.М. Олейник, С.Н. Данияров -1979.-№43.-С. 29-31.

3. Азаматов В.И., Свихнушин Н.М. Методы изучения неоднородных коллекторов в связи с оценкой запасов нефти и газа. М.: Недра, 1976. -216 с.

4. Авт. свид. СССР № 1656410, G01N 15/08, Кернодержатель / З.Р. Борсуцкий, О.В. Кочкин, Г.П. Хижняк. № 4483074; Заявлено 12.09.88; Опубл. 15.06.91, БИ № 22.

5. Авт. свид. СССР № 1310515, Способ определения коэффициента вытеснения нефти водой / JI.B. Сиротенко, В.П. Потапов. Заявлено 15.01.87; Опубл. 15.05.87, БИ№18.

6. Авт. свид. СССР № 706750, кл GO INI 5/08, Способ определения параметров тупиковых зон порового пространства пород-коллекторов, авторы: Тульбович Б.И., Хижняк Г.П., № заявки 2646947, заявлено 18.07.78, опубликовано 30.12.79 БИ № 48

7. Аманов С.Н. Структура порового пространства пород-коллекторов как критерий оценки коллекторских свойств // НТС Нефтегазовая геология и геофизика. 1969. - Вып. 4.

8. Амелин И.Д., Сургучев М.Л., Давыдов A.B. Прогноз разработки нефтяных залежей на поздней стадии. М.: Недра, 1994.

9. Амикс Дж., Басс Д., Уайтинг Р. Физика нефтяного пласта. М.: Гостоптехиздат, 1962. - 572 с.

10. Анализ разработки Кустовского нефтяного месторождения: отчет по дог.: № 05zl 140; рук. Седорако А.Б.; исполн.: Вахитова P.P. и др. / ЗАО «УфаНИПИнефть». Уфа, 2006.

11. Анализ форм и идентификации индикаторных диаграмм гидродинамического каротажа/ВНИИГИК. Калинин, 1986. - 12 с. - Деп. в ВИНИТИ 5.06.86, № 4136-В86.

12. Ахметов Р.Т. Методика прогноза коэффициента вытеснения нефти водой по данным промысловой геофизики // Геология, геофизика и разработка нефтяных и газовых месторождений. М.: ВНИИОЭНГ, 2011. -№ 3.С. 62-63.

13. Багринцева К.И. Карбонатные породы- коллекторы нефти и газа -М.: Недра, 1977.-219 с.

14. Баишев Б.Т. Функции распределения проницаемости и учет неоднородности пласта при проектировании разработки нефтяных месторождений: Тр. ВНИИ. М., 1980. - Вып. 28. - С. 39-67.

15. Baker L.E. Effects of Dispersion and dead-end Pore Volume in Miscible Flooding. Soc. of Petr. Eng. J., v. 17, no 3, p.p. 219-227 (1977).

16. Batycky J. P., Maini B.B., Fisher V.B. Simulation of Miscible Displacement in Full-Diameter Carbonate Cores. SPEJ, 1982 (October), p.p. 647-657.

17. Берлин A.B., Тульбович Б.И., Хижняк Г.П. Расчет капиллярного давления при центрифугировании // Э.И. Разработка нефтяных месторождений и методы повышения нефтеотдачи / ВНИИОЭНГ. М., 1990.-№8.- С. 7-10.

18. Блюм Р.Г., Путилов М.Ф. Удельная поверхность, нефте- и водонасыщенность карбонатных пород как факторы, определяющие эффективность солянокислотных обработок: Тр. ПФ Гипровостокнефть. -Пермь, 1969. Вып. 4.

19. Боганик В.Н., Медведев А.И., Медведев Ю.А., Пестрикова И.А. Определение продуктивности, дебита и проницаемости для неоднородного пласта по данным ГДИС, ГИС и керна. Изд-во АИС, Каротажник, Тверь, вып. 116-117, 2004. С. 69-77.

20. Богомолова А.Ф., Орлова H.A. К изучению структуры порового пространства // Геология нефти и газа. 1961. - № 12. - С. 46-49.

21. Боксерман A.A., Оганджанянц В.Г., Фазлыев Р.Т. Некоторые результаты вытеснения нефти водой из неоднородной пористой среды: Тр. / ВНИИ. М.: Недра, 1965. - Вып. 28.

22. Борисов Ю.П., Рябинина З.К., Воинов В.В. Особенности проектирования разработки нефтяных месторождений с учетом их неоднородности. М.: Недра, 1976. - 285 с.

23. Борковский A.A., Антонишин Г.И. Структура порового пространства пород-коллекторов Ледяковского нефтяного месторождения (УССР) // НТИ Нефтегазовая геология и геофизика / ВНИИОЭНГ. 1961. -№ 19.

24. Brigham W.E. Mixing Equation in Short Laboratory Cores. Soc. Pet. Eng. J., 1974, vol. 257, p. 91-99.

25. Brook C.S, Pursell W.R. Surface area measurements on sedimentary rock / Trans. AJME, 1952, v. 195, p. 289-296.

26. Брунауэр С. Адсорбция газов и паров. Т. 1. М.: Изд-во иностр. лит., 1948.

27. Венделыптейн Б.Ю. Исследование разрезов нефтяных и газовых скважин методом собственных потенциалов. М.: Недра, 1966. - 206 с.

28. Венделыптейн Б.Ю., Резванов P.A. Геофизические методы определения параметров нефтегазовых коллекторов (при подсчете запасов и проектировании разработки месторождений). М.: Недра, 1978. - 318 с.

29. Воларович М.П., Меклер Ю.Б. К вопросу о соотношении гидродинамического и электрического коэффициентов извилистости для пористых сред // Коллоидный журнал. 1967.- Т. 29. - 654 с.

30. Волобуев Г.П. Прогнозирование гидрогеологических обстановок в нефтегазодобывающих районах. М.: Недра, 1986. - 192 с.

31. Вольфкович Ю.М. Применение методов эталонной порометрии для исследования структурных, поверхностных и влагопоглотительных свойств пористых сред // Коллоидный журнал. 1979. - № 4. - С. 640-648.

32. Гавура В.Е. Геология и разработка нефтяных и газовых месторождений. М.: ВНИИОЭНГ, 1995. - 496 с.

33. Галкин C.B., Кошкин К.А., Поплаухина Т.Б. Анализ структуры фонда эксплуатационных объектов при оперативной оценке остаточных запасов нефти // Геология, геофизика и разработка нефтяных и газовых месторождений. М., 2009. № 10, С. 37-39.

34. Галкин C.B., Поплаухина Т.Б., Распопов A.B., Хижняк Г.П. Оценка коэффициентов извлечения нефти для месторождений Пермского края на основе статистических моделей // Нефтяное хозяйство. 2009. - № 4. С. 38-39.

35. Галкин C.B., Хижняк Г.П., Ефимов A.A. Влияние геолого-технологических показателей на коэффициенты извлечения нефти (на примере залежей Пермского края) // Научные исследования и инновации / Перм. гос. техн. ун-т. Пермь. 2008. - Том 2, - № 4 - С. 78-83.

36. Геофизические исследования скважин и интерпретация данных ГИС: учеб. пособие / В. Н. Косков, Б. В. Косков. Пермь: Изд-во Перм. гос. техн. ун-та, 2007. - 317 с.

37. Геофизические методы изучения под счетных параметров при определении запасов нефти и газа / Б.Ю. Венделыптейн, Г.М. Золоева, Н.В. Царева и др. М.: Недра, 1985. - 248 с.

38. Гиматудинов Ш.К. Нефтеотдача коллекторов. М.: Недра, 1970.

39. Гиматудинов Ш.К. Физика нефтяного и газового пласта. М.: Недра, 1971.-309 с.

40. Гусев В.М., Косков В.Н., Некрасов A.C. Оценка гранулометрического состава кварцевых песчаников по данным гамма-метода // Направления и методика поисково-разведочных работ в Пермском Прикамье: Сб. науч. тр. / ИГиРГИ. М., 1984. - С. 86-90.

41. Дахнов A.B. Исследование связей между физическими свойствами максимально влажных терригенных пород и их проницаемостью: Автореф. дис. канд. геол. минерал, наук. - М., 1975. -24 с.

42. Дахнов В.Н. Геофизические методы определения коллекторских свойств и нефтегазонасыщения горных пород. М.: Недра, 1985. - 310 с.

43. Дахнов В.Н. Интерпретация результатов геофизических исследований разрезов скважин. М.: Недра, 1985. - 365 с.

44. Dean H.A. Soc. Petrol. J., 3, 49 (1963).

45. Дементьев Л.Ф., Акбашев Ф.С., Леонтьев Е.И. О физической сущности параметра Рк // НТС Проблемы нефти и газа Тюмени. Тюмень, 1974.-Вып. 22.

46. Дементьев Л.Ф., Хитров Е.А., Шурубор Ю.В. Применение информационных мер в нефтепромысловой геологии / Тр. ПермНИПИнефть. Пермь, 1974. - Вып. 10. - 155 с.

47. Денк С.О. Карбонатные трещинные коллекторы в Пермском Приуралье // Геология нефти и газа. 1992. - № 11. С. 43-47.

48. Денк О.Э. Моделирование фильтрационных систем сложнопостроенных коллекторов для решения задач нефтегазопромысловой геологии: Автореф. дис. докт. геол.-минерал, наук . Пермь, 1998.

49. Денк С.О. Структура и состояния фильтрационной системы пласта-коллектора. Пермь, 1999.

50. Дерягин Б.В., Захаваева H.H., Талаев М.В. Прибор для определения удельной поверхности порошкообразных и дисперсных тел по сопротивлению течения разреженного воздуха (руководство). М.: Изд-во АН СССР, 1953.-21 с.

51. Деч Г. Руководство к практическому применению преобразования Лапласа. М.: Наука, 1965.

52. Дифференциация запасов нефти и газа в неоднородных коллекторах / М.А. Жданов, В.И. Азаматов, Е.П. Гудков, В.М. Гусев. М.: Недра, 1982. - 176 с.

53. Добрынин В. М., Венделыптейн Б.Ю., Резванов P.A., Африкян А.Н. Промысловая геофизика: Учеб. для вузов. Под ред. д.г.-м.н. В.М. Добрынина, к.т.н.Н.Е. Лазуткиной — М.: ФГУП Издательство «Нефть и газ» РГУ нефти и газа им. И.М. Губкина, 2004. — 400 с.

54. Dullien F., Batra V. Determinaition of the structure of porous média. Ind. and Eng. Chem., 1970, v. 62, p. 25-53.

55. Ефимов A.A., Хижняк Г.П. Методические подходы при обосновании коэффициента вытеснения нефти турне-фаменских отложений Сибирского месторождения // Научные исследования и инновации / Перм. гос. техн. ун-т. Пермь. 2010. - Том 4, - № 1 - С. 21-26.

56. Заболоцкая В. А., Коджаев Ш.Я., Кочетков А. А. Экспериментальные исследования капиллярных явлений в пористых средах / ВНИИ.- М.: Недра, 1965. (Сб. науч. тр.; Вып. 28).

57. Забродин П.И., Раковский H.JL, Розенберг М.Д. Вытеснение нефти из пласта растворителями. М.: Недра, 1968. - 224 с.

58. Иванова М.М. Динамика добычи нефти из залежей. М.: Недра, 1976. - 249 с.

59. Измерение кривых распределения пор по радиусам и по капиллярным давлениям методом эталонной контактной порометрии: Методика / Ю.М. Вольфкович, Е.И. Школьников, В.Е. Сосенкин. М.: Изд-во ин-та электрохимии АН СССР, 1977. - 28 с.

60. Изучение геофизическими методами нефтяных месторождений на поздней стадии разработки / В.А. Кошляк, А.Н. Филонов, В.Ф. Козяр и др. М.: Недра, 1983.

61. Изучение коллекторов нефти и газа месторождений Западной Сибири геофизическими методами / Леонтьев Е.И., Дорогоницкая Л.М., Кузнецов Г.С. и др. М.: Недра, 1974.

62. Исякаев В.А., Хачатурян. Э.В. Анализ форм и идентификации индикаторных диаграмм гидродинамического каротажа / ВНИИГИК. -Калинин, 1986. 12 с. - Деп. в ВИНИТИ 5.06.86, № 4136-В86.

63. Ишкин И.П., Каганер М.Г. Распределение пор по их размерам в пористых телах // Журн. физ. химии. 1950. - № 7. - С. 943-952.

64. Carman Р.С. Fluid flow through granular beds. Trans. Inst. Chem. Eng., v. 15, p. 150, 1937.

65. Klinkenberg L.J. Analogy between diffusion and electrical conductivity in porous rock. Bull. Geol. Amer., vol. 62, p. 559, 1951.

66. Klinkenberg L.J. Pore size distribution of porous media and displacement experiments with miscible liquids. J. Petrol. Technol., 1957, vol. 9, no 4, p. 63-66.

67. Coats K.H. and Smith B.D. Dead-End Pore Volume and Dispersion in Porous Media. Soc. Pet. Eng. J., 1964, vol. 231, p.73-84.

68. Колоскова М.И. Применение ртутного поромера для изучения структуры пород-коллекторов: Тр. / ВНИИГаз. М.: Недра, 1964.- Вып. 20/28.- С. 65-71.

69. Комаров B.JI. Петрофизические основы повышения эффективности геофизических исследований скважин на нефтяных месторождениях Восточной окраины Русской платформы: Автореф. дис. докт. геол. минерал, наук. - М., 1971. - 36 с.

70. Комплексная оценка состояния и работы нефтяных скважин промыслово-геофизическими методами: учеб. пособие / В. Н. Косков,

71. Б. В. Косков, И.Р. Юшков. Пермь: Изд-во Перм. гос. техн. ун-та, 2010. -226 с.

72. Контроль за разработкой залежей нефти и газа геофизическими методами: учеб. пособие / В.Н. Косков. Пермь: Изд-во Перм. гос. техн. ун-та, 2009. - 77 с.

73. Котяхов Ф.И. Количественная оценка структуры поровых каналов пород нефтяных и газовых залежей // НТЖ Нефтепромысловое дело / ВНИИОЭНГ. 1962 . - № 1.

74. Котяхов Ф.И. Основы физики нефтяного пласта. М.: Недра, 1956. -359 с.

75. Котяхов Ф.И. Результаты анализа методов определения среднего радиуса пор пород нефтяных залежей // Новости нефтяной и газовой техники. Сер. Нефтепромысловой дело. 1962. - № 2. - С. 37-42.

76. Кофф Т.П., Карша В.О., Капранов Л.Н. Электронномикроскопическое изучение тонкой пористости грунтов // Разведка и охрана недр. 1973. - № 5.

77. Коцеруба Л.А. О применении центрифуги для определения содержания связанной воды и измерения капиллярных давлений: Тр. / ВНИГНИ. М., 1970. - Вып. 90. - С. 193-203.

78. Коцонис А.Н. О некоторых результатах изучения структуры порового пространства горных пород: Тр. / ВНИИЯГГ. М.: Недра, 1972. -Вып. 12.

79. Кочкин О.В., Хижняк Г.П. Определение коэффициентов проницаемости и вытеснения методом центрифугирования //

80. Совершенствование методов поисков, разведки и разработки нефтяных и газовых месторождений: Тез. докл. науч.-техн. конф. / Перм. политехи, инт. Пермь, 1989. - С. 69.

81. Коэффициент вытеснения нефти водой по образцам керна продуктивных горизонтов месторождений Татарии: отчёт; исполн.: Коцюбинский В.Л., Юдинцев Е.А., Веденин C.B. / ТатНИПИнефть, -Бугульма, 1990.

82. Kramers H., Alberda G. Chem. Eng. Sei., 2, 173, (1953).

83. Куликова А.Г. Изучение структуры порового пространства карбонатных коллекторов // НТС Нефтегазовая геология и геофизика / ВНИИОЭНГ. 1965. - № 7.

84. Кусаков М.М., Мекеницкая Л.И. Исследование состояния связанной воды на моделях газовых и нефтяных коллекторов: Тр. МИНХиГП. М.: Гостоптехиздат, 1959. - Вып. 25.

85. Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. М.: Физматгиз, 1965.

86. Leverett M.С. Capillary behaviour in porous solids. Trans. AIME, vol. 142, p. 152, 1941.

87. Левич В.Г., Маркин B.C., Чизмаджев Ю.А. // ДАН СССР. 166, 1401 (1966).

88. Левич В.Г., Маркин B.C., Чизмаджев Ю.А. // ДАН СССР. 168, 1364 (1966).

89. Levich V.G., Markin V.S., Chismadzhev Yu.A. Chem. Eng. Sei., 22, 1357,(1967).

90. Левич В.Г., Письмен C.H., Кучанов С.И. // ДАН СССР. 168, 392 (1966).

91. Леворсен А.И. Геология нефти и газа. М. - Л.: Гостоптехиздат, 1953.-671 с.

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

93. Литология природных резервуаров с использованием фациально-циклического метода и промыслово-геофизических данных: учеб. пособие / В.И. Пахомов, В.Н. Косков. Пермь: Изд-во Перм. гос. техн. ун-та, 2011. - 168 с.

94. Лыков A.B. Теория теплопроводности. М.: Гос. изд. техн.-теорет. лит., 1952. - 392 с.

95. Мамедов Т.М. Добыча нефти с применением углеводородных растворителей. М.: Недра, 1984. - 152 с.

96. Марин А.Р. Влияние механических примесей на приемистость нагнетательных скважин и охват пласта заводнением // ЭИ Нефтепромысловое дело / ВНИИОЭНГ. 1973. - № 18.

97. Марморштейн Л.М. Коллекторские и экранирующие свойства осадочных пород при различных термодинамических условиях. Л.: Недра, 1975. - 158 с.

98. Марморштейн Л.М. Методические рекомендации по изучению физических свойств горных пород. Л.: НИИГА, 1976. - 43 с. (ротапринт)

99. Медер Г.О. О вычислении проницаемости по пористости и гранулометрическому составу. Пер. с нем., № 1914.- М.: ВНИИгаз, 1967.

100. Меклер Ю.Б., Марморштей JI.M. Изменение структуры порового пространства песчаных коллекторов под давлением // Изв. АН СССР. Физика земли. 1971. - № 8. - С. 82-87.

101. Messer E.S. Interstitial water determination by an evaporation method. Petrol. Trans. AIME, 1951, vol. 192, p. 457-463.

102. Методические рекомендации по изучению структуры порового пространства пород-коллекторов при фильтрации взаиморастворимых жидкостей / Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. Пермь, 1988.-22 с.

103. Методические рекомендации по определению коэффициента вытеснения нефти расчетным способом / В.Г. Михневич, Б.И. Тульбович / ПермНИПИнефть. Пермь, 1980. - 12 с.

104. Методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом / В.Г. Михневич, Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. Пермь, 1988. - 12 с.

105. Методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом для продуктивных отложений Пермского Приуралья / В.Г. Михневич, Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. Пермь, 1994. - 12 с.

106. Методические рекомендации по определению параметров застойных зон порового пространства пород-коллекторов / В.М. Рыжик, Б.И. Тульбович, Г.П. Хижняк / ПермНИПИнефть. Пермь, 1980. - 17 с.

107. Микерин Б.П. Нагнетание сжиженных углеводородных газов в нефтяные пласты. М.: Гостоптехиздат, 1961.

108. Митрофанов В.П. Особенности фильтрационно-емкостных свойств карбонатных коллекторов Соликамской депрессии. М.: ОАО «ВНИИОЭНГ». - 2002. - 116 с.

109. Митрофанов В.П. Оценка насыщенности с учетом керновых и геофизических данных // Особенности геологии и разработки нефтяных месторождений Пермского Приуралья: Сб. науч. тр. / ИГиРГИ. М., 1981. -С. 60-62.

110. Митрофанов В.П., Хижняк Г.П., Кочкин О.В. Дифференцированная оценка коэффициента вытеснения низкопроницаемых карбонатных коллекторов // Дифференциация запасов и ресурсов нефти: Сб. науч. тр. / ИГиРГИ. М., 1992. - С. 179-185.

111. Моченов C.B. Применение статистических методов исследования к анализу микроструктур: Тез. докл. II Всесоюзн. школы-семинара. М., 1971.

112. Некрасов A.C. Оценка продуктивности карбонатных коллекторов сложного строения по данным ГИС // Геология, разработка, бурение и эксплуатация нефтяных месторождений Пермского Прикамья. Пермь, 1999. С. 317-319.

113. Некрасов A.C. Методика выделения и оценки вторичной пористости карбонатных пород-коллекторов смешанного типа. Сб. науч. тр. / ПермНИПИнефть. Пермь, 2003. С. 17-30.

114. О кондициях фильтрационно-емкостных свойств нефтенасыщенных пород / В.П. Митрофанов, Л.И. Чернова, В.Б. Бейзман, Г.П. Хижняк // Геология, геофизика и разработка нефтяных месторождений. 1999. - № 10. - С. 14-18.

115. Определение флюидонасыщенности горных пород / В.М. Новоселицкий, A.A. Маловичко, М.С. Чадаев, A.B. Горожанцев // ЭИ Региональная, разведочная и промысловая геофизика / ВИЭМС. 1982. - № 9.

116. Особенности петрофизических свойств гидрофобных пород / В.П. Митрофанов, Б.И. Тульбович, О.В. Кочкин, Г.П. Хижняк // Геология нефти и газа. 1992. - № 7. - С. 25-28.

117. ОСТ 39-195-86. Нефть. Метод определения коэффициента вытеснения нефти водой в лабораторных условиях. Взамен ОСТ 39070-78; Введ. 01.01.87. УДК 665.61.001.4. Группа А29.

118. ОСТ 309-204-86. Нефть. Метод лабораторного определения остаточной водонасыщенности коллекторов нефти и газа по зависимости насыщенности от капиллярного давления. Введ. 01.01.87. УДК 553.98: 543.06. Группа А29.

119. Панченко Е.В. и др. Лаборатория металлографии. М.: Металлургиздат, 1957.

120. Parcalabesoh J.D., Olteann Doina. Structura spatinlui poros al colectorelor de hidrocarburi fluide observata la microscopul electronic eu baleaj. Mine, petrol si gaze, 1976, v. 27, p. 81-87.

121. Parcell W.R. Capillary pressures their measurement using mercury and the calculation of permeability therefrom. Trans. AIME, vol. 186, 1949.

122. Пат № 111664 РФ, кл. G01N15/08, Кернодержатель/Хижняк Г.П., Юшков В.В.; заявитель и патентообладатель Общество с ограниченной ответственностью "ЛУКОЙЛ-Инжиниринг" № 2011127983; заявл. 07.07.2011; опубл. 20.12.2011, Бюл. № 35.

123. Методические рекомендации по подсчету геологических запасов нефти и газа объемным методом. Под ред. В.И. Петерсилье, В.И. Пороскуна, Г.Г. Яценко. Москва-Тверь: ВНИГНИ, НПЦ «Тверьгеофизика», 2003.

124. Пирсон С.Д. Учение о нефтяном пласте. М.: Гостоптехиздат, 1961.-570 е.: ил.

125. Пирсон С.Д. Учение о нефтяном пласте. М.: Недра, 1971.

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

127. Порай-Кошиц Е.А. Рентгенографический метод изучения строения высокодисперсных веществ // Методы исследования структуры высоко дисперсных и пористых сред. М.: Изд. АН СССР, 1953.

128. Плаченов Т.Г., Карельская В.Ф., Пулеревич М.Я. Изучение структуры пористых сред методом вдавливания ртути // Методы исследования структуры высокодисперсных и пористых тел. М.: Изд-во АН СССР, 1958. - С. 251-258.

129. Постников Е.В., Жидовинов С.Н. Современные тенденции развития ресурсной базы углеводородного сырья и пути повышения результативности геолого-разведочных работ в Урало-Поволжском регионе // Геология нефти и газа. 2008. - № 5. - С. 2-10.

130. Привалов И.И. Введение в теорию функций комплексного переменного. М.: Физматгиз, 1960. - 444 с.

131. Проектирование разработки нефтяных месторождений / П.А. Крылов, П.М. Белаш, Ю.П. Борисов и др. М.: Гостоптехиздат, 1962. - 430 с.

132. Распопов A.B., Мезрин В.Н., Юшков И.Р. Применение метода аналогии при анализе и проектировании разработки залежей углеводородов // Нефтегазовое и горное дело. Вестник ПГТУ. Вып. 6 / Перм. гос. техн. ун-т. Пермь, 2005. - С. 73-76.

133. Распопов A.B., Хижняк Г.П. Обоснование коэффициента вытеснения нефти водой с привлечением результатов исследований объектов-аналогов // Геология, геофизика и разработка нефтяных и газовых месторождений. 2009. - № 6. С. 39-43.

134. РД 39-0147035-214-86. Методическое руководство по расчету коэффициентов извлечения нефти из недр.- М., 1986.

135. РД 39-0147276-002-ВНИИ-86. Руководство по применению метода воздействия на пласт растворами полимеров при вытеснении нефтиминерализованными водами на месторождениях Башкирии. Уфа, 1986. -46 с.

136. Регуш В. А. К методике изучения структуры порового пространства и определения проницаемости терригенных коллекторов по шлифам: Тр. / ТатНИИ. Л.: Недра, 1968. - Вып. 12.

137. Результаты лабораторных исследований и опытно-промысловых работ по закачке полимеров в пласт / Распопов A.B., Хижняк Г.П., Мордвинов В.А., Юшков И.Р. // Геология, геофизика и разработка нефтяных и газовых месторождений. 2007. - № 10. - С. 60-63.

138. Ritter H.L., Drake L.C. Pore-size distribution in some typical porous substances. Ind. Enq. Chem., 1945. v. 17. p. 781-791.

139. Rose W., Bruce W. Evaluation of capillary character of petroleum reservoir rock. Trans. AJME, 1949, v. 186, p. 127-142.

140. Рухин А.Б. Основы литологии. M. - Л.: Гостоптехиздат, 1953.671 с.

141. Рыжик В.М., Тульбович Б.И., Хижняк Г.П. Метод «метки» определения застойных зон порового пространства пород-коллекторов // Инженерно-физический журнал. 1980. - Т. 39. - № 5. - С. 882-886.

142. Сабиров И.Х. Обзор лабораторных исследований по изучению влияния физических и физико-химических свойств пластовых систем и условий вытеснения на механизм нефтеотдачи пласта: Тр. / УФНИИ. М.: Недра, 1969. - Вып. 18.

143. Савич А.И. Оценка коэффициента извлечения нефти на стадии промышленной разведки и первого подсчета запасов нефти залежей призаводнении // Вопросы освоения нефтяных залежей Пермского Приуралья / Тр. ВНИИОЭГН, М., 1990, с. 9-13.

144. Свихнушин Н.М., Азаматов В.И. Методика изучения неоднородных коллекторов при оценке кондиций и подсчете запасов. М.: Недра, 1970. - 140 с.

145. Сиротенко JI.B. Об оценке проницаемости терригенных пород по данным каротажа // Геофизические методы поисков и разведки нефти и газа: Межвед. сб. науч. тр. / ПТУ. Пермь, 1979. - С. 143-145.

146. Сиротенко JI.B. О выделении проницаемых пород по данным каротажа при помощи многомерного статистического анализа // Вопросы обработки и интерпретации геофизических аномалий: Межвед. сб. науч. тр. / ПТУ. Пермь, 1977.-С.133-135.

147. Сиротенко JI.B., Потапов В.П. Оценка коэффициента вытеснения нефти водой по образцам пород, отобранных сверлящими керноотборниками // Геология, поиски и разведка нефтяных и газовых месторождений: Межвуз. сб. науч. тр. / ППИ. Пермь, 1987. - С. 64-69.

148. Slobod R., Chambers A. Prehn W. Use of centrifuge for determining connate water, residual oil and capillary pressure curves of care samples. -Petrol. Technol., 1951, v. 192, p. 127-133.

149. Cottschlich C.F., Axial Dispersion in a Packed Bed. AJCBE Jour. 1963, p. 88-92.

150. Сургучев M.JI., Горбунов A.T., Жданов С.A. и др. / Геолого-физические условия эффективного применения методов увеличения нефтеотдачи пластов // Нефтяное хозяйство 1974. - №4. - С. 29-34.

151. Сургучев М.Л. / Вторичные и третичные методы увеличения нефтеотдачи пластов. М.: Недра, 1985. 308 с.

152. Сургучев М.Л. Маслянцев Ю.В. Характеристика процесса заводнения пластов с учётом микронеоднородности при избирательной фильтрации: Тр. / ВНИИ. М.: Недра, 1965. - Вып. 28.

153. Сургучев M.JI. Методы контроля и регулирования процесса разработки нефтяных месторождений. М.: Недра, 1968. - 300 с.

154. Таблицы вероятностных функций. М.: ВЦ АН СССР, 1970. - Т. 1. - Вып. 2.

155. Thomas A., Alba P., Conraund G. Contribution and methods detude du milien poreux. Rev. Inst, franc, petrol., 1961, v. 16, p. 36-47.

156. Тульбович Б.И., Берлин A.B., Гусев В.М. Оценка смачиваемости поверхности пород по данным измерения сопротивления // Особенности геологии и разработки нефтяных месторождений Пермского Приуралья: Сб. науч. тр. / ИГиРГИ. М, 1981. - С. 68-72.

157. Тульбович Б.И. Изучение структуры пористых сред методом рефрактометрии взаиморастворимых жидкостей // Журн. физ. химии, 1974. -Т. 48. -№ 11.-С. 2875-2876.

158. Тульбович Б.И. Коллекторские свойства и химия поверхности продуктивных пород. Пермь, 1975. - 194 с.

159. Тульбович Б.И. Методы изучения пород-коллекторов нефти и газа. М.: Недра, 1979. - 199 с.

160. Тульбович Б.И. Метод смесимого вытеснения неокрашенных красителем жидкостей для изучения структуры пород-коллекторов // Геология нефти и газа. 1975. - № 10. - С. 29-33.

161. Тульбович Б.И., Хижняк Г.П. Динамический метод изучения застойных зон порового пространства коллекторов // Геология, разработка, бурение и эксплуатация нефтяных месторождений Пермского Приуралья: Сб. науч. тр. / ПермНИПИнефть. Пермь, 1976. - С. 34-39.

162. Тульбович Б.И., Шустеф И.Н. К вопросу о зависимости коэффициента вытеснения от степени смачиваемости поверхности коллектора: Тр. / ПермНИПИнефть. Уфа,1973. - Вып. 8. - С. 91-95.

163. Turner G.A., The Flow Structure in Packed Beds. Chem. Eng. Sci. 1958, p. 156-165.

164. Turner J.C. Brit. Chem. Eng., v. 9, p. 376 (1964).

165. Тюрин В.И. и др. Разработка методики статистического исследования микроструктуры горных пород по фотографиям шлифов на ЭВМ: Тез. докл. 11 Всесоюзн. школы-семинара. М., 1971.

166. Фионов А.И. Особенности строения переходных зон газонефтяных залежей в физических полях электрического тока и пластовых давлений. Изд-во АИС, Каротажник, Тверь, вып. 6, 2008. С.42-59.

167. Хабибулин З.А. Регулирование свойств смолистой и парафинистой нефти в пласте для совершенствования процессов разработки месторождений при заводнении: Автореф. дисс. Уфа, 1988.

168. Ханин А.А. Петрофизика нефтяных и газовых пластов. М.: Недра, 1976. - 294 с.

169. Ханин А.А. Породы-коллекторы нефти и газа и их изучение. -М.: Недра, 1969. 366 с.

170. Хижняк Г.П. Использование материалов геофизических исследований скважин для оценки проницаемости терригенных пород // Нефтяное хозяйство. 2011. - № 6. С. 86-88.

171. Хижняк Г.П. К вопросу изучения застойных зон порового пространства пород-коллекторов // Тез. докл. V науч. техн. конф. молодых ученых и специалистов / ПермНИПИнефть. - Пермь, 1978. - С. 56-57.

172. Хижняк Г.П. Изучение структуры порового пространства пород-коллекторов методом фильтрации взаиморастворимых жидкостей //

173. Математические методы описания горных пород и расчета их эффективных свойств: Сб. науч. тр. / МОИП. М.: Наука, 1986. - С. 50-58.

174. Хижняк Г.П. Петрофизические исследования динамических особенностей структуры порового пространства пород-коллекторов в связи с проблемами нефтеизвлечения (на примере залежей Пермского Прикамья): Дисс. на соиск. степени канд. техн. наук. Пермь, 2000.

175. Хижняк Г.П. Определение параметров застойных зон порового пространства пород-коллекторов при фильтрации взаиморастворимых жидкостей // Тез. докл. IX науч. техн. конф. молодых ученых и специалистов / ПермНИПИнефть. - Пермь, 1983. - С. 44-45.

176. Хижняк Г.П. Изучение структуры порового пространства пород-коллекторов методом фильтрации взаиморастворимых жидкостей // Математические методы описания горных пород и расчета их эффективных свойств: Сб. науч. тр. / МОИП. М.: Наука, 1986. С. 50-58.

177. Хижняк Г.П., Лядова H.A. Определение коэффициента нефтевытеснения различными агентами для пород-коллекторов Пермского края // Геология, геофизика и разработка нефтяных и газовых месторождений. 2008. - № 9. - С. 49-54.

178. Хижняк Г.П., Распопов A.B., Ефимов A.A. Методические подходы при обосновании коэффициента вытеснения нефти в различных геолого-физических условиях // Геология, геофизика и разработка нефтяных и газовых месторождений. 2009. - № 10. - С. 32-35.

179. Хижняк Г.П., Распопов A.B., Ефимов A.A. Эффективность вытеснения нефти пластовыми водами по данным лабораторных исследований на кернах // Нефтяное хозяйство. 2011. - № 10. - С. 60-61.

180. Хижняк Г.П., Юркив Н.И. Применение водо- и нефтерастворимых оторочек для увеличения коэффициента нефтеизвлечения на поздней стадии разработки // Геология, геофизика и разработка нефтяных и газовых месторождений. 2009. - № 9. - С. 48-51.

181. Центробежный метод нахождения функции распределения пор по размерам в пористых телах / О.С. Ксенжек, Е.А. Калиновский, С.А. Петрова, В.Е. Литвинова // Журн. физ. химии. 1967. - № 7. - С. 1602-1607.

182. Чейз Ф. Количественно-минералогический анализ шлифов под микроскопом. М.: Изд-во ИЛ, 1963.

183. Черников O.A. Преобразование песчано-алевролитовых пород и их пористость. М.: Наука, 1969. - 120 с.

184. Чизмаджев Ю.А. и др. Макрокинетика процессов в пористых средах. М.: Недра, 1971. - 364 с.

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

186. Шейдеггер А.Э. Физика течения жидкостей через пористые среды. М.: Гостоптехиздат, 1960.

187. Ширковский А.И. Определение и использование физических параметров пористой среды при разработке газоконденсатных месторождений // НТО Геология, разведка и разработка газовых и газоконденсатных месторождений / ВНИИЭгазпром. М., 1971.

188. Шурубор Ю.В. Основы теории и методы геологического обеспечения разработки нефтяных месторождений с применением математического моделирования и ЭВМ: Автореф. дис.докт. геол.-минерал. наук. М., 1985. - 36 с.

189. Шурубор Ю.В. Учебная дисциплина «Геометризация нефтяных и газовых залежей»: Учебно-метод. пособие / Перм. гос. техн. ун-т. Пермь, 1999.-57 с.

190. Элланский М.М. Петрофизические связи и комплексная интерпретация данных промысловой геофизики. М., Недра, 1978. - 215 с.

191. Энгельгардт В. Поровое пространство осадочных пород. М.: Недра, 1964. - 220 с.

192. Юркив Н.И. Физико-химические основы нефтеизвлечения. М.: ОАО «ВНИИОЭНГ», 2005. - 336 с.