Анализ и оптимизация технологического режима работы добывающей скважины

Анализ и оптимизация технологического режима работы добывающей скважины

Министерство образования и науки РФ

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования

Пермский Национальный Исследовательский Политехнический Университет

Кафедра Нефтегазовые Технологии

Курсовой проект по дисциплине:

Скважинная добыча нефти в осложненных условиях

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

Пермь, 2015

Исходные данные для курсового проектирования по дисциплине «Скважинная добыча нефти в осложненных условиях»

Общие данные: =0,146 м; =1 МПа; =10 % (масс.);

= 800 м; =1 МПа; =50 % (масс.);

=200 м; =30 м; β=0,85 (УСШН);

α = 25 °; =6 °С; =0,1 (УЭЦН);

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

В связи с ухудшенным состоянием призабойной зоны пласта (rпзп=5м), о котором говорит скин-фактор S=3, рекомендуется провести кислотную обработку.

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

Глиняной кислотой называется смесь 3-5%-й фтористо-водородной (HF) и 8-10%-й соляной кислот. Терригенные коллекторы, как правило, содержат малое количество карбонатов, в среднем 1-5% по массе. Основная масса таких коллекторов представлена силикатными веществами (кварц) и алюмосиликатами (каолин).

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

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

SiF4 + 4H20 = Si(0Н)4 + 2Н2SiF6

Кремнефтористоводородная кислота Н2SiF6 остается в растворе, а кремниевая кислота Si(0Н)4 при понижении кислотности раствора может образовывать гель кремниевой кислоты, выпадающий в осадок и закупоривающий призабойную зону.

Реакция алюмосиликатов с HF такова:

H4Al2Si209 + 14HF = 2A1F3 + 2SiF4 + 9Н20

Соль фтористого алюминия A1F3 остается в растворе, а фтористый кремний SiF4, соединяясь с водой образует Si(0Н)4 и Н2SiF6.

Таким образом, общая реакция следующая:

H4Al2Si209 + 12HF = 2A1F3 + Si(0Н)4 + Н2SiF6 + 5Н20

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

CaCO3 + HF = CaF2 + H2O + CO2

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

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

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

Закачивают расчетный оббьем кислотного раствора в скважину.

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

Нейтрализация кислотного раствора за счет реагирования его с обрабатываемой породой. Время нейтрализации зависит от давления и температуры и изменяется от 1 ч. до 36 ч.

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

2. Техника, применяемая при кислотной обработке

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

Закачка кислотных растворов в скважину осуществляется специальными насосными агрегатами на автомобильном шасси, например, «Азинмаш 30 А». Насосный агрегат включает в себя гуммированную цистерну для кислотного раствора, насос высокого давления с приводом от коробки отбора мощности автомобиля. Давление, создаваемое таким насосом, меняется от единиц до десятков МПа. Наряду с насосным агрегатом, для кислотных обработок используется цементировочный агрегат (типа ЦА-320М), который играет роль подпорного насоса для основного агрегата, подавая технологические жидкости на прием силового насоса. Кроме того, агрегат ЦА-320М, оборудованный насосом низкого давления и емкостями, позволяющими перемешивать кислотный раствор с различными реагентами, добавляемыми в него на скважине.

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

Рис. 1 - Примерная схема обвязки наземного оборудования при глинокислотной обработке: 1 - передвижная емкость для кислоты; 2 - стационарная емкость для кислоты; 3 - емкость для нефти; 4 - цементировочный агрегат; 5 - установка насосная АзИНМАШ-ЗОА; 6 - бункеры; 7 - основной насос; 8 - водяной насос; 9 - резервуар; 10 - насос; 11 - скважина

Расчет распределения плотности ГЖС (), удельного газосодержания жидкости (βг) и температуры (tскв) в интервале «забой- прием скважинного насоса». Построение зависимости для , βг, tскв.

УСШН

а) =10 % (масс.)

Расчет производим «снизу-вверх».

)Для расчета f принимаем

А для расчета dp/dH,  принимаем:

Зададимся 10 интервалами, тогда

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

) Определим температуру на затрубе скважины

) Определяем температуру потока, соответствующую заданным давлениям

) Рассчитаем текущее равновесное давление насыщения при :

