Boletín de Geología

Vol. 32, N° 1, enero-junio de 2010

METODOLOGÍA PARA DETERMINAR DENSIDAD

DE MICROFRACTURAS EN ROCAS, A PARTIR
DE PARÁMETROS ANISÓTROPOS. UNA
APROXIMACIÓN EN NÚCLEOS SINTÉTICOS
Jenny Rueda1, Hernan Madero1, Zully Calderón1,
Néstor Saavedra , Germán Ojeda2, Alberto Ortiz2, Carlos Piedrahita2
2

RESUMEN
La anisotropía es una propiedad de las rocas la cual es generada por diversos factores. Comprender cuáles son
estos factores requiere un análisis extenso. Para contribuir en dicho análisis, en esta investigación se estudió
la densidad de microfracturas planteando que éstas son generadoras de anisotropía en las rocas. Para evaluar
el impacto de las microfracturas en la anisotropía, en este trabajo, se desarrolló una metodología experimental
de laboratorio, en la cual se utilizaron modelos sintéticos de rocas con microfracturas simuladas, variando la
densidad de microfracturas (inclusiones); y se construyó un dispositivo para realizar mediciones de ondas de
ultrasonido. La anisotropía se obtuvo a partir de las mediciones ultrasónicas, estableciendo cómo influye la
densidad espacial de las microfracturas en las velocidades de las ondas longitudinales o de compresión P y de
corte o cizalla S. A partir de los resultados obtenidos, se estableció una relación cuantitativa entre la densidad
de microfracturas y la anisotropía de las ondas S cuando se propagan a través de un medio fracturado de
composición sintética, aportando evidencia de laboratorio con respecto a (1) el papel de las microfracturas
como generadoras de anisotropía; (2) la relación directa entre la anisotropía de la onda S desde una densidad de
microfracturas de 0% al 4.96% y para valores mayores del 4.96% la anisotropía presenta una disminución; y (3)
la variación relativa de la anisotropía de las ondas P y ondas S.

Palabras claves: Anisotropía, densidad de microfracturas, parámetros anisótropos, simulación de microfracturas,

ondas ultrasonido.

METHODOLOGY TO DETERMINE DENSITY OF MICROFRACTURE

ON ROCKS FROM ANISOTROPICS PARAMETERS: AN
APPROXIMATION IN SYNTHETICS CORES
ABSTRACT
The anisotropy is a property of the rocks which is generated by diverse factors. To understand which are those
factors, is required an extensive analysis. To contribute in the analysis mentioned before, this investigation
studied the microfractures density planning that the microfractures anisotropy generators in the rocks. To evaluate
the impact of the microfractures in the anisotropy, this work develops an experimental laboratory methodology,
in which synthetic models of rocks were used. This models were elaborated with simulated microfractures
(inclusions), varying their density. The anisotropy was obtained thanks to ultrasound measurements (an
ultrasound device measurement was build up), finding how influences the space density of the microfractures in
the speeds of the longitudinal waves or compression (P) and shears (S). As of the obtained results was establish
a quantitative relationship between the microfractures density and the anisotropy from ultrasonic waves (S)
when they spread through a fractured layer of synthetic composition, contributing with laboratory evidence
due to the fact that (1) the microfractures role as anisotropy generators; (2) the direct relationship between the
anisotropy from wave S and the density of microfractures from 0% to 4.96% and for values bigger than 4.96%
the anisotropy will decrease; and (3) the relative variation of the anisotropy of the waves P and waves S.

KeyWords: Anisotropy, microfracture density, anisotropic parameters, simulation of microfracture, ultrasound

waves

1
Universidad Industrial de Santander, Escuela de Geología, Bucaramanga, Santander, Colombia. e-mail: jennypaolarueda@
yahoo.com; hernandariomadero@yahoo.com
2
Ecopetrol S.A – Instituto Colombiano Del Petróleo, A.A. 4185 Bucaramanga, Santander, Colombia
79
 Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

INTRODUCCIÓN cantidad de esferas dentro del cubo estas van a tener la

misma forma y no varían las propiedades del material.
Una de las necesidades que se presenta en la industria Los modelos isótropos son comúnmente usados para
del petróleo es la de reconocer la orientación y densidad la determinación de las propiedades mecánicas de las
de las microfracturas en los yacimientos, ya que, a partir rocas. Esta suposición de isotropía se asume para un
del conocimiento de este parámetro se aporta evidencia fácil entendimiento y para disminuir la complejidad
para la identificación de trenes de permeabilidad, zonas en los cálculos. Sin embargo, la mayoría de elementos
de acumulación de hidrocarburos, así como volúmenes que componen la corteza terrestre no presentan estas
porosos, además de evaluarse el riesgo de exploración propiedades, sino que por el contrario son “anisótropos”.
y optimizar el tratamiento para la extracción de A diferencia de los materiales isótropos, en materiales
hidrocarburos. anisótropos las propiedades físicas varían de acuerdo
con la orientación en la cual se realizan las medidas.
Varios autores han considerado que las microfracturas Como se observa en la Figura 1b, las esferas tienden
son generadoras de anisotropía en las rocas; basados en a alinearse en cierta dirección, de manera que, si por
esta hipótesis, en esta investigación se buscó determinar ejemplo, se estuviera midiendo su permeabilidad, el
la densidad de microfracturas. El objetivo planteado, flujo tendría una dirección preferencialmente horizontal.
evidenció la necesidad de desarrollar una metodología
propia que permitiera manipular las condiciones de Factores que generan anisotropía en las rocas
los experimentos y obtener información a partir de
parámetros controlados. Esta metodología consistió en Entre los factores que generan anisotropía en las rocas
determinar la anisotropía producida por microfracturas se encuentra la estratificación, las fracturas, el contenido
en núcleos sintéticos en los cuales se manipularon la de arcilla por su estructura laminar, y la foliación en una
forma y densidad de microfracturas. Dicha anisotropía roca metamórfica, donde los minerales se orientan en
se identificó a partir de la medición de ondas ultrasónicas una dirección preferencial. Esta investigación se enfocó
producidas en el laboratorio, con las cuales se estableció en el estudio de la anisotropía que pueden producir las
cómo influye la densidad espacial de las microfracturas microfracturas.
en las velocidades de las ondas longitudinales o de
compresión P y de corte o cizalla S. A partir de lo ¿Porqué trabajar con microfracturas?
anterior fue posible establecer una relación cuantitativa La importancia de determinar la densidad de
entre la densidad de microfracturas y la anisotropía que microfracturas en las rocas está dada por varios
estas producen en un medio sintético similar a las rocas. factores. Cuando se realiza el modelado de yacimientos
naturalmente fracturados se buscan datos que ayuden
ISOTROPÍA VS ANISOTROPÍA a representar las condiciones originales a las que se
encuentran las rocas en el yacimiento, el conocimiento
Los materiales son “isótropos” cuando sus propiedades de la densidad de microfracturas disminuye la
físicas se mantienen constantes a través de ellos sin incertidumbre en los yacimientos, por lo que estos
cambiar su magnitud, independientemente del ángulo datos afectan directamente la representatividad que
u orientación en el cual se estén midiendo. Como se los modelos tienen sobre la realidad, esto implica que
observa en la Figura 1a, los ejes x, y, y z representan las proyecciones económicas sean adecuadas y no se
las propiedades del cubo, de manera que si se tiene una presenten pérdidas

