El valor umbral optimo del RQD

sobre redes fractales de fracturas

Carlos Paredes*, Santiago de Vicente*, Jess Pedro de Vicente, Alberto Herrera, Maria
Lpez, Ismael Osuna, Juan Sordo, Ricardo Lan**.
* Departamento de Matemtica Aplicada y Mtodos Informticos
** Departamento de Explotacin de Recursos Minerales y Obras Subterrneas

ste trabajo constituye la aportacin realizada en la Adems, se sabe que, a partir de un nico conjunto
E asignatura Complejidad, Caos y Fractales por
cinco de sus alumnos matriculados. Su contenido
de espaciados, es posible obtener un rango de valores
posibles del RQD cambiando el valor de t y/o de L. Si bien,
aborda la determinacin apropiada del RQD (Rock Quality frente al aumento de L el RQD posee un comportamiento
Designation, Deere et al. 1967) como parmetro asinttico hacia el valor real. En los estudios
discriminante entre macizos segn su grado de geomecnicos se ha establecido t=0.1m como valor umbral
fracturacin a tenor del comportamiento anistropo de universal, si bien se trata de un valor elegido
ste; un problema frecuente que surge en la realizacin de arbitrariamente. Por ello, dado que el cambio en el valor
estudios geotcnicos de los medios fracturados. El tomado para t tiene efectos notables sobre el RQD, se
resultado obtenido ofrece una solucin ptima, hace necesario estudiar un procedimiento que permita
descartando la alternativa tradicional (de acuerdo con la establecer el valor ms apropiado para el umbral a utilizar
definicin de RQD se toma un espaciado mnimo de 0.1m basndose en las caractersticas de los espaciados entre
entre discontinuidades), para cada macizo fracturado y las discontinuidades que se encuentran en el macizo en
segn sea la caracterizacin fractal de los espaciados vez de usar un valor tpico.
entre las fracturas. Las discontinuidades existentes en la masa rocosa
le confieren a sta caractersticas de heterogeneidad y
Introduccin anisotropa, ambas de vital importancia ya que afectan, por
ejemplo, al comportamiento de un almacenamiento de
La utilizacin del RQD en los trabajos de ingeniera CO2, de residuos slidos urbanos o radiactivos, un acufero
geotcnica es una prctica muy comn, ya que se trata de subterrneo, un yacimiento de hidrocarburos, etc. que
un parmetro que permite caracterizar y categorizar pudiera encontrarse en su interior. Esta variabilidad
geomecnicamente los macizos rocosos. Al menos dos espacio-direccional en la geometra de la red de
razones lo fundamentan: a) resulta ser una medida de discontinuidades es la que se utiliza en este trabajo para
obtencin relativamente simple, de la misma forma que se seleccionar el valor apropiado de t para calcular el RQD.
hace sobre sondeos puede hacerse sobre afloramientos Aunque la heterogeneidad es probablemente uno de los
rocosos; y b) la definicin del RQD se encuentra aspectos ms importantes en la caracterizacin de los
universalmente reconocida, determinada con el mismo macizos rocosos, se ha tomado la anisotropa como
procedimiento, lo que permite efectuar comparaciones primera caracterstica a considerar en la definicin de t, al
entre diferentes ubicaciones rocosas y su clasificacin. ser sta una cualidad intrnseca a las discontinuidades.
Para el clculo del RQD, ya sea sobre sondeos o Se presenta una tcnica para la determinacin del
itinerarios lineales de longitud L entre dos umbral RQD apropiado para macizos rocosos en los que
discontinuidades, se evala el cociente entre el total de las discontinuidades se ubican espacialmente de forma
sumar las longitudes de tramos no fracturados mayores invariante al cambio de escala, lo que da lugar a que los
que un umbral (t) entre L. Es decir, el RQD esta definido en espaciados se distribuyan mediante una ley de potencia.
trminos del porcentaje de tramos de roca intacta Es frecuente encontrar este tipo de distribuciones habida
(espaciados entre discontinuidades) mayores que t. De cuenta de que los espaciados son el resultado de un
esta forma se facilita una cuantificacin porcentual del proceso de fragmentacin repetitivo y autoorganizado de
grado de roca intacta -no fracturada- respecto al valor los bloques litosfricos sobre un amplio rango de escalas.
mnimo considerado t. Una evidencia experimental de este fenmeno es la ley de
Sin embargo, con el fin de progresar y mejorar la Gutemberg-Richter. Adems se muestra cmo el uso del
aplicacin y comprensin del comportamiento del RQD es valor ptimo calculado mejora la utilidad del RQD ya que lo
preciso considerar algunos de los inconvenientes convierte en una medida mucho ms discriminatoria que
asociados a su uso. En la prctica, resulta que los valores considera no slo la heterogeneidad del macizo, en cuanto
del RQD entre 40% - 60% se encuentran mucho menos a la densidad de la fracturacin, sino tambin la
frecuentemente que los valores extremos superiores al anisotropa.
80% o inferiores al 10%. Si se supone que no hay una
razn fsica que justifique que puede encontrarse una Variacin del RQD en los macizos
geometra especfica de las discontinuidades a favor de rocosos
otra (y por lo tanto de los espaciados entre estas),
entonces la menor frecuencia de algunos valores del RQD Los macizos rocosos se encuentran y se han
puede deberse a cmo ha sido definido ste. encontrado, desde los tiempos de su formacin, sometidos
a una serie de esfuerzos que, en su mayora, son de

