Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-1

Capítulo 13

Información Geomecánica a
Partir de Registros Acústicos

General

En este capítulo se presentan las bases teóricas para la obtención de parámetros

geomecánicos a partir de los registros petrofísicos. Los registros petrofísicos están
compuestos por sondas que tienen varios emisores y receptores de ondas
acústicas que permiten evaluar el tiempo de tránsito de las ondas P y S en la
formación. Esto proviene de la teoría tradicional de propagación de ondas en
medios elásticos continuos y las modificaciones que se hace de ella utilizando
materiales poro-elásticos. Los módulos elásticos de la formación equivalen a
valores dinámicos debido al elevado intervalo de frecuencia al que trabajan los
emisores de ondas; sin embargo, el interés de la geomecánica generalmente está
asociado a problemas que son de naturaleza cuasi-estática. Afortunadamente, es
posible corregir los valores de los módulos dinámicos obtenidos a través de
registros acústicos mediante la calibración con ensayos geomecánicos en núcleos.
Estos resultados se presentan antes de la culminación del estudio para que se
tomen medidas correctivas lo más pronto posible.

Teoría de propagación de ondas en medios elásticos continuos

La determinación de las ecuaciones de movimiento para un medio elástico

comienza con el análisis de equilibrio de un pequeño elemento; esta derivación ha
sido desarrollada por muchos autores (Timoshenko y Goodier, 1951; Kolsky, 1963;
Richart y otros, 1970) y por lo tanto sólo se presentan los pasos más importantes.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez H.

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-2

En primer lugar se tiene que considerar la variación de esfuerzos en cada cara del
elemento tal como se muestra en la Figura 13.1. El equilibrio del elemento puede
ser expresado escribiendo la sumatoria de fuerzas que actúan paralelamente a
cada eje. Por ejemplo, la ecuación de equilibrio para la dirección x es la siguiente:

Figura 13. 1.- Sumatoria de esfuerzos normales y de corte de un elemento de material elástico.

 ∂σ x   ∂τ xy   ∂τ 
 σx + ∆x ∆y∆z − σ x ∆y∆z +  τ xy + ∆y ∆x∆z − τ xy ∆x∆z +  τ xz + xz ∆z ∆x∆y − τ xz∆x∆y = 0
 ∂x   ∂y   ∂z 

para la dirección y, se tiene:

 ∂σ y   ∂τ yx   ∂τ yz 
σy + ∆y ∆x∆z − σ y ∆x∆z +  τ yx + ∆x ∆y∆z − τ yx ∆y∆z +  τ yz + ∆z ∆y∆x − τ yz ∆y∆x = 0
 ∂y   ∂x   ∂z 

y para la dirección z, se tiene:

 ∂σ   ∂τ   ∂τ zy 
 σ z + z ∆z ∆x∆y − σ z∆x∆y +  τ zx + zx ∆x ∆z∆y − τ zx ∆z∆y +  τ zy + ∆y ∆z∆x − τ zy ∆z∆x = 0
 ∂z   ∂x   ∂y 
(10. 1)

Suponiendo que las fuerzas del cuerpo son despreciables y aplicando la segunda
ley de Newton en la dirección x se llega a:

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-3

 ∂σ x ∂τ xy ∂τ xz  ∂ 2u
 + +  ∆x∆y∆z = ρ( ∆x∆y∆z) 2 (10. 2)
 ∂x ∂y ∂z  ∂t

para la dirección y se tiene que:

 ∂σ y ∂τ yx ∂τ yz  ∂ 2v
 + +  ∆x∆y∆z = ρ( ∆x∆y∆z) 2 (10. 3)
 ∂y ∂x ∂z  ∂t

y para la dirección z se tiene que:

 ∂σ z ∂τ zy ∂τ z x  ∂ 2w
 + +  ∆x∆y∆z = ρ( ∆x∆y∆z) 2 (10. 4)
 ∂z ∂y ∂x  ∂t

Donde u, v, y w son los desplazamientos en las direcciones x, y y z

respectivamente. Las ecuaciones de movimiento en las tres direcciones pueden
entonces ser escritas de la siguiente manera:

∂ 2 u ∂σ x ∂τ xy ∂τ xz
ρ 2 = + + (10. 5)
∂t ∂x ∂y ∂z

∂ 2 v ∂τ yx ∂σ y ∂τ yz
ρ = + + (10. 6)
∂t 2 ∂x ∂y ∂z

∂ 2 w ∂τ zx ∂τ zy ∂σ z
ρ 2 = + + (10. 7)
∂t ∂x ∂y ∂z

Utilizando las siguientes relaciones para un medio elástico, se tiene:

σ x = λε vol + 2Gε x (10. 8)

τ xy = τ yx = G γ xy (10. 9)

σ y = λε vol + 2Gε x (10. 10)

τ yz = τ zy = Gγ yz (10. 12)

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-4

σ z = λε vol + 2Gε z (10. 12)

τ zx = τ xz = Gγ zx (10. 12)

y donde:

E
G=
2(1 + ν )
(10. 14)

νE
λ=
(1 + ν )(1 − 2ν )
(10. 15)

donde ν es la relación de Poisson, G es el módulo de corte o módulo de rigidez, λ

es la constante de Lamé y ε vol es la deformación volumétrica o dilatación dada por
ε vol = ε x + ε y + ε z . También se necesitan las siguientes relaciones para la
deformación y rotación en función de los desplazamientos (ω es la rotación en
torno a cada eje):

∂u ∂v ∂u ∂w ∂v
εx = γ xy = + 2ϖ x = −
∂x ∂x ∂y ∂y ∂z

∂v ∂w ∂v ∂u ∂w
εy = γ yz = + 2ϖ y = − (10. 16)
∂y ∂y ∂z ∂z ∂x

∂w ∂u ∂w ∂v ∂u
εz = γ zx = + 2ϖ z = −
∂z ∂z ∂x ∂x ∂y

Cuando estas relaciones son combinadas con las ecuaciones de movimiento,

entonces se tiene:

∂ 2u ∂ε vol  ∂ 2u ∂ 2u ∂ 2u
ρ 2 = (λ + G) + G 2 + 2 + 2  (10. 17)
∂t ∂x  ∂x ∂y ∂z 

∂ 2v ∂ε vol  ∂ 2 v ∂ 2v ∂ 2v 
ρ 2 = (λ + G) + G 2 + +  (10. 18)
∂t ∂y  ∂x ∂y 2 ∂z 2 

∂ 2w ∂ε vol  ∂ 2 w ∂ 2w ∂ 2w 
ρ 2 = ( λ + G) + G 2 + +  (10. 19)
∂t ∂z  ∂x ∂y 2 ∂z 2 

que son las ecuaciones de movimiento para un medio elástico, isótropo e infinito
homogéneo.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-5

Dos soluciones son posibles para esta ecuación de movimiento. La primera

solución describe la propagación de una onda con pura deformación volumétrica
(sin rotación, pero con cambio de volumen) mientras que la otra describe la
propagación de una onda de pura deformación rotacional (con rotación, pero sin
cambio de volumen).

La primera solución se obtiene diferenciando cada una de ellas con respecto a x, y

y z, respectivamente; sumando luego las expresiones se obtiene:

∂ 2 ε vol
ρ = (λ + 2G)∇ 2ε vol (10. 20)
∂t 2

donde ∇ 2 es el operador Laplaciano en coordenadas cartesianas que se define

como:

∂2 ∂2 ∂2
∇ = 2
+ +
∂x 2 ∂y 2 ∂z 2

La ecuación 10.20 tiene la forma de la ecuación de onda, si se expresa de la

siguiente manera

∂ 2 ε vol
ρ = v p2∇ 2 ε vol (10. 21)
∂t 2

donde vp es la velocidad de propagación de la onda de deformación volumétrica,

la cual es conocida por onda primaria, onda compresional o simplemente onda P y
viene dada por:

λ + 2G
vp = (10. 22)
ρ

La otra solución de la ecuación de movimiento puede ser obtenida diferenciando

las ecuaciones 10.18 y 10.19 con respecto de z y y, respectivamente, y eliminando
la deformación volumétrica por substracción; de esta operación resulta la siguiente
expresión:

∂ 2  ∂w ∂v   ∂w ∂v 
ρ 2 
−  = G∇ 2  −  (10. 23)
∂t  ∂y ∂z   ∂y ∂z 