FIGURA 1. Representación gráfica de los materiales isótropos (a) y anisótropos (b). En el primer caso (a), las propiedades
físicas del medio son las mismas independientemente de la dirección. Para el segundo caso (b), las propiedades físicas del medio
cambian de acuerdo con la dirección analizada.

80
 Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita

¿Por qué se requiere una nueva metodología para Ventajas de la metodología convencional
determinar la densidad de microfracturas? • Con el análisis de sección delgada de las rocas
Actualmente los datos de densidad de microfracturas en además de determinar la densidad de microfracturas,
las rocas, se estiman a partir de análisis petrográficos. se observan e identifican el contenido mineralógico,
En este análisis se utiliza una metodología la cual arreglo, tamaño, forma de los componentes,
consiste en, realizar una sección o lámina delgada de la condiciones de formación, clasificación y otras
roca. La elaboración de láminas o secciones delgadas se características que varían según el tipo de roca.
deben aplicar en rocas compactas con un mínimo de tres
cortes, uno para cada plano con el fin de tener un control
tridimensional durante el análisis (INEGI, 2000).
Desventajas de la metodología convencional
Luego de realizada la sección delgada, en un papel • En la elaboración de secciones delgadas se deben
transparente colocado sobre las secciones delgadas se realizar mínimo tres cortes uno para cada plano con el
trazan manualmente las microfracturas que se observan fin de tener un control tridimensional durante el análisis.
con claridad y se mide la longitud de cada una de ellas • Si se realiza solo un corte el dato es areal
(Figura 2) (Samaniego et al., 2003). Con los datos de • El tiempo en el que se obtiene el dato de densidad de
longitudes, se calculan la densidad de microfracturas microfracturas no es inmediato, la elaboración y análisis
(D) Ecuación 1: de la sección requiere un tiempo aproximado de días.
Siendo A el área de la imagen y Li la longitud de cada • Los costos de la elaboración son elevados, en
promedio una sección delgada tiene un precio de
microfractura.
aproximadamente 35 dólares y el análisis 125 dólares.
 1   Li  
D =   * Σ   * 2 (1) De acuerdo a las desventajas que presenta la
 A   2  
herramienta actual en la determinación de la densidad
de microfracturas en las rocas, se planteó emplear una
nueva herramienta, que optimizara estas desventajas.
Para esto se propuso analizar la anisotropía estudiada
desde la propagación de ondas ultrasónicas en un medio
fracturado de composición sintética.

FUNDAMENTOS TEÓRICOS PARA

DETERMINAR LA DENSIDAD DE
MICROFRACTURAS A PARTIR DE LA
ANISOTROPÍA (IDENTIFICADA DESDE
LA VELOCIDAD DE PROPAGACIÓN
DE ONDAS ULTRASÓNICAS)
Ondas en medios elásticos
El comportamiento de la propagación de las velocidades
ultrasónicas en medios porosos es un tema de interés por
su importancia en la industria del petróleo. Una onda
ultrasónica se puede definir como la propagación de una
perturbación en un medio elástico. Cuando una roca
es sometida a cualquier tipo de esfuerzo esta sufre una
deformación temporal. Dado que las rocas se comportan
elásticamente, los estudios de propagación de ondas en
medios rocosos han demostrado que el esfuerzo aplicado
sobre ellas causa una perturbación (deformación), la
FIGURA 2. Ejemplo de trazas de fracturas. Las trazas se cual se propaga a lo largo de toda la roca comportándose
miden manualmente y se calculan los parámetros indicados. como una onda con cierta velocidad. Basados en este
A= área de la imagen; D = densidad de las microfracturas. comportamiento elástico de las ondas en las rocas, se han
Tomado de Samaniego et al., (2003).
realizado investigaciones en el área de la geomecánica,

81
 Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

que han permitido establecer relaciones entre el Esta investigación se enfocó en la simetría hexagonal,
comportamiento elástico de las rocas y las velocidades la cual muestra como casos geofísicos sistemas de
sónicas del medio poroso. Las ondas que se propagan fracturas verticales y paralelas.
en medios elásticos, son principalmente de dos tipos:
ondas P (primarias) y ondas S (secundarias). Las ondas Isotropía hexagonal o transversal
P También llamadas compresivas, son ondas que se
Este tipo de simetría también es llamado anisotropía
propagan por compresión y tensión. El movimiento de
azimutal (Thomsen, 2002) y es el tipo más simple de
las partículas de roca por las ondas P consiste en una
simetría real porque tiene un solo eje de simetría de
alternancia de contracciones y distensiones, es decir,
rotación (aquí 2 de los 3 ejes son equivalentes entre sí, y
se acercan y alejan entre sí moviéndose en la dirección
diferentes del tercero). En el caso más común, el eje de
de propagación de la onda. Las ondas S llamadas de
simetría es vertical (Figura 3), así que puede ser llamado
cizalla ó transversales, se propagan por movimientos
isotropía transversal vertical ó “VTI”. La isotropía a la
perpendiculares a su dirección de propagación.
que se hace referencia está limitada al plano horizontal
(Transversal).
Para caracterizar la anisotropía desde la velocidad de La isotropía transversal vertical (VTI), proporciona una
propagación de estas ondas, es necesario entender el descripción adecuada de la propagación de onda en un
comportamiento del material que se está estudiando, en medio de capas delgadas con interfases horizontales
este caso las rocas, las cuales tienen de acuerdo a su (rocas estratificadas) y formaciones de esquistos
constitución un tipo de simetría específico. estratificados.

Simetría en medios elásticos