1
origen tectnico. Superada la capacidad competente de la asociados a las direcciones medias (Tabla 1) que han sido
roca, su efecto produce la ruptura, dando lugar a fracturas establecidas para cada familia a partir de datos de campo
y fallas. Comnmente, cada discontinuidad tiene un estructurales.
aspecto de superficie tridimensional alabeada quasi plana Dado que el RQD depende de la longitud total de
y suele agruparse en familias pseudoparalelas con espaciados mayores que un valor umbral, para una
aquellas que son singenticas para un campo de esfuerzos longitud de muestreo y umbral fijos, a medida que aumenta
particular. Cuando se analiza el conjunto de familias de el nmero de fracturas en el transecto, su frecuencia de
fracturas sobre un transecto (camino unidimensional que fracturacin aumenta, la longitud de tramos no fracturados
se traza sobre el macizo: itinerario, sondeo, calicata, disminuye y de la misma forma lo hace el RQD. Si los
galera, etc.) la importancia de la fracturacin se cuantifica espaciados entre fracturas (tramos de roca intacta) se
-1
mediante la frecuencia de fracturacin ( en m ), valor que distribuyen segn una funcin de densidad de probabilidad
se encuentra relacionado con el nmero de f(x), la probabilidad de encontrar un espaciado
discontinuidades de cada familia que se encuentran en el aleatoriamente elegido se encuentre entre x y x+dx es
(i, para la familia i-sima). En 1983, los profesores J. A. f(x)dx. Si sobre un transecto de longitud L la frecuencia de
Priest y J. A. Hudson mostraron que, para una masa discontinuidades es , se encontrarn un total de L
rocosa con n conjuntos de discontinuidades planas y discontinuidades o espaciados. Entonces, el nmero de
persistentemente pseudo paralelas, la frecuencia de espaciados entre x y x+dx es Lf(x), cubriendo una
fracturacin global a lo largo de un transecto es: longitud de Lf(x)dx. Lo que permite definir un valor terico
n del RQD mediante el % de longitud, frente al total L, de la
= i cos i (1) suma de todas las longitudes Lf(x)dx. Esta suma se
i =1 efecta para todos los espaciados x comprendidos entre el
donde i es el ngulo agudo que forma la direccin del valor umbral admisible t y la longitud mxima posible L:
L
RQD(, t , L) = 100 xf ( x )dx
transecto con la normal a la familia i-sima.
(2)
t
2
N Normalmente t>>L (t=0.1m y L10 m), por lo que en
la prctica el valor terico del RQD se determina a partir de
RQD(,t,). Esta clara dependencia del RQD con
provoca su comportamiento anistropo. Sin embargo, la
obtencin de los valores del RQD a partir de datos reales
max=7.017m-1 presenta leves variaciones frente a la direccin de
52.2
muestreo. Sus valores son normalmente prximos al 80%,
140
.04
por lo que su representacin polar sera prcticamente
istropa. Este hecho, lejos de la realidad anistropa del
medio fracturado, sugiere que el RQD debiera corregirse
para reflejar la anisotropa del medio analizado y tomar un
rango de valores ms amplio (entre 0% y 100%) y
discriminatorio (no slo el 20% de los posibles) de lo que
min=4.21m-1 ocurre para cada direccin. Como solucin a esta situacin
se propone examinar en (2) si una eleccin apropiada de t
pudiera maximizar el rango de valores que toma el RQD
para cada que se encuentre entre los valores extremos
de la frecuencia de fracturacin max y min. La bsqueda
del valor de t apropiado (que se denotar por t*) se
Figura 1. Representacin en coordenadas polares de la efectuar analticamente mediante la maximizacin de:
frecuencia de (1) para las 4 familias de fracturas de la
Tabla 1, para direcciones de muestreo en el plano RQD( min , max , t , ) =
(3)
horizontal. Los segmentos en gris denotan las direcciones
de las frecuencias mnimas locales.
RQD( min , t , ) RQD( max , t , )
es decir, determinando el t* que verifica:
Familia Direccin () Frecuencia (m-1) d
dir.buz./buz. RQD( min , max , t*, ) = 0 (4)
1 58 / 25 1.96 dt
2 186 / 81 2.87 Sustituyendo (2) en (3), derivando para obtener la
3 50 / 72 3.61 ecuacin (4), y teniendo en cuenta la regla de Leibnitz para
4 120 / 68 2.44 las integrales definidas se llega a la ecuacin general:

Tabla 1. Parmetros de las familias utilizadas para min t * f min ( t*) + xf min ( x )dx
t* t
representar la Figura 1.
(5)

Esta expresin muestra cmo, para una situacin max t * f max ( t*) + xf max ( x )dx = 0
t* t
general en la que cada familia de discontinuidades posee
un valor de su frecuencia de fracturacin diferente del que es vlida para aquellas funciones de distribucin de
resto, la frecuencia de fracturacin global del macizo espaciados definidas en t* y siendo fmin(x) y fmax(x) aquellas
(Figura 1) tiene un comportamiento anistropo (segn la asociadas a las direcciones en las que se encuentran min
direccin de muestreo).La anisotropa de en funcin de la y max, respectivamente.
direccin del transecto posee una serie de mnimos locales

2
Pareto como distribucin de espaciados log(b 2 n )
en el RQD =
log b
Aunque para la modelizacin estadstica de los Una variable aleatoria se distribuye segn una ley
espaciados entre discontinuidades rocosas se han de Pareto si la probabilidad es inversamente proporcional a
aplicado muchas funciones de distribucin distintas, entre una potencia del valor de la variable:
las ms frecuentemente utilizadas (exponencial negativa, Pr[X > x ] = Cx ; x> (6)
siendo el parmetro de la distribucin. A partir de (6) se
gamma, Weibull, log normal), una de las ms simples y
cuya aplicacin tiene un fundamento terico inspirado en el
determina la funcin de densidad de probabilidad:
proceso fsico de fracturacin, es la distribucin de Pareto.
El material rocoso puede fragmentarse por muy diversas f ( x ) = Cx1 ; x> (7)
causas: por meteorizacin o ataque directo de los y, teniendo en cuenta la condicin de normalizacin en f(x)
elementos atmosfricos, por explosin o variaciones
sbitas de la presin por impacto de una onda de choque,
o por exceso de los lmites de resistencia del material a la
(como
f ( x )dx = 1 ), la constante de proporcionalidad