где  и  - соответственно молярные доли метана и азота в попутном газе

Массовый состав смеси принимаем

) Рассчитываем относительную плотность газа по воздуху:

Молекулярная масса газа:

Вычисляем молекулярную массу данного газа:

) Определим объем выделившегося газа

, где

) Определим остаточную газонасыщенность нефти в процессе разгазирования (Р=1,5 МПа; Т=294 К)

) Определим относительную плотность выделившегося газа

, где

) Определим относительную плотность растворенного газа, остающегося в нефти при данных условиях ее разгазирования

) Определим удельное приращение объема нефти за счет единичного изменения ее газонасыщенности  и объемный коэффициент нефти

12) Определим коэффициент сверхсжимаемости газа

а) Определим Рпр и Тпр

Таблица 1 - Данные о критических давлениях и температуре различных компонентов

б) Определим коэффициент сверхсжимаемости для азота

в) Определим коэффициент сверхсжимаемости углеводородной составляющей газа

т.к Рпр находится в интервале 0≤Рпр≤3,8 ,а Тпр - в интервале 1,17≤Тпр≤2,0, то

г) Определим коэффициент сверхсжимаемости газовой смеси

) Определяем удельный объем газожидкостной смеси

) Определим удельную массу смеси при стандартных условиях

) Определяем корреляционный коэффициент необратимых потерь давления

) Определим идеальную плотность жидкостной смеси

) Определяем полный градиент давления в точках с заданными давлениями Р<Рнас

) Определяем dH/dP

) Проводим численное интегрирование зависимости dH/dP=f(P)

) Определим удельный объём газа

3. Расчет распределения давления в ЭК при обводненности n=10%. УСШН

Таблица 2

Рис. 2 - График распределения давления в ЭК

Рис. 3 - График распределения плотности ГЖС в ЭК

Рис. 4 - График распределения температуры ГЖС в ЭК

Рис. 5 - График распределения удельного газосодержания жидкости в ЭК

б)=50 % (масс.).

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

Дальнейший расчет см. пункт а).

4. Расчет распределения давления в ЭК при обводненности n=50%. УСШН

Таблица 3

Рис. 6 - График распределения давления в ЭК

Рис. 7 - График распределения плотности ГЖС в ЭК

Рис. 8 - График распределения температуры ГЖС в ЭК

Рис. 9 - График распределения удельного газосодержания в ЭК

УЭЦН

а) =10 % (масс.)

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

) Определим температуру на затрубе скважины

Дальнейший расчет см. пункт «УСШН а)».

5. Расчет распределения давления в ЭК при обводненности n=10% УЭЦН

Таблица 4

Рис. 10 - График распределения давления в ЭК

Рис. 11 - График распределения плотности ГЖС в ЭК

Рис. 12 - График распределения температуры ГЖС в ЭК

Рис. 13 - График распределения удельного газосодержания в ЭК

а) =50 % (масс.)

Где  - средний геотермический градиент скважины

тогда

Дальнейший расчет см. пункт «УСШН а)».

6. Расчет распределения давления в ЭК при обводненности n=50% УЭЦН

Таблица 5

Рис. 14 - График распределения давления в ЭК

Рис. 15 - График распределения плотности ГЖС в ЭК

Рис. 16 - График распределения температуры ГЖС в ЭК

Рис. 17 - График распределения удельного газосодержания в ЭК

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

УСШН

Условия: β=0,85, вынос воды с забоя.

а) =10 %

Рис. 18 - График распределения давления в НКТ и ЭК

.Проверяем выполнение условий выноса (ненакопления) воды в интервале «забой-прием».

Где - приведенная скорость нефти в колонне диаметром , м/с;

- кинематическая вязкость нефти, м2/с.

Условие полного выноса воды из интервала «забой-прием»:случай

Вне зависимости выполнения условия

,

где - необходимая глубина спуска лифта.случай

Численное значение  рассчитывается по зависимости

.

 на глубине

. По графику определим глубину погружения насоса

Давление на приеме насоса

9,56 МПа

Рассчитаем физические свойства продукции на приеме насоса

) Определим количество растворенного в нефти газа

Где С = 0,5 - эмпирический коэффициент

) Определим объемные коэффициенты нефти и жидкости

в) Определим расход жидкости

