1

METODOLOGIA PARA LA CARACTERIZACION DE YACIMIENTOS NO

CONVENCIONALES DE GAS SHALE EN COLOMBIA-VIM








RICARDO ANDRES PEREZ RIVERA
JUAN DIEGO PERDOMO RINCON







UNIVERSIDAD SURCOLOMBIANA
FACULTAD DE INGENIERA
PROGRAMA DE PETRLEOS
NEIVA
2013

2


METODOLOGIA PARA LA CARACTERIZACION DE YACIMIENTOS NO
CONVENCIONALES DE GAS SHALE EN COLOMBIA-VIM





RICARDO ANDRES PEREZ RIVERA
JUAN DIEGO PERDOMO RINCON



Trabajo de grado presentado como requisito acadmico para optar al titulo de
Ingeniero de Petrleos.

Director.
MSc. JAIRO ANTONIO SEPLVEDA GAONA
Ingeniero de Petrleos



UNIVERSIDAD SURCOLOMBIANA
FACULTAD DE INGENIERA
PROGRAMA DE PETRLEOS
NEIVA
2013




3

Nota de aceptacin


___________________________________________
___________________________________________
___________________________________________
___________________________________________
___________________________________________
___________________________________________







___________________________________________
Firma del Director




___________________________________________
Firma del Evaluador




___________________________________________
Firma del Evaluador







Neiva, Huila Enero 2013

4

DEDICATORIAS.


Con todas mis fuerzas dedico ste trabajo a todas las personas que nos
colaboraron, para bien o para mal, a terminar nuestro proyecto. Que los
quiero mucho, que los bendigo de todo corazn y que gracias.
RICARDO ANDRES PEREZ RIVERA


JUAN DIEGO PERDOMO RINCON.



















5

TABLA DE CONTENIDO

Pag.

