Permeability prediction belongs to the secondary interpretation which means that it is based on results of preliminary open-hole logs interpretation in terms of volumetric rock model (usually, Effective porosity and Shaliness) and normally does not involve direct inputs from open-hole logs.
Some open-hole logs (like variation between near and far zone resistivity) respond to permeability but can not be used in quantitative analysis because the value of response on other rock properties which are difficult to assess and account for.
Inputs & Outputs
Inputs |
| Outputs | ||
---|---|---|---|---|
\phi_e(l) | Effective porosity |
| Absolute permeability | |
V_{sh}(l) | Shaliness |
where l(x,y,z) is the along-hole depth and (x,y,z) – 3D spatial coordinates.
The model parameters are being fit to absolute permeability core data for each lithofacies separately LFI(l).
Models
Exponential model
Popular empirical model:
(1) | k = k_0 \, \exp \big[ \beta_{\phi} \, \phi \big] |
Power model
Popular physical model:
(2) | k = k_0 \, \phi^{m} |
In the early times this model was though to be empirical but later on it was shown that it has a clear physical meaning.
Dual-component power model
Если общий объем низкопоровых пропластков в пределах залежи велик и суммарные запасы внутри таких пропластков значительны, то пренебрегать проницаемостью в таких пропластках нельзя и необходимо строить двухкомпонентную модель, учитывающая различие в модели проницаемости для низкопроницаемых и высокопроницаемых пропластках.
Популярной моделью на этом пути является модель двухстепенной проницаемости:
(3) | k = k_1 \, \phi^{m_1} + k_2 \, \phi^{m_2} |
в котором компонента k_1 \, \phi^{m_1} определяет проницаемость для низкопроницаемых пропластков, а компонента k_2 \, \phi^{m_2} определяет проницаемость для высокопроницаемых пропластков.
Эта модель имеет физическое обоснование в рамках классической модели гидродинамического потока в поровом коллекторе с фрактальной структурой.
Kozeny-Carman model
Модель Козени-Кармана является исторически первой и наиболее популярной физической моделью проницаемости, основанной на гидродинамических уравнениях флюида в поровом коллекторе с упрощенной структурой независимых друг от друга капилляров разной длины, диаметра и формы и представляет собой линейную зависимость индекса качества пород от нормализованной пористости \phi_r:
(4) | RQI = FZI \phi_r |
где множитель FZI называется Динамическим Индексом Пород (Flow Zone Index) или сокращенно ДИП, пропорционален размеру породообразующих зерен D_g:
(5) | FZI \sim D_g |
и предполагается постоянным для каждой литофации
Явное выражение проницаемости в модели Кармана-Козени дается следующей формулой:
(6) | k = 1014.24 \, FZI^2 \, \frac{\phi^3}{( 1 - \phi )^2} |
General approach to permeability modelling
В отличие от простой теоретической модели Кармана-Козени в большинстве практических случаев динамический индекс пород не является постоянной величиной и варьируется внутри заданной литофации, в зависимости от конкретного значения пористости и глинистости пласта:
(7) | FZI = FZI(V_{sh}, \phi) |
которая учитывает вариабельность пористости и степени заглинизированности коллектора внутри одной литофации.
На практике динамический индекс пород удобнее коррелировать с относительной (нормализованной) пористостью \phi_r:
(8) | FZI = FZI(V_{sh}, \phi_r) |
После того, как зависимость FZI(\phi) построена, рассчитывается индекс качества коллектора:
(9) | RQI = FZI \, \phi_r |
а на его основе и проницаемость
(10) | k = 1014.24 \, \phi \, RQI^2 = 1014.24 \, \phi \, FZI(V_{sh}, \, \phi)^2 \, \phi_r^2 = 1014.24 \, FZI(V_{sh}, \, \phi)^2 \, \frac{\phi^3}{( 1 - \phi)^2} |
В частном случае, когда
FZI(\phi) = \rm const является постоянной величиной внутри каждой литофации, модель проницаемости
(10) сводится к модели Козени-Кармана
(6).
Dual-component Kozeny-Carman model
Для более тонкой настройки проницаемости можно использовать обобщенную двух-компонентную модель Кармана-Козени
(11) | FZI(V_{sh},\phi_r) = FZI_1(V_{sh}) \, \phi_r^{m_1} + FZI_2(V_{sh}) \, \phi_r^{m_2} |
где
(12) | FZI_1(V_{sh}) = FZI_{01} \, (1- V_{sh}/V_{sh1})^{g_1} |
(13) | FZI_2(V_{sh}) = FZI_{02} \, (1- V_{sh}/V_{sh2})^{g_2} |
Константы \{ FZI_{01}, \, FZI_{02} \} описывают предельно высокие значения динамического индекса заданной литофации для пропластков с низкой заглинизированностью.
Константы \{ V_{sh1}, \, V_{sh2} \} описывают критические значения глинистости при которых соответствующая компонента проницаемости исчезает.
Константы \{ g_1, \, g_2 \} описывают степень сцементированности глин и при низких значениях приводят к ослаблению зависимости динамического индекса пласта от коэффициента глинистости.
Эта модель описывает большое количество практических случаев в широких пределах изменения пористости и глинистости.
Главное отличие формулы (11) от (3) является то, что в координатах \{ FZI, \, \phi_r \} керновые данные, как правило, имеют гораздо лучшую кучность, чем в координатах \{ k, \, \phi \} и формула (11) описывает гораздо более слабую зависимость от пористости (то есть числа \{ m_1, \, m_2 \} очень малы), что позволяет строить более точные модели проницаемости.
Artificial Neural Network permeability model
Нейросетевая модель проницаемости основывается на универсальной нейросетевой регрессии динамического индекса пород FZI на глинистость V_{sh} и нормализованную пористость \phi_r:
(14) | \{ V_{sh}, \phi_r \} \rightarrow {\rm Neural Network} \rightarrow FZI |
Она в состоянии описать более сложные виды зависимости, чем двух-компонентная модель Кармана-Козени (11).
Однако, как и все виды виды универсальных регрессий, эта модель требует достаточного количества репрезентативного исходного материала по каждой литофации, для настройки коэфициентов регрессии (в данном случае для обучения нейронной сети).
Cut-off
Для абсолютного большинства коллекторов проницаемость низкопоровых пропластков имеет более слабую зависимость от пористости.
В первом приближении этой зависимостью можно пренебречь и считать проницаемость в зоне малых значений пористости равной нулю:
(15) | k = \left\{ \begin{aligned} k_0 \, (\phi - \phi_0 )^{m} , \, \phi> \phi_0 \\ 0 , \, \phi \le \phi_0 .\\ \end{aligned} \right. |
где \phi_0 – минимальная пористость, определяющая отсечку проницаемого коллектора.
See also
Petroleum Industry / Upstream / Subsurface E&P Disciplines / Petrophysics / Reservoir Data Logs (RDL) @model