DETERMINACIÓN DE LA PERMEABILIDAD AL LÍQUIDO A PARTIR DE PRUEBAS

DE DESPLAZAMIENTO

PRESENTADO POR GRUPO 01-01:

María Camila Ospina M. Cod. 20161146079
Ana María Garzón Puentes Cod. 20162152466
Juan José Soto Ramírez Cod. 20161149334

PRESENTADO A:
Ing. Javier Andrés Martínez Pérez
Docente de la asignatura “Análisis de núcleos”

FACULTADE DE INGENIERÍA
PROGRAMA DE PETRÓLEOS
20 DE NOVIEMBRE DEL 2019

1
 CONTENIDO

1. OBJETIVOS ........................................................................................................................... 3
1.1 OBJETIVO GENERAL ........................................................................................................ 3
1.2 OBJETIVOS ESPECIFICOS ............................................................................................... 3
2. MATERIALES Y EQUIPOS ................................................................................................. 4
3. MARCO TEÓRICO................................................................................................................ 5
4. PROCEDIMIENTO ................................................................................................................ 7
5. TABLAS DE DATOS ............................................................................................................ 8
6. MUESTRA DE CÁLCULOS ................................................................................................. 9
7. RESULTADOS..................................................................................................................... 11
8. ANÁLISIS DE RESULTADOS ........................................................................................... 11
9. CONCLUSIONES Y RECOMENDACIONES ................................................................... 12
10. CUESTIONARIO .............................................................................................................. 13
11. TALLER ............................................................................................................................ 21
12. BIBLIOGRAFIA ............................................................................................................... 33

2
 1. OBJETIVOS

1.1 OBJETIVO GENERAL

Determinar la permeabilidad absoluta a partir de pruebas de desplazamiento de fluidos a

través de núcleos a condiciones de presión y temperatura del yacimiento.
1.2 OBJETIVOS ESPECIFICOS

 Familiarizarse con el funcionamiento y manejo del equipo de desplazamiento

positivo, “PERMEÁMETRO LÍQUIDO” y de cada uno de sus subsistemas
constituyentes para llevar a cabo una prueba con resultados representativos.

 Identificar las ventajas y desventajas del equipo empleado para la realización de la

prueba.

 Aprender acerca de la relación existente entre la permeabilidad de una muestra y la

saturación.

 Identificar las diferencias que se presentan en los dos métodos empleados para hallar
la permeabilidad de una muestra mediante el permeámetro de gas y permeámetro de
líquidos.

 Afianzar el conocimiento en el empleo de la ecuación de Darcy para el cálculo de la

permeabilidad absoluta

 Reconocer el concepto de permeabilidad como un fundamento primordial para el

estudio y análisis de un yacimiento y principalmente en el método de recobro y
estimulación empleados en el mismo.

3
 2. MATERIALES Y EQUIPOS

ELEMENTOS EVIDENCIA FOTOGRÁFICA

PERMEÁMETRO LÍQUIDO

FLUIDOS DE DESPLAZAMIENTO Agua de formación, salmuera sintética, aceite

mineral, crudo.

MUESTRA DE NÚCLEO

PROBETA DE RECOBRO

4
 3. MARCO TEÓRICO

La permeabilidad, detonada por K, es la capacidad del medio poroso para dejar pasar los
fluidos a través de él. Matemáticamente se expresa por la ley de Darcy y es una medida del grado
y tamaño en que los espacios porosos están interconectados. Además, es una medida de la
conductividad de los fluidos y, por su analogía con los conductores de electricidad, también se
define como el reciproco de la resistencia que un medio poroso ofrece al flujo de fluidos (Paris
de ferrer, 2009).
Flujo Lineal del Fluidos incompresibles (Martínez Pérez, 2019)
Si el medio poroso es un sistema lineal horizontal, se supone que el flujo ocurre a través de un
área A seccional constante, de tal forma que los extremos del sistema sean planos paralelos
abiertos al flujo, y que la presión en cada extremo del sistema sea constante sobre la superficie.
Si un fluido incompresible fluye a través del elemento dx, entonces la velocidad del fluido v y la
tasa de flujo q son constantes en todos los puntos. El comportamiento de flujo en este sistema
puede expresarse en forma diferencial por medio de:

Resolviendo la integral anterior y expresándola en unidades de campo:

Donde:
k = Permeabilidad Absoluta, md
μ = Viscosidad del fluido, cp
A = Área Seccional, cm2
q = Tasa de Flujo, bbl/d
P = Presión, psi
L = Distancia, pies.
La diferencia de presión en la ecuación anterior no es la única fuerza de empuje en el yacimiento
También existe la gravitacional, que debe tomarse en cuenta para determinar la dirección y la
tasa de flujo, puesto que la fuerza que ejerce el gradiente líquido está siempre dirigida
verticalmente hacia abajo, mientras que los resultados de la que se genera debido a una caída de
presión pueden ser en cualquier dirección. La fuerza que causa el flujo sería entonces la suma de

5
estos dos vectores. La energía que hay disponible para el movimiento de fluidos en cada punto
de una región determinada se llama potencial de flujo, por lo tanto, para que haya flujo entre dos
puntos, se necesita una diferencia de potencial. El potencial en determinado punto se puede
expresar como:

Donde:
Ø= Potencial de flujo
P = Presión en el punto estudiado
g = Aceleración de la gravedad
ρ = Densidad del fluido
Δz = distancia vertical desde un punto de referencia en el yacimiento hasta un punto
seleccionado.
Cuando la gravedad se desprecia, el segundo término en la ecuación anterior se hace cero y se
habla de un delta de presión para que exista flujo en determinada dirección.

Tabla 1. Clasificación de los valores de permeabilidad1

1
Tomado de: https://es.slideshare.net/estefaniaeminem/curso-caracterizacin-dinmica-11-nov-12

6
4. PROCEDIMIENTO

7
 5. TABLAS DE DATOS

Permeabilidad T DP alta DP baja Caudal Viscosidad