INTRODUCCIN 8
1. GAS SHALE EN EL MUNDO Y EN COLOMBIA 9
2. PROPIEDADES DE LOS YACIMIENTOS DE GAS SHALE 15
2.1. Fractabilidad y Mineraloga 15
2.2. Espesor y Red Poral 20
2.3. Fisilidad y Encogimiento de la Matriz 21
2.4. Presin capilar 23
2.5. Contenido de gas 24
2.5.1. Isoterma de adsorcin 27
2.5.1.1. Volumen de Langmuir 28
2.5.1.2. Presin de Langmuir 28
2.6. Madurez trmica 29
2.6.1. Mtodo de reflectancia de la vitrinita 29
2.7. Carbn Orgnico Total 30
2.7.1. Rock-Eval Pirlisis 31
2.8.
Presin del yacimiento y transiente de presin en el
yacimiento
33
2.9. Descripcin de Corazones 36
2.9.1. Mtodos Convencionales para la medicin de porosidad 37
2.9.1.1. Picnmetro de Baja Presin (LPP 37
3.
EVALUACION DE PROPIEDADES DE YACIMIENTOS DE
GAS SHALE
41
3.1. F.A.S.T. CBM 42
3.1.1. Clculos Volumtricos 42
3.1.2.
Cambio de la Matriz y la Relacin Porosidad/Permeabilidad.
Porosidad de fractura y su relacin con el encogimiento de la
matriz.
44
3.1.2.1. Seiddle y Huitt 44
3.1.2.2. Palmer y Mansoori 47
3.2. Presin de flujo variable 50
3.3. Teora del anlisis del transiente curvas tipo 50
3.3.1.
Clculo de los parmetros: Curvas tipo radiales Pozos de
Gas
53
3.3.2.
Clculo de los parmetros: Curvas tipo para Fractura de
conductividad finita
53
3.4. Presin Capilares 55
3.4.1. Mtodo Analtico 55
3.4.2. Brooks Corey 55
3.5. Contenido de gas 56
3.5.1. Isoterma de Langmuir 56
3.5.1.1. Adsorcin 56
3.5.1.2. Formulacin de la isoterma de langmuir 56
6

3.5.2. Los clculos de yacimientos de gas 57
3.5.2.1. Ecuacin de gas libre 57
3.5.2.1.1. Ecuacin de gas libre Yacimientos de Shale 57
3.5.2.1.2. Ecuacin de gas adsorbido Yacimientos de Shale 57
3.6. Madurez trmica 58
3.7. Carbn Orgnico Total 58
3.8. Descripcin de Corazones 59
4.
EVALUACION ESPECIFICA DE LAS PROPIEDADES DEL
POZO USCO-69
60
4.1.1. Prediccin y evaluacin de un pozo colombiano USCO 69 68
4.2 Arreglo o disposicin de pozos 73
5. CONCLUSIONES 75
REFERENCIAS 77














INTRODUCCION

Diversas cuencas sedimentarias de Colombia, contienen secuencias litolgicas de
excelentes condiciones generadoras de hidrocarburo. La ventana de generacin
de los shales y los carbones incluidos en dichas secuencias, muestra que la fase
de generacin de gas biognico y termognico se ha venido presentando a lo
7

largo de varios perodos prolongados de tiempo, en donde las condiciones de
entrampamiento de dichos gases han podido generar grandes acumulaciones de
metano.
La aplicacin de un procedimiento detallado en donde se incluyen los principales
aspectos a tener en cuenta a la hora de verificar el potencial de un rea, para la
ocurrencia de estos gases ha mostrado que las posibilidades de identificar
grandes acumulaciones de metano en shale de varias cuencas colombianas son
bastante altas. Lo anterior, sumado a la estratgica posicin de nuestro pas y al
inters de empresas extranjeras con experiencia en el desarrollo de este tipo de
proyectos, nos convierte en un potencial productor fuerte de metano.













CAPITULO I GAS SHALE EN EL MUNDO Y EN COLOMBIA

El uso de gas natural se ha venido implementando desde hace miles de aos,
fueron los chinos quines transportaron el gas natural a travs de bambes para
evaporaciones en la salmuera para la obtencin de sal.

El primer pozo de gas natural en los EEUU ocurri en 1815 (accidentalmente),
cuando el capitn James Wilson perforando un domo salino encontr una
estructura de shale que contena grandes cantidades de gas. Treinta (30) aos
8

despus William Tompkins utiliz el gas natural en procesos de manufacturacin
de la sal llevndolo a travs de tuberas de madera para su utilizacin en un horno
de gas (Burning Springs). A partir del ao 1980, la industria petrolera ha buscado
muchas maneras de hacer ste proceso de extraccin ms eficiente, debido a las
grandes cantidades de gas que se han encontrado en estos yacimientos de
shale.
1

Desde los aos 70s y 80s, el temor por la disminucin de las reservas de gas
natural demand una investigacin para desarrollar mtodos y estimar el volumen
de gas en yacimientos no convencionales, como los de gas shale, areniscas
apretadas y mantos de carbn as como investigaciones en mtodos de
extraccin. Debido a que el proceso de extraccin del petrleo y el gas en estos
tipos de yacimientos fue muy desafiante, su recuperacin o recobro durante estas
dos (2) dcadas era ampliamente anti-econmica
2
. Gracias a los mtodos de
extraccin mejorados, los recursos recuperables de gas shale en USA estn ahora
siendo aprovechados para suministrar la demanda de energa domstica.


Figura 1.1 Tendencia del consumo de las reservas de hidrocarburos
Convencionales y no convencionales
En el ao 2000, la produccin de gas shale representaba el 1% del suministro del
gas natural estadounidense. Sin embargo desde el 2005 la cifra ha podido
alcanzar la cuarta parte del abastecimiento de la poblacin americana, que
contina en constante aumento debido en gran parte al desarrollo intensivo de la
extraccin del gas shale que hizo viable su produccin y que hoy en da cuenta
con 862 trillones de pies cbicos que son econmicamente factibles de extraer. Lo
que ha llevado al gobierno estadounidense a aumentar la proyeccin de consumo
nacional de gas shale hasta en un 45% y una disminucin de la dependencia por
gas extranjero del 11% al 1%.

9

Desde una perspectiva ambiental, el gas natural es ms amigable que los
productos derivados del petrleo debido a que es ms limpio: emite menos gases
de efecto invernadero y contaminantes durante la combustin (mercurio, dixido
de azufre entre otros provenientes del carbn). Sin embargo hay que tener
cuidado con respecto al impacto ambiental y la gestin de los procesos que
podran afectar a los acuferos (una mala cementacin y tratamientos de aguas
residuales deficientes) y la calidad del agua subterrnea.

Un estudio solicitado por la US Energy Information Administration
4
permiti revelar
que existen 6,622 trillones de pies cbicos tcnicamente recuperables, una cifra
prometedora dado los 6,609 trillones de pies cbicos de reservas probadas de gas
natural en el mundo. Estas cifras, sumadas a las proyecciones de consumo e
intercambio de gas en los diversos mercados mundiales, adems de las
condiciones de ciertos pases, han permitido establecer al gas shale como
alternativa real de inversin y donde existen casos emblemticos como el de
China y EE.UU., que en el muy largo plazo podran convertirse en exportadores de
gas natural.

Ahora en el caso de Suramrica, las reservas de gas shale para Chile se
estimaron en 64 trillones de pies cbicos que supera con creces a los 3.5 trillones
de pies cbicos de reservas probadas de gas natural, mientras que Argentina
posee reservas estimadas por 774 trillones de pies cbicos y probadas de 13.4
trillones de pies cbicos. Esto convierte a Argentina en un potencial exportador a
largo plazo, que la ubica como el tercer pas del mundo con recursos potenciales
de gas no convencional, detrs de China y muy cerca de Estados Unidos. Chile,
que obtiene el 52% de su consumo va importaciones tiene una buena oportunidad
de iniciar la explotacin y produccin de sus reservas de gas shale, o por el
contrario reiniciar importaciones desde Argentina. Ver figura 1.2

La consultora Advanced Resources International
3
tiene como fin reportar los
recursos mundiales del gas shale; esta compaa fue contratada por el gobierno
de EE.UU para realizar un estudio de las reservas de gas shale en los diferentes
pases a nivel mundial. Segn la investigacin realizada, los recursos de gas no
convencional de China suman 1275 TPC; los de Estados Unidos 862 TPC , y los
de Argentina 774 TPC. El mundo petrolero concibe esta posibilidad como la gran
oportunidad de prolongar la vida til de muchsimos yacimientos que declinan la
produccin con el paso del tiempo. Los petroleros estn contentos. Sucede que
las reservas de hidrocarburos evaluadas podran ser mucho ms importantes con
las nuevas tecnologas que permiten extraer combustibles fsiles de otros
reservorios del subsuelo.

10


Figura 1.2 Cuencas de gas shale en el Cono Sur
3

Desafortunadamente Colombia no cuenta con un buen horizonte en el
abastecimiento energtico a futuro por parte de los combustibles fsiles. Todos los
campos de mayor produccin en el pas se encuentran en la fase de recuperacin
asistida (recobro mejorado), campos maduros, y se necesita buscar una fuente
alterna para el abastecimiento energtico del pas. Como se ha venido diciendo es
el gas shale la fuente que posibilitara que el pas siga siendo autosuficiente, o por
el contrario pasaramos a ser importadores de GNL de gas natural venezolano.
Las reservas totales (probadas, probables y posibles) de gas natural en Colombia
son de apenas 8.45 TPC
6
y la mayora se encuentran hacia el norte del pas. La
produccin de gas natural se haba mantenido en alza hasta el 2010, desde
entonces se ha mantenido constante y se teme que para el 2015 si no hay un
cambio en el plan de abastecimiento nacional el pas ser un importador neto de
gas natural. Por el contrario el consumo tendr una tasa de crecimiento promedio
interanual (2011-2020) del 4.8%, un aumento de cerca de 400 MPCD en el
escenario ms alarmante. Por lo visto se requiere la colaboracin interinstitucional
tanto del sector pblico (MME, CREG, UPME, ECOPETROL) como del sector
privado (Pacific Rubiales) para implementar soluciones hacia la problemtica del
abastecimiento energtico del pas.
11


Figura 1.3 Reservas estimadas de gas no-convencional en el mundo (2010)

Desde el 2008, la Agencia Nacional de Hidrocarburos ha venido evaluando el
potencial de los yacimientos no convencionales en Colombia con un estimado de
28-32 Tpc
7
. Aqu las cuencas que presentan mayor prosperidad segn el informe
de Arthur D. Little
7
en yacimientos no convencionales de gas shale son el Valle
Medio del Magdalena, Csar y Ranchera; sin embargo se han encontrado indicios
de ste tipo de hidrocarburo en la cuenca del Valle Inferior del Magdalena, en el
bloque Guama en los pozos Pedernalito-1x y Cotorra-1x.
8
12


Figura 1.4 Cuencas potenciales de gas shale en Colombia










13

CAPITULO II PROPIEDADES DE LOS YACIMIENTOS DE GAS SHALE

A menudo la produccin atractiva inicial de un yacimiento de gas shale declina
rpidamente. Esto hace que sea econmicamente marginado y algunas veces
antiptico operacionalmente. As que para producir eficientemente ste tipo de
yacimientos, es importante entender los mecanismos de flujo y los parmetros
propiedades que controlan a los fluidos y a la roca, y que tienen un impacto
significativo en el desarrollo a un tiempo determinado de produccin de estos
recursos.

Es el primer paso para entender cmo funcionan estos tipos de yacimientos no
convencionales, porque va a permitir la viabilidad del proyecto al agregar un buen
nmero de datos necesarios para argumentar buenas estimaciones con los
procesos de simulacin y el incremento de la produccin en un futuro cercano. Por
eso es que se evalan las propiedades ms crticas en un escenario previo para
realizar proyecciones acerca de las condiciones ptimas y especificaciones
tcnicas de los equipos.
9


2.1 Fractabilidad y Mineraloga

Las compaas energticas estn agregando los recursos de hidrocarburos no
convencionales a sus portafolios. Usualmente se refieren a ellos como yacimientos
con una permeabilidad extremadamente baja y no pueden ser producidos a ratas
o volmenes econmicos sin haber realizado una estimulacin en las regiones
cerca al wellbore. La estimulacin recomendada para ste tipo de recursos es la
de aplicar el procedimiento de fracturamiento hidrulico que se utiliza para
incrementar el rea de drenaje en los pozos, un rea de contacto mayor entre
estos y la formacin, con lo que se crean caminos para que los hidrocarburos
puedan ser producidos.

Las consideraciones de los modelos petrofsicos son muy importantes para
determinar el proceso de fracturamiento que se va a llevar a cabo. Pues es ste el
que determina el tipo de fluido que va a acompaar al propante, la seleccin del
tipo de propante (tamao y volumen) y dnde tiene que ser iniciado el proceso de
fracturamiento. Adems de stas consideraciones geomecnicas tambin se
deben tener en cuenta las propiedades geoqumicas con las que cuente el pozo: la
mineraloga y la sensibilidad del fluido al fracturamiento, la utilizacin de cidos en
el proceso de estimulacin y si los surfactantes son benficos para el proceso,
entre otros.

Es de sta manera como se ve que la fractabilidad es una propiedad fundamental
en el diseo del fracturamiento, debido a que est directamente relacionada con la
composicin de la roca, se espera que el contenido de arcilla en ella no supere el
45%
10
(el contenido de arcilla es proporcional a la capacidad de absorcin de
impacto que posee sta) y que haya un buen contenido de cuarzo en la roca de
14

materiales no dctiles que permitan el buen desarrollo de sa estimulacin. La
fractabilidad puede ser evaluada de dos (2) maneras: Una es la ecuacin de Wang
y Gale
11
para calcular el Indice de fractabilidad y la otra es la utilizacin de la razn
de Poisson y el mdulo de Young
12
.

El Brittleness Index (BI) ndice de fractabilidad
11
es una herramienta sencilla que
permite evaluar la facilidad que posee la roca para fracturarse. Est en funcin de
la composicin mineralgica del shale y hay una ecuacin propuesta reportada en
la literatura que relaciona el contenido mineralgico y el TOC.

(2.1)


Donde:

BI : Indice de Fractabilidad
Q : Cuarzo
Cl : Arcilla
Dol : Dolomita
Lm : Caliza (Calcita)
TOC : Contenido Orgnico Total

En s, consiste en el clculo de un porcentaje o razn que determina la cantidad
fracturable y refractable (cuarzo y dolomita) de la roca lo que significa que sta
resistir altas temperaturas sin ablandarse y tendr una buena estimulacin en la
fractura haciendo el proceso ms adecuado.

Otra manera de evaluar los esfuerzos necesarios que determinan la fractura, es
interpretando algunas caractersticas de la roca como la de reflejar la habilidad de
sta a fallar bajo presin y mantener una fractura una vez la roca ha sido
fracturada, las cuales son evaluadas por la razn de Poisson y el mdulo de
Young respectivamente. La ductilidad juega un papel importante porque no
permite que el yacimiento sea propenso a fracturarse debido a que sta propiedad
intentar sanar la fractura. Por el contrario al ser el yacimiento ms fractable ser
ms propenso a ser naturalmente fracturado y responder ms favorablemente a
los tratamientos para el fracturamiento hidrulico.

En ste mtodo se cuantifica el factor de fractabilidad teniendo en cuenta lo
anteriormente mencionado, osea combinando ambas propiedades mecnicas de
la roca. Tambin difiere de otros mtodos porque no est basado enteramente en
mediciones de core
12
. Utiliza interpretaciones petrofsicas que cubren el shale as
como las fronteras de las capas de roca por lo que el intervalo que va a ser
fracturado necesita ser enteramente registrado.

Como se puede observar en la figura 2.1 entre menor sea el valor de la Razn de
15

Poisson y mayor sea el Mdulo de Young la roca es ms fractable, si stas
relaciones son inversas la roca tendra un comportamiento dctil. El concepto
general de la grfica es la divisin parcial de los cuadrantes entre puntos de shale
dctil y fractable que seran los cuadrantes noreste y suroeste respectivamente.
Pero qu son exactamente el Mdulo de Young y la Razn de Poisson?


Como se puede ver en la figura 2.2 al
aplicarse una fuerza F a una
cuerda/cilindro de masa m, de
compresin tensin se va a poder
apreciar una deformacin un cambio
de longitud. Esta deformacin va a ser
directamente proporcional a la fuerza
que se aplique y a la longitud de la
cuerda e inversamente proporcional al
nmero de cuerdas y al rea transversal
con las que se cuenten. Por lo que la
ecuacin del Mdulo de Young puede
expresarse de la siguiente manera:


Fig. 2.1 Cuerda/Cilindro

Por lo que se convertira en

Siendo Y la constante de proporcionalidad mdulo de Young.
La ecuacin que se encuentra en la literatura es la siguiente y puede ser
interpretada como una analoga a la ecuacin anteriormente mostrada.

(2.2)

Dnde:

E : Mdulo de Young (psia)
c : Esfuerzo de agotamiento - strain (adimensional)
o : Presin - stress (psia)
En la figura 2.2 se ilustra el comportamiento de la vs

16


Figura 2.2 Diagrama tensin - deformacin.

El mdulo de Young viene representado por la tangente a la curva en cada punto.
Para materiales como el acero (y valores de tensin dentro del rango de
reversibilidad completa de deformaciones) resulta aproximadamente constante
dentro del lmite elstico. Sin embargo si el material est representado por una
curva de tensin-deformacin que no tiene tramos lineales, no se puede usar la
expresin anteriormente mostrada (Eq 2.11). Por lo que se puede decir que la
razn entre el estiramiento y la deformacin no es directamente proporcional y se
define un Mdulo de Young aparente al definirlo como el incremento de esfuerzo
aplicado al material y el cambio correspondiente a la deformacin unitaria que
experimenta en la direccin de aplicacin del esfuerzo:

Como se puede ver el mdulo de Young se refiere a la razn entre la presin axial
y la deformacin. Describe la naturaleza elstica de una sustancia dada y puede
describir convenientemente la cantidad de deformacin de un objeto dado cuando
una presin es aplicada. A menudo se refiere como a la fuerza o rigidez del
material. Entre mayor sea el valor del mdulo de Young menos deformacin
ocurrir a una presin dada.

La razn de Poisson () relaciona el cambio en tamao de un objeto en diferentes
ejes. Cuando la fuerza compresiva es aplicada a un eje en particular de un
material, habr una deformacin por tensin a lo largo de los diferentes ejes. Esta
propiedad es la razn de los esfuerzos de contraccin a los esfuerzos de
extensin. Para darle un valor a la direccin, es positiva cuando el esfuerzo de
17

agotamiento ocurre en la direccin de una fuerza de extensin.



Figura 2.4 Deformacin de un cuerpo y Razn de Poisson


Se expresara, si la fuerza de compresin/tensin estuviera siendo aplicada sobre
el eje axial, de sta manera:

(2.3)


Siendo la tensin o el esfuerzo (strain).


Figura 2.3 Ejemplo de grafico de PR vs YM

Richman et al (2008)
12
us el Mdulo de Young y la Razn de Poisson
normalizados para estimar el ndice de fractabilidad de la siguiente manera:

(2.4)
18

(2.5)

(2.6)


Dnde:
YM_n : Mdulo de Young normalizado
YMS_C: Mdulo de Young
PR_n: Razn de Poisson normalizada
PR_C: Razn de Poisson


2.2 Espesor y Red Poral

El espesor de un yacimiento es muy importante debido a que la produccin de gas
se ve directamente influenciada por el volumen del espacio poroso con el que se
cuenta en conjunto con el gas adsorbido. Anteriormente slo se tena en cuenta
una (1) sola fase de gas contenida en el yacimiento, la del espacio poroso. Esto
haca que el volumen del yacimiento se viera sobreestimado, debido a que se
ignoraba la capa adsorbida que tomaba una parte del volumen poroso
adhirindose al rea poral. Lo que significaba que el volumen poroso se vea
reducido y se necesitaba reformular el balance de materia y en s todo el sistema
de transporte de los hidrocarburos dentro del yacimiento
13
.

Entonces entre mayor sea el espesor neto del yacimiento, mayores volmenes de
gas se van a encontrar. Como sabemos la porosidad es relativamente muy baja,
por lo tanto entre mayor espesor de formacin se tenga, mayor ser la produccin
de gas libre que se va a tener en el yacimiento. Este gas libre es el que se
encuentra en los poros y en las fracturas en la mayora de materia inorgnica.
Esta es la razn por la cual los yacimientos de gas se depletan tan rpidamente,
porque su deplecin est en funcin del volumen de gas libre contenido en la roca
(sin haber realizado un proceso de estimulacin.)

El contenido de gas y la capacidad de adsorcin son los parmetros claves para la
determinacin de gas original in-situ. El gas original in-situ es calculado
mediante
14
:

(2.7)


Donde:

: Gas original en sitio (in situ)
: rea, acres
19

: Espesor de la formacin, ft

: Densidad del shale en

: Contenido de gas promedio de la matriz

2.3 Fisilidad y Encogimiento de la Matriz

La fisilidad es la propiedad que poseen algunas rocas de partirse o dividirse por
planos (extenuados en lminas muy delgadas o sheets) separados entre s por
distancias pequeas. Esto es el resultado de la disposicin subparalela de los
minerales laminares que caractersticamente componen la roca como las arcillas,
micas, cloritas, etc. Si la roca presenta fisilidad se denomina lutita, sino se llama
fangolita.
19
La permeabilidad es la capacidad de un material para permitir que un
fluido fluya sin alterar su estructura interna. Se afirma que un material es
permeable si deja pasar a travs de l una cantidad apreciable de fluido en un
tiempo dado, e impermeable si la cantidad de fluido es nula. En nuestro caso la
permeabilidad es una propiedad muy importante debido a que controla la entrega
de energa del pozo. Esta propiedad est directamente relacionada con la fisilidad
y con la red poral del yacimiento porque en conjunto con la presin es la que
modifica el camino por el que va a fluir el gas.

La permeabilidad en los shales o en los mantos de carbn est definida por las
fracturas y la matriz de la roca. La matriz almacena el gas por adsorcin y en la
matriz se refleja el flujo de gas por el proceso de difusin. Pero qu tiene que ver
la fisilidad con la permeabilidad? En concreto se trata de un proceso que se ve
bien reflejado cuando se analiza la curva de permeabilidad absoluta vs presin y
que se conoce como encogimiento de la matriz (matrix shrinkage). Este proceso,
utilizando en los simuladores slo la permeabilidad de las fracturas, consiste en
que la curva de Ki vs K observada en el simulador presenta un comportamiento
normal al aplicar una presin hasta el punto de quiebre: La matriz se encoge hasta
un punto mximo de presin, hasta el punto de presin de fractura y luego la
relacin K/Ki aumenta ms de 1. Esto se puede explicar debido a que la roca se
fractura y que el espacio brindado por las lminas subparalelas es mayor, al
juntarse el espacio poroso de las fracturas logradas antes de la intervencin de la
presin y despus se ve un incremento considerable en la permeabilidad de la
roca.

El modelo modificado de Palmer y Mansoori mostrado en la figura 2.5 corresponde
a un conjunto de parmetros que se han tenido en cuenta para estudiar el proceso
de desorcin y su relacin con la fisilidad, la permeabilidad y consecuentemente
con la produccin de estos tipos de pozos como son:

- El factor de compactacin pues ste orienta las fracturas horizontal o
verticalmente dependiendo del valor que se tenga (rango de 0 a 1 con los
lmites correspondientes a una orientacin vertical u horizontal
respectivamente).
20

- La presin de desorcin, puesto que el encogimiento slo ocurre una vez el
gas ha sido desorbido (que slo ocurre por debajo de la presin de
desorcin).
- El mdulo dependiente de la porosidad, mdulo bruto Bulk, as como el
mdulo de restriccin Axial M.

Figura 2.5 Razn de Permeabilidades K/Ki vs Presin

La ecuacin de Mansoori y Palmer modificada es la siguiente:

(2.8)

Donde:

= Permeabilidad final
= Permeabilidad inicial
= porosidad final

= porosidad inicial
n = exponente (normalmente seteado en 3, aunque podra llegar a ser mayor a 12)

/i puede ser hallada al integrar la siguiente ecuacin en un rango de cambio de
presin:

Por encima de la presin de desorcin

21


Por debajo de la presin de desorcin

Dnde
a = constante del Mdulo bruto Bulk (K)
b = constante del Mdulo de Restriccin-axial (M)
dP = cambio de presin
d = cambio de porosidad
f = factor de acoplamiento (de 0 a 1), (adimensional)
p = exponente del Mdulo de Restriccin-axial (M)
P = presin del yacimiento (psia)
P
L
= constante de presin de Langmuir (psia)
q = exponente del Mdulo bruto Bulk (K)

= tension maxima que puede ocurrir cuando P se acerca a 0 (adimensional)
= compresibilidad del grano, slo de la matriz y no de las fracturas, (1/psi)


2.4 Presin capilar

Cuando el agua inyectada, durante un procedimiento de recuperacin mejorada,
se pone en contacto con la matriz a travs de las superficies de fractura, el
mecanismo de imbibicin de la fase agua empieza a desempear un papel
importante en el rendimiento de la produccin en el pozo. El efecto imbibicin es
causado por la presin capilar entre las fases de gas y de agua. Cuanto menor
sea la permeabilidad del yacimiento, mayor ser la presin capilar. En los
yacimientos de baja permeabilidad, la presin capilar puede ser de cientos psi
ms (Holditch, 1979), y por lo tanto el efecto de imbibicin que resulta en la
retencin de lquidos se considera como un problema particularmente grave
(Economides y Martin, 2007). El proceso de imbibicin est dominado por las
fuerzas de la presin capilar. En un estudio realizado, dos curvas de presin
capilar se generan en base a correlaciones empricas (Gdanski et. al. 2009) y en la
que se utilizo la matriz del shale, como se muestra en la figura 2.6.Para un caso
base la

se utiliza tanto para la matriz y la zona de invasin, mientras

se
utiliza para la zona de invasin, para el caso en el que la invasin del agua altere
la presin capilar. Para las fracturas con gran permeabilidad, la presin capilar es
muy pequea, esencialmente considerado cero.






22




Figura 2.6 Curvas de presin capilar utilizados en los modelos de simulacin
20


2.5 Contenido de gas

En los yacimientos de gas en mantos de carbn para determinar el contenido de
gas se indican los volmenes de gas expulsado (a partir de muestras de carbn
colocados en contenedores), gas residual (de carbn triturado), y gas perdido
(calculado).

En yacimientos de gas shale es habitual reportar el gas in situ en trminos de
scf/ton, la adopcin de la nomenclatura (CBM) metano en mantos de carbn. Para
el clculo del almacenamiento de gas total, los volmenes de gas shale in situ se
consideran generalmente por los siguientes trminos:

- Un componente volumtrico,

, involucra hidrocarburos almacenados en

espacios porosos en forma de gas libre, corregido para el volumen ocupado
por la fase sorbida. El volumen de gas libre se cuantifica mediante
modificaciones de mtodos de evaluacin estndar de yacimientos.
- Un componente sorbido,

, con el gas almacenado fsicamente a lo largo

del rea de superficie de los micro y mesoporos. la cantidad de gas
adsorbido generalmente se cuantifica a partir de mediciones de la isoterma
de absorcin, mediante el establecimiento de una isoterma de adsorcin de
equilibrio.
- Un componente volumtrico,

, con presencia de gas disuelto en el

hidrocarburo lquido. Este volumen se combina generalmente con gas
23

adsorbido en yacimientos que contienen gran fraccin de hidrocarburo
lquido en espacios porosos.
- Un componente volumtrico,

, con presencia de gas disuelto en agua de

formacin. La cantidad de gas disuelto se estima a partir de los clculos de
la solubilidad bruta de la roca. Aunque tradicionalmente no se ha
considerado importante, en un estudio reciente se discute el mejoramiento
significativo de la solubilidad de gases en agua formacin cuando est
confinada a pequeos poros (Daz-Campos et al., 2009)
17
.

Por lo tanto, tenemos

el total de gas in situ:

(2. 9)


Donde

se define:

= 32.0368
C

C

(2. 10)

Donde:

: Porosidad

: Saturacin de agua

: Densidad de la roca

: Factor volumtrico del gas

Donde para un sistema de un solo componente, introducimos:


C

= 1.318 x 10
-6
M

) (2. 11)