De acuerdo al tipo de simetría que puede tener un sólido
(las rocas), existe un número de constantes elásticas
necesarias para describir su comportamiento elástico.
Entre menos complejo sea el tipo de simetría, menos
constantes se necesitan para caracterizar. En total hay
seis tipos de simetría en los medios elásticos.

TABLA 1. Tipos de simetría en medios elásticos.

FIGURA 3. Representación de un medio con isotropía
Simetría en medios elásticos
transversal con eje de simetría vertical, esquema de medios
Ejemplos de casos estratificados. Modificado de Rex Brian et al., (2003).
Constantes
Tipo de geofísicos que
elásticas
Simetría contienen este tipo de
independientes
simetría El otro tipo de isotropía transversal es el que tiene un eje
Muy pocos casos de simetría horizontal, (Figura 4), o HTI, aquí el plano
Cúbica 3 geofísicos, ejemplo, isótropo es el vertical (Thomsen, 2002), un ejemplo
cristal de Halita
de esta es donde las rocas tienen fracturas verticales
Esquistos, secuencias alineadas.
de capas delgadas,
Hexagonal 5 o bien sistema de
fracturas paralelas y
verticales
Esquistos, secuencias
de capas delgadas
horizontales afectadas
Ortorrómbica 9
por un sistema de
fracturas paralelas y
verticales
Material afectado
por dos familias
Monoclínica 13
no ortogonales de
fracturas paralelas. FIGURA 4. Representación de los medios de interés para
Medios con múltiples reservorios fracturados. HTI, Medio con eje de simetría
Triclínica 21 horizontal. Modificado de Rex Brian et al., (2003).
sistemas de fracturas

82
 Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita

Como se observa en la Figura 5, el modelo de simetría Determinación cuantitativa de la anisotropía en

consta de un plano de eje de simetría (el plano x1, x3) y un medios transversalmente isótropos HTI
plano isótropo. En medios HTI las propiedades elásticas
son las mismas en cualquier dirección perpendicular Este tipo de simetría está caracterizada por cinco
al eje (eje de simetría), pero son diferentes paralelas al constantes elásticas independientes (Tabla 1), Ecuación
eje. Además la velocidad es más baja cuando se mide 2, (Thomsen 2002;).
paralela al eje de simetría y más alta cuando se mide
perpendicular al eje de simetría (Rueger, 1996). c11 c13 c13 
c c33 (c33 − 2c44 ) 
 13 
c (c33 − 2c44 ) c33 
CHTI =  13 
 c44 
 c55 
 
 c55  (2)

Estas constantes elásticas (Ecuaciones 3 a 8) se

determinan a partir de la densidad de la roca (r) y las
velocidades de propagación de las ondas medidas en
diferentes direcciones (Figura 6).

C11 = ρV 2 P1 (3)

c66 = ρV 2 s1a (4)

FIGURA 5. Medio transversalmente isótropo con eje de
simetría horizontal HTI, se presenta un plano isótropo y un
plano del eje de simetría, las propiedades medidas paralelas C33 = ρV 2 P 2 (5)
al eje de simetría son bajas comparadas con las medidas
perpendiculares a dicho eje. Tomado de Rueger, (1996).
C 44 = ρV 2 S 1b (6)
El medio de interés, de acuerdo al objetivo planteado es
el medio transversalmente isótropo con eje de simetría
c11 − c12
horizontal o HTI. En este caso la onda S se polariza en c66 = (7)
dos ondas: una onda horizontal (Figura 5, plano X1) la 2
cual es relativamente lenta debido a que viaja en forma
perpendicular a las fracturas, y una onda vertical (Figura 4 ρ 2V 4 P 3 − 2 ρV 2 P 3 (C11 + C 33 + 2C 44 )
5, plano X2), la cual es relativamente rápida, debido a C13 = −C 44 +
+ (C11 + C 44 )(C 33 + C 44 )
que viaja paralelamente a las fracturas. (8)

FIGURA 6. Configuración para la medición de la propagación de ondas para determinar anisotropía en

modelos HTI.

83
 Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

El cálculo de estas constantes permite caracterizar el las microfracturas el espaciamiento y las dimensiones
medio estudiado, y medir el grado de anisotropía de las de las microfracturas deben ser mucho más pequeñas
rocas a partir de los parámetros anisótropos, como se que la longitud de onda (Thomsen, 2002). Debido a
muestra en las ecuaciones 9, 10 y 11 desarrolladas por que las longitudes de onda con las que se trabaja en
Thomsen. laboratorio son menores a 2 centímetros de longitud, en
c11 − c33 esta investigación se estudiaron fracturas menores a 2
ε= cm de longitud y de aperturas menores a 1mm. En este
2c33 (9)
estudio, a las fracturas que cumplen estos requisitos de
c66 − c44 longitud y apertura se les denominan “microfracturas”.
γ = (10)
2c44
(c13 + c44 ) 2 − (c33 − c44 ) 2 (11)
δ =
2c33( c33 − c44 )
e: Anisotropía del material obtenida desde la onda P
(Dependencia de la velocidad de la onda P respecto a la
dirección de propagación).
g: Anisotropía del material obtenida desde la onda S FIGURA 7. Representación en forma gráfica de las
(Dependencia de la velocidad de la onda S respecto a la suposiciones de la teoría de Hudson Tomado de Pachano et
al., (2005).
dirección de propagación).
d: Relación de la anisotropía determinada desde la onda P y S.
Hudson (1981) determinó la densidad de microfracturas
(grietas, inclusiones) como se muestra en la Ecuación (12).
A partir del anterior análisis se identifica el grado de
anisotropía producida por microfracturas en el material nr 3
estudiado, para poder identificar si hay una relación entre = (12)
dicha anisotropía y la densidad de microfracturas en el v
medio, se debe tener conocimiento de este parámetro. ε: Densidad de microfracturas (inclusiones)
n: Número de inclusiones
Modelos de fracturas r: Radio de la inclusión
v: Volumen del cilindro o matriz donde se encuentran
Existen diferentes teorías con las cuales se describe el
dichas inclusiones.
efecto de las fracturas sobre la propagación de ondas.
Estas teorías se encuentran divididas principalmente
A partir de la ecuación 12, se estimará la densidad de
en dos tipos: el primero es conocido como teoría de
microfracturas sobre los núcleos sintéticos, lo cual
grietas (planteada principalmente por Hudson, 1980)
permite manipular los parámetros que incluye dicha
el cual considera la distribución de grietas pequeñas,
ecuación, posteriormente a estos núcleos sintéticos se
desconectadas, paralelas y subparalelas, las cuales
les medirá los parámetros anisótropos para establecer si
tienen forma de monedas; mientras que el segundo
hay alguna relación entre estos factores, la densidad de
tipo se conoce como teoría de fracturas (planteado por
microfracturas y la anisotropía del material.
Schoenberg, 1983) el cual trata las fracturas como un
sistema de planos paralelos infinitos. En este trabajo se
utilizó la teoría propuesta por Hudson (1980), por ser Procedimiento preliminar para desarrollar la
la que mejor se ajusta a la descripción de los modelos metodología propuesta
físicos. Además, por permitir la representación de la
A diferencia de la metodología convencional en la cual
estructura de un yacimiento (Pachano et al., 2005).
la densidad de microfracturas se determina a partir
de análisis petrográficos, donde el tiempo y costos de
Modelo De Grietas En Forma Elipsoidal (Penny obtención son elevados y además de no realizarse un
Shaped) Teoría De Hudson, 1980. análisis en tres planos se obtiene un análisis areal;
Suposiciones: en este trabajo se propone una metodología que
1. Las grietas son elipsoides con una pequeña relación permite determinar la densidad de microfracturas,
de aspecto. 2. No existe interconexión entre las grietas aproximadamente en dos hora y con un costo
(Figura 7). 3. Para leer la anisotropía producida por mínimo, inferior comparado con el de la metodología