(md) (°F) (psia) (psia) (cc/min) (cP)
47.87 30 0 11.54 3 1
45.09 30 0 3.62 3 1
50.35 30 0 6.49 3 1
51.69 30 182.59 5 3 1
42.9 30 0 8.6 3 1
60 30 0 8.86 3 1
34.27 30 168.83 8.4 3 1
48.71 30 106.9 8.4 3 1
42.06 30 0 3.87 3 1
55.06 30 210.11 13.03 3 1
55.53 30 0 8.37 3 1
50.22 30 611.49 2.68 3 1
44.27 30 159.66 5.53 3 1
56.18 30 0 9.15 3 1
45.51 30 150.48 0 3 1
42.61 30 0 12.67 3 1
45.4 30 95.44 5.82 3 1
31.39 30 175.71 7.8 3 1
59.81 30 10.57 10.26 3 1
40.57 30 221.58 9.64 3 1
48.71 30 0 9.25 3 1
57.85 30 0 11.77 3 1
38.17 30 242.22 6.15 3 1
35.33 30 28.92 7.11 3 1
51.15 30 40.39 5.3 3 1
51.42 30 0 6.74 3 1
49.33 30 299.56 4.38 3 1
42.15 30 0 2.77 3 1
47.75 30 86.26 4.65 3 1
43.47 30 0 2.63 3 1
42.15 30 26.63 4.33 3 1
42.9 30 418.83 8.56 3 1
54.91 30 0 4.19 3 1
Resultados con el permeámetro del liquido

L(cm) 6.7
D(cm) 3,75
A(cm2) 11,045
Viscosidad(cp) 1
Datos del núcleo y fluido de desplazamiento

8
 6. MUESTRA DE CÁLCULOS

PARAMETROS ESTADISTICOS A CALCULAR:

 PROMEDIO
∑𝑛𝑖=1 𝐾𝑖
̅=
𝐾
∑𝑛𝑖=1 𝑛𝑖

47.87 + 45.09 + 50.35 + 51.69 … … . +43.47 + 42.15 + 42.9 + 54.91

̅=
𝐾
33

1554.78
̅=
𝐾 = 47.114545 𝑚𝑑
33
 DESVIACIÓN ESTÁNDAR

∑𝑛𝑖=1(𝐾 − 𝐾̅ )2
𝜎=√
(𝑛 − 1)

(47.87 − 47.1145)2 + (45.09 − 47.1145)2 + ⋯ + (42.9 − 47.1145)2

𝜎=√
(33 − 1)

𝜎 = 7.16284 𝑚𝑑
 NÚMERO DE INTERVALOS:
1 + 3,322 ∗ 𝐿𝑂𝐺(33) = 6.044

 TAMAÑO DE INTERVALOS:
𝑉𝑎𝑙𝑜𝑟 𝑚á𝑥𝑖𝑚𝑜 − 𝑉𝑎𝑙𝑜𝑟 𝑚í𝑛𝑖𝑚𝑜 = 59.81 𝑚𝑑 − 31.39 𝑚𝑑 = 28.42 𝑚𝑑

 TAMAÑO DE CLASE:
𝑡𝑎𝑚𝑎ñ𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜𝑠 28.42
= = 4.7366667 𝑚𝑑
𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜𝑠 6
9
Intervalos de permeabilidad
Frecuencia Frecuencia
Límite Marca de Frecuencia Frecuencia
Límite superior absoluta relativa
inferior clase absoluta (fi) relativa
acumulada (Fi) acumulada
31.39 59.8 45.60 3 3 9.1% 9.0

59.81 88.23 74.02 2 5 6.1% 15.1

88.23 116.65 102.44 11 16 33.3% 48.4

116.65 145.07 130.86 6 22 18.2% 66.6

145.07 173.49 159.28 6 28 18.2% 84.8

173.49 201.91 187.70 5 33 15.2% 100.0

201.91 230.33 216.12 0 0.0%

SUMA 33 100%
DISTRIBUCION ESTADISTICA DE LA PERMEABILIDAD K (md)

12

10
Frecuencia absoluta

6 Histograma de
frecuencia
4 Poligono de
frecuencias
2

0
33.76 38.50 43.23 47.97 52.71 57.44 62.18
Marca de clase

10
 7. RESULTADOS

MEDIDA VALOR
MEDIA ARITMETICA 47.1145
DESVIACION ESTANDAR 7.16284
MAXIMO 59.81
MINIMO 31.39
TAMAÑO DE INTERVALOS 28.42
TAMAÑO DE CLASES 4.7366667
MODA 42.9
MEDIANA 47.75

8. ANÁLISIS DE RESULTADOS
El análisis de los resultados parte desde la graficacion de los mismo en un histograma de
permeabilidad, acá se puede observar la gran variación y distribución de los datos por lo que
inferir un comportamiento es bastante relativo. Sin embargo, un intervalo resalta sobre los demás
y es allí donde encontramos la media aritmética de 47.11. todas estas variaciones pueden deberse
principalmente a que muchos de los poros no están interconectados por lo cual la tortuosidad
(indirectamente la permeabilidad) se verá afectada dando resultado alejados de la absoluta u otro
error pudo ser que varias veces la prueba se corrió y se detuvo antes de que el diferencial de
presión fuera 0 o cercano a el, habiendo resultados de mas de 600 psi de diferencia.
En la grafica se puede observar la línea de polígono de frecuencia la cual ayuda a ver el
comportamiento de los datos, asimilando al comportamiento a la campana de gauss demostrando
que en el P50 se encuentra la media aritmética.
Al usar agua u algun otro liquido se corre el peligro de que este reaccione con la roca, por lo
cual la humectabilidad puede interactuar y hacer que este se pegue a la muestra lo que va a hacer
que el de desplazamiento se vea un poco interrumpido además de que al usar este metodo que es
muy invasivo se debe dejar un buen tiempo la muestra secar para llevarla a otra prueba o bien se
puede secar en horno para acelerar el proceso
Por este motivo se puede considerar que el metodo para hallar permeabilidad por gas es mas
acertado y que aparte de que no hay reacción con la roca y la invasión es cero se tiene la
correcion por efecto de klikenberg para tener mayor certeza del resultado.
Para finalizar, todos los datos fueron sumistrados por lo cual no se tiene certeza de que la prueba
se haya corrido de la mejor manera, si no hubieron problemas durante estas o que tipo de muestra
fue la que se uso, ya que si esta podría haber contenido un mínimo de arcillas estas podrían
hidratarse y hacer la que la permeabilidad varie aun mas, además la calibración del equipo, si

