Método de Mills

El método de Mills se desarrolló en la década de 1930 y es comúnmente utilizado en la

industria.
El método discutido aquí generalmente se atribuye a Mills, aunque algunas de las
características del método han evolucionado desde el desarrollo original.

Desarrollo del Método de Mills

Las ecuaciones para carga máxima de la varilla pulida (PPRL- peak polished rod load) y
carga mínima de varilla pulida (MPRL- minimum polished rod load) para unidades
(Clase I) con varillas de acero y convencionales son:

PPRL=W f +W r ( 1+∝ ) (1)

PPRL=máximacarga de la varilla pulida ,lbf
W f = peso del fluido en el área netadel émbolo ,lbf
W r = peso de las varillasen elaire , lbf
∝=factor de aceleración Mills

MPRL=W r + ( 1−∝−0.127 G ) (2)

MPRL=carga mínima de varilla pulida ,lbf
G=gravedad específica del fluido

Donde:
SN 2
∝= ( 3)
70500
S=longitud de carrera de la varilla pulida ,∈¿
N=velocidad de bombeo ,1 /min

Y
W f =0.433 GL ( A P −A r ) ( 4)
L=longitud de la varilla , ft
A p =área de la seccióntransversal del émbolo de la bomba ,¿ 2
Ar =área de la sección transversal de la varilla, ¿2

Estas ecuaciones suponen que solo las varillas experimentan aceleración y que la
aceleración que se produce es la máxima bajo movimiento armónico simple. La
fricción se desprecia en ambas cargas.
Se utiliza un conjunto similar de ecuaciones para las unidades de Clase III. Para
unidades de balance de aire:
PPRL=W f +W r ( 1+ 0.7 ∝) (5)
MPRL=W r + ( 1−1.3 ∝−0.127 G ) (6)

Para las unidades Mark II:

PPRL=W f +W r ( 1+ 0.6 ∝) (7)
 MPRL=W r + ( 1−1.4 ∝−0.127 G ) (8)

El efecto de contrapeso, la fuerza del contrapeso que coloca en la varilla pulida, puede
aproximarse mediante la siguiente ecuación:
PPRL+ MPRL
CBE= (9)
2
CBE=efecto de contrapeso , lbf

Como una aproximación, la siguiente ecuación puede usarse para obtener el torque
máximo para una unidad de balance convencional o balance de aire:

PT =0.283 S ( PPRL−MPRL) (10)

PT =torque máximo de la caja de cambios ,∈−lbf

Para unidades Mark II:

PT =0.25 S ( PPRL−1.3 MPRL ) (11)
PT =torque máximo d elacaja de cambios ,∈−lbf

La potencia de varilla pulida a menudo se aproxima por:

SN (PPRL+ MPRL)
PRHP= (12)
750.000
PRHP= potencia de la varilla pulida , hp

La potencia nominal de la placa de identificación del motor principal debe ser de 1,5 a 2
veces la potencia de la varilla pulida, el menor valor que se utiliza para motores de alto
deslizamiento y el valor más alto para motores de deslizamiento normales.

Una aproximación de la carrera efectiva del émbolo comúnmente utilizado con el

método de Mills es el siguiente:
S P=S+ e P−( e t + er ) (13)
S p=longitud de carrera de la bomba,∈¿
e p=alargamiento de varillas por aceleración ,∈¿
e r =alargamiento de varillas por gravedad ,∈¿
e t =alargamiento de latubería sin anclajes debido a la cargade fluido ,∈¿

Donde:
L2 ∝
e P =40.8 (14 )
Yr
Y r =módulo de young para el material de la varilla , psi

La constante (40.8) en la ecuación 14 se reemplaza comúnmente empíricamente con

46.5 al diseñar un sistema con varilla de acero cónica.
(W ¿ ¿ f +0.127 GW f ) L
e t =12 (15)¿
At Y t
At =área de la sección transversal de la pared de la tubería , ¿2
Y t =módulo de young para el material de latubería , psi

L1 Ln
e r =12W f
[ Y 1 Ar
+…+
]
Y n A rn
(16)