accin de fuerzas tectnicas. En cualquiera de estos C = > 0 . Una vez establecida f(x), puede utilizarse
casos, el macizo rompe con mayor frecuencia a travs de
ciertos planos de debilidad, diaclasas, juntas, fracturas o sta en (2) para determinar (5). Sin embargo, interesa
fallas que a travs de la roca sana. Una gran variedad de expresar f(x) incorporando antes de operar, para reducir
estos procesos de fragmentacin puede interpretarse en el nmero de parmetros en RQD(,t,). Para ello, si se
trminos de dimensin fractal de probabilidad o grado de tiene en cuenta que el espaciado medio <x> es el inverso
fracturacin . de la frecuencia de fracturacin , entonces:
Segn este comportanmiento, el nmero de
fragmentos N(x) cuyo tamao es mayor que x es < x >= x x1 dx = = 1 (8)
proporcional a x- (habitualmente 2<<3). Un sencillo
1
modelo sinttico, que utiliza la informacin de estos planos donde se deduce para que <x>>0m ha de ser >1. Esta
persistentes en el medio, y la invarianza al cambio de relacin permite reescribir (7) como:

escala del proceso de fragmentacin, permite generar una
1
distribucin de tamaos de bloques (Figura 2) que se rige f (x) = ; x> (9)
por una ley de potencia x-. El proceso, a cada escala, x x
mantiene sin romper un nmero de bloques n. El resto, y sustituyendo (9) en (2) e integrando cuando L:
rompe segn un patrn de fracturas ortogonales. Sobre
estos segundos se repite el proceso a una escala b veces 1 1
inferior, manteniendo n bloques sin romper y fracturando el 100 ; t
resto.
RQD(, t , ) = t ; >1 (10)
Si, por ejemplo, como se muestra tambin en la 100; t <
Figura 2, los cuatro primeros bloques se reubican
secuencialmente, y se muestrease segn un transecto La incorporacin de una distribucin de espaciados
lineal, el conjunto de espaciados entre fracturas tambin se de tipo fractal crea un comportamiento de tipo ley de
distribuir de la misma forma que el tamao de bloques. potencia en el RQD.
Para el modelo bidimensional, puede demostrarse que, en
general:

j=0 j=1 j=2 j=3 j=4

Figura 2. (arriba) Ilustracin del modelo de fragmentacin 2D hasta j=4 que genera una distribucin de bloques de tipo
fractal (=1.58) con n=1 y b=2. (abajo) Modelo simple de medio fracturado creado mediante redistribucin aleatoria de los
bloques generados para j=5 y colocacin lineal (hipottica traza de muestreo en rojo).

3
 Si en sta se sustituye directamente la expresin
(9), tras la correspondiente manipulacin aritmtica, se
tiene una ecuacin en t*:
RQD(t=0.1,min,max) max
max max max 1

t *max max max
(12)
min
RQD(t=0.2,min,max) min min min 1
=0
t * min min min
t=0.1 cuya solucin, teniendo en cuenta (8), puede obtenerse en
t=0.2 funcin de los extremos de la frecuencia de fracturacin:
t=0.5
( min 1)
t=1.5
t* = max (13)
min max
min ( max 1)
siendo:
Figura 3. Dependencia del RQD terico (10) con la
( max 1)( min 1)
frecuencia de fracturacin para diferentes valores del =
umbral de muestreo t>=0.08m. El valor de [0, ] ( max min )
-1

considerando que >1 y (8).