11
nunca se filtro aire a la prueba etc. Son muchos los errores que se pueden llevar al no haber
realizado la prueba por lo que la discreción con el manejo de datos debe ser alto.

9. CONCLUSIONES Y RECOMENDACIONES
Los errores alrededor de esta prueba anteriormente mencionados son muchos, las presiones no
fueron las que dictan en la guía, la reactividad por acción de la humectabilidad es casi seguro que
estuvo, la realización de la prueba y problemas al llevarla a cabo no se sabe como estuvo; todos
estos ayudan a que la permabilidad absoluta no sea la real de la roca, además de que no hay
forma de tener una correcion de dicho valor.
Los resultados que se obtienen con esta prueba no son tan deseados como lo pueden ser en otras
pruebas por lo que se recomiendo que para determinar la permeabilidad se lleve a cabo otra
prueba.
Según datos históricos de Colombia la permeabilidad cercana a 47mD es algo elevado al
promedio de los yacimientos colombianos por lo que si de la muestra se logra comprobar que
había hidrocarburo será un buen prospecto para la explotación.

12
 10. CUESTIONARIO

1. ¿Cuál es la definición de una arena que tiene una permeabilidad de un Darcy? (Paris
de ferrer, 2009)
Una arena que tiene permeabilidad de un Darcy es aquella en la cual un fluido con una
viscosidad de un 1centipoise, se moverá a una velocidad de un centímetro por segundo bajo un
gradiente de presión de una atmósfera por centímetro.
2. ¿Cuáles son las dimensiones de la permeabilidad? (Paris de ferrer, 2009)
La ley de Darcy para flujo lineal viene dada por la ecuación:
𝐾𝜌 𝑑
𝑣=−
𝜇 𝑑𝑙
La cual escrita con las dimensiones correspondientes son:
𝐿 [𝑀/𝐿3 ] [𝐿2 /𝑇 2 ]
= [𝐾] ( )
𝑇 [𝑀/𝐿𝑇] [𝐿]

Por lo tanto, 𝐾 = [𝐿2 ]

La unidad fundamental SI de permeabilidad es el m2, se define de la siguiente manera: una
permeabilidad de un metro cuadrado permitirá un flujo de 1 m3/s de fluido de una viscosidad de
1 Pa·s a través de un área de 1 m2 bajo una gradiente de presión de 1 Pa/m. Un darcy es igual a
0.986923 x 10-12 m2.
Debido al exponente tan pequeño, esta es una unidad difícil para uso común. La unidad de
permeabilidad preferida por la Sociedad de Ingenieros de Petróleos de AIME es el micrómetro
cuadrado (μm2). Un darcy (la unidad tradicional es igual a 0.986923 μm2.)
3. ¿Cuáles son los tipos de permeabilidad que existen? Explique. (Paris de ferrer, 2009)
Los tipos de permeabilidad son:
 Permeabilidad absoluta, k: es aquella que se obtiene cuando existe una sola fase que
satura 100% el medio poroso.
 Permeabilidad efectiva, ki: es una medida de la conductancia del medio poroso hacia un
fluido específico cuando existe más de una fase en él, las cuales fluyen simultáneamente
(petróleo, gas o agua). Esta es función de la saturación del fluido considerado y es
siempre menor que la permeabilidad absoluta.
 Permeabilidad relativa, kri: es la relación entre la permeabilidad efectiva de un fluido
(petróleo, gas o agua) a una saturación dada y la permeabilidad absoluta del medio
poroso. Esta también es función de la saturación del fluido y siempre será menor o igual a
la unidad, pero si en la roca sólo existe un fluido, su permeabilidad relativa es igual a 1.0.

13
 4. ¿Cuáles son los métodos matemáticos para determinar la permeabilidad promedia de
un yacimiento? (Paris de ferrer, 2009)
Debido a que no es homogéneo en toda su extensión, para determinar la permeabilidad promedio
se tienen tres métodos matemáticos:
1) Permeabilidad ponderada por espesor:
El yacimiento debe estar compuesto por estratos paralelos separados por barreras impermeables
que impidan el flujo a través de ellos, como se muestra en la figura:

Figura 1. Flujo lineal a través de estratos paralelos.

En este caso, el caudal total está representado por la sumatoria del caudal en cada estrato:

𝑘̅𝑎ℎ𝑡 ∆𝑃 𝑘1 𝑎ℎ1 ∆𝑃 𝑘2 𝑎ℎ2 ∆𝑃 𝑘3 𝑎ℎ3 ∆𝑃

𝑞𝑡 = = + +
𝜇𝐿 𝜇𝐿 𝜇𝐿 𝜇𝐿
Como se evidencia en la figura, la caída de presión (∆𝑃), longitud (L), viscosidad (𝜇) y ancho
(a) en cada estrato, es igual, por lo se puede simplificar la ecuación para obtener la permeabilidad
promedio (𝑘̅):
∑𝑛𝑖=1 𝑘𝑖 ℎ𝑖
𝑘̅ =
∑𝑛𝑖=1 ℎ𝑖
Donde:
𝑘𝑖 = Permeabilidad de cada estrato.
ℎ𝑖 = Espesor de cada estrato.
Si el ancho de cada estrato difiere, el promedio se obtiene reemplazando el espesor por el área de
cada estrato:

14
 ∑𝑛𝑖=1 𝑘𝑖 𝐴𝑖
𝑘̅ =
∑𝑛𝑖=1 𝐴𝑖

Las anteriores ecuaciones pueden ser aplicadas tanto para flujo lineal como para flujo radial
paralelo.
2) Permeabilidad promedio armónico
El yacimiento presenta variación de la permeabilidad lateralmente, por lo que los estratos se
encuentran dispuestos en serie, como lo muestra la siguiente figura:

Figura 2. Flujo lineal a través de estratos en serie.

A diferencia del caso anterior, permanece la tasa de flujo constante a través de las capas,
variando la caída de presión. La caída de presión total es igual a:

𝑞𝜇𝐿𝑇 𝑞𝜇𝐿1 𝑞𝜇𝐿2 𝑞𝜇𝐿2

∆𝑃𝑇 = = + +
𝐴𝑘̅ 𝐴𝑘1 𝐴𝑘2 𝐴𝑘2