si usamos la expresión para la rotación de ϖ x tenemos que:

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-6

∂2ϖx
ρ 2 = G∇2ϖx (10. 24)
∂t

donde esta expresión tiene la forma de la ecuación de onda si se expresa de la

siguiente manera:

∂ 2ϖ x
= v 2s∇2ϖx (10. 25)
∂t 2

Expresiones similares pueden ser obtenidas para ϖ y y ϖ z , lo cual implica que la

relación es propagada con velocidad:

G
vs = (10. 26)
ρ
donde Vs representa la velocidad de propagación de una onda de pura
deformación rotacional (distorsión sin compresión), la cual es conocida como onda
de corte, onda secundaria o simplemente onda S .

Este análisis demuestra que un medio infinito, elástico e isótropo puede sostener
solamente dos tipos de ondas las ondas P y S. Las velocidades de onda P para
barras es diferente (menos rápida) que la velocidad en medios infinitos; mientras
que la velocidad de corte en barras es igual a la velocidad en medios infinitos. Es
por esto que muchos métodos experimentales de laboratorio utilizan la onda de
corte para determinación de parámetros en el campo.

Si el medio no es infinito sino que contiene una frontera, entonces es posible que
existan otros tipos de ondas tales como las ondas Rayleigh y las ondas Love. En
geometrías complejas tales como cilindros llenos de fluidos en medios infinitos
(hoyos de perforación), es posible que existan otros tipos de ondas conocidas
como ondas de fluido (mud wave) y de tubo (tube wave), entre las cuales destacan
las ondas Stoneley.
En estos momentos las aplicaciones geomecánicas de ondas, trabajan
principalmente con las ondas compresionales y de corte debido a lo sencillo del
modelo infinito elástico e isótropo. Sin embargo, muchos centros de investigación
están desarrollando teorías que trabajan con otros tipos de ondas.

Propagación de ondas en perforaciones petroleras

Las herramientas acústicas están compuestas por uno o varios receptores de

ondas acústicas colocados a una distancia de uno o varios emisores de ondas (ver
Figura 13.2). Estas herramientas tienen circuitos electrónicos que permiten medir
tiempos de tránsito de las ondas a medida que se mueve la herramienta en el
pozo, y conociendo la distancia viajada por la onda es posible conocer la velocidad
de propagación de las ondas.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-7

Figura 13. 2.- Ejemplo de una herramienta petrofísica para medir ondas acústicas.

Se simulará el sistema real compuesto por una perforación petrolera llena de fluido
de perforación, perforada en una formación geológica de interés por medio de un
sistema idealizado, compuesto por un cilindro lleno de fluido situado en un medio
elástico y homogéneo (ver Figura 13.2). Un emisor de ondas hipotético dentro de
este cilindro lleno de fluido, sólo puede emitir ondas compresionales ya que el
fluido no permite ondas de corte; sin embargo estas ondas compresionales van a
chocar con las paredes del medio elástico sufriendo refracciones que pueden
crear otros tipos de ondas entre las cuales destacan las S por su interés
geomecánico.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-8

Figura 13.3.- Modos de propagación de ondas en un cilindro dentro de un medio

elástico, isótropo y homogéneo que simula una perforación petrolera en una
formación geológica.

Tal como se observa en la Figura 13.4, en la pared del pozo ocurre una interfase,
donde se sabe por física básica que las ondas incidentes crearán ondas reflejadas
y ondas refractadas. Cuando un rayo acústico en un medio M 1 de velocidad de
propagación de onda V1, choca oblicuamente con la frontera de otro medio M 2 de
velocidad de propagación V2; parte de la energía acústica es reflejada al medio M 1
como un rayo reflejado R y otra parte es transmitido dentro del medio M 2 como un
rayo refractado T. El ángulo de incidencia α y el de refracción θ están
relacionados por medio de:

senα senθ
=
V1 V2

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-9

Medio 1 Medio 2
I

Angulo incidencia

Angulo refraccion
Angulo refleccion
T

R
V2 V1

Figura 13. 4.- Esquema de rayos acústicos chocando en una interfase creando rayos reflejados y
refractados.