84
 Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita

convencional. Además la densidad de microfracturas se P y S, hallar las constantes elásticas que caracterizan el
estima en un volumen de roca. medio HTI (Ecuaciones 3 a la 8).
4. Desde las constantes elásticas, calcular los parámetros
Para entender el efecto de las microfracturas en la anisótropos (Ecuaciones 9 a la 11).
anisotropía del material desde las ondas ultrasónicas, 5. Hallar la relación entre la anisotropía del material
este trabajo desarrolló una metodología experimental obtenida desde los parámetros anisótropos e, g, d y la
de laboratorio. De acuerdo con el objetivo planteado, densidad de microfracturas obtenida desde Hudson
el cual es desarrollar una metodología para determinar 1981 (ecuación 12).
la densidad de microfracturas a partir de parámetros
anisótropos, fue necesario obtener una relación entre Obtención de ondas acústicas. La generación de
dichos parámetros y la densidad de microfracturas. las ondas ultrasónicas se basa, en que una unidad de
Para evaluar el impacto de las microfracturas en pulso envía una señal eléctrica al transductor, el cual
la anisotropía se elaboraron núcleos sintéticos de mediante un cristal piezoeléctrico interno genera una
roca para poder manipular en ellos los parámetros onda de tensión de baja energía y alta frecuencia. A su
de microfracturas. A estos modelos sintéticos se les vez los transductores o palpadores deben ser acoplados
hizo inclusiones simulando microfracturas las cuales a la superficie del material mediante medios acoplantes
tuvieron densidad variable en cada modelo. Además, para transmitir el máximo de energía posible, en este
cada modelo cumplió con las condiciones del modelo caso se empleó inicialmente epóxico de Plata. La onda
HTI (modelo transversalmente isótropo), donde la viaja entonces a través de la roca, siendo captada por
matriz de roca debe ser isótropa y homogénea y la el transductor receptor, que puede situarse en varias
anisotropía solo es entonces causada por microfracturas posiciones en función del método empleado (Figura 6) el
paralelas y verticales. cual, a su vez, convierte la energía de la onda en energía
eléctrica. El tiempo de la transmisión puede visualizarse
Se fabricaron 35 núcleos sintéticos en el Laboratorio de en el osciloscopio, mediante un display de lectura, en
Tecnología Operacional (TOP) del Instituto Colombiano microsegundos, a partir del cual se calcula, conocida
del Petróleo. Estos se elaboraron a partir de una mezcla de la distancia entre emisor y receptor, la velocidad de
arena fina más resina epóxica, cuyas medidas se adecuaron propagación de la onda y, posteriormente, estimar los
al tamaño de un núcleo de perforación (3.8 cm de diámetro parámetros anisótropos.
y 7 cm de longitud), apropiado para realizar experimentos
en laboratorio, con un volumen aproximado de 64 cm3. Estas mediciones se realizaron en el Laboratorio de
Las microfracturas son simuladas con pequeñas láminas Mecánica de Rocas del Instituto Colombiano del Petróleo.
de aluminio de forma circular de 6mm de diámetro y 0.07
mm de espesor, estas se iban ubicando de forma paralela y Montaje de cristales piezoeléctricos en rocas
vertical, una con respecto a otra, a medida que se construían En experimentos realizados previamente sobre rocas
los núcleos. De estos núcleos, 30 se elaboraron variando para medir los tiempos de propagación de las ondas se
el número de inclusiones. El primer núcleo sintético se fijaron directamente los piezoeléctricos a la roca con
fabrica sin inclusiones, al siguiente núcleo se le incrustan epóxico de Plata.
5 inclusiones, y así sucesivamente se van aumentando las
inclusiones de cinco en cinco, hasta llegar a 130 inclusiones, La configuración utilizada para el desarrollo de la
lo que corresponde al 5.38% de densidad de microfracturas metodología propuesta se muestra en la Figura 6, como
según Hudson, 1981; a partir de 130 inclusiones se insertaron se observa en esta Figura, para llevar a cabo las medidas
de 20 en 20 hasta llegar a 250 microfracturas equivalentes a de los parámetros anisótropos, se necesitan, dos pares de
10.55% en la densidad de microfracturas (según Hudson, cristales S y un par P dispuestos sobre las caras planas
1981). Los 5 núcleos restantes se utilizaron para realizar la del cilindro, y tres cristales P ubicados en las paredes
verificación de resultados. del cilindro de manera que un par queda enfrentado
horizontalmente y otro a 45º. Los cristales PZT se
METODOLOGÍA PROPUESTA cortaron antes de realizar el montaje en el laboratorio,
ya que estos se adquirieron en el mercado con
La metodología propuesta consiste en: dimensiones más grandes de las requeridas. La forma
1. Medir los tiempos de propagación de las ondas P y S de los cristales es inicialmente circular con un diámetro
2. Determinar la velocidad de propagación de las ondas de 2cm y un área de 3.14 cm2. Después de realizar el
PyS corte, el área de cada cristal fue de aproximadamente 25
3. A partir de la velocidad de propagación de las ondas

85
 Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