Где  - дебит дегазированной нефти; - планируемый дебит жидкости

г) Определим расход свободного газа

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

д) Определим расход газожидкостной смеси

3) По графику на рис. определим давление на выходе из насоса

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

Таблица 6

) Определим потери давления в клапанах насоса

Примем, что  и

Определим максимальную скорость движения жидкости в отверстии седла клапана

Для подобранного насоса диаметр клапанов d=20 и 14 мм

в) Определим число Re

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

г) По графику Степановой определяем коэффициент расхода клапана

д) Определим плотность дегазированной жидкости

е) Определим потери давления

) Определим давление в цилиндре насоса при всасывании и нагнетании

) Рассчитаем количество утечек в зазоре плунжерной пары

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

/с

 нет свободного газа

) Рассчитаем коэффициент наполнения

удовлетворяет условию.

б)=50 %

Рис. 19 - График распределения давления в НКТ и ЭК

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

Таблица 7

5) Определим потери давления в клапанах насоса

Примем, что  и

Определим максимальную скорость движения жидкости в отверстии седла клапана

Для подобранного насоса НСВ1-32 диаметр клапанов d=20 мм и d=14 мм

в) Определим число Re

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

г) По графику Степановой определяем коэффициент расхода клапана

д) Определим плотность дегазированной жидкости

е) Определим потери давления

6) Определим давление в цилиндре насоса при всасывании и нагнетании

) Рассчитаем количество утечек в зазоре плунжерной пары

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

/с

) Рассчитаем коэффициент наполнения

удовлетворяет условию.

УЭЦН

Условия: βвх=0,1, вынос воды с забоя.

а) =10 %

Рис. 20 - График распределения давления в НКТ и ЭК

.Проверяем выполнение условий выноса (ненакопления) воды в интервале «забой-прием».

Где - приведенная скорость нефти в колонне диаметром , м/с;

- кинематическая вязкость нефти, м2/с.

Условие полного выноса воды из интервала «забой-прием»:

Вне зависимости выполнения условия , где - необходимая глубина спуска лифта => происходит полный вынос воды с забоя, на любой глубине спуска лифта.

. Так как при любой глубине спуска насоса будет выполняться условие выноса воды, то мы спустим на глубину где =8 МПа.  МПа.

При таких условиях , что удовлетворяет условию =>.

б)=50 %

Рис. 21 - График распределения давления в НКТ и ЭК.

.Проверяем выполнение условий выноса (ненакопления) воды в интервале «забой-прием».

Где - приведенная скорость нефти в колонне диаметром , м/с;

- кинематическая вязкость нефти, м2/с.

Условие полного выноса воды из интервала «забой-прием»:

Вне зависимости выполнения условия

,

где - необходимая глубина спуска лифта => происходит полный вынос воды с забоя, на любой глубине спуска лифта.

. Так как при любой глубине спуска насоса будет выполняться условие выноса воды, то мы спустим на глубину где =8 МПа.

При таких условиях , что удовлетворяет условию =>.

Выбор скважинной насосной установки (УСШН и УЭЦН)

УСШН

а) =10 %

Для данной скважины подобран насос НСВ1-32 и по графику Адонина 5СК-6-1,5-1600

) Рассчитаем коэффициент усадки нефти

) Рассчитаем требуемую подачу насоса

С другой стороны

Где  - площадь поперечного сечения плунжера;  - длина хода полированного штока; N - число двойных ходов;

) Задаемся значением , тогда

) Определим величины упругих деформаций штанг и труб

Где  - длина колонны штанг; - площадь поперечного сечения труб; Е - модуль упругости материала штанг

Где  - доля штанг с площадью поперечного сечения  в общей длине штанговой колонны;

Тогда

) Рассчитаем коэффициент деформации

) Рассчитаем коэффициент утечек

7) Рассчитаем коэффициент подачи

б) =50 %

Для данной скважины подобран насос НСВ1-32 и по графику Адонина 5СК-6-1,5-1600

) Рассчитаем коэффициент усадки нефти

) Рассчитаем требуемую подачу насоса

С другой стороны

Где  - площадь поперечного сечения плунжера;  - длина хода полированного штока; N - число двойных ходов;

) Задаемся значением , тогда

) Определим величины упругих деформаций штанг и труб