Expresión que se puede simplificar si se considera que el caudal (q), viscosidad (𝜇) y área (A) en
cada uno de los estratos es igual. Por lo anterior, despejando la permeabilidad promedio se
obtiene que esta será igual a:

∑𝑛𝑖=1 𝐿𝑖
𝑘̅ =
∑𝑛𝑖=1(𝐿/𝑘)𝑖

15
Donde:
𝑘𝑖 = Permeabilidad de cada estrato. 𝐿𝑖 = Longitud de cada estrato.
Si el sistema es radial, como se ilustra en la siguiente figura, se debe seguir el mismo
procedimiento anterior de sumar las caídas de presión en cada estrato teniendo en cuenta los
radios (r):

Figura 3. Flujo radial a través de estratos en serie.

𝑟 𝑟 𝑟 𝑟
𝑞𝜇ln (𝑟𝑒 ) 𝑞𝜇ln (𝑟1 ) 𝑞𝜇ln (𝑟2 ) 𝑞𝜇ln (𝑟𝑒 )
𝑤 𝑤 1 2
∆𝑃𝑇 = = + +
2𝜋ℎ𝑘̅ ̅̅̅1
2𝜋ℎ𝑘 ̅̅̅2
2𝜋ℎ𝑘 ̅̅̅3
2𝜋ℎ𝑘
Simplificando, la permeabilidad promedio se puede obtener con la siguiente ecuación:
𝑟
ln (𝑟𝑒 )
𝑤
𝑘̅ = 𝑟
ln (𝑟 𝑖 )
𝑖−1
∑𝑛𝑖=1
𝑘𝑖

3) Permeabilidad por promedio geométrico:

Procedimiento que se define mediante la siguiente ecuación:
∑𝑛𝑖=1 ℎ𝑖 𝑙𝑛𝑘𝑖
̅
𝑘 = 𝑒𝑥𝑝 [ ]
∑𝑛𝑖=1 ℎ𝑖

Donde:
𝑘𝑖 = Permeabilidad de cada estrato.
ℎ𝑖 = Espesor de cada estrato.

16
Si los espesores son iguales, la ecuación se puede simplificar a:

𝑘̅ = (𝑘1 ∗ 𝑘2 … . .∗ 𝑘𝑛 )1/𝑛
Donde n es el número total de muestras.
5. Si la muestra no está completamente saturada, ¿el valor de la permeabilidad será alto o
bajo? ¿Por qué?
Teniendo en cuenta que para determinar el valor de permeabilidad usando la ecuación de Darcy
se debe cumplir que la muestra se encuentre saturada al 100% por una sola fase, al no
encontrarse completamente saturada el resultado obtenido será menor y además incorrecto, por
lo que no sería representativo de la muestra.
6. ¿La permeabilidad absoluta medida con agua es igual a la determinación de la misma
medida con otro líquido (por ejemplo el petróleo)? ¿Por qué? (Paris de ferrer, 2009)
La permeabilidad absoluta es la capacidad que tiene el medio poroso de dejar pasar un fluido a
través de él, por lo cual es una propiedad de la roca y no está influenciada por el tipo de fluido
que la satura y pasa por ella, por lo que sí es medida con agua o con otro fluido el valor obtenido
para dicha propiedad petrofísica debe ser el mismo. Sin embargo, para una medición correcta de
permeabilidad es requisito que el fluido a utilizar no reaccione con la roca y que sature 100% el
espacio poroso de la roca.
Los factores que la afectan la permeabilidad son las características de la roca (empaquetamiento,
tamaño y distribución de los granos, el contenido de arcilla, la compactación de la roca).
7. Investigar sobre la presión de confinamiento a que son sometidas las muestras en las
diferentes pruebas de laboratorio.
Existen diferentes ensayos de resistencia de laboratorio los cuales son:
* Ensayo Triaxial2:
Permite realizar ensayos en una variedad de procedimientos para determinar la resistencia al
cortante, la rigidez, características de deformación de las muestras, consolidación y
permeabilidad. Se ensayan muestras cilíndricas dentro de una membrana delgada de caucho,
colocándolas dentro de una celda triaxial con dos tapas rígidas y pistones arriba y debajo de la
muestra.
La celda se llena de un fluido especial, se aplica una presión determinada sobre el fluido (σ3), la
cual se transmite por éste a la muestra. Los esfuerzos de cortante se aplican mediante fuerzas de
compresión verticales accionadas por pistones. La presión de poros dentro de la muestra puede
medirse a través de un pequeño tubo o bureta en contacto con la muestra. Para cada presión de
confinamiento se obtiene el esfuerzo desviador (∆σ) que se requiere para hacer fallar la muestra.

2
Tomado de: http://www.kerwa.ucr.ac.cr/bitstream/handle/10669/13369/8391-11758-1-SM.pdf?sequence=1

17
El drenaje de la muestra se realiza a través de las piedras porosas y el cambio de volumen de
agua puede medirse. Alternativamente, si no se permite drenaje, se puede medir la presión de
poros. Realizando varias pruebas se puede obtener la envolvente de Mohr para un suelo
determinado. El comportamiento Esfuerzo–deformación es determinado por la presión de
confinamiento, la historia de esfuerzos y otros factores.
 Formas de realizar el ensayo Triaxial:
a. Ensayo Consolidado drenado: El ensayo se realiza lentamente para permitir el drenaje del agua
dentro de la muestra e impedir que se puedan generar presiones de poros. La velocidad de ensayo
debe ser tal, que las fluctuaciones en la presión de poros sean despreciables y en cualquier caso
no superiores al 5% de la presión efectiva de confinamiento
b. Ensayo Consolidado no drenado, con medición de presión de poros: Se permite el drenaje
durante la aplicación del esfuerzo de confinamiento colocándolo lentamente, pero se impide
durante la aplicación del esfuerzo desviador. Los ensayos no drenados deben realizarse a una rata
que no exceda una deformación unitaria del 2% por hora.
c. Ensayo No consolidado - No drenado o ensayo rápido: No se permite el drenaje durante la
aplicación de la presión de confinamiento y el esfuerzo desviador.
 Fuentes de error:
-Las muestras tienden a deformarse como un barril, lo cual conduce a sobreestimación de la
resistencia al cortante.
-En el proceso de saturación, la muestra tiende a alterarse por cambio de volumen, lo cual puede
determinar una pérdida de resistencia.
 Características de la muestra
-Razón Largo/ancho de la probeta debe ser 2.0 a 2.5:1
-Extremos deben ser paralelos y pulidos, sin grietas.
-Ancho de muestra debe ser >10 veces el tamaño medio del grano.

* Ensayo de Corte Directo3:

Es simple y económico de realizar, pero presenta los inconvenientes del poco control que se
tiene sobre las condiciones de drenaje, la dificultad para medir presiones de poro y algunos
problemas inherentes a los mecanismos de las máquinas que realizan los ensayos. Las ventajas
de los ensayos de Corte Directo son su facilidad de ejecución, la cual permite la realización de
una cantidad grande de pruebas en poco tiempo y la posibilidad de realizar ensayos sobre
superficies de discontinuidad. En este ensayo la resistencia al cortante puede medirse en un plano

3
Tomado de :229-3_esfuerzoyresistenciaalcortante.pdf

18
predeterminado, cortando la muestra con una determinada orientación. La superficie de falla es
predefinida y no depende de las propiedades del suelo.
La muestra se coloca en una caja compuesta por dos anillos, uno superior y otro inferior, los
cuales pueden desplazarse horizontalmente el uno con respecto al otro al aplicarse una fuerza de
cortante. Las muestras no pueden saturarse completamente pero un grado de saturación
relativamente alto se puede obtener sumergiendo la muestra en agua por un periodo largo de
tiempo, antes del ensayo. Sin embargo, debe tenerse mucho cuidado con los efectos de
saturación sobre algunos materiales, especialmente los suelos expansivos.
Se dibuja una curva esfuerzo-deformación para cada ensayo, en la cual se determinan los valores
de la resistencia máxima y la resistencia residual. Se realizan varias pruebas para el mismo tipo
de suelo con diferentes presiones normales y se dibuja la envolvente de falla para obtener
gráficamente los valores de cohesión y ángulo de fricción. Se recomienda un mínimo de cinco
pruebas para cada tipo de suelo.
*Ensayos con deformación controlada o con esfuerzo controlado:
El esfuerzo de corte puede ser aplicado incrementando los esfuerzos en forma gradual y
midiendo la deformación producida (Esfuerzo controlado) o moviendo las partes del equipo a un
desplazamiento dado y midiendo el esfuerzo resultante (deformación controlada). Este tipo de
ensayo no puede determinar el esfuerzo pico y la resistencia residual en forma precisa.
*Ensayo de compresión simple o uniaxial:
Es un ensayo de compresión de una muestra cilíndrica con una relación diámetro longitud 1:2.
La muestra es comprimida axialmente hasta que ocurre la falla. La resistencia al cortante se
asume que es igual a la mitad de la resistencia a la compresión. Este ensayo es utilizado con
frecuencia para conocer la resistencia no drenada de suelos cohesivos. Debe tenerse en cuenta
que los resultados son expresados en términos de esfuerzos totales, debido a que no se mide
ningún tipo de presión de poros y los ensayos en limoso arenas o materiales fisurados no tienen
ninguna validez. El ensayo es solo aplicable a suelos cohesivos que no permiten la salida de agua
durante el proceso de carga. Generalmente, el valor de la resistencia no drenada se supone igual a
la mitad del valor de la resistencia encofinada.
Su = ½ qu
En todos los casos, debido a las incertidumbres asociadas con el ensayo, el muestreo y su
preparación, esta prueba de laboratorio solamente puede utilizarse como un estimado
aproximado de la resistencia en el sitio.
*Ensayo de carga puntual4:
Se emplea para determinar la resistencia a la compresión simple de fragmentos irregulares de
roca, a partir del índice de resistencia a la carga puntual (Is), de tal forma que el stress aplicado
se convierte a valores aproximados de resistencia a la compresión uniaxial según el diámetro de

4
Tomado de: Lab_1_Ensayos_carga_puntual_uniaxial_y_triaxial.pdf

19
la muestra. Las ventajas de esta prueba son que se pueden usar muestras de roca irregulares sin
preparación previa y que la máquina es portátil.
El procedimiento consiste en romper una muestra entre dos puntas cónicas metálicas accionadas
por una prensa.

8. Explicar el funcionamiento de la contrapresión. ¿Si la contrapresión es mayor a la

presión de producción existirá flujo?
La contrapresión es aquella presión ejercida de forma opuesta a la dirección de flujo del fluido,
siendo una condición perjudicial para el flujo y simula la presión que ejercen factores como los
fluidos de la perforación o la formación hacia los fluidos de la formación. Si se presenta la
situación en la que la contrapresión sea igual a la presión de producción, no existirá flujo del
fluido de producción porque se habrá establecido un equilibrio entre las presiones y los fluidos
no se moverán a través del medio poroso. Caso similar se presenta cuando la contrapresión es
mayor a la presión de producción ya que tampoco existirá flujo ya que los fluidos se desplazan
hacia la región de menor presión, por lo que no fluirán al pozo donde la presión es mayor.
Por lo anterior, que la contrapresión sea igual o mayor a la presión de producción existirá es una
condición no adecuada para el eficiente flujo de fluidos de la formación hacia el pozo.

20
 11. TALLER

1. Analice el rango de variación de permeabilidades de los diferentes tipos de rocas según

los estudios realizados en Freeze. R. A & Charry. J. A. y preséntelos en la figura 5.26
en el libro fundamentos de ingeniería de yacimientos de Magdalena Paris de Ferrer.

Figura 1. Permeabilidades según litología. (Paris de ferrer, 2009)

En la figura se puede observar el rango de variación de permeabilidad de los diferentes tipos

de rocas. Como era de suponer, las rocas que presentan las mayores permeabilidades son las
rocas sedimentarias, las calizas y las rocas fracturadas; la permeabilidad es una propiedad muy
importante y de mayor interés en la industria ya que esta es la que me va a permitir el fluido de
hidrocarburos en el yacimiento.