mm2, obteniendo entonces 10 cristales a partir del cristal dispuso de una placa de acero sobre la cual se ubicaron
inicial. La frecuencia de oscilación de los cristales está los cristales adhiriéndolos con epóxico de plata. Esta
directamente relacionada con el área de estos, la cual placa se adaptó de tal manera que su acople a la roca
fue de aproximadamente 80 a 100 KHz para los cristales sea el mejor, ya que de esto depende que la calidad de
S y de 50 a 60 KHz para cristales P después del corte. la onda sea óptima. Debido a que el montaje de los
cristales sobre la roca tiene las disposiciones observadas
De acuerdo con lo anterior, con un cristal S, uno P y en la Figura 6, se tomaron dos pares de placas para
un sobre de epóxico de plata (3g), es posible realizar realizar las medidas longitudinales (Figura 8), es decir
solo dos montajes experimentales. Tanto el costo de los a lo largo del núcleo, ubicando las platinas sobre las
cristales como el epóxico de plata es elevado, con el fin partes planas del cilindro.
de reducir estos costos se trabajó en la implementación
de un dispositivo para el ahorro de cristales, ya que En la Figura 8a se observan las placas para realizar
el promedio de pruebas que se realizaron implicó un medidas longitudinales, en cada platina se observan 3
número elevado de cristales y epóxico de plata lo cual cristales, 1 P y 2 S polarizados a 90º entre ellos. En la
pudo haber limitado la aplicación de la metodología. Figura 8b se observa una estructura Metálica externa la
Este dispositivo se describe a continuación. cual se hizo para fijar las placas, esta es denominada
soporte para núcleos. En la Figura 8c se observan las
Dispositivo para ahorro de cristales placas para realizar medidas diametrales. Estas placas
Como se mencionó anteriormente, se diseñó un tienen una forma tal que en su parte interna son paralelas
dispositivo para optimizar el uso de los cristales. a la curvatura del cilindro (Figura 8c).
Basados en el funcionamiento de los piezoeléctricos se

FIGURA 8. Dispositivo para ahorro de cristales. (a) Placas para medidas longitudinales. (b) Soporte para plugs con el respectivo
montaje para las mediciones longitudinales. (c) Soporte para plugs con el respectivo montaje para las mediciones diametrales.

FIGURA 9. Fotografía de la señal recibida en el osciloscopio de una onda P y una onda S, mediciones realizadas empleando el
dispositivo para el ahorro de cristales.

86
 Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita

Este mecanismo consta de una estructura metálica la ondas, se procedió a realizar las medidas de los tiempos
cual soporta el núcleo y las platinas y se utiliza para de propagación de las ondas en los núcleos sintéticos
presionar las platinas contra el núcleo con ayuda de con microfracturas simuladas; a partir del tiempo de
tornillos para conseguir un mejor acople. Además, propagación se estimaron las velocidades P y S, para
para mejorar aun más el acople se aplicó grasa de baja posteriormente calcular los parámetros anisótropos
viscosidad o de vacío a las placas y a la roca (Figuras (Ecuaciones 9, 10 y 11), desde las constantes elásticas
8b y 8c). La señal que se logró con este dispositivo fue (Ecuaciones 3, 4, 5, 6, 7, y 8). Conociendo el número
óptima, es decir, la amplitud de la señal ocupó más del de inclusiones en cada núcleo sintético, el radio de las
60% de la pantalla del osciloscopio y la señal de ruido inclusiones así como el volumen del núcleo, fue posible
fue mínima, con lo cual se pudo identificar el arribo de establecer la densidad de microfracturas desde Hudson
la onda sin mayor dificultad (Figura 9). 1981, Ecuación 12. Teniendo tanto los parámetros
anisótropos, así como la densidad de microfracturas
Para validar las mediciones de la propagación de ondas (Ver anexo) se plantea examinar dichos datos y analizar
y determinar si la señal obtenida era la correcta, se usó si existe alguna relación entre estos.
como patrón el aluminio, ya que la velocidad de la onda
P y S son conocidas (Tabla 2). En el laboratorio las En las Figuras 10 a 13 se muestran estas relaciones
medidas experimentales obtenidas para el tiempo de gráficamente, así como la desviación estándar de los
propagación de la onda P y S fueron de 6 μs y 13,5 μs datos, en las gráficas la desviación estándar se presenta
respectivamente; los datos reportados se estiman en 5.9 como barras verticales.
μs y 12.5 μs respectivamente. Con lo anterior se muestra
que el montaje es aceptable para medir los parámetros El comportamiento de la onda S (Figura 10) muestra
anisótropos, ya que se observa el retraso normal de la una tendencia lineal clara, se observa una relación
onda S con respecto a la onda P (Tabla 2). directamente proporcional entre la densidad de
microfracturas y la anisotropía caracterizada desde la
TABLA 2. Comparación de los datos de velocidades onda S. Este comportamiento se mantiene hasta un
reportados en el Aluminio comparados con los obtenidos en el porcentaje de densidad de microfracturas del 4.96%,
dispositivo ahorrador de cristales. y en los valores de mayor porcentaje de densidad de
Datos Reportados
Datos de
microfracturas esta relación cambia, disminuyendo
Carbó Rafael et al, considerablemente la anisotropía desde 0,4 hasta 0,15.
Laboratorio
2000

Diámetro (m) 0.038 0.038

Tiempo
5,92E-6 6.00E-6
Longitudinal (s)

Velocidad
Longitudinal 6420.00 6333.33
(m/s)

Tiempo
12.5 E-6 13.5E-6
Transversal (s)

Velocidad
3040 2814
Transversal (m/s)

Después de verificar que la señal obtenida era la adecuada

se procedió a realizar la medición de la propagación de
las ondas P y S en los núcleos sintéticos. Con el fin de FIGURA 10. Relación de la Anisotropía caracterizada desde
obtener repetitividad en los resultados se realizaron la onda S promediada Vs la densidad de microfracturas.
alrededor de cinco mediciones en cada posición de las
microfracturas con respecto al cristal para cada núcleo. Por el contrario, en las gráficas obtenidas de las
relaciones entre los parámetros anisótropos e (Figura
11) y d (Figura 12) con la densidad de microfracturas,
RESULTADOS no se observan tendencias que permitan confirmar el
comportamiento esperado de aumento en la anisotropía
Una vez se obtuvo a partir del dispositivo que contenía a medida que se aumenta el porcentaje de densidad de
los cristales piezoeléctricos la señal correcta de las microfracturas. En la Figura 11 se observa una tendencia

87
 Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

