4.

Presión de Poro

4.1 Resumen

El siguiente capitulo, muestra los efectos de generación de presión de

formación anormal. Sus estudios muestran como influyen en la forma de
analizar las geopresiones. Se analizan los mecanismos que alteran las
presiones de formación normal por efectos de compactación (intervalos
someros) y por expansión de los fluidos (debido a zonas de baja
compresibilidad). La calibración de las líneas de tendencia de compactación
normal, se ajustan bajo estos criterios, definiendo los mejores match de tiempo
de transito y resistividad. Finalmente se presentan las metodologías para
determinar la presión de poro en función de la ecuación mas utilizada a nivel
mundial (Eaton. B.) así como los resultados obtenidos después de su
determinación.

4.2 Objetivo

Presentar una metodología para la determinación de la presión de poro de

acuerdo los conceptos fundamentales de presión de formación y línea de
compactación normal para tiempo de transito.

4.3 Introducción

Se ha escrito mucho referente a la predicción de la presión de la formación y el

gradiente de fractura. En la ingeniería de perforación, los métodos relacionados
a la determinación de la presión de poro anormal puede separarse en dos
categorías: métodos predictivos y métodos de verificación.

Los métodos de predicción utilizan la velocidad de una onda (sísmica) en los

intervalos de interés. Los métodos de detección utilizan parámetros físicos
como son el exponente “d” y pruebas de goteo.

Posiblemente no existen ecuaciones o métodos en los cuales se puedan

aplicar en todo el mundo, en todas ocasiones y bajo todas las condiciones,
porque todos los métodos están basados en un grado de variación a partir de
datos empíricos. Todos los métodos de predicción de presión de poro tienen la
base en el método origina de Therzagi: cuya ecuación es la siguiente:

SC  Ppn   v
Ecu. 4.1
En donde:

Pp = Presión de poro [Psi]

S = Presión de Sobrecarga [Psi]
σv = Esfuerzo vertical efectivo [Psi]

33
La presión de sobrecarga (S) es la presión por los pesos de los sedimentos
suprayacentes.

El esfuerzo vertical efectivo (σv) es el resultado del contacto de los granos, y

es una fuerza diamétricamente opuesta a la presión de sobrecarga resultado
de la matriz de la roca.

Los mejores métodos para la predicción de presión de poro, son los métodos
de Toman y Johnson, profundidad equivalente y Eaton.

Estos son métodos gráficos que utilizan los registros de densidad, velocidad y
valores eléctricos. La base esta fundamentada bajo la premisa de las
condiciones normales de valores de porosidad, que decrecen con la
profundidad a través de una tendencia normal.

Figura 4.1 Tendencia de Comparación normal con zonas anormales. (Adams, N, 1996)

En formaciones de presión anormal, la porosidad en el lugar tiene una

desviación de su línea normal. Durante los procesos de depositación la
sobrecarga incrementa tanto como los sedimentos son depositados. El
incremento de la sobrecarga compacta las formaciones, resultando en un
decremento de la porosidad. Los procesos de compactación ocurren mientras
el agua es expulsada de los espacios porosos y el esfuerzo de sobrecarga es
transferido de los fluidos de los poros a la matriz. Las presiones anormales
ocurren cuando a los fluidos de los poros no se les permite escapar. La
retención del agua, no permite la compactación, por lo tanto, la porosidad en el
lugar no decrece con la profundidad.

Las propiedades de las lutitas son medidas por varios registros y están
directamente relacionadas con la porosidad. Cuando estos valores son
graficados en función de la profundidad, la sección de presión normalmente
sigue una tendencia normal, mientras la porosidad en el lugar decrece. Una
desviación de los valores graficados de su tendencia normal es un indicador de
tendencia anormal. Los primeros métodos pueden utilizar también,

34
información sísmica para transformar la velocidad en intervalos de tiempo de
transito.
Las propiedades eléctricas y acústicas están directamente relacionadas con la
porosidad. Por eso una medición de las propiedades de las arcillas debe ser
representativa de la porosidad del lugar, esto es porque el incremento en la
presión de sobrecarga no se transfiere a la matriz, pero si a los fluidos de los
poros, lo que da como resultado una presión anormal.