Medio 1 Medio 2
I

Angulo incidencia
Angulo refraccion
Angulo refleccion
Onda S
refractada
Onda P
reflejada Onda P
refractada

V1 V 2s , V2P

Figura 13. 5.- Esquema de los ángulos de incidencia de la onda compresional en

el fluido de perforación y las ondas compresional y onda de corte creadas por
refracción en un medio elástico homogéneo que es la formación geológica de la
pared del pozo.

Como se explicó anteriormente, un emisor acústico en un fluido como el que se

encuentra en la perforación, sólo puede crear ondas compresionales ya que las
ondas de corte no pueden existir en fluidos. Como se puede observar en la Figura
13.5, esta onda compresional incidente va a crear sólamente una onda
Introducción a la Geomecánica Petrolera Andrés R. Vásquez
 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-10

compresional reflejada, y no una onda de corte, ya que la propagación sigue

ocurriendo en el fluido; pero en la formación que es un sólido esta onda incidente
creará una onda compresional refractada y una onda de corte refractada las cuales
se propagarán en la formación para ser luego detectadas por el receptor.

A medida que el ángulo de incidencia de la onda compresional aumenta, se llegará

a un ángulo crítico, donde el ángulo de refracción de la onda de corte se vuelve de
90°, lo que implica que la onda S se propaga paralela a la pared del pozo. Esto
también implica que una onda propagándose por la formación paralela a la pared
del hoyo, radiará energía hacia el fluido del hoyo a este mismo ángulo crítico.

Figura 13.6.- Propagación de ondas de una herramienta acústica monopolar en

una formación consolidada con alta rigidez.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-11

Figura 13.7.- Propagación de ondas de una herramienta acústica monopolar en

una formación no consolidada con baja rigidez.

Herramientas petrofísicas de mediciones acústicas

Las herramientas petrofísicas de mediciones acústicas han existido desde hace

muchos años, sin embargo, éstas sólo medían el tiempo de llegada de la onda P.
La primera generación de estas herramientas tenían un emisor y un receptor a una
distancia fija, por lo que el tiempo de tránsito permitía calcular la velocidad de
propagación en unidades de distancia sobre tiempo (o la rapidez, en unidades de
tiempo/longitud). Para corregir problemas de centralización y de calidad de hoyos
se modificaron estas herramientas para incluir varios receptores y emisores, por lo
que esta segunda generación de herramientas se conocieron como sónicos
compensados. Estas eran herramientas monopolares (sin polarización), con
lecturas analógicas donde se escogía el tiempo de llegada de la onda P viendo la
forma de la onda grabada, las cuales estuvieron en el mercado por mucho tiempo,
hasta que a mediados de los años 80 apareció la tercera generación de
herramientas acústicas.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-12

Esta tercera generación de herramientas acústicas fueron los sónicos digitales

monopolares que permiten procesar el tren de ondas completo, compuesto por las
ondas P, S y Stoneley y por lo tanto representan la primeras herramientas que
permitieron obtener parámetros geomecánicos a partir del registro acústico. Estas
herramientas permiten medir los tiempos de tránsito de las ondas P y S los cuales
permiten calcular las velocidades de propagación de dichas ondas para elaborar
un registro completo del pozo. El conocimiento de las velocidades de propagación
(o el tiempo de tránsito) de las ondas P y S permite calcular los módulos elásticos
dinámicos de una manera directa, utilizando la ecuación de onda y la teoría de
elasticidad.

Las herramientas monopolares permiten que el tiempo de llegada de las ondas P

sea fácil de determinar para casi todo tipo de formación consolidada ya que esta
es rigida (ver Figura 13.6). Sin embargo, el tiempo de llegada de la onda S resulta
mucho más complicado, especialmente si la formación es de una rigidez que se
aproxime al tiempo de tránsito de las ondas en el lodo del hoyo. Esta condición se
presenta en formaciones de arenas no consolidadas y por lo tanto se tuvo que
diseñar una cuarta generación de herramientas que funcionaran con ondas
polarizadas en forma dipolar.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-13

Las herramientas monopolares contienen un emisor omnidireccional de presión, el