La ecuación correspondiente para el desplazamiento de la bomba es:

PD=0.1484 A P S P N (17)
PD=desplazamiento teórico de la bomba, B /D

Las ecuaciones 1 a 17 forman la base de un sistema de bombeo diseñado por el método

de Mills.

MÉTODO API RP-11L

En 1954, en un intento de desarrollar métodos más precisos, un grupo fabricantes
comisiona un estudio en el “Midwest Research Institute” para entender más acerca del
complejo comportamiento de los sistemas de bombeo mecánico. La API publicó el
resultado de este estudio en 1967 como “Recommended Practice 11L”. Desde su
aparición, la API RP-11L se ha convertido en el método de diseño más popular, sin
embargo, el método tiene muchas limitaciones debido a las suposiciones cuando fue
desarrollado.

Consideraciones:
 Llenado completo de la bomba de subsuelo (sin interferencia de gas)
 Varillas de acero con diseño API.
 Poco deslizamiento del motor.
 No debe existir grandes efectos de fricción o aceleración del fluido.
 No hay efectos por aceleración del fluido.
 Profundidades mayores a 2000 pies.

Las tablas Norma API RP 11L hacen más fácil la selección y evaluación de equipos de
bombeo, sin necesidad de cálculos excesivos.

DATOS REQUERIDOS
 Nivel fluido, H=…………………………………………ft
 Profundidad de la bomba, L=…………………...ft
 Velocidad de bombeo, N=………………………..SPM
 Longitud de carrera, S=………………………….....in
 Diámetro del émbolo de la bomba, D=………in
 No de varilla=………………………………………………..
 Peso específico del fluido, G=……………………..
 Tamaño de la tubería=……………………………….in
 Tamaño y diseño de la varilla de bombeo=……………

Datos requeridos para el diseño de bombas, valores hallados en “TABLA 1.

NORMA API RP 11L” y dependen del número de varilla y diámetro del émbolo de
la bomba.
W r = peso de las varillasen elaire , lbs/ pie(1)
F c =factor de frecuencia(2)
Er =constante elástica de varillas , pulgs /(lbs−ft)(3)
Et =constante elástica de latubería de producción , pulgs /lb− pie (4)

DESARROLLO DEL MÉTODO API RP-11L

Carga de fluido sobre la bomba

La carga de fluido sobre la bomba (Fo) depende de la gravedad específica del fluido, la
profundidad de levantamiento y el diámetro del émbolo de la bomba.

F 0 (lbs)=0.34 ×G × D 2 × H (5)
G=gravedad específica del fluido , adimensional
D=diámetro del émbolo de la bomba , pulg
H=nivel del fluido , pies

Estiramiento de varillas
El cálculo del estiramiento de varillas adimensional, (Fo/Skr), es una de las relaciones
claves para determinar una carta dinagráfica parecida. La constante elástica de la sarta
de varillas (Er) se encuentra en la “TABLA 1. NORMA API RP 11L” de acuerdo al
número de varilla y diámetro del émbolo de la bomba.
1
La relación indica que a la profundidad total de las varillas (profundidad de la
Kr
bomba) se estirarán una cantidad de pulgadas por cada libra aplicada sobre la varilla.
1
( pulgs /lb )=Er × L(6)
Kr
1
=constante elástica de la sartade varillas , pulg /lb
Kr
Er =constante elástica de varillas , pulgs /(lbs−ft)
L= profundidad d elabomba , pies

F0
La relación indica que a la profundidad total de las varillas (profundidad de la
Sk r
bomba) se estirarán un porcentaje de la carrera de superficie cuando se levanta la carga
de fluido.

F0 F 0 × Er × L
= (8)
Sk r S
S=longitud de carrera , pulg

Longitud efectiva de la carrera del émbolo

Calcular la velocidad de sonido de la varilla

 a o (ft /seg)=19.7 / √¿ ¿
W r = peso de las varillasen elaire , lbs/ pie

El peso de las varillas en el aire (Wr) en lbs/pie se encuentra en la “TABLA 1.

NORMA API RP 11L” de acuerdo al número de varilla y diámetro del émbolo de la
bomba.