En las partes someras normalmente la presión de la formación, registra valores

de resistividad que son mayores con respecto a los valores que puedan
presentarse con la variación de la profundidad dando como resultado una
tendencia normal. El incremento en los valores de resistividad de las arcillas es
resultado de los valores de porosidad del lugar, que decrecen por efectos de
compactación normal. En secciones de presión anormal hay una desviación de
los valores de resistividad de las arcillas con respecto a su tendencia normal.
La desviación es debido a los bajos valores de resistividad de la tendencia
normal debido a un incremento de porosidad en el lugar. Cuando utilizamos
registros acústicos, los valores de tiempo de transito decrecen con la
profundidad en la sección normalmente represionada del pozo. En secciones
anormales de presión, los valores de tiempo de transito se desvían a valores
mayores con respecto a la porosidad del lugar.

Por lo tanto, dado este contexto, es importante conocer las propiedades de las
lutitas en toda la extensión de sus propiedades, que como ya describimos
arriba son: resistividad, tiempo de transito y porosidad. Sin embargo, es
importante para este estudio, conocer los valores de densidad total de las
lutitas dada la definición de sobrecarga.

En este estudio es importante hablar de las propiedades de las arcillas y su

impacto en las mediciones de las propiedades de interés, hablando en términos
de los registros. Por esta razón, hablaremos de las arcillas y su importancia.

4.4 Conceptos generales.

Presión de Formación.

La presión de formación, también conocida como presión de poro (Pp), es

aquella a la que se encuentran confinados los fluidos en el espacio poroso de
la formación, estos fluidos pueden ser gas, aceite y/o agua salada.

Para entender las fuerzas responsables que provocan la presión de poro se

deben de considerar los eventos geológicos en la zona de interés. Un ejemplo
de los mas simples y mas comunes ocurre en los sedimentos someros que
fueron depositados lentamente en un ambiente de sedimentación deltaico, el
materia detrítico es transportado por el río hacia el mar, los sedimentos
formados inicialmente no son consolidados y tienen porosidades y

35
permeabilidades relativamente altas, el agua de mar mezclada con estos
sedimentos permanece en comunicación ejerciendo una presión hidrostática.

Una vez que la sedimentación ocurre, el peso de las partículas sólidas es

soportado por el punto de contacto grano a grano; así que, la presión
hidrostática del fluido contenido en los espacios poroso de los sedimentos
depende solamente de la densidad del fluido.

Presión de Formación Normal.

Cuando la presión de formación es aproximadamente igual a la presión

hidrostática teórica a la profundidad de interés, se considera que es una
presión de formación normal, esta presión regularmente es expresada como un
gradiente hidrostático. En general, la presión de formación normal es la presión
hidrostática ejercida por una columna de agua de 80,000 ppm de NaCl, cuya
densidad es de 1.074 gr/cm3 desde la superficie hasta la profundidad de
interés.

Presión de Formación Anormal.

Este término se utiliza para describir las presiones de formación que son
mayores que la normal. Este proceso es atribuido a dos procesos:

3.1 Al incremento del esfuerzo aplicado a una roca compresible.

3.2 La expansión de fluidos dentro de un espacio poroso limitado.

Ambos procesos son más efectivos en litologías de grano fino que grueso.

Compactación.

La relación entre la compactación del fluido se puede entender analizando un

ejemplo expuesto por Hubbert y Rubey tomando como referencia la teoría de
Terzaghi. En la figura 4.2, un resorte simula la comunicación entre las
partículas de arcilla y los pistones del sistema. Los piezómetros son utilizados
para registrar la presión del fluido. Mientras se aplica una presión sobre el
pistón superior, la longitud del resorte entre los pistones disminuye tanto como
el agua no escape del sistema. Por lo tanto, en la etapa inicial la presión
aplicada es soportada enteramente por la presión opuesta generada por el
agua.

36
 Figura 4.2 Relación de esfuerzo de sobrecarga, confinado y matricial.