cual crea un pulso de onda compresional en el fluido de perforación, que a su vez
se propaga hacia la formación. Cuando este pulso llega a la pared del pozo, esta
perturbación radial a su vez excita ondas compresionales y de corte que se
propagan hacia dentro de la formación. A medida que estas ondas
compresionales y de corte se propagan en la formación y de regreso al pozo, las
mismas inducen un tren de ondas compresionales en el fluido, conocidas como
ondas de cabeceras (head waves), las cuales son las únicas ondas medidas por
los receptores de la herramienta. Las ondas compresionales en la formación
siempre son más rápidas que las ondas a través del fluido del hoyo; sin embargo,
las ondas de corte en algunas formaciones de poca rigidez pudieran ser más
lentas que las ondas del fluido del hoyo. Las ondas de cabecera son creadas
solamente cuando la velocidad de las ondas en la formación son más rápidas que
las del fluido y por lo tanto la onda compresional siempre puede ser medida por
herramientas monopolares. Luego de las ondas de cabecera, llegan las ondas de
tubo y las ondas Stoneley, las cuales son ondas dispersivas (su velocidad es
dependiente de la frecuencia). La onda de tubo es producida por las reflexiones de
las ondas emitidas reverberando en el hoyo. Las ondas Stoneley son ondas de
superficie que se forman alrededor del hoyo y que viajan más despacio que las
ondas del fluido del hoyo. Tal como se presenta en la Figura 13.7; no es entonces
posible medir las ondas de corte en formaciones de poca rigidez (formaciones
lentas).

Figura 13.8.- Principio de operación de un transmisor dipolar en las herramientas acústicas.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-14

Las herramientas acústicas dipolares fueron creadas para medir los tiempos de
tránsito de la onda de corte en formaciones de poca rigidez. Esta herramienta
utiliza emisores direccionales en forma de pistones que crean un aumento de
presión en un lado del pozo y una disminución de presión en el otro lado. Este
cambio de presión causa una perturbación de tipo dipolar (por ello su nombre), la
cual crea una pequeña flexión de las paredes del hoyo tal como se muestra en la
Figura 13.8. Esta pequeña flexión del hoyo se propaga como una onda flexural que
es dispersiva, ya que a bajas frecuencias tiene la misma velocidad de la onda de
corte y a frecuencias más altas tiene valores de velocidad mayores. Las
herramientas dipolares pueden medir estas ondas flexurales hasta en las
formaciones menos rígidas y por lo tanto su uso en este tipo de formación se hace
obligatorio.

Cálculo de los módulos elásticos dinámicos

En la derivación de la ecuación de onda para medios elásticos continuos se pudo

observar que las velocidades de propagación de las ondas P y S son función de los
módulos elásticos. Esto implica que si se conocen los tiempos de tránsito o las
velocidades de propagación de las ondas se pueden determinar los módulos
elásticos.
En materiales poroelásticos también se pueden determinar módulos elásticos por
medio de tiempos de tránsito de las ondas. Sin embargo, los módulos elásticos
calculados de esta manera (mediciones dinámicas), dan resultados diferentes de
aquellos determinados por ensayos triaxiales (mediciones cuasi-estáticas o
estáticas). Desafortunadamente las aplicaciones geomecánicas en la ingeniería de
petróleo son de procesos estáticos o cuasi-estáticos los cuales sólamente pueden
ser medidos en el laboratorio. Esto implica que cualquier módulo elástico
determinado por tiempos de tránsito de ondas serán valores dinámicos que no
pueden ser utilizados directamente en aplicaciones.

Los registros petrofísicos acústicos miden el tiempo de tránsito de las ondas

compresionales y de las ondas de corte entre un emisor y un receptor a una
distancia determinada. Utilizando las velocidades de propagación de la onda
compresional Vp y la onda de corte Vs se pueden entonces definir los módulos
elásticos dinámicos de la siguiente manera:

Modulo de corte:

G = ρVs2 (10. 27)

Relación de Poisson:

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-15

Vp2 − 2 Vs2
ν=
( )
(10. 28)
2 Vp2 − Vs2

Constante de Lamé:

λ = ρVp2 − 2ρVs2 (10. 29)

Módulo de compresibilidad

4
K = ρVp2 − ρVs2 (10. 30)
3