N
Calcular la relación
No
N N× L
= (10)
N o 15 × ao
N=velocidad de bombeo , spm
N o =frecuencia natural de la sarta simple de varillas , spm

N
Calcular la relación
No
N
N No
= (11)
N o FC
N o =frecuencia natural de lasarta de varillas telescópicas , spm
F C =factor de frecuencia
1
Calcular la relación
Kt
1
( pulgs/lb)=Et × L(12)
Kt
1
=constante elástica de la porción de la sarta de tubería no anclada , pulgs/lb
Kt
Et =constante elástica de latubería de producción , pulgs /lb− pie
La constante Et es un dato que facilita el fabricante.
S
Calcular la relación, P
S
SP
(13)
S
Desde la “GRÁFICA 1. RELACIÓN ADIMENSIONAL (Sp/S) PARA CALCULAR
N
LA CARRERA EFECTIVA DEL ÉMBOLO” con el valor de trazamos hasta
No
F0
encontrar la línea correspondiente al valor y podemos determinar el valor de la
Sk r
S
relación P .
S

Calcular la carrera efectiva del émbolo (Sp)

 SP
S P ( pulg)= [( ) ] [(
S
× S − Fo×
1
Kt
(14 )
)]
S P=longitud efectiva de lacarrera del émbolo , pulgs
S=longitud de carrera de la varilla pulida , pulgs

Desplazamiento de la bomba
PD=0.1166 × S P × N × D2 (15)
PD=desplazamiento de labomba , bbls/día

Peso de las varillas en el aire

El peso de las varillas en el aire (Wr) en lbs/pie se encuentra en la “TABLA 1.
NORMA API RP 11L” de acuerdo al número de varilla y diámetro del émbolo de la
bomba.

Peso total de la sarta de varillas en el aire

El peso total de la sarta (W) será:
W (lbs)=W r × L(16)
W =peso total de la sarta de varillas en el aire , lbs
W r = peso de las varillasen elaire , lbs/ pie
L=longitud de la sarta , pies

Peso total de la sarta de varillas en flotación

W rf =W [ 1−0.128 ( G ) ] (17)
W rf = peso total de la sarta de varillas en flotación , lbs
G=gravedad específica del fluido , adimensional

Torque máximo en la caja de engranajes

La “GRÁFICA 4. RELACIÓN ADIMENSIONAL(2T/(S^2*kr)) PARA CALCULAR
EL TORQUE MÁXIMO EN LA CAJA DE ENGRANAJES” permite calcular una
2T
relación adimensional de torque
S2 kr ( )
, utilizando los factores adimensionales base

N F0
conocidos; y
N o Sk r

2T
(21)
S2k r
Wrf
Para determinar Ta debemos calcular la relación .
Skr
W rf
(18)
Sk r
Condición:
Wrf
=0.3 → Ta=1
Skr
Wrf
≠ 0.3 → T se determina por gráfica
Skr
La “GRÁFICA 5. VALOR DE AJUSTE (Ta) PARA CORREGIR EL TORQUE
Wrf
MÁXIMO” determina Ta con la relación ≠ 0.3 . Utilizando los factores
Skr
N F0
adimensionales base conocidos; y determinamos un %x que sirve para corregir
N o Sk r
Ta.
Ta (23)
Wrf
T a=1+(%x)×
Skr
10 [
−0.3
]
T a=factor de correcón de torque

2T
PT ( C ) =0.5 ×
( )
S2 k r
× Sk r × S × Ta, pulgs−lbs(26)

2T
PT ( M )=0.48 ×
( )S2k r
× Sk r × S ×Ta ,, pulgs−lbs

PT ( A )= [ 0.233× PPRL ( A )−0.3 × MPRL ( A ) ] × S , , pulgs−lbs

Carga máxima en la varilla pulida

La “GRÁFICA 2. RELACIÓN ADIMENSIONAL (F1/Skr) PARA CALCULAR
CARGA MÁXIMA EN LA VARILLA” permite obtener una relación adimensional
F1
( )
Sk r
, para calcular la carga máxima en la varilla pulida, utilizando los factores

