Estimulacion de Pozos
Estimulacion de Pozos
Estimulacion de Pozos
MATERIA:
RECUPERACION MEJORADA
AUTOR:
ACOSTA EDISON PAUL
CARPIO BRYAN STEVEN
GAVILANES JOCELYNE HAILYS
QUIZHPE FREDY JOSE
TEMA:
INYECCION DE SURFACTANTE NO IONICO-ALCOHOL
DOCENTE:
ING. FIDEL CHUCHUCA
CARRERA:
INGENIERÍA EN PETRÓLEO
AULA:
8/1
CONCLUSIONES ................................................................................................................ 25
RECOMENDACIONES ...................................................................................................... 25
BIBLIOGRAFIA .................................................................................................................. 26
INTRODUCCION
Esta teoría sobre la inyección de alcohol presentada se basa en datos físicos conocidos y hace
uso de ecuaciones de Buckley- Leverett. Para explicar el uso de la teoría tenemos tres
características necesarias para el desplazamiento de alcohol al momento de inyectarlo, la
primera es un sistema de dos fluidos inmiscibles, es decir, el caso Buckley- Leverett; la
segunda es un sistema de dos fluidos ligeramente miscibles que implican una extensión de
la ecuación de Buckley- Leverett y la tercera es un sistema de los
tres componentes, incluyendo la inyección de fluido puro rico en alcohol como un caso
especial, que implica una extensión más general de la ecuación de Buckley- Leverett.
El uso del alcohol isopropílico está limitado porque inicialmente es miscible con el petróleo
y el agua connata y, por lo tanto, el contenido de alcohol de la zona de mezcla se diluye por
debajo del nivel necesario para mantener la miscibilidad, por lo tanto, reduce su propiedad
de desplazamiento haciendo que se necesite inyectar altas cantidades de alcohol; el uso de
cantidades grandes de alcohol hace posible la efectividad del mecanismo inyector, pero es
más costoso.
OBJETIVO GENERAL
OBJETIVOS ESPECIFICOS
Desplazamiento miscible
Movilización miscible
Si la formulación de la fase acuosa es tal que ésta pueda formar una sola fase con el aceite,
entonces las fuerzas capilares no entran en línea de cuenta ya que la interface desaparece. El
desplazamiento miscible puede obtenerse con soluciones alcohólicas, con soluciones
concentradas de surfactante o con dióxido de carbono supercrítico. Al obtenerse las
condiciones del desplazamiento miscible, la ley de Darcy se aplica; sin embargo, el fluido no
es homogéneo desde el punto de vista de la densidad y de la viscosidad y pueden producirse
segregaciones. El estudio de la distribución de tiempos de residencia permite caracterizar las
propiedades del flujo.
Una baja en la tensión interfacial ayuda a mejorar la producción del hidrocarburo confinado,
sabiendo que el número de capilaridad del fluido debe estar elevado para que se cumpla.
Fig. 1. Crudo recuperado vs número de capilaridad.
La teoría del desplazamiento inmiscible unidimensional fue desarrollada originalmente por
Buckley y Leverett, además nos habla del principio de conservación masiva, otorgando tres
clases:
La base de esta teoría es que la ley de Darcy es válida para cada uno de los dos fluidos que
fluyen simultáneamente a través de un medio poroso. El sujeto tiene varios supuestos, el más
importante es que la presión capilar es cero, la ecuación de Darcy se reduce a las de Buckley
Leverett
𝜕𝑆𝐷 𝑈 𝜕𝑓𝐷
( ) =− ( )
𝜕𝑡 𝑥 ∅𝐴 𝜕𝑥 𝑡
Esta ecuación diferencial se resolvió para obtener la ecuación de Buckley- Leverett, a partir
de la cual se deduce que cualquier saturación se produce a través del medio poroso con la
velocidad dada por la pendiente de la curva de flujo vs saturación.
Estas diferencias de velocidad en un perfil de saturación pueden generar una discontinuidad
de saturación, y la velocidad de esta discontinuidad se puede determinar mediante el uso de
la construcción de la ecuación de Buckley- Leverett. Por lo tanto, para concentraciones de
agua constantes, el perfil de saturación se puede determinar en diferentes momentos del
proceso de desplazamiento por inyección de alcohol.
Las características a tomar en cuenta para esta teoría de desplazamiento son tres:
La siguiente figura muestra un ejemplo de curva de flujo de la fase acuosa 𝑢𝑤 trazada contra
la saturación de la fase acuosa 𝑆𝑤 . El flujo se puede calcular fácilmente a partir de
viscosidades conocidas y permeabilidades relativas, la pendiente de la curva es cero para las
saturaciones entre 0 y 𝑆1, y entre 𝑆2 y 1, lo que indica la velocidad cero de estas saturaciones.
Entre 𝑆1 y 𝑆2 , la velocidad de saturación es finita y alcanza un máximo en la pendiente más
pronunciada de la curva 𝑢𝑤 : 𝑆𝑤 . El medio poroso en el que tiene lugar el proceso de
desplazamiento unidimensional se muestra a continuación, la saturación de la fase acuosa 𝑆𝑤
se mide verticalmente y la distancia x horizontalmente. Ahora se puede considerar un
elemento de medio poroso de longitud ∆𝑥 que contiene la interfaz A: B. Supongamos que la
saturación de la fase acuosa 𝑆𝑤 en este elemento es inicialmente como se muestra por la línea
de puntos.
El fluido Y es desplazado por el fluido X del medio poroso unidimensional como se muestra
en la siguiente figura. Estos fluidos no están en equilibrio entre sí y, por lo tanto, cuando al
comienzo del proceso de desplazamiento tienen una interfaz común, tienen lugar hasta que
se crea una zona de amortiguamiento que contiene los fluidos A y B.
Fig. 4. Sistema de dos fluidos ligeramente miscibles-comportamiento de interface.
Sistema de los tres componentes. –
𝜕𝑆𝐷 𝑈 𝜕𝑓𝐷
( ) =− ( )
𝜕𝑡 𝑥 ∅𝐴 𝜕𝑥 𝑡
𝐴 = 𝐶𝑤𝑤 𝑢𝑤 + 𝐶𝑤𝑜 𝑢𝑜
𝑎, 𝑏 = 𝑟𝑎𝑑𝑖𝑜𝑠 𝑑𝑒𝑓𝑖𝑛𝑖𝑑𝑜𝑠
𝐶 = 𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛
𝑘 = 𝑝𝑒𝑟𝑚𝑒𝑎𝑏𝑖𝑙𝑖𝑑𝑎𝑑 𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑎
𝑚 = 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒 𝑑𝑒𝑓𝑖𝑛𝑖𝑑𝑎
𝑡 = 𝑇𝑖𝑒𝑚𝑝𝑜
𝑢 = 𝐹𝑙𝑢𝑗𝑜 𝑣𝑜𝑙𝑢𝑚é𝑡𝑟𝑖𝑐𝑜
𝛼 = 𝐶𝑜𝑤 𝑢𝑤 + 𝐶𝑜𝑜 𝑢𝑜
𝛽 = 𝐶𝑜𝑤 𝑆 + 𝐶𝑜𝑜 (1 − 𝑆)
𝜇 = 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑑𝑎𝑑 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎
𝜌 = 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑
∅ = 𝑝𝑜𝑟𝑜𝑠𝑖𝑑𝑎𝑑
𝑤 = 𝑓𝑎𝑠𝑒 𝑎𝑐𝑢𝑜𝑠𝑎
𝑜 = 𝑓𝑎𝑠𝑒 𝑎𝑐𝑒𝑖𝑡𝑜𝑠𝑎
La conservación de las ecuaciones de masa puede escribirse para el aceite y para el agua
como:
𝑋𝑤 − 𝑋𝑤 𝑣1 = 𝑢𝑤 𝑤 − 𝑣1 𝑆𝑤 𝑤 + 𝑢𝑜 𝐶𝑤𝑜 − 𝑣1 𝑆𝑜 𝐶𝑤𝑜
𝑢𝑤 − 𝑣1 𝑆𝑤 𝑢𝑜 − 𝑣1 𝑆𝑜
𝑋𝑜 = ( ) 𝐶𝑜𝑤 + ( ) 𝐶𝑜𝑜
1 − 𝑣1 1 − 𝑣1
𝑢𝑤 − 𝑣1 𝑆𝑤
𝑚=
1 − 𝑣1
𝑋𝑜 = 𝑚 𝐶𝑜𝑤 + (1 − 𝑚) 𝐶𝑜𝑜
𝑢𝑤 − 𝑣1 𝑆𝑤 𝑢𝑜 − 𝑣1 𝑆𝑜
𝑋𝑤 = ( ) 𝑤+( ) 𝐶𝑤𝑜
1 − 𝑣1 1 − 𝑣1
𝑋𝑤 = 𝑚 𝑤 + (1 − 𝑚) 𝐶𝑤𝑜
𝐵𝑋 𝐵𝑌
𝑎 = 𝐴𝐵 𝑏 = 𝐴𝐵
Entonces
𝑢𝑤 − 𝑣1 𝑆𝑤
𝑎=𝑚=
1 − 𝑣1
La velocidad 𝑣1 de la discontinuidad de saturación es, por lo tanto:
𝑢𝑤1 − 𝑎
𝑣1 =
𝑆𝑤1 − 𝑎
La velocidad 𝑣2 de la segunda discontinuidad de saturación se obtiene de manera similar:
𝑢𝑤2 + 𝑏
𝑣2 =
𝑆𝑤2 + 𝑏
𝜕(𝜌 𝑣) 𝜕𝜌
+ =0
𝜕𝑥 𝜕𝑡
La conservación de las ecuaciones de masa se puede escribir para cada uno de los tres líquidos
constituyentes del sistema ternario considerado aquí, pero, como se asumió que los líquidos
eran incompresibles, la tercera ecuación es una continuación de las otras dos. La saturación
de la fase acuosa se puede definir como:
𝜕 𝜕
(𝐶𝑤𝑤 𝑢𝑤 + 𝐶𝑤𝑜 𝑢𝑜 ) + ∅𝐴 (𝐶𝑤𝑤 𝑆 + 𝐶𝑤𝑜 (1 − 𝑆)) = 0
𝜕𝑥 𝜕𝑡
Y de petróleo
𝜕 𝜕
(𝐶𝑜𝑤 𝑢𝑤 + 𝐶𝑜𝑜 𝑢𝑜 ) + ∅𝐴 (𝐶𝑜𝑤 𝑆 + 𝐶𝑜𝑜 (1 − 𝑆)) = 0
𝜕𝑥 𝜕𝑡
𝜕𝐴 𝜕𝐵
+ =0
𝜕𝑥 𝜕𝑡
𝜕𝛼 𝜕𝛽
+ =0
𝜕𝑥 𝜕𝑡
Los flujos de la fase acuosa 𝑢𝑤 y de la fase del petróleo 𝑢𝑜 se determinaron a partir de la ley
de Darcy; la ecuación está sujeta a las suposiciones habituales sobre la gravedad, la tensión
superficial, las fuerzas potenciales, etc.
Las viscosidades de las dos fases se escribieron como funciones lineales de la concentración
de agua en la fase acuosa w, a saber:
𝑆 − 𝑆1 4
𝑘𝑤 = ( ) 𝑆1 = 0,2 ≤ 𝑆 ≤ 1
1 − 𝑆1
y para la fase oleica
𝑆2 − 𝑆1
𝑘𝑜 = 𝑆1 = 0,2 ≤ 𝑆 ≤ 𝑆2 = 0,6
𝑆2 − 𝑆1
𝑘𝑜 = 0 𝑆2 = 0,6 < 𝑆 ≤ 1
Los flujos de las dos fases por lo tanto fueron funciones de la saturación de agua S y la
concentración de agua w. Para la fase acuosa,
𝜇𝑤
𝜇𝑜
𝑢𝑤 = 𝑢𝑤 (𝑆, 𝑤) = 𝑈
𝜇𝑤 𝑘𝑤
𝜇𝑜 + 𝑘𝑜
𝑘𝑤
𝑘𝑜
𝑢𝑜 = 𝑢𝑜 (𝑆, 𝑤) = 𝑈
𝜇𝑤 𝑘𝑤
𝜇𝑜 + 𝑘𝑜
𝑢𝑤 (𝑆, 𝑤) + 𝑢𝑜 (𝑆, 𝑤) = 𝑈
Las concentraciones de los componentes en las fases acuosa y oleica se expresaron como
funciones de la concentración de agua en la fase acuosa w = 𝐶𝑤𝑤 . Las constantes numéricas
se determinaron por computadora utilizando un programa de mínimos cuadrados modificado
obteniendo lo siguiente:
𝑢𝑤
𝑆𝑤 =
𝑈
Esta solución fue obtenida por Gardner y Warren y no se reproducirá aquí. Si este 𝑆𝑤 se
sustituye en las funciones A, B, 𝛼 y 𝛽 definidas en el apéndice B, se sigue que:
𝐴=𝑈𝐵 𝛼=𝑈𝛽
Por lo tanto
𝜕𝐵 ∅𝐴 𝜕𝐵
+ =0
𝜕𝑥 𝑈 𝜕𝑡
𝜕𝛽 ∅𝐴 𝜕𝛽
+ =0
𝜕𝑥 𝑈 𝜕𝑡
𝑈
𝐵(𝑆, 𝑤) = 𝑓1 (𝑥 − 𝑡)
∅𝐴
𝑈
𝛽(𝑆, 𝑤) = 𝑓2 (𝑥 − 𝑡)
∅𝐴
Donde𝑓1 y 𝑓2 denotan funciones arbitrarias.
La solución iso-tachykinetic implica que para un 𝑆𝑤 particular en el medio poroso hay una
w correspondiente. Para una w fija, se puede calcular y trazar una curva de flujo en el plano
𝑢𝑤 : 𝑆𝑤 . En general esta curva tiene una forma de S. El 𝑆𝑤 para el que se mantiene la solución
iso-tachykinetic es el identificado por el punto de intersección de la curva de flujo y la
diagonal a través del origen.
Inicialmente, t es cero, para cualquier otra t debe haber una X correspondiente, de manera
que para 𝛽(𝑆, 𝑤) fija, el argumento (𝑥 − 𝑈𝑡/∅𝐴) permanece sin cambios. Esto significa que
los planos de constante 𝑆𝑤 y constante correspondiente w viajan hacia adelante a través del
medio poroso con una velocidad 𝑈/∅𝐴.
Función I(w).
Se usaron valores particulares de I (w) en la construcción de Buckley-Leverett, donde se
señalaron sus propiedades geométricas. La función fue utilizada también en cálculos de
velocidades de planos de concentraciones constantes, obteniendo las siguientes variables:
∆𝐶𝑜𝑤 = 𝐶´𝑜𝑤 ∆𝑤
∆𝐶𝑤𝑜 = 𝐶´𝑤𝑜 ∆𝑤
∆𝐶𝑜𝑜 = 𝐶´𝑜𝑜 ∆𝑤
Donde los símbolos imprimados denotan derivados con respecto a w. Las ecuaciones de línea
recta para las dos líneas por lo tanto son:
𝑦 − 𝐶𝑜𝑜 𝑥 − 𝐶𝑤𝑜
=
𝐶𝑜𝑜 − 𝐶𝑜𝑤 𝐶𝑤𝑜 − 𝑤
𝑦 − 𝐶𝑜𝑜
𝐼(𝑤) =
𝐶𝑜𝑜 − 𝐶𝑜𝑤
Los símbolos imprimidos denotan derivación con respecto a w. De ahí nos lleva a la
interpretación geométrica de I(w)
𝑦 − 𝐶𝑜𝑜 𝐵𝐶
𝐼(𝑤) = =
𝐶𝑜𝑜 − 𝐶𝑜𝑤 𝐴𝐵
𝜕𝑤 𝜕(𝑆, 𝑤) 𝜕𝑥
( ) = ( )
𝜕𝑡 𝑥 𝜕(𝑥, 𝑡) 𝜕𝑆 𝑤
𝜕𝑤 𝜕(𝑆, 𝑤) 𝜕𝑡
( ) =− ( )
𝜕𝑡 𝑡 𝜕(𝑥, 𝑡) 𝜕𝑆 𝑤
𝜕𝑆 𝜕(𝑆, 𝑤) 𝜕𝑥
( ) =− ( )
𝜕𝑡 𝑥 𝜕(𝑥, 𝑡) 𝜕𝑤 𝑠
𝜕𝑆 𝜕(𝑆, 𝑤) 𝜕𝑡
( ) =− ( )
𝜕𝑡 𝑡 𝜕(𝑥, 𝑡) 𝜕𝑤 𝑠
𝜕𝐴 𝜕𝑡 𝜕𝐴 𝜕𝑡 𝜕𝐵 𝜕𝑥 𝜕𝐵 𝜕𝑥
− − + =0
𝜕𝑆 𝜕𝑤 𝜕𝑤 𝜕𝑆 𝜕𝑆 𝜕𝑤 𝜕𝑤 𝜕𝑆
𝜕𝛼 𝜕𝑡 𝜕𝛼 𝜕𝑡 𝜕𝛽 𝜕𝑥 𝜕𝛽 𝜕𝑥
− − + =0
𝜕𝑆 𝜕𝑤 𝜕𝑤 𝜕𝑆 𝜕𝑆 𝜕𝑤 𝜕𝑤 𝜕𝑆
Quedando:
𝜕𝐴 𝜕𝛽 𝜕𝐵 𝜕𝛼
− =0
𝜕𝑆 𝜕𝑆 𝜕𝑆 𝜕𝑆
𝜕𝑡 1
= … … … … … … … … … … … … . (1)
𝜕𝑆 𝐹(𝑆, 𝑤)
𝜕𝑡 1 𝜕𝑥 1 𝜕𝑥
= + … … … . . (2)
𝜕𝑤 𝐺(𝑆, 𝑤) 𝜕𝑤 𝐻(𝑆, 𝑤) 𝜕𝑆
Donde:
𝜕𝐴 𝜕𝛼 𝜕𝛼 𝜕𝐴
−
𝐹(𝑆, 𝑤) = 𝜕𝑤 𝜕𝑆 𝜕𝑤 𝜕𝑆
𝜕𝐵 𝜕𝛼 𝜕𝛽 𝜕𝐴
−
𝜕𝑤 𝜕𝑆 𝜕𝑤 𝜕𝑆
𝜕𝐴
𝐺(𝑆, 𝑤) = 𝜕𝑆
𝜕𝐵
𝜕𝑆
𝜕𝐴 𝜕𝛼 𝜕𝛼 𝜕𝐴
−
𝐻(𝑆, 𝑤) = 𝜕𝑤 𝜕𝑆 𝜕𝑤 𝜕𝑆
𝜕𝐵 𝜕𝛼 𝜕𝐴 𝜕𝛽
−
𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤
Sustituyendo 1 y 2:
1 𝜕𝑥 1 𝜕𝑥 1 𝜕𝑥
𝑑𝑡 = 𝑑𝑆 + 𝑑𝑤 + 𝑑𝑤 … … … … … … (5)
𝐹(𝑆, 𝑤) 𝜕𝑆 𝐺(𝑆, 𝑤) 𝜕𝑤 𝐻(𝑆, 𝑤) 𝜕𝑆
Las funciones F (S, w) y H (S, w) pueden representarse en términos de concentración w, flujo
de Cwo y C de la fase acuosa y saturación S para dar
𝑢𝑤 + 𝐼(𝑤)
𝐹(𝑆, 𝑤) =
𝑆 + 𝐼(𝑤)
𝜕𝑢 𝑤
{𝑢𝑤 + 𝐼(𝑤)}
𝐻(𝑆, 𝑤) = 𝜕𝑆
𝜕𝑢 𝑤
𝐽(𝑤)(𝑈𝑤 − 𝑆) + {𝑆 + 𝐼(𝑤)}
𝜕𝑤
𝑢𝑤 + 𝐼(𝑤) 𝜕𝑢 𝑤
𝐼(𝑤) 𝜕𝑆
𝐻(𝑆, 𝑤) = … … … … … … … … … (9)
𝐽(𝑤) 𝑆 + 𝐼(𝑤) 𝜕𝑢 𝑤
(𝑢𝑤 − 𝑆) +
𝐼(𝑤) 𝐼(𝑤) 𝜕𝑤
Las funciones F(S, w), G (S, w) y H (S, w) pueden ser examinadas más de cerca con la ayuda
de un flujo arbitrario 𝑢𝑤 : curva de saturación S.
𝑢𝑤 + 𝐼(𝑤)
𝐹(𝑆, 𝑤) = … … … … … … … … … … … … … . (10)
𝑆 + 𝐼(𝑤)
𝜕𝑢 𝑤
𝐺(𝑆, 𝑤) = … … … … … … … … … … … … . . (11)
𝜕𝑆
𝜕𝑢 𝑤
{𝑢𝑤 + 𝐼(𝑤)}
𝐻(𝑆, 𝑤) = 𝜕𝑆
𝐽(𝑤)(𝑢𝑤 − 𝑆)
𝑢𝑤 + 𝐼(𝑤) 𝜕𝑢 𝑤
𝐼(𝑤) 𝜕𝑆
𝐻(𝑆, 𝑤) = … … … … … … … … … … … … … . . (12)
𝐽(𝑤)
(𝑢 − 𝑆)
𝐼(𝑤) 𝑤
Fig. 8. Curva de flujo vs saturación de la fase acuosa.
RECOMENDACIONES
REED, J. J.-I.-R. (1961). Mechanism of Alcohol Displacement of Oil from Porous. SPE 1536-
G.
Wachmann, C. (s.f.). A Mathematical Theory for the Displacement of Oil. SPE 879.