Una manera de registrar esta presión en el resorte, es en términos de la razón
de la presión del fluido y de la presión total aplicada por el pistón superior.
Entonces este esfuerzo del resorte puede ser expresado por la siguiente
relación:

λ= p/S
Ecu. 4.2
En la misma figura 1, se observa que mientras no exista escape de agua del
sistema, la relación del resorte será de uno, λ=1, si existe una ligera expulsión
de agua del sistema, la relación del esfuerzo del resorte será menor que uno,
λ<1 y cuando el agua es expulsada totalmente del sistema, el agua residual y
el resorte estarán en total equilibrio. Este modelo se asemeja a la analogía de
las arcillas saturadas de agua y soportando una sobrecarga generada por los
estratos superiores de un punto de interés. Entonces dada esta analogía, se
puede definir la siguiente expresión:

σ=S-Pp
Ecu. 4.3
En zonas del terciario, existen zonas de lutitas que pueden contener
intercalaciones de arenas. Estas intercalaciones actúan como mecanismos que
permiten que el agua contenida en las arcillas, pueda circular a consecuencia
del valor de sobrecarga a la profundidad deseada. Finalmente, estos
fenómenos están relacionados con la porosidad y su compactación.

Mecanismos de generación de la presión de formación.

Glenn Bowers, utilizó la teoría para estudiar el comportamiento de

compactación de las lutitas a través de la relación del esfuerzo efectivo y la
velocidad. Durante su estudio determinó que la velocidad del sonido en lutitas
sistemáticamente dependía del criterio que genera la presión.

El esfuerzo efectivo continuamente incrementa con respecto a la profundidad,

por lo tanto el esfuerzo efectivo versus la velocidad, tendrían una tendencia
normal.

37
Los dos mecanismos considerados para estimar la presión de formación son
por efectos de baja compactación y por efectos de expansión.

Baja compactación.

La sobrepresión ocurre cuando los fluidos de la formación están entrampados

por la baja permeabilidad, lo que genera que el agua escape o sea expulsada
por el peso de los sedimentos nuevos. Este efecto se conoce como
desequilibrio de compactación.

La baja compactación no puede causar decremento al esfuerzo efectivo, por lo

tanto, la curva virgen también aplica para formaciones sobrepresionadas por
baja compactación. Esto generará que el esfuerzo efectivo sea constante con
respecto al tiempo, lo que generará que la velocidad sea constante en la zona
virgen.

Efectos de expansión

El calor, maduración de los hidrocarburos, migración, expulsión y expansión del

agua debido a la diagénesis de las arcillas, son mecanismos de generación de
presión por expansión.

El exceso de presión es resultado de la contracción de la matriz de la roca,

dado que el volumen del fluido intenta incrementar. Esto implica que el
esfuerzo efectivo disminuya mientras la depositación continua sepultando ese
sedimento, provocando así, un cambio en la velocidad del sonido con respecto
a la zona normal de compactación.

No existe una teoría que determine si el cambio de la velocidad es por la

subcompactación o expansión. Lo que si es seguro, es que la sobrepresión es
consecuencia de una relación de la compresibilidad del fluido y la matriz.
Ambos actúan como un resorte en paralelo. Si la compresibilidad de la matriz
es mayor que el esfuerzo, entonces el fluido carga o actúa contra la
sobrecarga, por otro lado si la compresibilidad de los fluidos es mayor que la
roca, entonces el esfuerzo matricial actuará contra la sobrecarga en mayor
proporción.

Por lo tanto, la subcompactación se considera típica para generar sobrepresión

a profundidades someras, en donde la formación es muy suave.

4.5 Metodología

38
La metodología utilizada para la determinación de la presión de poro es a
través del método de Glenn Bowers. Para este caso, el programa integra una
serie de parámetros A, B y Vo, en donde cada término son factores de ajuste
determinados a través de una serie de datos representativos dos zonas, la
Zona Virgen y la Zona de Descarga.

Este método es únicamente de ajuste para el tiempo de transito de los