(2.12)

Donde:

: Volumen de Langmuir
: Presin a la que se quiere evaluar el gas adsorbido en

: Presin de langmuir en
M: Peso molecular del gas

: Densidad de la fase sorbida

24

Donde gas disuelto en aceite es igual a:

(2. 13)

Donde:

: Porosidad

: Relacin de gas en solucin

: Saturacin de aceite

: Densidad de la roca

: Factor volumtrico del petrleo (Vol de yacimiento/Vol de superficie)

Donde el gas disuelto en agua es igual a:

(2.14)

Donde:

: Porosidad

: Relacin de gas en solucin

: Saturacin de agua

: Densidad de la roca

: Factor volumtrico del agua (volumen de yacimiento/volumen de superficie)

En los clculos estndar actuales de la industria, las ecuaciones (2.53) y (2. 14) no
se aplican (GRI, 1997). El gas en solucin en hidrocarburos y agua mvil, y el gas
adsorbido dentro de la materia orgnica se combinan en el anlisis de la isoterma
de adsorcin; por lo tanto, la ecuacin (2. 9) se reduce a:

(2.15)


Sustituyendo las ecuaciones (2. 10) - (2.12) en la ecuacin (2.15) da la ecuacin
del gas in situ para un yacimiento de un sistema para un solo componente:

= 32.0368
C

(2.16)


25

En esta ecuacin, el primer trmino al lado derecho representa la capa de
adsorcin corregido por el volumen del gas libre, mientras que el segundo trmino
corresponde al volumen del gas adsorbido. En esta ecuacin, el factor volumtrico
de formacin, Bg, el peso molecular, M, la densidad de la fase sorbida, s y los
parmetros de adsorcin de Langmuir

son dependientes de la
composicin y las reglas de mezcla simples se recomiendan a menudo para el
clculo de sus valores promedios.

Tambin con el fin de determinar la variabilidad del contenido de gas con la
profundidad se colectan muestras que corresponden en la columna estratigrfica
a la parte inferior, media y superior de cada seccin realizando un anlisis de
desgasificacin de muestras de carbn o de shale que pueden ser efectuadas
empleando canisters los cuales fueron diseados y construidos por el Centro
Internacional de Fsica (CIF). Estos canisters registran en forma automtica el
volumen de gas liberado por las muestras de carbn de igual manera esta
metodologa es utilizada para registrar el volumen de gas shale liberado;

el
propsito de este muestreo es obtener muestras tanto de carbn o shale
saturadas en gas. Igualmente este tipo de muestreo se realiza teniendo en cuenta
la posicin estratigrfica de los mantos de carbn y su grado de exposicin, ya
que es necesario que el frente de la mina se encuentre activo al momento de
muestrear. Se hace importante abarcar gran parte de la columna en cada una de
las reas
16
.


2.5.1 Isoterma de adsorcin

El shale es capaz de almacenar cantidades significativas de gas. El mecanismo
por el cual esto ocurre es llamado adsorcin. En la adsorcin las molculas de una
sustancia se adhieren a la superficie de otra. La adsorcin puede ser visualizada
como un imn adherido a una superficie de metal. Esta es diferente de la
absorcin donde una sustancia es atrapada dentro de otra, como la esponja
remojada por agua. La adsorcin es un proceso reversivo, porque este involucra
fuerzas de atraccin dbiles.

La isoterma de adsorcin de Langmuir asume que el gas adherido a la superficie
de shale o del carbn cubre la superficie como una mono capa de gas. Casi todo
el gas almacenado por adsorcin existe en un condensado, cerca al estado
liquido. A bajas presiones este estado denso permite que volmenes ms grandes
sean almacenados por sorcin tanto como sea posible por compresin.

La formulacin general de la isoterma de Langmuir es:

26

2.5.1.1 Volumen de Langmuir

Este es la mxima cantidad de gas que puede ser adsorbida en una pieza de
shale a presin infinita (Fig. 2.6). Este valor es asintticamente aproximado por la
isoterma a medida que la presin incrementa.

Generalmente, las unidades para el parmetro de volumen de Langmuir

son
scf/ton (volumen de gas por unidad de masa de roca) el volumen puede ser
convertido a scf /

(volumen de gas por volumen unitario de roca) multiplicndolo

por la densidad de la roca.





Figura 2.6 Isoterma de Langmuir tpica mostrando la relacin con el (VL)



2.5.1.2 Presin de Langmuir

Este parmetro afecta la forma de la isoterma. La presin de Langmuir es la
presin a la cual el volumen de Langmuir puede ser adsorbido.



Figura 2.7 Isoterma de Langmuir tpica mostrando la relacin con el (PL)


2.6 Madurez trmica

La madurez trmica proporciona una indicacin de la paleotemperatura mxima
alcanzada por una roca origen. Los hidrocarburos comienzan a ser generados por
Presin
Contenido
de gas V(P)
Presin
Contenido
de gas V(P)
27

encima del umbral de temperatura de 60 C. Vassoyevich et al. (1970) describe
este proceso, donde turbas y lignitos se deshidratan y pierden otros voltiles y el
kergeno se divide en los cuatro tipos distintivos como el "salto de carbonizacin.
La "ventana de petrleo" se encuentra entre las temperaturas de 60-120 C, la
ventana del gas entre 120-150 C. A temperaturas superiores a 150 C, se dice
que la materia orgnica esta posterior a la maduracin y ya no es reactiva para el
desarrollo de hidrocarburos. A temperaturas de 200 C, los compuestos orgnicos
se reducen al grafito y metano. La reflectancia de la vitrinita es el mtodo ms
comn para la determinacin de la madurez trmica, que se completa mediante el
examen microscpico de kergeno o con la montura de la totalidad de la roca y el
registro de la reflectividad de partculas a travs de un fotomultiplicador.


2.6.1 Mtodo de reflectancia de la vitrinita

La Reflectancia de la Vitrinita (VR) es el indicador de maduracin orgnica ms
comn utilizado en la industria del petrleo. McCarteny y Teichmller (1972) lo
describen simplemente como el reflejo de macerales de un haz incidente normal
de luz no polarizada. La reflectividad (R) puede ser registrado como Rv a un
mximo porcentaje (Rv max% ) o Ro en porcentaje (Ro%). Ambos son mediciones
de porcentaje de luz reflejada por la muestra, calibrado frente a un material que
muestra la reflectancia un 100% (es decir, un espejo). El anlisis VR tiene
aceptacin mundial como una tcnica capaz de producir una medicin precisa de
palaeotemperaturas mximas en cuencas con hidrocarburos.

No existe un lmite marcado entre madurez e inmadurez, debido a que materias
orgnicas de composicin diferente responden a velocidades diferentes a los
incrementos de temperatura.
<0.60 Inmaduro
0.60 -1.00 Ventana de Petrleo
1.00 1.40 Ventana de Condensado/Gas Hmedo
>1.40 Ventana de Gas Seco
Gas producible puede ser encontrado acerca de 1.0 %Ro


Figura 2.8 Cambio de Color debido a la Reflectancia de la Vitrinita


28

2.7 Carbn Orgnico Total

El carbono orgnico total (TOC) es una medida de la riqueza orgnica de las rocas
sedimentarias. El anlisis TOC es el primer paso en la seleccin de una formacin
con potencial de hidrocarburos. Los valores de TOC se reportan en porcentaje de
peso de carbono. Por ejemplo, TOC a un valor de un porciento (1 %) en peso
significa que en 100 g de muestra de sedimento, hay 1 g de carbono orgnico.
Hay tres componentes de un valor del TOC, a saber, la materia orgnica extrable
(MOE de carbono), carbono convertible y fraccin de carbono residual. MOE de
carbono es la fraccin de carbono en el aceite y el gas generado pero no
expulsado. El carbono convertible representa el potencial restante de una muestra
de sedimento para generar petrleo y gas. Bajo condiciones de presin y de
temperatura adecuadas, esta parte del TOC se transformara en betn luego en
aceite y finalmente en gas. La fraccin residual de carbono representa el carbono
orgnico que no tiene potencial para producir petrleo o gas debido a su
composicin qumica y estructura, ya que tiene muy poco hidrgeno por unidad de
carbono orgnico.

Los dos mtodos principales en el anlisis para la medicin del TOC son el
mtodo de combustin de Leco y el combinado de la pirolisis de oxidacin, mtodo
de Rock-Eval.

En el mtodo Leco solo se necesita 1 g de muestra. La muestra se tritura y luego
el carbono inorgnico (carbono de carbonato) se elimina por inmersin en HCl
durante 12-16 horas con agitacin intermitente. Posteriormente, la muestra se
enjuaga libre de HCl mediante el uso de agua, papel de filtro y un matraz de
filtracin. El analizador de carbono Leco se calibra con un acero tipo estndar de
contenido de carbono conocido. Despus, el carbono de la muestra se oxida a
dixido de carbono que es detectado por un receptor de infrarrojos (IR) o un
detector de conductividad trmica (TCD). El contenido de carbono total y el
contenido de carbono orgnico total se determinan y la diferencia es el carbono de
carbonato total.


2.7.1 Rock-Eval Pirlisis de la materia orgnica consiste en calentar las muestras
a una velocidad de 25 C / min entre 300 C a 600 C. Cerca de 100 mg de
muestra se mantiene en un horno de pirlisis. La temperatura del horno se
increment primero a 300 C y se mantiene constante. El detector de libre
ionizacin registra el primer pico (llamado S1) correspondiente a la liberacin de
hidrocarburos libres. Despus, la temperatura se aumenta gradualmente a 600 C
a 25 C / min. Esto conduce a la rotura de la materia orgnica o de carbono
convertible obtenindose otro pico (denominada S2). Los picos S1 y S2 se
normalizaron multiplicando por 0,083 para obtener MOE estimado y los valores
convertibles de carbono (0,083 se deriva del promedio en porcentaje en peso de
carbono en hidrocarburos y por la conversin de mg de hidrocarburos por g de
roca o de partes por mil a partes por cien (%)). Por oxidacin de la muestra en un
segundo horno, la fraccin de carbono residual se determina (denominado valor
29

S4 y se convierte en%).
15

La muestra se deja enfriar a temperatura ambiente despus de completada la
pirlisis. Durante el enfriamiento, el

atrapado se libera resultando un tercer

pico S3. El ndice de hidrgeno (HI) est dado por (S2 100/TOC). El ndice de
oxgeno (OI) est dada por (S3 100/TOC). El ndice de hidrgeno es una medida
de la riqueza del hidrgeno de la roca fuente y cuando se traza contra el ndice de
oxgeno (OI) en un diagrama modificado de Van Krevelen (Fig.2.7), el ndice de
hidrogeno (HI) se puede utilizar para determinar el tipo de kergeno presente.


Figura 2.9 Diagrama de Van Krevelen
15

El valor de Tmax es un parmetro normalizado, calculado a partir de la
temperatura a la que el pico S2 alcanza su mximo y proporciona una medida de
madurez trmica de la materia orgnica, Espitali et al. (1977). Cualquier muestra
con Tmax inferior a 430

C se considera trmicamente inmadura. Las muestras con