Módulo de Young

E=
(
ρVs2 3Vp2 − 4 Vs2 )
(V )
(10. 31)
2
p − Vs2

Tradicionalmente, los valores medidos por las compañías de servicio vienen dados
en unidades de rapidez, que se expresa generalmente en unidades de tiempo
sobre longitud y por lo tanto es el inverso de la velocidad. Las unidades
tradicionalmente usadas por las compañías de servicio son el microsegundo por
pié (µseg/pié).

Para ecuaciones de trabajo se utilizan directamente las unidades de los registros y

las relaciones de la teoría de elasticidad, se pueden entonces expresar los
módulos de la siguiente manera

Relación de Poisson:

2
1  ∆t s 
  −1
2  ∆t p 
ν= 2
(10. 32)
 ∆t s 
  − 1
 ∆t p 

Módulo de corte (libras/pulg2)

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-16

ρb
G=ρ ×a (10. 33)
∆t s

Módulo de Young (psi)

E = 2G ( 1 + ν ) (10. 34)

Módulo volumétrico (psi)

 1 4 
K b = ρb  2 − 2
 ×a (10. 35)
 ∆t p 3∆t s 

Compresibilidad

1
Cb = (10. 36)
Kb

donde a = 1.34 x 1010 si la densidad ρ d está expresada en gr/cc y los tiempos de

tránsito de la ondas ∆t están expresadas en µseg/pié (convención de las compañías
de servicio).

Módulos elásticos estáticos y dinámicos

Tal y como se presentó anteriormente, se pueden calcular los módulos elásticos a

partir de velocidades de propagación de ondas; y debido a que estas ondas son
generadas a frecuencias altas, a los módulos calculados se les conoce como
módulos elásticos dinámicos. Los resultados obtenidos en estas ecuaciones son
generalmente distintos de los módulos medidos en los ensayos mecánicos de
laboratorio (módulos estáticos), donde se utiliza la relación entre esfuerzos axiales
σ x y las deformaciones axiales ε x y laterales ε y tal como se muestra a continuación:

σx
Es =
εx

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-17

εy
νs = −
εx

donde la tasa de carga oscila alrededor de 103 segundos. Este tipo de carga
monotónica en el laboratorio es cuasi-estática, sin embargo, tradicionalmente se
conoce como estática.

Aunque las investigaciones análíticas y experimentales no han permitido entender

completamente el fenómeno; sí está establecido que los módulos dinámicos son
mayores o iguales a los módulos estáticos. El módulo de Young dinámico puede
ser hasta 4 veces mayor que el módulo de Young estático y la relación de Poisson
dinámica puede ser hasta 3 veces mayor que la reación de Poisson estática;
dependiendo del material.

El único lugar donde es posible medir módulos elásticos estáticos y dinámicos es

el laboratorio. Esto se logra sometiendo una muestra del material de formación a
un ensayo triaxial donde el equipo esté instrumentado para realizar mediciones
dinámicas. Durante el proceso de compresión triaxial se obtendrá la curva
esfuerzo-deformación de donde se pueden calcular los módulos elásticos estáticos
y haciendo mediciones de propagación de ondas en diferentes instantes durante el
proceso de compresión, se podrán calcular los módulos elásticos dinámicos. Esto
quiere decir que es fundamental tener ensayos de laboratorio para poder
determinar experimentalmente la relación entre módulos estáticos y dinámicos
para crear correlaciones predictivas. Por lo anteriormente expuesto, resulta
importante utilizar módulos estáticos o dinámicos cuando la situación lo amerite.

Existen varias razones para explicar la diferencia entre módulos elásticos

dinámicos y módulos dinámicos estáticos para materiales poroelásticos. Entre las
mas importantes se tienen:

a) Diferencias en frecuencias y tiempos de carga: El intervalo de frecuencia