Mientras menos sea el tamaño de los granos que conformen la roca, menor será su
permeabilidad, siendo las de menor permeabilidad las arcillas (de 10-7 a 10-4 Darcy). A pesar de
ello, son de importancia para la industria petrolera ya que actúan como rocas sello debido a que
su baja permeabilidad impide la circulación de los fluidos hidrocarburos, permitiendo la
formación de trampas.

Las morrenas presentan una de las variaciones de permeabilidad más amplias, por el contrario las
gravas presentan variaciones más estrechas.

De otro lado, las rocas ígneas y metamórficas presentan permeabilidades de un orden menor a
10-4 Darcys, debido a que estos tipos de rocas presentan una mala porosidad y son muy
consolidadas. Por ello, a causa de a su baja permeabilidad son poco de poco interés para la
industria. Sin embargo, cuando se encuentran fracturadas, debido a procesos como alteraciones
hidrotermales, aumentan su permeabilidad, pasando a estar entre 10-3 y 101 Darcys, lo que puede
generar caminos migratorios para los hidrocarburos.

21
El rango de variación de la permeabilidad es una aproximación para ciertas rocas ya que el
estudio se realizó par a un determinado número de rocas, por lo que otras rocas, que pertenezcan
a un mismo tipo de los ilustrados en la figura, pueden encontrarse fuera del rango de
permeabilidad mostrado.

2. Investigar sobre las correlaciones para determinar la permeabilidad de la roca según el

tipo de litología. (Paris de ferrer, 2009)

Se encontró que además de ecuaciones para determinar la permeabilidad de la roca según la

litología, 
 se han desarrollado correlaciones para calcular la permeabilidad absoluta teniendo en
cuenta el fluido que satura la roca. Entre ellas se tienen:

-Ecuación de Timur: Para estimar la permeabilidad a partir de la saturación de agua connata y la

porosidad:

∅4.4
𝑘 = 8.58102 ∗ 2
𝑆𝑤𝑐

-Ecuación de Morris y Biggs: Para estimar la permeabilidad en yacimientos de petróleo y gas:

2
∅3
𝑘 = 62,5 ∗ ( ) 𝑃𝑎𝑟𝑎 𝑢𝑛 𝑦𝑎𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑑𝑒 𝑝𝑒𝑡𝑟ó𝑙𝑒𝑜
𝑆𝑤𝑐
2
∅3
𝑘 = 2,5 ∗ ( ) 𝑃𝑎𝑟𝑎 𝑢𝑛 𝑦𝑎𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑑𝑒 𝑔𝑎𝑠
𝑆𝑤𝑐

Las correlaciones para determinar la permeabilidad relativa teniendo en cuenta la litología de la

roca son:

Tabla4. Correlación para determinar la permeabilidad según la litología.

Nombre de la Condiciones en las que se

Ecuación
correlación emplea:

22
 𝐾𝑟𝑔
= 𝐴 ∗ (0.0435 + 0.4556𝐴)
𝐾𝑟𝑜
Wahl y
Yacimiento de areniscas.
asociados
1 − 𝑆𝑔𝑐 − 𝑆𝑤𝑖 − 𝑆𝑜
𝐴=
𝑆𝑜 − 𝑆𝑜𝑟

𝑆 Arenas no consolidadas en
𝐾𝑟𝑤 = (1 − 𝑆)2 ; 𝐾𝑟𝑜 = 𝑆 3 ; 𝑆 = 1−𝑆𝑜
𝑤𝑖 procesos de drenaje.

𝑆𝑤 −𝑆𝑤𝑖 Arenas consolidadas en

𝐾𝑟𝑜 = (1 − 𝑆)2 ; 𝐾𝑟𝑤 = 𝑆 3 ; 𝑆 = 1−𝑆𝑤𝑖 procesos de imbibición.

𝐾𝑟𝑤 = (1 − 𝑆)4; 𝐾𝑟𝑜 = 𝑆 3 (2 − 𝑆);

Corey y Arenas consolidadas en
asociados 𝑆𝑜 procesos de drenaje.
𝑆=
1 − 𝑆𝑤𝑖

𝐾𝑟𝑤 = (1 − 𝑆)4; 𝐾𝑟𝑜 = 𝑆 3 (2 − 𝑆);

Procesos de drenaje en
𝑆𝑔 empuje de gas en solución.
𝑆=
1 − 𝑆𝑤𝑖

𝑆𝑤 −𝑆𝑤𝑖 Arenas no consolidadas

𝐾𝑟𝑜 = (1 − 𝑆)3 ; 𝐾𝑟𝑤 = 𝑆 3 ; 𝑆 = 1−𝑆𝑤𝑖 con grano bien clasificado.

𝐾𝑟𝑤 = (1 − 𝑆)2 ∗ (1 − 𝑆 1.5 ) Arenas no consolidadas

Willie con grano pobremente
𝐾𝑟𝑜 = 𝑆 3.5 clasificado.

𝐾𝑟𝑤 = (1 − 𝑆)2 ∗ (1 − 𝑆 2 )
Arenas cementadas y
formaciones limosas.
𝐾𝑟𝑜 = 𝑆 4

23
 𝑆𝑤 − 𝑆𝑤𝑐 3
𝐾𝑟𝑤 =( ) ;
1 − 𝑆𝑤𝑐
Areniscas y areniscas
Jones
arcillosas.
(𝑆𝑤 − 𝑆𝑤𝑐 )21
𝐾𝑟𝑜 =
(1 − 𝑆𝑤𝑐 )2

4 2
𝐾𝑟𝑜 = 𝑆𝑜𝑛 ; 𝐾𝑟𝑤 = 𝑆𝑤𝑛

𝑆𝑜 − 𝑆𝑜𝑟
𝑆𝑜𝑛 = ( )
Spivak 1 − 𝑆𝑜𝑟 − 𝑆𝑤𝑖

𝑆𝑤 − 𝑆𝑤𝑖
𝑆𝑤𝑛 = ( )
1 − 𝑆𝑜𝑟 − 𝑆𝑤𝑖

𝐾𝑟𝑤 = 𝑆 4

𝐾𝑟𝑜 = (1 − 2𝑆)1.25 ∗ [2 − √1 − 2𝑆] Yacimientos con empuje

Naar- de agua, procesos de
Henderson 𝑆𝑤 − 𝑆𝑤𝑖 imbibición en arenas
𝑆=
1 − 𝑆𝑤𝑖 consolidadas.