Tmax entre 430 - 450C se considerar trmicamente madura y en la ventana de
aceite o gas hmedo y las muestras con Tmax mayor de 450C se consideran en
la ventana de gas, Emeis y Kvenvolden (1986).

Liro et al. (1994) demostr variabilidad orgnica a travs de las diferencias en el
TOC y tipos de kergeno dentro de intervalos transgresivos y regresivos del shale.
El explico que estos shales son capaces de generar grandes cantidades de
hidrocarburo lquido y el carcter del kergeno tambin sugiere la presencia de
material propenso al crudo. Si el contenido del TOC vara desde 1 hasta 10% en
peso, implica una roca fuente de alta calidad segn Popov et al. (2000).

La Tabla. 2.1 propuesta por Jarvie (1991), muestra los valores de TOC
comnmente utilizados para shales y carbonatos.

30


INTERPRETACIONES CLASICAS DEL CONTENIDO DE TOC EN LOS
SEDIMENTOS
(BASADA EN LA MADUREZ TEMPRANA DE LA VENTANA DE ACEITE)
(a)
CUALIDAD
TOC EN SHALES (wt%) TOC EN CARBONATOS (wt%)

POBRE 0.0-0.5 0.0-0.2
MEDIO 0.5-1.0 0.2-0.5
BUENO 1.0-2.0 0.5-1.0
MUY BUENO 2.0-5.0 1.0-2.0
EXCELENTE >5.0 >2.0


(b)
POTENCIAL
TOC (wt%)

DE GENERACION

INSUFICIENTE 0.0-0.5
MARGINAL 0.5-1.0
SUFICIENTE >1.0
Tabla. 2.1 Valores comunes de TOC de shales y carbonatos


2.8 Presin del yacimiento y transiente de presin en el yacimiento

Los pozos de gas shale requieren una gran longitud lateral horizontal que
intercepte mltiples fracturas hidrulicas para ser econmicamente viables. Debido
a la naturaleza de una permeabilidad ultrabaja de la matriz del shale, se
necesitara un perodo significativo de tiempo para un pozo productor en
regmenes transitorios de flujo. Por lo tanto, las caractersticas del transiente de
presin para un pozo de gas son considerablemente importantes como en la
evaluacin de tratamiento de fractura mediante la estimacin de parmetros de
fractura y el yacimiento y la prediccin del comportamiento de la produccin a
largo plazo de los pozos en la recuperacin de gas.

La respuesta del transiente de presin para un pozo horizontal de gas shale con
mltiples fracturas transversales es controlada por una combinacin especfica de
propiedades del yacimiento y la configuracin del pozo. La investigacin de las
caractersticas de respuesta del transiente de presin es virtud de una serie de
factores y mecanismos de flujo. Los factores que se consideran que se incluyen
son la permeabilidad de la matriz, la conductividad de las fracturas hidrulicas, el
espaciamiento del clster, y el tamao / aumento de la permeabilidad en la zona
estimulada. Los impactos del gas de desorcin y conductividad de la fractura
dependiente del esfuerzo en el comportamiento del transiente de presin. Este tipo
de estudio proporciona informacin detallada sobre las caractersticas del
transiente presin durante la produccin de pozos de gas shale.
31


Las caractersticas del transiente de presin para los yacimientos de gas,
necesitan ser analizados en trminos de pseudopresin (Lee, et al. 2003). La
pseudopresin se define convencionalmente por (Al-Hussainy, et al 1966.):

= 2

(2.18)

Y la cada de pseudopresin, se calcula as:

A

(2.19)

Donde

es la pseudopresin de la presin inicial del yacimiento.

Basado en la presin del fondo de pozo simulada, la pseudopresin es la primera
que se calcula. Luego, la respuesta obtenida puede ser derivada y ambas
respuestas pseudopresin y funcin de derivada, estn en funcin del tiempo en
un grfico log-log para la identificacin de los regmenes de flujo.

De las respuestas del transiente de presin, cinco regmenes de flujo distintos
pueden ser identificados. Una imagen ilustra estos cinco regmenes (Fig.2.10) y
otra imagen presenta la respuesta del transiente de presin (Fig.2.11). La
pseudopresin est indicada con marcadores rellenados y la derivada de la
pseudopresin est indicada con marcadores abiertos.

La primera es una fractura radial / Formacin de flujo lineal (RF / LF),
caracterizada con las respuestas de derivadas planas (Larsen y Hegre, 1994).
Durante este rgimen de flujo, el flujo de gas en la fractura es convergente
radialmente dentro del pozo horizontal mientras el gas fluye dentro de la fractura
linealmente en la formacin (Fig. 2.10a). Este rgimen es corto, slo dura
alrededor de una hora. Luego resulta una pendiente de 1/4 de derivada, lo que
indica un flujo bilineal (Fig. 2.10b). Su duracin es de ms de un da (entre 1 y 30
horas), es interesante notar que mientras la conductividad de la fractura
adimensional es de aproximadamente 67 aos, el flujo bilineal est presente
durante un perodo de tiempo relativamente considerable. La conductividad de la
fractura adimensional se define como:

(2.20)

Donde k
f
: Permeabilidad de fractura w
f
: Ancho de fractura

: Permeabilidad de la formacin. X
f
: Longitud de fractura


Despus del BL, una formacin de flujo lineal se produce con una pendiente de
0.5 de derivada hasta 4x10
3
horas. El flujo va hacia la fractura linealmente al
interior del yacimiento estimulado entre las fracturas (Fig. 2.11c). En la figura, el
32

interior del yacimiento se esboza en el tablero rectangular como pesado y en el
exterior del yacimiento est entre los tableros rectangulares pesados y ligeros.
LF es seguido por un rgimen nico de flujo llamado Flujo cuasi -estable (Fig.
2.11d). Teniendo en cuenta que en el grfico de la derivada, exhibe una lnea
recta con una pendiente de 0.88, durando mucho tiempo desde 4x10
3
horas hasta
5x10
4
horas. Durante este rgimen de flujo, la interferencia de presin entre
fracturas domina mientras empieza el flujo a lo largo de la barrera entre el interior
y el exterior del yacimiento pero es insignificante. Este es un estado de flujo
pseudo estable aproximado, con el cual el interior del yacimiento es depletado con
contribucin limitada desde la parte exterior del yacimiento. Esta puede ser la
razn de porque la pendiente de la derivada no es uno pero se aproxima a este
valor.

El ltimo rgimen indefinible es el flujo linear a tiempo tardo exterior con una
pendiente de 0.5 como respuesta a la derivada, durante el cual un flujo linear
domina desde la parte exterior del yacimiento hasta la parte interior mientras la
parte interior del yacimiento ha sido depletada. Este empieza aproximadamente
alrededor de las 2x10
5
horas. Hay una larga transicin entre este flujo linear tardo
y el previo estado quasi estable de flujo. Es interesante notar que el flujo
dominante de frontera no se desarrolla al final de la simulacin (50 aos). Esto
resplada la interpretacin comn que para yacimientos de shale, el drenaje de la
parte exterior del yacimiento se limita y el mayor drenaje para un yacimiento
horizontal fracturado viene de la parte interior del yacimiento entre las fracturas
hidrulicas, tambin por lo general se refiere como el volumen de yacimiento
estimulado (Maxwell, 2006 y Ozkan, 2009.)
18


Figura 2.10 Cinco regmenes de flujo tpicos para pozos horizontales de gas shale
con fracturas hidrulicas.
33




Figura 2.11 Respuestas de la cada de presin y de la derivada.


2.9 Descripcin de Corazones

La composicin de la matriz de las shale es bastante variable haciendo
complicado obtener la densidad del grano.

Es casi imposible calcular con precisin la saturacin de gas y agua a partir de la
informacin de registros a menos de que sea calibrada con informacin de
ncleos. El efecto de arcilla en la porosidad neutrn usualmente elimina los
efectos de gas en el registro.

Es importante conocer con certeza la porosidad en este tipo de yacimientos, ya
que estos pueden afectar significativamente el clculo de reservas (OGIP).


34



Figura 2.12. Diagrama de Referencia de Niveles de Porosidad en un yacimiento de
Gas Shale.

2.9.1 Mtodos Convencionales para la medicin de porosidad

2.9.1.1 Picnmetro de Baja Presin (LPP)

Es uno de los aparatos ms confiables para la medicin de porosidad en lutitas.
Este equipo usa la ley de Boyle para medir el volumen de grano en las muestras
trituradas mediante la variacin de los cambios en la presin de la celda llena de
helio, este se expande en una celda de calibrada previamente que contiene la
muestra. La muestra es colocada en una celda de volumen conocido. El helio
tambin permanece en otra celda de referencia de volumen y presin conocidas
(20 psi). Las dos celdas estn conectadas, causando una cada de presin en una
celda a medida que se llenan los espacios vacos en la otra que contiene la
muestra. El volumen que no es llenado por el helio es aquel ocupado por el
volumen de los granos de la muestra y los poros no-interconectados. El Volumen
de los granos puede entonces ser calculado teniendo en cuenta las presiones
antes y despus de conectar las celdas. Este equipo se calibra con esferas de
acero de volumen conocido, cada vez que se va a realizar una medicin.

35



Figura 2.13 Cilindro para la trituracin de la muestra de Shale.




Figura 2.14 Picnmetro de Baja Presin LPP (Low Pressure Pycnometer)


El volumen total de la muestra con el volumen de grano obtenido del LPP, se usa
para calcular la porosidad, que generalmente es llamada como porosidad total o la
porosidad de la ley de Boyle Ecuacin (2.90).


b
g b
b
p
V
V V
V
V
= = |
(2.21)




36

CAPITULO III EVALUACION DE PROPIEDADES DE
YACIMIENTOS DE GAS SHALE

Para evaluar cada una de las propiedades cuantitativas del yacimiento, los
ingenieros se valen de ciertas herramientas que les permiten determinar la
viabilidad de ste. En este estudio, se utilizaron programas para simular el
comportamiento del yacimiento y la relacin de sus propiedades ms
fundamentales entre s y en conjunto con la produccin. El paquete de
simuladores usados, fue brindado por la empresa Fekete que es la encargada de
llevar a cabo muchos proyectos en Norteamrica desde 1973 (prcticos y fciles
de usar, con soluciones avanzadas a los proyectos de optimizacin de produccin
de recuperacin y de ingeniera de yacimientos.)

De los 8 productos que maneja ste software, ste estudio se construy con el
F.A.S.T CBM y el Fekete Harmony. Estos son los encargados del desarrollo y
anlisis de yacimientos no convencionales de mantos de carbn. Para ste caso,
un yacimiento de gas shale tiene caractersticas similares a las de un yacimiento
de mantos de carbn (un sistema de doble porosidad fractura-matriz, un sistema
geomecnico, un sistema de adsorcin, as como un sistema de red poral, entre
otros), lo que hizo que se utilizara un software de yacimientos de mantos de
carbn para yacimientos de gas shale (puesto que la universidad no cuenta con
sta herramienta).

El CBM incluye herramientas de anlisis de datos de produccin para modelos
analticos y numricos, puede realizar pronsticos de produccin de gas/agua un
cotejo histrico si es requerido y estimar ciertas propiedades del yacimiento; las
herramientas que incluye para realizar estos clculos son: volumtricas,
herramientas de declinacin convencional, balance de materia dinmico y esttico
y curvas tipo de Agarwal-Gardner (para mantos secos). Tambin realiza anlisis
de riesgo usando simulaciones de Monte Carlo. Su mdulo de monitoreo de
produccin puede ser usado para monitorear la produccin de nuevos pozos.
Cuenta as mismo con varios mtodos para calcular la permeabilidad relativa y los
efectos del encogimiento de la matriz y pueden ser fcilmente integrados al
anlisis.

El software de Fekete Harmony sirve para evaluaciones de reservas,
modelamiento y caracterizacin del yacimiento, adquisiciones y cesiones del
campo, desempeo del monitoreo y de la optimizacin, planeamiento del
desarrollo del campo, yacimientos no convencionales y convencionales.





37

3.1 F.A.S.T. CBM

Como se mencion anteriormente se adapt el simulador para yacimientos de
mantos de carbn a yacimientos de gas shale. La manera en que ste software
evala las propiedades permite conocer un poco ms del comportamiento del
yacimiento y de su relacin con las propiedades fundamentales que determinan el
inicio o no de un proyecto en un campo.

Ahora, se va a realizar un marco de la evaluacin. Cules son los pasos a seguir y
a tener en cuenta para realizar una buena aproximacin de las propiedades del
yacimiento de gas no convencional.


3.1.1 Clculos Volumtricos

Para los clculos volumtricos se tiene en cuenta una ecuacin ajustada por los
parmetros de adsorcin de Langmuir (VL y PL) de la ecuacin de gas para un
yacimiento convencional.

(3.1)


Si se ve con atencin la ecuacin anterior, se puede apreciar que la porcin de
volumen de gas libre del yacimiento convencional (el 2 trmino), se encuentra
presente en el clculo de volumen de gas almacenado en un yacimiento no
convencional. Es la fraccin de gas adsorbida (el 1
er
trmino) la que determina el
volumen de gas adicional, que se encuentra presente en los yacimientos con una
alta capacidad de sorcin (no convencionales). Los estudios
22
han mostrado que
la industria del gas shale estaba determinando incorrectamente los volmenes
para determinar el gas in situ en los yacimientos con una alta capacidad de
sorcin, al no contar apropiadamente con el volumen ocupado por la fase
adsorbida. Entonces una porcin del volumen total poroso sera ocupado por el
gas adsorbido y no estara disponible para molculas de gas libre.

(3.2)

(3.3)


El modelo de Langmuir es el ms ampliamente usado en la industria petrolera y en
varias ocasiones anteriores se compararon sus clculos con los de un modelo
termodinmicamente consistente
22
. Sin embargo, las conclusiones de los
38

experimentos fueron muy positivas para el modelo termodinmicamente dbil de
Langmuir pues el clculo de volumen de gas in situ variaba menos de un 1% con
respecto al modelo experimental.