Dado que estos extremos se encuentran
relacionados con el valor de los exponentes de la
Dado que las frecuencias para un medio fracturado:
-1 distribucin de Pareto, la solucin puede ponerse en
0<min<max< , las diferentes curvas trazadas en la Figura
funcin de los mismos:
3, tomando diferentes valores de t en (10), permiten
1
confirmar que el rango abarcado por los valores del RQD
1 max min
entre los extremos de las frecuencias vara con t. Para el
t* = max (14)
min 1
ejemplo de esta figura, el rango es mayor para t=0.2m que
para t=0.1m; por lo que, para la estimacin del RQD, el
segundo valor de t sera una mejor eleccin que el primero. De esta forma, se obtiene un valor ptimo del
De esta forma, para cada pareja de valores (max, min), el umbral que se encuentra explcitamente definido en
comportamiento del rango para cada valor de t, (Figura 4) trminos de los valores extremos de la frecuencia de
posee un valor de t=t* para el que el rango es mximo. fracturacin (13) o de los exponentes de Pareto en las
Ciertamente ste no es el 100% porque tambin depende mismas direcciones observadas en el macizo rocoso. Para
de los valores max y min que se hayan obtenido con (1), su determinacin, ser necesario llevar a cabo los
con lo que el mximo es para una configuracin particular oportunos estudios de campo geolgico-estructurales que
de discontinuidades. permitan obtener con fiabilidad el comportamiento
anistropo de la frecuencia de fracturacin global segn
(1), del que se extraern max y min. Sobre la direccin que
corresponda a cada uno se estimarn max y min,
respectivamente, as como (supuesto istropo) sobre un
conjunto de espaciados tomados en ambas direcciones. Si,
para una de ellas, la muestra de espaciados medida es
{xi}i=1,m (tras ser ordenados de mayor a menor:
xixj,1i<jm), puede utilizarse el estimador de mxima
verosimilitud de Hill:
1
r
= r 1 (log x i log x r +1 )
i=1
C = 1 = (r / m) x r+1
Figura 4. Dependencia del rango terico (3), utilizando el para la distribucin de Pareto (6), si tras su representacin
RQD (10), con los valores del umbral de muestreo grfica bilogartmica se encuentra un ajuste lineal a partir
t>=0.08m, para max=7.017m-1 y min=4.21m-1, de la Fig. 1. de x=xr+1. Para la aplicacin prctica del resultado obtenido
al ejemplo de la Tabla 1, los valores extremos de las
La determinacin analtica valor umbral t=t* para el frecuencias han de ser extrados de la Figura 1:
que se localiza el valor mximo del rango (3), max=7.017m-1 y min=4.21m-1, a los que corresponden
suponindose (9) como distribucin de espaciados, se ha max=2.279 y min=1.507, segn (8), respectivamente. El
realizado tomando la frmula de Leibnitz (5) en la que: valor ptimo de t*=0.264m a partir de (13) (14), el cual
proporciona los valores extremos del RQD de 54.3% y
xf ( x ) = 0 21.6%, es decir, un rango til del 32%. Comparando este
t valor con el rango del 14%, obtenido al utilizar el valor
con lo que la ecuacin algebraica no lineal a resolver es: tradicional de t=0.1m, como se aprecia en la Figura 4, se
max t * f max ( t*) min t * f min ( t*) = 0 (11) entiende cmo, en este ejemplo, el uso de t* ha
incrementado el rango al menos duplicndolo. En la Figura
5 se ha representado la variacin direccional del RQD,

4
para el conjunto de familias de la Tabla 1, calculado con el todas las direcciones posibles utilizando el umbral ptimo
valor de t*, el valor de t tradicional y un valor umbral de amplia el rango de variacin a ms del 70%, valor muy
t=1.0m, para la misma orientacin horizontal de la Figura 1. superior al 27.5% que se obtiene utilizando el valor umbral
Puede observarse cmo el uso de t* incrementa no solo el tradicional.
rango sino el grado de su anisotropa. Esto demuestra
cmo el uso de t* hace que la funcin RQD(,t*,) sea
mucho ms discriminatoria que con cualquier otro valor del
umbral.
N
0.1
RQD con t=0.1m

min 0.2

RQD con t=1m 0.4

0.6

5.0

max
RQD con tptimo =0.264m
=0.378m

Figura 5. Diagramas polares del RQD para diferentes

umbrales.

Dado que, segn la expresin (14), el valor de t*