(longitud de banda) de los emisores de las primeras herramientas sónicas para
medir el tren de ondas completa estaban entre 5 y 18 kHz con las frecuencias de
más energía cercano a los 12 kHz; es decir, 12000 ciclos por segundo, lo que
implica que el tiempo de carga está en el orden de 10-4 segundos. Los últimos
avances en electrónica y diseño de transductores han causado que las
herramientas dipolares más recientes ahora hayan reducido la frecuencia con
mayor energía a valores cercanos a los 3 kHz. Aunque esta reducción de
frecuencias es un logro importante, los problemas geomecánicos de importancia,
arenamiento, fracturamiento hidraúlico, estabilidad de hoyo, compactación, etc.,
son fenómenos cuasiestáticos que se pudieran asignar conservadoramente con
tiempo de carga de 103 segundos. Como puede observarse hay más de 7 órdenes
de magnitud entre la velocidad de carga en ensayos dinámicos y en ensayos
estáticos y como es de esperarse carga en materiales poroelásticos causará

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-18

diferencias entre los módulos elásticos estáticos y los módulos elásticos

dinámicos.

b) Diferencias en el intervalo de deformaciones: Es conocido que los módulos

dinámicos de materiales porosos pueden sufrir diferencias debido al nivel de
deformación al que es sometido. En particular, el módulo dinámico no cambia si
las deformaciones cíclicas se mantienen por debajo de un umbral mínimo que en
arenas no consolidadas se encuentra cercano a 10-4 (Dobry y Swiger, 1979; Dobry
y otros, 1982; Vásquez, 1988; Vásquez y Dobry, 1989). Una vez que este umbral
es excedido, empieza a ocurrir una degradación del valor de los módulos que
creará una diferencia entre valores dinámicos y estáticos. Esta degradación de los
módulos en función de la deformación puede llegar a ser muy importante (Seed y
Idriss, 1970; Seed y otros, 1984; Sun y otros, 1988).

Debido a que existen diferencias entre el cálculo de módulos elásticos por medio
de mediciones dinámicas a partir de propagación de ondas y métodos estáticos
basados en la curva esfuerzo deformación; se debe utilizar algún tipo de
procedimiento para corregir los valores dinámicos. Un procedimiento usado para
generar valores continuos de los módulos estáticos (generalmente los parámetros
módulo de Young y la relación de Poisson) a partir de los módulos elásticos
dinámicos de los registros es el siguiente:

1. Realizar una corrida de un perfil sónico en el pozo para determinar el módulo de

Young y la relación de Poisson dinámica a nivel de campo.

2. Realizar mediciones dinámicas en núcleos en el laboratorio del módulo de

Young y de la relación de Poisson . Estas mediciones son realizadas basadas
en el tiempo de viaje de la ondas P y S (mediciones acústicas) y serán llevadas
a cabo sobre núcleos que se encuentren a condiciones representativas del pozo
en cuanto a esfuerzos, temperatura y presión de poros.

3. Realizar medidas estáticas en núcleos en el laboratorio del módulo de Young y

de la relación de Poisson, mediante el uso de la curva esfuerzo-deformación de
los ensayos triaxiales. Esto también implica someter al núcleo a condiciones de
temperatura, presión de confinamiento y presión del yacimiento. El paso anterior
y éste se deben efectuar en una misma muestra y en un mismo ensayo con
equipo triaxial especializado capaz de realizar ensayos estáticos y dinámicos
simultáneamente.

4. Elaborar correlaciones entre el módulo de Young y la relación de Poisson

estática y dinámica de laboratorio. Esta correlación será usada para extrapolar
el módulo de Young y la relación de Poisson dinámica medida por perfiles y
poder obtener valores continuos de módulos estáticos en el campo. De esta
manera se podrá producir un perfil continuo y calibrado de la variación del
módulo de Young y de la relación de Poisson estática con profundidad.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-19

Coeficiente poroelástico de Biot

Utilizando el principio de esfuerzos efectivos se sabe que parte del esfuerzo total σ
sobre un elemento de material poroso es llevado por el esqueleto mineral (esfuerzo
efectivo, σ’ ) y parte por la presión de poros, tal y como lo demuestra la ecuación
siguiente:

σ = σ '+ αu

Por lo tanto, el coeficiente poroelástico de Biot α mide la eficiencia de la presión

de poros para soportar los esfuerzos totales. En sistemas porosos con buena
interconexión, como el caso de arenas con porosidades y permeabilidades altas,
el valor de α es 1. Por el contrario en sistemas sin porosidad como metales, el
valor de α es 0, y por lo tanto los esfuerzos efectivos y los esfuerzos totales son
iguales.

