Ramos Fernández Daniel Tarea 04 18/09/2019
Ecuación de Darcy[1]
En 1856 Darcy publicó su trabajo, en el cual se describı́an estudios experimentales de
flujo de agua a través de filtros de arena no consolidada, los cuales tenı́an como objetivo
procesar los requerimientos diarios de agua potable del pueblo de Dijon (Francia).
La Ecuación de Afluencia básica es la ley de Darcy en su forma diferencial, la ecuación
nos permite analizar el flujo de hidrocarburos desde el yacimiento hasta el pozo:
ka A dp
q=−
µ dx
A: área de flujo, pie2
ka : permeabilidad absoluta, mD
bl
q: gasto, D
µ: viscosidad, cP
lb
dp pg 2
dx : gradiente de presión, pie
Ecuación de Darcy para flujo lineal incluyendo efectos gravitacionales:
kA dp g dD
q=− − ρ
µ dx gc dx
Otros sistemas de unidades usados en la Ecuación de Darcy:
SI Inglés CGS Darcy Campo
K m2 f t2 cm2 Darcy mD
lbf dina
P Pa pg 2 cm2 atm psia
m3 pie3 cm3 cm3 ST B
Q s s s s D
lbf ·s
µ Pa · s pie2 cP cP cP
A m2 pie 2
cm2 cm2 pie2
Desarrollo de la Ecuación de Darcy en flujo lineal:
ka A dp
q=−
µ dx
qµ
dp = − dx
ka A
Z P2 Z L
qµ
dp = − dx
P1 ka A 0
qµ
P2 − P1 = − L
ka A
ka A(P2 − P1 )
q=−
µL
1
Ramos Fernández Daniel Tarea 04 18/09/2019
Ecuación de Darcy en flujo radial:
ka A dp
q=− ; A = 2πrh
µ dx
ka 2πrh dp
q=−
µ dr
qµ dr
dp = −
ka 2πh r
Z pe Z re
qµ 1
dp = − dr
pwf ka 2πh rw r
2πhka (pe − pwf )
q=−
µln rrwe
Para flujo de aceite:
0,00708hko (pws −pwf )
q= µo Bo ln( rrw
ko : permeabilidad del aceite, mD
e
)
h: espesor del yacimiento, pie
pws : presión del yacimiento, pglb2
lb
pwf : presión del fondo fluyendo, pg 2
re : radio de drene, pie
rw : radio del pozo, pie
Ecuación de continuidad en diferentes sistemas coordenados[2]
2
Ramos Fernández Daniel Tarea 04 18/09/2019
El desarrollo de la ecuación de continuidad (EC) parte de un balance de materia sobre un elemento de volumen
representativo (EVR) arbitrario con N fluidos. Para el caso de la figura anterior, donde se muestra una porción de
yacimiento en coordenadas cartesianas con flujo en las direcciones x, y y z; esta expresión puede obtenerse al realizar
un balance de materia en el medio poroso, para cada fluido f presente, como:
(mi − mo )f,x + (mi − mo )f,y + (mi − mo )f,z = mf,a (1)
donde mi , mo y ma representan la cantidad de materia que entra al sistema, la que sale y la que permanece acumulada.
A su vez, estos términos pueden expresarse para cualquiera de las n direcciones de flujo, como:
mf,n = (ṁA∆t)f,n = (ρvA∆t)f,n (2)
mac = mt+∆t − mt = [Vb φ(sρ)f ]t+∆t − [Vb φ(sρ)f ]t (3)
y dado que el fluido entra en las caras situadas en x, y y z y sale en las que están en x + ∆x, y + ∆y + ∆z, al considerar
que Vb es constante, la Ec. (1) se reescribe como:
h i h i
(ρuA)x − (ρuA)x+δx + (ρuA)y − (ρuA)y+δy + (ρuA)z − (ρuA)z+δz ∆t = Vb (φsρ)f,t+∆t − (φsρ)f,t (4)
f
que al dividir por ∆t y Vb resulta en:
1 h i 1 h i
(ρuA)x − (ρuA)x+δx + (ρuA)y − (ρuA)y+δy + (ρuA)z − (ρuA)z+δz ∆t = (φsρ)f,t+∆t − (φsρ)f,t (5)
Vb f ∆t
y como
Vb = ∆x∆y∆z, Ax = ∆y∆z, Ay = ∆x∆z, Az = ∆x∆y (6)
al aplicar lı́mites cuando ∆x∆y∆z → 0 y ∆t → 0 sobre la Ec. (5), se llega a la siguiente expresión:
∂ ∂ ∂ ∂
(ρux )f + (ρuy )f + (ρuz )f = − (φsρ)f (7)
∂x ∂y ∂z ∂t
La Ec. (7) puede ser representada a través de la divergencia (∇·) como:
∂
∇ · (ρu)f = − (φsρ)f (8)
∂t
Ecuación de continuidad en coordenadas cilı́ndricas
1 ∂ 1 ∂ ∂ ∂
(rρur )f + (ρuθ )f + (ρuz )f = − (φsρ)f (9)
r ∂r r ∂θ ∂z ∂t
3
Ramos Fernández Daniel Tarea 04 18/09/2019
Ecuación de continuidad en coordenadas esféricas
1 ∂ 2 1 ∂ 1 ∂ ∂
(r ρur )f + (ρuθ )f + (sin ΘρuΘ )f = − (φsρ)f (10)
r2 ∂r r sin Θ ∂θ r sin Θ ∂Θ ∂t
Ecuaciones de estado utilizadas en los problemas de flujo[3]
Ecuación de estado para fluido incompresible
ρ = constante
∂ρ
=0
∂t
Ecuación de estado para fluido ligeramente compresible
ρ = ρo (1 + CP )
Ecuaciones de estado, para fluidos compresibles
Para un gas ideal
m PM
ρ= =
V RT
Para un gas real
PM
ρ=
ZRT
Referencias
[1] Trejo Martı́nez J., C. (2019). Comportamineto de Yacimientos
[2] Gallardo Ferrrera, H. E., Becerra Zepeda, M. (2015). FUNDAMENTOS DE LA CARACTERIZACIÓN DINÁMI-
CA DE YACIMIENTOS. Ciudad Universitaria: Facultad de Ingenierı́a, UNAM.
[3] Hernández Garcı́a, M. A., Domı́nguez Vargas, G. C., (1984). APUNTES DE SIMULACIÓN MATEMÁTICA DE
YACIMIENTOS. Ciudad Universitaria: Facultad de Ingenierı́a, UNAM.