Где  - длина колонны штанг; - площадь поперечного сечения труб; Е - модуль упругости материала штанг

Где  - доля штанг с площадью поперечного сечения  в общей длине штанговой колонны;

Тогда

) Рассчитаем коэффициент деформации

) Рассчитаем коэффициент утечек

) Рассчитаем коэффициент подачи

УЭЦН

Необходимый напор ЭЦН можно определить из уравнения основной характеристики:

где  - динамический уровень, из графика распределения = 68 м.

 - напор, затрачиваемый на создание противодавления, на устье;

 - потери напора на трение при движении ГЖС в НКТ;

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

Напор, затрачиваемый на создание противодавления, на устье:

где  - средняя плотность жидкости в НКТ

Потери напора на трение при движении ГЖС в НКТ:

где  - скорость движения жидкости в НКТ, м/с:

Определим число Рейнольдса:

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

Тогда потери напора на трение при движении ГЖС в НКТ:

Положительная работа выделившегося из нефти газа по подъему жидкости, выраженная в метрах столба жидкости.

Тогда требуемый напор насоса:

Теперь по известному дебиту 43,8 м3/сут и по найденному необходимому напору 173,67 м выбираем насос ЭЦН5-80-1400. в=50%:= 939 м, Qж=110,4 м3/сут

Таблица 8 - Характеристики рабочей области насоса

Расчет распределения ,  и  в колонне НКТ (интервал «выкид насоса - устье скважины»);

УСШН

а) =10 % (масс.)

Расчет производим «сверху-вниз».

) Определим величину шага изменения давления

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

) Определим температуру на затрубе скважины

) Определяем температуру потока, соответствующую заданным давлениям

5) Рассчитаем текущее равновесное давление насыщения при :

где  и  - соответственно молярные доли метана и азота в попутном газе

Массовый состав смеси принимаем

) Рассчитываем относительную плотность газа по воздуху:

Молекулярная масса газа:

Вычисляем молекулярную массу данного газа:

) Определим объем выделившегося газа

) Определим остаточную газонасыщенность нефти в процессе разгазирования (Р=1,5 МПа; Т=294 К)

) Определим относительную плотность выделившегося газа

) Определим относительную плотность растворенного газа, остающегося в нефти при данных условиях ее разгазирования

) Определим удельное приращение объема нефти за счет единичного изменения ее газонасыщенности  и объемный коэффициент нефти

) Определим коэффициент сверхсжимаемости газа

а) Определим Рпр и Тпр

,

Таблица 9 - Данные о критических давлениях и температуре различных компонентов

б) Определим коэффициент сверхсжимаемости для азота

в) Определим коэффициент сверхсжимаемости углеводородной составляющей газа т.к Рпр находится в интервале 0≤Рпр≤3,8 ,а Тпр - в интервале 1,17≤Тпр≤2,0, то

г) Определим коэффициент сверхсжимаемости газовой смеси

) Определяем удельный объем газожидкостной смеси

) Определим удельную массу смеси при стандартных условиях

) Определяем корреляционный коэффициент необратимых потерь давления

) Определим идеальную плотность жидкостной смеси

) Определяем полный градиент давления в точках с заданными давлениями Р<Рнас

) Определяем dH/dP

) Проводим численное интегрирование зависимости dH/dP=f(P)

) Определим удельный объём газа

7. Расчет распределения давления в НКТ при обводненности n=10% УСШН

Таблица 10

Рис. 22 - График распределения давления в НКТ

Рис. 23 - График распределения плотности ГЖС в НКТ

Рис. 24 - График распределения температуры ГЖС в НКТ

Рис. 25 - График распределения удельного газосодержания жидкости в НКТ

б)=50 % (масс.).

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

Дальнейший расчет см. пункт а).

8. Расчет распределения давления в НКТ

Таблица 11

Рис. 26 - График распределения давления в НКТ

Рис. 27 - График распределения плотности ГЖС в НКТ

Рис. 28 - График распределения температуры ГЖС в НКТ

Рис. 29 - График распределения удельного газосодержания в НКТ

УЭЦН

а) =10 % (масс.)

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

) Определим температуру на затрубе скважины

Дальнейший расчет см. пункт «УСШН а)».

9. Расчет распределения давления в ЭК