El valor de α de una roca de interés debe ser medido en el laboratorio utilizando un

tipo de ensayo geomecánico especializado de carga isótropa que produce valores
estáticos. Sin embargo, el valor de α también puede ser estimado por el registro
acústico, pero desafortunadamente éste es un valor dinámico que siempre es
bastante menor que los valores estáticos calculados en el laboratorio. Por lo tanto
es importante tener valores de la constante poroelástica de Biot medida en el
laboratorio y no utilizar los valores calculados directamente por el registro acústico.
La diferencia entre los valores estáticos y dinámicos puede ser substancial, lo que
va a tener un gran efecto en el cómputo de los esfuerzos efectivos.

Sistema QLA2 para el análisis petrofísico geomecánico

El paquete de computación QLA2 fue desarrollado por las empresas

Schlumberger y Geographix para hacer análisis petrofísico. Este programa tiene la
capacidad de trabajar con modelos y/o funciones creadas por el usuario. A este
paquete V.V.A. Consultores, C.A. le ha añadido los módelos geomecánicos para
calcular los módulos elásticos dinámicos, la resistencia de la roca, los esfuerzos de
la formación alrededor del hoyo y el valor de “drawdown” máximo para que no
ocurra el arenamiento. Es importante notar que los resultados de estos
procesamientos ya estarán calibrados para la roca específica del area ya que
incluye los resultados de los ensayos de laboratorio sobre las rocas del área. Una
vez que se obtengan los tiempos de tránsito de las ondas P y S y la densidad total
en función de la profundidad de la compañia de servicio, se hacen, según se

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-20

requiera, los siguientes procesamientos independientes: Propiedades mecánicas

básicas, Estabilidad de hoyos y Resistencia de las arenas.

Procesamiento del Registro Sónico del Pozo LSE-5018

Se presentan los resultados de dos procesamientos de los registros del pozo LSE-
5018, calculados a partir del Software QLA-2. El uso de este software de
procesamiento precisa que los datos estén en formato binario (LBS). Las curvas
datos que van a ser procesados son:

GR
RHO (B)
Tiempos de viaje, Onda P y Onda S

Esta información es suminisrada en formato ASCII, por lo que se hace la

transformación a formato binario y luego se manipula gráficamente para controlar
la imagen (escala, curvas, color, etc). En base a estas curvas datos y del uso de la
opción del QLA2 para crear conjuntos de ecuaciones y modelos se hacen los
cálculos necesarios pa ra la interpretación.

En este trabajo las ecuaciones usadas corresponden al cálculo de la porosidad y el

cálculo de las constantes elásticas dinámicas reseñadas anteriormente
(ecuaciones 10.32, 10.33 y 10.34)

La porosidad se calculó con la siguiente ecuación:

PHID=(RHOM-RHOB)/(RHOM-RHOF)
Donde:

RHOM es la densidad de la matriz,

RHOF es la densidad del fluido,

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-21

Figura 13.9.- Gráfico del registro mecánico del pozo LSE-5018 mostrando velocidades de ondas compresional
y onda de corte y los módulos elásticos.

Un intervalo de los resultados se muestran en dos figuras. El la Figura 13.9 se tiene

5 “tracks” o pistas, en el primero está la curva GR y “VSHALE”, en el segundo se
tiene la calidad del hoyo en base al cáliper y el “BIT SIZE”, el tercer “track”, situado
después de la escala de la profundidad, se presenta la curva dato densidad
(RHOB) y la curva calculada para la porosidad. En el cuarto “track” se presentan los
datos de tiempo de tránsito incluidos: Tiempo de viaje de la onda P, tiempo de
viaje de la onda S y relación Vp/Vs para dejar en el último “track” las constantes
elásticas dinámicas. En la Figura 13.10 se presentan las curvas de un registro
típico de resistividad.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez

 Capítulo 13 Información Geomecánica a Partir de Registros Acústicos Pagina 13-22

Figura 13.10.- Gráfico del registro de resistividades del pozo LSE-5018 mostrando las resistividades llanas y
profundas.

Introducción a la Geomecánica Petrolera Andrés R. Vásquez