Análisis de producción en gas no convencional
Análisis de producción en gas no convencional
FACULTAD DE INGENIERÍA
TESIS
Que para
00 obtener el título de
Ingeniero Petrolero
P R E S E N T A (N)
Abimael Avila Torres
DIRECTOR(A) DE TESIS
Dr. Jorge Alberto Arévalo Villagrán
Ilustración 2-15 Rango de aplicación del método de presión al cuadrado a 300°F ... 24
Ilustración 2-17 Grafica de 1/p contra pseudo presión con viscosidad del gas ideal y
gas real............................................................................................................................................ 26
Ilustración 2-19 𝒑𝑫𝒕𝑫 𝒚 𝒎𝑫𝒕𝑫 𝒗𝒔 𝒕𝑫 para flujo radial para líquido, gas ideal y gases
reales. .............................................................................................................................................. 28
Ilustración 2-20 Relación Presión-tiempo de los distintos flujos de fluidos................ 30
Ilustración 3-2 Grafica de diagnóstico con sus distintos regímenes de flujo. ............ 39
Ilustración 3-9 Grafica especializada de flujo para régimen lineal. donde podemos
apreciar ambos casos: con intercepción en cero y con intercepción en la recta. ..... 46
Ilustración 3-14 Grafica especializada de flujo para régimen de flujo dominado por la
frontera externa a la condición de producción de presión de fondo fluyendo
constante,𝒑𝒘𝒇. .............................................................................................................................. 54
Ilustración 4-1 Grafica 𝒒𝑮𝒑𝒗𝒔 𝒕 para determinar los parámetros "a" y "m". ................ 62
Ilustración 4-2 Cuatro pasos propuestos para calcular el ajuste de producción. ...... 63
Ilustración 4-4 Caso de simulador(k=0.1 mD, A=200 acres, Fcd=200,xf=300 ft, h=40)
se muestra un EUR 8% menor , pero un desajuste en el perfil de producción, si
usamos el potencial de recuperación..................................................................................... 66
Ilustración 4-5 El mismo caso de simulación que la Fig. 4-4, el cual muestra un
excelente ajuste en ambos perfiles de producción, si únicamente usamos 𝟐𝒏𝒅 y 𝟑𝒏𝒅
año de los datos de producción............................................................................................... 67
Ilustración 4-6 Grafica Especializada YM-SEPD .................................................................. 68
Ilustración 4-7 Esquema conceptual de bloques de matriz del modelo lineal para
pozos hidráulicamente fracturados. Área transversal de la cara del pozo,𝑨𝒄𝒘 =
𝟐𝒙𝒆𝒉. ................................................................................................................................................ 72
Ilustración 4-8 Grafica de las cinco regiones de flujo del moldeo de Bello. ................ 73
Ilustración 5-13Grafica especializada del método YM-SEPD para el pozo CIGE 231. 86
Ilustración 5-16 Grafica especializada del método de Potencias para el pozo CIGE
231. ................................................................................................................................................... 88
Ilustración 5-20 Histórico de presión- producción del pozo CIGE 218. ......................... 92
Ilustración 5-54 Histórico de producción de gas del pozo Shale C. ............................. 119
Ilustración 5-63 Grafica especializada de gasto de gas contra t(a.m.),𝒒𝒈 𝒗𝒔 𝒕𝒂, 𝒎, del
pozo Shale C................................................................................................................................ 125
Ilustración 5-68 Grafica especializada del método de YM-SEPD del pozo Shale C. 129
Ilustración 5-75 Comparación de ajuste de la producción de gas del pozo CIGE 231,
mediante curvas de declinación ............................................................................................ 140
Ilustración 5-76 Grafica de comparación del ajuste de la producción de gas del pozo
CIGE 218. ...................................................................................................................................... 142
Ilustración 5-77 Grafica de ajuste de la producción de gas del pozo Shale A. ......... 144
Ilustración 0-1 Propiedades pseudo críticas de los gases naturales. ......................... 180
Ilustración 0-2 Corrección de las propiedades pseudo- criticas 𝑯𝟐𝑺 𝒚 𝑪𝑶𝟐. ............. 183
Ilustración 0-3 Factor z para gases naturales con 𝟎 < 𝒑𝒓 < 𝟏𝟎. .................................... 186
Ilustración 0-4 Factor z para gases naturales con 𝟗 < 𝒑𝒓 < 𝟐𝟎. .................................... 187
Ilustración 0-5 Variación de la 𝒄𝒓𝑻𝒓 para gases naturales con 𝟏. 𝟓 < 𝑻𝒓 ≤ 𝟏. 𝟒(esbozo
izquierdo) y variación en 𝒄𝒓𝑻𝒓 para gases naturales con 𝟏. 𝟒 < 𝑻𝒓 ≤ 𝟑(esbozo
derecho). ....................................................................................................................................... 190
Ilustración 0-2 Flujo lineal ideal dentro de las fracturas. ................................................ 213
Ilustración 0-6 Grafica normal del modelo tipo curva para régimen de flujo lineal con
intersección en cero y no intersección en cero de la línea recta, bajo la condición de
producción a gasto constante de gas,𝒒𝒈. ........................................................................... 218
Ilustración 0-7 Grafica normal del modelo de tipo curva para régimen de flujo lineal
con intersección en cero y no intersección en cero de la línea recta, bajo la
condición de presión de fondo fluyendo constante,𝒑𝒘𝒇. ............................................... 220
Ilustración 0-12 Grafica semi- log del modelo tipo curva para régimen radial con no
intersección en cero de la línea recta, bajo la condición de producción a gasto de
flujo constante de gas,𝒒𝒈. ....................................................................................................... 229
Ilustración 0-13 Grafica semi- log del modelo tipo curva para régimen radial con no
intersección en cero de la línea recta, bajo la condición de producción a presión de
fondo fluyendo constante,𝒑𝒘𝒇. .............................................................................................. 232
Ilustración 0-14 Ilustración de grafica normal del modelo de tipo curva para régimen
de flujo esférico para evaluar algunos valores de yacimientos de baja permeabilidad
con no intercepción en cero de la línea recta, bajo la condición de producción a
gasto de flujo constante de gas,𝒒𝒈. ...................................................................................... 235
Ilustración 0-17 Grafica normal del modelo de tipo curva para régimen de flujo
dominado por la frontera externa con no intersección en cero de la línea recta, bajo
la condición de producción a gasto de flujo constante,𝒒𝒈. ........................................... 239
Ilustración 0-18 Grafica de modelo semi- log del modelo tipo curva para régimen de
flujo dominado por la frontera externa con no intersección en cero de la línea recta,
bajo la condición de producción a presión de fondo fluyendo constante,𝒑𝒘𝒇. ....... 244
LISTA DE TABLAS
Tabla 3-4 Funciones específicas del tiempo para las distintas geometrías de flujo
.............................................................................................................................. 45
Tabla 5-9 Información General de los casos de capo de pozos de gas de baja
permeabilidad. ..................................................................................................... 135
Tabla 5-10 Información de los yacimientos de gas de baja permeabilidad ......... 135
Tabla 5-11 Datos de los fluidos obtenidos de los pozos productores de gas en
formaciones de baja permeabilidad. .................................................................... 136
Tabla 5-12 resultados calculados de los regímenes de flujo bilineal y lineal tardío.
............................................................................................................................ 136
Tabla 5-13 Resultados calculados del regimen de flujo dominado por la frontera
externa. ............................................................................................................... 137
Tabla 5-16 Resultados del ajuste de producción de gas de los modelos YM-
SEPD y Potencias. .............................................................................................. 138
Tabla 5-17 Resultados del ajuste de la producción de gas del modelo de
Crecimiento Logístico. ......................................................................................... 139
Tabla 5-18 Comparación del pronóstico de EUR- CIGE 231 .............................. 140
Tabla 11-1 Tabla de propiedades físicas de los gases a 14.7 psia y 60 °F. ........ 175
1. CAPITULO I INTRODUCCION
1.1. Resumen
1
datos de historia de producción, presión, porosidad, propiedades de los fluidos,
gasto de producción, espesor, presión inicial y temperatura. El análisis de los pozos
de formación de lutitas es más complejo, por ello se adiciono el análisis de
metodología de Bello (Bello, 2009), el cual presenta una estimación de la
permeabilidad y volumen original de gas,𝐺, mediante el uso de función
pseudotiempo y gasto normalizado utilizando las gráficas especializadas para
diferentes regímenes de flujo.
Mediante el análisis de los datos de producción a través de los diferentes modelos;
se presentan los resultados de los regímenes identificados en cada pozo, las
propiedades del yacimiento calculadas mediante gráficas especializadas y el
volumen de gas recuperado utilizando de los distintos métodos empíricos. Por lo
que, haciendo uso de todo lo anterior, se realizó un análisis completo del
comportamiento de cada pozo presentando las discusiones y resultados.
Adicionalmente, en el trabajo se incluyen tres apéndices. El apéndice A muestra la
derivación de la ecuación de difusión, la cual se deriva de la combinación de la
ecuación de continuidad, la ley de flujo y la ecuación de estado. Este modelo
representa el flujo de fluidos en el medio poroso a través del tiempo y espacio. En
el apéndice B se señalan las diferentes propiedades fisio-químicas de los gases
naturales, los cuales son fundamentales para entender la interacción del
movimiento de los fluidos en el medio poroso (sistema roca-fluidos), además se
presentan algunos métodos para estimar algunas de las propiedades físico-
químicas requeridas en el cálculo de la ingeniería de yacimientos y producción de
gas. Por último, el apéndice C presenta la derivación de las soluciones analíticas
para los pozos productores de gas con las condiciones de presión de fondo fluyendo
constante,𝑝𝑤𝑓 , y gasto de gas constante,𝑞𝑔 ,incluyendo, los modelos de flujo y el
desarrollo de las ecuaciones de interpretación de los distintos periodos de flujo para
calcular algunas propiedades del yacimiento.
2
2. Mostrar las ecuaciones de interpretación para los diferentes regímenes de
flujo presentados en yacimientos de baja permeabilidad, las cuales fueron
desarrolladas y adaptadas para flujo de gas real a partir de soluciones
analíticas, semi- analíticas y numéricas considerando diferentes regímenes
de flujo bajo las consideraciones a flujo constante de gas y presión de fondo
fluyendo constante.
3. Presentar los fundamentos de las diferentes metodologías empíricas y semi-
analíticas de análisis de presión-producción en yacimiento no
convencionales.
4. Desarrollar un software especializado en el análisis de la presión-producción
de yacimiento de gas de baja permeabilidad para el uso de la comunidad
estudiantil.
5. Proponer una metodología general que incluya tanto las modelos semi-
analíticos y empíricos mostrados en este trabajo, para su mejor rendimiento,
según sean las características del yacimiento analizado.
6. Realizar el análisis de distintos pozos productores con datos validados de
presión y producción, tanto de arenas compactas o de baja permeabilidad
(Tight), como de lutitas gasíferas (Shale) mostrando el procedimiento de
análisis de las metodologías usadas, además de mostrar el rendimiento de
cada método en los distintos pozos de campo usados.
Dentro de los resultados más relevantes obtenidos, se tienen:
1. Exposición didáctica de la teoría de los yacimientos de gas de baja
permeabilidad y las metodologías de análisis de presión-producción
empíricas y semi- analíticas.
2. Metodologías analíticas y empíricas fiables en el análisis de datos de
pozos productores de gas de baja permeabilidad y lutitas gasíferas.
3. Buen rendimiento del software programado con las metodologías de
análisis de presión producción de yacimientos no convencionales para el
uso didáctico en el proceso de enseñanza y aprendizaje.
4. Recomendaciones obtenidas en el análisis de los casos de campo
analizados en este trabajo.
3
(tight y shale), así como de conceptos importante del flujo de fluidos en el medio
poroso.
En el Capítulo III, se muestra las metodologías las diferentes gráficas de diagnóstico
y específicas de análisis para la identificación y caracterización (interpretación),
respectivamente, cada régimen de flujo detectado en los datos de producción de los
pozos productores en los yacimientos de gas, empleando ecuaciones analíticas
respectivas para el cálculo de algunos parámetros del yacimiento. También se
muestra la técnica de superposición del tiempo y los diferentes escenarios de flujo
que se presentan en los yacimientos de gas no convencionales.
El Capítulo IV , se muestra el desarrollo teórico y la metodología de uso de los
métodos de análisis de presión-producción empíricos: modelos de Duong (Duong,
2010), Modelo de declinación Exponencial Extendida (Valko & lee, 2010), Modelo
Modificado de Declinación Extendida (Yu, Lee, Miocevic, Li, & Harris, 2013), Modelo
de potencias (IIk, Rushing, A., & Blasingame, 2008a) y el Modelo de Crecimiento
Logístico (Clark, 2011).
En el Capítulo V, se aplican las metodologías, mostradas en los Capítulos III y IV,
en seis diferentes campos de formaciones de arenisca compactas (baja
permeabilidad o Tight) y de lutitas (Shale), donde se analizan los diferentes
regímenes de flujo mediante graficas de diagnóstico y posteriormente se hace un
análisis individual de cada régimen de flujo con graficas especializadas, después se
hace uso de las ecuaciones de modelos analíticos para el cálculo de algunas
propiedades del yacimiento. Por último, se muestra el cálculo de la recuperación de
gas de los pozos mediante el ajuste de curvas de declinación y se exponen las
conclusiones de los resultados que arrojaron los modelos de análisis de presión-
producción.
En el Capítulo VI, se presenta las conclusiones obtenidas, limitaciones del uso de
los modelos de análisis de presión- producción y algunas recomendaciones en el
uso de los modelos mostrados en el trabajo.
Finalmente se presenta la nomenclatura, referencias y algunos apéndices
desarrollados a lo largo del trabajo.
4
2. CAPÍTULO II ANTECEDENTES Y REVISIÓN BIBLIOGRÁFICA
2.1. Introducción
El gas se considera un combustible de transición que va a ir sustituyendo a otros
combustibles fósiles más contaminantes, por ese motivo se espera que su consuma
crezca hasta 2030, quizás incluso hasta [Link] tal motivo es necesario el
desarrollo de nuevas tecnologías para el óptimo desarrollo de los campos de gas
convencionales y no convencionales. Aunque la mayor parte de la tecnología para
los pozos de aceite se aplican a los pozos de gas, existen diferencias a resaltar, es
por ello conocer estas diferencias y tener un buen conocimiento de los fundamentos
de la ingeniería de yacimientos de gas y de los yacimientos no convencionales.
En la industria es fundamental el conocimiento y buen manejo de los conceptos
fundamentales de las propiedades que interactúan en el yacimiento, dado que un
buen manejo de estos podrá darnos el cómo se comporta nuestro yacimiento y
poder tomar las mejores decisiones en su optima explotación.
En el siguiente capítulo se presentan las definiciones, características y propiedades
del gas natural, así como los conceptos más importantes que se presentan en los
yacimientos de formaciones de armenias compactas (Tight) y de lutitas (Shale), así
como los conceptos referentes al fracturamiento Hidráulico, a los distintos
regímenes de flujo de fluidos y los conceptos del flujo de fluidos en el medio poroso.
2.2. Gas natural.
EL gas natural es definido como un fluido homogéneo de baja viscosidad y densidad
que no tiene volumen definido y toma la forma del recipiente que lo contiene.
Generalmente, el gas natural es una mezcla de hidrocarburos y gases no
hidrocarburos. Los gases hidrocarburos que normalmente son encontrados es
metano(𝐶𝐻4 ), etano(𝐶2 𝐻4 ), propano(𝐶3 𝐻8 ), butano(𝐶4 𝐻10 ), pentano(𝐶5 𝐻12 ) y
pequeñas cantidades de hexano(𝐶6 𝐻14 ) y fracciones pesadas(𝐶7+ ). Los gases no
hidrocarburos incluyen al Dióxido de carbono(𝐶𝑂2 ), sulfuro de hidrogeno(𝐻2 𝑆) y
nitrógeno(𝑁2 ).El gas natural se encuentra como gas húmedo o gas seco, lo cual
depende de la cantidad de hidrocarburos condensables presentes en la mezcla, el
pentano y más pesados son considerados como hidrocarburos condensables, a
temperatura y presión atmosférica el pentano es un líquido, mientras que el metano,
etano, propano y butano existen en fase gaseosa a condiciones atmosféricas.
5
De los componentes del gas el más común es el 𝐶1 que generalmente varia su
porcentaje entre 70% al 98% mol, seguido del 𝐶2 entre 1% al 10% y los que
+presentan menor porcentaje son el 𝐶6 (trazas a 0.5%) y el 𝐶7+ que generalmente
presenta trazas.
Los componentes no hidrocarburos son los siguientes: 𝑁2 , 𝐶𝑂2 , 𝐻2 𝑆 y Helio. El
sulfuro de hidrogeno y dióxido de carbono se denominan gases ácidos porque ellos
se disocian en agua para generar soluciones acidas. El sulfuro de Hidrogeno es muy
toxico y corrosivo, mientras el dióxido de carbono es únicamente corrosivo.
No Hidrocarburos
Nitrógeno Trazas-15
Bióxido de Carbono Trazas-2
Ácido Sulfhídrico Trazas-2
Helio Hasta el 5%, generalmente trazas o no
existe.
6
En el Apéndice B se presentan algunas propiedades fisicoquímicas de los gases
naturales y los métodos para estimar algunas de sus propiedades fisicoquímicas
requeridas para los cálculos en la ingeniería de yacimientos de gas.
7
3. Gas condensado. El gas condensado tiene alta cantidad de
hidrocarburos líquidos y puede ocurrir como gas en el yacimiento. La
fuente más deseable puede ser de tipo no asociado porque este puede
ser producido a alta presión. El gas asociado o disuelto es separado
del aceite crudo a baja presión en el separador y además exige más
gastos de compresión. El gas condensado representa una gran
cantidad de gas asociado con el líquido.
8
En la Fig. 2-1 podemos observar como las curvas de calidad tienden a
cercarse más a la curva de presión y burbuja, estando presente muy poco
liquido en la región de dos fases cerca de la curva de presión de roció. Al
describir el comportamiento de una muestra analizada en el yacimiento
comenzando por el punto 𝑃1 (donde la presión es mayor a la presión de rocio)
el hidrocarburo se encuentra en fase gaseosa, disminuyendo la presión llega
al punto 𝑃2 donde apensar de estar por debajo de la cricondenbara siguen
estando en fase gaseosa, pero en este punto entra en la tubería de
producción hasta llegar a la batería de separación(en este trayecto existe una
disminución en la presión y en la temperatura) que es el punto 𝑃3 , en este
punto ya se han condensado algunos de los componentes mas pesados de
la mezcla estando presentes las dos fases. Por lo tanto, no se formará liquido
en el yacimiento, el termino húmedo se deriva de que en la tubería de
producción o en la batería de separación se entra en las dos fases, la
formación de líquido es causada porque la energía cinética de las moléculas
disminuye lo suficiente como consecuencia de la disminución de la
temperatura, cambiando a liquido por las fuerzas atractivas entre las
moléculas.
El gas seco está compuesto en su mayor parte por metano y etano con
pequeños porcentajes de componentes pesados (en su composición también
pueden estar presentes no hidrocarburos como nitrógeno y dióxido de
carbono). El gas seco contiene vapor de agua, que sería un condensado.
Seco significa libre de hidrocarburos líquidos mas no libre de agua, la
temperatura de los yacimientos que contienen este tipo de fluidos es mucho
mayor que la cricondenterma, en un diagrama de fase. Las curvas de calidad
se acercan mucho más a la curva de burbuja cargándose las concentraciones
de líquidos a esta, no se presentan condensados en el yacimiento. La energía
cinética es tan alta comparada con las fuerzas de atracción por lo que no se
condensa la mezcla.
En la Fig. 2-2 el caso 𝑃1 ,𝑃2 y 𝑃3𝑎 el gas no entra en la región de dos fases,
por lo tanto, los hidrocarburos no se condensan, pero en el caso 𝑃1 ,𝑃2 y 𝑃3𝑏
si se entra en la región de dos fases, debido al vapor de agua presente en la
mezcla.
9
Ilustración 2-2 Diagrama de fase(presión-temperatura) para un yacimiento de gas
seco (Arévalo J. D., 2009).
10
Ilustración 2-3 Yacimientos Convencionales y No Convencionales (Administration,
2017).
11
El contenido de materia orgánica y el gas adsorbida son los factores que gobiernan
el comportamiento los yacimientos de gas no convencionales con baja
permeabilidad. (F., 2014). El comportamiento puede ser representado a través de
modelos conceptuales tomando los siguientes conceptos en consideración:1)
Mecanismo de almacenamiento en yacimientos de gas no convencionales,2)
transporte en yacimientos de gas no convencionales, y 3) adsorción física del gas y
efectos de desorción en yacimientos no convencionales.
2.3.1. Mecanismo de almacenamiento en yacimientos de gas no
convencionales.
Hay dos principales tipos de gas almacenado en yacimientos de gas no
convencionales: gas libre en los poros de la matriz y gas adsorbido en la superficie
de la matriz. Un modelo teórico de triple-porosidad incluye ambos, gas libre y gas
adsorbido el cual es usado para considerar todo el gas que esta almacenado en la
formación que contiene la materia orgánica. (Arevalo-Villagran J.A. C. P., 2014)
Un sistema de triple-porosidad es una combinación de un sistema de doble
porosidad, fracturas en la matriz, y gas adsorbido. Los microporos de la matriz son
considerados como porosidad 1, las fracturas naturales como porosidad 2, y el gas
adsorbido es la porosidad 3, Fig. 2-4, sin embargo, el almacenamiento del gas
adsorbido no está en la superficie de las partículas de la formación de la matriz.
12
Ilustración 2-4 Modelo de almacenamiento de triple porosidad para yacimientos de
gas no convencionales (Arévalo Villagran, Catellanos Páez, Martínez Romero, & Pumar
Martínez , 2014).
13
Estos mecanismos pueden actuar juntos o separados, ya sean las propiedades del
gas y de la roca. El mecanismo que domina el transporte en la porosidad primaria
de la matriz obedece la primera Ley de Fick (gradiente de concentración). El flujo
de Darcy no es alcanzado debido al hecho de que la permeabilidad es muy baja,
Ec. 2.1.
−𝐷 𝐴𝑧𝑠𝑐 𝑅 𝑇𝑠𝑐 𝑑𝐶 (2.1)
𝑞𝑔 =
𝑃𝑠𝑐 𝑑𝑥
14
2.3.3. Adsorción física de gas y desorción en yacimiento no
convencionales.
Los yacimientos no convencionales presentan un mecanismo de almacenamiento
diferente a los convencionales, el cual es el fenómeno de adsorción de las
moléculas de gas de las paredes de la roca orgánica. El proceso de este fenómeno
ocurre cuando la estructura electrónica del átomo o de las moléculas es apenas
perturbada sobre adsorción. En otras palabras, las moléculas adsorbidas conservan
su naturaleza química. (Arevalo-Villagran J.A. C. P., 2014)
El modelo de Langmiur ha sido el modelo más usado para describir el fenómeno de
la adsorción de gas en sólidos. El cual considera que las moléculas de gas son
adsorbidas en un único lugar y no afecta las moléculas vecinas. El modelo de
Langmiur puede ser representada por la Ec. 2.2:
𝑉𝐿 𝑝 ( 2.2)
𝑉𝑎 =
𝑝𝐿 + 𝑝
Donde 𝑉𝑎 es el volumen total del gas adsorbido por unidad de volumen que está en
equilibrio en el yacimiento a la presión 𝑝; 𝑉𝐿 es el volumen de Langmiur o el volumen
máximo adsorbido por unidad de volumen en el yacimiento a una presión infinita, y
𝑝𝐿 es la presión Langmiur que representa la presión a la cual el volumen adsorbido
𝑉𝑎 es igual a la mitad del volumen Langmiur, 𝑉𝐿 .
Las isotermas de Langmiur describe la máxima cantidad de gas en un yacimiento
no convencional que puede estar almacenado dentro ciertas condiciones con
contenido de materia orgánica, presión y temperatura, Fig 2-6.
15
2.3.4. Gas de lutita (Shale gas).
16
El concepto de triangulo de recursos fue usado por Masters (1979), el concepto es
que todos los recursos naturales están distribuidos normalmente en la naturaleza.
Como vas descendiendo en el triángulo de recursos-gas, los yacimientos son de
bajo nivel, lo que significa que la permeabilidad del yacimiento está decreciendo,
Fig. 2-7. Sin embargo, el tamaño es mayor que los yacimientos de mejor calidad.
El tema principal es que los depósitos de baja calidad de gas natural requieren
mejores tecnologías y adecuados precios de gas antes de ser económicamente
desarrollados y producidos
Ilustración 2-7 Triangulo de Recursos (Usman Ahmed). (Holditch 2006, Masters 1979
y Gray 1977).
17
2.4. Fractura miento hidráulico.
La capacidad de un pozo para producir hidrocarburos o recibir fluidos de inyección
es limitada por la permeabilidad natural del yacimiento y los cambios producidos en
la región vecina al pozo como resultado de las operaciones de perforación o de otro
tipo. El fracturamiento hidráulico, también conocido como estimulación hidráulica,
mejora el flujo de hidrocarburos mediante la creación de fracturas en la formación,
que conectan el yacimiento con el pozo. (Nolen-Hoeksema, 2013)
Una fractura hidráulica es una fractura inducida por presión, causada por la
inyección de fluido altamente viscoso en una formación rocosa objetivo. El fluido es
bombeado en a la formación a presiones que exceden la presión de fractura miento;
la presión a la cual se fracturan las rocas. Para ingresar a una zona para el
tratamiento de estimulación, los ingenieros disparan la tubería de revestimiento a
través del intervalo de interés y utilizan tapones recuperables para asilar dicho
intervalo de otras zonas abiertas. Luego, este intervalo se presuriza hasta alcanzar
la presión de ruptura de la formación, o presión de iniciación de la fractura, punto en
el cual la roca se rompe y se forma una fractura.
18
Los esfuerzos locales controlan la orientación y la dirección de propagación de las
fracturas hidráulicas, Fig. 2-8. Las fracturas hidráulicas se abren en la dirección de
menor resistencia. Los tres esfuerzos principales se incrementan con la
profundidad. El esfuerzo vertical principal, conocido generalmente como esfuerzo
de sobrecarga, es causado por el peso de la roca que subrayase un punto de
medición. Su gradiente vertical se denomina litostático. Los esfuerzos horizontales
mínimo y máximo son los otros dos esfuerzos principales. Sus gradientes verticales,
que varían debido a la cuenca y litología, son controlados por los esfuerzos locales
y regionales, principalmente a través de la litología.
19
El efecto de incremento de drene de fluidos decrece rápidamente con el tiempo.
Esto se debe a que la fisura se cierra y el pozo vuelve a sus condiciones casi
originales. Para evitar que el cierre de la fractura, los ingenieros ingresan una
lechada de fracturamiento compuesta de fluido y apuntalante; solido que mantiene
abierta la fractura después de que se interrumpe el bombeo, Fig 2-10.
20
ecuación, la cual se desarrolló originalmente para el estudio de transferencia de
calor. El flujo de fluidos a través de medios porosos es directamente análogo al flujo
de calor a través de sólidos. La siguiente expresión es la ecuación de difusividad
para coordenadas radiales.
1 𝜕𝑟 𝜕𝑝 𝜙 𝜇 𝑐𝑡 𝜕𝑝 (2. 3)
(𝑟 ) =
𝑟 𝜕 𝜕𝑟 𝑘 𝜕𝑡
Ahora para expresar la Ec. 2.5 en forma similar a la ecuación de difusividad para
flujo de líquidos, Ec.2. 4. Definimos la función pseudo- presión como:
𝑝
𝑝 (2. 6)
𝑝𝑝 (𝑝) = 2 ∫ 𝑑𝑝
𝑝𝑏 𝜇𝑧
Para reportar las variaciones en las propiedades en el gas con el tiempo, definimos
el pseudo- tiempo:
𝑡
𝑑𝑟 (2. 7)
𝑡𝑎𝑝 (𝑡) = ∫
0 𝜇𝑐𝑡
21
1𝜕 𝜕𝑝 𝜙 𝜇 𝑐𝑡 𝜕𝑝 (2. 8)
(𝑟 ) =
𝑟 𝜕𝑟 𝜕𝑟 0.0002637 𝑘 𝜕𝑡
Dado que la dependencia de las propiedades del gas por efecto de la presión, Ec.
2.10, es no-linear, es una ecuación diferencial parcial. Si asumimos que el
termino 𝑝⁄𝜇𝑧 es constante con respecto a la presión y el termino 𝜇𝑐𝑔 es evaluado
con una presión promedio y se trata como constante se resuelve la Ec. 2.5 en
términos de la presión (Ec. 2.6). Generalmente, esas suposiciones son válidas
únicamente por altas temperaturas y presiones, Fig. 2-12 a Fig. 2-14,las cuales
muestran la relación entre presión y 𝑝⁄𝜇𝑧 a diferentes temperaturas y varias
gravedades específicas de gas, como se observa en la ilustración 𝑝⁄𝜇𝑧 es constante
con la presion a temperaturas superiores a 3000 [psia] para 100°F,5000 [psia] para
200°F y 6500 [psia] para 300°[Link] conclusión, a valores altos en cuanto a la
gravedad especifica del gas y a la presión , los valores de 𝑝⁄𝜇𝑧 varían en función
con la presión. Las figuras implican que las soluciones para la ecuación de difusión
para el gas real, las cuales asumen constante a 𝑝⁄𝜇𝑧, se deben utilizar únicamente
para gases a muy altas presiones. (Wattenbarger J. L., 1996)
Considerando que el producto 𝜇𝑧 es constante con la presion y que el producto 𝜇𝑐𝑔
se evalúa con un promedio de la presión,𝑝̅ ,se resuelve la Ec. 2.5 en términos de la
presión al cuadrado Ec. 2.10. La hipótesis de que el producto 𝜇𝑧 es constante, es
válida solamente para bajas gravedades específicas y presiones con altas
temperaturas. El esbozo de la Fig.2-15 a Fig.2-17 ilustran la variación del producto
𝜇𝑧 con respecto a la presión para diferentes gravedades específicas del gas y
temperaturas de 100,200 y 300°F, respectivamente.
Se nota que el producto 𝜇𝑧 es constante en función de la presión a presiones
menores que 1200 psia para 100°F, 1750 psia para 200°F y 2200 psia para 300 °F.
Con valores altos para la gravedad especifica del gas, el producto 𝜇𝑧 varía en
función de la presión, entonces las soluciones para la ecuación de difusión deben
ser utilizadas en términos de la presión al cuadrado solamente a muy bajas
presiones y gravedades específicas de gas, con altas temperaturas.
22
Las consideraciones que se toman en cuenta para la derivación de la ecuación de
difusividad son razonables para altas presiones en gases (aquellos con presiones
mayores a 3000 psia).Para gases a bajas presiones (aquellos con menores a 3000
psia), la hipótesis de que la compresibilidad y viscosidad sea constante no es válida,
pero esto se maneja mediante el concepto de “pseudo- presión del gas real” la cual
es una ecuación para yacimientos de gas, con el objeto de obtener una mayor
precisión para todos los rangos de presiones. Este método riguroso para linealizar
la ecuación diferencial parcial Ec. 2.5 es con la función de la pseudo- presión del
gas real definida por Al-Hussainy.
𝑝
𝑝 (2. 10)
𝑚(𝑝) = 2 ∫ 𝑑𝑝
𝑝𝑜 𝜇𝑔 𝑧
La función de la pseudo- presión del gas real integra las variaciones de 𝑝, 𝑧 𝑦 𝜇 con
la presión, lo que permite que la Ec. 2.5 se resuelva sin la limitación de la suposición
de ciertas propiedades como constantes en función de ciertos valores de la presión.
Se introduce esta función dentro de la derivación de la ecuación de difusión, en
términos de la pseudo- presión y el tiempo, la ecuación de difusión para el gas
queda:
1𝜕 𝜕𝑚(𝑝) 𝜙 𝜇 𝑐𝑡 𝜕𝑚(𝑝) (2. 11)
(𝑟 )=
𝑟 𝜕𝑟 𝜕𝑟 0.0002637 𝑘 𝜕𝑡
23
Ilustración 2-14 Rango de aplicación 2-13 Rango de aplicación de los
de los métodos de presión a 300°F. métodos de presión al cuadrado a
(Wattenbarger J. L., 1996) 100°F. (Wattenbarger J. L., 1996)
Ilustración 2-16 Rango de aplicación del Ilustración 2-15 Rango de aplicación del
método de presión al cuadrado a 200°F. método de presión al cuadrado a 300°F.
(Wattenbarger J. L., 1996) (Wattenbarger J. L., 1996)
24
El estado de flujo estacionario raramente puede obtenerse en la realidad por que
los pozos de gas usualmente producen gas de un yacimiento finito limitado o de un
volumen de drene. No puede haber flujo a través de limite exterior. Así, la presión
declina como la producción continua. Un verdadero estado estacionario requeriría
presión para mantener constante en el límite exterior, lo que implica flujo a través
del límite exterior. (R. A.-h. , 1965)
25
𝜕𝑝 (2.18)
( ) = 𝑓(𝑖, 𝑡)
𝜕𝑡 𝑖
Aronofsky y Jenkins demostraron que la producción de un gas ideal de un sistema
cerrado radial puede ser aproximadamente cercano a la solución para flujo de
líquidos en estado transitorio de van Everdingen y Hurst. (R. A.-h. , 1965)
La ecuación para flujo radial de un gas real está dada por la Ec. 2.19:
Ilustración 2-17 Grafica de 1/p contra pseudo presión con viscosidad del gas ideal y
gas real (Wattenbarger J. L., 1996).
En vista del parecido entre (𝜇𝑐𝑔 )𝑣𝑠 𝑚(𝑝) para gas real, y 𝑝−1 𝑣𝑠 𝑝2 para gases
ideales. Es razonable esperar soluciones para flujo de gases reales para
correlacionar como funciones de tiempo adimensional basado en valores iniciales
de viscosidad y compresibilidad. Que nos deja:
𝑘𝑡 (2.20)
𝑡𝐷 = 2
𝜙 (𝜇𝑐𝑔 )𝑖 𝑟𝑠𝑐
26
Además, definimos la caída de presión adimensional para gas real 𝑚𝐷 (𝑟𝐷 , 𝑡𝐷 ):
𝜋 𝑘 ℎ 𝑇𝑠𝑐 [𝑚(𝑝𝑖 ) − 𝑚(𝑟, 𝑡)] (2.21)
𝑚𝐷 (𝑟𝐷 , 𝑡𝐷 ) =
𝑞𝑠𝑐 𝑝𝑠𝑐 𝑇
Donde 𝑟𝐷 = 𝑟⁄𝑟𝑠𝑐 .La caída de presión adimensional de la pseudo-presión del gas
real es semejante a la caída de presión adimensional de van Everdingen-Hurst (van
Everdingen, 1949),𝑝𝐷 (𝑟𝐷 , 𝑡𝐷 ).
27
Ilustración 2-19 𝒑𝑫 (𝒕𝑫 ) 𝒚 𝒎𝑫 (𝒕𝑫) 𝒗𝒔 𝒕𝑫 para flujo radial para líquido, gas ideal y gases
reales. (Wattenbarger J. L., 1996)
En los casos de gas ideal y gas real de la Fig. 2-20, conduce a caídas de presión
adimensional, las cuales son menos que el caso del líquido y hay dependencia de
la tasa de flujo. Otra diferencia importante es en el caso de 𝑄 = 0.05. La línea de
gas ideal termina en el punto donde la presión del pozo es cero. La solución del gas
real termina cuando la presión del pozo está al 10 por ciento de la presión inicial.
Sin embargo, la Fig. 2-20 no muestra que tiempo de producción para los casos de
gases reales para alcanzar la presión de producción limitada es cerca de dos y la
mitad de los tiempos requeridos para los casos de los gases ideales.
La línea de gas natural está cerca del caso de los líquidos, pero la línea de
condensados es muy debajo del caso de los líquidos. La presión de producción
terminal se alcanza más rápido en el caos de los condensados que en la línea del
caso de gas natural. Esto indica la importancia de la variación de las propiedades
de gas encima del resultado.
La correlación de la Fig. 2-19 y Fig. 2-20 incluyen la dependencia de la
permeabilidad a la presión. Si la una aproximación de los efectos de la caída de
presión sobre la condensación de líquidos y reducción en la permeabilidad cerca
del pozo puede ser hecha, el rendimiento puede estimarse de:
𝑟𝑑 𝜋 𝑘 ℎ 𝑇𝑠𝑐 (2.22)
ln = [𝑚´(𝑝̅ ) − 𝑚´(𝑝)]
𝑟𝑤 𝑞𝑠𝑐 𝑝𝑠𝑐 𝑇
28
Donde
𝑝
𝑘 𝑝 𝑑𝑝 (2. 23)
𝑚´(𝑝) = 2 ∫
𝑝𝑚 𝜇 𝑧
Y k está en función de la presión.
Eilerts determino el comportamiento con estado estacionario, en la región de flujo
no-Darcy cerca del pozo productor. Como resultado, el efecto skin del estado
estacionario puede ser introducido para darla siguiente aproximación para el flujo
radial de gases reales durante la producción:
𝜋𝑘ℎ𝑇𝑠𝑐 [𝑚(𝑝̅ ) − 𝑚(𝑝)] 𝑟𝑑 (2.24)
= ln + 𝑠 + 𝐷𝑞
𝑞𝑠𝑐 𝑝𝑠𝑐 𝑇 𝑟𝑤
Donde s es el efecto skin y D es el coeficiente de flujo no-Darcy.
29
Ilustración 2-20 Relación Presión-tiempo de los distintos flujos de fluidos (Ahmed,
Resevoir Engineering Handbook, 2006).
En el caso de los líquidos, una ecuación similar a la Ec. 2.26 puede ser derivada
usando el concepto de estado de flujo pseudo- estacionario. La cual es la condición
que es alcanzado por producción a flujo constante de líquido cuando la tasa de
presión declina se vuelve constante en cada parte del yacimiento. Esta condición
es expresada matemáticamente usando Laplaciano de la presión igual a la
constante. Sin embargo, el Laplaciano de la pseudo- presión del gas real no puede
ser igual a una constante, una condición de flujo similar al estado de flujo pseudo-
estacionario parece existir para flujo de gas ideal y flujo de gas real. La existencia
de tales como condición esta sugerida por la Ec. 2-26 y la Fig. 2-20, la cual presenta
una inspección del comportamiento de la presión durante el periodo que la Ec. 2-26
se adapta al caso de Eilerts.
Como se vio en la Fig. 2-19, 𝑚𝐷 (𝑡𝐷 ) no cambia a una tasa constante durante ese
periodo. Sin embargo, se iguala a la solución del líquido a tiempo tempranos. Lo
más interesante de la Fig. 2-20 es que 𝑚𝐷 (𝑟𝐷 , 𝑡𝐷 ) para todas las ubicaciones
radiales es esencialmente paralelo. Así, el perfil de 𝑚(𝑝) es independiente del
30
tiempo. Esta condición puede ser descrita por el Laplaciano de 𝑚(𝑝) igual a una
constante. En cualquier acontecimiento, la Ec. 2.27 describe el comportamiento de
flujo radial con un sistema radial a tiempo largos con extraordinaria exactitud.
Otra consecuencia de la Fig. 2-20 es que la distribución de 𝑚(𝑝) puede ser obtenida
fácilmente. Por ejemplo, la siguiente ecuación también describe razonablemente el
flujo:
𝑟𝑒 𝜋𝑘ℎ𝑇𝑠𝑐 (2.27)
ln (0.606 )= [𝑚(𝑝̅ ) − 𝑚(𝑝𝑠𝑐 )]
𝑟𝑤 𝑞𝑠𝑐 𝑝𝑠𝑐 𝑇
31
Ilustración 2-21 Distribución de la presión en la formación cercana del pozo productor
(Wattenbarger J. L., 1996).
El volumen original de gas se define como la cantidad de gas que se estima existen
inicialmente en un yacimiento y se encuentra en equilibrio a la temperatura y presión
prevaleciente en el yacimiento. Este volumen incluye a las acumulaciones
conocidas, económicas o no, recuperables o no, a la producción obtenida de los
campos explotados o en explotación también a las cantidades estimadas en los
yacimientos que podrían ser descubiertos. Este volumen original de gas, G, puede
inferirse por procedimientos determinísticos y probabilísticos, los primeros incluyen
principalmente a los volumétricos, balance de materia, análisis de producción y los
de simulación numérica. Los segundos modelan la incertidumbre de parámetros
como porosidad, saturación de agua, espesores netos, entre otros, como funciones
de probabilidad que producen en consecuencia una función de probabilidad para el
volumen original.
32
2.6. Análisis de tipo curva.
Para el análisis de los modelos de forma correcta, es necesario tomar en cuenta los
siguientes puntos:
A. Flujo dominado por la frontera.
B. Condiciones de operación estables (gasto de flujo y presión de fondo
fluyendo constante).
C. Desarrollo de ecuaciones de forma empírica.
D. Determinación de la recuperación final de forma rápida y sencilla.
E. No se requiere de datos de presión.
F. b depende del mecanismo de producción.
33
La Tabla 2-2 muestra los tres modelos de declinación y sus valores de b
correspondientes, además de las ecuaciones de gasto y de producción acumulada.
Tabla 2-2 Modelos de declinación
34
2.7. Yacimiento Naturalmente Fracturados
35
Donde 𝑘𝑚𝑎 = permeabilidad de la matriz y 𝑘𝑓 = permeabilidad de la fractura. El
parámetro 𝛼, el cual caracteriza la geometría del sistema, es definido como:
4𝑗(𝑗 + 2) (2.30)
𝛼= 2
𝐿
Donde 𝐿 =característica adimensional del bloque de la matriz y 𝑗 = numero de
conjuntos de planos que limitan el medio menos permeable.
El segundo parámetro de doble porosidad es la tasa de almacenamiento, el cual es
definido como:
(𝜙𝑉𝑐𝑡 )𝑓 (𝜙𝑉𝑐𝑡 )𝑓 (2.31)
𝜔= =
(𝜙𝑉𝑐𝑡 )𝑓+𝑚𝑎 (𝜙𝑉𝑐𝑡 )𝑓 + (𝜙𝑉𝑐𝑡 )𝑚𝑎
Donde 𝑉 =la razón de volumen total de un medio al BULK volumen total del sistema
y 𝜙 = tasa de PV de un medio al volumen total del medio. Los sufijos 𝑓 + 𝑚𝑎 se
refieren a la fractura y el sistema total (fractura más matriz).
Barenblatt y más autores, asumen un estado de flujo pseudo- estacionario; mientras
otros, como deSwaan, asumen un estado de flujo transitorio en la matriz. Gringarten
señala que, a pesar de que muchas pruebas de presión exhiben un comportamiento
de Inter porosidad de flujo en estado pseudo- estacionario; otras pruebas sugieren
comportamiento de Inter porosidad de flujo transitorio. Una posible explicación es
que dentro del flujo Inter poroso siempre ocurre bajo condiciones transitorias, pero
pueden exhibir un estado pseudo- estacionario si un impedimento de flujo existe
entre el medio menos permeable(matriz) y el medio más permeable(fractura).
36
3. CAPÍTULO III DIAGNÓSTICO Y CARACTERIZACIÓN DE REGÍMENES
DE FLUJOS: ECUACIONES Y GRÁFICAS
3.1. Introducción
En los yacimientos de baja permeabilidad de arenas compactas y lutitas gasífera se
han identificado que presentan un régimen de flujo transitorio tardío (Lineal, Bilineal,
radial y esférico) durante meses, años e incluso décadas de su vida productiva del
pozo, hasta alcanzar las fronteras (pseudo- estacionario). Estos patrones de flujo
transitorios se pueden llegar a presentar en diferente forma en el yacimiento y su
identificación es crucial para la caracterización correcta del yacimiento y poder hacer
uso eficiente de las distintas herramientas.
En este capítulo, se presentan las diferentes graficas de diagnóstico que podemos
′
∆𝑚(𝑝) ∆𝑚(𝑝)⁄
usar para su identificación en el yacimiento ( 𝑞 𝑣𝑠 𝑡, [ 𝑞𝑔 𝑣𝑠 𝑡 y 𝐺𝑝 𝑣𝑠 𝑡),
]
𝑔
sus modelos de flujo, las graficas especializadas de flujo con sus ecuaciones
analíticas para el cálculo de las distintas propiedades del yacimiento ,tales como, el
volumen original,𝐺, las reservas probadas o recuperables de gas, la conductividad
de la formación,𝑘ℎ, la permeabilidad de la formación,𝑘, el área de la sección
transversal,𝐴𝑐 , el daño de la formación,𝑠, el efecto del régimen de flujo, el área de
drene del pozo,𝐴, el volumen poroso asociado al área de drene,𝑉𝑝 . Además, se
muestra la técnica de superposición de tiempo, la cual tiene como finalidad de ser
utilizada cuando las condiciones de presión de fondo fluyendo constante,𝑃𝑤𝑓 ,y gasto
constante de gas,𝑞𝑔 , presentas variaciones sustanciales respecto al tiempo de
explotación.
Adicionalmente, el Apéndice A se presenta el desarrollo de las ecuaciones de
interpretación, las cuales se adaptaron para el flujo real de gases a partir de
diferentes soluciones analíticas y semi- analíticas presentadas en la literatura
técnica petrolera con la finalidad de analizar los datos de presión-producción
provenientes de pozos productores reales.
En la Fig. 3-1 tenemos una gráfica de diagnóstico de flujo doble logarítmica del
cambio de la presión y la derivada del cambio de la presión contra el tiempo. Se
tienen tres regiones en el tiempo (al comienzo del tiempo, tiempo intermedio y
tiempos tardíos) a la gráfica de diagnóstico de flujo.
Al comienzo del tiempo de producción, las respuestas del cambio de la presión y la
derivada del cambio de la presión son controladas por el almacenamiento del
agujero del pozo y por efectos cercanos al agujero de este. En esta etapa se
37
incluyen el almacenamiento o llenado del agujero del pozo productor, el factor de
daño, la penetración parcial (la cual es una forma geométrica del factor de daño), la
redistribución de fases y la finita e infinita conductividad hidráulica de las fracturas.
En la etapa intermedia el comportamiento del yacimiento es forma infinita, lo que
significa que los límites del yacimiento no incluyen en este lapso en particular. Los
datos que aparecen en esta región proporcionan la mejor estimación en cuanto a la
permeabilidad del yacimiento.
Finalmente, en la región de tiempo tardíos los efectos de la frontera son lo que
dominan las respuestas en la evaluación de datos. Los efectos que influyen en la
respuesta de la presión ocasionados por la frontera externa son: falla de sello,
yacimientos cerrados, diversos contactos (gas/agua, aceite/ agua y gas/aceite).
La Fig. 3-1 muestra una gráfica de diagnóstico de flujo en la cual se observa una de
las mayores ventajas de la gráfica de diagnóstico de flujo:” Identificar los regímenes
de flujo que dominan la respuesta de la presión durante un tiempo en la evaluación
de datos mediante el valor de sus pendientes. En el caso de la Fig. 3-1 se observa
que al comienzo del tiempo de explotación existe un efecto del almacenamiento del
fluido en el agujero del pozo productor, gráficamente se observa debido a una
pendiente con valor igual a la unidad. En el siguiente corte se observa un régimen
de flujo esférico, con la derivada del cambio de la presión se tiene una pendiente de
la línea recta de -1/2. Después es seguido por un flujo radial, el cual es identificado
por la línea de la derivada del cambio de la presión en forma horizontal. Por último,
en el periodo tardío se puede identificar que se tiene una pendiente igual a la unidad,
lo que indica existe algún tipo de energía externa, la cual es una fuente de “recarga”,
38
por ejemplo, un acuífero que provee de energía suficiente para recargar el
yacimiento.
Ilustración 3-2 Grafica de diagnóstico con sus distintos regímenes de flujo (Arévalo
J. D., 2009).
39
3.3. Graficas De Diagnostico Para Los Regímenes De Flujo
∆𝒎(𝒑)
Ilustración 3-3 Graficas de diagnóstico, 𝒕 𝒚 𝑮𝒑 𝒗𝒔 𝒕, para identificar un régimen
𝒒𝒈 𝒗𝒔
lineal (Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case
histories, 2001).
40
∆𝒎(𝒑)
Ilustración 3-4 Graficas de diagnóstico, 𝒒 𝒕 𝒚 𝑮𝒑 𝒗𝒔 𝒕, para identificar un régimen
𝒈 𝒗𝒔
bilineal (Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case
histories, 2001).
∆𝒎(𝒑)
Ilustración 3-5 Graficas de diagnóstico , 𝒒 𝒕 𝒚 𝑮𝒑 𝒗𝒔 𝒕, para identificar un régimen
𝒈 𝒗𝒔
radial (Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case
histories, 2001).
41
∆𝒎(𝒑)
Ilustración 3-6 Graficas de diagnóstico , 𝒕 𝒚 𝑮𝒑 𝒗𝒔 𝒕, para identificar un régimen
𝒒𝒈 𝒗𝒔
esférico (Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case
histories, 2001).
∆𝒎(𝒑)
Ilustración 3-7 Graficas de diagnóstico , 𝒒 𝒕 𝒚 𝑮𝒑 𝒗𝒔 𝒕, para identificar un régimen
𝒈 𝒗𝒔
esférico a la condición de flujo constante de gas,𝒒𝒈. (Arévalo Villagran, Analysis of
Long-term Behavior in tight gas resevoirs:case histories, 2001)
42
∆𝒎(𝒑)
Ilustración 3-8 Graficas de diagnóstico , 𝒒 𝒕 𝒚 𝑮𝒑 𝒗𝒔 𝒕, para identificar un régimen
𝒈 𝒗𝒔
esférico a la condición de presión de fondo fluyente constante,𝒑𝒘𝒇 . (Arévalo Villagran,
Analysis of Long-term Behavior in tight gas resevoirs:case histories, 2001)
Las Tablas 3-1,3-2 y 3-3 indican la pendiente de la línea recta que se encuentran
′
∆𝑚(𝑝)⁄ ∆𝑚(𝑝)⁄
en las gráficas de diagnóstico de 𝑡 [ 𝑞𝑔 𝑣𝑠 𝑡,
] 𝑞𝑔 𝑣𝑠 𝑡 y 𝐺𝑝 𝑣𝑠 𝑡,
respectivamente.
43
Tabla 3-1 Determinación del régimen de flujo mediante la gráfica de diagnóstico [∆m(p)/qg]'
vs t=ct^n.
44
3.4. Graficas Especializadas Para Cada Régimen De Flujo
Una vez que se caracteriza cada periodo de flujo detectado en las gráficas de
diagnóstico de flujo se elaboran las gráficas especializadas haciendo
[𝑚(𝑝𝑖 )−𝑚(𝑝𝑤𝑓 )]
𝑣𝑠 𝑢𝑛𝑎 𝑓𝑢𝑛𝑐𝑖𝑜𝑛 𝑒𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐𝑎 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜 𝑠𝑒𝑔𝑢𝑛 𝑙𝑎 𝑔𝑒𝑜𝑚𝑡𝑟𝑖𝑎 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜,
𝑞𝑔
estas funciones de tiempo se exponen en la Tabla.3-4, de la línea recta producida
en estas graficas se estima la pendiente, la no intersección en cero y el final de la
misma. Posteriormente se estiman algunos valores del yacimiento, utilizando las
ecuaciones de interpretación para cada régimen de flujo analizado.
Tabla 3-4 Funciones específicas del tiempo para las distintas geometrías de flujo
1
Flujo Lineal 𝑡 ⁄2
Flujo bilineal 1
𝑡 ⁄4
Flujo radial
log(𝑡)
Flujo esférico
Flujo dominado por la 1⁄𝑡 1⁄2
frontera externa 𝑡
45
Ilustración 3-9 Grafica especializada de flujo para régimen lineal. donde podemos
apreciar ambos casos: con intercepción en cero y con intercepción en la recta.
(Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories,
2001)
46
Tabla 3-5 Ecuaciones de interpretación para régimen de flujo lineal en un modelo
homogéneo e infinito para las condiciones de producción a presión de fondo fluyendo
constante,𝑝𝑤𝑓 y a gasto de flujo constante de gas,𝑞𝑔 .
47
Tabla 3-6 Calculo de Área de sección transversal,𝐴𝑐 , y distancia a la frontera,𝐿, según el
tipo de modelo.
Modelo 𝐴𝑐 𝐿
Bloque Lineal 𝑤ℎ 𝐿
Fractura Hidráulica 4 ℎ 𝑥𝑓 𝑦𝑒
Fractura Hidráulica (𝑥𝑒 = 𝑥𝑓 ) 4 ℎ 𝑥𝑒 𝑦𝑒
Pozo productor en un bloque del yacimiento 4 ℎ 𝑥𝑓 𝑦𝑒
Veta de alta permeabilidad con flujo lineal 𝜋 𝑟𝑒2 ℎ
sencillo 2 𝜋 𝑟𝑒2 ℎ⁄2
Veta de alta permeabilidad con flujo lineal doble 2 𝜋 𝑟𝑒2 𝑛𝑣𝑒𝑡𝑎𝑠 ℎ⁄(2 𝑛𝑣𝑒𝑡𝑎𝑠 )
Vetas de n-altas permeabilidad con flujo lineal
doble
Ilustración 3-10 Grafica especializada de flujo para régimen bilineal. donde podemos
apreciar ambos casos: con intercepción en cero y con intercepción en la recta.
(Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories,
2001)
48
Se calcula el termino √𝑘1 𝑤 utilizando la ecuación de interpretación de la Tabla 3-7
para flujo bilineal. Para el caso de flujo bilineal con no intersección en cero de la
línea recta, se evalúa el factor de daño de la formación. Pero la estimación del factor
de daño de la formación 𝑠1 requiere del conocimiento de los valores de k y 𝐴𝑐 .
⁄𝟒
𝟗𝟖𝟒 𝑨𝒄𝟏 𝑻
√𝒌𝟏 𝒘 = ( )
𝟒 𝒎̃ 𝑩𝑮𝑪
√𝒌(𝝓𝝁𝒈 𝒄𝒕 )𝒊
Ilustración 3-11 Grafica de especializada de flujo para régimen radial. donde podemos
apreciar ambos casos: con intercepción en cero y con intercepción en la recta.
(Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories,
2001)
49
Se calcula el producto 𝑘 ℎ (potencial de flujo), el factor de daño,𝑠, la longitud media
de la fractura,𝑥𝑓 , el radio de investigación,𝑟𝑖𝑛𝑣 y el volumen original de gas,𝐺 al radio
de investigación mediante el uso de las ecuaciones descritas en la Tabla 3-8.
El volumen original de gas, G y el radio de investigación,𝑟𝑖𝑛𝑣 se consideran valores
mínimos, si toda la historia de producción aun tiende sobre la línea recta en la gráfica
semilogarítmica. Si el comportamiento una sigue en línea recta, el último tiempo de
producción es utilizado como tiempo final sobre la línea recta,𝑡𝑓𝑙𝑟 ,en las ecuaciones
de la Tabla 3-9.
Tabla 3-8 Ecuaciones de interpretación para régimen de flujo radial en un modelo
homogéneo para la condición de producción a presión de fondo fluyendo constante,𝑝𝑤𝑓 .
𝟏𝟔𝟒𝟎 𝑻
𝒌𝒉 =
̃ 𝑹𝑷𝑪
𝒎
𝟏 𝒃𝑹𝑷𝑪 𝒌 𝒉 𝟎. 𝟎𝟎𝟔𝟑𝟑 𝒌
𝒔= [ − 𝐥𝐨𝐠 ( ) − 𝟎. 𝟑𝟓𝟐]
𝟎. 𝟖𝟔𝟗 𝟏𝟔𝟒𝟎 𝑻 (𝝓𝝁𝒈 𝒄𝒕 ) 𝒓𝟐𝒘
𝒊
𝟎. 𝟎𝟎𝟔𝟑𝟑 𝒌
𝒙𝒇 = √
𝒌𝒉 𝒃
𝐞𝐱𝐩 ( 𝟕𝟏𝟐𝑹𝑷𝑪
𝑻 − 𝟐. 𝟐) 𝝓(𝝁𝒈 𝒄𝒕 )𝒊
𝒌 𝒕𝒇𝒊𝒓
𝒓𝒊𝒏𝒗 = 𝟎. 𝟏𝟕𝟕𝟗 √
(𝝓𝝁𝒈 𝒄𝒕 )
𝒊
𝟏𝟔𝟑 𝑻 𝒕𝒇𝒍𝒓
𝑨= ( )
(𝝓𝝁𝒈 𝒄𝒕 ) 𝒉 𝒎 ̃ 𝑹𝑷𝑪
𝒊
𝟏𝟔𝟑 𝑻 𝒕𝒇𝒍𝒓
𝑽𝑷 = ( )
(𝝁𝒈 𝒄𝒕 ) 𝒎 ̃ 𝑹𝑷𝑪
𝒊
𝟏𝟔𝟑 𝑻𝒔𝒈𝒊 𝒕𝒇𝒍𝒓
𝑮= ( )
(𝝁𝒈 𝒄𝒕 𝑩𝒈 ) 𝒎 ̃ 𝑹𝑷𝑪
𝒊
50
Tabla 3-9 Ecuaciones de interpretación para régimen de flujo radial en un modelo
homogéneo para la condición de producción a gasto de flujo constante de gas,𝑞𝑔 .
𝟎. 𝟎𝟎𝟔𝟑𝟑 𝒌
𝒙𝒇 = √
𝒌𝒉 𝒃
𝐞𝐱𝐩 ( 𝟕𝟏𝟐𝑹𝑮𝑪
𝑻 − 𝟐. 𝟐) 𝝓(𝝁𝒈 𝒄𝒕 )𝒊
𝒌 𝒕𝒇𝒍𝒓
𝒓𝒊𝒏𝒗 = 𝟎. 𝟏𝟕𝟕𝟗 √
(𝝓𝝁𝒈 𝒄𝒕 )
𝒊
𝟏𝟔𝟑 𝑻 𝒕𝒇𝒍𝒓
𝑨= ( )
(𝝓𝝁𝒈 𝒄𝒕 ) 𝒉 𝒎 ̃ 𝑹𝑮𝑪
𝒊
𝟏𝟔𝟑 𝑻 𝒕𝒇𝒍𝒓
𝑽𝑷 = ( )
(𝝁𝒈 𝒄𝒕 ) 𝒎 ̃ 𝑹𝑮𝑪
𝒊
𝟏𝟔𝟑 𝑻𝒔𝒈𝒊 𝒕𝒇𝒍𝒓
𝑮= ( )
(𝝁𝒈 𝒄𝒕 𝑩𝒈 ) 𝒎 ̃ 𝑹𝑮𝑪
𝒊
51
Ilustración 3-12 Grafica especializada para régimen de flujo esférico. (Arévalo
Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories, 2001)
𝟐⁄
𝟏𝟎𝟎𝟗𝟖 𝑻 𝟑
𝒌 = −[ √(𝝓𝝁𝒈 𝒄𝒕 )𝒊 ]
𝒎̃ 𝑬𝑮𝑪
𝟏𝟒𝟐𝟒 𝑻
𝒓𝒆𝒔𝒇 =
𝒌 𝒃𝑬𝑮𝑪
52
Ilustración 3-13 Grafica especializada de flujo para régimen de flujo dominado por la
frontera externa a la condición de producción a gasto de flujo constante de gas,𝒒𝒈 .
(Arévalo Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories,
2001)
𝟓𝟕 𝑻 𝟏
𝑽𝒀 = ( )
(𝝓𝝁𝒈 𝒄𝒕 ) 𝒎 ̃ 𝑭𝑮𝑪
𝒊
𝑽𝒀
𝑨=
𝒉
𝟓𝟕 𝑻 𝟏 𝟏
𝑨= ( )( )
(𝝓𝝁𝒈 𝒄𝒕 ) 𝒎 ̃ 𝑭𝑮𝑪 𝒉
𝒊
𝟐. 𝟐𝟒𝟓𝟖 𝑨
𝑪𝑨 =
𝒃𝑭𝑮𝑪 𝒌 𝒉 𝟐
𝐞𝐱𝐩 [𝟐 ( 𝟏𝟒𝟐𝟒 𝑻 ) − 𝟐𝒔] 𝒓𝒘
𝟓𝟕 𝑻 𝟏
𝑽𝒑 = ( )
(𝝁𝒈 𝒄𝒕 ) 𝒎 ̃ 𝑭𝑮𝑪
𝒊
𝟓𝟕 𝑻 𝑺𝒈𝒊 𝟏
𝑮= ( )
(𝝁𝒈 𝒄𝒕 𝑩𝒈 ) 𝒎 ̃ 𝑭𝑮𝑪
𝒊
53
3.4.6. Régimen De Flujo Dominado Por La Frontera Externa Con La
Condición De Producción A Presión De Fondo Fluyendo
Constante,𝑷𝒘𝒇
Ilustración 3-14 Grafica especializada de flujo para régimen de flujo dominado por la
frontera externa a la condición de producción de presión de fondo fluyendo
constante,𝒑𝒘𝒇 . (Arévalo Villagran, Analysis of Long-term Behavior in tight gas
resevoirs:case histories, 2001)
54
Tabla 3-12 Ecuaciones de interpretación para régimen de flujo dominado por la frontera
externa en un modelo homogéneo para la condición a presión de fondo fluyendo
constante,𝑝𝑤𝑓 (declinación exponencial).
𝟐𝟓 𝑻 𝟏
𝑽𝒀 = ( )
(𝝓𝝁𝒈 𝒄𝒕 ) 𝒃𝑭𝑷𝑪 𝒎 ̃ 𝑭𝑮𝑪
𝒊
𝑽𝒀
𝑨=
𝒉
𝟐𝟓 𝑻 𝟏
𝑨= [ ]
(𝝓𝝁𝒈 𝒄𝒕 ) 𝒃𝑭𝑷𝑪 𝒉 𝒎 ̃ 𝑭𝑮𝑪
𝒊
𝟓𝟔 𝑻 𝟏 −𝟐𝒃𝑭𝑷𝑪 𝒌 𝒉
𝑪𝑨 = [ ] 𝐞𝐱𝐩 [ ]
𝒉 (𝝓𝝁𝒈 𝒄𝒕 ) 𝒃𝑭𝑷𝑪 𝒓𝒘 𝒎
𝟐 ̃ 𝑭𝑮𝑪 𝟏𝟒𝟐𝟒 𝑻
𝒊
𝟐𝟓 𝑻 𝟏
𝑽𝒑 = ( )
(𝝁𝒈 𝒄𝒕 ) 𝒃𝑭𝑷𝑪 𝒎 ̃ 𝑭𝑮𝑪
𝒊
𝟐𝟓 𝑻 𝑺𝒈𝒊 𝟏
𝑮= ( )
(𝝁𝒈 𝒄𝒕 𝑩𝒈 ) 𝒃𝑭𝑷𝑪 𝒎 ̃ 𝑭𝑮𝑪
𝒊
55
[𝑚(𝑝𝑖 )−𝑚(𝑝𝑤𝑓 )]
3.-Flujo Lineal tardío. a) Se crea la gráfica especializada de 𝑣𝑠 √𝑡 y
𝑞𝑔
se identifica la línea recta en los tiempos tardíos como se muestra en la Fig. 3-15.
Después, se calcula la pendiente y la no intersección en cero de la línea recta, si
esta existe, para cualquiera de las condiciones de producción, y el tiempo final de
la línea recta con raíz cuadrada de la gráfica del tiempo,𝑡𝑓𝑟𝑐 . b) Mediante las
expresiones de interpretación descritas en la Tabla 3-5, se calcula el producto √𝑘𝐴𝑐 ,
el volumen de poro,𝑉𝑝 y el volumen original de gas, Además, para el caos de flujo
lineal con no intersección en cero de la recta, se evalúa el daño de la formación,𝑠1
y/o 𝑏.
Tabla 3-13 Ecuación de interpretación para régimen de flujo lineal temprano en un modelo
homogéneo para las condiciones de producción a presión de fondo fluyendo constante,𝑝𝑤𝑓
y gasto de flujo constante de gas,𝑞𝑔 .
56
3.5. Construcción de la gráfica especializada de superposición de tiempo.
Superposición del tiempo para régimen de flujo bilineal con o sin intersección
en cero de la línea recta. a) Se construye la gráfica especializada de
[𝑚(𝑝𝑖 )−𝑚(𝑝𝑤𝑓 )] (𝑞𝑔𝑗 −𝑞𝑔𝑗−1 ) 4
𝑣𝑠 ∑𝑛𝑗=1 √𝑡𝑛 − 𝑡𝑗−1 y se detecta el comportamiento de la
𝑞𝑔 𝑞𝑔𝑛
línea recta, con su intersección en cero y no intersección en cero, como se muestra
en los esbozos de la Fig. [Link]és, se estima el valor de la pendiente,𝑚, y la
no intersección en cero de la línea recta,𝑏. b) Se evalua el termino √𝑘𝑓 𝑤 y el factor
57
de daño de la formación,𝑠1 ,utilizando las ecuaciones de interpretación de la Tabla
3-7.
58
3.5.3. Superposición de tiempo para régimen de flujo esférico con no
intersección en cero de la línea recta.
[𝑚(𝑝𝑖 )−𝑚(𝑝𝑤𝑓 )] (𝑞𝑔𝑗 −𝑞𝑔𝑗−1 ) 1
a) Se elabora la gráfica especializada de 𝑣𝑠 ∑𝑛𝑗=1 y
𝑞𝑔 𝑞𝑔𝑛 √𝑡𝑛 −𝑡𝑗−1
como se muestra en las Fig. 3-19. Después, se calcula la pendiente negativa y la
no intersección en cero de la línea recta. b) Se calcula la permeabilidad,𝑘, y el radio
equivalente de la esfera,𝑟𝑒𝑠𝑓 , utilizando las ecuaciones de interpretación de la Tabla
3-10.
59
Después, se calcula la pendiente de la línea recta y la no intersección en cero de
esta. b) Se calcula el volumen del yacimiento,𝑉𝑌, el factor de forma de Dietz´s,𝐶𝐴 ,
el área de drene,𝐴, y el volumen original de gas,𝐺, mediante las ecuaciónes de
interpretación presentadas en la Tabla 3-11.
60
4. CAPITULO IV MÉTODOS DE ANÁLISIS DE DATOS DE PRESIÓN Y
PRODUCCIÓN
4.1. Introducción.
El análisis de datos de presión y producción de yacimientos de gas y aceite de lutitas
ha evolucionado a nuevas necesidades como: geometrías de pozo complejas,
fracturamiento hidráulico, flujo transitorio durante la vida productiva de los pozos y
a un sin número de efectos no considerados en el análisis tradicional. Actualmente
se tienen tres técnicas de análisis: modelos empíricos, modelos semi- analíticos y
simulación numérica. Los métodos semi- analíticos son simples de usar, pero muy
a menudo dan un resultado aproximado debido a incertidumbre de propiedades del
yacimiento y la necesidad de datos de presión de confianza. Los métodos
numéricos, los cuales están basados en obtener una buena armonización de la
historia, a menudo producen una solución no única. Sin embargo, los métodos
empíricos debido a su simplicidad, sigue siendo el método más utilizado en la
industria.
Desde 2008, un nuevo número de métodos empíricos se han introducido en la
industria petrolera, específicamente para yacimientos de gas en yacimientos no
convencionales. Algunos de ellos son: declinación exponencial extendida (Stretched
exponential decline), declinación exponencial de ley de potencias (power law
exponential decline), modelo de Duong, modelo de crecimiento logístico (logistic
growth model), modelo para pozos verticales fracturados de Arévalo, modelo para
pozos horizontales multi fracturados de Bello. En este capítulo analizaremos los
métodos ya antes mencionados.
4.2. Modelo de Duong.
El modelo de Duong está basado en un flujo lineal a largo plazo. En un yacimiento
hidráulicamente fracturado, el flujo de fluidos se comporta como un régimen lineal
de larga duración, por lo tanto, la producción de gas va a ser:
𝑞 = 𝑞1 𝑡 −𝑛 (4.32)
Donde 𝑛 = 1⁄2 para flujo lineal, 𝑛 = 1⁄4 para flujo bilineal y 𝑞1 para la tasa de
producción de un día.
El gas acumulado se comporta así:
𝑡
(4.33)
𝐺𝑝 = ∫ 𝑞 𝑑𝑡
0
61
Lo que da:
𝑞 (1 − 𝑛) (4.35)
=
𝐺𝑝 𝑡
𝒒
Ilustración 4-1 Grafica 𝑮 𝒗𝒔 𝒕 para determinar los parámetros "a" y "m".
𝒑
(Duong, 2010)
En una gráfica log -log de 𝑞 ⁄𝐺𝑝 𝑣𝑠 𝑡, podemos encontrar representada una línea
recta con pendiente negativa (−𝑚) y una intercepción de 𝑎. La producción futura y
la producción acumulada puede ser predecida usando las siguientes dos relaciones:
𝑎 1−𝑚 −1)
𝑞 = 𝑞1 𝑡 −𝑚 𝑒 1−𝑚(𝑡 (4.37)
𝑞1 𝑎 (𝑡 1−𝑚−1) (4.38)
𝐺𝑝 = 𝑒 1−𝑚
𝑎
62
Ilustración 4-2 Cuatro pasos propuestos para calcular el ajuste de producción.
(Duong, 2010)
Procedimiento
Paso 1. Revisión de datos y corrección.
Se analizan los datos de la historia de producción del pozo realizando una gráfica
presión en la cabeza del pozo, flujo de gas y flujo de agua. Las correcciones del
gasto de gas se hacen basados en la presión de operación promedio.
Paso 2. Determinación de “a” y “m”.
𝑞
Se realiza una gráfica 𝑙𝑜𝑔 − 𝑙𝑜𝑔 𝑑𝑒 𝑣𝑠 𝑡𝑖𝑒𝑚𝑝𝑜mo se muestra en la Fig. 4-2. Se
𝐺𝑝
analiza la gráfica y se determinación la sección en donde se obtendrán los valores
representativos de “a” y “m” para el pozo.
Paso 3. Pronóstico de producción y estimación de reservas.
Determinación de q1. Para obtener q1 es necesario graficar gasto de gas contra t
(a, m). De la Ec. 4.8:
𝑎 1−𝑚 −1)
𝑞 = 𝑞1 𝑡(𝑎, 𝑚) = 𝑞1 𝑡 −𝑚 𝑒 1−𝑚(𝑡 (4.8)
La grafica 𝑞𝑔 vs t (a, m) debe darnos una línea recta con la pendiente de q1. Y de
la Ec. 4.9 obtenemos la producción de gas acumulado:
𝑞1 𝑎 (𝑡 1−𝑚−1) (4.9)
𝐺𝑝 = 𝑒 1−𝑚
𝑎
63
Ventajas
• Es un método sencillo de usar.
• Se obtienen buenos pronósticos de producción tanto en yacimiento de baja
permeabilidad como los de lutitas.
Desventajas
El modelo de Duong presenta problemas cuando existen alguno de estos casos:
• Historias de producción mayores a 2 años.
• Datos de producción con variaciones grandes.
• Mal ajuste de los datos en pozos que tienen un régimen dominado por la
frontera externa.
𝐺𝑝 1 −1 1 𝑞 (4.13)
𝑟𝑝 = 1 − = (Γ [ ]) × Γ [ , − ln ( )]
𝐸𝑈𝑅 𝑛 𝑛 𝑞0
64
Ritmo de producción adimensional:
𝑞 𝑡 𝑛 (4.14)
𝑞𝐷 = = exp [− ( ) ]
𝑞𝑡 𝜏
Producción acumulada adimensional
𝑄 𝜏 1 1 𝑡 𝑛 (4.15)
𝑄𝐷 = = {Γ [ ] × Γ [ , ( ) ]}
𝑞𝑡 𝑛 𝑛 𝑛 𝜏
𝑞0 𝜏 1
Donde 𝐸𝑈𝑅 = Γ [𝑛] es definida como la producción acumulada cuando el ritmo
𝑛
de producción se acerca a 0.
La ventaja más importante es la habilidad de calcular el potencial de recuperación.
Conociendo el ritmo de producción adimensional y el parámetro “𝑛”, podremos
fácilmente calcular el potencial de recuperación, sin conocer el parámetro 𝜏. (Valkó
, 2009) Para analizar los datos de producción se sugiere el siguiente procedimiento:
1. Preparar la serie de datos de 𝑞𝐷 y 𝑄.
2. Asumir un parámetro “𝑛”, calcular 𝑟𝑝 de la Ec. 4.13.
3. Graficar 𝑟𝑝 vs 𝑄. De la grafica debe aparecer una línea recta con dos
intercepciones:
• Intercepción en “Y” =1
• Intercepción en “X” =EUR
Si la intercepción en “Y” no es igual a 1, el parámetro “n” debe ser ajustado.
65
ajusta perfecto a los datos teóricos y tendrán una intercepción e igual a 1 y una
𝝉 𝟏
intercepción x igual a (𝒏) 𝜞 [𝒏]. La grafica también muestran un 𝒓𝒑 de los mismos
datos, si erróneamente asumimos un n=0.4(triángulos azules) o n=0.5 (triángulos
verdes). Se observa claramente que existe un error en el valor de n. (Valko & lee, 2010)
Ilustración 4-4 Caso de simulador(k=0.1 mD, A=200 acres, Fcd=200,xf=300 ft, h=40)
se muestra un EUR 8% menor , pero un desajuste en el perfil de producción, si
usamos el potencial de recuperación. (YU & Miocevic, An Improved Method to Obtain
Reliable Production and EUR Prediction for Wells with Short Production History in
Tight/Shale Reservoirs., 2013)
66
Sin embargo, si no usamos el primer año de producción y la metodología de
potencial de recuperación, en lugar definimos “n” & “𝜏” manualmente igualando
únicamente el 2𝑛𝑑 o/y 3𝑛𝑑 año de producción, podremos encontrar que toda la vida
de producción y el pronóstico de EUR final se igualan bastante bien con las Ecs.9 y
10, como se muestra en la Fig. 4-5.
Ilustración 4-5 El mismo caso de simulación que la Fig. 4-4, el cual muestra un
excelente ajuste en ambos perfiles de producción, si únicamente usamos 𝟐𝒏𝒅 y 𝟑𝒏𝒅
año de los datos de producción. (YU & Miocevic, An improved Method to Obtain Reliable
Production and EUR Prediction for Wells with Short Production History in Tight/shale
Resevoirs., 2013)
Pero existe un problema, en la vida real es difícil obtener un único valor de “n” & “𝜏”
por ajuste manual del 2𝑛𝑑 ó 3𝑛𝑑 año de producción debido al desconocido EUR. El
método de YM-SEPD resuelve ese problema.
De la Ec. 4.11, es fácil obtener lo siguiente:
𝑑𝑞 𝑡 𝑛𝑞 (4.16)
= −𝑛 ( )
𝑑𝑡 𝜏 𝑡
67
𝑞0 (4.18)
𝐿𝑛 [ ] = 𝜏 −𝑛 𝑡 𝑛
𝑞(𝑡)
Ilustración 4-6 Grafica Especializada YM-SEPD. (YU & Miocevic, An improved Method
to Obtain Reliable Production and EUR Prediction for Wells with Short Production History
in Tight/shale Resevoirs., 2013)
Ventajas
• Es un método fácil de aplicar que nos arroja pronóstico de producción
bastante buenos.
68
• El método de YM-SEPD es aplicable a para el pronóstico de producción y
EUR de yacimientos de gas seco, gas húmedo, gas retrogrado y aceites
negros, mientras su perfil de producción tenga la forma de “palo de hockey”.
• Se puede predecir razonables EUR de pozos que tienen 2 o 3 años de
historia de producción independientemente del porcentaje de EUR
remanente (Yu, Lee, Miocevic, Li, & Harris, 2013).
• A diferencia del método de Duong, el método de YM-SEPD no es sensible a
yacimientos con flujo dominado por la frontera.
Desventajas
• Con datos de producciones de 1 a 2 años de producción se obtienen errores
en el pronóstico de la producción de gas, pero se pueden calcular una pseudo
historia con el método de Duong para reducir el porcentaje de error (Yu, Lee,
Miocevic, Li, & Harris, 2013).
𝐷 = 𝐷∞ + 𝐷1 𝑡 −(1−𝑛) (4.23)
El parámetro b no es constante y está en función del tiempo. Sustituyendo la Ec.
4.23 en la Ec. 4.22, obtenemos:
̂𝑖 (1 − 𝑛)
𝑛𝐷 (4.24)
𝑏= 2 𝑡 −𝑛
̂𝑖 + 𝐷∞
[𝑛𝐷 𝑡 (1−𝑛) ]
Sustituyendo la Ec.4.23 en la Ec. 4.21 y resolviendo para el gasto(𝑞), obtenemos
la relación de gasto del “método de declinación exponencial de potencias”, la cual
es:
69
𝐷1 𝑛 (4.25)
𝑞 = 𝑞̂𝑖 exp [−𝐷∞ 𝑡 − 𝑡 ]
𝑛
La cual se reduce en la forma de relación de la tasa de declinación de ley de
potencias con “los-ratio” definida por IIk (IIk, Rushing, A., & Blasingame, 2008a).
̂𝑖 𝑡 𝑛 ]
𝑞 = 𝑞̂𝑖 exp[−𝐷∞ 𝑡 − 𝐷 (4.26)
Donde:
𝑞̂𝑖 = 𝐺𝑎𝑠𝑡𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑐𝑖𝑜𝑛 𝑑𝑒𝑓𝑖𝑛𝑖𝑑𝑎 𝑒𝑛 𝑙𝑎 𝐸𝑐. 4.25.
𝐷1 = 𝐼𝑛𝑡𝑒𝑟𝑝𝑒𝑐𝑖𝑜𝑛 𝑒𝑛 𝑑𝑒𝑐𝑙𝑖𝑛𝑎𝑐𝑖𝑜𝑛 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑎 1 𝑢𝑛𝑖𝑑𝑎𝑑 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜 𝑑𝑒𝑓𝑖𝑛𝑖𝑑𝑎 𝑒𝑛 𝑙𝑎 𝐸𝑐. 4.23.
70
El modelo de incremento logístico ha sido desarrollado para el pronóstico de
producción y la estimación de reservas en yacimientos de muy baja permeabilidad.
El modelo fue originalmente desarrollado por Tsoularis y Wallace en 2001 y es del
mismo tipo que el modelo de Blumberg. Las curvas de crecimiento de Blumberg son
hiperbólicas y es de esta forma en que los pozos no convencionales(tight) de aceite
y gas declinan.
En un principio el modelo se presentó de la siguiente forma:
𝐾(𝑡 + 𝑎)𝑛 (4.27)
𝑁(𝑡) =
𝑏 + (𝑡 + 𝑎)𝑛
Donde:
𝑁: 𝑃𝑜𝑏𝑙𝑎𝑐𝑖𝑜𝑛
𝐾: 𝐶𝑎𝑝𝑎𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑐𝑎𝑟𝑔𝑎
𝑎: 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝑏: 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝑛: 𝑃𝑎𝑟𝑎𝑚𝑒𝑡𝑟𝑜 𝐸𝑥𝑝𝑜𝑛𝑒𝑛𝑐𝑖𝑎𝑙
𝑡: 𝑇𝑖𝑒𝑚𝑝𝑜
El modelo de incremento logístico que se usa para el pronóstico de producción de
aceite y gas en yacimientos no convencionales(tight) es el siguiente:
𝐾𝑡 𝑛 (4.28)
𝑄(𝑡) =
𝑎 + 𝑡𝑛
Donde:
𝑄: 𝑃𝑟𝑜𝑑𝑢𝑐𝑐𝑖𝑜𝑛 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑑𝑎
𝐾: 𝐶𝑎𝑝𝑎𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑐𝑎𝑟𝑔𝑎
𝑛: 𝐸𝑥𝑝𝑜𝑛𝑒𝑛𝑡𝑒 ℎ𝑖𝑝𝑒𝑟𝑏ó𝑙𝑖𝑐𝑎
𝑎: 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
La capacidad de carga,𝑘, es prácticamente el aceite y gas recuperable del
yacimiento del mecanismo de producción. La capacidad de carga no toma en cuenta
las limitaciones económicas que podrían ser utilizadas para pronosticar la EUR. La
constante 𝑎 es el tiempo elevado a la potencia 𝑛 en el cual se ha recuperado la
mitad del aceite o gas existente.
71
El ritmo de producción se obtiene al derivar la Ec. 4.28, quedando como resultado
lo siguiente:
𝑑𝑄 𝐾𝑛𝑏𝑡 𝑛−1 (4.29)
𝑞(𝑡) = =
𝑑𝑡 (𝑎 + 𝑡 𝑛 )2
Donde:
𝑞: 𝑅𝑖𝑡𝑚𝑜 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖ó𝑛
El modelo de incremento logístico coincide con los datos de producción para los
pozos de aceite y gas no convencionales.
4.7. Modelo para pozos Horizontales multifracturados de Bello.
El modelo de Bello y cols (Wattenbarger R. O., 2008) consideran un pozo horizontal
con multifracturas, donde identificaron cinco regiones de flujo transitorio, a través de
modelos lineales de doble porosidad, para los cuales presento ecuaciones. En la
Fig.--- se muestra el modelo conceptual de los bloques de la matriz de la roca.
72
Ilustración 4-8 Grafica de las cinco regiones de flujo del moldeo de Bello.
(Bello, 2009)
En la región 1 identificamos un flujo lineal transitorio temprano, únicamente en el
sistema de fracturas. El drene de la matriz es despreciable:
1262 𝑇 1 (4.30)
√𝑘𝑓 𝐴𝑐𝑤 =
√𝜔(𝜙𝜇𝑐𝑡 )𝑓+𝑚 𝑚1
Donde 𝑚1 es la pendiente obtenida de la gráfica especializada ∆𝑚(𝑝)⁄𝑞𝑔 𝑣𝑠 √𝑡.
73
Donde 𝑚3 es la pendiente de una grafica de ∆𝑚(𝑝)⁄𝑞𝑔 𝑣𝑠 √𝑡.
74
5. CAPÍTULO V COMPARACION DE LOS METODOS EN EL ANALISIS DE
LOS CASOS DE CAMPO
5.1. Introducción
En este capítulo se muestra la aplicación de métodos semi- analíticos y métodos
empíricos para el análisis de la presión- producción de seis yacimientos de gas de
baja permeabilidad con diferentes características. El procedimiento que se presenta
es el siguiente: primero, se realizó un ajuste de los datos de presión mediante el
método de media móvil, para suavizar el histórico de producción. Después, se
identificaron los regímenes de flujo presentes en el yacimiento mediante graficas de
diagnóstico de pseudo- presiones normalizadas para así poder realizar un análisis
individual de cada régimen de flujo a través de graficas especializadas con
superposición de tiempo y expresiones matemáticas que se analizaron en el
capítulo 3. En el caso de los yacimientos de gas de lutitas se utilizó el método de
bello para identificar la permeabilidad y volumen original de gas. Posteriormente, se
realizó el análisis de la declinación de producción del yacimiento por medio de
curvas de declinación para poder obtener el mejor ajuste y un pronóstico de
producción.
En seguida se presenta un inventario y datos generales de los pozos analizados en
este trabajo. Después, se presenta el análisis de los primeros tres pozos, los cuales
son del tipo del tipo tight. Luego se muestra el procedimiento de los pozos de gas
de lutitas.
Finalmente, se muestra un resumen de este capítulo V.
75
5.2. Casos de Campo Tight Gas
76
Tabla 5-1 Datos general Pozo CIGE 231
Datos Generales
Fecha 27-ene-01
Campo Castlegate
Pozo CIGE 231
Tipo de gas Seco
Datos del Yacimiento
Presión inicial, pi 6700 psia
Presión de fondo Fluyendo, pwf 500.00 psia
Temperatura Inicial, Ti 220 oF
77
Ilustración 5-1 Histórico de Producción del pozo CIGE231
78
Ilustración 5-3 Grafica de diagnóstico de la pseudo presión normalizada contra
tiempo,∆𝒎(𝒑)⁄𝒒𝒈 𝒗𝒔 𝒕, del pozo CIGE 231.
Una vez detectado los regímenes de flujo se realizaron las gráficas especializadas
de cada régimen. La Fig. 5-4 es la gráfica especializada para flujo radial, como se
puede observar se tiene una recta que termina a los 12 días de producción. La recta
𝐷
tiene una pendiente 𝑚𝐶𝑅𝑅 = 695259 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] y una intercepción al eje de
𝐷
las ordenadas de 𝑏𝐶𝑅𝑅 = 320,000 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝].Mediante las ecuaciones del
Capítulo 3 para régimen radial se obtuvieron algunas propiedades del yacimiento
como: 𝑘 = 0.0139101𝑚𝑑 y 𝑠 = −4.74543.
La Fig. 5-5 muestra una gráfica especializada con superposición de tiempo para un
∆𝑚(𝑝)⁄
régimen lineal. Los datos de la caída de la pseudo- presión, 𝑞𝑔 , siguen una
tendencia de línea recta cuando son graficados contra la superposición del tiempo
de √𝑡.Lo que confirma que existe un régimen de flujo lineal tardío hasta los 165 días
𝐷 .5
de producción. La pendiente de la recta es de 𝑚𝐶𝑅𝐿 = 73,579 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] con
𝐷
una ordenada de 𝑏𝐶𝑅𝐿 = 300,000 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝]. Se obtuvo un valor de área de
drene de 𝐴 = 28.643[𝐴𝑐𝑟𝑒𝑠], distancia de la fractura hidráulica a la frontera de 𝑌𝑒 =
6684.17[𝑓𝑡] y un valor mínimo de volumen original de gas,𝐺 = 0.83315𝐵𝑠𝑐𝑓.
Como se puede observar en la Fig. 5-6, solo al inicio de la producción se tiene una
tendencia de lineal, lo que significa un flujo bilineal temprano. La pendiente que se
3
2 𝐷4
obtuvo de la gráfica especializada fue de 𝑚𝐶𝑅𝐵 = 444,438 [𝑝𝑠𝑖𝑎 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] y una
79
𝐷
ordenada de 𝑏𝐶𝑅𝐿 = 300,000 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝].Mediante las Ecuaciones del
Capítulo 3 se obtuvieron los siguientes parámetros del yacimiento: Factor de daño
de la formación 𝑆1 = 8.0366 y la raíz de la permeabilidad de mayor flujo por el ancho
de la fractura de √𝑘1 𝑤 = 2266.91.
80
Ilustración 5-6 Grafica especializada de pseudo presión normalizada contra
superposición de tiempo para flujo lineal del pozo CIGE 231.
Por último, se realizó la gráfica de diagnóstico del flujo dominado por la frontera,
Fig. 5-7, donde se obtuvo la tendencia de línea recta en los últimos días de
producción hasta los 152 días. Se obtuvo una pendiente 𝑚𝑆𝑆𝑃 = 1272.48 [𝑝𝑠𝑖𝑎2 −
𝐷 𝐷
− 𝑐𝑝] y una ordenada 𝑏𝑆𝑆𝑃 = 300,000 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝]. Se obtuvo un
𝑀𝑠𝑐𝑓
volumen poroso 𝑉𝑝 = 2.8[𝑀𝑀𝑟𝑐𝑓] y volumen original de gas de 𝐺 = 0.1462 𝐵𝑠𝑐𝑓.
81
Se procedió a analizar la declinación de pozo de gas mediante métodos empíricos
y ajustar su producción futura. En primer lugar, se analizó el pozo con el método de
Duong, se introdujeron los datos de la historia de producción en el software de
Matlab y se calcularon los parámetros “a” y “m” mediante una gráfica 𝑙𝑜𝑔 −
𝑞𝑔
log 𝑑𝑒 𝐺𝑝 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜g. 5.8. Donde obtuvimos unos valores de “a=-1.22” y
“m=1.924”. Se obtuvo un gasto inicial de 𝑞1 =1.19477 y un error de ajuste, Fig. 5-9
y 5.10, de producción de gas de 0.1822, esto se atribuye a que se cuenta con tan
solo 256 días de historia de producción y para un ajuste más confiable es necesario
mínimo 2 años de historia de producción.
82
Ilustración 5-9 Ajuste de producción de gas,𝒒𝒈, mediante el método de Duong del
pozo CIGE 231.
83
Para el análisis del método de declinación exponencial extendida de Valkó se
obtuvo:
𝑀𝑠𝑐𝑓
𝑞𝑖 = 3649.6
𝐷
Sustituyendo las relaciones entre producciones acumuladas para los años tres y
dos con respecto a la del uno, se desarrolló el sistema de ecuaciones siguiente.
1 1 23.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
= 𝑟21
1 1 11.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
1 1 35.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
= 𝑟31
1 1 11.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
84
Ilustración 5-11 Ajuste de producción de gas,𝒒𝒈, mediante el metodo de Decliacion
exponencial extendida del pozo CIGE 231.
85
El siguiente método que se analizo fue el método modificado de declinación
exponencial extendida de Yu y Miocevic (YU & Miocevic, An Improved Method to
Obtain Reliable Production and EUR Prediction for Wells with Short Production
History in Tight/Shale Reservoirs., 2013), que a diferencia del método de Valko es
un método más sencillo de usar. Primero se hizo una gráfica log-log log(qo/q(t)) vs
tiempo del cual se seleccionó el intervalo de 40 a 130 días, donde podemos
encontrar una tendencia lineal. De la Fig. 5-13 obtuvimos “n= 0.413771” de la
pendiente y 𝜏 se calculó con la Ec.5.1 y se obtuvo “𝜏 = 73.6025”, los cuales se
utilizaron para el cálculo del pronóstico de gas y producción de gas con las Ec. 5.2:
𝐿𝑛(𝐼𝑛𝑡) (5.1)
𝜏 = 𝐸𝑥𝑝 [− ]
𝑛
𝑡 𝑛 (5.2)
𝑞 = 𝑞0 × 𝐸𝑥𝑝 [− ( ) ]
𝜏
Ilustración 5-13Grafica especializada del método YM-SEPD para el pozo CIGE 231.
86
Ilustración 5-14 Ajuste de la producción de gas,𝒒𝒈, mediante el método YM-SEPD del
pozo CIGE 231.
87
El siguiente método por el cual se analizó el pozo CIGE 231 fue el método de
potencias presentado por IIk (IIk, Rushing, A., & Blasingame, 2008a) (IIK, Perego,
Rushing, & Blasingame, 2008b). Se calcularon los parámetros “D” y “b” mediante el
𝑑𝑞
método diferencia numérica de Bourdet. Después, se calculo mediante el método
𝑑𝑡
1 𝑞
de diferencia numérica de Bourdet y se grafico 𝐷 = − 𝑑𝑞 vs tiempo, Fig. 5-16, para
⁄𝑑𝑡
𝐷1
obtener los parámetros 𝑛 = 1 − 𝑚, 𝐷𝑖 = ,𝐷∞ = 0 y 𝑞𝑖 = 5.2. Por último, se calculó
𝑛
el pronóstico de producción de gas con la Ec. 5.4, obteniendo un Error del
0.0577903 como se puede observar en el ajuste de la Fig. 5-17.
1 𝑞 (5.3)
=−
𝐷 𝑑𝑞⁄
𝑑𝑡
̂𝑖 𝑡 𝑛 ]
𝑞 = 𝑞̂𝑖 exp[−𝐷∞ 𝑡 − 𝐷 (5.4)
Ilustración 5-16 Grafica especializada del método de Potencias para el pozo CIGE 231.
88
Ilustración 5-17 Ajuste de la producción de gas,𝒒𝒈, mediante el método de Potencias
del pozo CIGE 231.
89
Ilustración 5-18 Ajuste de la producción de gas,𝒒𝒈, mediante el método de
Crecimiento Logístico del pozo CIGE 231.
90
Pozo CIGE 218. En este caso se analizó un pozo seco tight que ha sido producido
por menos de un año,294 días. Como se puede observar en la Fig. 5-20 existen dos
fluctuaciones a los 270 días de producción, por lo que se suavizo la historia de
producción mediante el método de media móvil.
La Fig. 5-21 y 5-22 muestran las gráficas de diagnóstico doble logarítmica de
producción acumulada contra tiempo, 𝐺𝑝 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜, y la caída de pseudo presión
∆𝑚(𝑝)⁄
sobre gasto de gas contra tiempo, 𝑞𝑔 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜. Los datos de producción
en estas figuras muestran un flujo lineal temprano definido por una pendiente de un
medio, seguido de un comportamiento de pendiente de la recta de un cuarto, lo que
nos indica un régimen de flujo bilineal el cual termina a los 230 días. A los 250 días
se encontró otro comportamiento de la pendiente de un medio (régimen lineal tardío)
que termina a los 280 días.
Tabla 5-2 Datos generales pozo CIGE 218
Datos Generales
Fecha 27-ene-01
Campo Castlegate
Pozo CIGE 218
Tipo de gas Seco
Datos del Yacimiento
Presión inicial, pi 6700 psia
Presión de fondo Fluyendo, pwf 200.00 psia
Temperatura Inicial, Ti 220 oF
91
Ilustración 5-20 Histórico de presión- producción del pozo CIGE 218.
92
Ilustración 5-22 Grafica de diagnóstico de pseudo presión normalizada contra
tiempo,∆𝒎(𝒑)⁄𝒒𝒈 , del pozo CIGE 218.
Una vez definido el régimen de flujo se realizaron las gráficas especializadas para
los diferentes regímenes de flujo. En la Fig. 5-23 se muestra la gráfica especializada
para régimen radial donde observamos una recta lineal con pendiente 𝑚𝐶𝑅𝑅 =
𝐷
695259 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝].Una intercesión en cero y termino del régimen de flujo
radial a los 142 días. Mediante las ecuaciones del capítulo 3 se obtuvo una
permeabilidad, longitud media de la fractura y daño del pozo de: 𝑘 =
0.0220744 𝑚𝑑,𝑥𝑓 = 419.185[𝑓𝑡] 𝑠 = −4.93879.
93
Ilustración 5-23 Grafica especializada de pseudo presión normalizada contra
superposición de tiempo para flujo radial del pozo CIGE 218.
94
Ilustración 5-24 Grafica especializada de pseudo- presión de tiempo normalizada
contra superposición de tiempo para flujo Lineal del pozo CIGE 218.
95
Ilustración 5-25 Grafica especializada de pseudo- presión normalizada contra
superposición de tiempo para flujo Bilineal del pozo CIGE 218.
96
Ilustración 5-26 Grafica de ajuste de producción de gas mediante el método de Duong
del pozo CIGE 218.
97
Ilustración 5-28 Ajuste de la producción de gas mediante el método de SEPD del pozo
CIGE 218.
98
La Fig. 5-30 y 5-31 muestran el ajuste de la producción de gas y producción
acumulada, respectivamente, mediante el método modificado de Declinación
Exponencial Extendida. Mediante una gráfica log-log de la relación de
qo
log (q(t)) vs tiempo se obtuvieron los siguientes valores
n 0.2399
Int 0.3319
𝜏 99.0533
99
Ilustración 5-31 Ajuste de la producción acumulada de gas mediante el método de
YM-SEPD del pozo CIGE 218.
m -0.77846
D1 0.07805
n 0.221533
Di 0.3523
𝒒𝒊 2.3
100
Ilustración 5-32 Ajuste de la producción de gas mediante el método de Potencias del
pozo CIGE 218.
101
Ilustración 5-33 Ajuste de la producción de gas mediante el método de Crecimiento
logístico del pozo CIGE 218.
102
5.3. Casos de Campo Shale Gas
Pozo Shale A. En este caso, se describe paso a paso la metodología de análisis
del pozo Shale A. El pozo ha producido por 442 días, Fig. 5-35, pero existen
variaciones en la producción, se infiere que se hizo una estimulación a los 288 días.
Es por ello por lo que para el análisis del pozo se hizo para dos intervalos. El primer
intervalo fue del día 1 a los 222 días de los datos de producción de gas del pozo, el
segundo intervalo fue de 288 a 442 días. Para cada caso se consideró como un
pozo independiente, pero tomando los mismos datos del yacimiento y fluido de la
Tabla 5-5.
Primero se realizó el análisis de los regímenes de flujo que se encuentran en el
yacimiento mediante graficas de diagnóstico, utilizando el primero intervalo de los
datos del pozo. La Fig. 5-36 y 5-37 muestran las gráficas de diagnóstico doble
logarítmica de producción acumulada contra tiempo, 𝐺𝑝 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜, y la caída de
∆𝑚(𝑝)⁄
pseudo presión sobre gasto de gas contra tiempo, 𝑞𝑔 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜, para el
primer intervalo del pozo sale A. Los datos de producción en estas figuras muestran
un flujo Bilineal, la cual se identifica por una pendiente de un cuarto y termina a los
50 días de producción. Después, se presentó una tendencia lineal de un medio, Fig.
5-36, lo que significa un régimen de flujo lineal tardío.
Las Fig. 5-38 y 5-39, muestran las gráficas de diagnóstico de producción acumulada
contra tiempo, 𝐺𝑝 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜, y la caída de pseudo presión sobre gasto de gas contra
∆𝑚(𝑝)⁄
tiempo, 𝑞𝑔 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜, para el segundo intervalo del pozo Shale A. Como se
puede observar en la Fig. 5-38 se encontró una tendencia de un cuarto de
pendiente, lo que significa un régimen de flujo Bilineal, el cual se presenta duran
toda la segunda parte del pozo productor Shale A.
Una vez identificado los regímenes de flujo que dominan, se prosiguió a hacer el
análisis de cada régimen de flujo mediante graficas especializadas de la primera
parte del histórico de producción. La Fig. 5-40 muestra la gráfica especializado
régimen radial donde observamos una recta lineal con pendiente 𝑚𝐶𝑅𝑅 =
𝐷
345686 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝].Una intercesión en cero y termino del régimen de flujo
radial a los 57 días. Mediante las ecuaciones del capítulo 3 se obtuvo una
permeabilidad longitud media de la fractura y daño del pozo de: 𝑘 = 0.011[𝑚𝑑]𝑥𝑓 =
443.5 [𝑓𝑡] 𝑠 = −4.32.
103
Tabla 5-5 Datos general Pozo Shale A
Datos Generales
Fecha
Campo
Pozo Shale A
Tipo de gas Seco
Datos del Yacimiento
Presión inicial, pi 5100 psia
Presión de fondo Fluyendo, pwf 500.00 psia
Temperatura Inicial, Ti 207 oF
104
Ilustración 5-35 Histórico de presión- producción del pozo Shale A. Se puede
observar que existe bastante fluctuaciones de la presión- producción, es por ello por
lo que se optó por analizar en dos secciones la producción del pozo Shale A.
105
Ilustración 5-36 Grafica de diagnóstico de flujo de producción acumulada contra
tiempo,𝑮𝒑 𝒗𝒔 𝒕, de la primera sección del histórico de producción del pozo Shale A.
106
Ilustración 5-38 Grafica de Diagnostico de flujo de producción acumulada de gas
contra tiempo,𝑮𝒑 𝒗𝒔 𝒕, de la segunda sección del histórico de producción del pozo
Shale A.
107
Ilustración 5-40 Grafica especializada de pseudo- presión normalizada contra
superposición de tiempo para flujo radial de la primera sección del histórico de
producción del pozo Shale A.
108
Ilustración 5-42 Grafica especializada de pseudo- presión normalizada contra
superposición de tiempo para flujo Bilineal de la primera sección del pozo Shale A.
en cero.
Se procedió a hacer el análisis de la segunda parte de los datos del pozo mediante
graficas especializadas. La Fig. 5-43 muestra en análisis del régimen de flujo radial
con una gráfica especializada de superposición de tiempo contra pseudo- presión
normalizada. Debido a las grandes variaciones de gasto de gas se tiene una
distorsión en la gráfica especializada, debido a que muestra dos tendencias. Solo
se tomó en cuenta la primera tendencia en todas las gráficas especializadas. En la
Fig. 5-43 podemos observar que la primera tendencia tiene una pendiente de
𝐷
𝑚𝐶𝑅𝑅 = 767289 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] con una intersección en cero.
109
Mediante las ecuaciones vistas en el Capítulo 3 se calcularon los parámetros de la
Tabla 5-12 a 5-14. La Fig. 5-44 muestra la gráfica especializada del régimen
Bilineal, solo se tomó en cuenta la primera tendencia. Se obtuvo una pendiente de
𝐷 .5
𝑚𝐶𝑅𝐵 = 673,570 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] con una intersección en𝑏𝐶𝑅𝐿 = 8𝑥10ˆ5 [𝑝𝑠𝑖𝑎2 −
𝐷 .5
− 𝑐𝑝].
𝑀𝑠𝑐𝑓
110
Ilustración 5-44 Grafica especializada de pseudo- presión normalizada contra
superposición de tiempo para flujo Bilineal de la segunda sección del pozo Shale A.
111
Posteriormente, se realizó el análisis de producción de gas de pozo mediante cinco
métodos empíricos de curvas de declinación, únicamente para el primer intervalo
del pozo, debido a que las distorsiones de la segunda parte generan errores muy
grandes. La Fig. 5-45 y Fig. 5-46 muestra el análisis de producción de gas y
producción acumulada del yacimiento Shale A mediante el método de Duong.
Mediante graficas especializadas y un ajuste de regresión potencial se obtuvo un
𝑴𝑴𝒔𝒄𝒇
valor de “a=-1.129”, “m=1.334 ”y 𝒒𝒊 = 𝟐. 𝟕𝟓𝟏𝟔 .Se obtuvo un error de ajuste de
𝑫
0.2577 del análisis de la producción de gas, pero como se muestra en la Fig. 5-45,
se tiene mejor ajuste con la producción acumulada del gas. Esto se debe a que se
cuenta con muy pocos datos de producción de gas.
112
Al final se obtuvo un error de ajuste de 0.137032, el cual fue mejor que el obtenido
con el método de Duong.
El tercer método que se analizo fue el método modificado de declinación
exponencial extendida de Yu y Miocevic (YU & Miocevic, An Improved Method to
Obtain Reliable Production and EUR Prediction for Wells with Short Production
History in Tight/Shale Reservoirs., 2013), que a diferencia del método de Valko es
un método más sencillo de usar. Primero se hizo una gráfica log-log log(qo/q(t)) vs
tiempo del cual se seleccionó el intervalo de 9 a 150 días, dado que fue el rango de
datos donde se encontró una tendencia lineal. Posteriormente, de la gráfica anterior
se obtuvieron los valores de 𝑛 = 0.49455 y 𝜏 = 176.093, los cuales se utilizaron para
el pronóstico de la producción de gas y producción acumulada de gas, como se
muestra en la Fig. 5-49 y Fig. 5-50. Al final se obtuvo un error de ajuste de 0.12225.
Ilustración 5-47 Ajuste de la producción de gas mediante el método de SEPD del pozo
Shale A.
113
Ilustración 5-48 Ajuste de la producción acumulada de gas mediante el método de
SEPD del pozo Shale A.
114
Ilustración 5-50 Ajuste de la producción acumulada de gas mediante el método de
YM-SEPD del pozo Shale A.
115
Ilustración 5-51 Ajuste de la producción de gas mediante el método de Potencias del
pozo Shale A.
K 1.19061
n 0.961445
a 328.068
b 2.62409e+08
116
Ilustración 5-52 Ajuste de la producción de gas mediante el método de Crecimiento
Logístico del pozo Shale A.
117
Pozo Shale C. El tercer pozo que se analizo fue el denominado Shale C, el cual
cuenta con una historia de producción de 1274 días, pero tiene muchas
fluctuaciones en la producción, además se supusieron algunos datos del yacimiento
y fluidos, Tabla 5-7. Se hizo un filtrado de datos mediante el método de medial móvil
para suavizar la producción y obtener mejor resultados al momento de aplicar los
métodos empíricos, Fig. 5-55.
Las Fig. 5-56 muestra la gráfica de diagnóstico de flujo con ejes log-log de la
producción de gas acumulada contra el tiempo de producción,𝐺𝑝 𝑣𝑠 𝑡 para el pozo
Shale C. De igual manera, la Fig. 5-57 muestra una gráfica de diagnóstico de flujo
con ejes log-log de la caída de pseudo- presión y el gasto de flujo de gas contra el
tiempo de producción,∆𝑚(𝑝)⁄𝑞𝑔 𝑣𝑠 𝑡.En la Fig. 5-56 se muestra un flujo bilineal
temprano con pendiente de un cuarto, el cual se comienza a desviar a los 132 días
de producción. A los 132 días se comienza una tendencia de media pendiente, lo
cual confirma un régimen de flujo lineal que termina a los 560 días. En análisis de
flujo muestra el alcance de la primera frontera a los 211 días, Fig. 5-57. Al final se
identifica una tendencia de flujo pseudo- estacionario (PSS), ya que la curva de la
caída de la pseudo- presión, Fig. 5-57, se dobla desde la media pendiente hacia
abajo a los 1000 días.
Tabla 5-7 Datos generales pozo Shale C
Datos Generales
Fecha
Campo
Pozo Shale C
Tipo de gas Seco
Datos del Yacimiento
Presión inicial, pi 5100 psia
Presión de fondo Fluyendo, pwf 500.00 psia
Temperatura Inicial, Ti 210 oF
118
Nitrógeno, N2 0.001 Fracción
mol
Temperatura @ s.c. 60.00 oF
119
Ilustración 5-55 Histórico de producción suavizado mediante el método de media
móvil del pozo Shale C.
120
Ilustración 5-57 Grafica de diagnóstico de flujo de pseudo- presión normalizada
contra tiempo del pozo Shale C.
Para el análisis de los distintos regímenes de flujo del pozo Shale C, se utilizó la
técnica de superposición de tiempo, debido a que los datos de producción de gas
cuentan con muchas variaciones. La Fig.5-58 muestra la gráfica especializada de
121
∆𝑚(𝑝)
𝑣𝑠 𝑠𝑢𝑝𝑒𝑟𝑝𝑜𝑠𝑖𝑐𝑖𝑜𝑛 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜 𝐿𝑜𝑔(𝑡) para el pozo Shale C. Como se muestra
𝑞𝑔
en la Fig. 5-58, los datos de producción siguen una tendencia de línea recta con
𝐷
una pendiente de 𝑚𝐶𝑅𝑅 = 1.53 𝑒 + 06 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] y ordenada a cero, la línea
recta se comienza a desviar a los 145 dí[Link] el análisis de ecuaciones de
interpretación del capítulo 3 para flujo radial bajo la condición de presión de fondo
fluyendo constante,𝑃𝑤𝑓 , se estimaron los siguientes valores: 𝑘 = 0.00239𝑚𝑑 y 𝑠 =
−5.69816.
En la Fig. 5-59 los datos de producción se analizaron mediante la gráfica
∆𝑚(𝑝)
especializada de 𝑣𝑠 𝑠𝑢𝑝𝑒𝑟𝑝𝑜𝑠𝑖𝑐𝑖𝑜𝑛 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜 √𝑡. Se observa un
𝑞𝑔
𝐷 .5
comportamiento en línea recta con pendiente 𝑚𝐶𝑅𝐿 = 70,583.4 [𝑝𝑠𝑖𝑎2 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] y
una intersección a la ordenada en cero. Se utilizaron las ecuaciones de
interpretación, capitulo 3, para flujo línea bajo la condición de presión de fondo
fluyendo constante, 𝑃𝑤𝑓 , para estimar un área de drene de 𝐴 = 17.5186[𝐴𝑐𝑟𝑒𝑠],
distancia de la fractura hidráulica a la frontera de 𝑌𝑒 = 442.735[𝑓𝑡] y un valor mínimo
de volumen original de gas,𝐺 = 21.4364𝐵𝑠𝑐𝑓.
El análisis de los datos de producción de gas mediante la gráfica especializada de
la Fig. 5-60 muestran la tendencia línea, lo que confirma la existencia de un régimen
de flujo bilineal. La línea recta de la Fig. 5-60 tiene una pendiente de 𝑚𝐶𝑅𝐵 =
3
2 𝐷4
70583.4 [𝑝𝑠𝑖𝑎 − 𝑀𝑠𝑐𝑓 − 𝑐𝑝] y una ordenada en cero. La línea recta comienza a
122
realizaron las gráficas especializadas de pseudo presión normalizado contra pseudo
1
tiempo para flujo bilineal, [𝑚(𝑝𝑖 − 𝑚(𝑝𝑤𝑓 )]⁄𝑞𝑔 𝑣𝑠 (𝑡 ∗ ) ⁄4, y flujo
1
lineal, [𝑚(𝑝𝑖 − 𝑚(𝑝𝑤𝑓 )]⁄𝑞𝑔 𝑣𝑠 (𝑡 ∗ ) ⁄[Link] que el primero caso se usaron las
pendientes y se obtuvo una permeabilidad de matriz de 𝑘𝑚 = 0.0015824[𝑚𝑑] y una
permeabilidad de fractura de 𝑘𝑓 = 0.0000089311[𝑚𝑑].
123
Ilustración 5-60 Grafica especializada de pseudo- presión normalizada contra
superposición de tiempo para el flujo Bilineal del pozo Shale C.
124
Primero se analizó el pozo mediante el método de Duong, donde se grafico 𝑙𝑜𝑔 −
𝑞𝑔
𝑙𝑜𝑔 𝑑𝑒 𝐺𝑝 𝑣𝑠 𝑇𝑖𝑒𝑚𝑝𝑜 ediante una regresión potencial se calculo la pendiente “m=-
1.20362” y la ordenada el origen “a=1.70158”. Después, se grafico
𝑞𝑔 𝑣𝑠 𝑡(𝑎, 𝑚),Fig. 5-62, y se obtuvo un gasto inicial de gas de 𝑞𝑖 = 2.35932. Al final
se obtuvo un error de ajuste de producción de gas y gasto acumulado de gas del
0.148938, Fig. 5-64 y Fig. 5-65.
Ilustración 5-63 Grafica especializada de gasto de gas contra t(a.m.),𝒒𝒈 𝒗𝒔 𝒕(𝒂, 𝒎), del
pozo Shale C.
125
Ilustración 5-64 Ajuste de la producción de gas mediante el método de Duong del
pozo Shale C.
126
Para el análisis del método de declinación exponencial extendida de Valkó se
obtuvo:
𝑀𝑠𝑐𝑓
𝑞𝑖 = 20024.7
𝐷
Sustituyendo las relaciones entre producciones acumuladas para los años tres y
dos con respecto a la del uno, se desarrolló el sistema de ecuaciones siguiente.
1 1 23.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
= 𝑟21
1 1 11.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
1 1 35.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
= 𝑟31
1 1 11.5 𝑛
Γ [ 𝑛 ] − Γ [𝑛 , ( 𝜏 ) ]
Ilustración 5-66 Ajuste de la producción de gas mediante el método de SEPD del pozo
Shale C.
127
Ilustración 5-67 Ajuste de la producción acumulada de gas mediante el método SEPD
del pozo Shale C.
128
Ilustración 5-68 Grafica especializada del método de YM-SEPD del pozo Shale C.
129
Ilustración 5-70 Ajuste de la producción acumulada de gas mediante el método de
YM-SEPD del pozo Shale C.
130
Por medio de la Ec.5.8 se calculó el pronóstico de producción de gas, obteniendo
un error de ajuste de 0.114517, Fig. 5-71.
̂𝑖 𝑡 𝑛 ]
𝑞 = 𝑞̂𝑖 exp[−𝐷∞ 𝑡 − 𝐷 (5.8)
131
Tabla 5-8 Resultados modelo Crecimiento Logístico-Shale C
𝒌 1.28871
𝒏 0.817643
𝒂 189.614
𝒃 1.88974e+08
132
Ilustración 5-73 Ajuste de la producción acumulada de gas mediante el método de
Crecimiento Logístico del pozo Shale C.
133
La explicación geológica varia de caso en caso, pero la idea general es que existe
una trayectoria de flujo de alta permeabilidad rodeada por una formación de mucha
más baja permeabilidad. Una explicación es que existe un sistema de fracturas
naturales paralelas. EL flujo de fluido de la formación de baja permeabilidad dentro
de la trayectoria de flujo de la de mayor permeabilidad genera el flujo transitorio
bilineal.
En el caso de los Pozos Castlegate se puede observar que existe un régimen de
flujo lineal temprano, seguido de un flujo bilineal, un flujo lineal tardío y por último
se alcanza el régimen dominado por la frontera. El flujo lineal temprano es afectado
por √𝑘1 𝐴𝑐1 de la trayectoria de flujo de alta permeabilidad y 𝑘1 es la permeabilidad
de la trayectoria de flujo de la misma área. El periodo de flujo lineal tardío es
afectado por √𝑘𝐴𝑐 de la formación, en donde 𝐴𝑐 es el área de flujo en la trayectoria
provista de mayor permeabilidad para el flujo perpendicular proveniente de la
formación hacia el interior de la trayectoria de flujo de alta permeabilidad y 𝑘 es la
permeabilidad de la formación.
En el análisis de flujo transitorio de los pozos Tight 231 y Tight 218 parece indicar
algo sobre la naturaleza de este yacimiento de gas de baja permeabilidad. Como se
muestra en la Tablas 5-12 a 5-14, el yacimiento no se comporta como un yacimiento
homogéneo con valores sencillos de permeabilidad y de espesor. Además, este
comportamiento sugiere que la geometría de flujo de la fractura es más complicada
que el modelo común de un rectángulo con una longitud media de fractura,
permeabilidad y con la altura de la fractura que es equivalente al espesor del
yacimiento. Es probable que existan trayectorias naturales de alta permeabilidad,
como lo son vetas naturales de alta permeabilidad o tal vez fracturas naturales. Las
cuales se esperan que se encuentren en las proximidades del pozo productor.
134
Ilustración 5-74 Modelo lineal para flujo Horizontal (J.A., 2001).
Tabla 5-9 Información General de los casos de capo de pozos de gas de baja
permeabilidad.
135
Shale C 300 0.05 4 0.25
Tabla 5-11 Datos de los fluidos obtenidos de los pozos productores de gas en formaciones
de baja permeabilidad.
Pozo 𝑷𝒊 T 𝛾𝑔 𝑪𝒕 𝝁𝒈 𝑯𝟐 𝑺 𝑪𝑶𝟐 𝑵𝟐
[psia] [°F] [A] [1/psia [cp] [%] [%] [%]
E^-6]
Tight 6700 220 0.61 104 0.02414 0 0.002 0.1
231 7
Tight 6700 220 0.61 8.91 0.0287 0 0.002 0.01
218
Tight 6700 220 0.61 4.817 0.028 0 0.002 0.01
NBU
Shale A 5100 207 0.596 176.8 0.02395 0 0.04 0.001
1
Shale B 2821. 181. 0.621 310.3 0.019 0 0.015 0.01
5 4
Shale C 5100 210 0.6 189.2 0.004 0 0.02 0.001
Tabla 5-12 resultados calculados de los regímenes de flujo bilineal y lineal tardío.
136
Tabla 5-13 Resultados calculados del régimen de flujo dominado por la frontera externa.
𝑘𝑓 𝑘𝑚 𝐺 𝑘𝑓 𝑘𝑚 𝐺
137
Pozo Lim. Método de Modificado de Método de potencias
Econ Declinación Exponencial
ómic Extendida
o
[años] 𝑛 𝜏 int Error Gp 𝑚 𝐷1 𝑛 𝐷𝑖 𝑞𝑖 Error Gp
Fracci MMscf 1/D 1/D Mscf/ Fracció MMscf
ón D n
Tight 30 0.4 72.9 0.1 0.018 638.3 0.2 0.014 0.8 0.016 3.000 0.070 475.1
231 2 5
Tight 50 0.2 263.8 0.2 0.006 3421 0.5 0.060 0.4 0.137 2.500 0.019 596.3
218 6 7 3 6
Tight 20 0.5 116.4 0.0 0.007 198.5 4.1 57135 3.1 1829 0.320 0.020 2307.
NBU 3 9 7 2 0
Shale 20 0.4 175.6 0.0 0.032 1385. 0.1 0.003 1.1 0.003 3.000 0.037 668.3
A1 9 2 7
Shale 20 0.6 218.7 0.0 0.035 817.3 0.6 0.158 0.4 0.386 30.00 0.119 440.0
B 8 3 2
Shale 15 0.5 144 0.0 0.052 925.8 0.6 0.107 0.4 0.258 10.00 0.190 675.4
C 2 7
Tabla 5-15 Resultados del ajuste de la producción de gas de los modelos de Duong y
SEPD.
Tabla 5-16 Resultados del ajuste de producción de gas de los modelos YM-SEPD y
Pozo Limite
Económico Método de Duong Método de Declinación Exponencial
Extendida
[años] a m 𝑞𝑖 Error Gp 𝑞𝑖 𝜏 𝑛 Error Gp
Mscf/D [Fracción] MMscf Mscf/D [Fracción] MMscf
Tight 30 1.23 2.02 0.99 0.026 955.56 3652.65 0.36 47.8 0.019 780.63
231
Tight 50 1.02 0.84 2.03 0.012 5575.81 3632.02 0.18 15.4 0.011 3635.47
218
Tight 20 1.29 3.33 0.03 0.015 383.72 1170.03 0.18 100 0.017 1417.66
NBU
Shale 20 1.12 1.33 2.75 0.053 2385.83 4838.33 0.36 100 0.036 1929.13
A1
Shale 20 1.65 34.27 8.50002e^- 0.103 1474.82 9035.86 0.16 0.5 0.204 1481.82
B 19
Shale 15 1.20 1.80 2.06 0.069 1390.14 20206 0.20 0.5 0.159 935
C
Potencias.
138
Tabla 5-17 Resultados del ajuste de la producción de gas del modelo de Crecimiento
Logístico.
Haciendo el análisis de las curvas de declinación a los casos de campo de tipo tight
gas se observó que el ajuste de la curva depende de los siguientes factores: la
calidad de los datos de producción (no existan tantas variaciones o cierres en los
datos), la cantidad de datos de producción de gas y el régimen de flujo que domina
en el yacimiento.
En los casos de los pozos de tipo tight se observó un mejor ajuste de producción de
los métodos de declinación exponencial extendida (Valko & lee, 2010) y el modelo
modificado de declinación exponencial extendida (Yu, Lee, Miocevic, Li, & Harris,
2013), respecto a los modelos de Duong (Duong, 2010), Potencias Crecimiento
logístico (Clark, 2011). Como se puede observar en la Fig.5-75 y en la Tabla 5-18
del caso CIGE 231, se aprecia que el modelo modificado de declinación exponencial
extendida presenta un mejor ajuste con un ajuste de error de tan solo un 0.018, el
modelo de Declinación exponencial extendida muestra un buen ajuste con 0.019.
El modelo de Duong presenta una sobre estimación del pronóstico de producción,
esto debido a la poca cantidad de datos de producción de gas y dado que el pozo
ha alcanzo el régimen de flujo dominado por la frontera externa. En el caos de los
modelos de potencias y crecimiento logístico ofrecen una subestimación del
pronóstico de producción del pozo CIGE 231.
En el pozo CIGE 218 se muestra un mejor ajuste por parte del modelo modificado
de declinación exponencial con un error de ajuste del 0.0065. A diferencia del caso
anterior se muestra que los modelos de Declinación exponencial extendida y modelo
de potencias muestran un sobre estimación de la producción de gas, esto debido a
que solo se cuenta con 294 días de producción y esto es una limitante en el modelo
de SEPD, debido a que arroja datos muy conservadores como se ve en la Fig. 5-
76. Los modelos de Duong y Crecimiento logístico nos dan una subestimación de
la producción de gas. A pesar de los pocos datos de producción con que se cuentan
139
el modelo de Duong presenta un ajuste bastante bueno, además de que el pozo se
encuentra en un régimen de flujo transitorio de larga duración. En el último pozo,
NBU, el mejor ajuste se obtuvo del modelo modificado de declinación exponencial
y del modelo de Duong.
Tabla 5-18 Comparación del pronóstico de EUR- CIGE 231
Ilustración 5-75 Comparación de ajuste de la producción de gas del pozo CIGE 231,
mediante curvas de declinación
140
En general se obtuvieron mejores resultados con el modelo modificado de
declinación exponencial, debido a que los tres pozos analizados son similares en
cuestiones de la cantidad de datos de producción.
Yu (Yu, Best Practice of Using Empirical Method for Production Forecast and EUR
Estimation in Tight/Shale Gas Resevoirs, 2013) recomienda el uso del método
Modificado de Declinación exponencial para los yacimientos que cuentan con un
rango de permeabilidades de 0.1 - 0.001[md] , lo cual se comprobó en el análisis
de los pozos CIGE 231 y CIGE 218, los cuales cuentan con unas permeabilidades
de 𝑘 = 0.0139[md] y 𝑘 = 0.0220[md], a pesar de que solo se cuenta con unos
históricos de producción menores de un año. Además, Yu (Yu, Best Practice of
Using Empirical Method for Production Forecast and EUR Estimation in Tight/Shale
Gas Resevoirs, 2013) menciona que para rangos de permeabilidades menores de
0.001[md] se recomienda usar los modelos de Duong y modelo modificado de
Declinación exponencial con un histórico de producción de 3 a 4 años. Pero en el
caso del caso del pozo NBU, notamos que se aplica el mismo criterio, dado que
cuenta con una permeabilidad de 𝑘 = 0.0068[md] y con un histórico de producción
menor de un año.
En su artículo Duong (Duong, 2010) menciona que mediante la pendiente de la
gráfica especializada log-log de 𝑞 ⁄𝐺𝑝 𝑣𝑠 𝑑𝑖𝑎𝑠 podemos diferencias si es un pozo
profundo y una formación tight ,m<1, o un pozo shale, m>1. Si observamos los
datos obtenidos en la Tabla 5-15, podemos apreciar que solo un caso representa
este comportamiento y los otros pozos tipo tight tienen una pendiente mayor a uno,
lo cual no sigue el comportamiento descrito por Duong.
Tabla 5-19 Comparación del pronóstico de EUR-CIGE 218
141
Ilustración 5-76 Grafica de comparación del ajuste de la producción de gas del pozo
CIGE 218.
142
hidráulica). El fluido se encuentra fluyendo dentro de la trayectoria de flujo de alta
permeabilidad en una dirección perpendicular al flujo del yacimiento.
Los resultados obtenidos con el análisis de la de modelo de Arévalo-Castellanos
nos da unos resultados de permeabilidad bastante conservadores. Es por ello por
lo que se optó por hacer el análisis de la permeabilidad de forma individual,
mediante el método de Bello. Los resultados obtenidos son bastante aceptables. Se
realizo el análisis dos casos: arriba de presión de desorción y debajo de presión de
desorción.
En el caso del pozo Shale A se tiene un histórico de producción de 442 días. De los
207 días a los 284 días se da una etapa de cierres de pozo, donde al parecer se da
una fractura miento del pozo. Esto se comprobó al hacer el análisis de la
permeabilidad mediante el método de bello, donde se puede observar que la
permeabilidad de fractura del primer intervalo del histórico fue 𝑘𝑓 = 𝟎. 𝟎𝟎𝟎𝟎𝟓[𝐦𝐝] y
en el segundo intervalo, después del fracturamiento, resulto una permeabilidad de
fractura de 𝑘𝑓 = 𝟎. 𝟎𝟒𝟗𝟖𝟏 [𝐦𝐝]. Claramente existe un aumento en la permeabilidad
de fractura, lo que nos indica que hubo fracturamiento.
Tabla 5-20 Comparación del pronóstico de EUR-Shale A
Pozo Shale A
Método EUR %Error
143
Ilustración 5-77 Grafica de ajuste de la producción de gas del pozo Shale A.
Los métodos que nos dieron mejor ajuste en los pozos shale fueron: método
modificado de declinación exponencial y método de Duong. En el pozo Shale A se
observa que el método de declinación exponencial, método de potencias y el
método modificado de declinación exponencial muestran un buen ajuste. A
diferencia del método de Duong que a tiempo largo de producción muestra una
sobre estimación, al igual que el método de declinación exponencial. El método de
Crecimiento logístico nos da muy malos ajustes y una subestimación de la
producción de gas.
En el caso del pozo Shale C los métodos que mostraron mejor ajuste fueron los
métodos de Duong y el método modificado de declinación exponencial, pero como
en los casos anteriores muestra una sobre estimación de la producción. En este
caso el método de crecimiento muestra un ajuste bastante aceptable y a tiempo
largo de pronóstico es similar a los demás métodos.
En yacimientos con permeabilidad menor a 0.001[md] Yu (Yu, Best Practice of Using
Empirical Method for Production Forecast and EUR Estimation in Tight/Shale Gas
Resevoirs, 2013)recomienda usar el método de Duong y el método modificado de
144
declinación exponencial. En los casos analizados en este trabajo corroboramos las
recomendaciones, pero podemos agregar que el método de Duong genera sobre
estimaciones de la producción a tiempos largos.
Tabla 5-21Comparación del pronóstico de EUR-Shale C
Pozo Shale C
Método EUR %Error
Duong 1390.140 0.069
SEPD 935.000 0.159
YM-SEPD 925.840 0.052
Potencias 675.490 0.190
Crecimiento 1108.16 0.17
logístico
Ilustración 5-78 Grafica comparativa de ajuste de la producción de gas del pozo Shale
C.
145
6. CAPÌTULO VI CONCLUSIONES Y RECOMENDACIONES
6.1. Conclusiones
El objetivo de este trabajo se enfoca en el análisis del comportamiento de la
producción en los yacimientos no convencionales (baja permeabilidad y arenas
compactas, Tight y Shale) mediante varios métodos empíricos y semi- analíticos.
Después del estudio de los 6 campos de baja permeabilidad se obtuvieron las
siguientes conclusiones:
1. Se percibió un comportamiento similar en los pozos de arenas compactas y
lutitas. Donde presentan una gran declinación en las etapas tempranas de
producción, esto se debe a la producción rápida de gas, provenientes de la
fracturas hidráulicas y naturales. Esto resulta en la formación de flujos
transitorios de largo plazo bilineal y lineal, lo cual se da por el flujo de gas de
una formación de baja permeabilidad de matriz a una zona de alta
permeabilidad (Fracturas naturales o inducidas).
2. Se observo que los flujos transitorios de larga duración son una característica
predominante en los yacimientos de baja permeabilidad y pueden llegar a
durar meses e incluso años. Sin embargo, también se observó que algunos
pozos pueden llegar a mostrar un régimen de flujo dominado por la frontera
externa, lo que hace inferir que no existe flujo de fluidos fuera del volumen
estimulado de pozo para el caso de gas de lutitas, por ejemplo.
6. Los pozos de tipo Shale mostraron un únicamente un flujo lineal y/o bilineal
y no muestran flujo pseudo- radial como se espera debido a una fractura
146
miento hidráulico. Esto se debe principalmente las características asiladas
de los micro yacimientos en las formaciones de gas de lutita, la presión
capilar o las curva o las curvas de permeabilidades relativas deben ser
aplicadas a cada pequeño yacimiento, pero no puede ser aplicable a través
de todo el yacimiento. Es por ello por lo que no muestran un flujo radial tardío
o un régimen de flujo tardío.
6.2. Limitaciones
1. El método semi- analítico de Arévalo-castellano solo se aplican a flujo de gas.
2. Las fórmulas semi- analíticas aplicadas para el cálculo de las propiedades
mostradas (ver capítulo 3) se derivaron para las condiciones de flujo
constante de gas,𝑞𝑔 , y presión de fondo fluyendo constante,𝑝𝑤𝑓 .
3. La variación constante en los datos de producción y presión podrían
ocasionar grandes desajustes en la mayoría de los métodos aplicados en
este trabajo.
4. El método de Duong arroja errores muy grandes en yacimientos que ya han
alcanzado la frontera y presentara una sobre estimación de gas a pronósticos
grandes.
5. El método de SEPD muestra un desajuste en yacimientos con un rango de
permeabilidades de 0.01-0.001[md] y con datos de producción mayores de 2
años.
147
6. El método de Potencias calcula errores grandes en yacimientos con
permeabilidades menos a 0.001 [md].
7. Los modelos de curvas de declinación son modelos empíricos y tienen
algunas limitaciones:
• No involucran las restricciones naturales del yacimiento
• Dificultad de explicar el diseño del pozo.
• Los efectos geo mecánicos son ignorados
• No involucran datos geológicos del yacimiento.
• Las predicciones son cuestionables debido a todo lo anterior.
148
8. Complementar en software de Método de Análisis de presión y producción
en yacimientos no convencionales (MAPPYNC) con otras metodologías
empíricas y semi analíticas.
9. Realizar el análisis de más yacimientos No convencionales con diferentes
características y ver su comportamiento en los diferentes modelos de
análisis de presión y producción.
149
NOMENCLATURAS
𝐴 =área de drene del pozo,𝐿2 , [𝑓𝑡 2 , 𝑎𝑐𝑟𝑒𝑠].
150
𝐵 =factor de volumen de la formación.
̂1
𝐷 ̂1
=intercepción de la constante de declinación definida en la Ec. 4.23,[𝐷
= 𝐷1 ⁄𝑛].
151
𝐻2 𝑆 =ácido sulfhídrico, [𝑓𝑟𝑎𝑐𝑐𝑖𝑜𝑛 𝑚𝑜𝑙].
𝐾 =constante de acarreamiento.
𝑘ℎ =capacidad de flujo.
𝑚 =metros , 𝐿.
𝑚𝑑 =milidarcy.
𝑚(𝑝𝑖 ) =pseudo- presión del gas real a la presión inicial, 𝑚/𝐿𝑡 3 , [𝑝𝑠𝑖𝑎2 /𝑐𝑝].
152
𝑚(𝑝𝑤𝑓 ) =pseudo- presión del gas real a la presión de fondo fluyendo
constante, 𝑚/𝐿𝑡 3 , [𝑝𝑠𝑖𝑎2 /𝑐𝑝].
𝑚
̃ 𝐿𝑃𝐶 =pendiente para la constante 𝑝𝑤𝑓 evaluada del comienzo de la línea
recta de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 √𝑡, [𝑝𝑠𝑖𝑎2 − 𝐷1⁄2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝑅𝑃𝐶 =pendiente para la constante 𝑝𝑤𝑓 evaluada del comienzo de la línea
recta de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 log(𝑡), [𝑝𝑠𝑖𝑎2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐵𝑃𝐶 =pendiente para la constante 𝑝𝑤𝑓 evaluada del comienzo de la línea
4
recta de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 √𝑡 [𝑝𝑠𝑖𝑎2 − 𝐷1⁄4 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐵𝑃𝐶 = pendiente para la constante 𝑝𝑤𝑓 evaluada del comienzo de la línea recta
de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 1/√𝑡 [𝑝𝑠𝑖𝑎2 − 1/𝐷1⁄2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐹𝑃𝐶 =pendiente para la constante 𝑝𝑤𝑓 evaluada del comienzo de la línea
recta de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 𝑡, [𝑝𝑠𝑖𝑎2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐿𝐺𝐶 =pendiente para la constante 𝑞𝑔 evaluada del comienzo de la línea recta
de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 √𝑡, [𝑝𝑠𝑖𝑎2 − 𝐷1⁄2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝑅𝐺𝐶 =pendiente para la constante 𝑞𝑔 evaluada del comienzo de la línea recta
de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 log(𝑡), [𝑝𝑠𝑖𝑎2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐵𝐺𝐶 =pendiente para la constante 𝑞𝑔 evaluada del comienzo de la línea recta
4
de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 √𝑡 [𝑝𝑠𝑖𝑎2 − 𝐷3⁄4 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐵𝐺𝐶 = pendiente para la constante 𝑞𝑔 evaluada del comienzo de la línea recta
de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 1/√𝑡 [𝑝𝑠𝑖𝑎2 − 1/𝐷1⁄2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑚
̃ 𝐹𝐺𝐶 =pendiente para la constante𝑞𝑔 evaluada del comienzo de la línea recta
de la grafica ∆𝑚(𝑝)/𝑞𝑔 𝑣𝑠 𝑡, [𝑝𝑠𝑖𝑎2 /𝑀𝑓𝑡 3 − 𝑐𝑝].
𝑛𝑣𝑒𝑡𝑎𝑠 =número de vetas de permeabilidad.
𝑛 =parámetro exponencial,[𝑎𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙].
𝑁 =población,[𝑎𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙].
153
𝑝𝑐𝑒 = presión a condiciones estándar=14.596,𝑚/𝐿𝑡 2 , [𝑝𝑠𝑖𝑎].
154
𝑠 =factor de daño, [𝑎𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙 ].
𝑡 =tiempo, t, [𝑑𝑖𝑎𝑠].
𝑡𝑒𝑠𝑟 =Tiempo en el fin de la línea recta en una gráfica de raíz cuadrada del
tiempo.
𝑉 =Volumen molar,𝐿3 , [𝑓𝑡 3 @c. y. ]
155
𝑥𝑒 =distancia del pozo a la frontera externa, 𝐿, [𝑓𝑡]
𝜙 =porosidad, [𝑓𝑟𝑎𝑐𝑐𝑖ò𝑛].
𝜔 =factor acéntrico.
Subíndices
𝐴 =sistema A.
𝑎 =aparente.
𝑐𝑒 =condiciones estándar.
156
𝑃𝐶 =presión de fondo fluyendo constante
𝑓𝑙𝑟 =tiempo real del final de la línea recta en la gráfica semilogarítmica del tiempo.
𝑓𝑟𝑐 =tiempo real del final de la línea recta en la raíz cuadrada de la gráfica del
tiempo.
𝑓𝑙𝑏 =tiempo real del final de la línea recta en el régimen de flujo bilineal.
𝑓𝑙𝑟 =tiempo real del final de la línea recta en el régimen de flujo radial.
𝑓 =fractura.
𝑓 =formación.
𝑔 =gas.
𝑖 =condiciones iniciales.
𝑖𝑛𝑣 =investigación.
𝑥 =dirección x.
𝑦 =dirección y.
𝑧 =dirección z.
REFERENCIAS
157
1. Abbott, C. A. (1982). Application ofLlinear and spherical flow Analysis
Techniques to field Problems. (SPE 11088).
2. Administration, U. E. (2017). Energy Information Administration. Obtenido
de [Link]
3. Ahmed H. El-Banbi, a. R. (1998). Analysis of Linear Flow in gas Well
Production. (39972).
4. Ahmed H. El-banbi, R. A. (1996). Analysis of linear flow in Gas Well
Production. (SPE 39972).
5. Ahmed, T. (2006). Resevoir Engineering Hanbook. Elsevier.
6. Ahmed, T. (2006). Resevoir Engineering Handbook. GPP.
7. Al-Hussainy, R. (1965). The flow of real gases through poruous media. SPE.
8. Arévalo Villagran, J. A. (2001). Analysis of Long-term Behavior in tight gas
resevoirs:case histories. Texas A&M University.
9. Arévalo Villagran, J. A., Catellanos Páez, F., Martínez Romero, N., & Pumar
Martínez , F. (2014). Dynamic Characterization of Unconventional Gas
[Link] Cases. SPE, 17.
10. Arévalo, J. D. (2009). Analisis de la produccion en yacimientos de gas con
baja permeabilidad. Ciudad de México: UNAM.
11. Arévalo, J. F. (2012). Tecnicas de Modelado para la Simulacion Numerica
en Yacimientos y Pozos de Gas. Tesis Profesional.
12. Arevalo-Villagran J.A., C. P. (2014). Dynamic Characteriation of
Unconventional Gas [Link] Cases. SPE, 17.
13. Arevalo-Villagran J.A., C.-P. F.-R. (2014). Dynamic Characterization of
Unconventional gas Resevoir. FIeld Cases. SPE.
14. Arps, J. J. (1944). Analysis of Decline Curves. Houston : SPE.
15. Aviles, G. (2007). Propiedades de los Hidrocarburos e Identificacion de
Zonas con Fluidos en el Yacimiento Utilizando Registros Geofisicos de
Pozos. Tesis Profesional.
16. Bello, R. O. (2009). Rate Trasient Analysis in Shale Gas Resevoir with
Transient Linear Behaviour . PhD,Texas A&M U, College Station.
17. Castellanos Páez, F. (2015). Caracterización Dinámica de yacimientos no
convencionales de gas. CDMX: UNAM.
18. Clark, A. J. (2011). Decline Curve Analysis in Unconventional Resource
Play Using Logistic Growth Model. USA: University of Texas at Austin.
158
19. De la Cruz, G. (2007). Ingenieria en Produccion de Gas. Tesis profesional.
20. Duong, A. N. (2010). Rate-Decline Analysis for fractured Dominated Shale
Resevoirs. SPE.
21. F., A. V. (2014). Dynamic Characterization of Unconventional Gas
[Link] Cases. SPE.
22. Gringarten, A. R. (1975). Applied pressure Analysis for fractured wells. (SPE
4051).
23. Hagoort, J. (1988). Fundamental of Gas Resevoir Engineering. Elsevier.
24. Holditch, S. A. (2006). Tight Gas Sands. JPT, 7.
25. IIK, D., Perego, A., Rushing, J., & Blasingame, T. (2008b). Exponential vs.
Hyperbolic Decline in TIght Gas Sands-Understanding the Origin and
Implications for Reserve Estimates Using Arps´Decline Curves. (116731).
26. IIk, Rushing, A., J., & Blasingame. (2008a). Estimating Reserves Using the
Arps Hyperbolic Rate-Time Relation-Theoery,Practice and Pitfalls. (CIM
2008-108).
27. J.A. Arévalo-Villagrán, R. W.-V. (2001). Some History Cases of long-term
Linear Flow in tight gas wells. SPE.
28. J.A., A. V. (2001). Analysis of long-term behavior in tight gas resevoirs:case
histories. Texas A&M University.
29. James G. Speight, P. (2013). Shale Gas Production Processes. Elsevier.
30. Lee, W. a. (1981). Fracture Evaluation with Pressure Trasient Testing in
Low-Permeability Gas reservoir. JPT, 271.
31. McKinney, T. A. (2005). Advanced Resevoir Engineering. Burlington, USA:
Elsevier.
32. Nolem-Hoeksema, R. (2013). Elementos de fracturamiento hidraulico.
oilfield Review, 57-58.
33. Nolen-Hoeksema, R. (2013). Definición del concepto de fracturamiento
Hidráulico . Oilfield Review.
34. R., A. (s.f.). Naturally Fractured Resevoirs. PennwellBooks.
35. R., A.-h. (1965). The Flow of Real Gases Through Porous Media. SPE, 13.
36. Rodriguez, F. N. (2014). Sistema de Terminación Mutietapas en un Pozo
del Área de Chicontepec. Ciudad de Mexico : UNAM.
159
37. Samaniego-V., H. C.-L. (1981). Transient Pressure Analysis for Fractured
Wells. SPE.
38. Usman Ahmed, D. N. (s.f.). Unconventional oil and gas resources
exploitation and Development. CRC Press.
39. Valkó , P. P. (2009). Assingning Value to Stimulation in the Barnet Shale:A
Simultaneous Analysis of 7000 Plus Production Histories and Well
Completion Records. SPE, 19.
40. Valko, P. P., & lee, W. J. (2010). A Better Way to Forecast Production From
Unconventional Gas Well. SPE, 16.
41. van Everdingen, A. F. (1949). The aplication of LaPlace Transformation to
flow Problems In Resevoir.
42. Verduzco, D. F. (2014). Flujo de gas a traves de medios porosos.
43. Villagran, J. A. (2001). Analysis of long-term behavior in tight gas
resevoir:cases histories. Texas A&M University.
44. Wattenbarger, J. L. (1996). gas resevoir Engineering. SPE.
45. Wattenbarger, R. O. (2008). Rate Transient Analysis in Naturally Fractured
Shale Gas Rsevoir. SPE.
46. Yu, S. (2013). Best Practice of Using Empirical Method for Production
Forecast and EUR Estimation in Tight/Shale Gas Resevoirs.
Alberta,Canada: SPE.
47. YU, S., & Miocevic, D. J. (2013). An Improved Method to Obtain Reliable
Production and EUR Prediction for Wells with Short Production History in
Tight/Shale Reservoirs. SPE(SPE 168684), 11.
48. YU, S., & Miocevic, D. J. (2013). An improved Method to Obtain Reliable
Production and EUR Prediction for Wells with Short Production History in
Tight/shale Resevoirs. SPE, 11.
49. Yu, S., Lee, W. J., Miocevic, D. J., Li, D., & Harris, S. (2013). Estimating
Proved Reserves in Tight/shale Wells Using the Modified SEPD Method.
(166198).
160
En este apéndice presentaremos la derivación de la ecuación de difusividad para
un flujo de fluido ligeramente compresible en un yacimiento homogéneo. La
ecuación de difusión es un modelo matemático que un principio se desarrolló para
el estudio de transferencia de calor, pero dado que el flujo de fluidos a través de los
medios porosos es directamente análogo al flujo de calor a través de sólidos, se
utilizó para describir el comportamiento del flujo de fluidos en el medio poroso. La
ecuación de difusión se obtiene a partir de los siguientes principios físicos: (a) Ley
de la conservación de la masa (ecuación de continuidad), (b) ecuación de
movimiento (Ley de Darcy), (c) una ecuación de Estado y (d) Ley de la conservación
de la Energía. Para la combinación de estos principios físicos es necesario hacer
las siguientes consideraciones: (1) fluido ligeramente compresible, (2) gradiente de
presión pequeño, (3) flujo cilíndrico, (4) no hay variación vertical de la presión, (5)
son despreciados los efectos gravitacionales cuando se utiliza la ley de
movimiento(flujo isotérmico, régimen laminar),(6) la permeabilidad es constante, (7)
temperatura constante en el yacimiento,(8) la compresibilidad es constante y
pequeña,(9) viscosidad constante, (10) la porosidad del yacimiento cambia con la
presión, pero permanece constante con respecto a la posición, medio poroso
incompresible(espesores uniformes),(11) fluido en una sola fase y (12) la Ley de
Darcy es válida(flujo laminar).
Ecuación de continuidad (Verduzco, 2014) (Al-Hussainy, 1965)
El principio de conservación de masa para un fluido isotérmico a través del medio
poroso es expresado por la bien conocida ecuación de continuidad. Para flujo en un
sistema 𝑥, 𝑦, 𝑧 se tiene el volumen elemental representativo siguiente:
Ilustración 0-1 Volumen de control del medio poroso. (Wattenbarger J. L., 1996)
161
(𝑚𝑎𝑠𝑎 𝑞𝑢𝑒 𝑒𝑛𝑡𝑟𝑎) − (𝑚𝑎𝑠𝑎 𝑞𝑢𝑒 𝑠𝑎𝑙𝑒) = (𝐴𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖𝑜𝑛) (A-
40)
Para el caso de la Fig. 10-1 se expresa de la siguiente manera
La Ec. A-7 es la ecuación de continuidad para flujo de fluidos en una sola fase.
162
[Link] Ecuación de Estado
La compresibilidad isotérmica de un fluido se define:
1 𝜕𝑉 (A-
𝑐=− ( ) 47)
𝑉 𝜕𝑝 𝑇
163
Para las condiciones de flujo en un sistema 𝑥, 𝑦, 𝑧 la Ec. A-14 puede escribirse en
función de la presión empleando la Ec. A-15:
𝑘𝑥 𝜌 𝜕Φ 𝑘𝑥 𝜕p (A- 55)
𝑢𝑥 = − =−
𝜇 𝜕𝑥 𝜇 𝜕𝑥
𝑘𝑦 𝜌 𝜕Φ 𝑘𝑦 𝜕p
𝑢𝑦 = − =−
𝜇 𝜕𝑦 𝜇 𝜕𝑦
𝑘𝑧 𝜌 𝜕Φ 𝑘𝑧 𝜕p
𝑢𝑧 = − = − ( + 𝑔𝜌)
𝜇 𝜕𝑧 𝜇 𝜕𝑧
𝑝𝑀 𝑘 𝑝𝑀 𝜕 𝑝𝑀 (A- 59)
∇∙( ∗ ∆𝜌) − (𝑞 ∗ ) = − (𝜙 )
𝑧𝑅𝑇 𝜇 𝑧𝑅𝑇 𝐼 𝜕𝑡 𝑧𝑅𝑇
𝑝 𝑝 𝜕 𝑝 (A- 60)
∇ ∙ (𝑘 ∇ρ) − (𝑞 ∗ ) = − (𝜙 )
𝜇𝑧 𝑧 𝐼 𝜕𝑡 𝑧
𝑝 𝑞∗ 𝑝 𝜙 𝜕 𝑝 (A- 61)
∇ ∙ (𝑘 ∇ρ) − ( ) = ( )
𝜇𝑧 5.615 𝑧 𝑃 5.615 𝜕𝑡 𝑧
Donde:
𝑘 = 𝑝𝑒𝑟𝑚𝑒𝑎𝑏𝑖𝑙𝑖𝑑𝑎𝑑, 𝑝𝑒𝑟𝑚(1 𝑝𝑒𝑟𝑚 = 1.127 𝐷𝑎𝑟𝑐𝑦𝑠)
164
𝑝 = 𝑝𝑟𝑒𝑠𝑖𝑜𝑛, 𝑙𝑏⁄𝑝𝑔2
𝜇 = 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑑𝑎𝑑, 𝑐𝑝
3
𝑃𝑐𝑦
𝑞 = 𝑔𝑎𝑠𝑡𝑜, ⁄
𝐷
𝑡 = 𝑡𝑖𝑒𝑚𝑝𝑜, 𝑑𝑖𝑎𝑠
𝑦 5.615 𝑒𝑠 𝑙𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝑐𝑜𝑛𝑣𝑒𝑟𝑠𝑖𝑜𝑛𝑑𝑒 𝑝𝑖𝑒𝑠 3 𝑎 𝑏𝑎𝑟𝑟𝑖𝑙.
Entonces
𝑃𝑐𝑒 𝑧𝑇 ∗ (A- 67)
𝑞𝑝∗ = ( )𝑞
𝑇𝑐𝑒 𝑃 𝑐𝑒
𝑝 ∗
𝜙 𝜕 𝑝 (A- 69)
∇ ∙ (𝑘 ∇ρ) − 𝐶1 𝑇 𝑞𝑐𝑒 = ( )
𝜇𝑧 5.615 𝜕𝑡 𝑧
Donde
1 𝑃𝑐𝑒 1 14.7 𝑃𝑆𝐼𝐴 𝑃𝑆𝐼𝐴 (A- 70)
𝐶1 = = = 0.00503 [ 3 ]
5.615 𝑇𝑐𝑒 5.615 (460 + 60)°𝑅 𝑝 °𝑅
165
Flujo de Gases Ideales
Considerando momentáneamente flujo únicamente en la dirección x, el termino 𝑝∇ρ,
se puede expresar:
1 (A- 71)
𝑝∇ρ = 𝑑𝑝2 ó 𝑝∇ρ
2
Para estas condiciones de flujo, el primer término del lado izquierdo de la Ec. A-30
es:
𝑑𝑝 1 𝑑𝑝2 (A- 73)
=
𝑑𝑥 2𝜇𝑧 𝑑𝑥
1 2 2 1 𝐶1 𝑇 ∗ 𝜙 𝜕 𝑝 (A- 74)
∇ 𝑝 + ∆𝑝2 ∇ ∙ ( )− 𝑞𝑐𝑒 = ( )
2𝜇𝑧 2𝜇𝑧 𝑘 5.615 𝜕𝑡 𝑧
166
1. El producto 𝜇𝑧 es aproximadamente constante.
2. Los gradientes de presión al cuadrado son pequeños,∇p2 ≪< 1,y (∆𝑝2 )2 →
0.
3. El gas se comporta idealmente (z=1) y la viscosidad es constante e
independiente a la presión.
2 2
𝐶1 𝑇 ∗ 𝜙𝜇𝑐𝑔 (𝑝) 𝜕𝑝2 (A- 77)
∇ 𝑝 − 2𝜇𝑧 𝑞 =
𝑘 𝑐𝑒 5.615 𝑘 𝜕𝑡
167
𝑃
𝑘𝑝 (A- 81)
𝑃𝑝 (𝑝) = 2 ∫ 𝑑𝑝
𝑃0 𝜇 𝑍
Donde:
𝑝0 es una presión de referencia.
Si 𝑘 no depende de la presión, la Ec. A-41 puede simplificarse:
𝑃
𝑝 (A- 82)
𝑃𝑝 (𝑝) = 2 ∫ 𝑑𝑝
𝑃0 𝜇 𝑍
El potencial 𝑝𝑝 (𝑝) puede estimarse por medio de los dos métodos siguientes:
168
2𝑝2 Δ𝑝 (A- 84)
𝑃𝑝 (𝑝2 ) = ( ) ,𝑝 = 0
𝜇2 𝑧2 2 1
∗
𝜙 𝜇 (𝑝)𝑐𝑔 (𝑝) 𝜕𝑝𝑝 (𝑝) (A- 89)
∆ ∙ (𝑘∇𝑝𝑝 (𝑝)) − 2𝐶1 𝑇𝑞𝑐𝑒 =
5.615 𝜕𝑡
169
𝜙 𝜇(𝑝)𝑐𝑔 (𝑝) 𝜕𝑝𝑝 (𝑝) (A- 90)
∆2 𝑝𝑝 (𝑝) =
𝑘 𝜕𝑡
Considerando en forma rigurosa la situación física en que se tiene flujo de una sola
fase, gas en contacto con agua a condiciones de saturación congénita, la
compresibilidad total 𝑐𝑡 (𝑝) puede expresarse:
𝑐𝑡 (𝑝) = 𝑆𝑔 𝑐𝑔 (𝑝) + 𝑆𝑤 𝑐𝑤 + 𝑐𝑓1 (A- 93)
1 𝜕𝜙 (A- 94)
𝑐𝑓1 = = (1 − 𝜙)𝑐𝑓
𝜙 𝜕𝑝
170
componente en una mezcla es el número de libras-mol de ese componente dividido
por el número total de moles de todos los componentes en esa mezcla. Luego, para
un sistema de 𝑛 componentes, la fracción molar es:
𝑛𝑐
(B-
𝑦𝑖 = 𝑛𝑖 ⁄∑ 𝑛𝑗 95)
𝑗=1
Donde 𝑦𝑖 = fracción molar del componente i, 𝑛𝑖 = numero de libras-mole del
componente i, y 𝑛𝑐 es el numero de componentes en el sistema.
Ley de los gases Ideales.
Para comenzar la discusión del comportamiento de los gases reales, se considera
un gas hipotético llamado gas ideal. Un gas Ideal (Perfecto) presenta las siguientes
propiedades:
• El volumen ocupado por la molécula es insignificante en comparación con el
volumen total ocupado por el gas.
• No existen fuerzas de atracción y repulsión entre las moléculas y las paredes
del contenedor en donde se aloja el gas.
• Los choques entre las moléculas son perfectamente elásticas (no existiendo
perdida de energía interna durante los choques).
Una ecuación que describe la relación entre el volumen ocupado por el gas, la
presión y la temperatura se denomina ecuación de estado (EOS). La forma de la
ecuación de estado del gas ideal se desarrolló a partir de observaciones empíricas
que, para una determinada masa de gas en una temperatura constante, el producto
de la presión y el volumen,𝑝𝑉, es contante (Ley de Boyle) y, para una determinada
masa de gas en una presión constante, la relación 𝑣𝑜𝑙𝑢𝑚𝑒𝑛⁄𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎,𝑉⁄𝑇, es
constante(Ley de Charles). Combinando la Ley de Boyle y la de Charles, se obtiene
la ecuación de estado para un gas ideal:
𝑚 (B-
𝑝𝑉 = 𝑛𝑅𝑇 = ( ) 𝑅𝑇
𝑀 96)
171
temperatura(60 °F y 14.65 psia).Este concepto implica que, para las condiciones
estándares, el volumen molar es constante y puede ser usado para convertir masa
a volumen o bien, como en algunas ocasiones se requiere, para convertir un
volumen dado a condiciones estándares a masa.
172
Varias propiedades de los gases tienen valores iguales para gases similares (tales
como los hidrocarburos parafínicos) en cuanto a presión y temperatura reducida.
Presión reducida y temperatura reducida para componentes puros se definen como:
𝑝 (B-
𝑝𝑟 = ⁄𝑝𝑐
100)
Y
𝑇𝑟 = 𝑇⁄𝑇 (B-
𝑐 101)
Respectivamente. La presión pseudo- reducida y la temperatura pseudo- reducida
para mezclas se define como:
𝑝 (B-
𝑝𝑝𝑟 = ⁄𝑝𝑝𝑐
102)
Y
El punto critico (𝑝𝑐 , 𝑇𝑐 ) para una sustancia pura es la presión y temperatura a la cual
las propiedades de las fases de líquido y vapor serán idénticas. A presiones por
arriba de 𝑝𝑐 ,el líquido y el gas no pueden coexistir, sea cual sea la temperatura, a
temperaturas por encima de 𝑇𝑐 , la substancia no puede ser licuada, sea cual sea la
presión. Para sustancia puras, 𝑝𝑐 y 𝑇𝑐 están determinadas experimentalmente. Para
mezclas, 𝑝𝑝𝑐 y 𝑇𝑝𝑐 ambas son obtenidas con alguna consistencia de reglas de
mezclado o están estimadas por correlaciones. Estos valores computados de 𝑝𝑝𝑐 y
𝑇𝑝𝑐 no son verdaderos puntos críticos: por ejemplo, las propiedades de las fases de
líquido y vapor no serán idénticas en el punto (𝑝𝑐 , 𝑇𝑐 ).
La observación segura de las propiedades del gas, tal como el factor 𝑧, serán
aproximadamente las mismas a obtener la temperatura y presión reducidas para
una sustancia pura y forma similar de gases por medio de las bases del principio de
correspondencia de estados. Este comportamiento también será observado para
mezclas y gases químicamente similares; por tanto, correlaciones del factor 𝑧 para
gases puros y mezclas de gas están basadas en este principio.
Propiedades de los gases naturales. (Wattenbarger J. L., 1996)
Al analizar las características de un yacimiento se deben comprender todas las
propiedades que interactúan para el movimiento de fluidos (sistema roca fluidos),
173
debido a que, en la industria petrolera, este es un aspecto de gran importancia, ya
sean en medios porosos (para la extracción o inyección de estos al yacimiento) o
por tuberías (para el transporte de fluidos), en este caso se describen las
propiedades del gas y la manera en que se obtienen.
Una propiedad es cualquier característica medible en una región cerrada por limites
imaginarios, las propiedades dependen del tamaño y de la masa total de los límites
fijados, algunas propiedades se definen en términos de otras, las propiedades
pueden ser clasificadas en extensivas e intensivas, las propiedades extensivas
dependen de la cantidad de masa como volumen, la masa, etc. Las propiedades
intensivas son independientes de la cantidad de masa presente como la
temperatura, la presión, etc. Las propiedades específicas que están definidas por
unidad de masa son propiedades intensivas.
Ahora bien, un fluido es una sustancia que se deforma continuamente al ser
sometido a un esfuerzo cortante (esfuerzo tangencial) no importa cuán pequeño sea
este. Con esta definición queda entendido que un fluido puede ser liquido o gas.
El gas es una fase de la materia homogénea de baja densidad y viscosidad, que no
tiene volumen ni forma fija, toma la forma y el volumen del recipiente que lo contiene,
una de las propiedades más evidentes de los gases es que ejercen presión sobre
sus alrededores en el cuerpo que los contiene. Está conformado por varios
componentes, tales como el nitrógeno, bióxido de carbono, ácido sulfhídrico,
metano, etano, propano, isobutano, n-butano, isopentano, n-pentano, n-hexano, n-
heptano, n-octano, n-nonano, n-decano, oxígeno, hidrogeno, helio y vapor de agua.
Cada uno de estos componentes poseen propiedades físicas. Además, se
encuentran en un cierto porcentaje de la mezcla de gas; por tanto, las propiedades
de cada componente constituyen las propiedades del gas.
Las propiedades físicas de los compontes puros que se presentan en los gases,
evaluadas a condiciones estándar de 14.7 𝑙𝑏⁄𝑝𝑔2 𝑎𝑏𝑠 y 60 °𝐹, se enlistan en la
Tabla 11-1. Estas propiedades, incluyen la formula química, el peso molecular,
temperatura y presión critica, densidades del líquido y gas y viscosidad del gas
(Para componentes más ligeros que el pentano). Estas propiedades de
componentes puros son usadas en cálculos basados en las reglas de mezclado
para desarrollar pseudo- propiedades para mezclas de gas, incluyendo el peso
molecular aparente y gravedad especifica del gas.
Peso molecular aparente para una mezcla de gas.
Debido a que la mezcla de gas está compuesta de moléculas de varios tamaños y
de distintos pesos moleculares. Una mezcla de gas se comporta como si tuviera
definido peso molecular. Esta observación del peso molecular para una mezcla de
gas con 𝑛𝑐 componentes es llamado aparente o peso molecular promedio molar y
se determina por:
174
𝑛𝑐
(B-
𝑀 = ∑ 𝑦𝑖 𝑀𝑖 104)
𝑖=0
Donde 𝑀 = peso molecular aparente de la mezcla de gas, [𝑙𝑏⁄𝑙𝑏𝑚 − 𝑚𝑜𝑙]; 𝑀𝑖 =
peso molecular del componente i del gas, [ 𝑙𝑏⁄𝑙𝑏𝑚 − 𝑚𝑜𝑙 ]; y 𝑦𝑖 = fracción mol de la
fase de gas del componente i [Fracción].
Tabla 0-1 Tabla de propiedades físicas de los gases a 14.7 psia y 60 °F.
Peso
Temperatura Presión Densidad Densidad Viscosidad
Fórmula molecular
Componente crítica crítica del líquido del gas del gas
química (lbm/ lbm-
(°R) (lb/pg2abs) (lbm/ ft3) (lbm/ ft3) (cp)
mol)
Hidrógeno H2 2.109 59.36 187.5 4.432 0.005312 0.00871
Monóxido de
CO 28.010 239.26 507.5 49.231 0.073810 0.01725
carbono
Nitrógeno N2 28.013 227.16 493.1 50.479 0.073820 0.01735
Sulfuro de
H2S 34.08 672.35 1,306.0 49.982 0.089810 0.01240
hidrógeno
Dióxido de
CO2 44.010 547.58 1,071.0 51.016 0.116000 0.01439
carbono
Aire - 28.963 238.36 546.9 54.555 0.076320 0.01790
175
Se dan valores de densidad de los líquidos para estos componentes, que puedan existir como líquidos a 60 °F y 14.7
(lb/pg2abs), se estima la densidad del líquido para componentes que son gases naturales a estas condiciones.
Se define como la relación de la densidad del gas y la densidad del aire, ambas
están medidas a determinada presión y temperatura:
𝜌𝑔 (B-
𝛾𝑔 = ⁄𝜌𝑎
105)
176
desarrollados con anterioridad (tales como el procedimiento de kay (Wattenbarger
J. L., 1996)). Las reglas de mezclado de Stewart proveen regularmente más
precisión en sus resultados. El segundo método, desarrollado por Sutton, estima las
propiedades pseudo- criticas cuando la composición del gas no es conocida. El
método de Sutton requiere de manera considerable menores cálculos aritméticos
en comparación con el método de reglas de mezclado de Stewart, y es el método
de preferencia para cuando la rapidez en los cálculos es importante, sin importar
mucho la precisión en los mismos. También es usada la gravedad especifica del
gas para una constante en lugar de una detallada composición del hidrocarburo. El
método de Sutton es más preciso que el de las reglas de mezclado de kay.
También presentaremos correlaciones de corrección de la temperatura y presión
pseudo-criticas debido a la presencia de contaminantes comúnmente asociados con
la producción del gas natural. Tales contaminantes son el ácido sulfhídrico (𝐻2 𝑆),
bióxido de carbono (𝐶𝑂2), nitrógeno, y vapor de agua. Mas aun, se mostrará una
técnica de cálculo para estimar la gravedad especifica del gas-húmedo y gas y
condensado. Esta gravedad especifica también puede ser usada en el método de
Sutton para estimar las propiedades pseudo-criticas.
Estimación de las propiedades pseudo-criticas cuando la composición del
gas es conocida: reglas de mezclado de Stewart.
Stewart y otros autores (Wattenbarger J. L., 1996) trabajaron con 21 diferentes
reglas de mezclado concluyendo su trabajo mostrado las Ec. B-19 hasta las Ec. B-
24. Estas reglas de mezclado proveen los resultados más consistentes usando una
simple cubica regla de mezclado con datos experimentales, en comparación con los
resultados computarizados que también son muy confiables. El término “simple
cubica” se refiere a las ecuaciones de estado en forma cubica (Ecuaciones de
Estado de van ser Waals y Redlich-kwong). Debido a que estas reglas de mezclado
dan resultados más precisos, el método de Stewart será usado para estimar
presiones y temperaturas pseudo- críticas para la obtención de factor 𝑧, la
compresibilidad del gas y la viscosidad.
Con el método de Sutton también se puede hacer correcciones para componentes
de mayor peso molecular. Las ecuaciones de Kessler-Lee (Wattenbarger J. L.,
1996)(Ec. B-13 hasta B-15) son usadas para calcular las propiedades críticas de la
fracción del heptano plus (𝐶7+ ).
Procedimiento de cálculo con el método de Stewart.
1. Si una fracción significante de componentes pesados (𝐶7 𝑦 𝑝𝑒𝑠𝑎𝑑𝑜𝑠)está
presente en la mezcla de gas natural, mediciones en el laboratorio para la
obtención del peso molecular y de la gravedad específica para el 𝐶7+ son
necesarios para su utilización en las reglas de mezclado para el cálculo de la
gravedad de la mezcla y de las propiedades pseudo-criticas. Las ecuaciones
177
de Whitson y Kessler-Lee (Wattenbarger J. L., 1996) son recomendadas para
estimar las propiedades críticas de la fracción del 𝐶7+ .
𝑇𝑐 2 3 )
𝜉𝑘 = ( ) (0.3129 𝑦𝐶7+ − 4.8156 𝑦𝐶7+ + 27.3751𝑦𝐶7+
√𝑃𝑐 𝐶7+
3. Obtener las presiones críticas y las temperaturas de los
componentes remanentes de la Tabla 11-1.
4. Determinar las presiones y temperaturas pseudo- críticas del gas.
A. Calcular los parámetros J y K.
𝑛𝑐 𝑛𝑐 2 (B- 110)
1 𝑦𝑇𝑐 2 𝑇𝑐
𝑗 = ∑ ( ) + [∑ (𝑦√ ) ]
3 𝑃𝑐 𝑖 3 𝑃𝑐
𝑖=1 𝑖=1
𝑖
Y
178
𝑛𝑐
𝑦 𝑇𝑐 (B-
𝐾 = ∑( ) 111)
𝑖=1
√𝑃𝑐 𝑖
B. Los parámetros de correlaciones 𝐽´ 𝑦 𝐾´ para la fraccion de 𝐶7+ .
𝐽´ = 𝐽 − 𝜉𝑗 (B-
112)
Y
𝐾´ = 𝐾 − 𝜉𝑘 (B-
113)
C. Calcular la temperatura y presión pseudo- critica.
𝑇𝑝𝑐 = 𝐾 ´2 ⁄𝐽´ (B-
114)
Y
𝑝𝑝𝑐 = 𝑇𝑝𝑐 ⁄𝐽´ (B-
115)
Estimación de las propiedades pseudo- criticas cuando la composición del
gas no es conocida: correlación de Sutton.
El método propuesto por Stewart para calcular las propiedades pseudo- criticas
requiere información sobre la composición del gas, sin embargo, análisis de
laboratorio a menudo no están disponibles. Usando datos de 264 muestras de gas,
Sutton desarrollo una correlación para estimar la presión y temperatura pseudo-
critica como una función de la gravedad especifica del gas. Las curvas de
correlación de Sutton son mostradas en la Fig. 11-1, estas están basadas en una
mayor base de datos que las usadas por Standing (Wattenbarger J. L., 1996) y por
consecuencia difieren significativamente de las curvas de Standing. Sutton ajusto
los datos en bruto con ecuaciones cuadráticas y obtuvo las siguientes ecuaciones
empíricas relacionadas con las propiedades pseudo- críticas de los hidrocarburos a
la gravedad especifica del gas:
𝑝𝑝𝑐ℎ = 756.8 − 131.0𝛾ℎ − 3.6𝛾ℎ2 (B-
116)
Y
𝑇𝑝𝑐ℎ = 169.2 − 349.5𝛾ℎ − 74.0𝛾ℎ2 (B-
117)
179
Donde 𝑝𝑝𝑐ℎ =presión pseudo-critica de los componentes hidrocarburos, [psia];
𝑇𝑝𝑐ℎ = temperatura pseudo- critica de los componentes hidrocarburos, [°R]; y
𝛾ℎ =gravedad especifica del gas de los componentes hidrocarburos(aire=1.0).
Las Ecs. 24,25 y la Fig. 11-1 son aplicables para 0.57 < 𝛾ℎ < 1.68. Y si el gas
contiene < 12 % mol de CO2 ,< 3% mol de nitrógeno, y no contiene 𝐻2 𝑆, entonces 𝛾ℎ
puede ser determinada de la siguiente forma:
1. Si el gas es seco (no forma condensados), y si el separador de gravedad del
gas es usado, entonces 𝛾ℎ = 𝛾𝑔 .
2. Si la gravedad del fluido en la corriente del pozo es computarizada, entonces
𝛾ℎ = 𝛾𝑤 .
Sin embargo, si el gas contiene >12% mol de CO2 , > 3% mol de nitrógeno, o algo
de 𝐻2 𝑆, entonces la gravedad del gas de los hidrocarburos puede ser calculada
por:
180
𝛾𝑤 − 1.1767 𝑦𝐻2 𝑆 − 1.5196 𝑦𝐶𝑂2 − 0.9672 𝑦𝑁2 − 0.6220 𝑦𝐻2 𝑂 (B-
𝛾ℎ = 118)
1 − 𝑦𝐻2 𝑆 − 𝑦𝐶𝑂2 − 𝑦𝑁2 − 𝑦𝐻2 𝑂
Donde 𝛾𝑤 = 𝛾𝑔 si el separador de gravedad del gas ese está utilizando.
Una vez que la gravedad especifica del gas de los componentes hidrocarburos es
estimada, las propiedades pseudo- críticas de la mezcla de los hidrocarburos serán
calculadas con la correlación de Sutton dada por las Ecs. B-24 y B-25 o en la Fig.
11-1. Las propiedades pseudo- críticas de toda la mezcla, que incluyen los
contaminantes, se estiman con las siguientes ecuaciones.
𝑝𝑝𝑐 = (1 − 𝑦𝐻2 𝑆 − 𝑦𝐶𝑂2 − 𝑦𝑁2 − 𝑦𝐻2 𝑂 )𝑇𝑝𝑐ℎ + 1306 𝑦𝐻2 𝑆 + 1071 𝑦𝐶𝑂2 (B-
+ 493.1 𝑦𝑁2 + 3200.1 𝑦𝐻2 𝑂 119)
Y
𝑇𝑝𝑐 = (1 − 𝑦𝐻2 𝑆 − 𝑦𝐶𝑂2 − 𝑦𝑁2 − 𝑦𝐻2 𝑂 )𝑇𝑝𝑐ℎ + 672.35 𝑦𝐻2 𝑆 + 547.58 𝑦𝐶𝑂2 B-
+ 227.16 𝑦𝑁2 + 1164.9 𝑦𝐻2 𝑂 120
Donde los coeficientes de la presión critica (Ec. B-27) y temperatura critica (Ec. B-
28) son las fracciones mol de los contaminantes. Notar que la formas de las Ec. B-
26 hasta las Ec. B-28 inicialmente propuestas por Standing no contienen
correcciones para vapor de agua.
Notar también que la presión temperatura pseudo-criticas calculadas con las Ecs.
B-27 y B-28 no son las correctas si la mezcla del gas está contaminada con los
componentes no hidrocarburos. Correcciones para contaminantes del gas
natural, 𝐶𝑂2 , 𝐻2 𝑆,nitrógeno y vapor de agua son visualizadas más adelante.
181
Donde la temperatura 𝑇´𝑝𝑐 , y la presión 𝑝´𝑝𝑐 pseudo-criticas, ajustadas por
contaminación de 𝐻2 𝑆 𝑦 𝐶𝑂2 son:
𝑇´𝑝𝑐 = 𝑇𝑝𝑐 − 𝜉 (B-
122)
Y
𝑝𝑝𝑐 𝑇´𝑝𝑐 (B-
𝑝´𝑝𝑐 = 123)
[𝑇𝑝𝑐 + 𝐵(1 + 𝐵)𝜉]
182
Ilustración 0-2 Corrección de las propiedades pseudo- criticas 𝑯𝟐 𝑺 𝒚 𝑪𝑶𝟐 .
(Wattenbarger J. L., 1996)
183
La temperatura y presión pseudo-criticas corregidas son:
𝑇𝑝𝑐 ´ − 227.1𝑦𝑁2 − 1165𝑦𝐻2 𝑂 (B-
𝑇𝑝𝑐 " = + 𝑇𝑝𝑐,𝑐𝑜𝑟 126)
1 − 𝑦𝑁2 − 𝑦𝐻2 𝑂
𝑝𝑝𝑐 ´ − 493.1 𝑦𝑁2 − 3200 𝑦𝐻2 𝑂 (B-
𝑝𝑝𝑐 " = + 𝑝𝑝𝑐,𝑐𝑜𝑟 127)
1 − 𝑦𝑁2 − 𝑦𝐻2 𝑂
Donde 𝑇𝑝𝑐 ´ y presión 𝑝𝑝𝑐 ´ son la temperatura y presión pseudo-critica corregida por
𝐻2 𝑆 𝑦 𝐶𝑂2 con la correlación de Wichert y Aziz. Si no hay 𝐻2 𝑆 𝑜 𝐶𝑂2 en la mezcla de
gas, entonces 𝑇𝑝𝑐 ´ = 𝑇𝑝𝑐 , y presión 𝑝𝑝𝑐 ´ = 𝑝𝑝𝑐 .
184
A. Si el gas no contiene nitrógeno o vapor de agua, entonces 𝑝𝑝𝑐 ” = 𝑝𝑝𝑐 ´
y 𝑇𝑝𝑐 ” = 𝑇𝑝𝑐 ´.
B. Si el gas contiene nitrógeno y/o vapor de agua, entonces calcular
𝑝𝑝𝑐 ” 𝑦 𝑇𝑝𝑐 ” con los métodos presentados en la sección pasada.
4. 𝑃𝑝𝑐 y Tpc son valores apropiados para el uso de las correlaciones para obtener
el factor 𝑧, la compresibilidad del gas y la viscosidad del gas.
Correlación para obtener el factor 𝒛 por Dranchuk y Abou-kassem.
Standing y Katz (Wattenbarger J. L., 1996) presentaron una correlación grafica del
factor 𝑧 para los gases naturales como una función de la presión pseudo-
reducida,𝑝𝑟 , y temperatura pseudo-reducida,𝑇𝑟 . Dranchuk y Abou-Kassem
acomodaron 11 constates de las ecuaciones de estados con datos de Standing y
Katz y extrapolaron estas correlaciones para altas presiones reducidas. En la Fig.
11-3 y 11-4, están las gráficas del factor 𝑧 vs 𝑝𝑟 calculadas por Dranchuk y Abou-
Kassem con su ecuación de estado.
Dranchuk y Abou-Kassem desarrollaron su ecuación de estado principalmente para
estimar el factor 𝑧 con rutinas de equipo. Ellos trabajaron con la ecuación de estado
a 1500 puntos de datos con un promedio de error absoluto de 0.486%, cuando el
factor 𝑧 es una función de𝑇𝑟 𝑦 𝑝𝑟 . El factor 𝑧 basado sobre esta ecuación de estado
es exacto dentro de los estándares usuales de ingeniería por 0.2 ≤ 𝑝𝑟 ≤ 30 y 1.0 <
𝑇𝑟 ≤ 3.0, y 𝑝𝑟 < 1.0 y 0.7 < 𝑇𝑟 ≤ 1.0.
Se recomienda el siguiente procedimiento para estimar el factor 𝑧 con la correalcion
de Dranchuk y Abou-Kassem.
1. Calcular las propiedades pseudo-criticas corregidas para
𝐻2 𝑆, 𝐶𝑂2 , 𝑁2 𝑦 𝐻2 𝑂. 𝑝𝑝𝑐 ” 𝑦 𝑇𝑝𝑐 ” usando el procedimiento presentado en la
sección anterior.
2. Calcular las propiedades reducidas,𝑝𝑟 = 𝑝⁄𝑝𝑝𝑐 ” y 𝑇𝑟 = 𝑇⁄𝑇𝑝𝑐 ”.
3. Estimar el factor 𝑧 usando la Fig. 11-3 y 11-4.
𝐵𝑔 = 𝑉𝑅 ⁄𝑉𝑆𝐶 (B-
128)
185
𝑉𝑅 = 𝑍𝑛𝑅𝑇⁄𝑝 (B-
129)
Ilustración 0-3 Factor z para gases naturales con 𝟎 < 𝒑𝒓 < 𝟏𝟎. (Wattenbarger J. L.,
1996)
186
Ilustración 0-4 Factor z para gases naturales con 𝟗 < 𝒑𝒓 < 𝟐𝟎.
(Wattenbarger J. L., 1996)
187
Asumimos las condiciones estándares de 𝑇𝑆𝐶 = 60°𝐹 = 519.67°𝑅, 𝑝𝑆𝐶 =
14.65 𝑝𝑠𝑖𝑎, 𝑦 𝑧𝑆𝐶 = 1
𝑍𝑇(14.65 [𝑝𝑠𝑖𝑎]) 𝑧𝑇 𝑓𝑡 3 (B-
𝐵𝑔 = = 0.0282 132)
(1)(519.65 [°𝑅])𝑝 𝑝 𝑠𝑐𝑓
188
Compresibilidad del gas.
La definición del coeficiente isotérmico de compresibilidad, o simplemente
compresibilidad, es:
1 𝜕𝑉 1 𝜕𝐵𝑔 1 𝜕𝜌 (B-
𝑐=− ( ) =− ( ) =− ( ) 138)
𝑉 𝜕𝑃 𝑇 𝐵𝑔 𝜕𝑝 𝑇 𝜌 𝜕𝑝 𝑇
𝑐𝑟 𝑇𝑟 (B-
𝑐𝑟 =
𝑇𝑟 142)
Mattar y otro autores (Wattenbarger J. L., 1996) utilizaron 11 constantes de
ecuaciones de estado para generar los dos esbozos de la Fig.11-5, en donde el
producto de la compresibilidad pseudo- reducida y temperatura pseudo-
reducida,𝑐𝑟 𝑇𝑟 , se grafica como una función de la presión pseudo-reducida ,𝑝𝑟 , y la
temperatura pseudo-reducida,𝑇𝑟 . La temperatura pseudo-reducida y la presión
pseudo-reducida para una mezcla de gas están definidas por las ecuaciones B-8 y
B-9, respectivamente. Los cálculos son 𝑐𝑔 y el factor 𝑧 son aplicables sobre el rango
0.2 ≤ 𝑝𝑟 < 30 𝑦 1.0 < 𝑇𝑟 ≤ 3.0 𝑦 𝑝𝑟 < 1.0 𝑦 0.7 < 𝑇𝑟 < 1.0. Sin embargo, este
método es pobre en cuanto a resultados en t𝑇𝑟 = 1.0 y 𝑝𝑟 > 1.0.
189
Ilustración 0-5 Variación de la 𝒄𝒓 𝑻𝒓 para gases naturales con 𝟏. 𝟓 < 𝑻𝒓 ≤ 𝟏. 𝟒(esbozo
izquierdo) y variación en 𝒄𝒓 𝑻𝒓 para gases naturales con 𝟏. 𝟒 < 𝑻𝒓 ≤ 𝟑(esbozo
derecho). (Wattenbarger J. L., 1996)
190
986.4 (B-
𝑋 = 3.448 + + 0.01009𝑀
𝑇 146)
191
ecuación de Darcy. Conforme el número de Reynolds se incrementa, las fuerzas de
inercia se incrementan de manera significativa, y el flujo no se modela ya por la
ecuación de Darcy. A esto se le nombra efecto de flujo no-Darcy. Muchos
experimentos (Wattenbarger J. L., 1996) se han realizado para mediciones de 𝛽, en
cuanto a las propiedades de la roca. Utilizando 355 muestras de arena y 29
muestras de calizas, jones determino experimentalmente 𝛽 y desarrollo
correlaciones describiendo 𝛽 como función de la porosidad y de la permeabilidad.
La siguiente correlación se recomienda para estimar 𝛽.
𝛽 = 1.88𝑥1010 𝐾 −1.47 𝜙 −0.53 (B-
150)
El componente no-Darcy en la ecuación de flujo es significativo solamente en el
área de alta velocidad y de alta caída de presión, cercana al agujero del pozo, por
lo tanto, el efecto de flujo no-Darcy usualmente es incorporado dentro de las
ecuaciones de flujo de fluidos como un factor de daño adicional como gasto
dependiente. EL factor de daño total, el cual es el valor determinado en el análisis
de presión transitoria, es un valor aparente,𝑠´,que incluye ambos conceptos: el
factor de daño,𝑠, y el termino que representa los efectos de flujo no-Darcy,𝐷𝑞 .𝛽 es
incorporado dentro del termino 𝐷.
2.715𝑥10−15 𝛽𝑘𝑀𝑝𝑆𝐶 (B-
𝐷= 151)
ℎ𝑟𝑤 𝑇𝑆𝐶 𝜇𝑔 , 𝑤𝑓
En donde 𝐷 = coeficiente de flujo no-Darcy,[𝑀𝑓𝑡 3 @ c. e. ]−1;𝜇𝑔 , 𝑤𝑓 =presión
dependiente de la viscosidad del gas evaluada con la presión de fondo fluyendo,
[cp]; y 𝑀 =peso molecular del gas, [𝑙𝑏𝑚⁄𝑙𝑏 − 𝑚𝑜𝑙]. El coeficiente de flujo no-Darcy
no es contante, pero varía en función de la presión. Después del comienzo de la
explotación, el radio de la región de flujo con alta velocidad inicialmente se
incrementa con el tiempo y pronto se estabiliza.𝐷 es inversamente proporcional a la
viscosidad del gas evaluada en 𝑝𝑤𝑓 .La viscosidad es directamente relacionada a la
presión, lo cual significa que está en función del tiempo. Como la presión declina, la
viscosidad asa mismo declina, se causa por tanto un incremento en 𝐷. Debido a
que 𝐷 no es constante, el análisis en cuanto a pruebas de potencial del pozo con
métodos desarrollados para flujo de líquidos trae consigo resultados erróneos si el
flujo no-Darcy afecta la respuesta de la presión.
192
Apéndice C. DERIVACIÓN DE LAS ECUACIONES DE INTERPRETACIÓN
PARA LOS DIFERENTES REGÍMENES DE FLUJO DETECTADOS EN
YACIMIENTOS DE GAS CON BAJA PERMEABILIDAD.
193
(158)
𝑘𝑡
𝑥𝑒 = 0.1125√
(𝜙𝜇𝑐𝑡 )𝑖
La distancia a la frontera externa debe ser considerado el mínimo valor si todos los
datos aún están en la mitad de la pendiente de la tendencia. Otra vez, en estas
ecuaciones es necesario conocer el valor de la permeabilidad,𝑘. Este es un punto
débil en determinar la distancia a la frontera externa con este método.
Calculo del área de drene, A.
Para el cálculo del área de drene no es necesario conocer el valor de la
permeabilidad, 𝑘, o de 𝑥𝑒 . El área esta dado mediante la siguiente ecuación:
𝐴 = 4𝑥𝑒 𝑦𝑒 (159)
Pero 𝑥𝑒 = 𝑦𝑒
𝐴 = 4𝑥𝑒 2 (160)
El área de drene puede ser usada para calcular el área de drene de la pendiente y
el final del tiempo de 1⁄2 de la pendiente.
200.8107 𝑇 1 (161)
𝑚𝐿𝐺𝐶 = [ ]
√(𝜙𝜇𝑐𝑡 ) ℎ √𝑘 𝑥𝑒
Despejamos la permeabilidad
200.8107 𝑇 1 (162)
√𝑘 = [ ]
√(𝜙𝜇𝑐𝑡 ) ℎ 𝑚𝐿𝐺𝐶 𝑥𝑒
Sabemos que la distancia a la frontera externa para podemos obtener mediante la
siguiente expresión:
(163)
𝑘𝑡
𝑥𝑒 = 0.1125√
(𝜙𝜇𝑐𝑡 )𝑖
Reordenamos la ecuación anterior y sustituimos en la Ec. C-10.
𝑡 (C-13)
𝑥𝑒 = 0.1125√𝑘√
(𝜙𝜇𝑐𝑡 )𝑖
200.8107 𝑇 1 𝑡 (164)
𝑥𝑒 = 0.1125 [ ] √
√(𝜙𝜇𝑐𝑡 ) ℎ 𝑚𝐿𝐺𝐶 𝑥𝑒 (𝜙𝜇𝑐𝑡 )𝑖
Despejamos 𝑥𝑒 de la Ec. C-13:
194
200.8107 𝑇 1 𝑡 (165)
𝑥𝑒 𝑥𝑒 = 0.1125 [ ] √
√(𝜙𝜇𝑐𝑡 ) ℎ 𝑚𝐿𝐺𝐶 (𝜙𝜇𝑐𝑡 )𝑖
22.6 𝑇 √𝑡 (166)
𝑥𝑒 2 =
ℎ(𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐿𝐺𝐶
90.36 𝑇 √𝑡 (167)
𝐴=
ℎ(𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐿𝐺𝐶
Calculo del Volumen de poro,𝑉𝑝 .
90.36 𝑇 √𝑡 (169)
𝐴ℎ𝜙𝑖 =
(𝜇𝑐𝑡 )𝑖 𝑚𝐿𝐺𝐶
90.36 𝑇 √𝑡 (170)
𝑉𝑝 =
(𝜇𝑐𝑡 )𝑖 𝑚𝐿𝐺𝐶
90.36 𝑇 √𝑡 (172)
𝑉𝑝 = (1 − 𝑠𝑤𝑖 )
(𝜇𝑐𝑡 )𝑖 𝑚𝐿𝐺𝐶
𝐺=
𝐵𝑔𝑖
90.36 𝑇√𝑡(1 − 𝑠𝑤𝑖 ) (173)
𝐺=
(𝜇𝑐𝑡 )𝑖 𝐵𝑔𝑖 𝑚𝐿𝐺𝐶
195
Además, sabemos:
(1 − 𝑠𝑤𝑖 ) = 𝑠𝑔𝑖 (174)
Sustituimos la Ec. C-24 en la Ec. C-23 y reacomodamos.
90.36 𝑇𝑠𝑔𝑖 √𝑡 (175)
𝐺= ( )
(𝜇𝑐𝑡 𝐵𝑔 )𝑖 𝑚𝐿𝐺𝐶
Si la compresibilidad domina, el cálculo de gas se vuelve insensible para el valor
usado para 𝑠𝑤 .Por lo tanto tenemos
𝑐𝑡 = 𝑐𝑓 + 𝑠𝑔 𝑐𝑔 + 𝑠𝑤 𝑐𝑤 = 𝑠𝑔 𝑐𝑔 (176)
Sustituyendo en la Ec. C-25:
90.36 𝑇𝑠𝑔𝑖 √𝑡 (177)
𝐺= ( )
(𝜇 𝑠𝑔 𝑐𝑔 𝐵𝑔 )𝑖 𝑚𝐿𝐺𝐶
196
(C-31)
0.00633 𝑘𝑡
𝑥𝑒 = √ √
0.25 (𝜙𝜇𝑐𝑡 )𝑖
(C-32)
𝑘𝑡
𝑥𝑒 = 0.159122√
(𝜙𝜇𝑐𝑡 )𝑖
𝐴 = 4𝑥𝑓 𝑦𝑒 (C-34)
Pero consideramos que 𝑥𝑒 = 𝑥𝑓 = 𝑦𝑒 , por lo tanto
𝐴 = 4𝑥𝑒 2 (C-35)
𝑡 (C-38)
𝑥𝑒 = 0.159√𝑘√
(𝜙𝜇𝑐𝑡 )𝑖
315.43𝑇 1 𝑡 (C-39)
𝑥𝑒 = 0.159 √
√(𝜙𝜇𝑐𝑡 )𝑖 ℎ 𝑚𝐿𝑃𝐶 𝑥𝑒 (𝜙𝜇𝑐𝑡 )𝑖
Despejamos 𝑥𝑒 la Ec. C-39 y re arreglamos la ecuación
197
315.43𝑇 1 𝑡 (C-40)
𝑥𝑒 𝑥𝑒 = 0.159 √
(𝜙𝜇𝑐𝑡 )𝑖 ℎ 𝑚𝐿𝑃𝐶 (𝜙𝜇𝑐𝑡 )𝑖
50.153𝑇 1 (C-41)
𝑥𝑒 2 = √𝑡
(𝜙𝜇𝑐𝑡 )𝑖 ℎ 𝑚𝐿𝑃𝐶
200.612 𝑇 √𝑡 (C-45)
𝐴𝜙𝑖 ℎ =
(𝜇𝑐𝑡 )𝑖 𝑚𝐿𝑃𝐶
200.612 𝑇 √𝑡 (C-46)
𝑉𝑝 =
(𝜇𝑐𝑡 )𝑖 𝑚𝐿𝑃𝐶
198
200.612 𝑇 √𝑡 (C-48)
(1 − 𝑠𝑤𝑖 )
(𝜇𝑐𝑡 )𝑖 𝑚𝐿𝑃𝐶
𝐺=
𝐵𝑔𝑖
Reordenamos la Ec. C-48
𝑐2 = 𝑏𝐵𝐺𝐶 = 𝑋𝑠 𝑔 (C-54)
Donde X es:
1424𝑇 (C-
𝑋=
𝑘√𝐴𝑐 55179)
Remplazando X en la Ec. C-54
1424𝑇 (C-56)
𝑏𝐵𝐺𝐶 = 𝑠𝑔
𝑘√𝐴𝑐
199
1424𝑇 (C-57)
𝑏𝐵𝐺𝐶 = 𝑠1
𝑘√𝐴𝑐
Despejando 𝑠1 y el daño al comienzo del régimen de flujo puede ser calculado como:
𝑘√𝐴𝑐 (C-58)
𝑠1 = 𝑏𝐵𝐺𝐶
1424𝑇
2
0.00633𝑘𝑡𝑓𝑙𝑟 (C-
𝑟𝑖𝑛𝑣 = 61180)
(𝜙𝜇𝑐𝑡 )𝑖 0.20
(C-62)
0.00633 𝑘 𝑡𝑓𝑙𝑟
𝑟𝑖𝑛𝑣 =√
(𝜙𝜇𝑐𝑡 )𝑖 0.20
200
2
𝐴 = 𝜋 𝑟𝑖𝑛𝑣 (C-65)
𝑡𝑓𝑙𝑟 (C-69)
𝑟𝑖𝑛𝑣 = 0.1779√𝑘√
(𝜙𝜇𝑐𝑡 )𝑖
𝑡𝑓𝑙𝑟 (C-70)
1640 𝑇
𝑟𝑖𝑛𝑣 = 0.1779√ √
𝑚𝑅𝐺𝐶 ℎ (𝜙𝜇𝑐𝑡 )𝑖
Sustituimos la Ec. C-70 en la Ec. C-65
2 (C-71)
1640 𝑇 𝑡𝑓𝑙𝑟
𝐴 = 𝜋 (0.1779√ √ )
𝑚𝑅𝐺𝐶 ℎ (𝜙𝜇𝑐𝑡 )𝑖
Reacomodando términos
𝑇 𝑡𝑓𝑙𝑟 (C-72)
𝐴 = 𝜋 (51.9034 )
𝑚𝑅𝐺𝐶 ℎ (𝜙𝜇𝑐𝑡 )𝑖
Quedando finalmente:
163 𝑇 𝑡𝑓𝑙𝑟 (C-74)
𝐴=
𝑚𝑅𝐺𝐶 ℎ (𝜙𝜇𝑐𝑡 )𝑖
201
𝑉𝑝 = 𝐴(𝜙)𝑖 ℎ (C-75)
Sustituyendo la Ec. C-74 en la Ec. 75, queda:
163 𝑇 𝑡𝑓𝑙𝑟 (C-76)
𝑉𝑝 = (𝜙)𝑖 ℎ
𝑚𝑅𝐺𝐶 ℎ (𝜙𝜇𝑐𝑡 )𝑖
Eliminando 𝜙𝑖 𝑦 ℎ
163 𝑇 𝑡𝑓𝑙𝑟 (C-77)
𝑉𝑝 = ( )
(𝜇𝑐𝑡 )𝑖 𝑚𝑅𝐺𝐶
Tenemos lo siguiente:
(1 − 𝑠𝑤𝑖 ) = 𝑠𝑔𝑖 (C-81)
Sustituimos la Ec. C-81 en la Ec. C-80
163 𝑇 𝑡𝑓𝑙𝑟 𝑠𝑔𝑖 (C-82)
𝐺=
(𝜇𝑐𝑡 )𝑖 𝐵𝑔𝑖 𝑚𝑅𝐺𝐶
Reacomodando términos, finalmente queda:
163 𝑇 𝑠𝑔𝑖 𝑡𝑓𝑙𝑟 (C-83)
𝐺= ( )
(𝜇𝑐𝑡 𝐵𝑔 )𝑖 𝑚𝑅𝐺𝐶
202
Ecuaciones de interpretación para régimen de flujo radial en un modelo
homogéneo infinito a condiciones de presión de fondo fluyendo contante,𝑝𝑤𝑓 .
2
Sabemos que 𝑡 = 𝑡𝑓𝑙𝑟 , 𝑡𝐷𝐿 = 𝑡𝐷𝑟𝑖𝑛𝑣 , 𝐿2 = 𝑟𝑖𝑛𝑣 y 𝑡𝐷𝑟𝑖𝑛𝑣 = 0.20, por lo tanto:
0.00633𝑘𝑡𝑓𝑙𝑟 (C-85)
𝑡𝐷𝐿 = 2 = 0.20
(𝜙𝜇𝑐𝑡 )𝑖 𝑟𝑖𝑛𝑣
2
0.00633𝑘𝑡𝑓𝑙𝑟 (C-86)
𝑟𝑖𝑛𝑣 =
(𝜙𝜇𝑐𝑡 )𝑖 0.20
(C-87)
0.00633𝑘𝑡𝑓𝑙𝑟
𝑟𝑖𝑛𝑣 =√
(𝜙𝜇𝑐𝑡 )𝑖 0.20
1 (C-88)
0.0063 2 𝑘𝑡𝑓𝑙𝑟
𝑟𝑖𝑛𝑣 =( ) √
0.20 (𝜙𝜇𝑐𝑡 )𝑖
203
1640 𝑇 (C-91)
𝑘ℎ=
𝑚𝑅𝑃𝐶
1640 𝑇 (C-
𝑘 =
ℎ 𝑚𝑅𝑃𝐶 92181)
De la Ec. C-89 reordenamos términos
𝑡𝑓𝑙𝑟 (C-93)
𝑟𝑖𝑛𝑣 = 0.1779√𝑘√
(𝜙𝜇𝑐𝑡 )𝑖
𝑡𝑓𝑙𝑟 (C-
1640 𝑇
𝑟𝑖𝑛𝑣 = 0.1779√ √ 94182)
ℎ 𝑚𝑅𝑃𝐶 (𝜙𝜇𝑐𝑡 )𝑖
Reordenando términos
𝑇 𝑡𝑓𝑙𝑟 (C-96)
𝐴 = 𝜋 (51.9034 )
𝑚𝑅𝑃𝐶 ℎ (𝜙𝜇𝑐𝑡 )𝑖
Finalmente obtenemos:
163.0593 𝑇 𝑡𝑓𝑙𝑟 (C-
𝐴= 97184)
(𝜙𝜇𝑐𝑡 )𝑖 ℎ 𝑚𝑅𝑃𝐶
204
Reordenamos y eliminamos términos
163 𝑇 𝑡𝑓𝑙𝑟 (C-100)
𝑉𝑝 =
(𝜇𝑐𝑡 )𝑖 𝑚𝑅𝑃𝐶 𝑖
Reordenando términos:
163 𝑇 𝑠𝑔𝑖 𝑡𝑓𝑙𝑟 (C-
𝐺= 106187)
(𝜇𝑐𝑡 𝐵𝑔 )𝑖 𝐵𝑔𝑖 𝑚𝑅𝑃𝐶𝑖
205
Ecuaciones de interpretación para flujo esférico en un modelo homogéneo
infinito a condiciones de explotación a gasto contante constante,𝑞𝑔 .
Calculo de la permeabilidad,𝑘.
𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄ 1
En una grafica 𝑞𝑔 vs √𝑡 obtenemos una línea recta con la siguiente
pendiente:
2⁄
3 (C-109)
10098 𝑇 √(𝜙𝜇𝑐𝑡 )𝑖
𝑘=[ ]
𝑚𝐸𝐺𝐶
206
Ecuaciones de interpretación para régimen de flujo dominado por la frontera
externa en un yacimiento homogéneo infinito a las condiciones a gasto
constante de gas,𝒒𝒈 , en estado pseudo-estacionario (PSS).
207
2 𝑏𝐹𝐺𝐶 𝑘ℎ 𝐴 2.2458 (C-120)
exp ( − 2𝑠) = 2
1424 𝑇 𝑟𝑤 𝐶𝐴
2.2458 𝐴 (C-121)
𝐶𝐴 =
2 𝑏𝐹𝐺𝐶 𝑘ℎ 2
exp ( 1424 𝑇 − 2𝑠) 𝑟𝑤
Calculo del volumen de poro del yacimiento,𝑉𝑝 , y del área de drene del pozo,𝐴.
57 𝑇 1 (1 − 𝑠𝑤𝑖 ) (C-129)
𝐺= ( )
(𝜇𝑐𝑡 )𝑖 𝑚𝐹𝐺𝐶 𝐵𝑔𝑖
208
Sabemos que:
(1 − 𝑆𝑤𝑖 ) = 𝑆𝑔𝑖 (C-130)
Por lo tanto:
57 𝑇 1 𝑠𝑔𝑖 (C-131)
𝐺= ( )
(𝜇𝑐𝑡 )𝑖 𝑚𝐹𝐺𝐶 𝐵𝑔𝑖
Reacomodando términos, obtenemos:
57 𝑇 𝑠𝑔𝑖 1 (C-132)
𝐺= ( )
(𝜇𝑐𝑡 𝐵𝑔 )𝑖 𝑚𝐹𝐺𝐶
4 𝜋(0.00633)𝑘
𝐴=
2.2458 𝐴
2.303 𝑙𝑛 (𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐹𝑃𝐶
2.2458 𝐴
𝑟𝑤2
𝑏𝐹𝑃𝐶 𝑘 ℎ
exp [2 ( 1424𝑇 )] 𝑟𝑤2
( ( ) )
209
(C-137)
4 𝜋(0.00633)𝑘
𝐴=( )
exp( 2𝑏𝐹𝑃𝐶 𝑘ℎ)
2.303 𝑙𝑛 ( ) (𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐹𝑃𝐶
1424 𝑇
(C-138)
4 𝜋(0.00633)𝑘
𝐴=( )
2.303 (𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐹𝑃𝐶 2𝑏𝐹𝑃𝐶 𝑘ℎ
1424 𝑇
4 𝜋(0.00633) (1424)𝑇 (C-139)
𝐴=( )
2.303 (𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐹𝑃𝐶 2𝑏𝐹𝑃𝐶 ℎ
24.5923 𝑇 (C-140)
𝐴=( )
(𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐹𝑃𝐶 𝑏𝐹𝑃𝐶 ℎ
Finalmente obtenemos:
56 𝑇 1 −2 𝑏𝐹𝑃𝐶 𝑘 ℎ (C-144)
𝐶𝐴 = 2
( ) exp [( )]
ℎ (𝜙𝜇𝑐𝑡 )𝑖 𝑏𝐹𝑃𝐶 𝑟𝑤 𝑚𝐹𝑃𝐶 1424𝑇
Calculo del volumen del yacimiento, VY.
Se tiene el producto de:
𝐴ℎ = 𝑉𝑌 (C-145)
Tenemos la Ec. C-141 y se sustituye en la Ec. C-145, se obtiene:
25 𝑇 1 (C-146)
𝐴ℎ = 𝑉𝑌 = ( )ℎ
(𝜙𝜇𝑐𝑡 )𝑖 ℎ 𝑏𝐹𝑃𝐶 𝑚𝐹𝑃𝐶
O
25 𝑇 1 (C-147)
𝑉𝑌 = ( )
(𝜙𝜇𝑐𝑡 )𝑖 𝑏𝐹𝑃𝐶 𝑚𝐹𝑃𝐶
210
Calculo del Volumen de poro del yacimiento.
EL volumen poro está representando por la siguiente ecuación:
𝑉𝑝 = 𝐴(𝜙)𝑖 ℎ (C-148)
Sustituimos la Ec. C-141 en la Ec. C-148
25 𝑇 1 (C-149)
𝑉𝑝 = ( )( ) (𝜙)𝑖 ℎ
(𝜙𝜇𝑐𝑡 )𝑖 𝑏𝐹𝑃𝐶 ℎ 𝑚𝐹𝑃𝐶
25 𝑇 1 (C-150)
𝑉𝑝 = ( )( )
(𝜇𝑐𝑡 )𝑖 𝑏𝐹𝑃𝐶 𝑚𝐹𝑃𝐶
Calculo del volumen original de gas, G.
Se tiene la siguiente ecuación:
𝑉𝑝 (1 − 𝑠𝑤𝑖 ) (C-151)
𝐺=
𝐵𝑔𝑖
25 𝑇 1 (1 − 𝑠𝑤𝑖 ) (C-153)
𝐺=( )( )
(𝜇𝑐𝑡 )𝑖 𝑏𝐹𝑃𝐶 𝑚𝐹𝑃𝐶 𝐵𝑔𝑖
Sabemos que:
(1 − 𝑠𝑤𝑖 ) = 𝑠𝑔𝑖 (C-154)
Por lo tanto:
25 𝑇 1 𝑠𝑔𝑖 (C-155)
𝐺=( )( )
(𝜇𝑐𝑡 )𝑖 𝑏𝐹𝑃𝐶 𝑚𝐹𝑃𝐶 𝐵𝑔𝑖
Reacomodando términos, finalmente obtenemos:
25 𝑇𝑠𝑔𝑖 1 (C-156)
𝐺= ( )
(𝜇𝑐𝑡 𝐵𝑔 )𝑖 𝑏𝐹𝑃𝐶 𝑚𝐹𝑃𝐶
211
3.1.1 Régimen de flujo lineal en yacimiento de tight gas
El flujo lineal se produce cuando las trayectorias de flujo son paralelas y el fluido
fluye en una sola dirección. (McKinney, 2005) En Fig. 12-1 podemos observar una
representación ideal de un sistema de flujo lineal. Este tipo de flujo se presenta en
un gran número de pozos de baja permeabilidad. Algunas veces el flujo lineal
transitorio dura muchos años, debido a la extremadamente baja permeabilidad de
los yacimientos. El análisis de producción de flujo lineal a corto plazo se puede
caracterizar con tratamiento de fracturas, pero el flujo lineal a largo plazo puede ser
controlado en algunos casos por la geometría del yacimiento. En algunos casos el
flujo lineal de largo plazo en yacimientos de tight gas se puede presentar en
yacimientos lineales o largos, vetas de alta permeabilidad entre dos pozos sin
límites de flujo, comportamiento transitorio de doble porosidad para yacimientos
radiales, pozos interceptados por fracturas verticales y pozos horizontales con
fracturas. (J.A. Arévalo-Villagrán, 2001)
Al final del flujo lineal, la respuesta de la presión de los pozos comienza a disminuir
como el flujo entra desde afuera de las puntas de fracturas, lo que indica un efecto
de la frontera externa. (Gringarten, 1975) Para estos pozos, parece ser flujo lineal
durante flujo transitorio y dominado por la frontera externa.
212
Ilustración 0-2 Flujo lineal ideal dentro de las fracturas. (Arévalo J. D., 2009)
En varios pozos se ha observado la presencia de flujo lineal a largo plazo que dura
meses y hasta años. Una razón para el flujo lineal a largo plazo es debido a que
esos pozos han sido perforados y completados en yacimientos de gas de baja
permeabilidad con apropiado fracturamiento hidráulico, el cual la mitad de la longitud
de la fractura,𝑥𝑓 ,puede continuar y extenderse a los limites laterales en la dirección
𝑥,𝑥𝑒 ,lo que significa que 𝑥𝑓 = 𝑥𝑒 .
213
Well
2xe
L
Ilustración 0-3 Sistema lineal Ilustración 0-4 Sistema lineal con matriz
homogéneo. (Arévalo Villagran, Analysis paralela. (Arévalo Villagran, Analysis of
of Long-term Behavior in tight gas Long-term Behavior in tight gas
resevoirs:case histories, 2001) resevoirs:case histories, 2001)
Asumimos que el pozo se encuentra en el centro del are de drene rectangular y las
fracturas tienen una conductividad infinita. La distancia a la frontera en la dirección
perpendicular a la fractura es 𝑦𝑒 .El área de drene de la geometría del pozo de la
Fig. 12.5 es 4𝑥𝑒 𝑦𝑒 .
214
xf
LINEAR
FLOW ye
xe
Sin embargo, la geometría de flujo de la Fig 3.3 es más simple que la mayoría de
los casos de los pozos fracturados en la literatura petrolera. Generalmente, los
casos incluyen yacimientos infinitos, lo que significa que se acaba el flujo lineal y se
convierte en flujo pseudo- radial. Sin embargo, la fractura es tratada como si tuviera
una conductividad limitada.
Wattenbarger ha adaptado las soluciones de Miller and Nabor and Barham para
pozos fracturados en la geometría de la Fig 3.3. Incluyen ambas el gasto
contante,𝑞𝑔 , y la presión de fondo fluyendo constante,𝑝𝑤𝑓 ,para flujo lineal en
yacimientos rectangulares.
Dichos autores contribuyeron con las siguientes ecuaciones que nos ayudan al
análisis del comportamiento de flujo lineal a largo plazo en yacimientos de gas con
baja permeabilidad y con valores adimensionales de conductividad de la fractura,
𝐹𝑐𝑑 > 50.
Las Ec. C-157 y Ec. C-158 se aplican para flujo lineal dentro de la fractura. Además,
fueron adaptadas para yacimientos de gas con muy buena exactitud.
1 𝜋 (C-158)
= √𝜋𝑡𝐷𝐿
𝑞𝐷 2
215
Sustituimos la pseudo presión real para el gas,𝑚𝐷 , por la presión adimensional para
líquido,𝑝𝐷 , para utilizarlo en casos de comportamiento lineal para gasto de flujo
contante de gas,𝑞𝑔 , y a presión de fondo fluyendo contante, 𝑝𝑤𝑓 ,respectivamente,
en el comienzo de la explotación.
1 𝜋 (C-159)
= √𝜋𝑡𝐷𝐿
𝑚𝑤𝐷𝐿 2
Debemos notar que las soluciones para un yacimiento lineal infinito a gasto contante
𝜋
y presión contante solo difiere por un factor de 2 .
3.1.3 Modelo para régimen de Flujo lineal en un yacimiento Homogéneo infinito con
una producción a gasto constante de gas,𝒒𝒈 .
En la solución de líquidos puede ser adaptada para flujo de gas usando la pseudo-
presión real. Remplazamos 𝑝𝑤𝐷 por 𝑚𝑤𝐷 , el cual se defino como:
216
Sabemos que:
𝑚𝐷𝐿 = 𝑚𝑤𝐷𝐿 + 𝑆 𝑔 (C-164)
Sustituimos la Ec. C-160 y Ec. C-161 en la Ec. C-164 para poder llegar al modelo
de la ecuación de la línea recta 𝑦 = 𝑚𝑥 + 𝑏 => ∆𝑝 = 𝑐1 √𝑡 + 𝑐2 , la cual es el modelo
tipo curva para flujo lineal.
Δ𝑚(𝑝)⁄
Despejamos 𝑞𝑔
(C-167)
Δ𝑚(𝑝) 1424𝑇 𝑘𝑡
=[ ] [√𝜋 0.00633√ + 𝑆𝑔 ]
𝑞𝑔 𝑘ℎ (𝜙𝜇𝑐𝑡 )𝑖 𝐿2
Consideramos lo siguiente:
1424𝑇 (C-168)
𝑋=
𝑘ℎ
217
200.8107 𝑇 1 (C-170)
𝑐1 = 𝑚𝐿𝐺𝐶 = [ ]
√(𝜙𝜇𝑐𝑡 ) ℎ √𝑘 𝐿
Ilustración 0-6 Grafica normal del modelo tipo curva para régimen de flujo lineal con
intersección en cero y no intersección en cero de la línea recta, bajo la condición de
producción a gasto constante de gas,𝒒𝒈. (Arévalo Villagran, Analysis of Long-term
Behavior in tight gas resevoirs:case histories, 2001)
𝐴 (C-173)
𝑋𝑒 ℎ =
4
218
803.2427 𝑇 1 (C-175)
√𝑘𝐴 = [ ]
√(𝜙𝜇𝑐𝑡 ) 𝑚𝐿𝐺𝐶
1 𝜋 (C-177)
= √𝜋𝑡𝐷𝐿
𝑚𝑤𝐷𝐿 2
(C-181)
𝑘ℎ[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 𝜋 0.00633𝑘𝑡
= √𝜋 + 𝑆𝑔
1424𝑞𝑔 𝑇 2 (𝜙𝜇𝑐𝑡 )𝑖 𝑥𝑒2
Δ𝑚(𝑝)
Despejamos 𝑞𝑔
219
(C-182)
Δ𝑚(𝑝) 𝜋 𝑇 𝑘𝑡 1424 𝑇 𝑔
= 1424 √𝜋 0.00633 √ 2
+ 𝑆
𝑞𝑔 2 𝑘ℎ (𝜙𝜇𝑐𝑡 )𝑖 𝑥𝑒 𝑘ℎ
Hacemos lo siguiente:
1424 𝑇 (C-183)
𝑋=
𝑘ℎ
Haciendo las operaciones y reordenando la Ec. C-182:
Δ𝑚(𝑝) 315.43 𝑇 (C-184)
= √𝑡 + 𝑋 𝑆 𝑔
𝑞𝑔 √(𝜙𝜇𝑐𝑡 )𝑖 √𝑘ℎ𝑥𝑒
La Ec. C-184 nos representa el modelo de tipo curva para flujo lineal. Que podemos
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )]⁄
encontrar en una gráfica normal de 𝑞𝑔 vs √𝑡, donde podemos
encontrar la pendiente mediante:
315.43 𝑇 (C-185)
𝐶1 = 𝑚𝐿𝑃𝐶 =
√(𝜙𝜇𝑐𝑡 )𝑖 √𝑘ℎ𝑥𝑒
Ilustración 0-7 Grafica normal del modelo de tipo curva para régimen de flujo lineal
con intersección en cero y no intersección en cero de la línea recta, bajo la condición
de presión de fondo fluyendo constante,𝒑𝒘𝒇 . (Arévalo Villagran, Analysis of Long-term
Behavior in tight gas resevoirs:case histories, 2001)
La Ec. C-185 representa la pendiente de la línea recta con base en la distancia del
pozo a la frontera externa y la fractura hidráulica. Tomamos en cuenta que 𝑦𝑒 = 𝑥𝑒 =
𝑥𝑓 , por lo tanto:
220
315.43 𝑇 (C-186)
𝑚𝐿𝑃𝐶 =
√(𝜙𝜇𝑐𝑡 )𝑖 ℎ √𝑘𝑥𝑒
Despejando √𝑘𝑋𝑒 de la Ec. C-186:
315.43 𝑇 (C-187)
√𝑘𝑥𝑒 =
√(𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐿𝑃𝐶 ℎ
Sabemos que el área se calcula mediante la siguiente ecuación
𝐴𝑐 = 4𝑥𝑒 ℎ (C-188)
𝐴𝑐 (C-189)
𝑥𝑒 ℎ =
4
𝐴𝑐 315.43 𝑇 (C-191)
√𝑘 =
4 √(𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐿𝑃𝐶
1261.72 𝑇 (C-192)
√𝑘𝐴𝑐 =
√(𝜙𝜇𝑐𝑡 )𝑖 𝑚𝐿𝑃𝐶
221
Ilustración 0-8 Modelo de flujo Bilineal. (Arévalo J. D., 2009)
Un flujo bilineal existe cuando la mayoría del fluido que ingresa al pozo proviene de
la formación y los efectos de la punta de la fractura no afectan aun el
comportamiento del pozo. En una gráfica log-log de 𝑃𝑤𝐷 vs 𝑡𝐷𝑏 , el comportamiento
de la presión del flujo bilineal exhibe una línea recta cuya pendiente es igual a 1⁄4.
La presión adimensional del pozo para un periodo de flujo bilineal esta dado por:
𝜋 1⁄ (C-195)
𝑃𝑊𝐷𝐵 = 𝑡𝐷𝐵4 + 𝑠 𝑔
Γ(1.25)√2(𝑘𝑓 𝑏𝑓 )𝐷
o
2.45 1⁄ (C-196)
𝑃𝑊𝐷𝐵 = 𝑡𝐷𝐵4 + 𝑠 𝑔
√(𝑘𝑓 𝑏𝑓 )𝐷
222
Esta ecuación indica que la gráfica de 𝑃𝑊𝐷𝐵 𝑣𝑠 4√𝑡𝐷𝐵 genera una línea recta cuya
pendiente es 2.45 , interceptando al origen. Y modela la caída de presión
⁄√(𝑘 𝑏 )
𝑓 𝑓 𝐷
−1.53 (C-198)
𝑡𝐷𝐵𝐹 ≃ 0.025 [(𝑘𝑓 𝑏𝑓 ) − 1.5] 𝑝𝑎𝑟𝑎 1.6 ≤ (𝑘𝑓 𝑏𝑓 ) ≤ 3
𝐷 𝐷
223
−4 (C-199)
4.55
𝑡𝐷𝐵𝐹 ≃ − 2.25 𝑝𝑎𝑟𝑎 (𝑘𝑓 𝑏𝑓 ) ≤ 1.6
𝐷
√(𝑘𝑓 𝑏𝑓 )𝐷
[ ]
224
Para poder llegar al modelo de línea recta, el cual es el modelo de tipo curva para
flujo bilineal. Hacemos la sustitución de la Ec.C-200 y Ec. C-201 en la Ec. C-202.
𝑘ℎ[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 2.45 1⁄ (C-203)
= 𝑡𝐷𝐵4 + 𝑠 𝑔
1424𝑞𝑔 𝑇
√(𝑘𝑓 𝑏𝑓 )𝐷
El tiempo adimensional se representa como en la Ec. C-204.
0.00633𝑘𝑡 (C-204)
𝑡𝐷𝐿 =
(𝜙𝜇𝑐𝑡 )𝑖 𝐿2
Reacomodando términos
1⁄ 1⁄
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] (1424) ∗ 2.45 ∗ 0.00633 4 𝑇 𝑘𝑡 4 1424𝑇 (C-209)
= + 𝑠𝑔
𝑞𝑔 𝑘𝑓 𝑤𝑓 (𝜙𝜇𝑐𝑡 )𝑖 𝐴𝑐 𝑘√𝐴𝑐
𝑘√𝐴𝑐 √
𝑘√𝐴𝑐
Para mayor comodidad hacemos
1424𝑇
𝑋= (C-210)
𝑘√𝐴𝑐
225
1⁄
1⁄ 1⁄ 2 (C-211)
2
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 984.072 𝑇 𝑘 2 (𝐴𝑐 ) 1 1
𝑘 ⁄4 𝑡 ⁄4
= 1⁄ 1⁄ 1⁄
+ 𝑋𝑠 𝑔
𝑞𝑔 √𝑘𝑓 𝑤 (𝜙𝜇𝑐𝑡 )𝑖 4 𝐴𝑐 4
𝑘 𝐴𝑐 2
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 984.072 𝑇 4
(C-212)
= 1 1 √ 𝑡 + 𝑋𝑠 𝑔
𝑞𝑔 1 ⁄ ⁄
𝑘 ⁄4 𝐴𝑐 2 √𝑘𝑓 𝑤 (𝜙𝜇𝑐𝑡 )𝑖 4
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 984.072 𝑇 4 (C-213)
= 1 √ 𝑡 + 𝑋𝑠 𝑔
𝑞𝑔 ⁄
√𝑘(𝜙𝜇𝑐𝑡 )𝑖 𝐴𝑐 2 √𝑘𝑓 𝑤
4
La ecuación anterior es el modelo de tipo curva para flujo bilineal. Donde nos indica
que el cambio de presión es inversamente proporcional a la raíz de la conductividad
de la fractura y directamente proporcional a la raíz cuarta del tiempo. En una grafica
𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄ 4
𝑞𝑔 𝑣𝑠 √𝑡 podemos encontrar una línea recta que pasa por el origen,
donde la pendiente es dada por:
984.072 𝑇 (C-214)
𝑐1 = 𝑚𝐵𝐺𝐶 = 1⁄
√𝑘(𝜙𝜇𝑐𝑡 )𝑖 𝐴𝑐 2 √𝑘𝑓 𝑤
4
Despejamos el termino√𝑘𝑓 𝑤
226
984.072 𝑇 (C-215)
√𝑘𝑓 𝑤 = 1⁄ 𝑚
4
√𝑘(𝜙𝜇𝑐𝑡 )𝑖 𝐴𝑐 2 𝐵𝐺𝐶
[Link] Modelo para flujo radial en un modelo homogéneo infinito a las condiciones
de explotación a gasto constante de gas,𝒒𝒈 .
Si hacemos
227
𝑚𝐷𝐿 = 𝑚𝑊𝐷 (C-219)
Sustituimos la Ec. C-217 y la Ec. C-218 en la Ec. C-219
228
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] (C-227)
𝑞𝑔
1640 𝑇 1640 𝑇 0.00633𝑘
= log 𝑡 + (log + 0.3516 + 0.8686𝑠)
𝑘ℎ 𝑘ℎ (𝜙𝜇𝑐𝑡 )𝑖 𝑟𝑤2
La Ec. 129 es similar a la ecuación de línea recta𝑦 = 𝑚𝑥 + 𝑏 = 𝑐1 log 𝑡 + 𝑐2 . En una
𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄
gráfica de 𝑞𝑔 vs tiempo la pendiente está dada por:
1640 𝑇 (C-228)
𝑚𝑅𝐺𝐶 =
𝑘ℎ
Ilustración 0-12 Grafica semi- log del modelo tipo curva para régimen radial con no
intersección en cero de la línea recta, bajo la condición de producción a gasto de
flujo constante de gas,𝒒𝒈. (Arévalo Villagran, Analysis of Long-term Behavior in tight gas
resevoirs:case histories, 2001)
Despejando 𝑘 ℎ
1640 𝑇 (C-229)
𝑘ℎ=
𝑚𝑅𝐺𝐶
229
1 𝑏𝑘ℎ 0.0063𝑘 (C-231)
𝑠= ( − log − 0.352
0.8686 1640𝑇 (𝜙𝜇𝑐𝑡 )𝑖 𝑟𝑤2
Sabemos que la pseudo- presión real está dado por la Ec. C-233 y el tiempo
adimensional por las Ec. C-234
0.00633𝑘𝑡 (C-234)
𝑡𝐷𝐿 =
(𝜙𝜇𝑐𝑡 )𝑖 𝐿2
Se hace:
1 (C-235)
𝑚𝐷𝐿 =
𝑚𝑊𝐷
230
𝑘ℎ[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 1 0.00633𝑘𝑡 (C-238)
= (2.3025850929 log + 0.80907 + 2𝑠)
1424𝑞𝑔 𝑇 2 (𝜙𝜇𝑐𝑡 )𝑖 𝑟𝑤2
Simplificamos términos
𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄
Despejamos 𝑞𝑔
231
Ilustración 0-13 Grafica semi- log del modelo tipo curva para régimen radial con no
intersección en cero de la línea recta, bajo la condición de producción a presión de
fondo fluyendo constante,𝒑𝒘𝒇 . (Arévalo Villagran, Analysis of Long-term Behavior in tight
gas resevoirs:case histories, 2001)
Despejando 𝑘 ℎ
1640 𝑇 (C-244)
𝑘ℎ=
𝑚𝑅𝑃𝐶
232
1. Flujo dentro de un hundimiento esférico ficticio.
Sabemos que:
𝑟𝑒 (C-248)
𝑟𝐷 =
𝑟𝑤
La solución de líquidos puede ser adaptada para flujo de gas usando la pseudo-
presión real. Remplazamos 𝑝𝑊𝐷𝐸 por 𝑚𝑊𝐷𝐸 , el cual se define como:
1 (C-250)
𝑚𝑊𝐷𝐸 = 1 −
√𝜋𝑡𝐷𝐿
233
Se hace:
𝑚𝐷𝐿 = 𝑚𝑊𝐷𝐸 (C-251)
1
Para poder llevar al modelo de la ecuación de línea recta 𝑦 = 𝑚𝑥 + 𝑏 = 𝑐1 + 𝑐2, el
√𝑡
cual es el modelo de tipo curva para flujo esférico. Sustituimos las Ec.C-160 y la Ec.
C-247 en la Ec. C-251
𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄
Despejamos la caída de la pseudo- presión, 𝑞𝑔 y reacomodamos
términos:
1
La Ec. C-256 es similar a la ecuación de la línea recta 𝑦 = 𝑚𝑥 + 𝑏 = 𝑐1 + 𝑐2 , la
√𝑡
cual es el modelo de tipo curva para flujo esférico.
[Link] Calculo de la permeabilidad,𝑘.
𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄ 1
En una grafica 𝑞𝑔 vs √𝑡 obtenemos una línea recta con la siguiente
pendiente:
234
10097.95 𝑇 √(𝜙𝜇𝑐𝑡 )𝑖 (C-257)
𝑐1 = 𝑚𝐸𝐺𝐶 = 3⁄
𝑘 2
2⁄
3 (C-259)
10098 𝑇 √(𝜙𝜇𝑐𝑡 )𝑖
𝑘=[ ]
𝑚𝐸𝐺𝐶
235
el flujo se estabiliza. Este periodo o régimen de flujo dominado por la frontera
externa inicial al final del periodo transitorio cuando la condición de no-flujo más
alejada de la pared del pozo es alcanzada por el disturbio de la presión y el área
total de drene comienza a contribuir en la producción, es decir, el área de drene ha
alcanzado las fronteras externas de no-flujo. Así mismo, las condiciones en la
cercanía del pozo, (gasto de flujo y presión) tienden a estabilizarse durante este
régimen de flujo dominado por la frontera.
Las fronteras se generan alrededor del pozo como resultado del gasto de flujo de
producción y variación de las propiedades de la formación (permeabilidad, espesor
de la zona productora, etc.). El ritmo de declinación de la presión en el yacimiento
dependerá de que tan rápido los fluidos sean producidos, la expansión de los fluidos
en el yacimiento y compactación del volumen de poros, pero el efecto más
significativo es el deterioro del comportamiento de afluencia, que se ve reflejada en
la declinación de la presión media del yacimiento y el incremento de la resistencia
de flujo, él se puede evaluar mediante el balance de materia. Las Fig.12-15 y 12.15
muestran dos casos de abatimiento a las condiciones de producción a gasto de flujo
constante (estado pseudo-estacionario PSS) y a presión de fondo fluyendo
constante (declinación exponencial), respectivamente.
236
Ilustración 0-16Dibujos de las gráficas en las cuales se detecta un régimen de flujo
dominado por la frontera externa bajo la condición de producción a presión de fondo
fluyendo constante (declinación exponencial). (Arévalo Villagran, Analysis of Long-term
Behavior in tight gas resevoirs:case histories, 2001)
Las Ecs. C-260 y C-261 pueden ser adaptadas para el flujo de gas mediante la
implementación de la caída pseudo- presión del gas real. Remplazamos la solución
de 𝑝𝑊𝐷𝐹𝐺𝐶 por 𝑚𝑊𝐷𝐹𝐺𝐶 , con la condición de explotación a gasto constante de
gas, 𝑞𝑔 ,:
1 𝐴 1 2.2458 (C-262)
𝑝𝑊𝐷𝐹𝐺𝐶 = 𝑚𝑊𝐷𝐹𝐺𝐶 = 2𝜋𝑡𝐷𝐴 + ln ( 2 ) + ln ( )+𝑠
2 𝑟𝑤 2 𝐶𝐴
237
Para la condición de explotación a presión de fondo fluyendo constante, 𝑝𝑤𝑓 :
2 4𝜋 1 (C-263)
𝑞𝑊𝐷𝐹𝑃𝐶 = 𝑚𝑊𝐷𝐹𝐺𝐶 = exp (− 𝑡𝐷𝐴 ) +
ln 𝑐𝑓𝑝𝑐 ln 𝑐𝑓𝑝𝑐 2𝑠
Sabemos que la solución analítica para el caso de flujo de gas real para la condición
de explotación a gasto constante de gas, 𝑞𝑔 , es dado por la Ec. C-262:
1 𝐴 1 2.2458 (C-264)
𝑚𝑊𝐷𝐹𝐺𝐶 = 2𝜋𝑡𝐷𝐴 + ln ( 2 ) + ln ( )+𝑠
2 𝑟𝑤 2 𝐶𝐴
Ahora hacemos:
𝑚𝐷𝐿 = 𝑚𝑊𝐷𝐹𝐺𝐶 (C-265)
Sustituimos las Ecs. C-160 y C-264 en la Ec. C-265 para poder llegar el modelo de
la ecuación de línea recta 𝑦 = 𝑚𝑥 + 𝑏 ⟶ ∆𝑝 = 𝑐1 𝑡 + 𝑐2, el cual es el modelo de tipo
curva para flujo dominado por la frontera externa a condiciones de explotación a
gasto de gas constante, 𝑞𝑔 :
238
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] (C-269)
𝑞𝑔
1424 𝑇 0.00633𝑘𝑡 1 𝐴 1 2.2458
= (2𝜋 ( ) + ln ( 2 ) + ln ( ) + 𝑠)
𝑘ℎ (𝜙𝜇𝑐𝑡 )𝑖 𝐴 2 𝑟𝑤 2 𝐶𝐴
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] (C-270)
𝑞𝑔
2 𝜋 (1424) (0.00633)𝑇𝑡
=( )
(𝜙𝜇𝑐𝑡 )𝑖 𝐴 ℎ
1424 𝑇 1 𝐴 1 2.2458
+ [ ln ( 2 ) + ln ( ) + 𝑠]
𝑘ℎ 2 𝑟𝑤 2 𝐶𝐴
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 56.6361 𝑇 1424 𝑇 𝐴 2.2458 (C-271)
= 𝑡+ [ln ( 2 ) +2𝑠]
𝑞𝑔 (𝜙𝜇𝑐𝑡 )𝑖 𝐴 ℎ 2 𝑘ℎ 𝑟𝑤 𝐶𝐴
Ilustración 0-17 Grafica normal del modelo de tipo curva para régimen de flujo
dominado por la frontera externa con no intersección en cero de la línea recta, bajo
239
la condición de producción a gasto de flujo constante,𝒒𝒈. (Arévalo Villagran, Analysis
of Long-term Behavior in tight gas resevoirs:case histories, 2001)
2 4𝜋 1 (C-273)
𝑞𝑊𝐷𝐹𝑃𝐶 = exp (− 𝑡𝐷𝐴 ) +
ln 𝑐𝑓𝑝𝑐 ln 𝑐𝑓𝑝𝑐 2𝑠
La Ec. 243 es la pseudo- presión del gas real adimensional con la condición de
explotación a presión de fondo fluyendo constante, 𝑝𝑤𝑓 , para flujo dominado por la
frontera externa.
La solución analítica para el caso de flujo de gas real para la condición de
explotación a presión de fondo fluyendo constante, 𝑝𝑤𝑓 , es análoga a la solución
analítica para el caso de líquido con la misma condición:
2 4𝜋 1 (C-276)
𝑚𝑊𝐷𝐹𝐺𝐶 = exp (− 𝑡𝐷𝐴 ) +
ln 𝑐𝑓𝑝𝑐 ln 𝑐𝑓𝑝𝑐 2𝑠
En donde:
240
2.2458 𝐴 (C-277)
𝑐𝑓𝑝𝑐 =
𝑟𝑤2 𝐶𝐴
Reordenando términos:
241
1424𝑇 (C-285)
𝑋=( )
𝑘ℎ
Se aplica logaritmo natural en ambos lados de la Ec. C-285:
242
1 [𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] (C-291)
( ) 2.303 𝑙𝑜𝑔 ( )
2.303 𝑞𝑔
1 4 𝜋0.00633𝑘
=( )( )𝑡
2.303 𝑙𝑛 𝑐𝑓𝑝𝑐 (𝜙𝜇𝑐𝑡 )𝑖 𝐴
1 1424𝑇 𝑙𝑛 𝑐𝑓𝑝𝑐
+( ) 2.303 log ( )
2.303 𝑘ℎ
1
+( ) 2.303 𝑙𝑜𝑔 (2𝑠 𝑋)
2.303
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] (C-292)
𝑙𝑜𝑔 ( )
𝑞𝑔
4 𝜋(0.00633)𝑘 1424𝑇 𝑙𝑛 𝑐𝑓𝑝𝑐
=( ) 𝑡 + log ( )
2.303 𝑙𝑛 𝑐𝑓𝑝𝑐 (𝜙𝜇𝑐𝑡 )𝑖 𝐴 𝑘ℎ
+ 𝑙𝑜𝑔 (2𝑠 𝑋)
4 𝜋(0.00633)𝑘 (C-293)
𝑐1 = 𝑚𝐹𝑃𝐶 = ( )
2.303 𝑙𝑛 𝑐𝑓𝑝𝑐 (𝜙𝜇𝑐𝑡 )𝑖 𝐴
243
Ilustración 0-18 Grafica de modelo semi- log del modelo tipo curva para régimen de
flujo dominado por la frontera externa con no intersección en cero de la línea recta,
bajo la condición de producción a presión de fondo fluyendo constante,𝒑𝒘𝒇 . (Arévalo
Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories, 2001)
Además, la Ec. C-292 produce una línea recta con no intersección en cero igual a:
1424𝑇 𝑙𝑛 𝑐𝑓𝑝𝑐 (C-295)
log 𝑐2 = log(𝑏𝐹𝑃𝐶 ) = log ( )
𝑘ℎ
El termino 𝑙𝑜𝑔 (2𝑠 𝑋) se desprecia por no considerar este tipo de daño. Sustituimos
la Ec. C-277 en la Ec. C-295:
2.2458 𝐴 (C-296)
1424𝑇 𝑙𝑛 ( 2 )
𝑟𝑤 𝐶𝐴
log(𝑏𝐹𝑃𝐶 ) = log ( )
𝑘ℎ
244
2.2458 𝐴 1 (C-300)
𝐶𝐴 = 2
𝑟𝑤 𝑏 2𝑘ℎ
exp ( 𝐹𝑃𝐶
1424𝑇 )
2.2458 𝐴 (C-301)
𝐶𝐴 =
𝑏𝐹𝑃𝐶 𝑘 ℎ
exp [2 ( 1424𝑇 )] 𝑟𝑤2
Ilustración 0-19 Grafica doble logarítmica donde se detecta un régimen de flujo lineal
temprano con valor de pendiente de un medio, seguido de un régimen de flujo Bilineal
con valor de la pendiente de la línea recta de un cuarto. (Arévalo Villagran, Analysis of
Long-term Behavior in tight gas resevoirs:case histories, 2001)
245
3.6.1 Modelos para régimen de flujo lineal temprano en un yacimiento homogéneo
a condiciones de producción de gasto constante de gas,𝒒𝒈 , y a presión de
fondo fluyendo constante,𝒑𝒘𝒇 .
Se aplica la relación anterior con las Ecs. C-302 y C-303, de la siguiente forma:
1 (C-305)
𝑃𝑤𝐷𝐿 =
𝑞𝑤𝐷𝐿
𝑃𝑤𝐷𝐿 (C-306)
=1
1⁄
𝑞𝑤𝐷𝐿
Reordenando términos:
2 1 (C-308)
𝑃𝑤𝐷𝐿 = ( ) ( )
𝜋 𝑞𝑤𝐷𝐿
1 𝜋 (C-
= 𝑃𝑤𝐷𝐿 ( )
𝑞𝑤𝐷𝐿 2 3188)
246
Conocemos el modelo de flujo lineal de respuesta para la presión adimensional
obtenido de modelos lineales para trayectorias con alta permeabilidad con flujo
lineal de fractura al comienzo del tiempo de explotación para el caso de líquidos con
la condición de explotación de producción a gasto constante de líquido,𝑞𝐿 :
2 (C-310)
𝑃𝑤𝐷𝐿 = √𝜋 𝜂𝐹𝐷 𝑡𝐷𝐿 + 𝑆 𝑔
(𝑘𝑓 𝑤)𝐷
1 𝜋 2 (C-311)
= ( )[ √𝜋 𝜂𝐹𝐷 𝑡𝐷𝐿 + 𝑆 𝑔 ]
𝑞𝑤𝐷𝐿 2 (𝑘𝑓 𝑤)
𝐷
1 𝜋 𝜋 𝑔 (C-312)
= √𝜋 𝜂𝐹𝐷 𝑡𝐷𝐿 + 𝑆
𝑞𝑤𝐷𝐿 (𝑘𝑓 𝑤)𝐷 2
2 (C-313)
𝑚𝑊𝐷𝐿 = √𝜋 𝜂𝐹𝐷 𝑡𝐷𝐿 + 𝑆 𝑔
(𝑘𝑓 𝑤)𝐷
1 𝜋 (C-314)
= √𝜋 𝜂𝐹𝐷 𝑡𝐷𝐿 + 𝑆 𝑔
𝑚𝑊𝐷𝐿 (𝑘𝑓 𝑤)𝐷
Las Ec. C-313 y C-314 son las ecuaciones que modelan el flujo lineal de fractura o
en una trayectoria de flujo con alta permeabilidad para la condición de explotación
a gasto de flujo constante,𝑞𝑔 y a presión de fondo fluyendo constante,𝑝𝑤𝑓 , para el
caso de gas, respectivamente, al comienzo del tiempo de explotación.
247
3.6.2 Modelo para régimen de flujo lineal temprano en un yacimiento homogéneo
infinito a condiciones de producción de gasto constante de gas,𝒒𝒈 .
248
(C-323)
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 1424 𝑇 2𝑘 𝑥𝑓 𝑘𝑓 0.00633 𝑡
=( )[ √𝜋 ( ) + 𝑆𝑔 ]
𝑞𝑔 𝑘ℎ 𝑘𝑓 𝑤 𝜙𝑓 𝑐𝑓 𝑡𝑖 𝜇𝑔𝑖 𝑥𝑓 2
Hacemos:
1424 𝑇 (C-324)
𝑋=( )
𝑘ℎ
Entonces:
Simplificando términos:
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 401.62 𝑇 (C-326)
= ( )√𝑡 + 𝑋𝑆 𝑔
𝑞𝑔 √𝜇𝑔𝑖 (𝜙 )
𝑐𝑡𝑖 𝑓 ℎ 𝑤 𝑘
√ 𝑓
Ilustración 0-20 Graficas de régimen de flujo lineal temprano con intersección en cero
(esbozo izquierdo) y no intersección en cero de la línea recta (esbozo derecho), bajo
la condición de producción a gasto de flujo constante de gas. (Arévalo Villagran,
Analysis of Long-term Behavior in tight gas resevoirs:case histories, 2001)
249
En una gráfica de 𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )⁄𝑞𝑔 𝑣𝑠 √𝑡, la pendiente se puede obtener
mediante la siguiente expresión:
401.62 𝑇 1 (C-328)
𝑐1 = 𝑚𝐿𝑇𝐺𝐶 = ( )
√𝜇𝑔𝑖 (𝜙𝑐𝑡𝑖 )𝑓 √𝑘𝑓 𝐴𝑐1
250
Sabemos que la pseudo- presión del gas real con la condición de explotación a
presión de fondo fluyendo constante es la siguiente:
Se hace la igualdad:
1 1 (C-335)
= + 𝑆𝑔
𝑞𝐷𝐿 𝑚𝑊𝐷
Sustituimos la Ec. C-334 y la Ec. C-333 en la Ec. C-335 para poder llegar al modelo
de la línea recta 𝑦 = 𝑚𝑥 + 𝑏 ⟹ ∆𝑝 = 𝑐1 √𝑡 + 𝑐2, por lo tanto, obtenemos:
251
(C-342)
[𝑚(𝑝𝑖 ) − 𝑚(𝑝𝑤𝑓 )] 1424 𝑇 𝜋𝑘 𝑥𝑓 𝑘𝑓 𝜙𝑖 𝐶𝑡𝑖 0.00633𝑘𝑡
= [ √𝜋 ( )( ) + 𝑆𝑔 ]
𝑞𝑔 𝑘ℎ 𝑘𝑓 𝑤 𝑘𝜙𝑓 𝑐𝑓 𝑡𝑖 (𝜙𝜇𝑐𝑡 )𝑖 𝑥𝑓 2
Hacemos:
1424 𝑇 (C-343)
𝑋=
𝑘ℎ
Por lo tanto:
252
Ilustración 0-21 Graficas de régimen de flujo lineal temprano con intersección en
cero(esbozo izquierdo) y no intersección en cero de la línea recta(esbozo derecho),
bajo la condición de producción de presión de fondo fluyendo constante,𝒑𝒘𝒇 . (Arévalo
Villagran, Analysis of Long-term Behavior in tight gas resevoirs:case histories, 2001)
253