El modelo de Langmuir aunque no es muy preciso para los clculos de gas in situ
(como el termodinmicamente riguroso IAS) es el empleado por el software,
Qu hacemos con la correccin de Ambrose?. De la misma manera, no siempre
se tiene un sistema de gas adsorbido de metano puro. El shale, as como el
carbn, pueden tambin adsorber cantidades apreciables de CO
2
, N
2
e
hidrocarburos ms pesados (etano, propano, entre otros). Cada gas no es sorbido
independientemente sino que compite por los mismos sitios de sorcin. Debido a
sta competencia la sorcin de cada componente es menor que cuando los
gases actan independientemente. Es de sta manera cmo una isoterma multi-
componente es necesaria para predecir la composicin de gas producido, el gas in
situ y las reservas. Esta isoterma multi-composicional es requerida para el recobro
primario al disminuir la presin y especialmente necesaria para recobro secundario
al inyectar CO
2
/ N
2
.

(3.4)

(3.5)

El software cuenta con unas suposiciones para un sistema binario en el clculo de
gas in situ. Para implementar el modelo de Langmuir se dice que el volumen de
gas in situ que sale es metano puro una composicin metano-CO
2
, que el gas
adsorbido est siempre en equilibrio con la fase de gas libre, la presin inicial y la
fraccin molar del gas sorbido es conocida y la isoterma composicional de
Langmuir es usada para calcular la composicin del gas desorbido. Sin embargo
los estudios realizados por Ambrose et al. 2010, nos muestra que el porcentaje
sorbido del gas libre cambia significativamente dependiendo de qu tipo de
isoterma es usada (composicional o simple). La diferencia entre la isoterma de gas
libre de un slo componente a la de n gas multicomposicional cambia cerca de un
20%, lo que lo hace muy importante si se tienen que determinar las reservas
recuperables en yacimientos donde una gran cantidad de gas contenida en el
hidrocarburo con mayor peso molecular (gas mojado). De la misma manera, un
gran impacto es observado, debido a que el gas sorbido requiere presiones de
abandono menores, para que puedan ser recuperados.


39

3.1.2 Cambio de la Matriz y la Relacin Porosidad / Permeabilidad. Porosidad
de fractura y su relacin con el encogimiento de la matriz.

As como el carbn, el shale cuenta con una matriz y un rea estimulada
fracturada llamada cleat, a travs de los cuales se mueve el gas. El gas es
almacenado por adsorcin dentro de la matriz y el rea estimulada, y el flujo de
ste se realiza slo cuando el agua ha salido de la roca y se alcanza la presin de
desorcin (proceso de difusin). Ahora se va a hablar acerca de la relacin de la
permeabilidad con la porosidad y el encogimiento e hinchamiento de la matriz.

3.1.2.1 Seiddle y Huitt

- El hinchamiento es inducido en la roca por el proceso de adsorcin, sin
embargo adems de ste fenmeno la roca experimenta otra deformacin por las
fuerzas mecnicas a la que es sometida (m y p respectivamente). Se ve cmo
se relacionan el coeficiente de hinchamiento de la matriz (Cm) con el volumen y la
presin de Langmuir (VL y PL). La ecuacin para el coeficiente de acoplamiento
mecnico (Cp que representa al comportamiento de la matriz sobre cargas) es
anloga a la de Cm.




(3.6)

(3.7)

(3.8)

mi y pi tienen que ser multiplicado por 10
6
para llegar a ser una cantidad
significativa, por lo que a condiciones normales (atmosfricas) es igual a 0.
Entonces:

(3.9)

(3.10)

Si el coeficiente de deformacin neto estuviera determinado por una relacin entre
los coeficientes mecnicos y de adsorcin sera:

40

(3.11)

Para la relacin de la porosidad con la permeabilidad, se tiene que:

(3.12)

Y luego de operar se tiene:

(3.13)

La porosidad se relaciona de una manera similar con la presin por medio de:

(3.14)

Que luego se convierte en

(3.15)

Al aplicar logaritmo natural y dividir las ecuaciones 3.13 y 3.15 se tiene que:

(3.16)

Si

= n, entonces:

(3.17)


- Como se dijo anteriormente, los cambios en la matriz estn relacionados
directamente con la capacidad contenida de gas adsorbido de sta y no con la
presin del gas. La geometra de una roca no convencional puede ser
representada como se muestra en la figura 3.1, los esfuerzos correspondientes a
los cambios en la matriz () estn representados por:



41


Figura 3.1 Geometra de la roca: matriz y fractura
(a y w respectivamente)

y si se iguala con la ecuacin 3.7 se tiene:

(3.18)

(3.19)

Para la porosidad de la fractura se tiene
23
:

(3.20)

Donde la amplitud de fractura y de matriz a una presin dada estn denotadas por:

Entonces una relacin entre la porosidad de la fractura actual y la inicial sera:

, al aplicar la ley de la oreja y separacin de trminos se tiene:

42

, se reemplaza la ecuacin 3.18 y:

, el ltimo trmino es muy

pequeo ltimo y reagrupando se tiene:

, ahora se reemplaza la ecuacin 3.20:

, Por ltimo se utiliza la ecuacin 3.19 para

despejar y se obtiene:

(3.21)


3.1.2.2 Palmer y Mansoori

Durante la produccin primaria de un yacimiento, los esfuerzos efectivos
incrementan y la permeabilidad disminuye debido a la compresin de las fracturas.
Sin embargo, en los yacimientos no convencionales (CBM y GS) vemos que el
drawdown lleva a produccin de metano, acompaado con encogimiento de la
matriz, que abre las fracturas e incrementa la permeabilidad. Los 2 efectos que
ocurren durante una cada de presin (compresin de fracturas y encogimiento de
la matriz), actan en direcciones opuestas a la permeabilidad.

La permeabilidad est en funcin de los esfuerzos y de la presin del poro, como
se menciona en diversos estudios
24
. De sta manera se parte de la ecuacin de
elasticidad lineal para cambios de esfuerzos en rocas porosas y se obtiene la
siguiente ecuacin:

Donde podemos observar cmo el esfuerzo incremental del volumen poroso
(adimensional) se relaciona con los esfuerzos de la roca y del grano,
respectivamente. Como se dijo que se iba a partir de la ecuacin de elasticidad
43

lineal (con un cambio de esfuerzo muy pequeo) y bajas condiciones de esfuerzos
en una direccin, el trmino

(Higgs et al 1994) lleva a la siguiente ecuacin:

(3.22)

Donde dS, dP y dT cambio en esfuerzo de overburden, presin de poro y de
temperatura respectivamente. Tambin f, , M, K y se refieren a una fraccin
que vara de 0 a 1dependiendo del tipo de roca
3
, la compresibilidad del grano, el
mdulo axial, el mdulo bruto y la expansin termal del grano respectivamente.

El mdulo de Young (E) y la razn de Poisson (v) estn directamente relacionados
con M y K a travs de la teora de elasticidad isotrpica:

(3.23)

(3.24)

Como en nuestro sistema de roca (gas shale o mantos de carbn) la porosidad es
tan pequea y el esfuerzo de overburden no cambia significativamente, la
ecuacin anterior puede ser reacomodada como:

Como dT es un trmino de expansin y contraccin, si disminuye, el tejido se
encoge y las fracturas aumentan en amplitud, lo que la hace anloga al proceso
de desorcin y al encogimiento de la matriz por lo que se puede escribir de la
siguiente manera:

Como el trmino entre parntesis puede expresarse por medio de mediciones
volumtricas de expansin en:

En la literatura sta b tiene un valor experimental de 0.02 psi
-1
, sin embargo tanto
el trmino b como el trmino deben ser determinados en un laboratorio
mediante pruebas mecnicas y de desorcin del metano en corazones
2
.

44

Como M y K son independientes de la presin y haciendo el reemplace de los
trminos anteriores, la ecuacin 3.22 se convierte en:

(3.25)

Donde Cm se refiere a:

Ahora al integrar la ecuacin 3.23 y al dividirla por

se obtiene finalmente:

(3.26)

La roca tiende gradualmente a volverse ms dura durante el drawdown (E suele
incrementar) debido a que las fracturas no tienen la capacidad de cerrarse en
rugosidades, rocas finas o mineralizacin. La falla en estas rocas no
convencionales (carbn o gas shale) puede incrementar la recuperacin de la
permeabilidad. Cuando las condiciones son favorables, el drawdown puede inducir
a una falla debido a los cambios en esfuerzos producidos por el encogimiento de
la matriz. La falla puede llevar a un incremento de la permeabilidad debido al
fenmeno conocido como dilactancia.

Este modelo puede ser comparado con el de Seidle y Huitt . Este (Seidle) tiene
una constante igual a 2 en vez de (K/M 1) en el trmino de encogimiento. Esto
significa que su trmino de encogimiento es cerca de 2.6
3
veces mayor que el de
la ecuacin de Palmer y Mansoori. Seidle y Huitt calculan la relacin

separadamente para los efectos de los esfuerzos y de encogimiento, luego son
multiplicados para obtener el efecto combinado. En cambio en el modelo de P&M
son agregados y no multiplicados.

Finalmente se debe tener en cuenta que la teora de Seidle asume que la
compresibilidad del volumen poroso es constante con los cambios de la presin de
poro, pero como se ha demostrado
3
esto no es cierto.

Sin embargo la nueva teora se aplica slo para cambios pequeos en los
esfuerzos como muestra la ecuacin linear de elasticidad. Debido a la porosidad
de la roca (que es muy pequea en nuestro caso), los cambios en la porosidad se
ven reflejados hasta en un factor de 2, en cambio los cambios correspondientes de
permeabilidad cambian menos de un factor de 10.

3.2 Presin de flujo variable

En el desarrollo inicial de esta tcnica de interpretacin, la presin del fluido se
asume constante. Un caso adicional de presin de flujo variable tambin se puede
45

considerar cuando al interpretar el grfico de la raz cuadrada del tiempo. El eje
estndar Y de tasa recproca es simplemente sustituir con la presin normalizada
para notar los cambios en la presin de flujo durante la parte del flujo lineal de los
datos. El Flujo lineal transitorio que parece ser una lnea recta en la presin
normalizada vs el grfico de la raz cuadrada del tiempo. La pendiente de la lnea
recta, m, se utiliza para calcular x
f
, as:

(3.27)

Como se dijo anteriormente, la observacin ha sido que la mayora de los pozos
no convencionales producen una alta declinacin para maximizar el recobro.
Debido a esta prctica comn, el enfoque del pronstico en el mdulo de presin
variable utiliza una disminucin hiperblica con la suposicin de una presin
constante para la duracin del pronstico.


3.3 Teora del anlisis del transiente curvas tipo

El mtodo de anlisis transitorio curvas tipo ofrece una perspectiva alterna que es
ideal para el anlisis de perodos muy cortos de produccin (temprano), y / o el
anlisis de yacimientos de muy baja permeabilidad.

En los mtodos Blasingame, Agarwal-Gardner y la integral de presin
normalizada, las curvas tipo se escalan de manera que hay convergencia para una
pendiente unitaria (pendiente unitaria). Esto se logra mediante el uso de un tiempo
adimensional que se basa en el rea (T
DA
o T
DD
). Una consecuencia de este tipo
de escala es que hay numerosos tramos de transientes, Si un tiempo
adimensional (T
D
) basado en radio de un pozo es escogido en su lugar, habr una
solo tramo de transiente con una serie de curvas de lmite dominada. Al ver una
nica cuerva tipo (por ejemplo: r
eD
= 28) no hay diferencia entre estos dos
formatos de escala (ver figura abajo). Sin embargo, durante la visualizacin de
todos las curvas tipo, la presentacin del transiente proporciona una base ms
conveniente para el anlisis de los datos del transiente. Esto se deduce del hecho
de que la seleccin de un tramo de transiente apropiada no es necesaria.

Considere la posibilidad de un conjunto de datos de produccin transitoria
(permeabilidad muy baja o vida de produccin limitada). Una comparacin de la
ajuste de la curva tipo utilizando el anlisis tasa-tiempo de Agarwal-Gardner con la
del formato del transiente se muestra a continuacin.

En este ejemplo la presentacin del transiente las curvas tipo proporciona un
mayor ajuste de las curvas tipo. Por la misma razn que los datos del transiente
funciona mejor con el formato del transiente (Q
D
vs t
D
), hay que sealar que lmite
dominado por anlisis de flujo no es aconsejable utilizarlo en este mtodo
25
.

46


Figura 3.2 Formatos de escala de cuerva tipo para la seleccin del mejor tramo


Figura 3.3 Comparacin de la ajuste de la curva tipo utilizando el anlisis tasa-
tiempo de Agarwal-Gardner con la del formato del transiente

Preparacin de datos

En el eje horizontal se da el tiempo de balance de materiale (pseudo-tiempo para
el gas) y en el eje vertical la Tasa Normalizada, estas dos variables se
representan:

Tasa Normalizada vs Tiempo Balance de Materia

Tasa Normalizada

Pozos de Gas

(3.28)

47

Pozos de Gas

|
c
A

c (

)
|
(3.29)

En el modificado (Anlisis de la Tasa del transiente RTA) grafico, el 1 / DER
i
(el
inverso de la derivada integral de presin) tambin se incluye.


Anlisis

La tasa normalizada se representa frente al tiempo de balance de materia en una
escala log-log del mismo tamao que las curvas tipo. Este grafico se llama
"grafico de datos". Cualquier unidad conveniente puede ser usada para normalizar
la rata o el tiempo debido a que un cambio en unidades simplemente causa una
variante a la escala logartmica de los datos sin procesar. Se recomienda graficar
tasas operadas diarias y no tasas mensuales; especialmente cuando se analizan
datos transitorios.