registros. La mejora en este método es que se considera la generación de la
presión por efectos de compactación y en segundo termino por efectos de
expansión de los fluidos.

La metodología que se utiliza es que, después de la selección de los intervalos

de interés y el promediado de estos, los datos de tiempo de transito son
utilizados para:

1) De los datos de tiempo de transito, se determinan las velocidades.

2) Posteriormente de los datos de densidad volumetrica, se determina la
densidad de la matriz.
3) Posteriormente se determina el esfuerzo de efectivo en unidades de
[Psi].
4) Se utilizan los parámetros de ajuste A, B y Vo para determinar una
velocidad normal Vn.

Para la curva virgen:

B
Vn  5000  Aσ v Ecu. 4.4

En donde
V=[Pies/seg]
σv= [Lb/Pulg^2]
A y B son parámetros que se muestran en el apéndice A.

Para Curvas de carga:

1
 1 B
  σz U 
V  5000  A σ zmax *  
σ 
  zmax  
 
Ecu. 4.5

En donde:
σv =((Vmax-5000)/A)1/B

5) Una vez determinados los parámetros de ajuste, se grafican las líneas

de tendencia de compactación normal de Tiempo de Transito vs
Profundidad, tanto para los datos del registro como para los de ajuste,
como en la figura 4.3.

39
 VELOCIDAD VS PROFUNDIDAD

Velocidad [Mpies/seg]

0 2 4 6 8 10 12 14 16 18 20
0

1000

2000

Profundidad [m]
3000

4000

5000

6000

Fig. 4.3 Línea de Compactación Normal

6) Posteriormente se utilizan los datos de la velocidad para regresarlos a

tiempo de transito, bajo la premisa de que la velocidad sónica es el
inverso del tiempo de transito. Se obtienen los nuevos valores para un
Tiempo de Transito normal y un Tiempo de Transito observado. Ver
tabla 4.1.

TIEMPO DE TRANSITO

Tiempo de Transito [mseg-pie]

10 100 1000
0

500

1000

1500
PROFUNDIDAD [PIES]

2000

2500

3000

3500

4000

4500

5000

TIEMP O DE TRANSITO
LCN "Bo wers " Tiempo deTransito

Figura 4.4 Tiempo de Transito Normal VS Tiempo de Transito Observado.

40
 Tabla 4.1 Resultados obtenidos después de la selección de los intervalos, promediado y ajuste con los
parametroa A,B y Vo de Bowers.

4.5.1 Cálculo de Presión de Poro

Algunas investigaciones realizadas por Mouchet y Mitchell6 en 1989 y por Yoshida7 y

otros en 1996, concluyeron que las ecuaciones de Eaton para la predicción de la
presión de formación son las más ampliamente utilizadas para calcular gradientes de
presión de formación en aguas profundas. Dichas ecuaciones son:

Para la Resistividad

1.2
P S S P   Ro 
     Ecu. 4.6
D D  D DN   RN 

Donde:
P  psi 
Gradiente de presión de formación,  
D  pie 

S  psi 
Gradiente de presión de sobrecarga,  
D  pie 

P  psi 
Gradiente de presión normal de formación,  
DN  pie 

Ro Resistividad de la lutita, de la lectura del registro,  ohm  m 

R N Resistividad de la lutita, de la línea de tendencia normal de compactación,  ohm  m 

Para la conductividad:

41
 1 .2
P S S P C N  Ecu. 4.7
    
D D  D DN  Co 

Donde:
 miliohms 
CN Conductividad de la lutita, de la línea de tendencia normal de compactación,  
 m 

 miliohms 
C o Conductividad de la lutita, de la lectura del registro,  
 m 

Si se utiliza el registro sónico de porosidad o los tiempos de tránsito del registro

sísmico, también se puede usar la ecuación:

3
P S S P   t n 
     Ecu. 4.8
D D  D D N   t o 

Donde:

 seg 
t n Tiempo de tránsito de la lutita, de la línea de tendencia normal de compactación,  pie 
 

t o Tiempo de tránsito de la lutita, de la lectura del registro  pie

seg 

 

Hottman y Johnson