para espaciados distribuidos segn Pareto depende 0.2
fundamentalmente de los valores extremos de la
frecuencia global, puede ser interesante examinar cmo t* min
vara con stos. Segn la Figura 6 puede observarse cmo
para valores de las frecuencias que son usuales en la 0.4
prctica de campo (normalmente inferiores a 10m-1) los 0.6
0.8
valores de t* son muy superiores a 0.1m. La aplicacin de 1.0
2.0
ste valor popular requerira trabajar con rocas 5.0
extremadamente fracturadas, con valores extremos altos
(superiores a 10m-1).

Aplicacin a la variacin del RQD en 3D max

Figura 6. Valores del umbral ptimo t* (m) de RQD, segn
Los anlisis presentados anteriormente se las frecuencias de fracturacin (izda.) y segn los
encuentran todos referidos a una seccin plana horizontal exponentes de la funcin de distribucin de Pareto
de muestreo. Sin embargo, a travs de la testificacin (=0.08m) en las mismas direcciones que las frecuencias
sobre sondeos, el muestreo de las estructuras del macizo extremas.
rocoso fracturado puede hacerse en cualquier direccin del
espacio, desde un punto. As, los principios de clculo Es muy frecuente considerar que los espaciados se
aplicados en 2D pueden extenderse a 3D mediante el uso comporten siguiendo las funciones de distribucin
de los valores apropiados de frecuencias extremas exponencial negativa, Gamma o Weibull. Algunos trabajos
(mxima y mnima) de fracturacin extrados de la publicados presentan un valor de t* para la primera
frecuencia global en 3D. distribucin que depende solamente de max y min. Para
La identificacin de estos valores extremos ha de Gamma y Weibull, el resultado en t* es el mismo que para
hacerse sobre la representacin de la densidad sobre una la exponencial negativa:
proyeccin semiesfrica equiangular como la de la Figura
2
7. Para el ejemplo tratado en este trabajo se han obtenido
los valores de max=7.82m (N234E, 25.2) y min=1.72m
-1 -1 t* = ln max
max min min
(N108E, 59.4) a los que corresponden max=2.279 y
min=1.507, segn (8), respectivamente. El umbral ptimo si se supone que el parmetro de escala que aparece en
es t*=0.378m, si =0.08m, calculado mediante las stas es la unidad.
expresiones (13) o (14).
El resultado del calculo en 3D para el RQD La comparacin de estos valores (representados en
proporciona unos valores mximo del 78% y mnimo del la Figura 8 drcha) con los que se han obtenido con la
7.5% lo cual permite afirmar que el anlisis extendido a ecuacin (14) muestra cmo los valores de t* son (Figura 8
izda.) mucho mayores para la distribucin de Pareto que

5
para la exponencial cuando las frecuencias de fracturacin las discontinuidades que atraviesan la masa rocosa sigan
-1
son muy bajas (<4m ) y contrariamente cuando estas una distribucin de Pareto, dependiente de: (a) las
frecuencias son elevadas (>8m-1). Para el ejemplo utilizado frecuencias mxima y mnima de fracturacin (max y min);
en este trabajo, t*=0.49m si los espaciados se comportan (b) los parmetros de la distribucin de Pareto: espaciado
segn una exponencial negativa. Este valor es ligeramente mnimo () y exponentes en las direcciones de los
superior al obtenido por Pareto. Esto implicara que, en la extremos estas las frecuencias (max y min).
determinacin prctica del RQD (2), si se considera la Mediante su anlisis grfico anistropo se han
distribucin de Pareto, entren a formar parte tramos ms establecido las direcciones en las que se dan los valores
cortos que si se aplica la distribucin exponencial. En extremos del RQD, coincidiendo el mximo RQD con la
consecuencia se modifican los extremos del RQD: 79.6% y direccin de menor frecuencia de fracturacin y el mnimo
9.6%, y el rango de valores abarcado 70%. As pues, para RQD con la de mxima frecuencia. La tcnica aqu
una correcta aplicacin prctica del procedimiento es presentada para espaciados anistropos que se
requisito indispensable un calibrado correcto de la funcin distribuyen segn Pareto es tambin apropiada para
de distribucin de los espaciados. macizos rocosos no homogneos, simplemente estudiando
cmo se distribuye espacialmente la frecuencia global,
N localizando sus extremos, y utilizando estos valores para el
6.
clculo de t*. En este procedimiento se han obtenido dos
0
expresiones en las que t* es proporcional a , una
dependiente de (max y min) y la otra de (max y min),
5 .0