de aumento en los intervalos de 0,0 y 5% de densidad aumentando linealmente la anisotropía a medida que
de microfracturas, pero a diferencia de la Figura 10 aumenta la densidad de microfracturas hasta el 4,96%,
esta muestra fluctuaciones elevadas en la desviación se realizó una regresión lineal obteniendo la ecuación de
estándar; como es el caso en el punto 8% de densidad la recta (Figura 13).
de microfracturas los valores que se obtuvo estaban en
un rango de 0,2 y 0,98 en la Anisotropía obtenida desde
la Onda P, menores fluctuaciones se presentan en los
intervalos de 3 a 5 % y de 8 a 10 % en la densidad de
microfracturas. El mismo comportamiento se presenta
en la Figura 12, con los datos en los intervalos 4-5 % y
9-10 % en densidad de microfracturas.

Las fluctuaciones que se observan con los parámetros

e y d, indican que no hay repetitividad en los datos,
dando incertidumbre a los resultados. Por lo tanto, no es
adecuado apoyarse en las relaciones de estos parámetros
para la determinación de la densidad de microfracturas.
FIGURA 13. Relación de la Anisotropía obtenida desde
la onda S promediada Vs la densidad de microfracturas, se
muestra la relación lineal entre estos parámetros.

Y=0,0742X + 0,0122
Y: Anisotropía caracterizada desde la Onda S (γ)
X: Densidad de microfracturas (%)

Con esta ecuación se determinará la densidad de

microfracturas a partir de parámetros anisótropos en
los núcleos reales, es decir, de acuerdo a la metodología
propuesta, una vez se realicen las mediciones de las
ondas P y S, se hallen las velocidades, las constantes de
rigidez y la anisotropía desde la onda S se determinará
la densidad de microfracturas.
FIGURA 11. Relación de la Anisotropía obtenida desde la
onda P promediada Vs la densidad de microfracturas.
Assad´d et al., (1992), realizaron un experimento
similar llegando hasta un porcentaje de densidad de
microfracturas del 10%, observando también una
tendencia de aumento de anisotropía de la onda S
hasta una densidad de microfracturamiento del 7%,
y por encima de este valor una disminución de la
anisotropía, ellos sugirieron que cuando se supera un
valor crítico de densidad de fracturas (inclusiones) estas
se interconectan entre sí, de manera que las dimensiones
de las inclusiones dejan de ser menores que la longitud
de onda, rompiéndose así las condiciones de anisotropía
en función de la escala.

FIGURA 12. Variación de la relación de ondas P y S (d) Vs la De acuerdo a los resultados obtenidos, además se puede
densidad de microfracturas. sugerir, que el comportamiento observado puede ser
el reflejo de que el medio fracturado alcanza un valor
A partir del análisis realizado a la relación entre la crítico de densidad de microfracturas (inclusiones),
anisotropía caracterizada desde la onda S y la densidad como consecuencia este se hace espacialmente más
de microfracturas (Figura 10), en donde se observa una homogéneo y su comportamiento tiende a ser menos
tendencia directamente proporcional entre la Anisotropía anisótropo.
desde la Onda S y la densidad de microfracturas,

88
 Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita

Aplicabilidad de la curva obtenida de la como consecuencia este se hace espacialmente más

anisotropía del material desde la onda S vs la homogéneo y su comportamiento tiende a ser menos
densidad de microfracturas anisótropo.

Según Crampin (1994), las rocas se pueden clasificar De acuerdo a los resultados obtenidos en las relaciones
según su grado de fracturamiento (Figura 14). La de la anisotropía de los núcleos sintéticos desde las onda
escala de Crampin (1994) varía entre un valor mínimo P (ε) y la relación entre la onda P y S (δ) con la densidad
de fracturamiento de 1.5% (e=0.015) y un máximo de de microfracturas, se descarta el uso de estos parámetros
4.5% (e=0.045). de anisotropía para emplearlos en la determinación de la
densidad de microfracturas, por sus fluctuaciones en la
desviación estándar, ya que no presentaron repetitividad
en los datos.

NOMENCLATURA
A: Área de la imagen para determinar densidad de
FIGURA 14. Esquema de la clasificación del grado de microfracturas, método convencional
fracturamiento de las rocas, según Crampin, (1994). Li: Longitud de cada microfractura
D: Densidad de microfracturas
Crampin (1994) además sugirió que hay un valor de e : Densidad de microfracturas (inclusiones)
densidad de fracturas que corresponde al “fracturamiento n: Número de inclusiones
crítico” de la roca. Como se observa en la Figura 15 el r: Radio de la inclusión
fracturamiento crítico se alcanza en 4.5%. ρ: Densidad de la roca
n: Volumen del cilindro o matriz donde se encuentran
“De acuerdo con la clasificación que planteó Crampin dichas inclusiones.
(1994) para el fracturamiento de las rocas, la curva Vp1: Velocidad P Paralela a las microfracturas
obtenida en esta investigación de la relación entre Vp2: Velocidad P Perpendicular a las microfracturas
la anisotropía desde la onda S (g) y la densidad de Vp3: Velocidad P a 45 grados de las microfracturas
microfracturas (Figura 13) tiene aplicabilidad en Vs1a: Velocidad de propagación S, paralela a las
rocas cuyo grado de fracturamiento sea inferior al microfracturas
fracturamiento crítico, es decir, en un rango de rocas Vs1b: Velocidad de propagación S, perpendicular a las
intactas a moderadamente fracturadas”. microfracturas
ts1a: Tiempo de propagación S, paralela a las
Esta metodología está en proceso de validación con microfracturas
núcleos reales. ts1b: Tiempo de propagación S, perpendicular a las
microfracturas
c11,c12, c13, c33, c44, c66: Constantes elásticas independientes
CONCLUSIONES e : Anisotropía de la onda P
g : Anisotropía de la onda S
La relación obtenida entre la anisotropía caracterizada d : Relación de la anisotropía de la onda P y S
desde la onda S y la densidad de microfracturas demuestra
que la anisotropía está relacionada directamente con la AGRADECIMIENTOS
densidad de microfracturas paralelas, verticales, y que
es posible determinar una relación entre estas dos. El Los autores queremos expresar nuestro agradecimiento
comportamiento observado muestra que a medida que a: ECOPETROL S.A. – ICP y a la UNIVERSIDAD
aumenta la densidad de microfracturas (inclusiones), INDUSTRIAL DE SANTANDER, por patrocinar
aumenta la anisotropía de la onda S, con una tendencia el desarrollo de esta investigación. Al Grupo de
lineal. Este comportamiento se presentó hasta el Investigación Estabilidad de Pozo, A Hernando
porcentaje de densidad de microfracturas del 4.96%, en Altamar, Karen Pachano, Herlin Gonzalez por sus
valores superiores, la anisotropía disminuye, esto puede aportes y asesoría en el entendimiento de la anisotropía.
ser el reflejo de que el medio fracturado alcanza un valor Al personal del laboratorio de Mecánica de Rocas
crítico de densidad de microfracturas (inclusiones), ICP, Jenny Carvajal, Carime Valera, Jorge Galvis,