Para este método se desarrollo la ecuación que mejor ajustara a partir de una
grafica en donde se utilizaron valores de resistividad. Este método en el eje de
las x, relaciona gradientes de todas las zonas del golfo de México para cada
valor de resistividad. Por lo tanto, es un método estadístico más que de
aplicación numérica.

La grafica se presenta en la figura siguiente, en donde los valores de Y,

representan Kg/cm3/m. Finalmente después de utilizar los método de LCN, se
aplican estas metodologías en las ventanas de operación.

42
4.5.2 Gráfica de la Presión de Poro

En la figura 4.5 se puede observar el resultado de la presión de poro obtenida

después de calibrar el tiempo de transito con los valores en los coeficientes:
A=2.6740, B=0.8617 y Vo=7200.

Posteriormente se muestran los valores para el tiempo de transito normal y

observado, (mostrados en la tabla 2.2) junto con los valores del tiempo de
transito se muestran también los de resistividad normal y observada.

En la tabla 2.3, se muestran los valores obtenidos para la presión de poro y

presión de fractura.

VENTANA OPERACIONAL SAMA 1

Presión de Poro

0.50 0.80 1.10 1.40 1.70 2.00 2.30

0.00
200.00
400.00
600.00
800.00
1000.00
1200.00
1400.00
1600.00
1800.00
2000.00
2200.00 Presión de Poro [gr/cm3]
Profundidad

2400.00 Presión Poro F.S [gr/cm3]

2600.00 Presión de Fractura [gr/cm3]
2800.00 Presión Fractura F.S. [gr/cm2]
3000.00
Esfuerzo de Sobrecarga [gr/cm3]
3200.00
EFECTO DE BROTE
3400.00
3600.00
3800.00
4000.00
4200.00
4400.00
4600.00
4800.00
5000.00
5200.00

Figura 4.5 Ventana Operacional, Calculo de Presión de Poro.

43
 Tabla 2.3 La columna morada muestra los valores de presión de poro sin factor de seguridad y
con factor de seguridad por efectos de suaveo y surgencia.

4.6 Conclusión

La determinación de geopresiones a través de promediado de datos es un

método mas simple y exacto.

La evaluación de geopresiones a través de la separación de la zona virgen

y la zona de descarga, permite obtener menor error de calculo para este
método.

La selección de intervalos es fundamental para obtener una ventana de

operación mejor definida. Sobre todo en aguas profundas, en donde el
efecto de la sobrecarga influye en que esta ventana sea mas estrecha,
podemos hablar de que este método con un tirante de 900 mts, permite el
asentamiento con mayor forma de operación.

4.7 Referencia

44
1) Mouchet J.P., A Mitchell. / ABNORMAL PRESSURES WHILE DRILLING Elf Aquitaine
Manuels Techniques 1989.

2) Yoshida, C. Ikeda, S. Y Eaton B.A: “AN INVESTIGATIVE STUDY OF RECENT

TECHNOLOGIES FOR PREDICTION, AND EVALUATION OF ABNORMAL
FORMATION PRESSURE AND FRACTURE PRESSURE IN NORTH AND SOUTH
AMERICA” IADC/SPE 36381,1996

3) Hottman, C.E, Jonson R.K. ESTIMATION OF FORMATION PRESSURE FROM LOG –

DERIVED SHALE PROPERTIES JPT, Agosto 1965.

4) Pennebaker, E.S./ AN ENGINNERING INTERPRETATION OF SEISMIC DATA SPE

2165

5) Terzaghi K, THEORICAL SOIL MECHANICS John Wiley & Sons, Inc. New York 1943

6) Gardner,G.H.F., Gardner L.W. / FORMATION VELOCITY AND DENSITY ,

THE DIAGNOSTIC BASIS FOR STRATIGRAPHICS TRAPS Geophysics Vol
39 Num 39.

7) Barker, J.W., Wood T.D./ ESTIMATING SHALLOW BELOW MUDLINE DEEPWATER GULF OF
MÉXICO FRACTURE GRADIENTS Annual Technical Forum Abril 2 – 3 1997

45