N F0
adimensionales base conocidos; y .
N o Sk r
F1
(19)
Sk r
Calcular Sk r
S
Sk r= (7)
1
kr

F1
PPRL ( C ) =W rf +
[( ) ]
Sk r
× Sk r lbs(24)

F1
PPRL ( M )=W rf +0.75
[( )
Sk r ]
× Sk r +0.25 × F o lbs

F1
PPRL ( A )=W rf + 0.85
[( )
Sk r ]
× Sk r + 0.15× F o lbs
W rf = peso total de lasarta de varillas en flotación , lbs
F o=carga de fluido sobre labomba ,lbs
C= convencional
M= Mark II
A= balance de aire

Carga mínima en la varilla pulida

La “GRÁFICA 3. RELACIÓN ADIMENSIONAL (F2/Skr) PARA CALCULAR LA
CARGA MÍNIMA EN LA VARILLA” permite obtener una relación adimensional
F2
( )
Sk r
, para calcular la carga mínima en la varilla pulida, utilizando los factores

N F0
adimensionales base conocidos; y .
N o Sk r
F2
(20)
Sk r

F2
MPRL ( C )=W rf −
[( ) ]
Sk r
× Sk r (25)

F2 + F 1
MPRL ( M )=PPRL ( M )−
[( ) ]Sk r
× Sk r

F2+ F1
MPRL ( A )=PPRL ( A )−
[( ) ]
Sk r
× Sk r

El rango entre las cargas máximas y mínimas en la barra pulida, gobiernan los límites de
esfuerzos impuestos sobre la sarta de varillas, y son factores claves en la fatiga y su vida
útil.

Potencia en la varilla pulida

La “GRÁFICA 6. RELACIÓN ADIMENSIONAL(F3/Skr) PARA CALCULAR LA
F3
POTENCIA DEL MOTOR” permite obtener una relación adimensional ( )
Sk r
, para

calcular la potencia en la varilla pulida, utilizando los factores adimensionales base

N F0
conocidos; y .
N o Sk r
F3
(22)
Sk r
−6 F3
PRHP(hp)=2.53× 10 × × Sk r × S × N (27)
Sk r

Peso para el contrabalance de la unidad de bombeo

CBE ( C )=1.06 × ( W rf +0.5 × F o ) , lbs(28)
CBE ( M )=0.53× ( PPRL ( M ) + MPRL ( M ) ) , lbs
CBE ( A )=0.52 × ( PPRL ( A ) +1.25 × MPRL ( A ) ) , lbs

Procedimiento de cálculo para un formulario típico basado en API RP 11L

1. Recolección de datos, éstos pueden ser de una instalación existente o de datos

calculados.
2. Calculo de los parámetros adimensionales independientes
3. Utilizando las gráficas de diseño API, obtener los parámetros adimensionales
dependientes
4. A partir de los parámetros adimensionales dependientes determinar los parámetros
operacionales del sistema de bombeo mecánico.

FORMULARIO TÍPICO BASADO EN API RP 11L

CÁLCULOS DE DISEÑO PARA SISTEMAS DE BOMBEO MECÁNICO
DATOS REQUERIDOS
 Nivel fluido, H=…………………………………………ft
 Profundidad de la bomba, L=…………………...ft
 Velocidad de bombeo, N=………………………..SPM
 Longitud de carrera, S=………………………….....in
 Diámetro del émbolo de la bomba, D=………in
 No de varilla=………………………………………………..
 Peso específico del fluido, G=……………………..
 Tamaño de la tubería=……………………………….in
 Tamaño y diseño de la varilla de bombeo=……………

FACTORES DESDE “TABLA 1. NORMA API RP 11L”

W r =( 1)
F c =(2)
Er =( 3)
Et =( 3)

CALCULAR VARIABLES ADIMENSIONALES.