𝑆𝑤𝑖 < 𝑆𝑤 < 1 − 𝑆𝑜𝑟

𝐾𝑟𝑔 (1 − 𝑆 ∗ )2 ∗ [1 − (𝑆 ∗ )2 ]
=
𝐾𝑟𝑜 (𝑆 ∗ )4
Torcaso y Procesos de drenaje en
Willie areniscas no consolidadas.
𝑆𝑜
𝑆∗ =
1 − 𝑆𝑤𝑖

𝐾𝑟𝑤 = √𝑆𝑤𝑛 ∗ 𝑆𝑤3

Rocas con porosidad
𝐾𝑟𝑔 = (1 − 𝑆𝑤𝑛 ) ∗ [𝑆𝑤 ∗ √1 − 𝑆𝑤𝑛 ] intergranular y para flujo
Pirson
de gas en condiciones de
𝑆𝑤 − 𝑆𝑤𝑖 drenaje.
𝑆𝑤𝑛 =
1 − 𝑆𝑤𝑖

24
 Rocas con porosidad
𝐾𝑟𝑤 = √𝑆𝑤𝑛 ∗ 𝑆𝑤3
intergranular y para flujo
simultaneo de agua y
1 − 𝑆𝑤 − 𝑆𝑤𝑖
𝐾𝑟𝑜 = petróleo en condiciones de
1 − 𝑆𝑜𝑟 − 𝑆𝑤𝑖 imbibición.

𝑆𝑔3
𝐾𝑟𝑔 =
(1 − 𝑆𝑤𝑖 )3
Sistemas trifásicos
preferiblemente mojados
𝑆𝑜3
𝐾𝑟𝑜 = por agua en condiciones de
(1 − 𝑆𝑤𝑖 )3
drenaje para areniscas no
consolidadas.
(𝑆𝑤 − 𝑆𝑤𝑖 )3
𝐾𝑟𝑤 =
(1 − 𝑆𝑤𝑖 )3

𝐾𝑟𝑔
Willie y 𝑆𝑔2 ∗ (1 − 𝑆𝑤𝑖 )2 − (𝑆𝑤 + 𝑆𝑜 − 𝑆𝑤𝑖 )2
Gardner =
(1 − 𝑆𝑤𝑖 )4 Sistemas trifásicos
preferiblemente mojados
𝑆𝑜3 ∗ (2𝑆𝑤 + 𝑆𝑜 − 2𝑆𝑤𝑖 ) por agua en condiciones de
𝐾𝑟𝑜 = drenaje para areniscas
(1 − 𝑆𝑤𝑖 )4
consolidadas.
(𝑆𝑤 − 𝑆𝑤𝑖 )4
𝐾𝑟𝑤 =
(1 − 𝑆𝑤𝑖 )4

Las correlaciones de Willie y Gardner pueden aplicarse a sistemas

mojados por petróleo si la saturación de agua Sw se cambia por la
saturación de petróleo.

Estimación estadística de la
permeabilidad relativa al
𝐾𝑟𝑜 =
petróleo en sistemas
Stone trifásicos con datos de
(𝐾𝑟𝑜𝑤 + 𝐾𝑟𝑤 ) ∗ (𝐾𝑟𝑜𝑔 + 𝐾𝑟𝑔) − (𝐾𝑟𝑤
permeabilidades relativas
+ 𝐾𝑟𝑔 )
para sistemas agua petróleo
y gas

25
3. Al yacimiento San Francisco se le corazonó la formación productora con un lodo base
agua y se solicitó al ICP, determinar la porosidad y la permeabilidad del intervalo
productor de caballos superior del SF-25. La tabla resume los resultados de la
porosidad y permeabilidad pie a pie de dicho intervalo. Calcular la distribución de la
permeabilidad del intervalo estadísticamente y cuánto vale la permeabilidad al 50%.
También halle el promedio geométrico de la permeabilidad. 


Tabla5. Datos de porosidad y permeabilidad para el pozo SF-25

Profundidad (ft) Porosidad (%) Permeabil

idad (mD)

2600 2601 0,1790 295

2601 2602 0,1740 262

2602 2603 0,1570 88

2603 2604 0,2000 87

2604 2605 0,1502 168

2605 2606 0,1913 71

2606 2607 0,1259 62

2607 2608 0,1338 187

2608 2609 0,1683 469

2609 2610 0,1873 77

2610 2611 0,1733 127

2611 2612 0,1600 161

2612 2613 0,1172 50

2613 2614 0,2015 58

2614 2615 0,1424 109

2615 2616 0,1869 228

2616 2617 0,2068 282

2617 2618 0,1963 776

26
 2618 2619 0,1678 87

2619 2620 0,1727 47

2620 2621 0,1690 16

2621 2622 0,1711 35

2622 2623 0,1397 47

2623 2624 0,1606 54

2624 2625 0,1858 273

2625 2626 0,1590 454

2626 2627 0,1864 308

2627 2628 0,1641 159

2628 2629 0,1871 178

2629 2630 0,1294 209

2630 2631 0,1654 216

2631 2632 0,1770 241

2632 2633 0,1745 333

2633 2634 0,1720 386

2634 2635 0,1675 429

∑ 5.903

-Permeabilidad promedio por distribución estadística:

los necesarios son 7 y tienen una amplitud de 108 mD. Los resultados son:

Tabla6. Permeabilidad promedio estadístico.

27
 Intervalos de porosidad
Marca de Frecuencia FRECUENCIA FRECEUNCIA RELATIVA
Lim inferior Lim superior
clase absoluta (fi) REALATIVA ACUMULADA
16 124.5714286 70.2857143 14 40.000 40
124.571429 233.1428571 178.857143 9 25.714 65.71428571 Media 200.828571
233.142857 341.7142857 287.428571 7 20.000 85.71428571 Desviacion estandar 161.973761
341.714286 450.2857143 396 2 5.714 91.42857143 maximo 776
450.285714 558.8571429 504.571429 2 5.714 97.14285714 minimo 16
558.857143 667.4285714 613.142857 0 0.000 97.14285714 tamaño de intervalos 760
667.428571 776 721.714286 1 2.857 100 numero de intervalos 7
suma 35 100 tamaño de clases 108.571429