La grafica de los datos se mueve sobre el grafico curva tipo, mientras que los ejes
de los dos grficos se mantienen paralelos, hasta obtener un buen ajuste. Existen
distintos curvas tipo al que se le debe intentar obtener el mejor ajuste para todos
los datos. La curva tipo que mejor se ajuste a los datos debe ser seleccionada y su
"r
e
/ r
wa
" (r
e
/ x
f
para la fractura) es un valor incluido.

El anlisis de las curvas tipo es hecho al seleccionar el punto de partida, y leyendo
sus coordenadas del grfico (q / p y tc)
ajustada
, y del grfico de curvas tipo (q
D
y
t
D
)
ajustada
. Al mismo tiempo, el valor principal "r
e
/r
wa
" (r
e
/ x
f
para la fractura de las
curvas tipo) de la curva ajustada es mostrado

Para crear un pronstico, la curva tipo seleccionada se remonta al grafico de
datos, y extrapolada. La tasa de futuro se lee del grafico de datos, de la curva tipo
trazada.

3.3.1 Clculo de los parmetros: Curvas tipo radiales
Pozos de Gas

La permeabilidad se obtiene reordenando la definicin de tasa adimensional:

(3.30)

(3.31)

48

(3.32)

La solucin para radio aparente del pozo viene de las definiciones del tiempo
adimensional y permeabilidad:

)(

(3.33)

Factor de dao se calcula as:

S = ln(

) (3.34)

Los parmetros de volumen y del rea se calculan as:

G =
t

)

(

(3.35)

A =

(3.36)

G : Volumen, Bcf
: rea, acres

3.3.2 Clculo de los parmetros: Curvas tipo para Fractura de conductividad
finita

Las curvas tipo para fractura de conductividad finita se basan en un yacimiento de
forma cuadrada con una fractura hidrulica en el centro. stos se representan
usando el formato de tiempo adimensional, y por lo tanto se clasifican como
curvas de formato transiente. Los tramos de transiente representan tres diferentes
conductividades fractura adimensional (FCD) valores (0,5, 5, y 500). Las
derivaciones de la barrera dominada diverge en cada una de las 4 longitudes del
yacimiento adimensional (1, 2, 5, 25) x
eD

49


Figura 3.4 Curvas tipo para fractura de conductividad finita



Conductividad de la fractura adimensional:

(3.37)

Longitud/ anchura del yacimiento adimensional:

(3.38)

Donde 2X
e
= dimensin principal del yacimiento


Pozos de Gas

La permeabilidad se obtiene reordenando la definicin de tasa adimensional:

(3.39)

(3.40)

(3.41)


50

La solucin para longitud media de la fractura viene de las definiciones del tiempo
adimensional y permeabilidad:

)(

(3.42)

Los parmetros de volumen y del rea se calculan as:

G =

)

