Особенности измерения абсолютной проницаемости горных пород
- Авторы: Песков А.В.1
-
Учреждения:
- Самарский государственный технический университет
- Выпуск: Том 28, № 2 (2020)
- Страницы: 73-83
- Раздел: Приборостроение, метрология и информационно-измерительные приборы и системы
- URL: https://journals.eco-vector.com/1991-8542/article/view/41714
- DOI: https://doi.org/10.14498/tech.2020.2.5
- ID: 41714
Цитировать
Полный текст
Аннотация
Рассмотрены аспекты измерения газопроницаемости на образцах терригенных и карбонатных пород коллекторов нефти и газа, а также искусственных образцах на отечественной установке «Дарсиметр». Для исследования использовались следующие типы пород коллекторов: поровые, трещиноватые, кавернозные. Уточнялась область применения закона Дарси для пород порового типа в интервале абсолютной проницаемости от 0.015 до 100 мД. Определялись коэффициенты проницаемости с учетом эффекта проскальзывания газа по уравнениям регрессии. Определялась кажущаяся проницаемость при малых перепадах давления. Для ряда образцов с низкой проницаемостью рассчитывался размер пор для связи со значением кажущей проницаемости. Расчет производился по полученным величинам структурных коэффициентов пород методом электрического удельного сопротивления и по значениям пористости, определенным с использованием метода Преображенского.
Для ряда образцов с естественной и искусственными трещинами и капилляром зависимость между градиентом давления и скоростью фильтрации становилась нелинейной и были определены виды законов фильтрации. Установление линейного или нелинейного закона контролировалось построением индикаторных кривых и расчетом числа Рейнольдса. Для терригенных пород высокой проницаемости определялись погрешности измерения коэффициентов газопроницаемости в различных интервалах перепада давления: подсчитывались дисперсия и коэффициент вариации, показавшие низкие значения.
Полный текст
Изучение абсолютной проницаемости горных пород на образцах пород является основным и начальным этапом для получения фильтрационно-емкостной характеристики пластов-коллекторов и покрышек [11]. Эти данные применяются при подсчете запасов месторождений и эксплуатации подземных хранилищ газа. Абсолютная проницаемость определяет качество флюидоупоров и дебиты скважин. Коэффициент абсолютной проницаемости определяется в петрофизических лабораториях, как правило, согласно ГОСТ 26450.2-85 «Породы горные. Метод определения коэффициента абсолютной проницаемости при стационарной и нестационарной фильтрации» или по сведениям, изложенным в отечественной и зарубежной литературе по петрофизическим исследованиям [2, 5, 6]. Проницаемость пород изменяется широко от 10-18 до 10-11 м2. Фильтрация флюида может осуществляться через поры, естественные и искусственно созданные трещины, каверны. Структура пустотного пространства различна у терригенных и карбонатных пород. У терригенных пород диаметры пор и соединяющие их каналы изменяются ненамного. В карбонатных коллекторах диаметры поровых каналов на 12 порядка меньше диаметра микрокаверн, составляющих основную емкость коллекторов. Сложный тип карбонатных пород коллекторов включает в пустотное пространство поры, каверны и трещины. Ширина эффективных трещин варьирует от 10 до 20 мкм, а длина от 2 мм до максимальных значений [16]. Трещины могут быть извилистыми и прямолинейными.
Абсолютная газопроницаемость определяется по эмпирическому закону Дарси.
Закон Дарси для определения газопроницаемости определяется выражением
(1)
где
Q − обьемный расход газа через образец;
PА, PВ − давление на входе и выходе образца;
F − площадь поперечного сечения образца;
k − проницаемость;
μ − вязкость газа;
l − длина образца.
Впервые исследованиями Л. Клинкенберга [1] было обнаружено, что проницаемость пород по газу kг может быть выше проницаемости по жидкости kж, что объясняется эффектом проскальзывания газа в условиях, когда длина среднего пробега молекул газа соизмерима с диаметром капилляра.
Значение проницаемости зависит от среднего внутрипорового давления и определяется выражением [9]
(2)
где
b − константа Клинкенберга (коэффициент скольжения);
pср − среднее давление при определении газопроницаемости.
Закон Дарси с учетом эффекта Клинкенберга определяется зависимостью
(3)
Приближенно константу Клинкенберга можно определить согласно уравнению [3]
b = 0,77 kж-0,39. (4)
Эффект проскальзывания газа (воздуха) тем больше, чем меньше коэффициент проницаемости горных пород и перепад давления, что было уточнено в работе [5].
Согласно [14] эффект Клинкенберга может проявляться только в лабораторных экспериментах при фильтрации разреженных газов в образцах керна низкой газопроницаемости меньше 1 мД; в пластовых условиях месторождений (при давлениях больше 11,5 МПа) эффект Клинкенберга ничтожно мал.
Закон Дарси при фильтрации газа через трещину определяется из зависимости [8]
(5)
где −проницаемость трещин.
Различают также изменение проницаемости трещины, связанное с изменением неровностей, шероховатостей берегов трещин введением коэффициента относительной шероховатости [10].
В ряде случаев линейная зависимость между скоростью течения флюида и давлением нарушается за счет проявления инерционных сопротивлений при увеличении скорости фильтрации. Связь между расходом флюида и градиентом давления при нелинейном законе фильтрации можно установить в виде [4, 6, 12]
(6)
где а, b − коэффициенты, учитывающие влияние вязкости и плотности флюида.
Границы применимости линейного закона фильтрации можно определить исходя из числа Re. В.Н. Щелкачев [7] получил выражение для числа Re:
(7)
где
V − cкорость фильтрации;
k − проницаемость;
кинематическая вязкость флюида;
m − пористость.
Критические значения числа Re в этой зависимости 112. При превышении этого значения определение фильтрационных параметров возможно проводить по зависимости (6).
Cоответствие фильтрации линейному закону согласно [6] можно контролировать по параметру Дарси:
где
V − cкорость фильтрации;
k − проницаемость;
динамическая вязкость флюида;
Δp − перепад давления;
L − длина образца.
Исходя из выражения (6) получают коэффициенты, характеризующие пропускную способность пористой среды [4]: Kμ вязкостной коэффициент, характеризующий пропускную способность породы для вязкого флюида, который измеряется в дарси; Kρ плотностной коэффициент, характеризующий пропускную способность породы для невязкого флюида. Плотностной коэффициент имеет размерность длины:
(9)
Чем больше площадь сечения поровых каналов, тем выше значение Kμ. C уменьшением извилистости поровых каналов возрастает Kρ. Еще один способ контроля линейности потока флюида с помощью графика Q=f(Δp). При линейном характере фильтрации расход флюида Q прямо пропорционален перепаду давления, под действием которого происходит фильтрация, и все точки лежат на прямой, выходящей из начала координат. Точки, выпадающие из этой зависимости, указывают на отклонение от линейного закона.
Целью данных экспериментальных исследований явилось изучение газопроницаемости пород в широком интервале абсолютной проницаемости от 0,15 до 100 мД при наличии в породах трех типов пустот: пор, каверн, трещин. Определялось соответствие фильтрации линейному закону и погрешность измерения проницаемости. Исследования проводились как на природном керне, так и на искусственных образцах, содержащих капилляры и трещины. Исследования проводились на отечественном приборе «Дарсиметр» в лаборатории физики пласта СамГТУ. Определялись погрешности определения газопроницаемости для нескольких интервалов перепадов давления на образцах. Для серии образцов рассчитывались величины kж по графику Клинкенберга, отношение максимального значения проницаемости Kмин, измеренного при минимальном давлении Pмин, к проницаемости по жидкости kмак/kж, константа Клинкенберга bр по зависимости (4) и константа Клинкенберга по приближенной формуле bпр (табл. 1). Стоит уточнить, что kж это истинная газопроницаемость, рассчитанная по графику Клинкенберга, а kмак кажущаяся повышенная проницаемость, измеренная при минимальном перепаде давления. Исходя из полученных данных по табл. 1 максимальным значением kмак/kж характеризуются образцы с минимальной абсолютной проницаемостью (образцы 59, 69, 63, 10), bр константа Клинкенберга была максимальной, хотя для образцов 85, 9 (содержащего каверны), искусственной щели, капилляра значения bр имеют минимальные значения при относительно средних значениях абсолютной проницаемости в рассматриваемой выборке. Для капилляра 1 и искусственной щели не наблюдается отклонения в законе Дарси в связи с незначительным значением сечения пустот. Для образца естественной трещины 3, капилляра 2, искусственных трещин значения bр не рассчитывались в связи с тем, что фильтрация газа через образцы осуществлялась по нелинейному закону.
Рассмотрим измерения параметров фильтрации через образцы 31, 19. Согласно табл. 2 число Re не превышает 1,1, что соответствует линейному закону фильтрации. Это также подтверждает индикаторная диаграмма, являющаяся линейной при всех градиентах давления. Подсчитанные средние выборочные проницаемости при близких давлениях характеризуются незначительными дисперсиями D =0.273 и коэффициентом вариации не более v = 0.85 %, что характеризует оптимальные условия регистрации параметров фильтрации и качество прибора.
Таблица 1 Результаты измерения проницаемости kж, bр, bпр изученных образцов
№ | Описание | pмин | pмак | Kмин | Kмак | Кж | Kмак/kж | bр | bпр |
– | Эталон | 0.286 | 1.6 | 1.275 | 1.5 | 0.806 | 1.86 | 0.98 | 0.65 |
69 | Известняк | 0.367 | 1.55 | 0.126 | 0.17 | 0.057 | 3.05 | 2.4 | 1.53 |
59 | Известняк | 0.24 | 1.47 | 0.028 | 0.04 | 0.006 | 7.1 | 7.3 | 3.35 |
85 | Известняк | 0.32 | 1 | 0.2 | 0.2 | 0.19 | 1.05 | 0.06 | 1.1 |
76 | Известняк | 0.43 | 1.54 | 0.49 | 0.54 | 0.42 | 1.27 | 0.31 | 1.02 |
60 | Известняк | 0.43 | 1.7 | 0.048 | 0.058 | 0.32 | 1.82 | 0.97 | 2.46 |
63 | Известняк | 0.38 | 1.7 | 0.047 | 0.067 | 0.016 | 4.07 | 3.67 | 2.59 |
10 | Песчаник | 0.35 | 1.6 | 0.13 | 0.18 | 0.035 | 5.2 | 4.74 | 1.5 |
49 | Известняк каверны | 0.27 | 1.47 | 3.84 | 4.31 | 1.16 | 3.72 | 2.85 | 0.43 |
56 | Известняк каверны | 0.285 | 1.74 | 12.08 | 13.4 | 11.1 | 1.2 | 0.23 | 0.28 |
15 | Искусственная микрощель | 0.26 | 1.68 | 11.49 | 12.09 | 10.7 | 1.11 | 0.07 | 0.308 |
11 | Капилляр 1 | 0.407 | 1.86 | 2.64 | 2.7 | 2.63 | 1.08 | 0.11 | 0.52 |
31 | Песчаник | 0.32 | 1.88 | 51.2 | 61.1 | 28.68 | 2.1 | 0.17 | 0.17 |
19 | Песчаник | 0.28 | 1.82 | 101 | 90 | 78.3 | 1.36 | 0.39 | 0.13 |
1 | Песчаник | 0.39 | 1.7 | 9.9 | 8.54 | 6.7 | 1.48 | 0.5 | 0.33 |
3 | Трещина известняк | 0.26 | 1.64 | 229 | 111 | – | – | – | – |
Погрешность измерения снижается при увеличении Δp от 0.25 до 1.3 атм и максимальна при Δp = 0.25 атм. Cогласно графику Клинкенберга kж = 28.68 мД, а Kмак/kж = 2.13 для образца 31. Для образца 19 Kмак/kж = 1.36. Значения Kмак/kж для рассматриваемых образцов 19, 31 имеют высокие величины, которые были ожидаемы так же, как для малопроницаемых образцов, например 76 и 1. Для образца 19 индикаторная диаграмма и число Re также соответствуют линейному закону фильтрации. Параметр Дарси Da оказался малоэффективным для определения характера фильтрации. Согласно графикам Клинкенберга истинная проницаемость для образца 19 составила 76.67 мД. Для изученных образцов песчаника 1931, содержащих поры без трещин и капилляров с абсолютной проницаемостью 50100 мД, с увеличением перепада давлений Δp до 1,65 атм расход Q, cм3/с возрастает до 38,50 см3/c и проявляется эффект Клинкенберга.
Таблица 2 Результаты измерения проницаемости k при различных давлениях Δp образца 31
№ | Q, cм3/c | P, атм | K, атм | Da | Re | К СР, мД | D, мД2 | V, % |
1 | 6.18 | 0.59 | 61.1 | 1.01 | 0.48 | 61.0
| 0.273
| 0.85
|
2 | 6.09 | 0.59 | 60.2 | 1.01 | 0.48 | |||
3 | 6.23 | 0.59 | 61.6 | 1.01 | 0.49 | |||
4 | 6.17 | 0.59 | 61.3 | 1.01 | 0.49 | |||
5 | 6.73 | 0.64 | 60.4 | 1.01 | 0.52 | |||
6 | 11.46 | 1.02 | 56.2 | 1.01 | 0.77 | 55.7
| 0.055
| 0.420
|
7 | 11.67 | 1.04 | 55.8 | 1.01 | 0.75 | |||
8 | 10.82 | 0.99 | 55.6 | 1.01 | 0.71 | |||
9 | 10.78 | 0.98 | 55.5 | 1.01 | 0.70 | |||
10 | 10.81 | 0.98 | 55.7 | 1.01 | 0.70 | |||
11 | 11.62 | 1.04 | 55.8 | 1.01 | 0.75 | |||
12 | 10.76 | 0.98 | 55.4 | 1.01 | 0.70 | |||
13 | 10.76 | 0.98 | 55.7 | 1.01 | 0.70 | |||
14 | 11.54 | 1.03 | 55.9 | 1.01 | 0.74 | |||
15 | 15.60 | 1.32 | 53.6 | 1.01 | 0.92 | 53.6
| 0.045
| 0.395
|
16 | 15.03 | 1.28 | 53.9 | 1.01 | 0.87 | |||
17 | 14.87 | 1.28 | 53.3 | 1.01 | 0.87 | |||
18 | 15.80 | 1.34 | 53.6 | 1.01 | 0.90 | |||
19 | 14.93 | 1.28 | 53.8 | 1.01 | 0.87 | |||
20 | 22.05 | 1.74 | 51.4 | 1.01 | 1.13 | 51.3 | 0.045 | 0.41 |
21 | 21.27 | 1.70 | 51.1 | 1.01 | 1.08 | |||
22 | 22.40 | 1.75 | 51.6 | 1.01 | 1.11 | |||
23 | 21.22 | 1.70 | 51.0 | 1.01 | 1.08 | |||
24 | 21.30 | 1.70 | 51.2 | 1.01 | 1.07 |
Рассмотрим фильтрацию через образец песчаника, содержащего трещину. Индикаторная диаграмма изгибается, и это свидетельствует о нелинейном законе фильтрации. Рассчитанные по формуле Дарси коэффициенты проницаемости K* искажены за счет нелинейной фильтрации и не являются истинными (табл. 3). Далее был получен график PQ/l от Q и рассчитаны коэффициенты а и b в уравнении (6).
Таблица 3 Результаты измерения проницаемости k при различных давлениях Δp образца с трещиной
Q, cм3/с | ΔP, атм | K*, мД |
8.9 | 0.259 | 229.9 |
17.6 | 0.6 | 170.5 |
27.7 | 0.97 | 145.3 |
40.4 | 1.484 | 118.1 |
44 | 1.644 | 111.0 |
Нелинейная зависимость для образца с трещиной имеет вид , если давление измеряется в паскалях, расход в м3/с. При этом Kμ = 510-13, Kϱ = 3.7110-10. Параметры трещины определялись по периметру, выходящему на торец образца; раскрытость составила b = 0.126 мм. Пористость трещины m = 4.410-4.
При этом проницаемость K, определенная по зависимости K = 85000b2m [8], cоставила 44710-15м2.
Рассмотрим фильтрацию через искусственный образец, содержащий капилляр 2 с диаметром 200 мкм (табл. 4). Индикаторная диаграмма показывает нелинейный характер фильтрации. Коэффициенты проницаемости через образец подсчитывались по закону Дарси с учетом зависимости проницаемости [8], где m пористость образца с капилляром; R − радиус капилляра; − cтруктурный коэффициент. В связи с нелинейной зависимостью фильтрации коэффициенты проницаемости, рассчитанные по закону Дарси в табл. 4, не являются истинными. Для образца, содержащего капилляр, при перепаде давлений Δp до 1,65 атм расход через образец снижается почти в 5 раз относительно образцов песчаника 1931, однако индикаторная кривая указывает на нарушение закона Дарси и фильтрация идет по нелинейному закону.
Нелинейная зависимость для образца с капилляром в диапазоне давлений до 1 атм имеет вид , если давление измеряется в паскалях, расход в м3/с. При этом Kμ = 127*10-15, Kρ = 8.84*10-12.
До давления 1.76 атм зависимость от Q аппроксимируется многочленом пятой степени при величине достоверности аппроксимации R2 = 0.996:
.
Таблица 4 Результаты измерения проницаемости k при различных давлениях Δp образца с капилляром № 2
№ | Q, cм3/c | P, атм | K, мД |
1 | 7.14 | 1.76 | 16.3 |
2 | 7.03 | 1.76 | 16.1 |
3 | 6.05 | 1.35 | 20.3 |
4 | 6.02 | 1.34 | 20.3 |
5 | 4.79 | 0.96 | 25.4 |
6 | 4.78 | 0.96 | 25.5 |
7 | 4.77 | 0.96 | 25.6 |
8 | 2.90 | 0.51 | 34.6 |
9 | 3.09 | 0.54 | 34.1 |
10 | 3.09 | 0.55 | 33.4 |
11 | 1.63 | 0.24 | 45.8 |
12 | 1.74 | 0.26 | 44.8 |
Проницаемость, рассчитанная по [8] по зависимости при структурном коэффициенте = 1, cоставила 79.510-15 м2. Это значение сопоставимо с найденным значением из графика Kμ = 127·10-15.
Рассмотрим фильтрацию через искусственный образец трещины раскрытостью 30 мкм. Это значение характерно для интервала средних величин раскрытости трещиноватых сложных карбонатных пород. Рассчитанная трещинная пористость составила m = 0,00064, рассчитанный коэффициент проницаемости k = 48мД, рассчитанный по закону Дарси коэффициент проницаемости k = 49 мД. Получена линейная зависимость расхода газа q от перепада давления Δp. Полученные расчетным и экспериментальным способами близкие значения коэффициента проницаемости указывают на достоверность получаемых результатов на искусственных моделях трещин, моделирующих природные трещины, и открывает возможности для определения проницаемости на искусственно созданных трещинах для гидроразрыва пласта одного из эффективных методов интенсификации добычи нефти.
Искуственный образец трещины с просветностью b = 47 мкм имеет расчетную пористость m = 0.00106, расчетную проницаемость k = 225мД. Экспериментальное значение проницаемости, полученное на установке «Дарсиметр» при Δp = 1.72 атм, составило k = 211мД. Получена нелинейная зависимость расхода газа q от перепада давления Δp. Полученные расчетным и экспериментальным способами значения коэффициента проницаемости для Δp = 1.72 близки. Для Δp = 0.4 атм измеренное значение на установке «Дарсиметр» cоставило k = 281 мД. Таким образом, нелинейный закон фильтрации газа осуществляется при раскрытости щели более 47 мкм.
По образцу с кавернами число Re ниже 1 и закон фильтрации является линейным, прослеживается увеличение проницаемости с уменьшением перепада давления, характерное для эффекта проскальзывния газа и Kмак/kж = 1.2. Параметр Дарси Da не отразил изменения характера фильтрации и при данных небольших параметрах расхода оказался малоэффективным.
Таблица 5 Результаты измерения и расчета удельного сопротивления пород ρп, параметра пористости Pп , структурного коэффициента ϕ, пористости m, проницаемости k, радиуса поровых каналов R
№ образца | ρп, ом·м | Pп | m, % | ϕ | ϕрасч | K, мд | R, мкм |
59 | 67.7 | 451.0 | 0.01 | 4.06 | 89.61 | 0.006 | 1.89 |
85 | 44.3 | 295.5 | 0.05 | 14.66 | 13.71 | 0.190 | 20.19 |
66 | 49.0 | 326.6 | 0.02 | 6.63 | 36.62 | 0.070 | 16.33 |
61 | 49.0 | 326.6 | 0.03 | 9.70 | 24.10 | 0.088 | 20.87 |
60 | 98.0 | 653.2 | 0.03 | 18.62 | 25.22 | 0.320 | 37.71 |
63 | 88.7 | 591.0 | 0.03 | 15.96 | 26.76 | 0.016 | 8.37 |
76 | 42.0 | 280.0 | 0.01 | 3.64 | 59.79 | 0.420 | 32.82 |
69 | 21.0 | 140.0 | 0.05 | 7.07 | 13.44 | 0.057 | 11.09 |
В ходе исследований пород методом удельного электрического сопротивления [9, 15] (табл. 5) определялись удельные сопротивления пород, насыщенных пластовой водой с минерализацией 50 г/л, рассчитывался параметр пористости, определялась пористость пород методом Преображенского, и по этим данным рассчитывался структурный коэффициент, а далее по зависимости И.Ф. Котяхова [13] размер пор.
Отклонение от закона Дарси, связанное с эффектом проскальзывания газа по [4], наблюдается для размеров пор менее 10 мкм, поэтому для ряда малопроницаемых образцов определялись средние размеры пор по данным пористости с учетом структурного коэффициента, измеренного по данным метода удельного электрического сопротивления (см. табл. 5). Для образцов 59, 69, 63 c размером пор менее 11 мкм характерны значения Kмак/kж = 3.57. Вероятно, в образцах 61, 66, 31 присутствует большая доля пор с размером менее 10 мкм.
Заключение
Измерения на аппаратуре «Дарсиметр» показали низкие погрешности при определении газопроницаемости не более (коэффициент вариации) V = 0,45 %. Измерения показали для ряда образцов возрастание проницаемости в несколько раз в области малых поровых давлений. Эффект проскальзывания газа для образцов, содержащих поры, тем больше, чем меньше коэффициент абсолютной проницаемости пород и перепад давления в образце. Применимость закона Дарси определяется размером поровых каналов. При среднем размере пор до 10 мкм наблюдается эффект проскальзывания Клинкенберга. Методом удельного электрического сопротивления определены удельные сопротивления образцов пород, рассчитан параметр пористости пород, пористость пород методом Преображенского, структурный коэффициент и размер пор. Выделены образцы, для которых размер пор менее 11 мкм имеет значительные величины Kмак/kж = 3.57. Для образцов, представленных трещинами или капиллярами, характерно отсутствие извилистости и значительные сечения пустотного пространства. Фильтрация по таким пустотам происходит по нелинейному закону. Таким образом, изучены образцы терригенных и карбонатных пород и искусственные образцы, содержащие капилляры и трещины. Определен вид зависимости для нелинейного закона фильтрации через капилляр, естественную тещину и ряд искусственных трещин. Для трещиноватых горных пород отклонение газопроницаемости от истинных значений может оказаться существенным, что в итоге сказывается на общей оценке фильтрационно-емкостных свойств пластовых систем. Измерения проницаемости при малых перепадах давления показали на исследованных образцах приращение значений проницаемости, даже если абсолютная проницаемость имеет значения 50100 мД.
Об авторах
Александр Владимирович Песков
Самарский государственный технический университет
Автор, ответственный за переписку.
Email: journal@eco-vector.com
(к.г.-м.н.), доцент кафедры «Разработка и эксплуатация нефтяных и газовых месторождений»
Россия, 443100, г. Самара, ул. Молодогвардейская, 244Список литературы
- Klinkenberg L.J. The permeability of porous media to liquids and gases // Am. Petroleum Inst. Drilling and Production Practice, 1941. Pp. 200–211.
- Методическое руководство по исследованию пород коллекторов нефти и газа физическими и петрографическими методами для производственных лабораторий. М.: ВНИГНИ, 1973.
- Виноградов В.Г., Дахнов А.В., Пацевич С.Л. Практикум по петрофизике. М.: Недра, 1990. 227 с.
- Ханин А.А. Петрофизика нефтяных и газовых пластов. М.: Недра, 1976. 295 с.
- Песков А.В., Ольховская В.А. Определение газопроницаемости пород с учетом эффекта скольжения газа // Нефтепромысловое дело. 2010. № 3. С. 10–12.
- Гудок Н.C., Богданович Н.Н., Мартынов В.Г. Определение физических свойств нефтеводосодержащих пород. M.: Недра-Бизнесцентр, 2007. 563 c.
- Щелкачев В.Н., Лапук Б.Б. Подземная гидравлика. М.: НИЦ, 2001. 735 c.
- Гиматудинов Ш.К., Ширковский А.И. Физика нефтяного и газового пласта. М.: Недра, 1982. 309 c.
- Добрынин В.М., Вендельштейн Б.Ю., Кожевников Д.А. Петрофизика. М.: Недра, 1991. 367 c.
- Михайлов Н.Н. Физика нефтяного и газового пласта. М.: Макспресс, 2008. 445 c.
- Орлов Л.И., Карпов Е.Н., Топорков В.Г. Петрофизические исследования коллекторов нефти и газа. М.: Недра, 1987. 214 c.
- Мирзаджанзаде А.Х., Аметов И.М., Ковалев А.Г. Физика нефтяного и газового пласта. М.: Недра, 1992. 268 c.
- Котяхов Ф.И. Физика нефтяных и газовых коллекторов. М.: Недра, 1972. 286 c.
- Троицкий В.М. Феноменологический подход к анализу экспериментальных данных о газопроницаемости в пористых средах. Истинная причина эффекта Клинкенберга // Вести газовой науки. 2017. № 2 (30). С. 111–121.
- Тульбович Б.И. Методы изучения пород коллекторов нефти и газа. М.: Недра, 1979. 368 c.
- Гмидт Л.П., Леви С.Ш. Атлас карбонатных пород коллекторов. М.: Недра, 1972. 81 c.