Posteriormente, se realizó la gráfica de frecuencia acumulada Vs. marca de clase:

Figura 5. Histograma: distribución estadística de la permeabilidad.

En la gráfica de frecuencia acumulada se lee el valor de la permeabilidad promedio, la cual

corresponde al 50%, siendo de 159 mD.

28
- Permeabilidad por promedio geométrico:

Para el presente promedio se emplea la siguiente fórmula:

∑𝑛 ℎ ∗ln 𝑘𝑖
[ 𝑖=1𝑛 𝑖 ]
̅=𝑒
𝐾 ∑𝑖=1 ℎ𝑖

Como el espesor de cada uno de los intervalos es de 1 ft, la fórmula se reduce a:

1
̅ = (𝑘1 ∗ 𝑘2 ∗ 𝑘3 … ∗ 𝑘𝑛 )𝑛
𝐾

1
̅ = (3.58987 ∗ 1075 )35
𝐾

̅ = 144.12 𝑚𝐷
𝐾

4. Dados los siguientes datos de permeabilidades obtenidas de un análisis de núcleos, se

desea estimar la permeabilidad promedio del yacimiento.

Tabla7. Datos de permeabilidad y profundidad del yacimiento.

Profundidad (ft) Permeabilidad (mD)

3998-4002 200
4002-4004 130
40004-4008 170
4008-4010 180
4010-4013 140

La permeabilidad promedio por peso es igual a:

∑𝑛𝑖=1 𝐾𝑖 ∗ ℎ𝑖
̅=
𝐾
∑𝑛𝑖=1 ℎ𝑖
Dónde:
𝐾𝑖 : 𝑃𝑒𝑟𝑚𝑒𝑎𝑏𝑖𝑙𝑖𝑑𝑎𝑑 𝑒𝑛 𝑐𝑎𝑑𝑎 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜
ℎ𝑖 : 𝐸𝑠𝑝𝑒𝑠𝑜𝑟 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑝𝑎

El espesor es igual a la diferencia entre los intervalos:

29
 ℎ1 = 4002 𝑓𝑡 − 3998 𝑓𝑡
ℎ1 = 4 𝑓𝑡
Luego, se multiplica la permeabilidad del intervalo por el espesor del intervalo:

𝐾𝑖 ∗ ℎ𝑖 = 200 ∗ 4𝑓𝑡 = 800 𝑚𝐷 𝑓𝑡

Posteriormente, se suma la multiplicación de la permeabilidad de cada intervalo con el espesor y

se suma el espesor total del núcleo. Los resultados se muestran a continuación:
Tabla8. Permeabilidad promedio por peso.

Profundidad Permeabilidad hi (ft) Ki*hi (mD*ft)

(ft) (mD)
3998-4002 200 4 800
4002-4004 130 2 260
40004-4008 170 4 680
4008-4010 180 2 360
4010-4013 140 3 420
∑ 15 2520

Para finalizar, se reemplazan las sumatorias en la fórmula de la permeabilidad promedio por

peso:
2520 𝑚𝐷 ∗ 𝑓𝑡
̅=
𝐾 ̅ = 168 𝑚𝐷
𝐾
15 𝑓𝑡

5. Sea un yacimiento formado por
 5 estratos de igual espesor se encuentran conectados

en serie. La longitud de cada segmento y sus respectivas permeabilidades son:

Tabla9. Datos de permeabilidad y longitud de los estratos productores del yacimiento.

Longitud(ft) Permeabilidad
(mD)
150 10
200 20
300 30
500 50
200 80

30
Calcule la permeabilidad promedio del yacimiento, suponiendo: a) sistema de flujo lineal, b)
sistema de flujo radial.

Si el sistema de flujo es lineal, la permeabilidad promedio es igual a:

∑𝑛𝑖=1 𝐿𝑖
̅=
𝐾
𝐿
∑𝑛𝑖=1 ( )
𝐾 𝑖

Por otro lado, suponiendo sistema de flujo radial, la permeabilidad promedio es igual a:
𝑟
ln (𝑟𝑒 )
̅= 𝑤
𝐾 𝑟
ln (𝑟 𝑖 )
𝑖−1
∑𝑛𝑖=1 [ ]
𝑘𝑖

Suponiendo un radio del pozo (𝑟𝑤 ) de 0.25 ft, en la siguiente tabla se muestran los resultados
obtenidos:

Tabla10. Permeabilidad promedio: flujo lineal y radial.

Longitud Permeabilidad 𝑳 𝒓𝒊 𝐥𝐧(𝒓𝒊 𝐥𝐧(𝒓𝒊 /𝒓𝒊−𝟏 )

( )
(ft) (mD) 𝑲 𝒊 /𝒓𝒊−𝟏 ) 𝒌𝒊

150 10 15 150 6.3969 0.63969

200 20 10 350 0.8473 0.042365
300 30 10 650 0.6190 0.02063
500 50 10 1150 0.5705 0.01141
200 80 2.5 1350 0.1603 0.000200375
∑ 1350 ∑ 47.5 ∑ 0.7161

Luego, se reemplazan las sumatorias obtenidas en las fórmulas correspondientes a flujo radial y
lineal.

La permeabilidad promedio para flujo lineal es igual:

1350 𝑓𝑡
̅=
𝐾 ̅ = 28.421 𝑚𝐷
𝐾
47.5 𝑓𝑡/𝑚𝐷

31
Ahora la permeabilidad promedio para flujo radial es igual:

1350 𝑓𝑡
ln ( )
̅= 0.25 𝑓𝑡 ̅ = 12 𝑚𝐷
𝐾 𝐾
0.7161 𝑓𝑡/𝑚𝐷−1

32
 12. BIBLIOGRAFÍA

 Escobar Macualo, F. H. (2005). Fundamentos de ingenieria de yacimientos. Neiva:

Universidad Surcolombiana.
 Martínez Pérez, J. A. (2019). guias de laboratorio, análisis de núcleos. Neiva.
 martinez, j. (2019). talleres Analisis de nucleos. neiva: usco.
 Paris de ferrer, M. (2009). fundamentos de ingenieria de yacimientos.
 Parra , R. (2011). propiedades fisicas de los fluidos de yacimiento. Neiva: Universidad
Surcolombiana.

33