Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

06-Modelo Radial

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 21

MODELO RADIAL

Universidad Nacional Autónoma de México.


Facultad de Ingeniería.
Curso intersemestral.
Introducción a la Simulación de Yacimientos.
Castañeda Nazario Héctor Iván.
Ayala García José Antonio.
¿QUÉ ES UN MODELO RADIAL?

Un modelo radial es una representación basada en la


forma de un cilindro, este nos ayuda a la interpretación del
flujo de fluidos en los núcleos y en las cercanías de un
pozo, este modelo solo cuenta con un pozo en el centro de
la malla.

2
1 DIMENSIÓN

1 Dimensión

3
2 DIMENSIONES

2 Dimensiones

4
3 DIMENSIONES

3 Dimensiones

5
FRONTERAS

Este modelo se discretiza de igual forma que un modelo


en coordenadas cartesianas.

Cuenta con 2 fronteras:


• Interna “Pozo”
• Externa “Límite del yacimiento”

6
ESPACIAMIENTO LOGARÍTMICO

El modelo se rige mediante un espaciamiento logarítmico “𝑎𝑙𝑔 ”,


y se interpreta de la siguiente manera:

1
𝑟𝑒 𝑛
𝑎𝑙𝑔 =
𝑟𝑤
re: Radio externo [ft]
rw: Radio del pozo [ft]
n: Número total de celdas
7
RADIO DE LA PRIMERA CELDA

La distancia de la primera celda se denomina “r1” y se


obtiene de la ecuación:

𝑎𝑙𝑔 (ln |𝑎𝑙𝑔 |)


𝑟1 = 𝑟𝑤
𝑎𝑙𝑔 − 1

r1: Radio de la primera celda [ft]


rw: Radio del pozo [ft]

8
RADIOS DE LA CELDA 2 A N

Para los radios de la celda 2 hasta la n se ocupa la ecuación:

𝑟𝑖 = 𝑎𝑙𝑔 𝑖−1 (𝑟1 )

i: Numero de celda
ri: Radio de la celda i [ft]
r1: Radio de la celda 1 [ft]

9
VB PARA LAS CELDAS 1 A N-1

El volumen bruto para las celdas que van desde la 1 hasta la n-1
se calcula de la forma:

2 2
𝑎𝑙𝑔 − 1 2 1
𝑉𝑏𝑖 = 𝑟 ∆𝜃 ∆𝑧 = [𝑓𝑡 3 ]
𝑎𝑙𝑔 2 ln |𝑎𝑙𝑔 2 | 𝑖 2

10
VB PARA LA CELDA N

El volumen bruto para las celda n se calcula de la forma:

2
ln |𝑎𝑙𝑔 | 2
𝑎𝑙𝑔 2 − 1
𝑎𝑙𝑔 − 1 1
𝑉𝑏𝑛 = 1 − 𝑟𝑒2 ∆𝜃 ∆𝑧 = [𝑓𝑡 3 ]
𝑎𝑙𝑔 2 ln |𝑎𝑙𝑔 2 | 2

11
TRANSMISIBILIDADES

Las transmisibilidades C2 y C5 se calculan de la forma:

𝛽𝑐 2𝜋𝑘𝑟𝑖 ∆𝑧𝑖
𝐶2𝑖 = 𝐶5𝑖 =
𝜇𝐵𝑜 ln |𝑎𝑙𝑔 |
Como en el caso de coordenadas cartesianas, es necesario el
calculo de C0.
𝐶0𝑖 = −(𝐶2𝑖 + 𝐶5𝑖 + 𝑔𝑡𝑖 )

12
LA ACUMULACIÓN

El termino de acumulación “gt” se obtiene de la misma forma


que en coordenadas cartesianas.

𝑉𝑏𝑖 ∅𝐶
𝑔𝑡𝑖 =
𝛼𝐵0 ∆𝑡

13
EJERCICIO

Propiedades Celda 1 Celda 2 Celda 3 Celda 4 Celda 5


Δz [ft] 30 30 30 30 30
μ [cp] 0.5 0.5 0.5 0.5 0.5
kr [mD] 150 150 150 150 150
∅ [%] 23 23 23 23 23

14
EJERCICIO

P=4000 psia
rw=0.25 ft
Tenemos condiciones de
Area=43560 acres
frontera.
B=1 rbl/stb
En la frontera izquierda
B°=1 rbl/stb
q=2000 BPD
Δt=1 días.
En la frontera derecha q=0
Tf=15 días.
BPD
Compresibilidad 1E-5
Δθ= 2π
15
RESULTADOS

re= 24576 [ft]

alg= 9.9659 [1]

16
RESULTADOS

Resultados Celda 1 Celda 2 Celda 3 Celda 4 Celda 5


r [ft] 0.6389 6.3672 63.4544 632.3775 6302.2
Vb [ft^3] 814.2843 8.0873E4 8.0322E6 7.9774E8 5.6118E10
C2 0 27.7189 27.7189 27.7189 27.7189
C5 27.7189 27.7189 27.7189 27.7189 0
gt 3.3354E-4 0.0331 3.2901 326.7691 2.2987E4
C0 -27.7192 -55.471 -58.728 -382.207 -2.3015E4
b 1998.7 -132.5083 -1.3160E4 -1.3071E6 -9.1948E7
17
MATRIZ SOLUCIÓN

18
RESULTADO

Resultados Celda 1 Celda 2 Celda 3 Celda 4 Celda 5


P a 1 día 3787.1 3859.2 3931.2 3995 4000
P a 15 días 3733.2 3805.3 3877.5 3949.4 3999.4

19
RESULTADOS
Δt=1 días.

20
FIN
21

También podría gustarte