(

(3.43)

A =

(3.36)

Donde:
G : Volumen, Bcf
: rea, acres

3.4 Presin Capilares
3.4.1 Mtodo Analtico - Sistema gas-crudo

(3.44)

(3.45)

Donde:

3.4.2 Brooks Corey

El modelo de Brooks Corey es ampliamente aceptado para el clculo de la presin
capilar para los clculos de la permeabilidad relativa. Este modelo es una
representacin modificada del modelo de Corey y representa la presin capilar en
una forma ms general.

Sistema gas-crudo

(3.46)

(3.47)

(3.48)

51

Donde:
= ndice de distribucin del tamao del poro

= saturacin de agua irreducible

= saturacin de petrleo residual a gas

3.5 Contenido de gas

3.5.1 Isoterma de adsorcin de Langmuir

La matriz del los sistemas de yacimientos de Shale tienen una inmensa capacidad
para el almacenamiento de metano. El mecanismo por el cual ocurre se denomina
adsorcin. En la adsorcin, las molculas de gas se adhieren a la superficie de
carbn o de material orgnico en el caso del Shale.

Tpicamente, los depsitos de Shale puede almacenar mucho mas gas por
adsorcin que lo que depsitos convencionales puede contener por compresin a
presiones inferiores a 1000 psia. Dado que el volumen de un sistema de fractura
es pequeo en comparacin con el volumen del yacimiento, el gas libre slo
corresponde a una pequea parte del gas almacenado en el Shale. Como
resultado, la relacin entre el volumen y la presin se describe a menudo
solamente por la isoterma de desorcin.

3.5.1.1 Formulacin de la isoterma de Langmuir

La liberacin de gas adsorbido se describe comnmente por una relacin de la
presin llamada la isoterma de Langmuir. La isoterma de adsorcin Langmuir
asume que el gas se adhiere a la superficie del Shale, y cubre la superficie como
una capa nica de gas (una monocapa). A bajas presiones, este estado denso
permite mayores volmenes de almacenar por sorcin que los que son posibles
por compresin.

La formulacin tpica de isoterma de Langmuir es:

(3.49)

Donde Cg es el contenido de gas medido en scf / ton de Shale.


3.5.2 Clculos de yacimientos de gas

Histricamente, en un yacimiento de gas, solo se consideraba gas in situ libre.
Debido a esto, slo un nombre se requera: OGIP. Sin embargo, con el creciente
uso de yacimientos de gas adsorbido en la industria, Fekete ha adoptado el
nombre de "OGIP
F
" para definir el gas in situ para un yacimiento de gas libre. Del
52

mismo modo, el nombre de "OGIP
A
" se utiliza para definir el gas in situ en un
yacimiento de gas adsorbido. El nombre OGIP ha sido mantenido para describir el
total original de gas in situ.

3.5.2.1 Ecuacin de gas libre

Clculo de Gas in situ libre original (OGIP
F
)

C (

(3.50)
Donde:
A : rea, acres
: Espesor de la formacin, ft
C : Porosidad, %

: Saturacin de gas inicial, %

: Factor volumtrico,

3.5.2.2 Ecuacin de gas adsorbido Yacimientos de Shale

Los yacimientos de gas shale suelen contener mucho ms gas adsorbido que gas
libre. Por lo tanto, los clculos OGIP de yacimientos de shale tambin se deben
tener en cuenta para la adsorcin. Para los yacimientos de shale, la siguiente
ecuacin se usa para calcular Gas in situ adsorbido original (OGIP
A
).

Clculo de Gas in situ adsorbido original (OGIPA).

(3.51)

Donde:
A : rea, acres
: Espesor de la formacin, ft

: Densidad del shale en

p : Presin a la que se quiere evaluar el gas adsorbido en

: Presin de langmuir en

: Volumen de Langmuir


3.6 Madurez trmica

La madurez trmica es proceso qumico dependiente de la Temperatura, que
ocurre en las cuencas sedimentarias como consecuencia natural de la alteracin
trmica de la materia orgnica depositada con los sedimentos. La extensin de la
reaccin (cantidad de petrleo generado) es funcin de la temperatura a la cual
ocurre la reaccin y el tiempo, en este caso vamos a evaluar el grado de madurez
de forma cualitativa ya que el software Fekete Harmony no predice el grado de
53

madurez trmica teniendo en cuenta los rangos de temperatura para la generacin
de la ventana de petrleo y gas; La "ventana de petrleo" se encuentra entre las
temperaturas de 60-120 C, la "ventana del gas" entre 120-150 C. A
temperaturas superiores a 150 C, se dice que la materia orgnica esta posterior
a la maduracin y ya no es reactiva para el desarrollo de hidrocarburos. A
temperaturas de 200 C, los compuestos orgnicos se reducen al grafito y
metano.

Como se dijo en el captulo II la reflectancia de la vitrinita es el mtodo ms
comn para la determinacin de la madurez trmica, que se completa mediante el
examen microscpico de kergeno o con la montura de la totalidad de la roca y el
registro de la reflectividad de partculas a travs de un fotomultiplicador

Debido que no existe un lmite marcado entre madurez e inmadurez el software
Fekete Harmony no lo evala, debido a que materias orgnicas de composicin
diferente responden a velocidades diferentes a los incrementos de temperatura es
por eso que se sigue el criterio de una evaluacin cualitativa para esta propiedad
del shale.

Grado de Madurez
<0.60 Inmaduro
0.60 -1.00 Ventana de Petrleo
1.00 1.40 Ventana de Condensado/Gas Hmedo
>1.40 Ventana de Gas Seco
Gas producible puede ser encontrado acerca de 1.0 %Ro
Tabla 3.1 Grado de madurez


3.7 Carbn Orgnico Total - TOC

En el estudio del potencial de generacin de las rocas generadoras de
hidrocarburos, se evala el contenido total de carbono orgnico (TOC) como
tambin otros parmetros evaluados en esta tesis

El contenido total de carbono orgnico (TOC) debe ser como mnimo de 2% en
peso para que una roca pueda generar cantidades comerciales de
hidrocarburos.
26
El TOC se evala en el laboratorio mediante pirolisis es por esto
que la descripcin del contenido total de carbono orgnico (TOC) en este
documento se hizo de acuerdo con la clasificacin propuesta por Jarvie (1991),
que muestra los valores de TOC comnmente utilizados para shales (tabla 2.1)

3.8 Descripcin de Corazones

La Descripcin de Corazones capta visualmente el brillo del carbn, las bandas, el
espaciamiento de las fracturas, la mineraloga, el espesor de shale, entre otros
54

factores. Proporciona datos acerca de la composicin, la permeabilidad y la
heterogeneidad de una secuencia de inters. Los datos de Corazones de shale
son obtenidos de estudios realizados a pozos perforados en Colombia y de la
formacin Marcellus en Estados Unidos clasificndolos cualitativamente que tan
buenos y malos puede ser estos yacimientos, teniendo tambin como ayuda el
software FEKETE que nos la opcin de variar ciertos parmetros para estimar
reservas y rentabilidad en un yacimiento de Gas Shale.




































55

CAPITULO IV - EVALUACION ESPECIFICA DE LAS
PROPIEDADES DEL POZO USCO-69

El inters en las nuevas oportunidades de desarrollo en Yacimientos No-
Convencionales, muestran a Colombia con un gran potencial de gas shale de lo
cual hace atractivo en nuestro pas una mayor investigacin y exploracin de este
recurso. Estos yacimientos de gas no convencional presentan inconvenientes ya
que estas formaciones de shale tienen una mineraloga muy variable que puede
estar dominada por minerales frgiles, como el cuarzo, los carbonatos y
feldespatos. La porosidad de la matriz de shale es muy variable. La permeabilidad
del shale es una funcin del tamao de grano fino y la porosidad (grado de
compactacin) y es extremadamente baja (por debajo de la gama de micro y
nanodarcy), lo que hace necesario caracterizarlos.

La necesidad de evaluar los parmetros fsicos de yacimientos de gas shale en
Colombia, se debe a su potencial que ha sido estimado aproximadamente en 82
Tcf de volmenes recuperables en cuencas que presentan prospectividad tales
como el Valle Medio del Magdalena, Cordillera Oriental, Cesar-Ranchera.
Tomando como base para este estudio, datos provenientes de la 7 Semana
Tcnica Internacional de la Universidad Surcolombiana y de documentos
realizados tanto a nivel nacional (UPME) como internacional (DLittle).

Se procede a determinar el potencial de reservas del pozo USCO-69 mediante la
siguiente tabla y con todas las referencias del captulo 2 y 3, con la suposicin que
el rea de drene sea de 120 acres.

Datos
porosidad total 0.21
saturacin agua 0.35
densidad roca - g/cm
3
2.26
factor vol. Gas - bbs/scf 0.0042
area acres 120
espesor - ft 2500
Metano
Peso Molecular 16
PL, psia 977.29
VL, scf/ton 15.88
densidad g/cm3 0.372
Tabla 4.1 Datos del Pozo USCO-69

- Resultados obtenidos de Gas almacenado, donde se puede apreciar la
sobreestimacin de las reservas al pasar por alto una reduccin en la
porosidad.
56



Fig. 4.1 Gas almacenado estimado por el
Mtodo convencional y el de Ambrose


- Resultados del perfil de Adsorcin


Fig. 4.2 Contenido de Gas y Presin de Desorcin del pozo USCO-69

A continuacin se muestran los resultados obtenidos al comparar 5 de las
formaciones ms importantes de EEUU y el pozo USCO-69. Usando las
referencias del captulo 3 se tiene que:



458
460
462
464
466
468
470
472
474
476
0 2000 4000 6000 8000
C
o
n
t
e
n
i
d
o

d
e

G
a
s

s
c
f
/
t
o
n

Presin, psia
Gas total corregido
Gast total - Industria
0
2
4
6
8
10
12
14
16
0 2000 4000 6000 8000
C
o
n
t
e
n
i
d
o

d
e

G
a
s

s
c
f
/
t
o
n

Presin, psia
Gas adsorbido
Presin de desorcin
57


DATOS
USCO-69 BARNETT FAYETTEVILLE
v 0.242 v 0.21 v 0.21
E 3.6E+05 E 5.2E+06 E 7.3E+06

MARCELLUS WOODFORD

v 0.24 v 0.2
E 3.7E+06 E 3.90E+06
Tabla 4.2 Mdulo de Young y Razn de Poisson

Los resultados obtenidos se produjeron con la tabla 4.2 y algunas suposiciones:
Los datos de la fraccin de la roca f, la compresibilidad del grano , y la relacin
del esfuerzo de langmuir con la constante de langmuir l/B, se tomaron como
constante para todos los ejemplos y son respectivamente 1, 1.30e-06 psi
-1
, y 8.0;
Adems, se utiliz la ecuacin de Palmer&Mansoori propuesta en el captulo 3.


Fig. 4.3 Presin vs Razn de Permeabilidad

Como se puede observar el mdulo de Young del pozo USCO-69 es muy bajo
comparado con el de los yacimientos en EEUU, lo que significa que a un esfuerzo
aplicado a la superficie de la roca va a haber una deformacin mayor; por lo que
se deduce que la roca colombiana cuenta con un alto contenido de arcillas y
material dctil (mayor que el de las formaciones americanas), corroborado por un
anlisis DRX realizado a unas muestras de la formacin del valle inferior a la cual
pertenece el pozo. Sin embargo tambin se ve un comportamiento anmalo en la
relacin k/ki y se cree que es mayormente debido a que las muestras no fueron
probadas en el laboratorio, mediante pruebas de desorcin y de elasticidad, para
determinar el valor verdadero de l/B, que se toma en ste ejemplo igual a 8.
Tambin, es obvio que, la suposicin de la compresibilidad del grano es algo
0.94
0.96
0.98
1
1.02
1.04
1.06
1.08
1.1
0 500 1000 1500 2000 2500 3000 3500 4000

K
/
K
i

Presin, psia
USCO-69
BARNETT
FAYETTEVILLE
MARCELLUS
WOODFORD
58

drstica debido a la variabilidad en la mineraloga mostrada en la tabla 4.2 que
deja entre dicha una gran posibilidad composicional para las rocas encontradas a
lo largo del Valle Inferior del Magdalena.

PROPIEDADES WOODFORD BARNETT USCO-69
% 4.7 7.17 16-21
Arcilla Total % 34 27 44-54
Cuarzo Total % 33 33 37-40
Carbonato Total % 8 20 1-50
TOC % 7.52 4.8 0.7-1.2
INDX DE FRACT. W&G Fraccin 0.443 0.509 0.420
v 0.2 0.21 0.242
E 3.90E+06 5.18E+06 3.60E+05
YM_n
41.450 59.684 -9.147
PR_n 80 76 63.2
INDX DE FRACT. RICHMAN % 60.725 67.842 27.0264
Tabla 4.3 Propiedades promedio de las formaciones ms relevantes

Se ve que se utilizaron las 2 ecuaciones propuestas para determinar el ndice de
fractabilidad de la roca. Richman utiliz los parmetros geomecnicos de la roca
para determinar un valor aproximado de las fuerzas de deformacin que actuaran
sobre varias rocas al aplicarles un mismo esfuerzo. Por su parte Wang&Gale,
decidieron hacerlo mediante parmetros mineralgicos. Como se ve en la tabla 4.3
los valores mineralgicos del pozo USCO-69. Se determin en ste estudio, que
debido a una falta de mineraloga por parte de las muestras y a lo cerca que se
encuentra ste valor del pozo de Woodford con USCO-69 que era ms confiable
usar la ecuacin de Richman normalizada, debido a que con menos parmetros
(2) se logr una mayor precisin (vase la Figura 4.3) en la determinacin del
ndice de fractabilidad.

Con diferentes datos de Mdulo de Young, tomados de distintos pozos de la
misma cuenca, se encuentra una variedad de relaciones, como se propuso
anteriormente (Siendo el caso 2 el mostrado anteriormente). Ver Figura 4.4

Siendo consecuente con el estudio, se grafic la Razn de Poisson contra el
Mdulo de Young para comprobar que la figura 2.3 de Richman se acercaba de
alguna manera a las conclusiones anteriormente mostradas en ste captulo. Ver
Figura 4.5


59


Fig. 4.4 P vs k/ki proporcionada por diferentes pozos de una misma cuenca



Fig. 4.5 Razn de Poisson vs Mdulo de Young.

La figura 4.5 demuestra que el grfico propuesto por Richman, donde dice que las
rocas que se encuentren en el 3er cuadrante (esquina izquierda inferior) tienen
afinidad a la fractura y que las del 1er cuadrante (esquina derecha superior) tienen
afinidad a absorber el impacto de los esfuerzos tiene relacin con la figura 4.3.



0.88
0.9
0.92
0.94
0.96
0.98
1
1.02
0 1000 2000 3000 4000
R
e
l
a
c
i

n

k
/
k
i

Presin, psia
1
2
3
0.0E+00
1.0E+06
2.0E+06
3.0E+06
4.0E+06
5.0E+06
6.0E+06
7.0E+06
8.0E+06
9.0E+06
1.0E+07
0.1 0.15 0.2 0.25 0.3 0.35 0.4
M

d
u
l
o

d
e

Y
o
u
n
g

Razn de Poisson
USCO-69
BARNETT
FAYETTEVILLE
MARCELLUS
WOORDFORD
60


4.1 Prediccin y evaluacin del comportamiento del transiente de presin en
el yacimiento
Para realizar esta evaluacin se tomaron los datos de la historia de produccin de
un pozo del Shale Marcellus en Estados Unidos, a partir de un pronstico a 5 aos
(60 meses) que tiene como propiedades del yacimiento y del pozo las siguientes:


Lex = 4175 ft Ye = 1120 ft
Nmero de Fracturas = 13 Xf = 150 ft
Pi = 3705 psia FCD = 20
T = 115 F K1 = 0.0007 mD
h = 150 ft K2 = 0.0007 mD
C = 7.1 Sgi = 66.2%
Cf = 5.6049E-06 1/psi Swi = 33.8%
VL = 85 scf/ton b = 2.6 g/cc
PL = 468 psia Rw = 0.35 ft
Ct = 1.3044E-04 1/psi
Tabla 4.4 Datos iniciales para la evaluacin del transiente de presin en el
yacimiento campo Marcellus.

Figura 4.6 Comportamiento de la presin y la produccin del pozo Marcellus.
Podemos observar que el pozo marcellus produjo a una rata inicial mxima de
5.31 MMpcsd de gas durante los 5 aos (60 meses).

Se procede a realizar el ajuste histrico del pozo Marcellus, una vez realizado el
ajuste se tiene los siguientes parmetros con las estimaciones de OGIP,
OGIPSVR y el gas recuperable para el pozo Marcellus, mostrados en la tabla 4.5
y se realiza el comportamiento de la presin y la produccin del pozo Marcellus
ajustado (ver Figura ).
61


Lex = 4175 ft Ye = 1120 ft
Nmero de Fracturas = 13 Xf = 45 ft
Pi = 3705 psia FCD = 32.061
T = 115 F K1 = 0.0037 mD
h = 150 ft K2 = 0.0011 mD
C = 7.1 Sgi = 66.2%
Cf = 5.6049E-06 1/psi Swi = 33.8%
VL = 85 scf/ton b = 2.6 g/cc
PL = 468 psia Rw = 0.35 ft
Ct = 1.3044E-04 1/psi
Resultados de Produccion
OGIPF = 8821 MMscf OGIPA = 4296 MMscf
OGIP = 13117 MMscf OGIPSRV = 1056 MMscf
EURG = 2373 MMscf RRG = 1774MMscf
Tabla 4.5 Resultados de la estimacione de OGIP, OGIPSVR y el gas recuperable
para el pozo Marcellus




Figura 4.7 Comportamiento ajustado de la presin y la produccin del pozo
Marcellus.

De acuerdo a la Tabla 4.5, se observa que la mayora de las propiedades se
mantienen constantes, slo vara la permeabilidad de rea estimulada, la
permeabilidad del rea no estimulada (matriz), la longitud de la fractura y la
conductividad adimensional.
Inicialmente se toman ambas permeabilidades como 0.0007 md porque la
informacin suministrada para ste caso slo habla de una permeabilidad, por lo
que se entiende que sea una permeabilidad promedio. Luego a la hora de simular
lo que se logr fue encontrar de manera coherente una k1> k2 como resultado del
fracturamiento hidrulico.
62

El valor de la longitud de la fractura cambi con respecto al real, lo que indica que
el simulador encontr este valor para realizar el ajuste, ya que, de haber utilizado
250 pies no se hubiera logrado. Vale la pena aclarar que el cambio de algunas
propiedades a la hora de hacer el ajuste histrico, es slo para hacer predicciones
o pronsticos futuros.
Por ltimo se nota que el valor de FCD se conserva como una conductividad finita.

El modelamiento se realiz utilizando el simulador FEKETE Harmony 2012 v1, con
el cual se construy el esquema del pozo con mltiples fracturas y el esquema de
las respuestas de la cada de presin y de la derivada del pozo caso base que se
muestra en la figura 4.8 y 4.9.


Figura 4.8 Respuestas de la cada de presin y de la derivada del pozo Marcellus.




Figura 4.9 Esquema del pozo multifracturado para el caso base realizado con
el simulador.

63

Despues de a ver realizado el esquema de las respuestas de la cada de presin
y de la derivada del pozo Marcellus (figura 4.8), se procede con la ayuda de
F.A.S.T welltest a identificar los regmenes de flujo tpicos para pozos
horizontales de gas shale con fracturas hidrulicas del pozo Marcellus (ver figura
4.10).

Figura 4.10 Derivada del pozo y los regmenes de flujo tpicos para pozos
horizontales de gas shale con fracturas hidrulicas, campo Marcellus.

Ahora con los valores de VL=85 y PL=468 se construy la isoterma de adsorcin
para este caso base (ver figura 4.11).

64


Figura 4.11 Isoterma de adsorcin para el pozo Marcellus.
Teniendo en cuenta los datos de produccin, del yacimiento y el esquema de pozo
y consideraron efectos de adsorcin incluyendo las constantes de Langmuir y
efectos geomecnicos se realizo el comportamiento del transiente curvas tipo
radiales con un Gp normalizada, como lo indica los datos de la tabla 4.6 para
tiempos de radial y de balance de materia y esquema del transiente se muestra en
las figura 4.12 y 4.13.

rED 7

rED 12
Rata Normalizada Rata Normalizada
Adsorcin Adsorcin
GeoMecnica GeoMecnica
TIME MBT TIME Radial
rED 7 rED 12
s -5.566 s -5.566
A 30 acres A 87 acres
Gp 599 MMscf Gp 599 MMscf
Pab 319.15 psi (a) Pab 319.15 psi
RFg 80% RFg 80%
Pavg 2713.17 psi (a) Pavg 3331.31 psi (a)
Tabla 4.6 Resultados del anlisis del transiente curvas tipo radial para el pozo
Marcellus.

65

Figura 4.12 Comportamiento del transiente para tiempos de balance de materia
curvas tipo radiales para el pozo Marcellus.
Figura 4.13 Comportamiento del transiente para tiempos de pseudo-superposicin
radial curvas tipo radiales para el pozo Marcellus.
Tambin se realizo el comportamiento del transiente curvas tipo para Fractura de
conductividad finita que es la que maneja el pozo Marcellus, con un Gp
normalizada consideraron efectos de adsorcin incluyendo las constantes de
Langmuir y efectos geomecnicos, como lo indica los datos de la tabla 4.7 para
tiempos de radial y de balance de materia y esquema del transiente se muestra en
las figura 4.14 y 4.15.


66

FCD 500

FCD 500
xED 5 xED 25
Rata Normalizada Rata Normalizada
Adsorcin Adsorcin
GeoMecnico GeoMecnico
TIME MBT TIME Radial
FCD 500 FCD 500
xED 5 xED 25
KfW 1994 md.ft KfW 1968 md.ft
A 65 acres A 827 acres
Gp 599 MMscf Gp 599 MMscf
Pab 800.7 psi (a) Pab 800.7 psi (a)
RFg 80% RFg 80%
Pavg 3181.14 psi (a) Pavg 3660.51 psi (a)
Tabla 4.7 Resultados del anlisis del transiente curvas tipo para fractura de
conductividad finita para el pozo Marcellus.

Figura 4.14 Comportamiento del transiente para tiempos de balance de materia
curvas tipo para fractura de conductividad finita para el pozo Marcellus.


67

Figura 4.15 Comportamiento del transiente para tiempos de pseudo-superposicin
radial curvas tipo para fractura de conductividad finita para el pozo Marcellus.

Por ltimo y tambin con la ayuda del simulador FEKETE Harmony 2012 v1 se
realizo con base al modelo de Brooks Corey que es ampliamente aceptado, el
comportamiento y clculo de la presin capilar con un

=20% y un

=
20% , mostrado en las figuras 4.16 y 4.17.

Figura 4.16 Curva de presin capilar agua-crudo y agua-gas.

68



Figura 4.17 Curva de presin capilar gas-crudo.

4.1.1 Prediccin y evaluacin de un pozo colombiano USCO 69

Para ste caso Colombiano se tomaron los datos de la historia de produccin de
un pozo USCO 69 en Colombia, a partir de un pronstico de un mes que tiene
como propiedades del yacimiento y del pozo las siguientes:

Pi = 3745 psia C = 21
K1 = 0.0010 mD h = 2500 ft
VL = 15.88 scf/ton Swi = 35%
PL = 977.29 psia Sgi = 65%
Cf = 3.5737E-06 1/psi b = 2.26 g/cc
Ct = 1.2368E-04 1/psi Nmero de Fracturas = 6
Contenido organico total (TOC)= 13%
RO = 0.63%
Predonima Kerogeno = Tipo III Y II
Tabla 4.8 Datos iniciales para la evaluacin del yacimiento y del pozo USCO 69.


Una vez realizado el ajuste se tiene los siguientes parmetros con las
estimaciones de OGIP, OGIPSVR y el gas recuperable para el pozo USCO 69,
mostrados en la tabla 4.9


69

Lex = 4175 ft Ye = 800 ft
Nmero de Fracturas = 6 Xf = 250 ft
Pi = 3546.74 psia FCD = 10
T = 120 F K1 = 0.0020 mD
h = 2500 ft K2 = 0.0006 mD
C = 21 Sgi = 65%
Cf = 3.5737E-06 1/psi Swi = 35%
VL = 15.88 scf/ton b = 2.26 g/cc
PL = 977.29 psia Rw = 0.35 ft
Ct = 1.2368E-04 1/psi
Resultados de Produccion
OGIPF = 291349 MMscf OGIPA = 8437 MMscf
OGIP = 299786 MMscf OGIPSRV = 187366 MMscf
EURG = 8354 MMscf RRG = 7755 MMscf
Tabla 4.9 Resultados de la estimacin de OGIP, OGIPSVR y el gas recuperable del
pozo USCO 69.


Figura 4.18 Comportamiento de produccin y presin del pozo USCO 69

Se realiza el comportamiento de la presin y la produccin del pozo USCO 69
ajustado (ver Figura 4.18), donde, podemos tambin observar que el pozo USCO
69 produjo a una rata inicial mxima de 1.4 MMpcsd de gas durante un meses.

70

De acuerdo a la Tabla 4.9, se observa que la mayora de las propiedades se
mantienen constantes, slo vara la permeabilidad de rea estimulada, la
permeabilidad del rea no estimulada (matriz), la longitud de la fractura, la
conductividad adimensional y pronosticando un valor del ancho del yacimiento de
800ft , segn el ajuste histrico del simulador al yacimiento Colombiano.
,
El simulador encontr el valor de la longitud de la fractura de 250 ft para realizar el
ajuste, vale la pena aclarar nuevamente que el cambio de algunas propiedades a
la hora de hacer el ajuste histrico, es slo para hacer predicciones o pronsticos
futuros.

Inicialmente teniendo la permeabilidad del rea estimulada igual a 0.0010 md,
informacin suministrada de tablas de datos del pozo USCO 69, se simulo
logrando encontrar de manera coherente una k1> k2 como resultado del
fracturamiento hidrulico.Por ltimo se nota que el valor de FCD se conserva
como una conductividad finita.


El modelamiento se realiz utilizando el simulador FEKETE Harmony 2012 v1, con
el cual se construy el esquema del pozo con mltiples fracturas y el esquema de
las respuestas de la cada de presin y de la derivada del pozo caso base que se
muestra en la figura 4.19 y 4.20.


Figura 4.19 Esquema del pozo multifracturado para el pozo USCO 69 realizado con
el simulador.

71


Figura 4.20 Respuestas de la cada de presin y de la derivada del pozo USCO 69.

Ahora con los valores de VL=15.88 y PL=977.29 se construy la isoterma de
adsorcin para el pozo USCO 69 (ver figura 4.21) y el comportamiento de la
presin en el SRV, como lo indica la figura 4.22.


Figura 4.21 Isoterma de adsorcin para el pozo USCO 69.

72


Figura 4.22 Comportamiento de la presin en el SRV para el pozo USCO 69.


Despus de a ver realizado el esquema de las respuestas de la cada de presin
y de la derivada del pozo USCO 69 (figura 4.20), fue imposible identificar los
regmenes de flujo tpicos para pozos horizontales de gas shale con fracturas
hidrulicas con la ayuda de F.A.S.T welltest ya que la prueba se hizo a un
intervalo muy corto de tiempo.

FCD 500
xED 5
Rata Normalizada
Adsorcin
GeoMecnico
TIME MBT
FCD 500
xED 5
KfW 0.1 md.ft
A 0.003 acres
Gp 1 MMscf
Pab 763.15 psi (a)
RFg 80%
Pavg 3472.25 psi (a)
Tabla 4.10 Resultados del anlisis del transiente curvas tipo para fractura de
conductividad finita para el pozo USCO 69.

Tambin se realizo el comportamiento del transiente curvas tipo para Fractura de
conductividad finita que es la que maneja el pozo Marcellus, con un Gp
normalizada consideraron efectos de adsorcin incluyendo las constantes de
Langmuir y efectos geomecnicos, como lo indica los datos de la tabla 4.10 para
Pseudo tiempo de balance de materia y esquema del transiente se muestra en las
figura 4.23

73


Figura 4.23 Comportamiento del transiente para tiempos de balance de materia
curvas tipo para fractura de conductividad finita para el pozo USCO69.



4.2 Arreglo o disposicin de pozos


Apartir de los datos de pozo y del yacimiento correspondientes a las condiciones
iniciales del pozo Marcelllus y del pozo USCO 69, con el fin de determinar la mejor
estrategia de perforar nuevos pozos en un yacimiento se utilizara el software
especializado F.A.S.T. Evolution se realizando 2 arreglos de pozos horizontales
multifracturado anteriormente nombrados. A partir de los cuales, se obtuvieron
datos de presin promedio, presin de pozo fluyendo y produccin acumulada por
pozo para cada arreglo, adems se evalu el mejor escenario econmico.

En cada arreglo se llev a cabo un pronstico a 5 aos (60 meses) y presiones de
fondo en cada pozo de 1000 psia. La longitud lateral de todos los pozos es de
4175 pies, slo se cambi el nmero de fracturas en cada arreglo.

El valor de G estimado y consultado del caso base y caso Colombiano.
Los parmetros o variables econmicas involucradas en el anlisis, para todos los
arreglos son mostrados a continuacin:

- Costo del pozo (US$/Pozo): Costo capital por pozo (perforacin,
completamiento, estimulacin, etc.), ste valor cambia para cada arreglo
debido a la diferencia en el nmero de fracturas.
- Costo de Compresin (3500 US$/hp): Costo de adquisicin de un
compresor. Se calcula con base a la presin de entrada y la presin de
descarga.
74

- Costos Operativos fijos (2500 US$/Pozo)/mes: Costos de labores
operativas, alquiler, impuestos de propiedad, etc.
- Costos Operativos variables (0.25 US$/Mpcs): Costos variables de
procesamiento, recoleccin y transporte de gas.
- Regalas (23%)
- Eficiencia de la Compresin (85%): Eficiencia del compresor
- Prdidas en superficie (5%): Porcentaje de prdidas de volumen de gas en
superficie.
- Precio del Gas (4 US$/Mpcs)
- Tasa de Descuento (10%): Porcentaje de descuento por ao Los valores de
cada parmetro fueron tomados como valores promedio de las principales
cuencas productoras de Estados Unidos.

Figura 4.24 Valor presente neto del proyecto para pozo Marcellus.

Figura 4.25 Valor presente neto del proyecto para pozo USCO 69

75

CAPITULO V - CONCLUSIONES

- La fractabilidad y la mineraloga juegan un papel fundamental para
determinar la capacidad que tiene el play de ser productivo. Puesto que al
ser la roca que lo conforma, ms plstica (contenido de arcilla alto)
convierte el proceso de fracturamiento en inviable, se necesitar mucha
ms energa para poder fracturar la roca y sta no va a responder
favorablemente debido a su plasticidad (capacidad de adsorcin y
dispersin de energa).

- La fisilidad y el cambio en la disposicin de la matriz, por la sorcin del gas
y por las caractersticas geomecnicas que sta posee, pasan a ser
parmetros claves al evaluar el potencial productivo del play de gas shale,
porque pueden restringir o aumentar significativamente el flujo de gas y el
factor de recobro del yacimiento.

- El Indice de fractabilidad de Richman y la ecuacin de Wang y Gale se
correlacionaron, sin embargo (ver la tabla 4.3) los valores del pozo
Woodford y USCO-69 son muy cercanos al determinarlos por el 2 mtodo,
y cuentan con una mineraloga muy diferente. Por lo que se determin que
el mejor mtodo para evaluar ste ndice era la ecuacin de Richman.

- Dos (2) modelos de encogimiento de la matriz fueron mostrados. El primer
modelo (Seiddle y Huitt) sustituy con una constante una relacin de los
mdulos axial y bruto en el trmino de encogimiento en el segundo modelo,
por sta razn se determin que el mejor modelo para evaluar el cambio de
la porosidad y de la permeabilidad era el de Palmer y Mansoori.

- El pozo USCO-69 cuenta con el repunte de permeabilidad, pero ste ocurre
muy tarde debido al alto contenido de arcillas y a su baja fisilidad.

- Es necesario corregir el volumen de gas almacenado, el cual puede
alcanzar un cambio de al menos un 0.3% en el pozo USCO-69. Este
parmetro es clave para determinar correctamente el gas in situ del
yacimiento.

- Se evalu estrategia tcnico-econmica, para realizar el anlisis de
rentabilidad de las fracturas, que tienen cada uno de los pozos horizontales
perforados en estudio en un yacimiento de una extensin dada, donde se
demuestra que el pozo USCO 69 necesita de una menor inversin es la del
arreglo generando un mayor valor presente neto.


76

- La simulacin numrica permiti en un determinado tiempo desarrollar un
ajuste histrico ms preciso y a partir de ste se llevaron a cabo
pronsticos en la produccin y presin de cada uno de los casos.

- Los resultados obtenidos de las simulaciones del pozo Marcellus y del pozo
USCO 69 demuestra que en Colombia existe un gran potencial y un futuro
comprometedor para la explotacin de gas shale.








































77





REFERENCIAS

1. WVExp.com. First Natural Gas Well [En Lnea]
http://www.wvexp.com/index.php/First_Natural_Gas_Well [Citado en Agosto, 2012]

2. Cornell University Cooperative Extension. Shale gas: a short history from NETL
[En Lnea]
http://www.cce.cornell.edu/EnergyClimateChange/NaturalGasDev/Documents/PDF
s/Shale%20Gas-%20a%20short%20history%20from%20NETL.pdf [Citado en
Agosto, 2012]

3. Diego Cabot - La Nacin (Argentina), REPORTS OF GAS SHALE IN
ARGENTINA. 15 de Mayo de 2011

4. Pontificia Univ. Catlica de Chile, LA REVOLUCIN DEL SHALE GAS. Mayo de
2011

5. Hernn F. Pacheco, EnerDossier, ANLISIS SOBRE EL MERCADO
ENERGTICO MUNDIAL
19 de Agosto de 2011

6. Unidad de Pleaneacin Minero Energtica y Ministerio de Minas y Energa.
PLAN DE ABASTECIMIENTO PARA EL SUMINISTRO Y TRANSPORTE DE
GAS NATURAL Noviembre 2010 [Citado en Agosto, 2012]

7. Julio Csar Vera Daz, PERSPECTIVAS DE LOS HIDROCARBUROS NO
CONVENCIONALES EN COLOMBIA XIII Congreso Nacional y IV Internacional de
Servicios Pblicos y TIC en Medelln, 22 de Junio de 2011, [Citado en Agosto,
2012]

8. Jesus David Pieros y Maria Adelaida Velsquez. InterBolsa, Pacific Rubiales
Energy: Predescubrimiento Cotorra [En Lnea]
http://www.interbolsa.com/c/document_library/get_file?uuid=db969bbb-f65a-4523-
93fb-e93198686a5b&groupId=81085

9. Metodologa para el modelamiento de produccin post fractura en proyectos de
shale gas. J.C, Crdenas Natfrac-ICP, R.D, Castillo Ecopetrol, E.R, Prez
Ecopetrol. Presentado en el XIV Congreso Colombiano del Petroleo organizado
por ACIPET. Bogota 22-25 de Nov 2001.

10. Shale Gas & Shale Oil: Solutions for your Shale Gas and Liquid-Rich
Challenges. Halliburton. 2011
78


11. Screening Criteria for Shale-Gas Systems. Fred P. Wang and Julia F. W. Gale.
Bureau of Economic Geology, Jackson School of Geosciences, The University of
Texas at Austin, University Station, Box X, Austin, Texas 78713

12. A practical use of Shale Petrophysics for stimulation design optimization: All
shale plays are not clones of the Barnett Shale. SPE 115258 Rick Rickman, Mike
Mullen, Erik Petre, Bill Grieser and Donald Kundert, SPE, Halliburton 2008.

13. Multi-component Sobed-phase Considerations for Shale Gas-in-place
Calculations. Ray J. Ambrose, The University of Oklahoma, Robert C. Hartman,
Weatherford Labs, and I. Yucel Akkutly, the University of Oklahoma.

14. E. Shtepani, L.A. Noll, L.W. Elrod and P.M. Jacobs A New Regression Based
Method for Accurate Measurement of Coal and Shale Gas Content 2010

15. NAMRITA SONDHI, Norman, Oklahoma. PETROPHYSICAL
CHARACTERIZATION OF EAGLE FORD SHALE 2011

16. http://www.agenciadenoticias.unal.edu.co/nc/detalle/article/crean-medidor-de-
gas-metano-para-minas-de-
carbon.html?cHash=ce5b363ef06f67e4183cdeda2592affc

17. Ray J. Ambrose, The University of Oklahoma, Robert C. Hartman, Weatherford
Labs, and I. Yucel Akkutlu, The University of Oklahoma, EEUU SPE 141416
MULTI-COMPONENT SORBED-PHASE CONSIDERATIONS FOR SHALE GAS-
IN-PLACE CALCULATIONS 2729 de marzo del 2011

18. PRESSURE TRANSIENT CHARACTERISTICS OF HYDRAULICALLY
FRACTURED HORIZONTAL SHALE GAS WELLS. Yueming Cheng, SPE, West
Virginia University. 17 19 de agosto del 2011

19. Bates, R.L., y Jackson, J.A. (eds.). 1987. Glosario de geologa, 3a. ed.
American Geological Institute, Alejandra, VA. 788 p.

20. Impact of Water Dynamics in Fractures on the Performance of Hydraulically
Fractured Wells in Gas Shale Reservoirs SPE 127863

21. Application of Thermogravimetric Fourier Transform Infrared Spectroscopy
Analysis in Determining Heating Protocol for Porosity Measurements in Shales,
Leonel Prez Solano, Neiva, Febrero 2011

22. Calculos de Shale gas in situ Parte II Efectos de la Adsorcin
Multicomponente del Gas. Ambrose. 2010

79

23. Experimental Measurement of Coal Matrix Shrinkage Due to Gas Desorption
and Implications for Cleat Permeability Increases John P. Seidle and L. G. Huitt,
Amoco E&P Technology Group. 1995

24. How Permeability Depends on Stress and Pore Pressure in Coalbeds: A New
Model Ian Palmer, SPE, and John Mansoori, Amoco Production Co. 1996 SPE
Annual Technical Conference & Exhibition, Denver, 69 October

25. Fekete Harmony/Reference Materials/General Concepts/Rate Transient
Analysis

26. Parmetros de Rocas Generadoras, Universidad Surcolombiana, Jorge Arturo
Camargo Puerto.