Таблица 12

Рис. 30 - График распределения давления в НКТ

Рис. 31 - График распределения плотности ГЖС в НКТ

Рис. 32 - График распределения температуры ГЖС в НКТ

Рис. 33 - График распределения удельного газосодержания в НКТ

а) =50 % (масс.)

) Определим температурный градиент потока

Где  - средний геотермический градиент скважины

тогда

) Определим температуру на затрубе скважины

Дальнейший расчет см. пункт «УСШН а)».

. Расчет распределения давления в НКТ

Таблица 13

Рис. 34 - График распределения давления в НКТ

Рис. 35 - График распределения плотности ГЖС в НКТ

Рис. 36 - График распределения температуры ГЖС в НКТ

с

Рис. 37 - График распределения удельного газосодержания в НКТ

Определение глубины начала интенсивного образования АСПО;

УСШН

а) =10 %

Рис. 38 - График распределения для Тскв и Таспо

Глубина интенсивного образования АСПО равна 1290м.

б) =50 %

Рис. 39 - График распределения для Тскв и Таспо

Глубина интенсивного образования АСПО равна 1200м.

УЭЦН

а) =10 %

Рис. 40 - График распределения для Тскв и Таспо

Глубина интенсивного образования АСПО равна 1330м.

б) =50 %

Рис. 41 - График распределения для Тскв и Таспо

пласт скважина газ насос

Глубина интенсивного образования АСПО равна 1180м.

Выбор и обоснование технологий предупреждения образования АСПО в колонне НКТ и удаления отложений;

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

Механический способ борьбы с отложениями АСПО - это способ очистки НКТ от парафина при помощи скребков. На скважинах эксплуатирующихся фонтанным способом или ЭЦН применяются установки УДС-1М, ПАДУ, ПАДУ-3 и т.д. Борьбу с отложениями парафина в затрубном пространстве проводят механическим путем с помощью скребков, закрепленных к штангам. Эти скребки автоматически поворачиваются на определенный угол при каждом ходе колонны штанг вниз.

Физический способ.

В настоящее время для предупреждения отложений АСПО на подземном оборудовании скважин месторождений «ЛУКОЙЛ - ПЕРМЬ» успешно применяются магнитные аппараты марки МАС-2,5. В настоящее время, с технической точки зрения, установка магнитных аппаратов возможна на фонтанных, оборудованных установками ШГН и ШВН скважинах. Перед применением магнитных аппаратов рекомендуется провести лабораторные испытания скважинной жидкости на магнитную восприимчивость.

В случае интенсивного отложения АСПО в колонне НКТ возможно более рентабельным будет применять нагревательных линий марки НКЛ-16-400-8-800 по ТУ 3542-001-001-48369-2003 (ООО «ПермНИПИнефть»).

Из химических методов предотвращения отложений асфальто-смоло-парафиновых веществ (АСПВ) применяют жидкие и твердые ингибиторы. В случае применения жидких ингибиторов, дозировка ингибиторов осуществляется с помощью поверхностных и глубинных дозаторов, выпускаемых отечественной промышленностью. Кроме того, на месторождениях применяются твердые ингибиторы серии ИКД, разработанные в ООО «ПермНИПИнефть», в которых происходит самодозирование за счет постепенного растворения ингибитора в скважинной жидкости. Установка контейнера возможна при всех способах эксплуатации скважины, действие контейнера рассчитано на 1 год. Твердые ингибиторы серии ИКД помещаются в добывающие скважины в перфорированных контейнерах, изготовленных из 4 или 8 насосно-компрессорных труб (НКТ) длиной 5 или 2,5 метра, соответственно, каждая камера контейнера отделена друг от друга перфорированной шайбой. Верхний и нижний концы контейнера заглушены такими же перфорированными шайбами. Твердый ингибитор ИКД помещается в камеры контейнера. После спуска глубинного оборудования и запуска скважины в работу, добываемые флюиды омывают реагент через перфорированные отверстия.

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

Список литературы

1.      Сборник задач по технологии и технике нефтедобычи. Мищенко И. Т.

.        Расчеты в добыче нефти. Мищенко И. Т.

.        Технология и техника добычи нефти. Щуров В. И.

.        Скважинная добыча нефти. И. Т. Мищенко.