F 0=0.34 xGx D2 xH (5)
1
=Er xL=(6)
Kr
1
Sk r=S+ =(7)
Kr
F0
=(8)
Sk r
a o=19.7+ √¿ ¿
N
=NxL+ ( 15 x a0 ) =(10)
No
N N
= + F C =(11)
No No
RESOLVER PARA SP Y PD
SP
=(13)
S
SP
S P=[( ) ] [(
S
× S − Fo×
1
Kt
=(14)
)]
PD=0.1166 × S P × N × D 2=(15)

DETERMINAR PARÁMETROS ADIMENSIONALES

W =W r × L=(16)
W rf =W [ 1−0.128 ( G ) ] =( 17 )
W rf
=(18)
Sk r

REGISTRE LOS VALORES NO DIMENSIONALES DESDE LAS GRÁFICAS

API 11L
F1
=(19)
Sk r
F2
=(20)
Sk r
2T
=( 21)
S2k r
F3
=(22)
Sk r
Ta=( 23)

RESOLVER LAS SIGUIENTES OPERACIONES

C= convencional
M= Mark II
A= balance de aire
F1
PPRL ( C ) =W rf +
[( )
Sk r ]
× Sk r =(24)

F1
PPRL ( M )=W rf +0.75
[(
Sk r ) ]
× Sk r +0.25 × F o=¿

F1
PPRL ( A )=W rf + 0.85
[(Sk r ) ]
× Sk r + 0.15× F o =¿

F2
MPRL ( C )=W rf −
[( )
Sk r
× Sk r =(25)
]
F2 + F 1
MPRL ( M )=PPRL ( M )−
[( ) ]
Sk r
× Sk r =¿
 F2+ F1
MPRL ( A )=PPRL ( A )−
[( ) ]
Sk r
× Sk r =¿

2T
PT ( C ) =0.5 ×
( )
S2 k r
× Sk r × S × Ta=(26)

2T
PT ( M )=0.48 ×
( )
S2k r
× Sk r × S ×Ta=¿

PT ( A )= [ 0.233× PPRL ( A )−0.3 × MPRL ( A ) ] × S=¿

F
PRHP=2.53 × 10−6 × 3 × Sk r × S × N =(27)
Sk r
CBE ( C )=1.06 × ( W rf +0.5 × F o )=(28)
CBE ( M )=0.53× ( PPRL ( M ) + MPRL ( M ) ) =¿
CBE ( A )=0.52 × ( PPRL ( A ) +1.25 × MPRL ( A ) )=¿

ANEXOS
DATOS DE BOMBA Y VARILLA API
TABLA 1. NORMA API RP 11L
NORMA API RP 11L
NORMA API RP 11L
GRÁFICA 1. RELACIÓN ADIMENSIONAL (Sp/S) PARA CALCULAR LA CARRERA EFECTIVA DEL PISTÓN
GRÁFICA 2. RELACIÓN ADIMENSIONAL (F1/Skr) PARA CALCULAR CARGA MÁXIMA EN LA VARILLA

GRÁFICA 3. RELACIÓN ADIMENSIONAL (F2/Skr) PARA CALCULAR LA CARGA MÍNIMA

EN LA VARILLA
GRÁFICA 4. RELACIÓN ADIMENSIONAL(2T/(S^2*kr)) PARA CALCULAR EL TORQUE MÁXIMO EN LA
CAJA DE ENGRANAJES

GRÁFICA 5. VALOR DE AJUSTE (Ta) PARA CORREGIR EL TORQUE MÁXIMO

GRÁFICA 6. RELACIÓN ADIMENSIONAL(F3/Skr) PARA CALCULAR LA POTENCIA DEL MOTOR
Bibliografía
JENNINGS, J. W. (1989). THE DESIGN OF SUCKER ROD PUMP SYSTEMS. TEXAS.

PANCHI ZAPATA CRISTIAN RAUL, C. M. (2011). ESTUDIO PARA EL CAMBIO DE SISTEMA DE

LEVANTAMIENTO ARTIFICIAL A BOMBEO MECÁNICO EN EL CAMPO LAGO AGRIO.
QUITO.

PDVSA. (2001). DISEÑO DE INSTALACIONES DE LEVANTAMIENTO POR BOMBEO MECÁNICO.

VENEZUELA.