89
 Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

Tecnología Operacional ICP, Laboratorio de Rocas Samaniego A. F., Álvarez S., Tolzon G. 2003. Estimación
UIS, Laboratorio De Química UIS, por la colaboración de densidades, distribuciones de longitud y longitud
técnica en el desarrollo de los experimentos. total de fracturas; un caso de estudio en la Falla de Los
Planes, La Paz, B.C.S. Boletín de la sociedad geológica
mexicana tomo LVI, ( 1), pp. 1-9
REFERENCIAS
Schoenberg, M. y Gilbert, K.E. 1983. Elastic waves
Ass’ad J. M., Robert H. T., y John A. M. 1992. A Physical in periodically layered media. Technical Programme
Model Study of Microcracks-induced anisotropy, and Abstracts of Papers, 45th Meeting of the European
Geophysics, 57, (12) 1562-1570. Association of Exploration Geophysicists (Oslo): 58-58.

Carbó, R. 2000. Difusión de ondas acústicas por una Thomsen L. 2002. Entendimiento de la Anisotropía Sísmica
esfera hueca con pared de muy poco espesor inmersa en Exploración y Explotación, Boletín AMGE, 41 ( 4).
en agua y llena de aire; Revista de Acústica, 37 (3, 4).

Crampin S. 1994. The Fracture Criticallity of Crustal Trabajo recibido: Marzo 24 de 2010
Trabajo aceptado: Junio 11 de 2010
Rocks, Geophys. J. Int. 118, 428-438.

Crampin S. 1998. Shear-Wave Splitting in a Critical

Crust: The Next Step, University of Edinburgh,
Departatment og Geology and Geophysics, Grant
Institute, West Mains Road, Revue de l’institut Francais
Du Petrole, 53, (5).

Hudson, J. A. 1980. “Overall properties of a cracked

solid”, Math. Proc. Camb. Phil. Soc., 88, 371-384.

Hudson, J.A. 1981. Wave speeds and attenuation of

elastic-waves in material containing cracks, Geophys. J.
R. astr. Soc., 64, 133–150.

INEGI, 2000. Instituto Nacional de Estadística

Geográfica e Informática, La Información Petrográfica y
Paleontológica en la Cartografía Geológica del INEGI,
Guía Normativo-Metodológica, México.

Pachano K. H., Altamar. M., Piedrahita C., Salinas T.,

Calderón Z. 2005. Propagación de Ondas en Medios
Efectivos - Mediciones en Modelos a Escala, Asociación
Colombiana De Ingenieros De Petróleos (ACIPET).

Rex B., Goodway B., Martin C., Gordon U. 2003. The

Impact of Azimuthal Anisotropy on Seismic AVO and
Petrophysical Response in a Fractured Wabamun Gas
Reservoir.

Rueger A. 1996. Reflexion Coefficients and azimutal

AVO analysis in anisotropy media: Ph.D. Thesis,
Colorado School of Mines.

90
 ANEXO. Datos promediados de los núcleos sintéticos con láminas de Aluminio como inclusiones.
Número de Densidad
Long. Diam. Densidad VS1b VS1a C33 C11 C66 C44 C12 C13
Microfract. de tS1b tS1a e d g
(m) (m) (Kg/m3) (m/s) (m/s) (Pa) (Pa) (Pa) (Pa) (Pa) (Pa)
Inclusiones Microfract(%). (ms) (ms)

37,00 36,40 0,018

0 0,00 0,0561 0,0381 1664,42 1515,36 1541,82 10457,4E+6 12245,7E+6 3958,1E+6 3823,2E+6 4329,5E+6 10173,7E+6 0,086 1,095
±0,73 ±0,78 ±0,013

36,00 34,40 0,050

5 0,21 0,0558 0,0381 1713,85 1549,12 1623,83 11671,2E+6 13340,8E+6 4519,5E+6 4114,4E+6 4301,7E+6 8275,4E+6 0,072 0,547
±0,83 ±0,42 ±0,015

36,50 34,50 0,061

10 0,42 0,0563 0,0381 1817,81 1542,51 1634,29 11424,1E+6 14264,9E+6 4855,3E+6 4325,3E+6 4554,3E+6 8320,8E+6 0,124 0,676
±0,23 ±0,19 ±0,011

35,50 33,30 0,068

15 0,63 0,0560 0,0381 1642,57 1577,51 1681,74 9802,5E+6 11250,9E+6 4645,8E+6 4087,7E+6 1959,3E+6 5440,7E+6 0,075 0,521
±0,28 ±0,28 ±0,001

37,70 35,10 0,077

20 0,86 0,0552 0,0381 1780,46 1464,22 1573,05 12704,8E+6 14518,5E+6 4406,9E+6 3817,3E+6 5704,7E+6 10473,8E+6 0,072 0,556
±0,20 ±0,69 ±0,017

91
36,90 34,10 0,085
25 1,06 0,0559 0,0381 1683,16 1514,91 1639,30 12465,8E+6 13148,9E+6 4523,2E+6 3862,8E+6 4102,6E+6 9608,0E+6 0,027 0,501
±0,10 ±0,10 ±0,001

38,00 34,70 0,098

30 1,27 0,0561 0,0381 1818,35 1478,65 1616,71 10472,7E+6 14379,4E+6 4752,7E+6 3976,7E+6 4874,0E+6 9743,0E+6 0,189 1,102
±0,70 ±0,10 ±0,022

40,40 36,40 0,116

35 1,49 0,0556 0,0381 1670,62 1374,98 1525,81 8949,8E+6 11065,3E+6 3889,4E+6 3158,7E+6 3286,5E+6 7914,7E+6 0,119 0,866
±0,46 ±0,12 ±0,016

39,00 34,40 0,141

40 1,69 0,0561 0,0381 1692,97 1440,38 1629,87 10651,6E+6 12807,6E+6 4498,2E+6 3513,5E+6 3811,2E+6 8366,4E+6 0,102 0,598
±0,83 ±0,61 ±0,045

39,10 34,30 0,150

45 1,91 0,0557 0,0381 1702,90 1424,30 1623,62 9656,0E+6 11319,7E+6 4489,1E+6 3454,5E+6 2341,6E+6 7064,5E+6 0,086 0,603
±0,10 ±0,10 ±0,001
Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita

Densidad
Número de Long. Diam. Densidad VS1b VS1a C33 C11 C66 C44 C12 C13
de tS1b tS1a e d g
Microfract. (m) (m) (Kg/m3) (m/s) (m/s) (Pa) (Pa) (Pa) (Pa) (Pa) (Pa)
Microfract. (ms) (ms)
 37,10 32,00 0,172
60 2,54 0,0560 0,0381 1848,55 1509,16 1749,70 10754,2E+6 14487,5E+6 5659,3E+6 4210,2E+6 3168,8E+6 76,12E+8 0,174 0,691
±0,10 ±0,14 ±0,006

38,40 32,50 0,196

65 2,75 0,0559 0,0381 1661,87 1457,48 1718,84 8547,3E+6 12126,8E+6 4910,0E+6 3530,8E+6 2306,8E+6 6679,4E+6 0,209 0,922
±0,50 ±0,17 ±0,021

38,50 32,40 0,206

70 2,96 0,0560 0,0381 1735,77 1455,16 1729,05 10353,6E+6 13153,6E+6 5190,3E+6 3676,4E+6 2773,0E+6 7822,4E+6 0,135 0,634
±0,60 ±0,53 ±0,026

40,20 32,90 0,245

75 3,17 0,0561 0,0381 1672,42 1396,29 1704,01 10087,1E+6 11372,6E+6 4856,3E+6 3260,6E+6 1660,0E+6 8279,6E+6 0,073 0,659
±0,19 ±0,21 ±0,011

40,60 33,00 0,255

80 3,37 0,0562 0,0381 1679,47 1385,80 1702,77 8559,9E+6 11365,4E+6 4869,7E+6 3225,6E+6 1626,0E+6 6672,4E+6 0,166 0,767
±0,41 ±0,20 ±0,011

40,60 3,26 0,276

85 3,59 0,0560 0,0381 1734,99 1381,30 1720,70 8796,7E+6 11674,3E+6 5137,6E+6 3310,6E+6 1399,0E+6 6271,2E+6 0,165 0,640
±0,46 ±0,46 ±0,004

38,90 30,70 0,305

92
95 4,00 0,0563 0,0381 1717,78 1445,70 1834,75 9760,6E+6 14764,8E+6 5786,2E+6 3590,6E+6 3192,4E+6 8935,6E+6 0,304 1,196
±0,43 ±0,89 ±0,035

40,70 31,90 0,314

100 4,12 0,0575 0,0381 1669,68 1411,79 1801,25 7158,9E+6 11815,6E+6 5417,3E+6 3327,9E+6 0981,0E+6 8052,6E+6 0,325 2,094
±0,10 ±0,10 ±0,001

40,20 30,60 0,363

105 4,43 0,0561 0,0381 1660,14 1397,05 1835,68 11621,7E+6 13066,7E+6 5594,2E+6 3240,5E+6 1878,2E+6 9489,9E+6 0,062 0,471
±0,21 ±0,12 ±0,013

40,80 30,60 0,389

115 4,86 0,0560 0,0381 1803,58 1373,68 1832,08 11952,6E+6 14723,1E+6 6053,8E+6 3403,4E+6 2615,5E+6 10226,3E+6 0,116 0,551
±0,21 ±0,12 ±0,002

46,20 34,40 0,402

120 4,96 0,0573 0,0381 1686,83 1238,29 1662,68 10598,3E+6 13821,7E+6 4663,4E+6 2586,5E+6 4494,9E+6 10007,3E+6 0,152 0,556
±0,21 ±0,23 ±0,016

Densidad
Número de Long. Diam. Densidad VS1b VS1a C33 C11 C66 C44 C12 C13
de tS1b tS1a e d g
Microfract. (m) (m) (Kg/m3) (m/s) (m/s) (Pa) (Pa) (Pa) (Pa) (Pa) (Pa)
Microfract. (ms) (ms)

39,70 34,60 0,160

Metodología para determinar densidad de microfracturas en rocas, a partir de parámetros anisótropos. Una aproximación en núcleos sintéticos

125 5,03 0,0589 0,0381 1668,01 1483,40 1703,68 8199,8E+6 13370,7E+6 4841,5E+6 3670,5E+6 3687,8E+6 5610,8E+6 0,351 1,034
±0,21 ±0,12 ±0,004
 37,00 35,60 0,040
130 5,48 0,0561 0,0381 1832,22 1515,67 1575,25 7214,8E+6 14431,4E+6 4546,6E+6 4209,1E+6 5338,1E+6 5356,4E+6 0,500 1,901
±0,21 ±0,23 ±0,006

38,30 36,30 0,057

150 6,25 0,0569 0,0381 1660,01 1484,86 1566,67 10253,7E+6 13422,0E+6 4074,4E+6 3660,0E+6 5273,2E+6 7729,9E+6 0,155 0,639
±0,21 ±0,12 ±0,001

33,60 31,50 0,070

170 7,81 0,0516 0,0381 1854,91 1533,65 1638,00 5563,2E+6 15931,7E+6 4978,4E+6 4363,0E+6 5974,8E+6 6307,1E+6 0,932 8,426
±0,21 ±0,69 ±0,023

37,50 35,80 0,050

190 7,87 0,0572 0,0381 1687,49 1524,38 1598,18 7167,7E+6 14162,6E+6 4310,3E+6 3921,6E+6 5542,0E+6 2568,1E+6 0,594 0,738
±0,21 ±0,23 ±0,007

35,80 35,40 0,012

210 8,92 0,0558 0,0381 1702,98 1556,66 1574,44 10431,0E+6 14991,7E+6 4222,0E+6 4126,7E+6 6547,7E+6 10301,7E+6 0,245 1,911
±0,21 ±0,50 ±0,012

35,50 35,10 0,011

230 9,92 0,0549 0,0381 1711,75 1546,48 1564,10 10484,6E+6 14030,8E+6 4187,6E+6 4093,8E+6 5655,5E+6 7800,2E+6 0,195 0,821
±0,21 ±0,10 ±0,001

35,10 33,10 0,062

250 10,55 0,0561 0,0381 1677,32 1598,29 1694,86 10538,5E+6 15592,2E+6 4818,2E+6 4284,8E+6 5955,7E+6 6926,1E+6 0,240 0,657
±0,10 ±0,10 ±0,001

93
Jenny Rueda, Hernan Madero, Zully Calderón, Néstor Saavedra, Germán Ojeda, Alberto Ortiz y Carlos Piedrahita