5.0
suponiendo que los espaciados poseen un istropo.

4.0
3.0

6.0
2.0
7.0

0.2

(m-1)
min
0.0 4.0 8.0 0.4

N 0.6
40 0.8
.0 1.0

5.0
30.

.0
40
0

30
.0 max
40.0

70.0
.0
20
10.0

A
RQD (%)
min
0.0 50 100

Figura 7. Grficos de la proyeccin semiesfrica

equiangular de la frecuencia global 3D de discontinuidades
(izda.) y del RQD obtenido con un umbral ptimo de 0.37m B
(drcha.), para las familias de discontinuidades de la Tabla
1, (trazadas segn sus crculos mayores) y de las que se
supone tienen un =0.08m istropo. max
Figura 8. Valores del umbral ptimo t* (m) de RQD, segn
Conclusiones las frecuencias de fracturacin, para la funcin de
distribucin exponencial negativa (drcha.); Zonificacin de
Se ha demostrado que la utilizacin de un valor t* frente a las funciones de distribucin analizadas, A:
ptimo t* frente al tradicionalmente recomendado como t*expon>t*Pareto, B: t*expon<t*Pareto.
longitud umbral en la determinacin del RQD permite hacer
ms rigurosa y discriminante la clasificacin geomecnica La correcta caracterizacin de resulta ser crtica en
de un macizo fracturado. El valor t* es, para el caso en que el tratamiento anteriormente expuesto y determinacin de

6
t*. El valor de es la cota inferior del umbral calculado. La
estimacin de depende de la eficacia, calidad y detalle de
los estudios efectuados en campo y de la poblacin
muestreada de espaciados. Por ello se recomienda no solo
asignar un valor estimado de , sino tambin conocer la
incertidumbre asociada a ste y, por ende, al RQD.
Adems, considerando el criterio establecido para el
clculo del umbral ptimo, el rango de valores
comprendido entre el mximo y mnimo RQD es el mayor
que se puede obtener, lo cual amplia la diversidad de
valores en la clasificacin y, por ende, su naturaleza
discriminatoria.

Referencias
Aban, I. B., Meerschaert, M. M., Panorska, A. K. 2006.
Parameter estimation of the truncated Pareto distribution.
J. Am. Stat. Ass., 101,473, 270-277(8).
Choi, S. Y., Park, H. D. 2004. Variation of rock quality
designation (RQD) with scanline orientation and length : a
case study in Korea. Int. J. Rock Mech. Min. Sci. 41, 207-
221.
Bodau, F. K., Long, L. T. 1994. The fractal character of
fracture spacing and RQD. Int. J. Rock Mech. Min. Sci. &
Geomech. Abstr. 31, 2, 127-134.
Deere, D. U. Hendron, A. J., Patton, F. D., Cording, E. J.
1967. Design of surface and near surface rock construction
th
in rock. In Failure and breakage of rock, proc. 8 U. S.
symp. rock mech., (ed. C. Fairhurst), 237-302, New York:
Soc. Min. Engrs, Am. Inst. Min. Metall. Pet. Engrs.
Harrison, J. P. 1999. Selection of the threshold value in
RQD assessments. Int. J. Rock Mech. Min. Sci. Technical
Note, 36, 673-685.
Hill, B. 1975. A simple general approach to inference about
the tail of a distribution. Ann. Statist. 3, 5, 1163-1173.
Newman, M. E. J. 2005. Power laws, Pareto distributions
and Zips law. Contemporary Physics, 46, 323-351.
Priest, S. D. 1993. Discontinuity analysis for rock
engineering. Chapman Hall. London. 473 p.