WWW - Ingenieria.unam - MX/ Revistafi: ISSN 1405-7743
WWW - Ingenieria.unam - MX/ Revistafi: ISSN 1405-7743
PP. 1-126
1-16 Seguridad estructural de plataformas autoelevables petroleras en operación Universidad Nacional Autónoma de México w Facultad de I ngeniería w México 2017
Campos Dante, Cabrera-Miranda José Manuel, Martínez-Mayorga Joan Manuel,
•
López-Cornejo Marcelino, Espinosa-Zamudio Esteban Ernesto
eneo-marzo 2017
17-29 Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo,
Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
31-42 Avances en investigación sobre el compostaje de biorresiduos en municipios menores de países en desarrollo. Lecciones
desde Colombia
Oviedo-Ocaña Edgar Ricardo, Marmolejo-Rebellon Luis Fernando, Torres-Lozada Patricia
•
43-53 Nuevo modelo para la expansión de sistemas eléctricos de distribución con generación distribuida considerando un
planeamiento multi-etapa coordinado
Núm. 1
Hincapié-Isaza Ricardo Alberto, Home-Ortiz Juan Manuel, Gallego-Rendón Ramón Alfonso
55-65 Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales and Marketing for Soft Drinks
Bottling Industry
Vol. XV I II
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen, De la Rosa-Hernández Romina Marisol
•
75-86 Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de Jesús, Pérez-Ruiz Juan Luis
111-126 Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species
as Bioindicators
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana,
Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
…
Ing. Leopoldo Silva Gutiérrez Dr. Alberto Ken Oyama Nakagawa
Secretario Administrativo Secretario de Desarrollo Institucional
Comité Editorial
Álvaro Aldama Rodríguez Dante Jaime Morán Zenteno
Consultor Instituto de Geología, UNAM
…
Revista Ingeniería Investigación y Tecnología
Gerardo Ruíz Solorio
Director Técnico y Editor
Eva Myriam Soroa Zaragoza Cuitláhuac Ernesto Sánchez Basilio Yuritzi Renee Olvera Rangel
Coordinación Editorial Coordinación de Arbitraje Coordinación de Diseño y Edición
Comité Científico Asesor
55-65 Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales
and Marketing for Soft Drinks Bottling Industry
Uso de mapas de dirección de valor para la selección de proyectos de Six Sigma: caso
de estudio en ventas y mercadotécnica para la industria embotelladora de bebidas
carbonatadas
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen,
De la Rosa-Hernández Romina Marisol
75-86 Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
Prueba de algoritmo de diagnóstico de turbinas de gas usando el Software ProDiMES
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de Jesús, Pérez-Ruiz Juan Luis
111-126 Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato
City: Use of Saxicolous Lichen Species as Bioindicators
Estudio de la distribución de metales pesados en la atmosfera de la ciudad de Guanajuato:
uso de especies de líquenes saxícolas como bioindicadores
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana,
Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
Ingeniería. Investigación y Tecnología es una publicación trimestral de la Facultad de Ingeniería de la
Universidad Nacional Autónoma de México.
Ingeniería. Investigación y Tecnología is a journal published quarterly by the Facultad de Ingeniería, UNAM.
Certificado de Reserva al Uso Exclusivo del Título Núm. 04-1998-060218085600-102, expedido por el
Instituto Nacional del Derecho de Autor el 2 de junio de 1998.
Certificado de Licitud de Título Núm. 10273 y de Contenido Núm. 7202, expedidos por la Comisión
Calificadora de Publicaciones y Revistas Ilustradas el 7 de enero de 1998.
Distribución: Dirección General de Publicaciones y Fomento Editorial, UNAM. Avenida del Imán
Núm. 5, Ciudad Universitaria, Coyoacán 04510, México, D.F. Tel. y fax (0155) 56 22 65 82.
Oficinas: Revista Ingeniería. Investigación y Tecnología: Edificio Principal de la Facultad de Ingeniería,
UNAM. Biblioteca Antonio Dovalí J. Circuito Interior, Ciudad Universitaria, CP. 04510, Delegación
Coyoacán, México DF o Apartado Postal 70-440, Administración de Correos ”Ciudad Universitaria”.
Tel. 5622 0876, Fax. 5616 1899. Lada Nacional (0155), Lada Internacional (5255).
Sitio Web: [Link]/~revistafi
E-mail: [Link]@[Link]
ISSN 1405-7743
Campos Dante
Instituto Mexicano del Petróleo (IMP)
Dirección de Servicios de Ingeniería
Correo: dcampos@[Link]
López-Cornejo Marcelino
Instituto Mexicano del Petróleo
Dirección de Servicios de Ingeniería
Correo: mrclopez@[Link]
Resumen
Se presentan un modelo de riesgo y otro de confiabilidad que definen la seguridad a exigir y la disponible, respecti-
vamente, de una plataforma marina autoelevable de piernas independientes ante el peligro meteorológico y oceano-
gráfico (metoceánico). Las exigencias se definen en términos del índice de confiabilidad óptimo anual, del factor de
reserva de resistencia (RSR) y de solicitaciones metoceánicas. La seguridad disponible se define en términos del ín-
dice de confiabilidad calculado a partir de tres estados límite de falla: el vuelco de la plataforma, el hundimiento
vertical y deslizamiento lateral de alguna de sus zapatas (spud-can). Si la seguridad disponible es mayor que la que
se exige, entonces, el empleo de la autoelevable será adecuado; de lo contrario, la autoelevable y su entorno estarían
en una situación muy vulnerable. Además estos modelos serán una herramienta para decidir rentar o no una au-
toelevable. Al final se realiza una aplicación numérica para dos lugares en el Golfo de México que mostraron un buen
comportamiento del sistema empleado.
Descriptores: seguridad requerida, seguridad disponible, elevación del casco, volteo de la autoelevable, costos de las
autoelevables.
Seguridad estructural de plataformas autoelevables petroleras en operación
Abstract
Algorithms and calculations are derived to define the target safety level and real safety level for an independent legs offshore
jack-up platform before meteorological and oceanographic hazards in a specific location. The risk model helps to define the re-
quired safety level which is defined by the annual optimal reliability index, the reserve strength ratio (RSR) and the maximum
environmental loads which must survive. Furthermore, the reliability analysis allows to know the installed (real) safety level of
the jack-up in terms of reliability index, for which three failure limit states are considered as follows: overturning of the platform,
vertical penetration of spud-cans and horizontal sliding of the foundation. If real safety level is higher than the target safety
level, then the performance of the jack-up platform is adequate; on the contrary, the jack-up platform would be under a very
vulnerable situation. This work is a tool for decision making on the rent of a jack-up rig. A numerical application for two loca-
tions in the Gulf of Mexico, which considers different water depths, shows a good safety behavior of this kind of offshore units for
the studied conditions.
Keywords: safety required, safety available, elevation the hull, Jack-up overturning, Jack-up costs.
a) b)
Figura 1. Plataforma autoelevable, a) La autoelevable GSF Adriatic VII prestando servicio a una plataforma
fija (Sharples, 2008), b) Isométrico del modelo empleado de la autoelevable
2 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 3
Seguridad estructural de plataformas autoelevables petroleras en operación
y su correspondiente índice de confiabilidad óptimo es Pérdidas por reposición, por retiro y abandono,
del pozo y porcentaje de la plataforma fija
bot = − Ф−1 [Pfot] (4)
La pérdida total de la autoelevable se considera me-
donde Ф(.) es la función de distribución normal están- diante el costo de reposición (CR). Se manejan valores
dar. constantes para el retiro y abandono de escombros (CA),
La figura 2, que representa a la Ec. 2, muestra que la pérdida del pozo en ejecución (CPZ) y la fracción del
para probabilidades de falla bajas, el costo inicial con- costo por perder la plataforma fija a la que se presta
trola el costo total; mientras que para probabilidades de servicio (γ CP), donde γ es un factor de reducción entre
falla altas, domina el costo futuro. El punto de costo to- cero y uno. Como hipótesis de análisis se considera
tal esperado mínimo define la Pf ot. que: a) la falla y el remplazo de la estructura ocurren
Una fracción (η) de la Pf ot es la probabilidad de falla una sola vez en el tiempo L, y b) la falla es igualmente
ambiental óptima anual (Pfo), asociada a los eventos me- probable en cualquier tiempo. Así, el valor esperado de
toceánicos extremos, así Pfo = η Pf ot y el índice de confia- los costos mencionados es
bilidad ambiental óptimo a exigir es
E[CR] + E[CA] + E[CPZ] + E[CP] =
βo = ‒Φ‒1[Pf o] (5) Pf (CR + CA + CPZ + γCP) PVF (r,L) (8)
∫ Pf CU (τ ) e − rrτdτ ,
L
E[CF] = E[CH] + E[CL] + E[CR] + E[CA] + E[CU] + E[CP] (6) E [CU] = 0
4 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
E [CF] = ζPf N (C1H + C1L )PP PVF(r , L) + Pf (C R + C A + C PZ + vación mínima de la cara inferior del casco de la au-
toelevable (ECMIN). Se estima que ECMIN = k HU, con
γ C P )PVF(r , L) + l Pf U d PVF(r , L) k = 0.726 aproximadamente, que considera la elevación
de la cresta y las mareas astronómica y de tormenta, tal
Entonces que, la ola última no impacte el casco.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 5
Seguridad estructural de plataformas autoelevables petroleras en operación
Tabla 2. Altura de ola última y elevación del casco asociados a los índices de confiabilidad mínimos
βo MIN = 3.15 βo MIN = 2.94 βo MIN = 2.81
Profundidad H100
Localización RSR = 1.18 RSR = 1.01 RSR = 0.91
(m) (m)
HU (m) ECMIN (m) HU (m) ECMIN (m) HU (m) ECMIN (m)
L-1 50.8 13.51 14.69 10.66 13.56 9.84 12.91 9.37
L-2 90.0 15.00 16.31 11.84 15.05 10.93 14.33 10.40
Tabla 3. Altura de ola última y elevación del casco asociados a los índices de confiabilidad máximos
βo MAX = 4.16 βo MAX = 4.00 βo MAX = 3.9
Localización RSR = 2.57 RSR = 2.27 RSR = 2.11
HU (m) ECMIN (m) HU (m) ECMIN (m) HU (m) ECMIN (m)
L-1 21.66 15.72 20.35 14.78 19.61 14.24
L-2 24.05 17.46 22.60 16.41 21.77 15.81
Estimación de los índices de confiabilidad Los resultados indican que en L-2 (tabla 3), la au-
anual por riesgo toelevable debe tener, al menos, un índice de confiabili-
dad ambiental anual de 4.16 y un RSR de 2.57, con una
Las figuras 3 y 4 muestran resultados para la Ec. 3. Se elevación mínima del casco de 17.46 m.
observa que al aumentar el costo o la utilidad por la
renta de la autoelevable, el índice de confiabilidad total Modelo de confiabilidad estructural
crece; pero, si el ΔCi aumenta, dicho índice se reduce. Se
encuentra que el índice de confiabilidad total óptimo Se emplea el índice de confiabilidad de Cornell para ob-
anual tiene un valor mínimo βot MIN = 2.72 (figura 3a) y tener el índice de confiabilidad estructural, como capa-
un máximo βot MAX = 3.83 (figura 4b). cidad disponible, de las autoelevables ante huracanes.
Para estimar la probabilidad de falla ambiental se Se consideran tres estados límite: el volteo de la au-
emplean factores de reducción (η) de 0.25, 0.50 y 0.75 toelevable, el hundimiento o deslizamiento horizontal
(Bea, 1997), con la Ec. 5 se calculan sus correspondien- de algunas de sus zapatas.
tes índices de confiabilidad y con la Ec. 11 el RSR a exi- Así, si Z es la función de estado límite que relaciona
gir (tabla 1). la resistencia R y la solicitación S, tal que
La tabla 2 muestra la profundidad, la H100 y la ECMIN
para el caso de los índices de confiabilidad mínimos z = r − s (13)
asociados a las dos localizaciones. La tabla 3 hace lo
propio para los índices de confiabilidad máximos. el índice de confiabilidad (β) se calcula como
6 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
4.0 4.0
ΔCi = 2 ΔCi = 2
3.5 ΔCi = 4 3.5 ΔCi = 4
ΔCi = 6 ΔCi = 6
b ot b ot
ΔCi = 8 ΔCi = 8
2.5 2.5
100 150 200 250 300 350 100 150 200 250 300 350
Costo de Jackup (MMUSD) Costo de Jackup (MMUSD)
a) ϑ= 0.05% b) ϑ= 0.30%
Figura 3. βot sin pérdida de la plataforma fija (γ= 0)
4.0 4.0
ΔCi = 2 ΔCi = 2
3.5 ΔCi = 4 3.5 ΔCi = 4
ΔCi = 6 ΔCi = 6
b ot ΔCi = 8 bot ΔCi = 8
ΔCi = 10 ΔCi = 10
3.0 3.0
ΔCi = 12 ΔCi = 12
2.5 2.5
100 150 200 250 300 350 100 150 200 250 300 350
Costo de Jackup (MMUSD) Costo de Jackup (MMUSD)
a) ϑ= 0.05% b) ϑ= 0.30%
Figura 4. βot con pérdida total de la plataforma fija (γ= 1)
PF = Φ (‒β) (15) Este colapso ocurre cuando una de las zapatas empieza
a trabajar a tensión, lo que se interpreta como un levan-
Como la ocurrencia de huracanes sigue un proceso de tamiento y separación de la zapata del suelo. La figura
Poisson con una tasa media de ocurrencia v, entonces la 5 ilustra esta falla con la gráfica que muestra inicial-
probabilidad de falla anual (Pf) de la plataforma se esti- mente a las tres piernas con una reacción de compre-
ma como sión similar, pero luego, al incrementarse las cargas
laterales de oleaje, dos piernas, las designadas como A
Pf = v PF (16) y B, incrementan su reacción a compresión, mientras
que en la otra, la C, la reacción vertical se reduce hasta
y su índice de confiabilidad anual ambiental disponible llegar a cero.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 7
Seguridad estructural de plataformas autoelevables petroleras en operación
Zapata C
Zapata B
Zapata A
Factor 2.22
C
A
B
Figura 5. Falla por volteo de plataforma autoevaluable
Una alternativa para el cálculo del momento de volteo tado de mar (Oj) es igual a la varianza del momento de
resistente la proporciona Det Norske Veritas (DNV) volteo (s s2 )
(2012), con la que resultan valores conservadores, los
cuales no se emplearon aquí. n
s s2 = ∑ (O j s s2, j )
La función de estado límite por volteo se aplica por j =1 (20)
8 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 9
Seguridad estructural de plataformas autoelevables petroleras en operación
Respecto a la evaluación de la media y desviación es- identificar un número suficiente de tormentas en una
tándar de la carga se presentan los dos enfoques que se base de datos de largo plazo. Para cada tormenta, el
explican a continuación. mayor valor de altura de ola significante (HS) es el valor
pico sobre el umbral, útil para describir su distribución
probabilista. La figura 7 muestra un registro de huraca-
Cargas ambientales en cada instante de una tormenta
nes correspondiente a 29 tormentas en 101 años (OWI,
En este caso, la media de la respuesta provocada por 2006), que señala el umbral inferior, valores de picos
oleaje es nula. Esto se observa en la figura 6, donde la locales y valores extremos seleccionados.
media de la reacción vertical por oleaje para un apoyo, La teoría extremos indica que la distribución de un
tiende a cero. máximo en una secuencia de valores arriba de un um-
Entonces, la media se puede calcular de la respuesta bral se generaliza mediante una distribución de Pareto
estática determinista de la estructura ante cargas de (Holthuijsen, 2007). DNV (2010) sugiere emplear una
viento (RV) y corriente (RC) distribución Weibull de dos parámetros o una exponen-
cial. Sin embargo, en el presente modelo se emplea la
R = RV + RC (27) distribución de Weibull de tres parámetros, cuya acu-
mulada se muestra por
La mayor incertidumbre en la carga la provoca el olea-
je, por ello se incluye únicamente la desviación estándar h − A C
de la respuesta por oleaje. Para un estado de mar la 1 − exp − S
FHS ( hs) = (32)
desviación estándar de la respuesta es B
10 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
5
4
3
2
1
Threshold 0.1
Figura 7. Valores pico de altura de ola
0
0 5000 10000 15000 20000 25000 30000 35000
significante (m) para huracanes en L-2
Date (OWI, 2006)
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 11
Seguridad estructural de plataformas autoelevables petroleras en operación
0.75 0.3
FHs(hS)
fHs(hS)
0.5 0.2
L-1 L-1
0.25 L-2 0.1 L-2
0 0
0 5 10 15 0 5 10 15
hS (m) hS (m)
Figura 9. Distribución de probabilidades de los huracanes en ambos lugares
a) Reacción vertical en zapatas para la dirección 240° b) Momento de volteo resistente (ft-kips)
momento resistente [ft-kips]
270°
300° 240°
2,038,098.43
1,212,274.33
330° 210°
1,060,345.87
1,449,487.03 1,239,549.25
3 0° 180°
1,341,203.83
A 2,217,321.46
Factor de carga
1,430,797.83
2 A
B 1,262,521.39
B
C 30°
1,090,282.86
150°
1 factor 2.22
1,235,253.09
2,070,385.65
0 60° 120°
-5000 0 5000 10000 15000 20000
Reacción vertical (kips) 90°
Cálculo de los índices de confiabilidad estructural Para estas aplicaciones el modo de falla dominante
es el deslizamiento de una zapata. Además, se aprecia
La confiabilidad para cada estado límite se calcula para que a mayor profundidad menor será el índice de con-
la dirección indicada en la tabla 5. Para cada estado de fiabilidad para cada estado límite, apreciación que se
mar se calcula el valor medio de la carga con análisis manifiesta en la localización L-2, con mayor profundi-
estático y la deviación estándar mediante análisis diná- dad, donde la confiabilidad de la plataforma autoeleva-
mico. Luego, se obtuvieron los índices de confiabilidad ble es menor para todos los casos analizados.
con la Ec. 14, cuyos resultados para la L-2 se muestran El índice de confiabilidad estructural debe ser anual
en las tablas 6 y 7; el índice más bajo fue de 5.6, corres- para compararse con las exigencias del modelo de ries-
pondiente a la falla por deslizamiento horizontal de go. Mediante las Ecs. 16 y 17 se calculan la probabilidad
una zapata. El mismo proceso se realiza para la L-1-. La de falla anual e índice de confiabilidad anual respecti-
figura 13 muestra los índices de confiabilidad de las au- vamente, para los índices de confiabilidad más bajos
toelevables para ambas localizaciones. presentados en la figura 13, asociados al deslizamiento
12 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
horizontal. Para la ocurrencia de huracanes se empleó ambas localizaciones con valores mayores respecto a
la tasa de excedencia v = 0.1 asociada al periodo de re- las exigidas por el modelo de riesgo, por ello es adecua-
torno mínimo de 10 años. La figura 14 muestra la con- do el empleo de la autoelevable.
fiabilidad anual, en capacidad de las autoelevables en
300° 240°
máxima reacción horizontal [kips]
12,537.49 270°
13,044.61
300° 240°
330° 210°
785.42
801.24
11,142.28
12,419.86 12,307.28 330° 210°
781.15
C
802.40
790.89
12,829.30
0° A 180° C
10,959.44
828.50
0° 180°
824.05
A
B
12,015.23 12,350.19
11,055.22
791.63 B
30° 150° 770.83
813.21
30° 150°
784.03
12,190.39
12,660.84
823.47
60° 120° 60° 120°
90° 90°
Figura 11. Reacción vertical (kips) de zapatas en L-1 Figura 12. Reacción horizontal (kips) de zapatas en L-1
Tabla 6. Confiabilidad en L-2 para las cargas ambientales en todo instante de las tormentas
Modo de falla Variable Distribución m s b
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 13
Seguridad estructural de plataformas autoelevables petroleras en operación
Tabla 7. Confiabilidad en L-2 para las cargas ambientales asociadas a las crestas de las tormentas
14 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
Campos Dante, Cabrera-Miranda J. Manuel, Martínez-Mayorga J. Manuel, López-Cornejo Marcelino, Espinosa-Zamudio E. Ernesto
Modelo de confiabilidad estructural metoceánico partment of Civil and Environmental Engineering, University
of California at Berkeley, March, 1997, pp. 96.
Se mostró un modelo de confiabilidad que permite el cál- Bentley. SACS: Structural Analysis Computer System, Bentley Insti-
culo del índice de confiabilidad estructural de la autoele- tute, 2012.
vable, como capacidad, dada su selección. Se requiere Campos D., Martínez-Mayorga J.M., Cabrera-Miranda J.M., Ló-
tomar en cuenta todas las características estructurales de pez-Cornejo M., Tenorio-Jaramillo C.A., Espinosa-Zamudio
la autoelevable y solicitaciones metoceánicas. E. Seguridad estructural de jackups durante su operación con
Se consideraron tres estados límite: el hundimiento base en estudios de riesgo y confiabilidad, XVIII Congreso
de una zapata, el deslizamiento lateral de un apoyo y el Nacional de Ingeniería Estructural, SMIE, Acapulco, 01-03 no-
volteo de la autoelevable. La carga actuante fue la me- viembre, 2012.
toceánica y la capacidad de la autoelevable se aplica por Campos D., Soriano A., Alamilla J.L., Ortega C. Aplicaciones del
la cimentación. Los análisis estructurales proporciona- riesgo y confiabilidad estructural en la norma mexicana de
ron la información necesaria para la aplicación del mo- diseño y evaluación de plataformas marinas fijas, 1er Simpo-
delo de confiabilidad. sio Análisis, diseño, evaluación, inspección y mantenimiento
Se realizó una aplicación con un tipo de plataforma de plataformas marinas, Aspectos relevantes, IMP y PEP,
autoelevable en dos localizaciones, lo que permitió ob- México, 18-19 octubre, 2010.
servar que a mayor profundidad el índice de confia- DNV. Self-Elevating Units, DNV-RP-C104. Norway, 2012.
bilidad disminuye y que el modo de falla dominante es DNV. Environmental Conditions and Environmental Loads,
el deslizamiento lateral de la cimentación, con un índi- DNV-RP-C105. Norway, 2010.
ce de confiabilidad anual de 5.99. Holthuijsen L.H. Waves in oceanic and coastal waters, Cambridge
University Press, 2007.
De ambos Kaiser M J, Snyder B F. Reviewing rig construction cost factors.
Offshore, Vol. 72, Issue 7, 2012.
Al aplicar ambos modelos se puede definir si una plata- Lind N.C., Davenport A.G. Towards practical applications of
forma autoelevable se puede emplear satisfactoriamen- structural reliability theory, en: Probabilistic design of rein-
te en determinado lugar con condiciones metoceánicas forced concrete buildings. ACI Publication, SP-31, 1972: 63-110.
intensas. Noble Asset México LLC. Noble Eddie Paul – Operations Manual.
En los casos de aplicación numérica, el máximo va- Revision 3. January 19, 2007, pp. 373.
lor de la confiabilidad anual del estudio de riesgo es Offshore Magazine. Worldwide survey of deepwater jack-up rigs.
4.16 (tabla 1), que es menor al índice de confiabilidad Drilling Rigs Capable of Operating in ≥ 350’ WD, 2011.
anual más bajo, obtenido por capacidad, que es 5.99. OWI. Update of metocean design data for zona norte and sonda
Esto significa que la plataforma autoelevable elegida de Campeche, final Report for IMP (Revised), Julio 18, 2006.
tendrá un buen desempeño en ambas ubicaciones. PEMEX. Diseño y evaluación de plataformas marinas fijas en el
Esto mismo significa que las capacidades en térmi- Golfo de México, NRF-003-PEMEX-2007, México, 2008, 82 p.
nos del factor de reserva de resistencia a obtener a par- PEMEX. Diseño y Evaluación de Plataformas Marinas Fijas en el
tir de los análisis de confiabilidad también estarán por Bahía de Campeche, NRF-003-PEMEX-2000, México, 2000,
arriba de los mínimos exigidos por el estudio de riesgo. 58 p.
Es importante mencionar que también la altura de ola PEMEX. Plan de respuesta a emergencias en la sonda de Campe-
que podrá soportar la plataforma será mayor que la mí- che (PRESC), México, 1997.
nima exigida. Sharples M. Post Mortem Failure Assessment of jack-up MODUs
During Hurricanes Katrina & Rita. Minerals Management
Referencias Service Report, 2008.
Stahl B. Reliability engineering and risk analysis, en: McClelland
Bea R.G. Risk based oceanographic criteria for design and requalifica- B., Reifel M.D. (Eds.), Planning and design of fixed offshore plat-
tion of platforms in the bay of Campeche. Report to PEMEX and forms, Van Nostrand Reinhold, Nueva York, 1986, pp. 59-98
IMP by Marine Technology & Management Group, De- (Chapter 5).
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM 15
Seguridad estructural de plataformas autoelevables petroleras en operación
16 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 1-16 ISSN 1405-7743 FI-UNAM
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 17-29
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: mayo de 2015, reevaluado: junio de 2016, aceptado: septiembre de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Jiménez-Aviles Gerardo
Universidad Autónoma del Estado de México
Facultad de Química
Departamento de Posgrado
Correo: gerjimav@[Link]
Valdivia-Barrientos Ricardo
Instituto Nacional de Investigaciones Nucleares
Departamento de Estudios del Ambiente
Laboratorio de Aplicaciones de Plasmas
Correo: [Link]@[Link]
Colín-Cruz Arturo
Universidad Autónoma del Estado de México
Facultad de Química
Departamento de Posgrado
Correo: acolinc@[Link]
Pacheco-Pacheco Marquidia
Instituto Nacional de Investigaciones Nucleares
Departamento de Estudios del Ambiente
Laboratorio de Aplicaciones de Plasmas
Correo: [Link]@[Link]
Pacheco-Sotelo Joel
Instituto Nacional de Investigaciones Nucleares
Departamento de Estudios del Ambiente
Laboratorio de Aplicaciones de Plasmas
Correo: [Link]@[Link]
Resumen
Este trabajo describe el desarrollo e implementación de un sistema de monitoreo y control de los parámetros físicos
involucrados en las etapas de inicio, sostenimiento y apagado de descargas de barrera dieléctrica generadas dentro
de un reactor de geometría coaxial. En el proceso se empleó el diseño gráfico de sistemas (GSD, graphical system de-
sign) para conjuntar eficientemente una interfaz gráfica de usuario (GUI, graphical user interface), la tarjeta de adqui-
sición de datos DAQ-6259 con canal de comunicación USB, una fuente de voltaje conmutado push-pull, controladores
de flujo gaseoso y equipo de diagnóstico químico y electrónico. El núcleo del enfoque GSD se encuentra en el código
de programación G de LabVIEW™; de esta manera, se recopilaron datos precisos y en tiempo real de las variables
(vgr.: flujo de aire y potencia aplicada), mismos que permitieron realizar un análisis del desempeño de las descargas.
Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
Durante las pruebas se emplearon diferentes concentraciones de una mezcla de helio contaminada con óxidos de
nitrógeno, con la finalidad de cuantificar el porcentaje de degradación/destrucción de los últimos. Los resultados
indican un valor mayor a 88%. Por seguridad, a la GUI se accede a través de la intranet institucional, debido a que
los compuestos generados pueden representar un riesgo para la salud humana.
Descriptores: contaminación atmosférica, descargas de barrera dieléctrica, diseño gráfico de sistemas, fuentes móvi-
les, óxidos de nitrógeno, plasma frío.
Abstract
This paper describes the development and implementation of a control and monitoring system for the physical parameters in-
volved in the stages of initiation, maintenance and shutdown of dielectric barrier discharges generated within a coaxial geometry
reactor. In the process, the graphical system design (GSD) was employed to efficiently combine a graphical user interface (GUI),
the data acquisition card DAQ-6259 with USB communication channel, a switch-mode power supply in push-pull configura-
tion, drivers for gas flow, chemical and electronic equipment for diagnosis. The core of the GSD approach is located in the G code
from LabVIEW ™, in this way, accurate data was collected in real time (e.g., airflow and power applied), which allowed to per-
form an analysis of the discharges performance. During testing, different concentrations of a mixture of helium contaminated
with nitrogen oxides, in order to quantify the percentage of its degradation. The results indicate a value greater than 88%. For
safety reasons, the GUI can be accessed remotely through institutional intranet, because of the generated compounds that may
represent a risk to human health.
Keywords: atmospheric pollution, dielectric barrier discharges, graphical system design, mobile sources, nitrogen oxides, cold plasma.
18 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo, Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
(gas oxidante), pues necesita funcionar a temperatura ganan electrón, X‒ átomo que gana electrón, X+ átomo
baja porque las moléculas de O3 decrecen rápido a tem- que cede electrón X2− y moléculas diatómicas excitadas.
peraturas elevadas y, al mismo tiempo, requiere una Las DBD, descargas de barrera o silenciosas, son un
presión relativamente alta (Gottschalk et al., 2010). Del conjunto de canales de conducción (microdescargas) en
mismo modo, sistemas en serie NTP-catalizador se han el espacio existente entre dos electrodos, separados por
combinado para tratar óxidos de nitrógeno (NOX) pro- una o varias capas dieléctricas (vidrio, sílice, polímeros o
venientes de automóviles: comprenden una etapa oxi- cerámicos), sobre los cuales se aplica voltaje alterno para
dativa (NTP) que convierte monóxido de nitrógeno iniciar y sostener las microdescargas. Las geometrías tí-
(NO) a dióxido de nitrógeno (NO2), y una fase de re- picas en los reactores son: de placas paralelas (planar o
ducción y almacenamiento (catalizador), que convierte coplanar) o cilíndricas (coaxial) (Kogelschatz, 2003).
el NO2 a nitrógeno (N2), dióxido de carbono (CO2), hi- Pacheco et al. (2008; 2012; 2013) reportan investiga-
drógeno (H2) y oxígeno (O) (Penetrante et al., 2002). ciones teórico-experimentales de tratamiento para ga-
Dado que las DBD operan a niveles altos de potencia y ses de combustión mediante DBD y se puntualiza el
tratan flujos grandes de gas a presión atmosférica con interés de emplearse a bordo de fuentes móviles (vehí-
caída de presión insignificante (basado en la tecnología culos). En su trabajo de 2008 realizaron el análisis de
madura de O3), potenciales aplicaciones en el control de cinética química, obtuvieron el mecanismo de degrada-
la contaminación atmosférica se han examinado siste- ción de los NOX y determinaron el comportamiento de
máticamente y en el presente trabajo se presenta la ins- las especies en la descarga. Los reactores utilizados fue-
trumentación basada en el GSD, de un reactor para ron coaxiales, de una y dos paredes de pírex, con bre-
DBD aplicado al tratamiento de NOX. cha de 1.4 cm y volumen para descarga de 65 µm3. En la
investigación de Estrada et al. (2011), con un consumo
Descargas de barrera dieléctrica de 20 W, se trató la mezcla de He (gas de acarreo), aire,
vapor de agua y NOX en un reactor de pírex, logrando
En el NTP la Te es el factor determinante al iniciar las obtener valores de remoción de 99.74% y 96.66% para
reacciones químicas, que en técnicas activadas térmica- niveles de concentración inicial de NOX de 50 y 300
mente resultarían ineficientes o imposibles de realizar. ppm, respectivamente.
Una vez iniciada la ionización, la distribución de ener- En las diversas aplicaciones de las DBD es importan-
gía la ejecutan los electrones (e‒), más ligeros que te conocer su comportamiento para poder cuantificar y
cualquier partícula y con alta energía cinética adquiri- analizar los diversos fenómenos que ocurren durante el
da, por lo que pueden acelerarse eficazmente para re- proceso por acción de las variables involucradas (flujo,
partir la energía al resto del gas mediante colisiones, voltaje aplicado o corriente demandada). Existen diver-
dentro de las cuales se encuentran: las elásticas, que sos dispositivos para realizar esta tarea, que van desde
cambian la energía cinética de las especies neutras, y las sensores individuales hasta controladores automatiza-
inelásticas, que convierten a las especies neutras en es- dos muy sofisticados (Jiménez et al., 2012).
pecies reactivas e inducen una serie de reacciones quí-
micas (tabla 1) que producen simultáneamente
Diseño gráfico de sistemas
portadores de carga y especies altamente reactivas en el
NTP (Meichsner, 2010). En (tabla 1) X2 representa las El GSD es una técnica de diseño (embebido), prototipo,
moléculas diatómicas, X átomos, X2+ moléculas diatómi- prueba, control y despliegue de aplicaciones, que com-
cas que ceden electrón, X2− moléculas diatómicas que bina programación gráfica en paralelo con hardware
COTS (commercial off-the-shelf) para obtener la mayor
Tabla 1. Reacciones en la Región Activa del NTP
cantidad de cálculo útil con un costo comercial bajo y
Ionización X2+ e− → X2+ + 2 e − simplificar el desarrollo de una aplicación, dando como
Captura de electrones X2+ e− → X2− resultado diseños robustos personalizados. El GSD per-
Captura disociativa de electrones X2+ e −
→ X− + X mite diseñar sin la necesidad de hacer uso del hardware
Disociación X2+ e− → X + X + e− final, posibilitando la realización de simulaciones y
− acortando los ciclos de integración e implementación
Excitación X2+ e → X2 + e −
(Bishop, 2007). En este proceso de instrumentación se
Ionización disociativa X2+ e− → X + + X + 2e −
emplearon:
Liberación de electrones X + 2e
−
2
−
→ X2 + 2 e −
Recombinación X + 2e
+
2
−
→ X2 1) LabVIEW™ 2013 para programar una GUI de un
Datos extraídos de Meichsner (2010) instrumento virtual (VI, virtual instrument);
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM 19
Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
A partir de los estudios realizados por Pacheco et al. GF = factor de corrección de la válvula
(2012), que establecen la composición y concentracio-
nes de una muestra sintética de gas de combustión, en Fuente de voltaje
este trabajo se usaron dos gases adquiridos en INFRA®:
la mezcla CAM 96 Middle Range (CAM), compuesta de Para generar y mantener las DBD es necesario alimentar
propano (C3H8), monóxido de carbono (CO), CO2 y NO, continuamente el sistema con una fuente bipolar, en caso
balanceados en N2 y Helio Industrial (He). Para sumi- contrario, el medio retorna a la fase gaseosa. La amplitud
nistrar la mezcla gaseosa se instalaron dos válvulas mínima de voltaje para que sucedan las descargas se de-
controladoras de flujo másico (MFC, mass flow contro- nomina Voltaje de Rompimiento (Vr), cuando su valor es
ller), una MKS™ MFC-M1000B (MFC1) para la CAM, y inferior al requerido la ionización es prácticamente nula,
una Omega™ FMA-A2315 (MFC2) para el gas de aca- una vez superado ligeramente ese nivel comienza a exis-
rreo, en este caso el He; ambas válvulas fueron conecta- tir un flujo de corriente entre los electrodos (plasma par-
das a un controlador MKS™ 247D. Debido a que los cialmente ionizado). El valor de Vr depende del gas, la
MFC se encuentran calibrados para flujos de N2, se ajus- presión de trabajo (P) y de la distancia entre los electro-
taron las escalas para el gas correspondiente con ayuda dos (d), parámetros vinculados con la Fórmula de Pas-
de las Ecuaciones 1-4 (MKS, 2013a; MKS, 2013b; Ome- chen (Ecuación 5, en la que δ y σ son constantes
ga, 2013). Los ajustes realizados se agrupan en la tabla específicas del gas inyectado) (Fridman, 2008).
2, donde al sumar los valores de la última columna se
determina que el flujo máximo de la mezcla He-CAM
δ pd
es de 67.75 l min-1, en una relación porcentual 85.84- Vr ( pd) = (5)
14.16%. σ + ln( pd)
20 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo, Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
1 (6)
Figura 1. Circuito eléctrico del inversor push-pull fR
2π C PP LPP
Si al C-DBD2 se acopla una fuente de voltaje alterno de
alta frecuencia, tanto la eficiencia como la conductivi- Reactor
dad eléctrica aumentan y el tamaño del equipo final se
reduce notablemente. La importancia de la frecuencia De acuerdo con la figura 2, las dimensiones y caracterís-
de excitación radica en su influencia en el comporta- ticas físicas del C-DBD2 son: 15.40 mm de radio exterior
miento general de la descarga, donde en el orden de (R1), 12 mm de radio interno (R2), 19.30 cm de longitud
kHz los electrones y los iones siguen las oscilaciones (D1), 4.85 mm de radio del aislante interno (R3); el volta-
del campo eléctrico, mientras que en un orden superior, je de salida de PPI (VPPI) se conecta a HV, que es un tubo
MHz o GHz, solo los electrones siguen dicho campo hueco de acero inoxidable de 13.62 cm de longitud, 1.57
(Roth, 1995). Una fuente conmutada (SMPS, switch-mo- mm de grosor y radio de 3.16 mm (R4); el electrodo
de power supply) permite trabajar con frecuencias de GND es una malla de acero inoxidable de 0.80 mm de
kHz mediante el cambio de estado de transistores. En grosor y 11.37 cm de longitud (D2). El volumen de las
este trabajo se ocupa como fuente principal de voltaje microdescargas es de 50.60 µm3. El reactor se elaboró
(VPP) la SMPS B&K® PRECISION XLN10014, que pue- con dos capas aislantes de alúmina, ya que propician
de proporcionar una potencia máxima de 1.44 kW (100 una descarga más homogénea respecto a reactores con
V @ 14.4 A), el diseño del sistema de excitación se enfo- una barrera.
ca en desarrollar un inversor con topología push-pull
(PPI, push-pull inverter) con base en la teoría expuesta Tarjeta de adquisición de datos
por Bersani (2009) y Mohammad (2007) (figura 1).
En el diagrama eléctrico del PPI se ubican los MOS- La tarjeta 6259 cuenta con 128 conectores de entrada/sa-
FET Q1 y Q2, matrícula IXFN80N50 (VDS(MÁX) = 500 V, lida en tiempo real, que incluyen 36 puertos analógicos:
IDS(MÁX) = 80 A y RDS(MÁX) = 55 mΩ), CPP es un condensador 4 salidas de 16 bits (AO, ±10 V) @ 2.00 MS s-1, y 32 entra-
de poliéster de 3.2 nF y forma el circuito taque resonan- das (AI, ±10 V) @ 1.00 MS s-1, que pueden referenciarse a
te con el transformador TPP, compuesto por las bobinas una tierra común (RSE, referenced single-ended) o estar ais-
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM 21
Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
22 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo, Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
3 2 4 5 6 11 3
10
Figura 4. Panel principal de la interfaz gráfica de usuario para control del reactor C-DBD2
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM 23
Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
24 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo, Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
Además, el voltaje visto por los transistores se eleva has- valores de la densidad volumétrica de potencia (DVP: po-
ta 4.5 veces el suministrado (figura 8), cuando idealmen- tencia distribuida por unidad de volumen disponible
te debería ser de 2.6 veces. La potencia de entrada se dentro del reactor), y de la Densidad de Energía de En-
obtuvo al multiplicar los valores de la corriente y el vol- trada (DEE: cantidad de energía aplicada por volumen
taje obtenidos con la interfaz virtual, mientras que la po- de gas). Con ambos conceptos, DVP y DEE, se constru-
tencia de salida se calculó mediante la Fórmula de ye el concepto Valor Fuente que permite visualizar las
Potencia para Ozonizadores de Manley (1943) a través cantidades de potencia y energía necesarias por unidad
de un algoritmo en MathCAD®; con la razón de ambas de volumen (µm3 y l) para cada dato de degradación
potencias, se obtuvo que el índice de eficiencia del siste- obtenido.
ma de potencia es en promedio del 0.90, mismo que varía Al aplicar 25 V y hacer fluir He se identifican oscila-
al filamentarse las descargas. ciones súbitas de 0 ‒ t1 en la gráfica de la potencia (debi-
Las condiciones mínimas para que inicien las descar- das a la corriente, ya que el voltaje se mantiene
gas son: He = 1.236 l min-1, VPP = 19 V, VPPI(pk) = 1.88 kV @ constante), suavizadas inmediatamente al incluir CAM
3.7 W; con ello, el sistema puede trabajar en estado esta- en 0.38 %vol; en este primer caso se consideraron las
cionario con un flujo reducido de 0.15 l min-1 de He y un muestras de 0.38 y 1.14 %vol. de CAM, que trabajaron
consumo de 3.2 W. En la figura 9 se observan las formas en estado estacionario con un consumo promedio de
de onda del voltaje aplicado, la corriente y la potencia de- 7 W. Posteriormente se aumentó la concentración de
mandadas; al incrementar gradualmente el flujo decae la CAM y la potencia decayó instantáneamente con la
corriente y las descargas se extinguen (zona de standby). consecuente extinción total de las microdescargas. Para
Posteriormente se hizo fluir una mezcla de He-CAM el caso de 35 V, se repiten las concentraciones anteriores
a las concentraciones y flujos enlistados en la tabla 4, y se nota un ascenso de la potencia sin que esto refleje
con potenciales de 25 a 55 V; cada una de las concentra- un aumento proporcional de degradación.
ciones para la fila de 55 V corresponde a los tiempos tx El procedimiento se repite para 45 y 55 V. En la figu-
identificados en la figura 10, también se entregan los ra 10, de t1 ‒ t2 la rampa de potencia aumenta 8 W y la
6 25 6 25 7 30
20 25
20
4 15 4 20
Corriente de 15 Corriente de 15
Corriente de 3.5
10 10 Resonancia
Corriente (mA)
Corriente (mA)
Corriente (mA)
2 2 Resonancia 10
Resonancia 5 5 5
VPPI (kV)
VPPI (kV)
VPPI (kV)
0 0 0 0 0 0
-5 -5 -5
-2 -10 -2 -10
-10
-3.5 -15
-15 -15
-4 -4 -20
-20 Fase -20 Fase -25
Fase
-6 -25 -6 -25 -7 -30
0 10 5 15 20 25 5 10 15 20 0 0 5 10 15 25 20 25
Tiempo (ms) Tiempo (ms) Tiempo (ms)
Figura 7. Variación de fPWM y DPWM para determinar la transferencia máxima de potencia, a) valores correctos, b) fPWM= 56 kHz y
DPWM=25.5%, c) fPWM=56 kHz y DPWM= 20%
Figura 8. Detección de valores pico para estimar el voltaje Figura 9. Identificación de las condiciones mínimas de potencia
máximo de operación del MOSFET para iniciar las descargas
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM 25
Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
descarga se vuelve filamental y heterogénea, fenómenos co no refleja una degradación mayor para VPP = 35 y
debidos a la elevada DVP y la baja concentración de 45V, a pesar de que la DVP y la DEE se intensifican 2 y
CAM. En el lapso t11 ‒ t13 se notan oscilaciones incorpora- 4 veces su valor inicial, respectivamente. Se identifica
das en la gráfica, derivadas de la filamentación homogé- que el porcentaje de degradación aumenta en cada ex-
nea de las microdescargas. En t12 la CAM está presente a perimento para la misma muestra conforme se incre-
13.95%vol. y se eleva a 15.61%vol., el sistema permanece menta VPP para flujos mayores a 0.14 l min-1 de CAM
por un instante hasta que la potencia decae naturalmente (1.14 %vol., 24 ppm de NOX); en estos casos sí es con-
en t13. Más allá de esta última concentración, las DBD gruente el aumento de la DVP y la DEE. El valor máxi-
tienden a localizarse dentro del reactor y a extinguirse. mo es de 74.07% de degradación cuando VPP = 55 V; sin
El valor máximo de degradación de NOX es de embargo, cuando el flujo de CAM es máximo (2.02 l
88.24% para una concentración inicial (CI) del 0.38% (9 min-1) únicamente alcanza un 19.27%.
ppm), en un flujo total de 12.48 l min-1 y con VPP = 25 V Con lo anterior, se culminaron las series de pruebas y
(DVP=141.50 mW µm-3 @ DEE=34.42 J l-1); este específi- se dispuso el equipo para su adecuación y mejoramiento.
26 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo, Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM 27
Control de un sistema de descargas de barrera dieléctrica basado en el diseño gráfico de sistemas
Conrads H. y Schmidt M. Plasma generation and plasma sources. Meichsner J. Nonthermal reactive plasmas, chapter 11, 267-298,
Plasma Sources Science and Technology, volumen 9 (número 4), en: Bonitz M. et al. eds., Introduction to Complex Plasmas, Sprin-
2000: 441. ger, 2010, ISBN 9783642105913.
Estrada N. et al. Chemical and electrical diagnosis of two configu- MKS (3 de enero de 2013). MKS Type 247D Four-channel readout,
rations of dielectric barrier discharges applied to nitric oxides Instruction manual [en línea]. Disponible en: [Link].
degradation. Superficies y Vacío, 2011, DOI: ND. com/docs/R/[Link].
Ferroxcube (3 de enero de 2013). U126/91/20. U, I Cores and Ac- MKS (3 de enero de 2013). MKS Type M100B Mass-flo controller
cessories, Datasheet [en línea]. Disponible en: [Link] and m10b mass-flo meter, Instruction Manual [en línea]. Dispo-
[Link]/FerroxcubeCorporateReception/datasheet/ nible en: [Link]/docs/R/[Link].
[Link]. Mohammad, K. Switch Mode Power Supply (SMPS) Topologies (Part
Fridman A. Plasma chemistry, Cambridge University Press, 2008, I), Microchip Technology, 2007, Inc., ID: AN1114.
ISBN 9780511398575. NI (2013, 5 de febrero). NI 625x Specifications [en línea]. Disponi-
Gottschalk C. et al. Ozonation of water and waste water, A practical ble en: [Link]
guide to understanding ozone and its applications, Wiley, 2010, Omega (3 de enero de 2013). FMA-A2100’s, FMA-A2300’s
ISBN 9783527319626. Massflow meters, user´s manual [en línea]. Disponible en:
Jiménez G. et al. Instrumenting and programming a virtual instru- [Link]/Manuals/manualpdf/[Link].
ment for an open loop system to control gliding centrifugal Pacheco M.J. et al. Application of non-thermal plasma on gas
thermal plasma, 9th CCE, 2012, DOI: 10.1109/ cleansing. Physica Scripta, 2008, DOI: 10.1088/0031-8949/2008/
ICEEE.2012.6421110. T131/014017.
Kogelschatz U. Dielectric-barrier discharges: Their history, dis- Pacheco M. et al. Removal of main exhaust gases of vehicles by a
charge physics, and industrial applications. Plasma Chemistry double dielectric barrier discharge, 14th LAWPP, 2012, DOI:
and Plasma Processing, 2003, DOI: 10.1023/A:1022470901385. 10.1088/1742-6596/370/1/012023.
Kunhardt E.E. Generation of large-volume, atmospheric-pressu- Pacheco M.J. et al. Tratamiento de emisiones de motores diesel por
re, nonequilibrium plasmas. IEEE Transactions on Plasma plasma fuera de equilibrio. Revista Internacional de Contamina-
Science, 2000, DOI: 10.1109/27.842901. ción Ambiental, 2013, DOI: ND.
Manley T.C. The electric characteristics of the ozonator discharge. Penetrante B.M. et al. Plasma-assisted catalytic storage reduction
Journal of the Electrochemical Society, 1943, DOI: 10.1149/1.3071556. system, Patente, 2002, US6038853 A.
Roth J.R. Industrial plasma engineering, volume 1, principles, Insti-
tute of Physics Publishing, 1995, ISBN 0750303174.
28 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM
Jiménez-Aviles Gerardo, Valdivia-Barrientos Ricardo, Colín-Cruz Arturo, Pacheco-Pacheco Marquidia, Pacheco-Sotelo Joel
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 17-29 ISSN 1405-7743 FI-UNAM 29
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 31-42
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: julio de 2015, reevaluado: febrero de 2016, aceptado: septiembre de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Torres-Lozada Patricia
Universidad del Valle
Facultad de Ingeniería
Correo: [Link]@[Link]
Resumen
El compostaje es una de las tecnologías más empleadas para el aprovechamiento de biorresiduos, no obstante, su implementación en paí-
ses en desarrollo no ha sido efectiva, debido, entre otros aspectos, a la limitada investigación para su aplicación. Este artículo presenta
aportes en torno a los avances en la investigación del compostaje de biorresiduos en municipios menores de países en desarrollo, con base
en resultados de seis años de estudios realizados por los autores. Se abordan: i) el análisis de la calidad fisicoquímica de los sustratos, ii) la
evaluación de opciones para mejorar el proceso y la calidad del producto, iii) el desarrollo de herramientas para la planeación y operación
de las instalaciones de compostaje. Estos avances han permitido reducir los tiempos de proceso, mejorar las condiciones para la higieniza-
ción del material, cumplir con los estándares de calidad del producto, identificar pruebas para mejorar el control y monitoreo in situ de la
estabilidad y madurez del producto y desarrollar herramientas para la selección de materiales de enmienda y de soporte. La reflexión
ayuda a elucidar que la investigación, acompañada de proyectos a escala piloto, se constituye en una estrategia para posicionar el compos-
taje de biorresiduos como una opción adecuada de manejo de la fracción orgánica de los residuos sólidos municipales (RSM).
Abstract
Composting is one of the most used technologies for biowaste recycling. However, its adoption has been ineffective in developing countries; one of the
causes being the limited research for implementation in these contexts. This paper presents a contribution around progress in research on biowaste
composting in small municipalities from developing countries, based on results from six years of studies conducted by the authors. Three topics are ad-
dressed: i) analysis of the physicochemical quality of substrates; ii) assessment of alternatives for improving the process and product quality; and iii)
development of tools for planning and operation of composting installations. This progress has allowed reducing processing times, improve the condi-
tions for material sanitization, fulfil the standards of product quality, identify tests for improving in situ control and monitoring of product maturity
and stability, and developing tools to select bulking and amendment materials. Research together with pilot-scale projects is a significant strategy for
positioning biowaste composting as an appropriate alternative for managing the organic fraction of municipal solid waste.
32 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM
Oviedo-Ocaña Edgar Ricardo, Marmolejo-Rebellon Luis Fernando, Torres-Lozada Patricia
nico y nutrientes), porosidad del sustrato, humedad, La localización, estilos de vida, hábitos de preparación
pH, oxígeno y temperatura del proceso, los cuales están de alimentos y condiciones socio-económicas, son factores
ampliamente documentados (Haug, 1993; Insam y de que afectan la cantidad y la composición de los biorresi-
Bertoldi, 2007; Krogmann et al., 2010). duos. Estos sustratos se caracterizan por una alta propor-
El producto del proceso se denomina compost y se ción de residuos de alimentos, tal como lo reportan
caracteriza por tener materia orgánica estabilizada, que Oviedo et al. (2015a) quienes encontraron en una instala-
es compatible con el crecimiento de plantas (Stentiford ción de compostaje en un municipio menor a Colombia,
y de Bertoldi, 2010) y tiene una incompleta humifica- que 92.9% (±4.3%) corresponde a alimentos sin procesar,
ción (Insam y de Bertoldi, 2007). Presenta propiedades entre los que se encontraban en mayor proporción pláta-
fisicoquímicas que determinan sus características y nos y tubérculos (32.6%) y frutas cítricas (15.5%). Estos
usos, entre las que se encuentran el carbono orgánico total residuos son altamente variables por su composición de
(COT), capacidad de intercambio catiónico (CIC), NTotal materia orgánica rápidamente degradable que lo hacen
y Ninorgánico, pH, conductividad eléctrica (CE), macro y mi- complejo de manejar (Iacovidou et al., 2012).
cronutrientes, densidad aparente, humedad, capacidad La evaluación fisicoquímica de los biorresiduos en
de retención de agua (CRA) y tamaño de partícula. Una esta población muestra que presentan condiciones críti-
descripción de estos parámetros se encuentra en Sulli- cas para el compostaje (tabla 1), tales como pH ácido,
van y Miller (2001) y Hargreaves et al. (2008). exceso de humedad, deficiencias en COT y PTotal y una
relación C/N inferior a la recomendada para el inicio
Compostaje de biorresiduos en municipios del compostaje. Estas características son similares a las
menores de países en desarrollo encontradas en otras poblaciones de Colombia (Mar-
molejo, 2011) y son relativamente distintas a las repor-
Identificación de las características de los biorresiduos tadas para biorresiduos en países desarrollados como
España y Francia, tal como lo muestran López et al.
Los biorresiduos son la fracción orgánica biodegrada- (2010), Huerta et al., (2011) y Francou et al. (2006) (e.g.
ble de los RSM y denotan los materiales separados en la contenido COT, NTotal, PTotal) (tabla 1).
fuente y recolectados selectivamente antes de introdu-
cirlos a la cadena de manejo de RSM. Provienen de resi- Evaluación de opciones para mejorar la calidad
duos de jardín y de alimentos de las viviendas y de de los sustratos y su efecto en el proceso
establecimientos comerciales. Se excluyen residuos y en la calidad del producto
agrícolas y forestales, estiércoles animales, biosólidos u
otros biodegradables (textiles, papel y madera) (EC, Debido a las limitaciones que presentan los biorresi-
2008). Mientras que en países desarrollados la propor- duos para un compostaje efectivo, es necesario evaluar
ción de biorresiduos en los RSM es del orden de 30 a opciones de acondicionamiento del sustrato.
40%, en países en desarrollo oscila entre 50 y 70% (Tros- Con el propósito de analizar la transformación de la
chinetz y Mihelcic, 2009). Estudios realizados en seis materia orgánica de este tipo de sustratos, los autores
municipios menores que Colombia, muestran que la desarrollaron un experimento en el que evaluaron la
proporción varía entre 62 y 65% (Marmolejo, 2011). influencia del pH y los nutrientes sobre la biodegrada-
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM 33
Avances en investigación sobre el compostaje de biorresiduos en municipios menores de países en desarrollo . Lecciones desde Colombia
bilidad aerobia de los biorresiduos (Soto et al., 2015). El Por otro lado, el pH ácido de los sustratos requiere
experimento se realizó a escala de laboratorio, em- controles operacionales como el incremento en la airea-
pleando técnicas respirométricas. Las condiciones más ción y la adición de material alcalino (Wong et al., 2009)
favorables para la cinética de transformación se presen- para minimizar el efecto de la acidez en el proceso; es-
taron cuando el sustrato se ajustó a pH entre 7 a 8 (i.e. tas opciones son efectivas en diversos estudios (Kuba et
se añadió NaOH a 50% de pureza) y conjuntamente se al., 2008; Wong et al., 2009; Kurola et al., 2011).
adicionaron nutrientes (i.e. se añadió 1 mL por cada 1 L La selección de materiales en instalaciones de com-
del reactor en forma de una solución de KH2PO4 – 8.5g, postaje requiere el desarrollo de criterios que conside-
K2HPO4 – 21.7g, Na2HPO4H2O – 33.4g, NH4Cl – 2.0 g), ren las características de los sustratos y condiciones
evidenciándose mayores tasas de biodegradación y de socioeconómicas y técnicas, tal como propuso Oviedo
crecimiento de biomasa. Esto reitera la necesidad de (2015), quien consideró los siguientes criterios:
acondicionar los sustratos para mejorar el desempeño
del proceso en las instalaciones de compostaje. i) Calidad del material: debe contribuir a mejorar las
El mejoramiento de la calidad de los sustratos en las condiciones de C/N, humedad, pH, otros nutrientes
instalaciones de compostaje, requiere propuestas de ope- y porosidad de los biorresiduos
ración, como la incorporación de materiales de soporte ii) Cantidad del material: debe estar disponible en las
(controlan la humedad, suministran porosidad y mejo- cantidades requeridas para emplearse en la instala-
ran la aireación) o de enmienda (mejoran las característi- ción de compostaje
cas químicas de sustratos) (Epstein, 2011). El co-com- iii) Costo de adquisición del material: debe tener un
postaje es otra medida que involucra dos sustratos con costo mínimo, de manera que no genere mayores
el propósito de mejorar la eficiencia del proceso y la cali- costos operativos en la instalación de compostaje
dad del producto, y el aditivo es una mezcla de microor- iv) Acceso para consecución del material: debe tenerse
ganismos, nutrientes minerales, formas de C fácilmente en cuenta el tipo de vía de acceso para recolectarlo
disponibles, enzimas y compuestos para balancear el v) Distancia para recolección del material: debe consi-
pH, que se utilizan para mejorar la actividad microbiana derarse la distancia desde la instalación de compos-
en los sustratos (Gabhane et al., 2012). taje hasta el centro de generación
En general, distintos materiales se han empleado vi) Disposición para entrega de material: debe estar dis-
para optimizar el compostaje de biorresiduos: Adhikari ponible la mayor cantidad de tiempo posible, para in-
et al. (2008) con paja de trigo y pasto picado, Francou et corporarse en el procesamiento de los biorresiduos
al. (2008) con residuos verdes y papel y cartón, Chang y vii) Facilidades de manejo del material: no deben gene-
Chen (2010) con cáscara de arroz y aserrín, y Kumar et rar mayores cambios operativos en la instalación de
al. (2010) con residuos de poda de jardín y hojas. Los compostaje (i.e. acondicionamiento del material en
materiales vegetales como la poda son fuente de N, actividades de selección y clasificación, secado y tri-
ayudan a regular la humedad y mejoran la estructura y turado). La tabla 2 presenta los criterios propuestos
aireación (López et al., 2010). El bagazo de caña de azú- para la selección de los materiales de enmienda y
car es un material lignocelulósico utilizado para mejo- soporte.
rar el COT en los sustratos (Iqbal et al., 2010).
En condiciones de acceso limitado a materiales de En el estudio de Oviedo (2015) se realizó un inventario
soporte o enmienda, el compost o el material en proce- de materiales de enmienda y soporte en un municipio
so (precompost) se emplean para la neutralización de menor que Colombia, identificando los materiales: resi-
ácidos y el suministro de especies microbiales que ace- duos de cosecha de frutas, bagazo de caña (BC), pasto
leran el proceso (Kato y Miura, 2008; Li et al., 2013). Ex- estrella (PE), material en proceso (MP) (precompost),
periencias en el compostaje de estiércol bovino (Kato y cenizas de madera (CM), estiércol bovino, estiércol por-
Miura, 2008) o de biorresiduos (Lu y Guo, 2009) han cino y pollinaza. La aplicación de la herramienta de se-
encontrado reducción en los tiempos de proceso. Aun- lección permitió escoger de entre los ocho materiales,
que la aplicación de precompost o compost en las insta- los cuatro siguientes: CM, MP, BC y PE, los cuales se
laciones de compostaje puede restringirse debido a la evaluaron en experimentos a escala piloto de compos-
necesidad de almacenar producto, se podría recircular taje de biorresiduos, como se indica a continuación:
material que se encuentra en proceso (precompost)
para mejorar el inicio del proceso de compostaje de bio- • Oviedo et al. (2014a) evaluaron la adición de CM
rresiduos, tal como se ha implementado en algunas ins- (proporción de 2, 4 y 8% en base húmeda definidas
talaciones de compostaje en el Valle del Cauca. con base en Kuba et al. (2008); Wong et al. (2009) y
34 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM
Oviedo-Ocaña Edgar Ricardo, Marmolejo-Rebellon Luis Fernando, Torres-Lozada Patricia
Kurola et al. (2011), la cual resultó favorable para masa alrededor de la humedad), identificando que
incrementar el pH inicial, mejorar el contenido nu- este contribuyó a mejorar la humedad inicial (i.e.
tricional de los sustratos (i.e. se incrementó conte- 72% en biorresiduos y 63.7% en la mezcla de biorre-
nido de KTotal y PTotal con 4% y 8%) y amortiguar los siduos con MP) y el COT de los sustratos (i.e. 31.2%
ácidos generados en la primera fase del proceso en biorresiduos y 36.9% en la mezcla de biorresi-
(i.e. tuvo valores al inicio de pH al inicio del pro- duos con MP), facilitando la obtención de mayores
ceso entre 7.2 y 9.0, mientras que sin adición de temperaturas en el rango termofílico respecto al
cenizas se presentaron 5.0 y 5.9); sin embargo, no compostaje de biorresiduos; no obstante, se presen-
tuvo efecto en acelerar el arranque del proceso y taron similares tiempos de duración de las etapas
se observó un excesivo incremento del pH con la termofílica y de enfriamiento y la obtención de un
adición de 4 y 8% de CM, que asociado con tempe- producto con características fisicoquímicas simila-
raturas termofílicas pudieron propiciar mayores res, reportándose únicamente diferencias significa-
pérdidas de NTotal. Las pruebas de germinación en tivas en la C/N y las cenizas (p<0.05).
los productos con CM mostraron presencia de fito- • Oviedo et al. (2015b) evaluaron en dos experimentos
tóxicos que pueden comprometer el uso del la incorporación del PE, en el primero (proporción
producto. de 34% base húmeda definida teniendo en cuenta
• Oviedo (2015) evaluó la incorporación de MP (pro- experiencias previas (Oviedo et al., 2013)) y del BC
porción de 25% estimada con base en balances de en el segundo (proporción de 22% base húmeda de-
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM 35
Avances en investigación sobre el compostaje de biorresiduos en municipios menores de países en desarrollo . Lecciones desde Colombia
finida a través de un balance de masa alrededor de buyen con el cumplimiento de estándares de calidad
la relación C/N), encontrando que en comparación del producto e incrementan la competitividad del com-
con las pilas con 100% de biorresiduos (control), se post de biorresiduos frente a otros productos de fertili-
aceleró el arranque del proceso, redujo la duración zación.
de la etapa termofílica y propiciaron condiciones Teniendo en cuenta que además de los biorresiduos,
adecuadas para la higienización del material (i.e. se las instalaciones de compostaje reciben material prove-
presentaron temperaturas en el rango termofílico niente del corte de césped y poda de árboles, se propu-
por más de dos semanas, tal como lo recomienda sieron otros experimentos que evaluaron a escala piloto
Böhm (2007) donde se reflejó en valores de colifor- el compostaje de residuos de alimentos no procesados
mes fecales de 0 NMP/g y 23 NMP/g en el primero y (RANP) y de poda de árboles y corte (RPAC) de césped.
segundo experimento, respectivamente); no obstan- Para tal efecto, se analizaron dos tratamientos con su
te, compuestos de difícil degradación en estos mate- réplica: tratamiento A (100% peso húmedo de RPAC) y
riales, pueden ser la causa de una mayor duración tratamiento B (40% peso húmedo de RANP y 60% de
en las fases de enfriamiento y maduración. La incor- RPAC). El sustrato del tratamiento B presentó mejores
poración del PE y el BC contribuyó a mejorar la cali- condiciones para el inicio del proceso (i.e. C/N: 20.1, hu-
dad de los productos en relación con el control medad: 74% y PTotal: 0.39%) respecto del tratamiento A
(tabla 3), obteniéndose niveles de COT, CIC, densi- (C/N: 14.8, humedad: 25.4% y PTotal: 0.28%). Esto se re-
dad y CRA que en su uso, pueden mejorar las pro- flejó en mejores condiciones de proceso, tales como me-
piedades del suelo, tales como la retención de nor tiempo para alcanzar las temperaturas termofílicas,
nutrientes y agua e incremento en la actividad mi- temperaturas más altas y de mayor duración para la
crobiana. Se observa además, un mayor cumpli- higienización del producto y menor duración de la fase
miento de la norma técnica colombiana. de enfriamiento (7 días) respecto del tratamiento A.
Asimismo, la calidad del producto se favoreció con
Los resultados en la calidad del producto guardan rela- contenidos superiores de nutrientes (Tratamiento B:
ción con lo reportado en estudios de países en desarro- NTotal: 2.1%, PTotal: 0.54% y KTotal: 3.1%; Tratamiento A:
llo (Barreira et al., 2006 en Sâo Paulo-Brasil; Adekunle et NTotal: 1.6%, PTotal: 0.45% y KTotal: 2.7%), que coinciden
al., 2010 en Abeokuta-Nigeria; Saha et al., 2010 en ciuda- con las mejoras en la calidad del producto reportadas
des de India), en los que el compostaje de biorresiduos en los estudios previos.
sin la incorporación de materiales (enmienda o sopor- Por otro lado, debido al alto contenido de humedad
te), se caracterizan por presentar altos valores de pH y identificado en los biorresiduos, Oviedo et al. (2014b)
CE, limitado contenido de COT o alto material inerte, evaluaron el efecto del incremento en la frecuencia de
restringiendo la comercialización del producto. Las es- volteo en el proceso de compostaje. En este estudio, re-
trategias de control presentadas en este estudio, contri- portaron la disminución en el tiempo del proceso con
36 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM
Oviedo-Ocaña Edgar Ricardo, Marmolejo-Rebellon Luis Fernando, Torres-Lozada Patricia
20% en las fases mesofílica, termofílica y de enfriamien- frecuencia en compost de biorresiduos son el Cd, Cr,
to, no obstante, se observó una disminución en el conte- Ni, Pb y Zn (Farrell y Jones, 2009). En cuatro experi-
nido nutricional en el producto. Los resultados de este mentos de compostaje de biorresiduos (Marmolejo,
estudio evidencian que los esquemas de control y mo- 2011; Oviedo et al., 2013; Oviedo et al., 2015b) se observó
nitoreo acordes con las condiciones locales, pueden ser que los productos cumplieron los límites de metales pe-
efectivos para mejorar el desempeño del proceso y sados de normas de Colombia, India, Reino Unido e
plantea nuevos retos en la investigación, tales como la Italia, lo cual se asocia al efecto de prácticas como la
evaluación conjunta del efecto del incremento de la fre- separación en la fuente.
cuencia del volteo y de la adición de materiales de so- La madurez y la estabilidad son otros de los crite-
porte (o enmienda) sobre el proceso y la calidad del rios para evaluar la calidad del producto en el uso de
compostaje. actividades agrícolas o mejorador de suelo; aunque
existen diferentes pruebas para evaluarlas, la aplicación
Consideraciones sobre la calidad del producto de pruebas de laboratorio en instalaciones de compos-
y potenciales usos
taje de municipios menores de países en desarrollo
puede ser limitada debido a costos y a la complejidad
Una de las condiciones críticas para determinar la cali- de su aplicación. En estas instalaciones se aplican prue-
dad de los productos y potenciar su uso y comerciali- bas de monitoreo in situ, con las que se toman decisio-
zación, está relacionada con la presencia de conta- nes operativas y de finalización del proceso. Oviedo et
minantes; un compost de calidad no se logra si los sus- al. (2015c) evaluaron pruebas empleadas para determi-
tratos contienen impurezas o componentes indesea- nar la estabilidad y madurez in situ (temperatura, pH,
bles, como metales pesados o compuestos orgánicos humedad, CE, olor y color) con pruebas de mayores
persistentes (Hargreaves et al., 2008; López et al., 2010). requerimientos de laboratorio (sólidos volátiles-SV, re-
El producto generado a partir de residuos mezclados, lación C/N, autocalentamiento, índice respirométrico-
puede ser de baja calidad y afectar la aceptación por IR e índice de germinación-IG). Los resultados
los agricultores, restringiendo su comercialización mostraron que las pruebas de monitoreo empleadas
(Zurbrügg et al., 2004). tradicionalmente in situ para evaluar la estabilidad
Esto se corroboró por los autores en el contexto de (temperatura, color y humedad) no fueron consistentes
un municipio menor de países en desarrollo, en el que con pruebas más complejas de laboratorio (IR, autoca-
se evaluó la influencia de la separación en la fuente so- lentamiento y SV). En relación con la madurez, se en-
bre el proceso y la calidad del producto del compostaje contraron resultados consistentes entre la prueba de
de biorresiduos (Marmolejo et al., 2010; Marmolejo, laboratorio (IG) con las pruebas in situ (pH y CE). Esto
2011). Este estudio mostró que a pesar de ser factible el muestra que en las instalaciones de compostaje se re-
proceso de compostaje de sustratos separados en la quiere el uso simultáneo de pruebas in situ con las de
fuente (SSF) y separados en planta (SSP), sin inhibición laboratorio. Estrategias de regionalización de laborato-
del proceso, el compostaje con los sustratos separados rios pueden ser efectivas para que varias instalaciones
en la fuente presentó menor tiempo para la transforma- de compostaje compartan infraestructura y equipos
ción de la materia orgánica (entre 10 y 19 días) y facilitó que les permitan el desarrollo de pruebas con cierta
la obtención de un producto con mayores contenidos complejidad.
de COT (9.7% respecto a 5.6% de SSP) y PTotal (0.65% res- De acuerdo con el uso del compost de biorresiduos,
pecto a 0.36% de SSP), mejores condiciones en cuanto a se ha identificado que mejora la fertilidad del suelo, es
CRA (173.1% respecto a 107.9% de SSP) y densidad fuente de nutrientes y materia orgánica, contribuye con
aparente (0.44 g/cm3 respecto a 0.65 g/cm3 de SSP), me- la estructura y el mantenimiento de las propiedades del
nor proporción de impurezas (3.47% respecto de 7.4% suelo (Odlare et al., 2011) y con la actividad biológica
de SSP) y mejor calidad microbiológica (menor en dos (Farrell y Jones, 2009). No obstante, los efectos depen-
órdenes de magnitud (1x102) respecto de SSP), factores den de las características del suelo, donde es necesario
clave de la calidad del producto final. considerar sus requerimientos para la aplicación del
Entre los contaminantes más comunes están los me- producto.
tales pesados; su presencia en el producto depende de En el experimento de Rosero (2013) (desarrollado en
su concentración en los biorresiduos (i.e. asociado a de- el marco de las investigaciones de los autores) se evaluó
ficiente separación en la fuente y recolección selectiva el crecimiento de maíz en un suelo andisol, con la apli-
de los residuos sólidos) y de la posible contaminación cación del producto mejorado provenientes del experi-
en el proceso; los metales pesados presentes con mayor mento con la adición de PE (tabla 3), con y sin
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM 37
Avances en investigación sobre el compostaje de biorresiduos en municipios menores de países en desarrollo . Lecciones desde Colombia
38 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM
Oviedo-Ocaña Edgar Ricardo, Marmolejo-Rebellon Luis Fernando, Torres-Lozada Patricia
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM 39
Avances en investigación sobre el compostaje de biorresiduos en municipios menores de países en desarrollo . Lecciones desde Colombia
40 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM
Oviedo-Ocaña Edgar Ricardo, Marmolejo-Rebellon Luis Fernando, Torres-Lozada Patricia
consulta: mayo 25, 2015]. Disponible en: [Link] International Conference, (3rd, 2009, Beijing, China), Bioinfor-
[Link]/urban/solid_wm/erm/CWG%20folder/[Link] matics and Biomedical Engineering, IEEE Editors, 2009, pp.
Hoornweg D., Bhada-Tata P. What a waste: A global review of 1-4.
solid waste management. Washington DC, Estados Unidos, Marmolejo L.F., Oviedo-Ocaña E.R., Jaimes J.C., Torres-Lozada P.
Urban Development Series Knowledge Papers No. 15. The Influencia de la separación en la fuente sobre el compostaje de
World Bank, 2012, 98 p. residuos sólidos municipales. Agronomía Colombiana, volu-
Huerta-Pujol O., Gallart M., Soliva M., Martínez-Farré F.X., López men 28 (número 2), 2010: 319-328.
M. Effect of collection system on mineral content of biowaste. Marmolejo L.F. Marco conceptual para la sostenibilidad de los sistemas
Resource Conservation and Recycling, volumen 55 (número 11), de aprovechamiento de residuos sólidos en cabeceras municipales
2011:1095-1099. menores a 20.000 habitantes del Valle del Cauca, (tesis doctorado
Iacovidou E., Ohandja D., Gronow J., Voulvoulis N. The House- en ingeniería), Colombia, Universidad del Valle, 2011, 190 p.
hold use of food waste disposal units as a waste management Odlare M., Arthurson V., Pell M., Svensson K., Nehrenheim E.,
option: A review. Critical reviews in environmental science and Abubaker J. Land application of organic waste – Effects on the
technology, volumen 42 (número 14), 2012:1485-1508. soil ecosystem. Applied Energy, volumen 88 (número 6), 2011:
Insam H., de Bertoldi M. Microbiology of Composting Process, 2210-2218.
en: Diaz L.F., de Bertoldi M., Bidlingmaier W., Stentiford E., OPS- Organización Panamericana de la Salud, BID- Banco Intera-
Waste Management Series 8, Compost Science and Technolo- mericano de Desarrollo, AIDIS- Asociación interamericana de
gy, United Kingdom, Elsevier Ltd., 2007, pp. 25-48. Ingeniería Sanitaria y Ambiental. Informe Regional del Pro-
Iqbal M.K., y Shafiq T., Ahmed K. Characterization of bulking yecto Evaluación Regional del Manejo de Residuos Sólidos
agents and its effects on physical properties of compost. Biore- Urbanos en América Latina y el Caribe, Estados Unidos, OPS,
source Technology, volumen 101 (número 15), 2010:1913-1919. 2010, 120 p.
Kato K. y Miura N. Effect of matured compost as a bulking and Oviedo-Ocaña E.R. Estrategias para la optimización del proceso y la
inoculating agent on the microbial community and maturity calidad del producto del compostaje de biorresiduos en municipios
of cattle manure compost. Bioresource Technology, volumen 99 menores de países en desarrollo, (tesis doctorado en ingeniería),
(número 9), 2008:3372-3380. Colombia, Universidad del Valle, 2015, 141 p.
Kim J., Park J., In B., Kim D., Namkoong W. Evaluation of pilot- Oviedo-Ocaña E.R., Daza-Torres M., Marmolejo-Rebellón L., Oso-
scale in-vessel composting for food waste treatment. Journal of rio-Ortiz A., Torres-Lozada P. Influencia de la incorporación
hazardous materials, volumen 154 (números 1-3), 2008: 272-277. de pasto estrella como material de soporte (Cynodon Plectos-
Krogmann U., Körner I., Diaz L.F. Composting: Technology, en: tachyus) en el compostaje de biorresiduos de origen munici-
Christensen T., Solid Waste Technology & Management, United pal. Ingeniería y Desarrollo, volumen 31 (número 2), 2013:
Kingdom, Blackwell Publishing Ltd., 2010, pp. 533-568. 251-271.
Kuba T., Tschöll A., Partl C., Meyer K., Insam H. Wood Ash Ad- Oviedo-Ocaña E.R., Marmolejo-Rebellon L.F., Torres-Lozada P.
mixture to Organic Wastes Improves Compost and its Perfor- Evaluation of the addition of wood ash to control the pH of
mance. Agriculture Ecosystems and Environment, volumen 127 substrates in municipal biowaste composting. Ingeniería, In-
(números 1-2), 2008: 43-49. vestigación y Tecnología, volumen 15 (número 3), 2014a:469-
Kumar M., Ou Y., Lin J. Co-composting of green and food waste at 478.
low C/N ratio. Waste Management, volumen 30 (número 4), Oviedo-Ocaña E.R., Marmolejo-Rebellon L.F., Torres-Lozada P.
2010: 602-609. Influencia de la frecuencia de volteo para el control de la hu-
Kurian J. Experiencias de tratamiento de residuos sólidos urbanos medad de los sustratos en el compostaje de biorresiduos de
en India. Residuos, volumen 97 (número 3), 2007:24-30. origen municipal. Revista internacional de contaminación am-
Kurola J., Mona A., Kontro M., Talves M., Romantschuk, M. Wood biental, volumen 30 (número 1), 2014b: 91-100.
Ash for Application in Municipal Biowaste Composting. Bio- Oviedo-Ocaña E.R., Torres-Lozada P., Marmolejo-Rebellon L.F.,
resourceTechnology, volumen 102 (número 8), 2011:5214-5220. Daza-Torres M., Torres-López W.A., Andrade M., Abonia-
Li Z., Lu H., Ren L., He L. Experimental and modeling approaches Gonzalez R. Evaluation of the variability of biowastes from
for food waste composting: A review. Chemosphere, volume 93 small towns in developing countries. Starting point for impro-
(número 7), 2013:1247-1257. ving their quality and composting process conditions. Waste
López M., Soliva M., Martínez-Farré F.X., Fernández M., Huerta- Management, 2015a, en evaluación.
Pujol O. Evaluation of MSW organic fraction for composting: Oviedo-Ocaña E.R., Marmolejo-Rebellón L.F, Torres-Lozada P,
Separate collection or mechanical sorting. Resources, Conserva- Daza M., Andrade M., Torres-López W.A. Abonia-Gonzalez
tion and Recycling, volumen 45 (número 4), 2010: 222-228. R. Effect of adding bulking materials over the composting
Lu Y. y Guo J. Effect of matured compost as a bulking agent and process of municipal solid biowastes. Chilean Journal of Agri-
inoculating agent on composting of municipal solid waste, en: culture Research, 2015b, en prensa.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM 41
Avances en investigación sobre el compostaje de biorresiduos en municipios menores de países en desarrollo . Lecciones desde Colombia
Oviedo-Ocaña E.R., Torres-Lozada P., Marmolejo-Rebellon L.F., Utilization in Horticultural Cropping Systems, 1ra ed., Estados
Hoyos L.V., Gonzales S. Barrena R., Komilis D., Sanchez A. Unidos, Boca Raton, Lewis Publishers, 2001, pp. 95-120.
Stability and maturity of biowaste composts derived by small Sundberg C., Navia T. Is there still a role for composting? Waste
municipalities: Correlation among physical, chemical and Management & Research, volumen 32 (número 6), 2014: 459-460.
biological indices. Waste Management, 2015c [en línea]. Dispo- Troschinetz A., Mihelcic J. Sustainable recycling of municipal so-
nible en: [Link] lid waste in developing countries. Waste Management, volu-
Rosero E.P. Evaluación de la aplicación de compost proveniente de bio- men 29 (número 2), 2009:915- 923.
rresiduos del municipio de Versalles y residuos de poda de jardín en Vergara S.E., Tchobanoglous G. Municipal solid waste and the en-
las propiedades químicas y biológicas de un suelo andisol y en la vironment: a global perspective. Annual Review of Environment
producción de maiz (zea mays), (tesis pregrado), Colombia, Uni- and Resources, volumen 37, 2012: 277-309.
versidad del Valle, 2013, 119 p. Wong J., Fung S., Selvam A. Coal fly ash and lime addition enhan-
Saha J.K., Panwar N., Singh M.V. An assessment of municipal so- ces the rate and efficiency of decomposition of food waste du-
lid waste compost quality produced in different cities of India ring composting. Bioresource Technology, volumen 100 (número
in the perspective of developing quality control indices. Waste 13), 2009:3324-3331.
Management, volumen 30 (número 2), 2010:192-201. Zurbrügg C., Drescher S., Patel A., Sharatchandra H. Decentrali-
Shekdar A.V. Sustainable solid waste management: An integrated sed composting of urban waste – an overview of community
approach for Asian countries. Waste Management, volumen 29 and private initiatives in Indian cities. Waste Management, vo-
(número 4), 2009:1438-1448. lumen 24 (número 7), 2004: 655-662.
Slater R.A., Frederickson J. Composting municipal waste in the Zurbrügg C., Gfrerer M., Ashadi H., Brenner W., Küper D. Deter-
UK: some lessons from Europe. Resources, Conservation and Re- minants of sustainability in solid waste management–The
cycling, volumen 32 (número 3), 2001: 359-374. Gianyar Waste Recovery Project in Indonesia. Waste Manage-
Smårs S., Gustafsson L., Beck-Friis B., Jönsson H. Improvement of ment, volumen 32 (número 11), 2012: 2126-2133.
the composting time for household waste during an initial
low pH phase by mesophilic temperature control. Bioresource
Technology, volumen 84 (número 3), 2002: 237-241.
Este artículo se cita:
SSPD- Superintendencia de Servicios Públicos Domiciliarios.
Diagnóstico sectorial. Plantas de aprovechamiento de resi- Citación estilo Chicago
duos sólidos, Colombia, SSPD, 2008, 45 p. Oviedo-Ocaña, Edgar Ricardo, Luis Fernando Marmolejo-Rebellon,
Soto J., Torres-Lozada P., Oviedo-Ocaña E.R., Marmolejo-Rebe- Patricia Torres-Lozada. Avances en investigación sobre el compostaje
llón L.F., Zambrano P. Influencia del pH y nutrientes sobre de biorresiduos en municipios menores de países en desarrollo. Lec-
la biodegradabilidad aerobia de biorresiduos de origen ciones desde Colombia. Ingeniería Investigación y Tecnología, XVIII,
01 (2017): 31-42.
municipal. Afinidad, volumen 72 (número 572), 2015, en
prensa. Citación estilo ISO 690
Stentiford E., De Bertoldi M. Composting process, en: Christensen Oviedo-Ocaña E.R., Marmolejo-Rebellon L.F., Torres-Lozada P. Avan-
T., Solid Waste Technology & Management, 1ra ed., United King- ces en investigación sobre el compostaje de biorresiduos en munici-
dom, Blackwell Publishing Ltd., 2010, pp 515-532. pios menores de países en desarrollo. Lecciones desde Colombia.
Sullivan D.M., Miller, R.O. Compost quality attributes, measure- Ingeniería Investigación y Tecnología, volumen XVIII (número 1), ene-
ments and variability, en: Stoffella P.J., Kahn B.A. Compost ro-marzo 2017: 31-42.
42 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 31-42 ISSN 1405-7743 FI-UNAM
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 43-53
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: septiembre de 2015, aceptado: julio de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Resumen
Este artículo presenta un nuevo modelo para la expansión de sistemas eléctricos de distribución con penetración de
generación distribuida, considerando un planeamiento multi-etapa coordinado. El problema se formuló como un
modelo de programación no lineal entero mixto y se solucionó empleando un algoritmo de búsqueda tabú. Para cada
etapa del horizonte de planeamiento se consideró la instalación de nuevos elementos (alimentadores, subestaciones
y generadores distribuidos), aumento de la capacidad de elementos existentes (alimentadores y subestaciones) y
costos operativos asociados a las pérdidas técnicas de energía en alimentadores. Los resultados obtenidos validan la
metodología propuesta al encontrar menores costos frente a dos escenarios de la literatura: planeamiento multi-
etapa no coordinado y coordinado (ambos sin generación distribuida).
Descriptores: algoritmo búsqueda tabú, generación distribuida, planeamiento multi-etapa no coordinado y coordi-
nado, sistemas eléctricos de distribución.
Abstract
This paper presents a new model for expansion of electric distribution system with penetration of distributed generation, consid-
ering coordinated multi-stage planning. The problem is formulated as a mixed integer non-linear programming model and is
solved using a tabu search algorithm. For each stage of the planning horizon are considered the installation of new elements
(feeders, substations and distributed generation), upgrading of existing elements (feeders and substations), and operating costs
associated to energy technical losses in feeders. The obtained results validate the proposed methodology because total costs are
lower when they are compared with two scenarios used in the literature: non-coordinated and coordinated multi-stage planning
(both without distributed generation).
Keywords: tabu search algorithm, distributed generation, non-coordinated and coordinated multi-stage planning, electricity
distribution system.
Nuevo modelo para la expansión de sistemas eléctricos de distribución con generación distribuida considerando un planeamiento ...
44 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM
Hincapié-Isaza Ricardo Alberto, Home-Ortiz Juan Manuel, Gallego-Rendón Ramón Alfonso
a) b)
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM 45
Nuevo modelo para la expansión de sistemas eléctricos de distribución con generación distribuida considerando un planeamiento ...
ij∈Ω ∑ c∈Ω (
∑ CFNTij ,c ,t σ NTij ,c ,t ) + ∑ ∑ (CF
ij∈Ω ET c∈ΩTC
ETij ,c ,t σ ET
ij ,c ,t
)+
NT TC
nE
∑ ∑ CFNSi ,s ,t σ NSi ,s ,t
i∈Ω NS s∈ΩTS
( ) j∈Ω ES s∈ΩTS
(
+ ∑ ∑ CFES j ,s ,t σ ES j ,s ,t + )
∑ β (t −1) DE (1)
t =1 ∑ ∑ CFNGi ,g ,t σ NGi ,g ,t
i∈Ω NG g∈ΩTG
( ) (
+ ∑ ∑ CFEG j ,g ,t σ EG j ,g ,t +
j∈Ω EG g∈ΩTG
)
( )
nN
β DE ∑ ∑ ∑ 3 PEt DN l DE I ij ,c ,t Rij ,c σ NTij ,c ,t + σ ETij ,c ,t
2
=l 1 ij∈ΩT c∈ΩTC
Note que los primeros seis términos de la ecuación (1) segunda sumatoria en estos seis términos se asocia a la
tiene una doble sumatoria, las cuales representan la lo- capacidad (tipo) de un tramo de red, subestación o GD
calización y el tipo de un elemento, respectivamente. (c, s o g).
Respecto a la primera sumatoria, los términos tres al A diferencia de los primeros seis términos, el sépti-
seis tienen solamente un subíndice, debido a que se mo término tiene tres sumatorias. La primera sumato-
asocia a la ubicación de subestaciones o generadores ria hace referencia al nivel de la curva de duración de
distribuidos (i o j). Por otro lado, los términos uno y dos carga, la segunda sumatoria se asocia a un tramo de red
tienen dos subíndices, debido a que se encuentran rela- i-j y la tercera sumatoria al calibre asociado al conduc-
cionados con los dos nodos de un tramo de red (ij). La tor de dicho tramo de red.
∀i ∈ Ω N
j∈Ωi c∈ΩTC
(
S Si ,l ,t = S Di ,l ,t − SGi ,l ,t + ∑ Vi ,l ,t ∑ σ NTij ,c ,t + σ ETij ,c ,t I ij ,c ,t *
) ∀l ∈ {1, 2,..., nN } (2)
∀t ∈ {1, 2,..., nE}
∀ij ∈ ΩT
I ij ,c ,t ≤ ∑
c∈ΩTC
I max
c (σ NTij ,c ,t + σ ETij ,c ,t ) ∀l ∈ {1, 2,..., nN } (3)
∀t ∈ {1, 2,..., nE}
∀i ∈ Ω S
(4)
S Si ,l ,t ≤ ∑ S (σ
s∈ΩTS
max
Ss NSi ,s ,t + σ ESi ,s ,t ) ∀l ∈ {1, 2,..., nN }
∀t ∈ {1, 2,..., nE}
∀i ∈ ΩG
SGi ,l ,t ≤ ∑
g∈ΩTG
SGmax
g
(
σ NGi ,g ,t + σ EGi ,g ,t ) ∀l ∈ {1, 2,..., nN } (5)
∀t ∈ {1, 2,..., nE}
∀i ∈ Ω N
Vi min
≤ Vi ,l ,t ≤ Vi max
∀l ∈ {1, 2,..., nN } (6)
∀t ∈ {1, 2,..., nE}
∑ (σ
c∈ΩTC
NTij ,c ,t )
+ σ ETij ,c ,t ≤ 1 ∀ij ∈ ΩT ∀t ∈ {1, 2,..., nE} (7)
∑ (σ
s∈ΩTS
NSi ,s ,t + σ ESi ,s ,t ≤ 1 ) ∀i ∈ Ω S ∀t ∈ {1, 2,..., nE} (8)
∑
g∈ΩTG
(σ NG i , g ,t + σ EGi ,g ,t ≤ 1 ) ∀i ∈ ΩG ∀t ∈ {1, 2,..., nE}
(9)
46 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM
Hincapié-Isaza Ricardo Alberto, Home-Ortiz Juan Manuel, Gallego-Rendón Ramón Alfonso
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM 47
Nuevo modelo para la expansión de sistemas eléctricos de distribución con generación distribuida considerando un planeamiento ...
Sin embargo, es importante tener en cuenta que al co- sis comparativo bajo tres escenarios de planeación. El
nectar un GD al sistema de distribución se pierde la radia- primer escenario considera un modelo de planeación
lidad de la red, al circular corrientes en ambos sentidos. multi-etapa no coordinado sin generación distribuida
De acuerdo con esto, es necesario utilizar un flujo de carga (pseudo-dinámico sin GD - caso 1). El segundo escena-
que se adapte a estas características operativas. En este rio emplea un planeamiento multi-etapa coordinado
artículo se emplea un flujo de carga pa-ra redes levemen- sin generación distribuida (dinámico sin GD - caso 2).
te enmalladas (Shirmohammadi y Hong, 1988). El tercer escenario resuelve el problema bajo un planea-
miento multi-etapa coordinado con generación distri-
buida (dinámico con GD - caso 3). Para los tres casos se
Integración del algoritmo bt a la solución
empleó el mismo horizonte de planeamiento.
del problema
Después de aplicar la metodología propuesta se ob-
La metodología propuesta en este artículo emplea una tuvieron los resultados presentados en la tabla 4. Las
estrategia maestro-esclavo. La estrategia maestro (pro- configuraciones finales de cada caso (después de la eta-
blema de inversión) se soluciona con el algoritmo BT. pa 3) se presentan en las figuras 5, 6 y 7. Los calibres
La estrategia esclavo (problema operativo) se resuelve seleccionados se observan entre paréntesis al lado de
con el flujo de carga. Cada configuración propuesta por cada tramo de red respectivo.
la técnica de solución (configuración del vecindario ge-
nerado) se evalúa por medio de la aplicación del flujo
de carga, el cual permite hallar el valor de la función
objetivo y validar la factibilidad. Una vez determinados
estos valores, el algoritmo BT continúa su proceso, has-
ta que el algoritmo alcanza el criterio de parada.
Pruebas y resultados
La metodología propuesta se aplica a un sistema de
prueba de 54 nodos (Miranda et al., 1994), el cual se mo-
difica para permitir la localización de generadores distri-
buidos (figura 3). En este sistema los cuadrados negros
corresponden a subestaciones existentes (nodos 51 y 52)
y los blancos a propuestas (nodos 53 y 54). Las líneas
continuas son tramos de red existentes y las líneas pun-
teadas son tramos de red propuestos. Los puntos negros Figura 3. Sistema de prueba de 54 nodos
son nodos de paso y de carga. Los lugares candidatos
para ubicar generadores distribuidos son los nodos 6, 18,
31, 36 y 45. Los datos de conductores, subestaciones y
generadores distribuidos se encuentran en las tablas 1, 2
y 3, respectivamente. La base de datos de la topología del
sistema se puede consultar con los autores.
El horizonte de planeamiento se dividió en tres eta-
pas de cinco años, en las cuales se consideró una tasa de
crecimiento anual de la demanda de 2%. Se considera Figura 4. Curva de duración de carga empleada
una curva de duración de carga de 100%, 70% y 40% del
valor nominal de la demanda, con una duración de 1000, De los tres casos implementados, el caso 1 es el más fácil
5760 y 2000 horas, respectivamente (figura 4). La tasa de de implementar. En este las etapas se resuelven de manera
descuento es de 10% y el precio de la energía $250 [kWh]. consecutiva, donde su complejidad matemática es de una
La metodología se implementó en Fortran 90 en un sola etapa, por lo tanto, el tiempo de cómputo es el más
computador con procesador Intel Core i7-3612QM bajo de los tres casos. Sin embargo, la calidad de la res-
CPU2.10GHz x 8, con 6GB de memoria RAM bajo siste- puesta es de peor calidad, respecto a los obtenidos en los
ma operativo UNIX. otros dos casos. En el caso 2 las etapas se estudian de ma-
Con el fin de verificar el comportamiento de la me- nera simultánea, incrementando ostensiblemente su com-
todología propuesta, en esta sección se realiza el análi- plejidad matemática y por ende, el tiempo de cómputo
48 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM
Hincapié-Isaza Ricardo Alberto, Home-Ortiz Juan Manuel, Gallego-Rendón Ramón Alfonso
requerido para su solución. A diferencia del caso anterior 52.51%, respectivamente. En el caso 1 se propone la ins-
(caso 1), se obtienen redes de menor costo, ya que requiere talación de las dos nuevas subestaciones (tipo 1), ya que
menos elementos al efectuar una adecuada distribución las dos subestaciones existentes no son capaces de ali-
de los mismos en el tiempo. En relación con el caso 3, ade- mentar el crecimiento de la demanda. Adicionalmente
más de considerar un planeamiento coordinado, conside- una subestación existente aumenta de capacidad al tipo
ra en su estudio la generación distribuida, incrementando 2. En el caso 2, también se construyen las dos nuevas
su complejidad matemática y con ello el tiempo de cóm- subestaciones; sin embargo, no aumentan capacidad las
puto. Por otro lado, la generación distribuida impacta po- otras subestaciones. En el caso 3 solamente se construye
sitivamente la distribución de flujos en la red, llevándolos una subestación tipo 1 en el nodo 53.
a valores que resultan favorables tanto en lo relacionado Respecto a los generadores distribuidos, en el caso 3
con las pérdidas de energía, como en la carga por los con- se instalan tres generadores de tipo 1 en los nodos 18, 31
ductores. Estos aspectos favorecen positivamente los cos- y 45, con un costo de M$975.
tos en inversión y operación, como se puede observar en Adicionalmente se observa que el caso 3 tiene el ma-
la tabla 4, en los casos 1, 2 y 3, respectivamente. yor costo de inversión (M$8043,5) en tramos de red. Los
De la tabla 4 se observa que el mayor costo de inver- casos 1 y 2 son menores un 5.09% y 2.75%, respectiva-
sión en subestaciones se presenta en el caso 1 mente. Esto se debe a que en el caso 3 se instalan conduc-
(M$8424.10), que es mayor a los casos 2 y 3 un 5.03% y tores de mayor capacidad, respecto a los otros dos casos.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM 49
Nuevo modelo para la expansión de sistemas eléctricos de distribución con generación distribuida considerando un planeamiento ...
Conclusiones
En este artículo se presenta una metodología para solu-
cionar el problema de la expansión de sistemas de dis-
tribución de energía eléctrica, usando como técnica de
solución un algoritmo de búsqueda tabú. El problema
se soluciona usando un planeamiento multi-etapa coor-
dinado, considerando la penetración de generación dis-
tribuida. Como análisis comparativo se realizan otros
Figura 6. Configuración final – Caso 2
dos casos de prueba, los cuales consideran modelos de
planeamiento multi-etapa no coordinado y coordinado
(ambos sin GD). Los resultados obtenidos muestran
que al emplear un modelo de planeamiento como el
propuesto en este artículo (multi-etapa coordinado con
GD), se reducen los costos totales de inversión al com-
pararse con los otros dos modelos de planeamiento
(planeamiento multi-etapa no coordinado y coordina-
do sin GD).
La metodología propuesta es útil para las empresas
distribuidoras de energía eléctrica debido a que la inver-
sión realizada en elementos de red (tramos de red, sub-
estaciones y generadores distribuidos), se compensa en
el tiempo con la disminución en el valor de las pérdidas
de energía, lo cual garantiza una mejor relación costo/
beneficio. Adicionalmente esta metodología se puede
emplear para sistemas con diferentes condiciones opera-
tivas, ya que permite modelar los GD como nodos PV o
Figura 7. Configuración final – Caso 3 PQ, dependiendo de la tecnología considerada.
50 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM
Hincapié-Isaza Ricardo Alberto, Home-Ortiz Juan Manuel, Gallego-Rendón Ramón Alfonso
Apéndice
En esta sección se presenta la nomenclatura empleada en este artículo.
Conjuntos
ΩEG Conjunto formado por GD existentes
ΩES Conjunto formado por subestaciones existentes
ΩET Conjunto formado por tramos de red existentes
ΩG Conjunto formado por GD nuevos y existentes
Ωi Conjunto formado por nodos conectados al nodo i
ΩN Conjunto de todos los nodos del sistema
ΩNG Conjunto formado por GD nuevos
ΩNS Conjunto formado por subestaciones nuevas
ΩNT Conjunto formado por tramos de red nuevos
ΩS Conjunto formado por subestaciones nuevas y existentes
ΩT Conjunto formado por tramos de red nuevos y existentes
Ω TC Conjunto formado por tipos de conductores
ΩTG Conjunto formado por tipos de GD
Ω TS Conjunto formado por tipos de subestaciones
Variables
σ EG i , g ,t
Variable binaria para repotenciar un GD existente en el nodo i, tipo g, en la etapa t
σ ESi ,s ,t Variable binaria para repotenciar una subestación existente en el nodo i, tipo s, en la etapa t
σ ETij ,c ,t Variable binaria para reconductorizar un tramo de red existente entre los nodos i-j, tipo c, en la etapa t
σ NG i , g ,t
Variable binaria para instalar un GD nuevo en el nodo i, tipo g, en la etapa t
σ NS i ,s ,t
Variable binaria para instalar una subestación nueva en el nodo i, tipo s, en la etapa t
σ NTij ,c ,t Variable binaria para instalar un tramo de red nuevo entre los nodos i-j, tipo c, en la etapa t
I ij ,c ,t Corriente entre los nodos i-j, en la etapa t (A)
SGi ,l ,t Potencia inyectada por un GD en el nivel l de la curva de carga, en el nodo i, en la etapa t (MVA)
S S i ,l ,t Potencia inyectada por una subestación en el nivel l de la curva de carga, en el nodo i, en la etapa t (MVA)
Vi ,l ,t Magnitud de la tensión en el nivel l de la curva de carga, en el nodo i, en la etapa t (kV)
Parámetros
β Tasa de descuento (%)
CFEGi ,g ,t Costo fijo de repotenciar un GD existente en el nodo i, tipo g, en la etapa t ($)
CFESi ,s ,t Costo fijo de repotenciar una subestación existente en el nodo i, tipo s, en la etapa t ($)
CFETij ,c ,t Costo fijo de reconductorizar un tramo de red existente entre los nodos i-j, tipo c, en la etapa t ($)
CFNGi ,g ,t Costo fijo de instalar un GD nuevo en el nodo i, tipo g, en la etapa t ($)
CFNSi ,s ,t Costo fijo de instalar una subestación nueva en el nodo i, tipo s, en la etapa t ($)
CFNTij ,c ,t Costo fijo de instalar un tramo de red nuevo entre los nodos i-j, tipo c, en la etapa t ($)
DE Duración de la etapa de planeamiento en años
DNl Duración del nivel l de la curva de duración de carga en un año (horas)
max
I c Capacidad máxima del conductor tipo c (A)
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM 51
Nuevo modelo para la expansión de sistemas eléctricos de distribución con generación distribuida considerando un planeamiento ...
max
SGg Capacidad máxima del generador distribuido tipo g (MVA)
max
S Ss Capacidad máxima de la subestación tipo s (MVA)
min
Vi Máxima tensión permitida en el nodo i (kV)
Vimax Mínima tensión permitida en el nodo i (kV)
Debido a las características de la metodología pro- Miranda V, Ranito J.V., Proença L.M. Genetic algorithms in opti-
puesta se puede aplicar en sistemas de distribución de mal multistage distribution network planning. IEEE Transac-
diferentes tamaños, características operativas, con re- tions on Power Systems, volumen 9 (número 4), 1994: 1927-1933.
des levemente enmalladas o inclusión de generación Nayanatara C., Baskaran J., Kothari D.P. Optimal location of dis-
distribuida, lo cual la hace atractiva para las empresas tributed generation using micro-genetic algorithm, en: 3rd
del sector eléctrico al permitir su adaptación a sistemas International Conference on Computation of Power, Energy,
de distribución bajo un ambiente de redes inteligentes Information and Communication (ICCPEIC), (Chenai, India),
(Smart Grids). 2014, pp. 525-530.
Paul Scherrer Institute (PSI). World Energy Scenarios: Composing
energy futures to 2050, World Energy Council (WEC), 2013
Referencias
[en línea] [fecha de consulta: 10 de agosto de 2015]. Disponible
Bagheri A., Monsef H., Lesani H. Integrated distribution net- en: [Link]
work expansion planning incorporating distributed genera- energy-scenarios-composing-energy-futures-to-2050/
tion considering uncertainties, reliability, and operational Popovic Z.N., Kerleta [Link]., Popovic D.S. Hybrid simulated an-
conditions. Electrical Power and Energy Systems, volumen 73, nealing and mixed integer linear programming algorithm for
2015: 56-70. optimal planning of radial distribution networks with distri-
El-Khattam W., Salama M.M.A. Distribution system planning buted generation. Electric Power Systems Research, volumen
using distributed generation, en: 16th Canadian Conference 108, 2014: 211-222.
on Electrical and Computer Engineering, (Montreal, Canadá), Ramírez I.J., Bernal J.L. Genetic algorithms applied to the design
2003, pp. 579-582. of large power distribution systems. IEEE Transactions on
Gallego R.A., Escobar A.H., Toro E. Técnicas Metaheurísticas de Op- Power Systems, volumen 13 (número 2), 1998: 696-703.
timización, segunda edición, Universidad Tecnológica de Pe- Rodrigues B., Sanches J.R., Cossi A.M., Contreras J. A multiobjec-
reira, Pereira, Colombia, 2008, pp. 167-216. tive model for distribution system planning based on tabu
Glover F. Y Laguna M. Tabu search, Boston, USA, Kluwer Acade- search, en: Proceedings 7th IEEE PowerTech, (Grenoble, Fran-
mic Publishers, 1997, pp. 59-91. cia), 2013, pp. 1-6.
Gopiya S.N., Khatod D.K., Sharma M.P. Analytical approach for Shaaban M.F., Atwa Y.M., El-Saadany E.F. A multi-objective ap-
optimal siting and sizing of distributed generation in radial proach for optimal DG allocation, en: 2nd International Con-
distribution networks. IET Generation. Transmission & Distri- ference on Electric Power and Energy Conversion Systems,
bution, volumen 9 (número 3), 2015: 209-220. (Sharjah, Emiratos), 2011, pp. 1-7.
Greatbanks J., Popovic D., Begovic M., Pregelj A., Green T.C. On Shirmohammadi D., Hong H.W. A Compensation-based power
optimization for security and reliability of power systems flow method for weakly meshed distribution and transmis-
with distributed generation, en: Proceedings 2nd IEEE Power- sion networks. IEEE Transactions on Power Systems, volumen 3
Tech, (Bologna, Italy), 2003, pp. 1-8. (número 2) 1988: 753-762.
Lavorato M., Rider M.J., García A.V., Romero R. A constructive Tabares A., Lavorato M., Franco J.F., Rider M.J. A mixed-binary
heuristic algorithm for distribution system planning. IEEE linear formulation for the distribution system expansion plan-
Transactions on Power Systems, volumen 25 (número 3), 2010: ning problem, en: 7th Transmission & Distribution Conferen-
1734-1742. ce and Exposition-Latin America, (Medellín, Colombia), 2014,
Mantway A.H., Al-Muhaini M.M. Multi-objective BPSO algorithm pp. 1-6.
for distribution system expansion planning including distri- Zou K., Agalgaonkar A.P., Muttaqi K.M., Perera S. Distribution
buted generation, en: IEEE Transmission and Distribution system planning with incorporating GD reactive capability
Conference and Exposition, (Chicago, USA), 2008, pp. 1-8. and system uncertainties. IEEE Transactions on Sustainable
Energy, volumen 3 (número 1), 2012: 112-123.
52 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM
Hincapié-Isaza Ricardo Alberto, Home-Ortiz Juan Manuel, Gallego-Rendón Ramón Alfonso
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 43-53 ISSN 1405-7743 FI-UNAM 53
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 55-65
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: enero de 2016, aceptado: agosto de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Sangabriel-Guillen Omar
Instituto Tecnológico y de Estudios Superiores de Monterrey, Campus Monterrey
Centro de Calidad y Manufactura
E-mail: osangabriel@[Link]
Abstract
The soft drinks bottling industry is one of the most successful industries in Mexico due to the high volume of con-
sumption of these drinks. However, changes on the market such as the introduction of new competitors and state
taxes have been producing a reduction on sales. In order to counter this effect, the sales and marketing functions of
these organizations must improve their results and to do that they need to improve their processes. Six Sigma is a
project driven methodology for process improvement that has been giving excellent results in the manufacturing
environment, but its use for sale and marketing processes is seldom. This case study will provide an example of how
it is possible to recognize Six Sigma projects by means of value driver maps for sales and marketing projects.
Keywords: Six Sigma, sales, marketing, project selection, value driver map, recognize.
Resumen
La industria embotelladora de bebidas carbonatadas es una de las industrias más exitosas en México debido al alto volumen de
consumo. Sin embargo, los cambios en el mercado, tales como la introducción de nuevos competidores y los impuestos federales
produjeron una reducción en las ventas. A fin de disminuir el efecto, las funciones de ventas y mercadotecnia de estas orga-
nizaciones deben mejorar sus resultados y hacer lo necesario para optimizar sus procesos. Six Sigma es una metodología para
mejorar procesos, se ejecuta mediante proyectos que han dado excelentes resultados en manufactura, pero su uso en procesos de
ventas y mercadotecnia es muy precario. Este caso de estudio proveerá un ejemplo de cómo es posible reconocer proyectos de Six
Sigma mediante Mapas de Direción de Valor para proyectos de mercadotecnia y ventas.
Descriptores: Six Sigma, ventas, mercadotecnia, selección de proyectos, mapa de dirección de valor, reconocer.
Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales and Marketing for Soft Drinks Bottling Industry
56 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen, De la Rosa-Hernández Romina Marisol
phases Define, Measure, Analyze, Improvement and a science. One starting point is to select projects that top
Control. Some variations were made to the DMAIC management would “buy”, however, there is no “ideal”
process in order to integrate the organizational goals procedure on how to select projects for SS, while there
(Recognize) and to assure the continuity of the impro- have been several attempts to make a process to select
vement achieved with projects (Replicate/Review), projects they fall short on how to do it constantly
transforming the process from DMAIC to R-DMAIC-R (Thompson, 2007).
(Temblador y Beruvides, 2009). This paper will present a framework for project se-
The structure of the DMAIC methodology allows lection in the sales and marketing environment and its
the use at any business function. Since the early days of results when applied to a soft drinks bottling company.
SS, there has been an unfortunate but common percep-
tion that SS is only oriented to manufacturing processes Proposed framework
and that a fact based problem solving methodology
does not transfer well to transactional processes, such Our proposed framework consists of the following sta-
as sales and marketing (Pestorius, 2007). ges:
The fundamental problem in these business
functions is due to inefficiency or ineffectiveness of at I. Establish commitment and goals with Chief Marke-
least one component of the sales process or lack of in- ting and Sales Officer (CMO) and his/hers staff.
formation within the process (Webb, 2006). II. Elaboration of a value driver map.
Pestorius (2007) mentions several reasons for which III. Project selection.
people working outside manufacturing may be uneasy
about using SS as a problem solving method, he lists Stage I. Establish commitment and with CMO
the following:
SS as methodology could be used to improve business
1) In manufacturing it is believed that almost every va- profitability, to drive out waste, to reduce cost of poor
riable can be controlled. quality and to improve the effectiveness and efficiency
2) Professional backgrounds of SS consultants are of all operations so as to meet or even exceed custo-
usually from manufacturing areas, that may repre- mers’ needs and expectations (Banuelas and Antony,
sent a barrier to the understanding of sales and mar- 2002). In the past, this methodology has been predomi-
keting processes. nantly used to improve manufacturing processes, but it
3) Quality is not an attribute typically measured in sa- is now increasingly applied to a wide variety of nonma-
les and marketing. nufacturing operations also known as transactional
4) Sales and marketing do not easily adopt standardi- processes.
zed processes. Most of the relevant transactional processes can be
found in the sales and marketing functions. Selling and
However, Guaspari (2001) mentions that it is possible marketing exist because someone must persuade so-
to implement SS on sales and marketing as long as the meone else to listen carefully, speak frankly, make a
tools used on the methodology are tailored to these decision and typically part with his or her money
processes. (Webb, 2006). These functions (S&M) must add value
In Table 1, the differences between manufacturing from a customer perspective.
and non-manufacturing on each of the DMAIC phases As any business improvement strategy, SS requires
are presented. commitment from top management also known as stee-
SS is a project driven methodology. An adequate ring committee. Typically, each member of the steering
project selection is one of the most important activities committee will take the role of sponsor of the program,
because it guides the identification and allocation of re- and therefore for a successful implementation they
sources in order to obtain more relevant benefits. One must be involved in the allocation of proper resources
addendum to the DMAIC structure was the introduc- for the program (Halliday, 2001). One major key for this
tion of the Recognize before Define phase. At the Re- to happen is to have an executive training of the
cognize phase, the customer, process and business methodology prior to the beginning of the initiative.
voice are identified and along with the strategic plan- Commitment often comes from the benefits (fore-
ning they are used to align the company’s objectives casted or previously obtained) of the implementation.
with the project list (Thompson, 2007). Pyzdek (2003) Antony (2006) mentions the following benefits from
says that project selection for SS is as much as an art as applying Six Sigma on transactional processes:
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM 57
Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales and Marketing for Soft Drinks Bottling Industry
• Effective decisions must be taken based on facts ins- • Discard systematically all the non-value operational
tead of emotions. Costs can be considerable reduced activities in order to improve delivery of service.
using well-known techniques for problem resolu-
tion. • The level of service will be reliable by minimizing
variability in service performance.
• Identify and appreciate the customers’ requirements
and expectations, it will influence the customer satis- • Improve preventive process instead of fire-fighting
faction and loyalty. approach.
• The shareholders satisfaction will increase through • Encourage teamwork in all levels of the company
trusted and capable internal operations.
Once the management is convinced from the potential
• Problems could be efficiently solved implementing benefits, it is time to establish how to tackle sales and
correct problem solving techniques and training. marketing problems. In order to address these pro-
‒ Accountability falls under a specific manager ‒ Accountability is frequently vague and therefore
and is clear and explicit to achieve the objectives of the quality’s goals are commonly unimportant
manufacturing
Measure ‒ Usually process maps or flowcharts exist before than a ‒ It is usual the absence of process maps or
SS project is begun flowcharts
‒ Measurements are typically repeatable and it is ‒ Measurements are typically not repeatable.
common in the measurement phase to perform a gauge Sometimes the attribute analysis is required and
repeatability and reproducibility analysis some others a measurement analysis is oriented more
in the reliability of the information rather than the
‒ External customers often specify limits for quality measurement process itself
characteristics
‒ Specification limits are defined by internal goals
rather than external customers
Analyze ‒ Many factors are divided into two categories: ‒ Uncontrollable or non-quantifiable factors are
controllable (temperature, speed, etc.) or quantifiable typical
noise factors (shift-to-shift or machine-to-machine
differences) that can be classified as blocks
Improve ‒ Fractional factorial experiments to determine the ‒ Experimentation is oriented to longitudinal studies
effects of factors on the quality characteristics are used rather than transversal. It is required a fine-tuning on
experimental design for transactional processes
‒ Improvements are made using direct changes on
controllable factors ‒ Improvements are typically performed by
mitigating the impact of uncontrollable factors on the
‒ A mathematical transfer function is often obtained process
‒ FMEA is used as a validation tool on the effectiveness ‒ FMEA has more value in transactional processes
and efficiency of the activities defined in measure/ than in manufacturing due to the activities defined
analysis phase for the highest scores of RPN are typically the inputs
for policy and procedure definition
58 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen, De la Rosa-Hernández Romina Marisol
blems, a value driver map is required to deploy, from the Value Driver Map and it should be place at the
the strategic down to the operative, performance indi- left of the sheet.
cators to be used that present an overview of S&M 2. A participative session with stakeholders is requi-
functions. red to define the components or elements that have
impact on the “Y” indicator. In order to do this, it
should be asked how does Y is calculated? or what
Stage II. Elaboration of a value driver map
elements are required to estimate Y? This process
The goal of stage II is to identify the operative indica- would result on having on the right of the indicator
tors (X’s) that have an impact on the Key Business Indi- another level for process performance indicators
cators (KBI) that represents the process performance that have a relationship with the indicators on the
main result (Y) through a value driver map (VDM). A left.
VDM is a tree diagram that shows a relationship among 3. In order to ensure all indicators required are on the
strategic and operational performance indicators. An value driver map, a validation process is necessary
equivalent of the VDM is the CTQ tree that is used to before continuing to the next level. This should be
translate needs considered vital by the customer into done through a necessity and sufficiency analysis.
product and services characteristics and links these Necessity means that the indicator is required in or-
characteristics to organizational processes (Six Sigma der to go to the next level and sufficiency would be
Academy, 2002). The difference between VDM and when it is not needed another indicator to go to the
CTQ tree is the focus on the processes and its result in next level. Once those two conditions are achieved it
the first one against the customer’s needs in the second is safe to assume that is possible to go to the next
one. In the context of project selection for SS it is possi- level.
ble to use a value driver map to find which strategic 4. Step 2 and 3 repeat themselves until the last practi-
indicators are the results (needs) of certain performance cal level is achieved (this is assured by step 3). The
indicators (drivers) and which are the key issues to re- last level indicators are the ones that will be selected
solve in order to improve that performance (require- for the SS projects.
ments) (Figure 1).
The process to elaborate a value driver map is the The deliverable for stage II will be a tree-like diagram
following: (Figure 2) that establishes the relationship of the results
(Y) with what it has been done on the processes (X).
1. At first, it is necessary to establish which is the main
indicator to be improved. Usually a Key Business Stage III. Project selection
Indicator is selected, this indicator will be the “Y” of
Project selection is a critical factor on any continuous im-
provement program, the correct identification of SS pro-
Maintenance
downtime jects leads to significant benefits. Managers and employees
X1 gain This can induce confidence in management as well as
Changeover time
in employees about the efficacy of the Six Sigma initiative,
Equipment
utilization X2 and as a result it will promote future investment and
efforts in the initiative (Kumar et al., 2009).
Crew availability. There is research about the existing relationship of
X3
project selection and the stalling of such projects. There
Line speed are two main reasons identified; 1) errors of problem
80% effective X4 definition traceable to the absence of certain objectively
utilization Rate/cycle time verifiable characteristics in the project itself, 2) errors in
Y Bottlenecks
judgment during the decision-making process (Bacda-
X5
yan, 2001).
Scrap The selection of six sigma projects is considered a
X6 multi-criteria decision where most of the information re-
First pass yield
levant to the problem is difficult and conflicting in natu-
Rework
X7 re. Selection criteria need to be prioritized so that those
that are most critical to the overall success of the business
Figure 1. Indicators Value Driver Map will have the most impact on the project selection.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM 59
Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales and Marketing for Soft Drinks Bottling Industry
Case study
NO Is it necessary?
In order to implement this framework, we worked with
a bottling company that was interested on doing pro-
YES
cess improvement projects for the sales and marketing
function, so we were hired as External coaches.
At stage I, we met with the Chief Executive Officer
(CEO) of the organization and its Sales and Marketing
Is it sufficient?
Director to establish the objectives of the implementa-
tion and who would be the team that would work on
this particular SS implementation. The organization is
YES
structured as the figure 3 shown below.
Top management decided that for the first exercise
Select Project to elaborate the value driver map a team formed by the
Marketing Manager (1), Revenue Management Mana-
ger (1) and the Regional Sales Managers (4, one for each
sales region) would be appropriate.
On the next meeting, the CMO and his staff were in-
End
formed about the initiative of implementing SS on the Sa-
les and Marketing functions, we discussed topics such as:
Figure [Link] chart for Stage II
60 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen, De la Rosa-Hernández Romina Marisol
Table 2. Critical success factors for SS implementation. (Adapted from Banuelas et al., 2006)
Top Connected
Customer Measurable and Learning and
Financial impact management to business
Critical criteria impact feasible growth
commitment strategy and
core competence
Harry and
Schroeder (2000)
Pande et al.
(2000)
Snee
(2001)
Goldstein
(2001)
Breyfogle et al.
(2001)
Pyzdek
(2000, 2003)
Lynch and
Soloy (2003)
Antony
(2004)
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM 61
Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales and Marketing for Soft Drinks Bottling Industry
Closing CEO
CMO
Emerging drinks
Volume/Total Volume by
Category
Light Products
Product mix Volume/Total Carbonated
drinks Portfolio Volume
Sales forecast
compliance(%)
%Coverage(%)
Asset productivity
(Sales Income minus Sales Expenses)
Major Sales indicator
Prospecting
Volume
Y
Man-hours
Maintenance
Fuel efficiency
Claims
Expense on Marketing /
Expenses Marketing expenses/unit
Incomes (%)
Absentism
Rotation
Extra time
62 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM
13
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen, De la Rosa-Hernández Romina Marisol
After the projects were selected, Regional Managers people believe that measurement on these areas is hard
would have to choose a start date for training on Six to achieve.
Sigma for each one the project leaders. The value driver map is a straightforward tool that
makes clear the relationship between processes and fi-
Conclusions nances. In this way the proposed framework helps top
management to have a notion of how an indicator has a
One of the hardest challenges that any continuous im- certain impact over the bottom-line and therefore it is
provement program faces it is the lack of support from easier for them to continue supporting these initiatives.
top management. Even today many organizations con- As the framework requires input from the main stake-
sider continuous improvement as a spend instead of an holders it ensures that their expertise is rested on the
investment, on this type of organization a program is relationships that are formed at each level of the map.
likely to fail on deliver results or simply fade away as The use of these framework on the case study ensu-
another management initiative. red that the top management main concern, the major
Six Sigma has proven in its many years of existence sales indicator, was always the goal of the SS initiative
that is a project driven methodology that works and de- and that all efforts were align to achieve such goal.
livers real results, however much of its success depends The outputs from the case study in this paper show
on management support (time, people, financial resou- that the value driver map approach provides an alter-
rces, etc.) and the adequate project selection. Project se- native way for selecting Six Sigma projects on sales and
lection is one of the most critical and challenging marketing functions. One advantage of the proposed
activities faced by Six Sigma companies. Most organi- framework is that it can be used on other key business
zations are able to identify a host of project opportuni- areas of the organization without making any changes.
ties, but the difficulty arises in sizing and packaging Furthermore, the proposed framework establishes a
those opportunities to create meaningful projects. guideline to quantify and select projects based on im-
Top management often manages in the language of pact on the bottom-line.
currency, while the processes measure its results on It is proposed as a further research to explore the
performance indicators. An effective communication possibility of creating mathematical models that could
between these two parts is crucial to achieve results establish the relationship of indicators at its correspon-
that are on the benefit of the organization. dent levels. This way an organization could predict
On the sales and marketing functions it is even more what would be the impact of increasing or decreasing a
challenging to translate these requirements as most performance indicator on the bottom-line.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM 63
Use of Value Driver Maps for Six Sigma Project Selection: A Case Study on Sales and Marketing for Soft Drinks Bottling Industry
References Webb M. Sales and marketing the Six Sigma way, Kaplan Publishing,
2006.
Antony-Jiju. Six sigma in the UK service organisations: results
from a pilot survey. Managerial Auditing Journal, volume 19 Bibliography
(issue 8), 2004: 1006-1012.
Antony-Jiju. Six sigma for service processes. Business Processes Banker-Rajiv D., Potter-Gordon, Srinivasan-Dhinu. Association of
Management Journal, volume 12 (issue 2), 2006: 234-248. Nonfinancial Performance Measures with the Financial Per-
Bacdayan P. Quality improvement teams that stall due to poor formance of a Lodging Chain. Cornell Hospitality Quarterly,
project selection: An exploration of contributing factors. Total volume 46 (issue 4), 2005: 394-412.
Quality Management, volume 12 (issue 5), 2001: 589-598. Breyfogle F.W., Cupello J. & Meadows B. Managing Six Sigma, 1st
Banuelas R., Tennant Ch., Tuersley I., Tang S. Selection of six sig- Edition, New York, USA, John Wiley & Sons, Inc., 2001, ISBN
ma projects in the UK. The TQM Magazine, volume18 (issue 5), 9780471396734
2006: 514-527. Dahlgaard-Jens J., & Dahlgaard-Park-Su Mi. Lean production, six
Banuelas Coronado R. & Antony-Jiju. Critical success factors for sigma quality, TQM and company culture. The TQM Magazi-
the successful implementation of six sigma projects in organi- ne, volume 18 (issue 3), 2006: 263-281.
zations. The TQM Magazine, volume 14 (issue 2), 2002: 92-99. Does-Ronald, Heuvel-Edwin, Mast-Jeroen, & Bisgaard-Soren.
CNN Expansión [on line]. México, tercero en consumo de refres- Comparing non-manufacturing with traditional applications
cos. Mexico City, April 8th 2013. [ref July 19th, 2014]. Availa- of Six Sigma. Quality Engineering, 15(1), 2002, 177-182.
ble on: [Link] Falconi V. TQC Control de Calidad Total (al estilo Japonés). Belo
al-lidera-el-consumo-mundial-de-refresco Horizonte, Fundação Christiano Ottoni, 1992, pp.222.
Espinosa E. Baja compra de refrescos en el país. Dinero en imagen, Goldstein M. Six sigma program success factors. Six Sigma Forum
May 10th, 2014 [ref: July 19, 2014]. Available on: [Link] Magazine, volumen 1 (número 1), 2001: 36-45.
[Link]/2014-05-10/36986 Harry M. & Schroeder R. Six Sigma: The breakthrough management
Guaspari J. Six Sigma on sales. Across the Board, volume 38 (issue strategy revolutionizing the world’s top corporations, 1st ed., New
5), 2001: 67-68. York, USA, Doubleday, 2000, ISBN 9788129117731
Halliday S. So what exactly is Six Sigma? Works Management, vo- Ittner C.D. & Larcker D.F. Innovations in performance measure-
lume 54 (issue 1), 2001: 15. ment: Trends and research implications. Journal of Manage-
Kubiak T.M. & Benbow D.W. The certified Six Sigma black belt hand- ment Accounting Research, volume 10, 1998: 205-238.
book, 2nd Edition, Wisconsin, Quality Press, 2009, ISBN Ittner C.D. & Larcker D.F. Coming up short on nonfinancial per-
9780873897327 formance measurement. Harvard Business Review, volume 81
Kumar M., Antony J., Chong-Byung R. Project selection and its (issue 11), 2003: 88-95
impact on the successful deployment of Six Sigma. Business Kennedy F.A. & Brewer P.C. The Lean Enterprise and Traditional
Process Management Journal, volume 15 (issue 5), 2009: 669-686. Accounting-Is the Honeymoon Over? Journal of Corporate Ac-
Pestorius M.J. Applying the science of Six Sigma to the art of sales and counting & Finance, volume 17 (issue 6), 2006: 63-74.
marketing, 1st Edition, Milwaukee, USA, ASQ Quality Press, Lynch D. & Soloy B. Improving the effectiveness of six sigma pro-
2006, ISBN 9780873896962 ject champions, ASQ’s Six Sigma Conference, 2003.
Pestorius M.S. Six Sigma outside of manufacturing: Apply Six Sig- Maynard-Harold B., & Zandin-Kjell B. Maynard’s industrial engi-
ma to sales and marketing. Quality Progress, 2007, 19-24. neering handbook, 5th Edition, New York, McGraw-Hill, 2001,
Pyzdek T. The Six Sigma project planner, 1st Edition, USA, McGraw- ISBN 9780070411029
Hill, 2003, ISBN 9780071411837 Pande P., Neuman R., Cavanach R. The six sigma way: How GE,
Six Sigma Academy. The Black Belt Memory Jogger: A Pocket Motorola, and other op companies are honing their performance,
Guide to Six Sigma Success, Salem, NH: GOAL/QPC, 2002. 1st ed., New York, USA, McGraw-Hill, 2000, ISBN
Temblador-Perez M.D. & Beruvides M.G. Systems approach to pro- 9780071358064
blem solving using Six Sigma. American Society For Enginee- Parmenter D. Key performance indicators: Developing, implementing
ring Management Conference Proceedings, 2009, pp. 14-17. and using winning KPI’s, 2nd Edition, Hoboken, New Jersey,
Thompson-Steven J. Improving the performance of Six Sigma; A case John Wiley & Sons, 2010, ISBN 978470545157
study of the Six Sigma process at Ford Motor Company, Uni- Pyzdek T. Selecting six sigma projects. Quality Digest, 2000, [on
veristy of Bedfordshire, 2007. line] [Ref July 19, 2014]. Available on: [Link]
Tohidi-Hamid. & KhedriLiraviasl-Kourosh. Six Sigma methodo- [Link]/sept00/html/[Link]
logy and its relationship with lean manufacturing system. Ad- Rother M. & Shook J. Learning to see: value stream mapping to add
vances in Environmental Biology, volume 6 (issue 2), 2012: value and eliminate muda, 1st edition, Salem, New Hampshire
895-906. USA, GOAL/QPC, 2002, ISBN 1576810437.
64 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM
Sangabriel-Guillen Omar, Temblador-Pérez María del Carmen, De la Rosa-Hernández Romina Marisol
Sangabriel O. Diseño y aplicación de una herramienta para el reconoci- Citation for this article:
miento de oportunidades de mejora para la implementación de exce-
Chicago style citation
lencia operacional, tesis (maestría en ciencias con especialidad
Sangabriel-Guillen, Omar, María del Carmen Temblador-Pérez, Ro-
en sistemas de calidad y productividad, MCP) - Tecnológico
mina Marisol de la Rosa-Hernández. Use of Value Driver Maps for Six
de Monterrey, Monterrey, México, 2012, 156 pp.
Sigma Project Selection: A Case Study on Sales and Marketing for Soft
Snee R. Dealing with the Achilles’ Heel of six sigma initiatives – Drinks Bottling Industry. Ingeniería Investigación y Tecnología, XVIII,
project selection. Quality Progress, volume 34 (issue 3), 2001: 01 (2017): 55-65.
66-69.
ISO 690 citation style
Sangabriel-Guillen O., Temblador-Pérez M.C., De la Rosa-Hernandez
R.M. Use of Value Driver Maps for Six Sigma Project Selection: A Case
Study on Sales and Marketing for Soft Drinks Bottling Industry. Inge-
niería Investigación y Tecnología, volume XVIII (issue 1), January-
March 2017: 55-65.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 55-65 ISSN 1405-7743 FI-UNAM 65
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 67-74
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: enero de 2016, aceptado: septiembre de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Ramírez-Méndez Esmeralda
Corporación Mexicana de Investigación en Materiales S.A. de C.V.
Correo: esmeraldaramirez358@[Link]
Cantú-Sifuentes Mario
Universidad Autónoma Agraria Antonio Narro
Departamento de Estadística y Cálculo
Correo: mcansif@[Link]
Resumen
El control estadístico de procesos es un conjunto de técnicas donde su empleo puede contribuir a mejorar la calidad
de cualquier proceso industrial. En particular, las cartas de control o monitoreo son herramientas poderosas para
detectar deterioro en la calidad de un proceso. Sin embargo, cuando las variables o características de calidad del
proceso son discretas, el uso de una carta de control puede implicar altos costos de muestreo. Aunque muchos es-
quemas existentes no restringen el tamaño de muestra a un valor determinado, se sabe que la elección del tamaño
de muestra tiene un efecto en el desempeño de una carta de control: usar un tamaño de muestra grande implica
mejor precisión en la estimación, por lo tanto, la carta tendrá mejor desempeño, y viceversa, usar un tamaño de
muestra pequeño implica peor precisión, por lo tanto, la carta tendrá un desempeño pobre. En este artículo se estudia
el desempeño de la carta de control D2 para variables discretas y se analiza el efecto que tiene el tamaño de muestra.
Adicionalmente, se obtienen diseños óptimos de la carta. La medida de desempeño que se emplea es la longitud
promedio de corrida (ARL, Average Run Length).
Descriptores: ARL, cartas de control, multinomial, tamaño de muestra.
Abstract
The application of the Statistical Process Control or SPC is useful to improve any industrial process. For example, the control
charts are powerful tools to detect quality deterioration. However, high sampling costs can be involved when the quality charac-
teristics are discrete. It is known that the sample size affects the control charts performance: a large sample size implies high
performance and vice versa, a small sample implies worst performance. In this article, the performance of D2 control chart is
studied and the sample size effect is analyzed. Additionally, optimal designs are obtained. The analysis is made in ARL terms.
Keywords: ARL, control charts, multinomial, sample size.
Carta de control para procesos multinomiales : estudio de desempeño
68 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM
Ramírez-Méndez Esmeralda, Cantú-Sifuentes Mario
La carta de control D2, propuesta por Mukhopadh- objetivo se fija para representar un proceso de alta ca-
yay (2008), no se basa en la distribución Chi cuadrada. lidad: proporciones bajas en las categorías no confor-
Esto implica que, en teoría, no se requieran tamaños mes y proporción alta en la categoría conforme.
grandes de muestra para funcionar de manera eficiente. Para el caso multinomial, una distancia de Mahala-
Para obtener una buena aproximación a la Chi cuadra- nobis generalizada se define como
da es necesario que el número esperado de piezas por
atributo sea de al menos 5. De otra manera, se hace una k N i ( pij − p j )2
Di2 = ( pi − p )T ∑ i−1 ( pi − p ) = ∑ (1)
corrección por finitud, por ejemplo la propuesta por Ya- j =1 pj
tes (1934). Sin embargo, resulta necesario comprobar
que la carta propuesta por Mukhopadhyay es eficiente Donde, Di2 mide, en un punto particular de tiempo, la
con muestras pequeñas y determinar si es una alterna- distancia entre el vector de proporciones observadas y
tiva adecuada para situaciones como las descritas en el el vector de proporciones objetivo. El límite de control
párrafo anterior. Con este fin, en este trabajo se estudia superior (LCS) es
el desempeño de la carta con diferentes tamaños de
muestra para investigar el efecto de reducir el tamaño N ( K − 1)
de muestra sobre la eficiencia de la carta. Adicional- LCS = i FK −1, Ni − K + 2 ,α (2)
mente se presentan diseños óptimos de carta para dife- Ni − K + 2
rentes escenarios. El desempeño se midió en términos
de la longitud promedio de corrida (ARL, Average Run La cantidad, FK‒1, Ni‒K+2, α es el cuantil α de la distribución
Length). Los resultados sugieren que la carta de control F con (K‒1) y Ni ‒K+2 y grados de libertad. El límite de
no puede ser eficiente cuando el tamaño de muestra es control inferior es cero. Ver Mukhopadhyay (2008) para
pequeño, por tal motivo se buscaron diseños óptimos más detalles.
que maximizan la eficiencia de la carta para esos casos. Una medida de desempeño comúnmente utillizada
Este artículo se organiza de la siguiente manera: la in- para evaluar la eficiencia de una carta de monitoreo es
troducción se presenta en la primera sección, la segunda el ARL. El ARL mide la rapidez con que la carta detecta
sección contiene el desarrollo de la investigación, siguien- un cambio en la calidad. Se supone que el proceso
do con el análisis y discusión de resultados en la tercera arranca dentro de control, es decir, p = p0 ; y en algún
sección. Finalmente, se presentan las conclusiones. tiempo después la calidad del proceso decrece ( p = p1 ) .
A continuación se presenta una breve descripción
Desarrollo del ARL.
La longitud promedio de corrida (ARL) se usa para
Enseguida se expone una descripción del gráfico D2. medir la rapidez con la que una carta de control detecta
Este gráfico se propuso por Mukhopadhyay (2008) y es un cambio o una perturbación del proceso que afecta la
capaz de monitorear más de un atributo simultánea- calidad del producto. El ARL se define como el número
mente. Se basa en una generalización de la distancia de promedio de muestras que se requieren para que la car-
Mahalanobis y en la distribución multinomial. Este grá- ta muestre una señal fuera de control. Si las observacio-
fico se adecua para monitorear procesos que puedan nes graficadas en la carta de control son independientes,
describirse de la siguiente manera: entonces el número de puntos graficados hasta que un
Durante un proceso de producción, un artículo ter- punto excede un límite de control es una variable alea-
minado se clasifica en una y solo una de k categorías toria geométrica como parámetro p. La media de la dis-
de defectos incluyendo la categoría de piezas confor- tribución geométrica es 1/p.
mes, es decir, hay k = 1 defectos en total. Sea pij la pro- La longitud promedio de corrida (ARL) se usa para
porción observada de artículos clasificados en la ca- medir la rapidez con la que una carta de control detec-
tegoría j, j = 1, 2, ..., k en una muestra de tamaño Ni, y ta un cambio o una perturbación del proceso que afec-
(pi1, pi2, pi3, ..., pik) = piT es el vector de proporciones aso- ta la calidad del producto. El ARL se define como el
ciadas a cada categoría. Ahora suponga que número promedio de muestras que se requieren para
p T = [ p1 , p2 ,..., pk ] denota el vector de categorías no con- que la carta muestre una señal fuera de control. Si las
formes en un proceso bajo control. p T puede estimarse observaciones graficadas en la carta de control son in-
mediante una base histórica de datos o especificarse dependientes, entonces el número de puntos grafica-
arbitrariamente. Entonces, piT tiene una distribución dos hasta que un punto excede un límite de control es
multinomial con parámetro p T. En lo consecutivo, p T se una variable aleatoria geométrica como parámetro p.
referirá como vector objetivo. Es claro que el vector La media de la distribución geométrica es 1/p.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM 69
Carta de control para procesos multinomiales : estudio de desempeño
Una carta de control es una secuencia de pruebas de Pr(E) = ∫ ∞FK−1, N −K+2 ,α fλ , K −1, Ni − K + 2 ( x) (5)
i
hipótesis. Entonces existen errores asociados conocidos
como error tipo I y error tipo II, que consisten en recha- Donde fλ , K −1, Ni − K + 2 (⋅) es la función de densidad de una
zar la hipótesis nula si es cierta y aceptar la hipótesis nula distribución F no central con parámetro de no centrali-
si es falsa, respectivamente. En este contexto, la hipótesis dad λ y con K − 1 y Ni − K + 2 grados de libertad. α, la
nula establece que el proceso se encuentra dentro de con- probabilidad de un error tipo I, es fijo y la cantidad
trol. Entonces, el error tipo I ocurre cuando el proceso se FK −1, N −K +2 ,α es el punto percentil (1 – α) de una distri-
i
declara fuera de control cuando no es así; a esto se le cono- bución F central con K − 1 y Ni − K + 2 grados de liber-
ce como una falsa alarma. En este caso, el ARL se denota tad. De hecho, α se fijó en 0.01. C= N (K − 1) . i
Ni − K + 2
por ARL0 y se desea que sea grande, de forma que se
tengan pocas falsas alarmas. El ARL0 se expresa en (3). Resultados y discusión
Por otra parte, cuando el proceso está fuera de control, el
ARL debe ser pequeño para proveer una rápida detec- Obtención de las curvas de potencia
ción de perturbaciones en la calidad del proceso. En este
caso, se denota por ARL1 y se calcula mediante (4). Para estudiar el desempeño de la carta, se obtuvieron
con (7) las curvas de potencia, suponiendo número de
1 categorías K = {3,6} y con tamaños de muestra
ARL0 = (3)
n, n = {10, 20, 30, 40, 50, 60, 70, 100}. Los resultados se
P
muestran en las figuras 1 y 2, respectivamente. Como el
Donde P es la probabilidad de cometer error tipo I; ARL está en función de la potencia, por razones de apa-
cuando esta probabilidad se especifica se denota con la riencia se optó por mostrar las curvas de potencia en
letra α. lugar de las de ARL.
Como puede observarse en las figuras 1 y 2, el tama-
1 ño de muestra tiene una influencia importante en el
ARL1 = (4)
P desempeño de la carta. Por ejemplo, la carta tiene el me-
jor desempeño cuando n = 100 y el peor cuando n = 10.
Donde, P es la probabilidad de concluir que el proceso
Aunque en Mukhopadhyay (2008) se sugiere que la car-
está fuera de control y que la hipótesis nula (el proceso
ta no debería presentar problemas relacionados al ta-
está dentro de control) es falsa. Se calcula como el com-
maño de muestra, con los resultados obtenidos, se
plemento de β, que denota la probabilidad de cometer
demuestra que la carta no puede ser eficiente cuando el
el error tipo II. Es fácil ver que el valor mínimo de ARL
tamaño de muestra es pequeño. Las curvas con el peor
es 1. Hay dos formas de calcular el ARL, una es obtener
desempeño para K = 3 son las que corresponden a
la expresión analítica de P y la otra es calcular un valor
n = {10, 20, 30} y para K = 6 son n = {10, 20, 30, 40, 50}.
aproximado del ARL mediante simulación. En este tra-
Además, puede observarse que el efecto es mucho ma-
bajo se obtuvo la forma analítica de P, considerando las
yor si el número de categorías K se acerca a n, el tamaño
propiedades distribucionales del estadístico de prueba.
de muestra. Así por ejemplo, la potencia de la carta
El evento “la carta señala un estado fuera de con-
para n = 10 es mucho menor para K = 6 que para K = 3,
trol” ocurre si= E { Di2 ≥ UCL}, entonces la probabilidad de
recordando que es deseable que la potencia, la probabi-
que la carta señale un estado fuera de control está dada
lidad de declarar que el proceso está fuera de control
por
sea grande, cuando λ > 0, que representa un proceso
P = Pr (E) multinomial fuera de control estadístico.
70 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM
Ramírez-Méndez Esmeralda, Cantú-Sifuentes Mario
proceso. Se usó un método de optimización, basado en En este modelo se necesita la solución x* = (xi|i = 1,
un algoritmo genético (AG), para encontrar los mejores 2, ..., r) que minimiza ARL1, así que el conjunto de pará-
parámetros de diseño x* = [N, LCS], el valor del tamaño metros x*∈χ debe determinarse. El algoritmo genético
de la muestra N, del límite de control LCS que minimi- (GA, Genetic Algorithm) usa una representación binaria
cen el ARL1 con ARL0, s y nmax fijos. Formalmente de una posible solución. Esta codificación se necesita,
ya que el GA manipula bits. En el modelo se empleó
1
Minimizar ARL1= una representación binaria de dieciséis por cada solu-
P
ción. La matriz de población inicial fue de 100 filas, que
donde P se calcula con la expresión (4), con λ > 0. representan cada bit de la solución binaria codificada,
es decir, cada individuo en la población representa un
Sea X, el espacio de búsqueda, donde x*∈χ. Las res- conjunto de (r × r) soluciones codificadas de 16 bits. Se
tricciones son: empleó selección por torneo de tamaño 2 con cruce de
un solo punto y una mutación simple. Las probabilida-
K ≤ N +1 (6) des de cruce y mutación fueron de 0.9 y 0.1, respectiva-
mente. El ARL1 definido en la ecuación (4) se usó como
K ≥ 2 (7) una función de evaluación con 200 generaciones en el
GA.
N ≤ nmax (8) Los parámetros de diseño que se muestran en la ta-
bla 1, se obtuvieron suponiendo ARL0min = 200, numero
ARL0 ≥ ARL0min (9) de categorías C = {3, 6, 9}, nmax = {20, 30}, y se considera-
ron magnitudes de corrimiento s = {0.5, 1.0, 1.5, 2.0, 2.5,
LCS ≥ 0 (10) 3.0}. Para cada combinación de valores, se encontró un
diseño óptimo mediante un GA (optimización local, ver
Donde (6) resulta de (ni – C + 2), (7) de (C – 1) los grados tabla 1). Por ejemplo, se obtuvo un diseño óptimo para
de libertad de la distribución F en (5), esas cantidades K = 3 y s= 0.5. Lo mismo se hizo para K = 6 y s = 1, así
deben ser iguales o mayores que 1, entonces sucesivamente.
K ≤ N + 1. La expresión (8) surge de la idea que el tama- El valor ARL0min= 200 se eligió porque se consideró
ño de muestra debe limitarse para acotar la búsqueda a una búsqueda de diseños con tamaños de muestra pe-
regiones de interés, de acuerdo con restricciones prác- queños. Seleccionar un valor grande (el valor estándar
ticas. La ecuación (9) es una restricción técnica. Por otra es 370) puede llevar a tamaños de muestra muy gran-
parte, el estadístico de prueba es una distancia, por lo des para asegurar una potencia de detección razonable
tanto, no es negativa, así el límite de control no debe ser (De Araújo et al., 2011).
negativo, lo cual se especifica en (10). Los valores K, Los diseños de la tabla 1 se pueden utilizar para mo-
nmax, ARL0min y s son entradas para el problema, en la nitorear procesos con características similares a los es-
práctica esta información se obtiene del proceso. cenarios probados, únicamente se elige la magnitud de
1 1
0.9 0.9
0.8 0.8
0.7 0.7
n=(10)
0.6 0.6 n=(10)
n=(20)
n=(20)
Potencia
Potencia
n=(30)
0.5 0.5 n=(30)
n=(40)
n=(40)
0.4 n=(50) 0.4 n=(50)
n=(60)
n=(60)
0.3 n=(70) 0.3 n=(70)
n=(80)
n=(80)
0.2 n=(90) 0.2 n=(90)
n=(100) n=(100)
0.1 0.1
0 0
0 10 20 30 40 50 60 70 0 10 20 30 40 50 60 70
Lambda Lambda
Figura 1. Curvas de potencia de la carta D2 para K = 3 Figura 2. Curvas de potencia de la carta D2 para K = 6
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM 71
Carta de control para procesos multinomiales : estudio de desempeño
72 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM
Ramírez-Méndez Esmeralda, Cantú-Sifuentes Mario
muestra, lo que evitará que la potencia de la carta tenga Jolayemi J.K. A Statistical model for the design of multiattribute
fluctuaciones. control charts. The indian journal of statistics, Series B, volumen
Cuando se emplea una carta de carta de control 61 (número 2),1999: 351-365.
multivariable o multiatributo, y la carta da una señal Khoo M.B.C., How-Chinh L., Zhang W., Chung-Ho C., Castaglio-
fuera de control, esta representa la información conte- la P. A synthetic double sampling control chart for the process
nida en las variables o atributos que se monitorean. Sin mean. IIE Transactions, volumen 43 (número 1), 2010: 23-38.
embargo, es común que no sea posible establecer, de doi:10.1080/0740817X.2010.491503.
manera específica, cuál o cuáles variables son las que Li J., Tsung F., Zou C. Multivariate binomial/multinomial con-
provocan que el proceso salga de control. De tal manera trol chart. IIE Transactions, volumen 46 (número 5), 2014:
que resulta necesario un método para identificar el atri- 526-524.
buto que con mayor probabilidad provoque una situa- Lu X.S., Xie M., Goh T.N., Lai C.D. Control chart for multivariate
ción fuera de control. Este método tendría un impacto attribute processes. International Journal of Production Research,
práctico, debido a que brindaría información útil para volumen 36, 1998: 3477-3489.
encontrar la causa raíz del problema. Investigar un mé- Mahadik S.B., y Shirke D.T. A Special variable sample size and
todo como el descrito aplicable a la carta de control es- sampling interval hotelling’s T 2 chart. The International Jour-
tudiada se deja como trabajo futuro. nal of Advanced Manufacturing Technology, volume 53 (números
1-4): 2011: 379-384. doi:10.1007/s00170-010-2819-8.
Mason R.L. y Joung J.C. Multivariate statistical process control with
Referencias
industrial applications, Society for industrial and applied
Bersimis S., Psarakis S., Panaretos J. Control charts : An overview. mathematics american statistical association, Philadelphia,
Quality and Reliability Engineering International, volumen 23 PA., 2002, pp. 4-5.
(noviembre 2006), 2007; 517-43, doi:10.1002/qre. Mukhopadhyay A.R. Multivariate attribute control chart using
Costa A., Branco F., Aparecida M., Machado G. Variable parame- mahalanobis D2 statistic. Journal of Applied Statistics, volumen
ter and double sampling charts in the presence of correlation: 35, 2008: 421-429. doi:10.1080/02664760701834980.
the markov chain approach. International Journal of Production Niaki-Seyed T.A. y Nasaji S.A. A hybrid method of artificial neu-
Economics, volumen 130 (número 2), Elsevier, 2011: 224-229. ral networks and simulated annealing in monitoring auto-
doi:10.1016/[Link].2010.12.021. correlated multi-attribute processes. The International Journal
De Araújo-Rodrigues, Aparecida A., Kahn-Epprecht E., Sacra- of Advanced Manufacturing Technology, volumen 56 (números
mento De Magalhães M. Double-sampling control charts for 5-8), 2011: 777-788. doi:10.1007/s00170-011-3199-4.
attributes. Journal of Applied Statistics, volumen 38 (número 1), Seif A., Moghadam M.B. Faraz A., Heuchenne C. Statistical me-
2011: 87-112. doi:10.1080/02664760903266007. rits and economic evaluation of t2 control charts with the
Faraz-Alireza R., Kazemzadeh B., Ahmad P., Moghadam M.B. On VSSC scheme. Arabian Journal for Science and Engineering, vo-
the advantages of economically designed the hotelling’s T2 lumen 36 (número 7), 2011: 1461-1470. doi:10.1007/s13369-
control chart with variable sample sizes and sampling inter- 011-0124.
vals. Quality & Quantity, volumen 46 (número 1), 2012: 39-53. Taleb H. y Limam M. Fuzzy Multinomial Control Charts. Bandini
doi:10.1007/s11135-010-9325-9. and S. Manzoni (Eds.): AI*IA 2005, LNAI 3673, 2005, pp. 553-
Ghute V.B. y Shirke D.T. A multivariate synthetic control chart for 563.
process dispersion. Quality Technology & Quantitative Manage- Taleb H. y Limam M. Multivariate fuzzy multinomial control
ment, volumen 5 (número 3), 2008: 271-288. charts. Quality Technology & Quantitative Management, volu-
Irianto D. y Juliani A. A two control limits double sampling con- men 3 (número 4): 2006: 437-453.
trol chart by optimizing producer and customer risks. ITB Yates F. Contingency Tables Involving Small Numbers and the χ2
Journal of Engineering Science, volumen 42 (número 2): 2010: Test. Supplement to the Journal of the Royal Statistical Society, vo-
165-178. lumen 1 (número 2), 1934: 217-235.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM 73
Carta de control para procesos multinomiales : estudio de desempeño
74 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 67-74 ISSN 1405-7743 FI-UNAM
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 75-86
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: febrero de 2016, aceptado: septiembre 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Felipe-Altamirano Salvador
National Polytechnic Institute
School of Mechanical and Electrical Engineering
Department of Graduate Studies and Research, Unit Culhuacan
Email: [Link]@[Link]
Loboda Igor
National Polytechnic Institute
School of Mechanical and Electrical Engineering
Department of Graduate Studies and Research, Unit Culhuacan
Email: iloboda@[Link]
Abstract
Since gas turbines are very complex and potentially unreliable machines, the improvement of their monitoring sys-
tems becomes an essential part. Considering this necessity, the present paper performs a gas turbine diagnostic al-
gorithm testing. The methodology proposed is formed by three stages. In the first stage, the commercial software
ProDiMES (Propulsion Diagnostic Method Evaluation Strategy) is used to simulate an engine fleet and generate data
with fault and no-fault conditions. In the second stage, a baseline model testing is implemented to improve the
healthy engine performance approximation. Finally, a fault recognition stage based on a pattern recognition techni-
que (Multi-Layer Perceptron) performs the diagnosis and calculates the probability of correct diagnostic decisions.
The results obtained show that: a) the software ProDiMES is an easy and convenient tool to evaluate gas turbine
diagnostic methods, b) the baseline model testing is a key step because it allows reducing the errors that can negati-
vely influence the diagnostic process and c) the algorithm correctly performs the fault recognition task.
Keywords: gas turbine diagnostics, baseline model, fault recognition, ProDiMES, multi-layer perceptron.
Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
Resumen
Debido a que las turbinas de gas son máquinas muy complejas y potencialmente no fiables, el mejoramiento de sus sistemas de
monitoreo se vuelve una parte esencial. Considerando esta necesidad, el presente trabajo realiza una prueba de algoritmo de di-
agnóstico de turbinas de gas. La metodología propuesta se forma por tres etapas. En la primera, el Software comercial ProDiMES
(Propulsion Diagnostic Method Evaluation Strategy) se emplea para simular una flota de motores y generar datos con condicio-
nes de falla y sin falla. En la segunda etapa, se implementa una prueba de modelos de referencia para mejorar la aproximación de
rendimiento de motor sano. Finalmente, una etapa de reconocimiento de fallas basada en una técnica de reconocimiento de pa-
trones (Perceptrón Multicapa) que diagnostica y calcula la probabilidad de decisiones diagnósticas correctas. Los resultados ob-
tenidos muestran que: a) el software ProDiMES es una herramienta fácil y conveniente para evaluar métodos de diagnóstico de
turbinas de gas, b) la prueba de modelos de referencia es un paso clave porque permite reducir los errores que influyen negativa-
mente al proceso de diagnóstico y c) el algoritmo realiza correctamente la tarea de reconocimiento de fallas.
Descriptores: diagnóstico de turbinas de gas, modelo de referencia, reconocimiento de fallas, ProDiMES, perceptrón multicapa.
76 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de J esús , Pérez-Ruiz Juan Luis
based on pattern recognition techniques such as Radial MacIsaac, 1983). This work uses the software ProDi-
Basis Network (Loboda et al., 2010), Probabilistic Neural MES as a tool to simulate a fleet of engines and study
Network (Tsalavoutas et al., 2000), Bayesian Network the algorithms for gas turbine fault diagnostics. ProDi-
(Romessis and Mathioudakis, 2004) and Multi-Layer Per- MES was developed by the NASA Glenn Research Cen-
ceptron (MLP) (Roemer and Kacprzynski, 2000; Volponi ter originally intended for benchmarking international
et al., 2000; Sampath and Singh, 2004). The latter techni- gas path diagnostic methods (Simon et al., 2008, 2013).
que has been applied widely in the past years and has ProDiMES includes an Engine Fleet Simulation
shown that is it not inferior to other methods (Loboda et (EFS) that generates simulated measurement parame-
al., 2010). Therefore, this work uses the MLP to calculate ter histories for each engine of the fleet. To simulate real
the probability of correct diagnosis. The methodology engine behavior, it works with a deterioration profile,
utilized is implemented in Matlab. Its neural network noise level and operation mode (takeoff and cruise)
toolbox (Beale et al., 2014) assists in an efficient develop- unique for each engine. The EFS is implemented in
ment of the algorithm. Matlab and consists of a Graphical User Interface (GUI),
The paper is organized as follows. Section 2 gives a Case Generator and a C-MAPPS Steady State engine
the description of the methodology proposed. Section 3 model (Simon, 2010). Each element is described below.
describes the pattern recognition method used. Section
4 presents the results obtained for the baseline model
Graphical User Interface (GUI)
testing and the gas turbine fault recognition.
In the Graphical User Interface (GUI), the user can con-
Methodology proposed for gas turbine diagnostics trol the type and the number of faults occurring in the
engine fleet. The maximal number of possible simula-
The methodology proposed (Figure 1) includes the fo- ted faults is 18 plus a no-fault case. Table 1 shows these
llowing main stages: fault cases and their magnitudes, within which the
faults are distributed uniformly. Module faults (ID 1-5)
1) Gas turbine data simulation using the software Pro- corresponding to Fan, LPC, HPC, HPT and LPT are si-
DiMES (Propulsion Diagnostic Method Evaluation mulated by adjusting at the same time efficiency and
Strategy), flow capacity parameters. Actuator faults (ID 6-7) co-
2) Baseline model testing to enhance the healthy engi- rresponding to VSV and VBV result from a mis-schedu-
ne performance approximation and ling between the commanded and current actuator
3) Gas turbine fault recognition and the calculation of position. There are also 11 different sensor faults (ID
diagnosis accuracy. The stages are described in the 8-18) whose magnitudes are in units of average measu-
below subsections. rement noise standard deviation σ. It is important to
mention that in EFS, each individual engine only expe-
Gas turbine data simulation using the software riences a single fault type. Other aspects to consider
ProDiMES are: the number of flights to generate the output data
(the maximal number is 5000 flights per engine); the
Most of the real gas turbine faults are not so severe or fault evolution type (abrupt or rapid); the flight of fault
rarely occur in practice to be sufficient for a complete initiation (the 11th flight is the lowest value) and the sen-
fault description. Also, physical experimentation is not sor noise (on or off).
a feasible option because of the very high costs. Instead, The EFS works with eleven sensed variables. The se-
mathematical models are applied (Saravanamuttoo and ven measured variables shown in Table 2 are available
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM 77
Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
for monitoring and are represented by a vector Y. Table Table 2. Monitored variables (Simon, 2010)
3 shows four measured variables used as operating con- ID Variable Symbol
ditions represented by a vector U . All the variables co- 1 Physical core speed Nc
rrespond to an engine standard measurement system. 2 Total pressure at LPC outlet P24
3 Static pressure at HPC outlet Ps30
4 Total temperature at LPC outlet T24
Case generator
5 Total temperature at HPC outlet T30
The case generator produces parameter histories after 6 Total temperature at HPT outlet T48
the user has selected the number and type of faults in 7 Fuel flow Wf
the GUI stage. One important characteristic is the ran-
Table 3. Operating conditions (Simon, 2010)
dom generation of unique faults, degradation profiles
ID Variable Symbol
and operating history for each engine in the fleet. This
1 Physical fan speed Nf
also includes the following assignations: the date
2 Total pressure at fan inlet P2
when the collection of the engine data starts, the city
3 Total temperature at fan inlet T2
pairs for the takeoffs of the engine, the ambient pres-
4 Ambient pressure Pamb
sure, the atmospheric temperature, the Mach number
and power setting parameters for takeoff and cruise.
The level and rate of gradual performance deteriora- C-MAPSS Steady-State engine model
tion for each engine are also considered by the Case
The Commercial Modular Aero-Propulsion System Simula-
Generator. They emulate the degradation that an air-
tion (C-MAPSS) Steady State is a high-bypass turbofan
craft engine experiences during its lifetime due to di-
engine model created for diagnostics research (Simon,
fferent effects such as fouling, erosion, and corrosion
2010). This model is run inside the EFS and receives the
of blades and vanes. The gradual deterioration is not
outputs from the Case Generator to produce the simu-
considered a fault and its development is much slower
lated measurement parameter histories for each engine,
than the produced by faults.
at takeoff and cruise of each flight. C-MAPSS works
Table 1. Simulated fault cases (Simon, 2010) with two spool speeds (fan and core speed). Figure 2
shows the station numbers, the modules and the simu-
ID Fault description Fault magnitude lated sensor variables of the C-MAPSS Steady-State
0 No-fault --- model.
1 Fan fault 1 a 7%
2 LPC fault 1 a 7% Baseline model testing
3 HPC fault 1 a 7%
4 HPT fault 1 a 7%
In order to know the current gas turbine condition by
means of measured gas path variables, it is necessary to
5 LPT fault 1 a 7%
describe correctly its healthy state. According to Loboda
6 VSV fault 1 a 7%
et al. (2004), a good approximation of healthy engine per-
7 VBV fault 1 a 19%
formance, also called baseline model, can be given by
8 Nf sensor fault ± 1 a 10 σ complete second order polynomials. Also, polynomials
9 Nc sensor fault ± 1 a 10 σ have shown to be better than other techniques (Loboda
10 P24 sensor fault ± 1 a 10 σ and Feldshteyn, 2010). Considering one monitored gas
11 Ps30 sensor fault ± 1 a 10 σ path variable as function of four operating condition ar-
12 T24 sensor fault ± 1 a 10 σ guments , the baseline model can be expressed as
13 T30 sensor fault ± 1 a 10 σ
14 T48 sensor fault ± 1 a 10 σ Y0 (U) =a1 + a2u1 + a3u2 + a4 u3 + a5u4 + a6 u1u2 + a7 u1u3 + a8 u1u4 + a0u2 u3 +
15 Wf sensor fault ± 1 a 10 σ
a8u1u4 + a9u2u3 + a10u2u4 + a11u3u4 + a12u12 + (1)
16 P2 sensor fault ± 1 a 10 σ
2 2 2
17 T2 sensor fault ± 1 a 10 σ a u +a u +a u
13 2 14 3 15 4
18 Pamb sensor fault ± 1 a 19 σ
*LPC=Low Pressure Compressor, HPC=High Pressure Compressor,
where a1, … , a15 are the model coefficients calculated
HPT=High Pressure Turbine, LPT=Low Pressure Turbine, VSV= Variable using the least squares method for all monitored
Stator Vane, VBV=Variable Bleed Valve variables.
78 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de J esús , Pérez-Ruiz Juan Luis
(1)
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM 79
Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
Gas turbine fault recognition dual model vector for all monitored variables is given
by
Deviation computation
Yr= Y φ(1+ δ Y φ) (7)
Due to the variation of gas turbine operation condi-
tions, absolute gas path monitored parameters change where δ Y φ and Yφ are vectors of correction deviations
as well. Since these changes are greater than the produ- and general model values respectively.
ced by faults, the latter remain hidden. Therefore, a
diagnostic process requires an important step of devia- 3. Final deviation calculation using individual models.
tion computation to reveal deterioration and faults Using individual baseline values Yr for a monitored
effects (Loboda et al., 2004). Three steps are proposed to variable, we obtain
calculate deviations needed for the gas turbine diag-
nostic algorithm: Y * −Yr
δ Yr = (8)
Yr
1) Initial deviation calculation using a general model,
2) Creation of individual models, where δYr is a final deviation. These deviations are the
3) Final deviation calculation using individual models. base of fault class formation.
Using the mean of these deviation vectors for the first D1, D2, ..., Dq (10)
n=10 simulated flights per engine (before the fault ini-
tiation) and considering one engine, we have Thus,
δ Yφ =
∑ nj=1 DYj (6) ∑ q
n −1
P( Dn ) = 1 and Pn≠ l ( Dn / Dl ) = 0 (11)
n
In ProDiMES, each fault class is constructed from pat-
where δ Y φ is a correction deviation vector and j is the terns with the change of only one fault parameter (sin-
index for flights. gular fault class).
80 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de J esús , Pérez-Ruiz Juan Luis
is formed by simulating all the 19 fault cases available ficial neural network intended for classification pro-
in ProDiMES for better verification of the diagnostic al- blems. It uses a back-propagation algorithm that pro-
gorithm (18 faults + 1 no-fault case), a determined num- pagates a signal for a given input vector, producing an
ber of engines per class and flights per engine for cruise output and adapting unknown coefficients based on
operating point. The flight fault initiation selected is 11 the error between a target and the network output. Fi-
with the option “fixed”, this means that the first 10 gure 6 shows the general structure of the MLP. The in-
flights will not experience any fault. The fault evolution put for each hidden
layer neuron is given by the sum of
rate is selected as “rapid”. Since accuracy of fault clas- an input vector Z multiplied by weights in a matrix W1
ses’ description depends on the number of simulated and a bias value (from a vector b1 ). The result is passed
patterns, 100 engines per class and 50 total flights per through a differentiable transfer function f1 (tansig)
engine are considered. However, the number of flights producing the neuron output (contained in a vector a1)
per engine for rapid fault evolution is 34. For this rea- within an interval of [-1, 1]. For the output layer neu-
son, 3400 patterns per class are employed. Thus, the to- rons, the computation is repeated in the same way
tal size of the training set is 64600 patterns (19 fault using the vector
a1 as input to the layer, a matrix W2, a
cases × 100 engines × 34 flights). bias vector b2, and a transfer function f2 (logsig), produ-
The validation set ZV is created to verify that the net- cing the network output contained in a vector a2 within
work can generalize the fault classes correctly. It is for- an interval of [0, 1]. The MLP is trained on known pair’s
med in the same manner as the training set ZT; however, vectors: the input vector and the target vector formed
its size is ten times smaller because it works with 340 by 1 and 0 representing the membership of a class. The
patterns per class. Therefore, the total size of the valida- number of input layer neurons is seven, which repre-
tion set is 6460 patterns (19 fault cases × 10 engines × 34 sent the seven monitored variables. There are nineteen
flights). Every pattern in the validation set belongs to a output layer neurons corresponding to the nineteen fault
known class. cases. Since the variation of the hidden neuron number
Figures 3-5 exemplify the fault class formation using produces different results, it is necessary to find the opti-
patterns of ZT in the space of two normalized devia- mal setting by trial and error that yields the maximal
tions: Figure 3 shows class 1 (no fault) and class 4 (HPC diagnosis accuracy. Another parameter to set is the num-
fault); Figure 4 shows faults 2-8 including component ber of epochs, which is the number of iterations to upda-
and actuator faults; Figure 5 shows sensor faults 10-14. te weights and biases. The network training also needs
some additional parameters such as the type of activa-
tion function and the back-propagation method.
Diagnosis accuracy
The fault recognition method selected classifies each Diagnostic algorithm results
pattern of the set ZV, producing the diagnosis dl. Com-
paring dl with a known class Dn for all validation set Errors of model variations
patterns, a confusion matrix is generated. Its diagonal is
formed by correct pattern classification probabilities As shown in (Cisneros et al., 2015), Variations 1 and 3
per class. A mean number P is obtained from these pro- contain displacements in their plots of errors δY at cer-
babilities representing the total diagnosis accuracy of tain intervals due to the influence of each engine. This
the gas turbine fault recognition. situation was addressed in Subsection deviation com-
putation to correct individualities of all engines. Figure
Fault recognition method 7 shows the errors of Variation 2 for 5000 flights and
one monitored variable. Here, the engine degradation
The fault recognition method chosen for this work is effect is observable through all flights. Each engine in
the Multi-Layer Perceptron (MLP). The MLP is an arti- the fleet experiences this inevitable situation.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM 81
Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
Figure 3. Class 1 (no fault) and class 4 Figure 4. Component and actuator faults Figure 5. Sensor faults
(HPC fault)
The total error e for each of the three variations
was more execution time. For this reason Variation 1 is the
calculated using the RMS of the vectors δ Y for all engi- model selected for the gas turbine diagnostic algorithm.
nes and flights. Table 6 shows these results. Variation 2
has the greater error (a mean for all monitored variables Fault diagnosis accuracy
of 0.0114). The reason of this could be because the de-
gradation is greater using 5000 flights. Variation 1 and In order to have the highest diagnosis accuracy, two
Variation 3 yield the lowest errors (a mean of 0.0021 and MLP parameters were tuned: the number of hidden la-
0.0023 respectively). Since the errors remain very close, yer neurons and the number of epochs. After perfor-
any variation of the two can be chosen. However, Varia- ming different computations, 54 neurons and 2000
tion 3 works with more data (300 engines x 270 flights) epochs produced the maximal validation probability.
than Variation 1 (100 engines x 90 flights) resulting in Figure 8 shows an example of this tuning. Figure 9 pre-
82 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de J esús , Pérez-Ruiz Juan Luis
Table 6. Total error e for all model variations and monitored variables
Baseline model Y1 Y2 Y3 Y4 Y5 Y6 Y7 Mean
Variation 1 0.0008 0.0010 0.0038 0.0003 0.0012 0.0024 0.0054 0.0021
Variation 2 0.0012 0.0026 0.0276 0.0014 0.0054 0.0296 0.0123 0.0114
Variation 3 0.0008 0.0010 0.0040 0.0004 0.0013 0.0031 0.0052 0.0023
sents the diagnosis probability of each fault class for 58.25% of recognition respectively) while the lowest
training and validation sets. Table 7 shows these values probabilities are produced by classes 8, 10 and 19
as well. For validation, the higher probabilities are ob- (9.25 %, 11.75 % and 1.25 % respectively). The differen-
tained from classes 2, 6 and 7 (60.50 %, 61.00% and ces between the probability values of both sets are ex-
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM 83
Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
plained by the limited pattern number of the validation The results showed that this network correctly per-
set. The increase of this number can produce more ac- forms this task; however, the diagnosis accuracy for
curate and closer results. The total diagnosis accuracy both sets (training and validation) seems to be relati-
for training is 40.66% and 38.29% for validation. P vely low. Some objective reasons of this could be: the
great number of fault classes and the low fault severi-
Conclusions ty randomly assigned to them. The other possible ex-
planation is that we did not take into consideration all
The data generation using the software ProDiMES peculiarities of the engine simulation in ProDiMES
allowed the simulation of healthy and faulty condi- because we are not authors of this simulator.
tions of an engine fleet in an appropriate environment Future works can consider working with different
facilitating the diagnostic process. This enabled the recognition methods such as Radial Basis Network,
test of our diagnostic approach on the data of a new Probabilistic Neural Network or Support Vector Machi-
engine. The stage of baseline model testing permitted nes in order to ensure that the fault recognition process
finding the best healthy engine performance approxi- was carried out correctly. The major issue is to increase
mation. This stage in an important part of the diagno- the diagnosis accuracy.
sis process because it directly affects the diagnosis
accuracy. After performing the calculations, the con- Acknowledgments
clusion is that the baseline model selected has low le-
vel of errors and the deviations computed with this The work has been carried out with the support of
model adequately reflect engine health conditions. In the National Polytechnic Institute of Mexico (IPN)
the gas turbine fault recognition stage, the Multi-La- and the National Council of Science and Technology
yer Perceptron was used to classify the fault patterns. (CONACYT).
Nomenclature
C-MAPSS Commercial Modular Aero-Propulsion Y Vector of monitored gas path variables
System Simulation
EFS Engine Fleet Simulation Yr, Yr Scalar and vector of individual baseline values
GUI Graphical User Interface Yϕ, Y ϕ Scalar and vector of general baseline values
HPC High Pressure Compressor Z, Z Scalar and vector of normalized deviations /
neuron input vector
HPT High Pressure Turbine ZT, ZV Training and validation sets
LPC Low Pressure Compressor γ Flow capacity
LPT Low Pressure Turbine ΔY, Δ Y Scalar and vector of initial deviation values
MLP Multi-Layer Perceptron δY, δ Y Scalar and vector of model errors
VBV Variable Bleed Valve δYr, δYr Scalar and vector of final deviations
VSV Variable Stator Vane δY ϕ Deviation correction vector
a ,a Layer output vectors h Efficiency
1 2
b1 , b2 Vectors of bias values s Standard deviation
C Matrix of model coefficients sY Mean deviation error
84 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM
Felipe-Altamirano Salvador, Loboda Igor, Cisneros-Azuara Felipe de J esús , Pérez-Ruiz Juan Luis
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM 85
Gas Turbine Diagnostic Algorithm Testing Using the Software ProDiMES
86 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 75-86 ISSN 1405-7743 FI-UNAM
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 87-99
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: febrero de 2016, aceptado: septiembre de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Ponce-Ceja Nohemí
Airbus Helicopters México Querétaro, S.A. de C.V.
Correo: noponce@[Link]
Resumen
La infraestructura ferroviaria es un elemento crítico para la efectividad de las cadenas de suministro de las empresas
automotrices, así como para la competitividad de los países en donde se localizan. En México, el sector automotriz
es uno de los sectores industriales estratégicos de la economía. A pesar de ello, hasta ahora no existe evidencia con-
sistente sobre cómo los flujos de carga automotrices presionan a la infraestructura ferroviaria de México como resul-
tado de la fuerte competencia global. En este sentido, la modelación espacial de los flujos de carga es una parte
esencial del proceso de planeación del transporte. El objetivo de este artículo es presentar un análisis espacial de los
flujos de carga automotrices desde una perspectiva de fluidez en la cadena de suministro. La evidencia revela que a
pesar de que México es actualmente un elemento clave del sistema de producción automotriz global, los flujos de
carga automotriz están extremadamente delimitados a las interacciones entre las plantas de manufactura localizadas
dentro de la región del Tratado de Libre Comercio para América del Norte (TLCAN). Como resultado, se direcciona
el proceso de planeación del transporte en México hacia un fortalecimiento de los corredores de carga del TLCAN.
Flujos de carga automotriz y su impacto en la infraestructura ferroviaria en México: un enfoque de fluidez en la cadena de suministro
De esta forma, el proceso de planeación deja de lado su rol crítico como soporte a la configuración de agrupamientos
industriales locales con mayor valor agregado. Este caso representa claramente el riesgoso dilema de las economías
emergentes cuando buscan mejorar el desempeño de las cadenas de suministro de las empresas globales instaladas
en sus territorios y, al mismo tiempo, integrar mayor valor agregado local.
Descriptores: infraestructura, transporte de carga, industria automotriz, organización industrial, fluidez en la cade-
na de suministro.
Abstract
Railway infrastructure is a critical element for automotive companies’ supply chain effectiveness as well as for countries’ com-
petitiveness where they are located. In Mexico, automotive sector is a strategic industrial sector of its economy. Even though,
until now there is not robust evidence about how automotive freight flows stress railway infrastructure in Mexico as a result of
its fierce global competition. In that sense, spatial modelling of freight flows is an essential part of the transportation planning
process. Thus, the objective of this paper is to present a spatial analysis of automotive freight flows based on the fluidity supply
chain approach. Evidence revealed that even if Mexico is currently a key member of the global automotive production system,
automotive freight flows are extremely circumscribed to interactions between manufacturing facilities located inside the NAFTA
(North American Trade Agreement) region. Consequently, they are driving the Mexican transportation planning process to
reinforce NAFTA transport freight corridors. Hence, this planning process is leaving aside its critical role in supporting con-
figuration of regional industrial clusters with more local value added. This case clearly represents the risky dilemma of emerging
economies when looking for improving global supply chains performance and at the same time, integrating local value added.
Keywords: infrastructure, freight transport, automotive industry, industrial organization, supply chain fluidity.
88 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM
Cedillo-Campos Miguel G., García-Ortega M.G., Martner-Peyrelongue C.D., Saucedo-Martínez Jania A., Ponce-Ceja N.
ma de organización de un sector clave para la competi- dad del sector, a través de la incorporación de más em-
tividad de México, así como la dinámica logística de sus presas mexicanas a las cadenas globales de valor. En
flujos. En la sección 3 se expone el enfoque metodológi- primer lugar, fortalecer el mercado interno mediante
co a través del cual se analizaron los flujos y sus impac- acciones como el ordenamiento del mercado de autos
tos en la infraestructura de transporte. En la sección 4 se usados importados y el mejoramiento del financia-
presentan los análisis que muestran que las cadenas de miento para la adquisición de vehículos nuevos. En se-
suministro automotriz tienen diferente grado de inte- gundo lugar, construir un entorno de negocios
gración con el contexto local y, en consecuencia, de uso alentando la productividad a través del desarrollo de
de la infraestructura de transporte. Ello representa no proveeduría local para la industria terminal y de auto-
solo un reto cada vez más crítico para el correcto diseño partes por medio del programa ProAuto. Dicho progra-
de corredores de transporte de carga, sino también para ma tiene por objeto el desarrollo de una red nacional de
su desempeño logístico. Finalmente, en la sección 5 se proveedores para aumentar el contenido nacional. En
expone evidencia revelando que a pesar de que México tercer lugar, impulsar el desarrollo de esquemas “triple
es actualmente un elemento clave del sistema de pro- hélice” vinculando al sistema de investigación y desa-
ducción automotriz global, los flujos automotrices de rrollo con los tomadores de decisión de gobierno y em-
carga están extremadamente delimitados a las interac- presas automotrices. Finalmente, en cuarto lugar,
ciones entre las plantas de manufactura localizadas diversificar los mercados de exportación, donde la mira
dentro de la región del Tratado de Libre Comercio para está en los países que han presentado las tasas de creci-
América del Norte (TLCAN). Como resultado, se direc- miento más dinámicas, como los integrantes de la
ciona el proceso de planeación del transporte en Méxi- Alianza del Pacífico y aquellos miembros del Acuerdo
co hacia el fortalecimiento de corredores de carga, que de Negociación Transpacífico.
son ejes del comercio del TLCAN, en detrimento de una En realidad, desde la conformación del TLCAN, la
integración regional-local enlazada a esos corredores producción anual de vehículos en México va en aumen-
de producción y transporte. to, con excepción de 2009, año de crisis económica. Para
2014, debido a que la producción automotriz de Brasil
Antecedentes estuvo por debajo de la producción que alcanzó en
2013, México terminó el año como líder en América La-
Para la mayor parte de las organizaciones industriales tina con una producción de alrededor de 3,300,000 uni-
localizadas en países emergentes, más que el impulso a dades, de la cual, 79% se destinó al mercado de
la integración de valor agregado local, el control de los exportación (figura 1). De hecho, entre enero y octubre
costos operativos es el principal objetivo. La transferen- de 2014, el país fue el séptimo mayor fabricante de ve-
cia de los conceptos de la “Manufactura Esbelta” al área hículos ligeros, el primero en América Latina, el cuarto
de logística y cadena de suministro ha hecho posible el exportador en el mundo, el sexto mayor exportador de
desarrollo de innovaciones basadas en la eficiencia de autopartes y el segundo proveedor de automóviles li-
las operaciones. Sin embargo, el incremento actual de la geros a Estados Unidos, después de Canadá (OICA,
volatilidad de los mercados, así como de la demanda, 2015).
aunado al incremento de los riesgos de interrupciones A través del tiempo, el patrón de localización de la
en los flujos de suministro (ya sea por causa de fenóme- industria automotriz ha variado según la política eco-
nos naturales o intencionalmente por el hombre) guían nómica seguida en el país. De este modo, a partir de la
el diseño de cadenas de suministro hacia la búsqueda firma del TLCAN, la industria automotriz adoptó una
de una mayor robustez que mejore su fluidez (Gaonkar nueva distribución espacial en donde las empresas ar-
y Viswanadham, 2004; Giunipero y Eltantawy, 2004; madoras se concentran tanto en la zona central del país,
The White House, 2012; Bueno y Cedillo, 2014; Cedillo como en la zona norte. Por otro lado, la industria de
et al., 2014a; Cedillo et al., 2014b). Donde fluidez en la autopartes, muestra una distribución más extendida,
cadena de suministro se entiende como: “la capacidad donde su mayor presencia se observa en torno a las
para alcanzar de forma consistente un flujo de procesos con- plantas armadoras. En la conformación de esta geogra-
fiable, seguro y preciso, apoyando con efectividad los objetivos fía productiva, un factor de localización fundamental es
de la cadena de suministro” (Cedillo et. al., 2015a). la disponibilidad de infraestructura ferroviaria y de
Dentro del contexto automotriz mexicano, en octu- servicios que atiendan las necesidades logísticas de la
bre de 2013 se anunció a través del Secretario de Econo- industria automotriz. Al respecto, el presidente de la
mía, la implementación de cuatro estrategias orientadas Asociación Mexicana de Ferrocarriles (AMF) expresó
a mantener la tendencia de crecimiento y competitivi- en febrero de 2014 que: “El servicio ferroviario está íntima-
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM 89
Flujos de carga automotriz y su impacto en la infraestructura ferroviaria en México: un enfoque de fluidez en la cadena de suministro
90 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM
Cedillo-Campos Miguel G., García-Ortega M.G., Martner-Peyrelongue C.D., Saucedo-Martínez Jania A., Ponce-Ceja N.
2014e; Avetisyan et al., 2015). Sin embargo, todo ello tie- La esencia de la metodología se sitúa en la identifi-
ne aún importantes retos por resolver. De hecho, a tra- cación y representación espacial de los flujos de carga
vés de un modelo de dinámica de sistemas, Cedillo et al. en la red ferroviaria del país. En general, las actuales
(2014a), exponen los críticos efectos de retraso en los aportaciones conceptuales y metodológicas en el tema
cruces de frontera, confirmando la importancia que tie- de los movimientos de carga se basan en modelos eco-
nen programas como el C-TPAT para hacer más eficien- nométricos que estiman volúmenes y tipos de mercan-
te la operación de las llamadas “NAFTA Supply Chains”. cías intercambiadas entre regiones, pero sin capacidad
para territorializar los flujos en espacios concretos, en
Enfoque metodológico redes y nodos específicos, o bajo determinadas combi-
naciones multimodales. La incorporación de la dimen-
Desde la perspectiva global de fluidez en la cadena de sión territorial implica la consideración de la morfología
suministro previamente definido y con base en una am- del espacio, la conectividad de la infraestructura y el
plia revisión de información bibliográfica, se desarro- origen y demanda de la carga.
llaron análisis espaciales. Lo anterior, a partir de un De manera específica, los insumos de la metodolo-
enfoque metodológico derivado de la necesidad de ge- gía fueron dos:
nerar herramientas y procedimientos que permitieran
sistematizar y representar espacialmente información, a) Bases de datos anuales de los movimientos ferrovia-
como por ejemplo, la matriz origen-destino de la carga rios de carga;
ferroviaria (García y Martner, 2011). b) Información georreferenciada de la red férrea del
La carencia de información para el desarrollo de análi- país (vías y estaciones).
sis geo-espaciales fue el elemento detonador para el ini-
cio del presente proyecto. Como consecuencia de la La construcción del soporte geográfico fundamental
desregulación del transporte ferroviario que inició en para la aplicación del proceso metodológico, constituyó
1996, se dejó de contar con información detallada por el conjunto de tareas más laborioso. Por un lado, inclu-
alrededor de 14 años1. Tiempo durante el que tampoco yó la incorporación de un importante número de esta-
se disponía de alguna herramienta metodológica que ciones no incluidas en la cobertura original, que a su
permitiera visualizar el movimiento de la carga en la vez, implicó la inclusión anual de nuevos nodos, ade-
red ferroviaria con fines de administración, planeación más de la depuración de la base de datos que se asocia
y diseño de políticas públicas. Por ello, el desarrollo de a la información georreferenciada. Por otro lado, invo-
un análisis orientado a sistematizar y representar terri- lucró la construcción de las condiciones de conectivi-
torialmente las bases de datos origen–destino de la car- dad de la red ferroviaria, lo cual es indispensable para
ga ferroviaria representó un importante avance en el que los trazos lineales que la representan se comporten
análisis sistemático del comportamiento, evolución y precisamente como una red que permita cumplir con
distribución de los flujos de carga en la red ferroviaria los propósitos de representar territorialmente los inter-
de México.2 cambios de carga entre pares de estaciones y los flujos
de carga que circulan a lo largo de los distintos tramos
1 Como consecuencia de la política de desregulación del transporte de la red ferroviaria nacional.
llevada a cabo por el gobierno mexicano en los años 90´s, las empre-
De este modo, una vez asociados los registros de
sas privadas escudadas en las ausencias y vacíos de la legislación, de-
origen y de destino de la carga a la base geográfica, se
jaron de generar o de reportar información al Estado. Este problema
parece no ser atípico, ni excepcional, de hecho lo advierte con preci-
realizaron en TransCAD3 tres procesos de trabajo:
sión el Sistema Nacional de Información Estadística y Geográfica de
INEGI quien reconoce cómo el problema de obligatoriedad y de nor-
malización en la generación de estadísticas en el país. En consecuen-
cia, la Secretaría de Comunicaciones y Transportes (SCT), no cuenta
con una base de datos estandarizada sobre el transporte ferroviario Comunicaciones y Transportes (SCT) y el Instituto Nacional de Esta-
de carga, vacío que da lugar a la dificultad de formular metodologías dística y Geografía (INEG), y en este marco se obtuvo la información
para el análisis estadístico de la información para generar indicadores para iniciar la integración y la representación geográfica de la matriz
que permitan la identificación de necesidades y la realización de pro- origen – destino de la carga ferroviaria.
nósticos que permitan orientar las labores de gestión en materia de 3 TransCAD, es un software diseñado para apoyar análisis diversos en
transporte ferroviario por parte de la autoridad del país. materia de transporte, combinando, capacidades de administración
2 Dicha necesidad se identificó en el seno del Comité Técnico Espe- de bases de datos, de mapeo digital y de presentación gráfica, con
cializado de Información Económica y Operativa del Sector Transpor- herramientas específicas para el análisis y modelado de sistemas de
te (CTE-IEOST), encabezado y conducido por la Secretaría de transporte, investigación de operaciones y de modelos estadísticos.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM 91
Flujos de carga automotriz y su impacto en la infraestructura ferroviaria en México: un enfoque de fluidez en la cadena de suministro
92 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM
Cedillo-Campos Miguel G., García-Ortega M.G., Martner-Peyrelongue C.D., Saucedo-Martínez Jania A., Ponce-Ceja N.
y sistemas de transporte y logística más aptos para res- Uno de los grandes aspectos que influyen en este tipo
ponder a estas necesidades, están dentro de los desafíos de comportamiento es la falta de confianza hacia la pro-
de mediano plazo de la industria automotriz y por su- veeduría local. No confían en que las empresas locales
puesto, de los países como México que acogen a esta cumplirán los altos estándares con que el sector auto-
industria. motriz mundial rige sus operaciones.
Asimismo, otro avance tecnológico que puede trans- En el caso específico del clúster automotriz localiza-
formar paradigmas en la industria automotriz es el uso do en la ciudad de Hermosillo, se encontró que peque-
de la impresión en tercera dimensión para la produc- ñas y medianas empresas comenzaron a participar en
ción de vehículos (Weller et al., 2015). Esta nueva tecno- una relación más cercana a los proveedores Tier 1 y la
logía conduciría a cambios hacia una producción armadora automotriz localizada en la ciudad, solo
localizada más cerca de su mercado final, dismuniyen- como resultado de que los dueños de las empresas par-
do el costo total de la producción (Gebler et al., 2014). ticipantes, habían trabajado previamente para los pro-
Pero además, las tecnologías pro ambientales están to- veedores y la armadora. En consecuencia, se tenía una
mando fuerza y empiezan a ser parte de las decisiones previa relación de confianza y conocimiento de los altos
estratégicas en la cadena de suministro automotriz, que estándares a los que dichas empresas proveedoras de-
buscan reducir su huella de carbono (Lee, 2011). La in- bían responder. Así también, las redes de confianza hi-
dustria automotriz actual trata de conformar un mode- cieron posible el intercambio y aprendizaje de ideas y
lo económico viable que, a su vez, tenga menores tecnología (Contreras et al., 2012). Sin embargo, se trata
impactos tanto en el medio ambiente como en la socie- de un esquema que lleva tiempo desarrollar.
dad (Nieuwenhuis y Katsifou, 2015). Queda claro que para mejorar la competitividad del
Por otro lado, la forma de operar de la industria au- sistema automotriz en México, es actualmente necesa-
tomotriz varía dependiendo su origen. Existen estudios ria la identificación de otros mecanismos de integración
que comparan la forma de trabajar de las armadoras de valor agregado local más dinámicos.
alemanas y japonesas, concluyendo que mientras las
plantas de origen japonés fabrican 3.8 distintos mode-
Impacto en la infraestructura de transporte
los por línea de ensamble, las alemanas, solo hacen 1.6
(Staeblein y Aoki, 2015). Con ello, la exigencia hacia los Tomando como base el contexto que ha favorecido la
sistemas logísticos locales dependiendo del origen na- localización de un importante número de empresas de
cional de la empresa, es también diferente. Asimismo, ensamble automotriz en México, se procedió a repre-
existen estudios comparativos sobre armadoras asiáti- sentar espacialmente la matriz O-D de los intercambios
cas, específicamente, entre Toyota y Hyundai en los ferroviarios de vehículos automotores. Lo anterior per-
que se reporta, que mientras Toyota se enfoca en la pla- mitió observar, por medio de las líneas de deseo que
neación y el sistema de trabajo para evitar que factores enlazan los nodos de emisión con los de recepción, que
externos influyan en sus procesos, Hyundai trabaja to- en primer lugar actualmente todas las zonas de produc-
mando en cuenta la incertidumbre externa en sus ope- ción automotriz se encuentran atendidas por el ferroca-
raciones para mantenerse flexible (Shim y Steers, 2012). rril. Las mismas también se encuentran conectadas por
Los aspectos anteriores forman parte del conjunto de ferrocarril con los puertos fronterizos y marítimos vin-
elementos en juego al momento de seleccionar un país culados con el comercio exterior de México (figura 2).
una nueva planta de ensamble automotriz. Además de que la mayoría de estas registran intercam-
Una vez seleccionado el país de implantación, los bios ferroviarios mayores a 50 mil toneladas (García, y
procesos de decisión se enfocan en la selección de la Martner, 2014). Como se puede observar en la figura 2,
región-local en donde las operaciones tendrán lugar. los principales intercambios representados por medio
Para el caso de México, diferentes estudios reportan de las líneas de deseo son poco densos. Esto se debe a la
que la relación que mantiene gran parte de las armado- especializada localización de la industria automotriz,
ras con el medio local es de “oportunidad”. Es decir, una además de la ubicación también muy específica de los
búsqueda por la reducción constante de costos. Cuando centros de mercado, de los centros de distribución na-
las armadoras se instalan en el país, atraen a sus pro- cional y de los puertos de comercio exterior, que son los
pios proveedores por lo que la relación con las empre- tres tipos de espacios involucrados en los movimientos
sas locales es regularmente muy distante. Se trata de este tipo de carga.
fundamentalmente de una relación de oportunidad en Como testimonios de la interdependencia económica
dónde la mejora operativa se centra en la disminución entre el ferrocarril y las armadoras establecidas en Méxi-
de costos (Cedillo y Pérez, 2010; Medeiros et al., 2011). co se puede citar que el segmento automotriz se convir-
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM 93
Flujos de carga automotriz y su impacto en la infraestructura ferroviaria en México: un enfoque de fluidez en la cadena de suministro
94 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM
Cedillo-Campos Miguel G., García-Ortega M.G., Martner-Peyrelongue C.D., Saucedo-Martínez Jania A., Ponce-Ceja N.
1,600,000
1,400,000
1,200,000
1,000,000
800,000
600,000
400,000
200,000
0
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM 95
Flujos de carga automotriz y su impacto en la infraestructura ferroviaria en México: un enfoque de fluidez en la cadena de suministro
La relación entre la industria automotriz y el trans- sólo se concluye el producto en territorio nacional, ge-
porte ferroviario es altamente benéfica para la produc- nerando empleo pero cuya derrama económica no al-
tividad de ambas partes. El crecimiento anual de la canza a crear mayores espacios internos en cascada
producción automotriz en el país ha impulsado las in- para el desarrollo económico y social.
versiones de las empresas ferroviarias en equipos de La expresión espacial de la industria automotriz es
arrastre especializados para vehículos automotores y el precisamente la de un espacio transnacionalizado, cuya
incremento en la oferta de servicios logísticos de estas y geografía la establece la propia cadena de suministro y
de los operadores intermodales. adquiere, por tanto, la forma de una red determinada
por los movimientos de intercambio. Se tienen flujos de
Conclusiones carga automotriz organizados a lo largo de grandes co-
rredores de carga, donde la capilaridad con el espacio
Debido a la importancia de los flujos de carga automo- local queda en entredicho, formando espacios sub-na-
trices en México, se está direccionando el proceso de cionales en la economía global. Su impacto para la in-
planeación del transporte en México hacia un fortaleci- fraestructura de transporte en México es regionalmente
miento de los corredores de carga del TLCAN para in- desigual y altamente enfocado a la interacción con el
crementar la fluidez en las cadenas globales de mercado de consumo ubicado en Estados Unidos.
suministro automotriz. Sin embargo, al no contar Méxi- En consecuencia, si bien se trata de un sector estraté-
co con una política pública exitosa de impulso al desa- gico (más si se considera que es uno de los pilares del
rrollo e integración de proveedores locales, este proceso desarrollo económico de la región del TLCAN), se ha
de planeación corre el riesgo de dejar de lado su rol crí- identificado como crítico el establecimiento de los pará-
tico como soporte a la configuración de agrupamientos metros que favorezcan la integración de las necesidades
industriales con mayor valor agregado local. Este caso de conexión de los grandes corredores de carga con el
representa claramente el riesgoso dilema de las econo- espacio regional-local. Todo ello con el fin de robustecer
mías emergentes cuando buscan mejorar el desempeño la red de infraestructura de transporte a fin de que se
de las cadenas de suministro de las empresas globales convierta en una palanca para el desarrollo de proveedo-
instaladas en sus territorios y al mismo tiempo, integrar res locales y de otros sectores clave de la economía.
mayor valor agregado local. Como parte del trabajo futuro, se desarrollará un
Las características de este caso se resumen en la for- análisis comparativo con los requerimientos de otros
mación de espacios geográficos puntuales, especializa- sectores productivos más vinculados con el desarrollo
dos y de alto desarrollo tecnológico que aprovechan las regional local, como es el caso del sector agroalimenta-
ventajas comparativas de vecindad entre México y Es- rio. El objetivo será determinar las características de sus
tados Unidos; la posibilidad de acceso de productos a impactos y necesidades logísticas de forma que sean
este mismo país, vía el TLCAN, así como el costo favo- considerados como parte de las políticas públicas que
rable de la mano de obra. busquen diseñar infraestructura con alta conectividad
Desde luego, también existen otras ventajas compe- logística, estableciendo un equilibrio con los requeri-
titivas en desarrollo, entre ellas, la creciente calificación mientos de otros sectores de la economía y en conse-
de la mano de obra mexicana, el aumento en la disponi- cuencia, favoreciendo un desarrollo regional equilibrado.
bilidad de servicios ferroviarios diseñados a la medida
de la industria automotriz, así como la mejora en la co- Agradecimientos
nexión ferroviaria con puertos marítimos y pasos fron-
terizos en donde la red férrea de México mejora su Los autores agradecen al Laboratorio Nacional en Siste-
conectividad con la red férrea de Estados Unidos. mas de Transporte y Logística del Instituto Mexicano
Sin embargo, esta industria y sus cadenas de sumi- del Transporte (IMT) el amplio apoyo brindado para el
nistro operan hasta ahora en México, más como encla- desarrollo de la presente investigación.
ves de producción que como motores de desarrollo
regional. La producción de autopartes hasta ahora es Referencias
mayoritariamente externa, el contenido nacional de los
vehículos hechos en México es de 33%. Se trata primor- Ahmed M. y Hendry L. Supplier development literature review and
dialmente de flujos compuestos por flujos de ingreso de key future research areas. International Journal of Engineering and
autopartes (importaciones) y flujos de salida de vehícu- Technology Innovation, volumen 2 (número 4), 2012: 293-303.
los terminados (exportación). En otras palabras, habla- Ascencio L., González-Ramírez R., Bearzotti L., Smith N., Cama-
mos de un proceso de producción exógeno en donde cho-Vallejo J. A collaborative supply chain management system
96 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM
Cedillo-Campos Miguel G., García-Ortega M.G., Martner-Peyrelongue C.D., Saucedo-Martínez Jania A., Ponce-Ceja N.
for a maritime port logistics chain. Journal of Applied Research 2015]. Disponible en: [Link]
and Technology, volumen 12 (número 3), 2014: 444-458. 02/10/ferrocarriles-clave-atraer-inversion-automotriz
Avetisyan M., Heatwole N., Rose A., Roberts B. Competitiveness De la O.G. y Matis T. Supply chain strategies for the international
and macroeconomic impacts of reduced wait times at US land interoceanic Mazatlan-Houston logistics corridor. Journal of
freight border crossings. Transportation Research Part A: Policy Applied Research and Technology, volumen 12 (número 4), 2014:
and Practice, volumen 78, 2015: 84-101. 167-179.
Azevedo S., Govindan K., Carvalho H., C.-M. V. An integrated Duran-Fernández R. y Santos G. Introduction to the special issue
model to assess the leanness and agility of the automotive in- regional development and transport infrastructure in Mexico.
dustry. Resources, Conservation and Recycling, volumen 66, Research in Transportation Economics, volumen 46, 2014a: 1-2.
2012: 85-94. Duran-Fernandez R. y Santos G. Regional convergence, road in-
Becerril-Torres O., Álvarez-Ayuso I., Del Moral-Barrera L. Do in- frastructure, and industrial diversity in Mexico. Research in
frastructure influence the convergence of efficiency in Mexico. Transportation Economics, volumen 46, 2014b: 103-110.
Journal of Policy Modeling, volumen 32 (número 1), 2010: 120- Duran-Fernandez R. y Santos G. A regional model of road acces-
137. sibility in Mexico: accessibility surfaces and robustness analy-
Brock G., y German-Soto V. Regional industrial growth in Mexico: sis. Research in Transportation Economics, volumen 46, 2014c:
do human capital and infrastructure matter? Journal of Policy 55-69.
Modeling, volumen 35 (número 2), 2013: 228-242. Duran-Fernandez R. y Santos G. Road infrastructure spillovers on
Bronk C., y González-Aréchica B. Mexico-United States border the manufacturing sector in Mexico. Research in Transportation
security: From a bilateral to a truly binational policy process. Economics, volumen 46, 2014d: 17-29.
Latin American Policy, volumen 2 (número 2), 2011: 152-181. Duran-Fernandez R. Infrastructure policy in the USA-Mexico bor-
Bueno A. y Cedillo-Campos M. Dynamic impact on global supply der: Evaluation and policy perspectives. Research in Transpor-
chains performance of disruptions propagation produced by tation Economics, volumen 46, 2014e: 70-102.
terrorist acts. Transportation Research Part E: Logistics and Francois J., y Manchin M. Institutions, infrastructure, and trade.
Transportation, volumen 61, 2014: 1-12. World Development, volumen 46, 2013: 165-175, 65-175.
Cedillo-Campos M. y Pérez P. Hybrid supply chains in emerging Fogliatto F., Da Silveira G., Borestein D. The mass customization de-
markets: the case of the Mexican auto industry. South African cade: an updated review of the literature. International Journal of
Journal of Industrial Engineering, volumen 21 (número 1), 2010: Production Economics, volumen 138 (número 1), 2012: 14-25.
193-206. Gaonkar R. y Viswanadham N. A conceptual and analytical fra-
Cedillo-Campos M., Sanchez C., Vadali S., Villa J., Menezes M. mework for the management of risk in supply chains. Procee-
Supply chain dynamics and the “cross-border effect”: The dings, ICRA, IEEE, International Conference on Robotics and
U.S.-Mexican border’s case. Computers and Industrial Enginee- Automation, volumen 3, 2004, pp. 2699-2704.
ring, volumen 72, junio 2014a: 261-273. García G. y Martner C. Metodología para integrar una matriz ori-
Cedillo-Campos M., Bueno A., González R., Jiménez E., Pérez G. gen-destino de la carga ferroviaria. Fase I: carga conteneriza-
Dynamic impact on global supply chains performance caused da. Reporte de Investigación ME 14/11, Instituto Mexicano del
by criminal disruptions. Journal of Applied Research & Technolo- Transporte, México, 2011.
gy, volumen 12 (número 4), 2014b: 684-694. García G. y Martner C. Movimientos ferroviarios por puertos y
Cedillo-Campos M. y Cedillo-Campos H. w@reRISK Method: a fronteras en México. Reporte de Investigación TE 22/14, Insti-
tool for security risk level classification of products. Safety tuto Mexicano del Transporte, México, 2014.
Science, volumen 79, 2015a: 358-368. Giunipero L.C. y Eltantawy R.A. Securing the upstream supply
Cedillo-Campos M., Lizarraga-Lizarraga G., Martner-Peyrelon- chain: a risk management approach. International Journal of
gue C. MiF3 method: Measurement of intermodal fluidity Physical Distribution & Logistics Management, volumen 9 (nú-
freight flows. Reporte técnico, Instituto Mexicano del mero 34), 2004: 698-713.
Transporte, México, 2015b. Gebler M., Schoot A., Visser C. A global sustainability perspective on
Contreras O., Carrillo J., Jorge A. Local entrepreneurship within 3D printing technologies. Energy Policy, volumen 74, 2014: 158-167.
global value chains: a case study in the Mexican automotive González A. Amarra México liderato de producción automotriz
industry. World Development, volumen 40 (número 5), 2012: en América Latina, El Economista, 2014, [fecha de consulta: 14
1013-1023. de noviembre de 2015] [en línea]. Disponible en: [Link]
Cravioto J., Yamasue E., Okumura H., Ishihara K. Road transport [Link]/industrias/2014/12/04/mexico-sera-primer-
externalities in Mexico: estimates and international compari- productor-autos-2014
son. Transport Policy, volumen 30, 2013: 63-76. Hockman G., Tabakis C., Zilberman D. The impact of internatio-
Cruz O. Ferrocarriles, clave para atraer inversión automotriz, Re- nal trade on institutions and infrastructure. Journal of Compa-
vista T21, 2014 [en línea] [fecha de consulta: 9 de julio de rative Economics, volumen 41 (número 1), 2013: 126-140.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM 97
Flujos de carga automotriz y su impacto en la infraestructura ferroviaria en México: un enfoque de fluidez en la cadena de suministro
Hu S. Evolving paradigms of manufacturing: from mass produc- de 2015] [en línea]. Disponible en: [Link]
tion to mass customization and personalization, Forty sixth gory/production-statistics/
CIRP Conference on manufacturing systems, Procedia CIRP Romero I. Impacto asimétrico de la crisis global sobre la industria auto-
7, 2013, pp. 3-8. motriz: Canadá y México comparados. Perspectivas para el futuro,
INA-Industria Nacional de Autopartes, A.C. Evolución y perspec- Comisión Económica para América Latina y el Caribe
tivas del sector automotriz, 2015 [fecha de consulta: 9 de julio (CEPAL), Organización de las Naciones Unidas, México, 2011.
de 2015] [en línea]. Disponible en: [Link] Shim W. y Steers R. Symmetric and asymmetric leadership cultu-
FTP/INA-Espa%C3%B1ol-completa..pptx res: a comparative study of leadership and organizational cul-
Johansson M. y Johnsson A. The package logic: a study on value ture at Hyundai and Toyota. Journal of World Business,
creation and knowledge flows. European Management Journal, volumen 47 (número 4), 2012: 581-591.
volumen 30 (número 6), 2012: 535-551. Staeblein T. y Aoki K. Planning and scheduling in the automotive
Keeling D. Transport research challenges in Latin America. Jour- industry: a comparison of industrial practice at german and
nal of Transport Geography, volumen 29, 2013: 103-104. japanese makers. International Journal of Production Economics,
Lee K.H. Integrating carbon footprint into supply chain manage- volumen 162, 2015: 258-272.
ment: the case of Hyundai Motor Company (HMC) in the au- The White House. National Strategy for Global Supply Chain Securi-
tomobile industry. Journal of Cleaner Production, volumen 19 ty, Washington, U.S.A., 2012.
(número 11), 2011: 1216-1223. Tomlin B. Managing supply-demand risk in global production:
Lim L., Alpan G., Penz B. Reconciling sales and operations mana- creating cost-effective flexible networks. Business horizon, vo-
gement with distant suppliers in the automotive industry: a lumen 57 (número 4), 2014: 509-519.
simulation approach. International Journal of Production Econo- Tseng M., Jiao R., Wang C. Design for mass personalization. CIRP
mics, volumen 151, 2014: 20-36. Annals-Manufacturing Technology, 2010: 175-178.
Martínez-Jurado P.J. y Moyano-Fuentes J. Lean management, Weller C., Kleer R., Piller F. Economic implications of 3D printing:
supply chain management and sustainability: a literature re- Market structure models in light of additive manufacturing
view. Journal of Cleaner Production, volumen 85, 2014: 134-150. revisited. International Journal of Production Economics, volu-
Medeiros G., Selitto M., Borchard M., Albert G. Procurement cost men 164, 2015: 43-56.
reduction for customized non-critical items in an automotive Wilsmeier G., Monios J., y Rodríguez J.P. Driver for outside-in
supply chain: an action research project. Industrial Marketing port hinterland integration in Latin America: the case of Vera-
Management, volumen 40 (número 1), 2011: 28-35. cruz, Mexico. Research in transportation Business & Manage-
Nieuwenhuis P. y Katsifou E. More sustainable automotive pro- ment, volumen 14, 2015: 34-43.
duction through understanding decoupling points in leagile Villa J. y Sacristán-Roy E. Privatization of Mexican railroads: fif-
manufacturing. Journal of Cleaner Production, volumen 95, teen years later. Research in Transportation Business & Manage-
2015: 232-241. ment, volumen 6, 2013: 45-50.
OICA-International Organization of Motor Vehicle Manufactu-
rers, Production statistics, 2015 [fecha de consulta: 9 de julio
98 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM
Cedillo-Campos Miguel G., García-Ortega M.G., Martner-Peyrelongue C.D., Saucedo-Martínez Jania A., Ponce-Ceja N.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 87-99 ISSN 1405-7743 FI-UNAM 99
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 101-109
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: junio de 2016, aceptado: octubre de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Sakhavand Navid
Rice University, Houston, Texas
Department of Civil and Environmental Engineering
E-mail: navids@[Link]
Kochandra Raji
National Institute of Technology, Calicut, India
School of Nano Science and Technology
E-mail: [Link]@[Link]
Ajayan Pulickel M.
Rice University, Houston, Texas
Department of Materials Science and NanoEngineering
E-mail: ajayan@[Link]
Shahsavari Rouzbeh
Rice University, Houston, Texas
Department of Civil and Environmental Engineering
E-mail: rouzbeh@[Link]
Abstract
Nanofluids have been recently studied for efficient and high impact dielectric and tribological fluids, which possess
superb thermal transport properties, compared to conventional fluids or lubricants. The present research theoretica-
lly analyzes, through Molecular Dynamic (MD) simulations, the physical performance of hexagonal boron nitride
(h-BN) nanosheets reinforced mineral oil (MO). This analysis leads to the prediction of nanofluids performance, uti-
lizing the Lennard-Jones potential. The computational results of viscosity performance are benchmarked at various
temperature levels and nanofiller concentrations. The theoretical results indicate that the simulations presented can
predict the behavior of nanoreinforced fluids.
Keywords: molecular dynamics, nanofluids, hexagonal boron nitride, viscosity.
Resumen
Los nanofluidos recientemente se han estudiado para fluidos dieléctricos y tribológicos eficientes y de alto impacto, los cuales
poseen excelentes propiedades de transporte térmico, comparado con fluidos o lubricantes convencionales. La presente investi-
gación analiza teóricamente, a través de simulaciones de Dinámica Molecular (MD), el desempeño físico de aceite mineral (MO)
reforzado con nanohojuelas de nitruro de boro hexagonal (h-BN). Este análisis permite la predicción del desempeño de los nano-
fluidos, utilizando el potencial Lennard-Jones. Los resultados computacionales del desempeño en viscosidad se han referenciado
en varios niveles de temperatura y concentraciones de los nanofluidos. Los resultados teóricos indican que las simulaciones
presentadas pueden predecir el comportamiento de los fluidos nano-reforzados.
Descriptores: dinámica molecular; nanofluidos; nitruro de boro hexagonal; viscosidad.
Theoretical Prediction of Physical Properties (Viscosity) on 2D-based Nanofluids
Introduction man et al., 2001; Choi et al., 2008; Gao et al., 2009; Xie and
Chen, 2011; Baby and Sundara, 2011; Botha et al., 2011;
Materials performance has become more relevant in Balandin, 2011, Choi, 2008). Other factors influencing
today’s industrial applications. A prime aspect that has the thermal transport performance of nanofluids are
been critical through years is the importance on redu- the liquid layering and viscosity. The viscosity of na-
cing the inherently losses from friction and wear, nofluids is another important transport property for
which reduce the performance of metal-mechanic pro- applications of nanofluids as a new class of HTFs in
cesses (Tzanakis et al., 2012). Friction plays an impor- thermal devices or systems such as heat exchangers or
tant role in diverse processes such as drilling, cutting, cooling systems. The improvement in thermal conduc-
working pair components and mechanisms, among tivity cannot be achieved by increasing the solid filler
others; causing wear and energy dissipation, which amount beyond a limit, since increase in filler concen-
subsequently impact on economic aspect (Jost, 1990). tration will increase the viscosity which will adversely
Wear is the major cause of energy and material loss in affect the fluid properties.
mechanical processes, as components are in constant Recently, two dimensional (2D) nanomaterials got
friction. However, with utilization of lubricants, com- remarkable scientific attention due to their unique phy-
ponents friction and wear could be minimized, resul- sical properties and high surface area (Coleman et al.,
ting into considerable energy and tooling savings 2011). The bulk counter parts of some of these 2D mate-
(Hernández et al., 2008). Nevertheless, while two or rials such as graphene and hexagonal boron nitride
more moving components are in contact, friction gene- (h-BN) are already known for their lubrication proper-
rates heat, which could degrade these lubricants; hen- ties for a long time. h-BN is an effective solid lubricant,
ce, this heat should be dissipated rapidly. it is suitable for diverse applications such as metal-wor-
With aid of nanotechnology, recent advances in tri- king processes where lubrication at high-temperatures
bology and heat transport have been observed. For is required and is widely used in high-temperature
instance, diverse studies have demonstrated that wear sealing materials of aerospace engines as well as
addition of nanoparticles to conventional lubricants anti-corrosive compound (Taha et al., 2013; Kimura et
can improve their tribological properties (Hernández al., 1999; Pakdel et al., 2011; Husain et al., 2013; Yi et al.,
et al., 2008; Lin, 2007; Taha, 2013; Chang et al., 2010; 2014; Cho et al., 2013; Golberg et al., 2010; Haubner et al.,
Hernández et al., 2010; Lee et al., 2012). The addition of 2002; Lin and Connell, 2012; Eichler and Lesniak, 2008;
nanoparticles to lubricants and heat-transfer fluids Mohan et al., 2013).
(HTFs) has shown superb results in tribological tests, Diverse investigations on nanofluids have been con-
as these could be deposited on the contact surfaces ducted, for instance, Hernández et al. (2010) studied the
and contribute to friction and wear reduction, also im- friction behavior of NiCrBSi coatings lubricated by CuO
proving load carrying capacity (Peña et al., 2015, 2014; nanoparticles suspended in poly-alphaolephin (PAO6).
Eswaraiah et al., 2011; Yu et al., 2008; Wu et al., 2007; Block-on-ring tests showed a decrease of up to 100% in
Peña et al., 2016). Similarly, several researchers have the friction coefficient with 2.0wt% CuO. Similarly, Yu et
also shown significant enhancements of reinforced al. (2008) reported improved lubricating properties by
HTFs with incorporation of nanoparticles in thermo- adding 0.2wt.% Cu nanoparticles to lubricant oil; in their
physical properties, such as thermal conductivity and study Cu formed a soft film by friction-shearing and
viscosity changes involving nanoparticle concentra- high pressure reducing the coefficient of friction up to
tion and temperature changes (Taha et al., 2013, 2012, 20%. Wu et al. (2007) examined the tribological proper-
2014; Sudeep et al., 2014).. ties of nanofluids of API-SF oil (SAE30 LB51153) and
Most of these lubricants and HTFs have the disad- base oil, with sphere-like CuO, TiO2, and nano-diamond
vantage of having low thermal conductivity. But when nanoparticles. Their anti-wear behavior is attributed to
reinforced by these suspended ultrafine structures such the sphere-like morphology of the nanoparticles resul-
as ceramics, metals, semiconductors, carbon nanotubes ting in a rolling effect between the surfaces; and the de-
(CNTs), and composite materials with varied shapes position of the nanoparticles on the worn surface, de-
and sizes (<100nm) they possess great potential of creasing the shearing stress and friction coefficient.
application to improve energy efficiency and heat Timofeeva et al. (2007, 2009) studied the thermal
transfer, which could be used in diverse areas and conductivity and viscosity of Al2O3 nanoparticles (11nm
applications such as in power generation, biomedical, to 40nm in diameter) dispersed within water and
tribological, cooling, nuclear, defense, microelectronics, ethylene glycol (EG). It is shown that nanoparticles are
among others (Taha et al., 2014; Sudeep et al., 2014; East- highly agglomerated; agglomeration state varies in
102 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM
Taha-Tijerina José, Sakhavand Navid, Kochandra Raji, M. Ajayan Pulickel, Shahsavari Rouzbeh
time, as the sample ages as well. The authors stated that ty as well. Lubricants thermo-oxidative stability is es-
the main parameters for controlling nanofluids’ ther- sential to minimize its degradation during service or
mal conductivity enhancement are the geometry, storing. Also, changes in fluids viscosity and auto-oxi-
agglomeration state, and surface resistance of nanopar- dation are promoted due to the high temperatures rea-
ticles. On the other hand, Wang et al. (1999) investigated ched. These issues have negative effect in lubricant
diverse fluids (water, pump fluid, engine oil and EG) properties; hence, the thermo-stability of lubricants is
with the addition of Al2O3 and CuO with 28nm and an important property.
23nm in diameter, respectively. Wang et al. measured Improvement in thermal transport of nanofluids
the viscosity of these systems and found an increment cannot be achieved by just increasing the solid filler
as nanoparticles agglomerate, also thermal conductivi- amounts beyond a limit, since increase in filler concen-
ty performance is observed to decrease; due to an effect tration will increase the viscosity which will adversely
of the agglomeration of the nanoparticles. Moreover, affect the fluid properties. With the advancement of
particles agglomeration is exacerbated by the size of the computational power, simulation tools are increasingly
reinforced fillers. Furthermore, Nguyen et al. (2007) in- used for prediction of material (mechanical (Sakhavand
vestigated the particle size effect for Al2O3 aqueous ba- et al., 2013; Sakhavand and Shahsavari, 2014), thermal
sed nanofluids and observed that the particle size (Raji et al., 2012), electrical (Yuge et al., 2005), etc.) pro-
effects on viscosity are more significant for high parti- perties. The atoms are modeled individually in this
cles concentration. method. The dynamic evolution of the system is gover-
On thermal transport, Taha et al. (2012) investigated ned by the interaction of atoms and calculation of the
MO reinforced with 2D-nanostructures of h-BN and total potential energy of the system from bonded and
graphene at low filler fraction (<0.10wt.%), it is obser- non-bonded interactions.
ved that the viscosity of the nanofluids decreases signi- Molecular dynamic (MD) simulations yield infor-
ficantly with temperature, while the enhancement in mation at the atomic or molecular level behaviors,
viscosity with the addition of 2D-nanofillers is very which contribute to the averaged transport properties
small (< 2% at 313K). This is an additional advantage of of the entire fluid-nanoparticle system. The computa-
the low filler fractions since the increase in viscosity tional domain consists of h-BN nanosheets homoge-
will decrease the effective thermal conductivity values neously dispersed within mineral oil (MO), which
as well as flow characteristics of the fluid. Moreover, consists of paraffinic, isoparaffinic, naphthenic, aro-
the relatively small increase in viscosity (<30%) at matic and poly-aromatic compounds such as cyclo-
0.35wt.% of h-BN is evidence that the solution is not hydrocarbons, benzene and pyrenes (Table 1, Figure
flocculating (Taha et al., 2012; Larson, 1998). 1). The MO is a standard electrical transformer oil,
In addition to the nanoparticle size, geometry and Nytro 10XN, from Nynas (2016). The present study
other liquid layering and Brownian motion aspects of aims at the prediction of the viscosity performance of
the reinforcements for HTFs, there is an important phe- h-BN/MO nanofluid using MD simulations (Sakha-
nomenon: the thermophoresis effect, which is impor- vand and Shahsavari, 2014) and comparison with ex-
tant in heat transport. This effect includes the increase perimental work.
in effective viscosity of nanofluids due to the presence
of nanoparticles and fluid density variation due to va- Materials and methods
riable volume fraction. Kandasamy et al. (2013) obser-
ved that the combined effect of thermophoresis and In our work we performed a liquid exfoliation of micro-
Brownian motion play a very dominant role on heat meter-sized layered h-BN crystals (Sigma Aldrich),
transfer in the presence of thermal stratification, mainly where we obtained 2D h-BN nanosheets (ranging from
due to the size of the nanoparticles. 100 to 500 nm). This process allows us to obtain thin
Lubricants can be used to minimize contact friction layer nanostructures containing a few atomic layers (~5
between components, resulting into considerable layers) that can be stabilized within MO via molecular
energy and tooling savings (Stachowiak and Batche- interactions as well as Brownian motion mechanism
lor, 2005). However, frictional heat generated, when (Taha et al., 2012). MD simulations are used to calculate
two or more moving surfaces are in contact, could the temperature dependency of viscosity and tribologi-
reach high temperatures and needs to be drawn out cal performance of MO at various filler fraction of 2D
from the system in an efficient manner, since it can hexagonal boron nitride (h-BN) nanosheets. Experi-
also degrade lubricants or oxidize them; thus, good mental part of the investigation was previously perfor-
lubricants must possess adequate thermal conductivi- med by Taha et al. (2012).
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM 103
Theoretical Prediction of Physical Properties (Viscosity) on 2D-based Nanofluids
Force fields and McDonald, 1990; Kang and Hwang, 2004). Len-
nard-Jones (LJ) potential accounts for the cross inte-
The consistent valance force field (CVFF) (Rafiee et al., ractions between h-BN atoms, and MO atoms. The
2013; Albe and Möller, 1998) is used for the interaction necessary LJ parameters are calculated using standard
of the MO molecules. CVFF has been successfully Lorentz-Berthelot combination rules (Plimpton, 1995)
used to simulate hydrocarbon molecules of crude oil from LJ parameters of MO determined by CVFF and
(Sekkal et al., 1998). The in-plane interactions of B and LJ parameters of B and N (Taha et al., 2012). This com-
N atoms in h-BN are determined via a bond-order bination has been used in h-BN and crude oil simula-
Tersoff-like potential parameterized for BN (Moon tions (Sekkal et al., 1998). A scalable parallel molecular
and Hwang, 2004), previously used to simulate BN dynamics code, LAMMPS (Plimpton, 1995) was used
sheets (Vaccarini et al., 2000; Verma et al., 2007; Hansen for the molecular dynamics simulations.
104 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM
Taha-Tijerina José, Sakhavand Navid, Kochandra Raji, M. Ajayan Pulickel, Shahsavari Rouzbeh
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM 105
Theoretical Prediction of Physical Properties (Viscosity) on 2D-based Nanofluids
Acknowledgements
Authors acknowledge the support from CONACYT
(213780). The financial support by Rice University is
greatly acknowledged. The supercomputer machines
utilized in this work were supported in part by NIH
award NCRR S10RR02950 and an IBM Shared Univer-
sity Research (SUR) Award in partnership with CISCO,
Qlogic and Adaptive Computing, and in part by the
Data Analysis and Visualization Cyber infrastructure
funded by NSF under grant OCI-0959097.
106 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM
Taha-Tijerina José, Sakhavand Navid, Kochandra Raji, M. Ajayan Pulickel, Shahsavari Rouzbeh
References Haubner R., Wilhelm M., Weissenbacher R., Lux B. Boron Nitri-
des-Properties, Synthesis and Applications. In: Jansen PDM,
Albe K., Möller W. Modelling of boron nitride: Atomic scale simu- Editor. High Performance Non-Oxide Ceramics II. Berlin,
lations on thin film growth. Computational Materials Science, Heidelberg: Springer Berlin Heidelberg, 2002, pp. 1-45.
volume 10 (issues 1-4), 1998; 111-115. Hernández-Battez A., González R., Viesca J.L., Fernández J.E.,
Baby T.T., Sundara R. Synthesis and Transport Properties of Metal Díaz Fernández J.M., Machado A., Chou R., Riba J. CuO, ZrO2
Oxide Decorated Graphene Dispersed Nanofluids. The Journal and ZnO nanoparticles as antiwear additive in oil lubricants.
of Physical Chemistry C, volume 115 (issue 17), 2011: 8527-8533. Wear, volume 265 (issues 3-4), 2008: 422-428.
Balandin A.A. Thermal Properties of Graphene, Carbon Nanotu- Hernández-Battez A., Viesca J.L., González R., Blanco D.,
bes and Nanostructured Carbon Materials. Nature Materials, Asedegbega E., Osorio A. Friction reduction properties of a
volume 10, 2011: 569-581. CuO nanolubricant used as lubricant for a NiCrBSi coating.
Botha S.S., Ndungu P., Bladergroen B.J. Physicochemical Proper- Wear, volume 268 (issues 1-2), 2010: 325-328.
ties of Oil-Based Nanofluids Containing Hybrid Structures of Husain E., Narayanan T.N., Taha-Tijerina J.J., Vinod S., Vajtai R.,
Silver Nanoparticles Supported on Silica. Industrial & Enginee- Ajayan P.M. Marine Corrosion Protective Coatings of Hexago-
ring Chemistry Research, volume 50 (issue 6), 2011: 3071-3077. nal Boron Nitride Thin Films on Stainless Steel. ACS Applied
Coleman J.N. et al. Two-dimensional nanosheets produced by li- Materials & Interfaces,volume 5 (issue 10), 2013: 4129-4135.
quid exfoliation of layered materials. Science, volume 331 (is- Jost H.P. Tribology-Origin and future. Wear, volume 136 (issue 1),
sue 6017), 2011: 568-571. 1990: 1-17.
Chang H., Li Z.Y., Kao M.J., Huang K.D., Wu H.M. Tribological Kandasamy R., Muhaimin I., Mohamad R. Thermophoresis and
property of TiO2 nanolubricant on piston and cylinder surfa- Brownian motion effects on MHD boundary-layer flow of a
ces. Journal of Alloys and Compounds, volume 495 (issue 2), nanofluid in the presence of thermal stratification due to solar
2010: 481-484. radiation. International Journal of Mechanical Sciences, volume
Choi C., Yoo H.S., Oh J.M. Preparation and heat transfer proper- 70, 2013: 146-154.
ties of nanoparticle-in-transformer oil dispersions as advan- Kang J.W., Hwang H.J. Comparison of C60 encapsulations into car-
ced energy-efficient coolants. Current Applied Physics, volume bon and boron nitride nanotubes. Journal of Physics: Condensed
8 (issue 6), 2008: 710-712. Matter, volume 16 (issue 23), 2004: 3901-3908.
Choi S.U.S. Nanofluids: A New Field of Scientific Research and Kimura Y., Wakabayashi T., Okada K., Wada T., Nishikawa H. Bo-
Innovative Applications. Heat Transfer Engineering, volume 29 ron nitride as a lubricant additive. Wear, volume 232 (issue 2),
(issue 5), 2008: 429-431. 1999: 199-206.
Cho D.H., Kim J.S., Kwon S.H., Lee C., Lee Y.Z. Evaluation of Larson R.G. The Structure and Rheology of Complex Fluids, USA,
hexagonal boron nitride nano-sheets as a lubricant additive in Oxford University Press, 1998.
water. Wear, volume 302 (issues 1-2), 2013: 981-986. Lee P.H., Nam J.S., Li C., Lee S.W. An experimental study on mi-
Eastman J.A., Choi S.U.S., Li S., Yu W., Thompson L.J. Anoma- cro-grinding process with nanofluid minimum quantity lubri-
lously increased effective thermal conductivities of ethylene cation (MQL). International Journal of Precision Engineering and
glycol-based nanofluids containing copper nanoparticles. Manufacturing, volume 13 (issue 3), 2012: 331-338.
Applied Physics Letters, volume 78 (issue 6), 2001: 718. Lin Y., Connell J.W. Advances in 2D boron nitride nanostructures:
Eichler J, Lesniak C. Boron nitride (BN) and BN composites for nanosheets, nanoribbons, nanomeshes, and hybrids with gra-
high-temperature applications. Journal of the European Ceramic phene. Nanoscale, volume 4 (issue 22), 2012: 6908.
Society, volume 28 (issue 5), 2008: 1105-1109. Lin J.C. Compression and wear behavior of composites filled with
Eswaraiah V., Sankaranarayanan V., Ramaprabhu S. Graphene- various nanoparticles. Composites Part B: Engineering, volume
based engine oil nanofluids for tribological applications. ACS 38 (issue 1), 2007: 79-85.
Applied Materials & Interfaces, volume 3 (issue 11), 2011: 4221- Mohan M., Thomas S., Taha-Tijerina J., Narayanan T.N., Sobhan C.B.,
4227. Ajayan P.M. Heat Transfer Studies in Thermally Conducting
Gao J.W., Zheng R.T., Ohtani H., Zhu D.S., Chen G. Experimental and Electrically Insulating Nano-Oils in a Natural Circulation
investigation of heat conduction mechanisms in nanofluids. Loop. In: Proceedings of ASME 2013 International Mechanical
Clue on clustering. Nano Letters, volume 9 (issue 12),2009: Engineering Congress & Exposition IMECE 2013. San Diego, Ca-
4128-4132. lifornia, USA, Volume 6B: Energy. p. V06BT07A040
Golberg D., Bando Y., Huang Y., Terao T., Mitome M., Tang C., Zhi Moon W.H., Hwang H.J. Molecular-dynamics simulation of struc-
C. Boron nitride nanotubes and nanosheets. ACS Nano, volu- ture and thermal behaviour of boron nitride nanotubes. Nano-
me 4 (issue 6),2010: 2979-2993. technology, volume 15 (issue 5), 2004: 431-434.
Hansen J.P., McDonald I.R. Theory of Simple Liquids, 2nd ed., Lon- Nguyen C.T., Desgranges F., Roy G., Galanis N., Maré T., Boucher
don, 1990. S., Angue Mintsa H. Temperature and particle-size dependent
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM 107
Theoretical Prediction of Physical Properties (Viscosity) on 2D-based Nanofluids
viscosity data for water-based nanofluids-Hysteresis pheno- phene oxide for efficient thermal management. RSC Advances,
menon. International Journal of Heat and Fluid Flow, volume 28 volume 4 (issue 47), 2014: 24887.
(issue 6), 2007: 1492-1506. Taha-Tijerina J., Peña-Paras L., Narayanan T.N., Garza L., Lapray
Nynas. NYTRO 10XN - Nynas Transformer Oil [Internet]. General C., Gonzalez J., Palacios E., Molina D., García A., Maldonado
information, 2016 [cited 2016 September]. D., Ajayan P.M. Multifunctional nanofluids with 2D nanos-
Pakdel A., Zhi C., Bando Y., Nakayama T., Golberg D. Boron nitri- heets for thermal and tribological management. Wear, volume
de nanosheet coatings with controllable water repellency. 302 (issues 1-2), 2013: 1241-1248.
ACS Nano, volume 5 (issue 8), 2011: 6507-6515. Taha-Tijerina J., Narayanan T.N., Gao G., Rohde M., Tsentalovich
Peña-Parás L., Taha-Tijerina J., Garza L., Maldonado-Cortés D., D.A., Pasquali M., Ajayan P.M. Electrically insulating thermal
Michalczewski R., Lapray C. Effect of CuO and Al2O3 nano- nano-oils using 2D fillers. ACS Nano, volume 6 (issue 2), 2012:
particle additives on the tribological behavior of fully formu- 1214-1220.
lated oils. Wear, 2015: 332-333, 1256-1261. Taha-Tijerina J.J., Narayanan T.N., Tiwary C.S., Lozano K., Chipa-
Peña-Parás L., Taha-Tijerina J., García A., Maldonado D., Gonzá- ra M., Ajayan P.M. Nanodiamond based thermal fluids. ACS
lez J.A., Molina D., Palacios E., Cantú P. Antiwear and Extre- Applied Materials & Interfaces, volume 6, 2014: 4778-4785.
me Pressure Properties of Nanofluids for Industrial Taha-tijerina J., Narayanan T.N., Avali S., Ajayan P.M. 2D Structu-
Applications. Tribology Transactions, volume 57 (issue 6), 2014: res-based Energy Management Nanofluids. In: Proceedings
1072-1076. of ASME 2012 International Mechanical Engineering Con-
Peña-Parás L., Maldonado-Cortés D., Taha-Tijerina J., García-Pi- gress & Exposition IMECE 2012, Houston, Texas, p. IMECE
neda P., Garza G.T., Irigoyen M., Gutiérrez J., Sánchez D. Ex- 2012-87890.
treme pressure properties of nanolubricants for metal-forming Taha-Tijerina J., Narayanan T.N., Gao G., Rohde M., Tsentalovich
applications. Industrial Lubrication and Tribology, volume 68 D.A., Pasquali M., Ajayan P.M. Electrically Insulating Ther-
(issue 1), 2016: 30-34. mal Nano-Oils Using 2D Fillers. ACS Nano, volume 6 (issue 2),
Plimpton S. Fast Parallel Algorithms for Short-Range Molecular 2012: 1214-1220.
Dynamics. Journal of Computational Physics, volume 117 (issue Timofeeva E.V., Gavrilov A.N., McCloskey J.M., Tolmachev Y.V.,
1), 1995: 1-19. Sprunt S., Lopatina L.M., Selinger J.V. Thermal conductivity
Raji K., Sobhan C.B., Taha-Tijerina J., Narayanan T.N., Ajayan and particle agglomeration in alumina nanofluids: experi-
P.M. Molecular dynamic simulation of thermal conductivity ment and theory. Physical Review E, Statistical, Nonlinear, and
of electrically insulating thermal nano-oil. In: Proceedings of Soft Matter Physics, volume 76 (issue 6 Pt 1), 2007: 61203.
ASME 2012 International Mechanical Engineering Congress Timofeeva E.V., Routbort J.L., Singh D. Particle shape effects on
& Exposition IMECE, 2012. Houston, Texas, pp. 1565-1571. thermophysical properties of alumina nanofluids. Journal of
Rafiee M.A., Narayanan T.N., Hashim D.P., Sakhavand N., Shahs- Applied Physics, volume 106 (issue 1), 2009: 14304.
avari R., Vajtai R., Ajayan P.M. Hexagonal Boron Nitride and Tzanakis I., Hadfield M., Thomas B., Noya S.M., Henshaw I., Aus-
Graphite Oxide Reinforced Multifunctional Porous Cement ten S. Future perspectives on sustainable tribology. Renewable
Composites. Advanced Functional Materials, volume 23 (issue & Sustainable Energy Reviews, volume 16 (issue 6), 2012: 4126-
45), 2013: 56245630. 4140.
Sakhavand N., Muthuramalingam P., Shahsavari R. Toughness Vaccarini L., Goze C., Henrard L., Hernández E., Bernier P., Ru-
governs the rupture of the interfacial H-bond assemblies at a bio A. Mechanical and electronic properties of carbon and bo-
critical length scale in hybrid materials. Langmuir, volume 29 ron–nitride nanotubes. Carbon, volume 38 (issues 11-12), 2000:
(issue 25), 2013: 8154-8163. 1681-1690.
Stachowiak G., Batchelor A.W. Engineering Tribology, 3rd ed., Verma V., Jindal V.K., Dharamvir K. Elastic moduli of a boron ni-
Oxford, UK, Elsevier Butterworth-Heinemann, 2005. tride nanotube. Nanotechnology, volume 18 (issue 43),2007:
Sakhavand N., Shahsavari R. Synergistic Behavior of Tubes, 435711.
Junctions, and Sheets Imparts Mechano-Mutable Functionali- Wang X., Xu X., Choi S.U.S. Thermal Conductivity of Nanoparti-
ty in 3D Porous Boron Nitride Nanostructures. Journal of Phy- cle-Fluid Mixture. Journal of Thermophysics and Heat Transfer,
sical Chemistry C: Nanomaterials and Interfaces, volume 118 volume 13 (issue 4), 1999: 474-480.
(issue 39),2014: 22730-22738. Wu Y.Y., Tsui W.C., Liu T.C. Experimental analysis of tribological
Sekkal W., Bouhafs B., Aourag H., Certier M. Molecular-dynamics properties of lubricating oils with nanoparticle additives.
simulation of structural and thermodynamic properties of bo- Wear, volume 262 (issues 7-8), 2007: 819-825.
ron nitride. Journal of Physics: Condensed Matter, volume 10 (is- Xie H., Chen L. Review on the Preparation and Thermal Perfor-
sue 23), 1998: 4975-4984. mances of Carbon Nanotube Contained Nanofluids. Journal of
Sudeep P.M., Taha-Tijerina J., Ajayan P.M., Narayanan T.N., Chemical & Engineering Data, volume 56 (issue 4), 2011: 1030-
Anantharaman M.R. Nanofluids based on fluorinated gra- 1041.
108 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM
Taha-Tijerina José, Sakhavand Navid, Kochandra Raji, M. Ajayan Pulickel, Shahsavari Rouzbeh
Yi M., Shen Z., Zhao X., Liang S., Liu L. Boron nitride nanosheets Citation for this article:
as oxygen-atom corrosion protective coatings. Applied Physics
Chicago style citation
Letters, volume 104 (issue 14), 2014: 143101.
Taha-Tijerina, José, Navid Sakhavand, Raji Kochandra, Pulickel M.
Yu H., Xu Y., Shi P., Xu B., Wang X., Liu Q. Tribological properties
Ajayan, Rouzbeh Shahsavari. Theoretical Prediction of Physical Pro-
and lubricating mechanisms of Cu nanoparticles in lubricant.
perties (Viscosity) on 2D-based Nanofluids. Ingeniería Investigación y
Transactions of Nonferrous Metals Society of China, volume 18 Tecnología, XVIII, 01 (2017): 101-109.
(issue 3): 2008: 636-641.
Yuge T., Ito N., Shimizu A. Nonequilibrium Molecular Dynamics ISO 690 citation style
Simulation of Electric Conduction. Journal of the Physical Socie- Taha-Tijerina J.J., Sakhavand N., Kochandra K., Ajayan P.M., Shahsava-
ty of Japan, volume 74 (issue 7), 2005: 1895-1898. ri R. Theoretical Prediction of Physical Properties (Viscosity) on 2D-
based Nanofluids. Ingeniería Investigación y Tecnología, volume XVIII
(issue 1), January-March 2017: 101-109.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 101-109 ISSN 1405-7743 FI-UNAM 109
I ngeniería I nvestigación y T ecnología
volumen XVIII (número 1), enero-marzo 2017 111-126
ISSN 1405-7743 FI-UNAM artículo arbitrado
Información del artículo: recibido: junio de 2016 aceptado: octubre de 2016
Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license
Miranda-Aviles Raúl
Universidad de Guanajuato, México
Departamento de Minas, Metalurgia y Geología
E-mail: rmiranda@[Link]
Abstract
The atmospheric deposition of some heavy metals was investigated using saxicolous lichen species (Xanthoparmelia
mexicana (Gyeln.) Hale, Xanthoparmelia tasmanica (Hook. f. & Taylor) Hale, Caloplaca aff. brouardii ([Link].) Zahlbr,
Caloplaca aff. ludificans Arup, and Aspicilia sp), samples were collected from three zones (rural, suburban and urban)
along the Guanajuato city, during the months of October-November 2012, April, July, and October 2013 and January
2014. Lichen samples were analyzed using the Inductively Coupled Plasma Mass Spectrometry technique. The con-
centrations of heavy metals in lichen samples from the Xanthoparmelia species ranged from (96.21 μg g -1) for lead
(Pb), (95.10 μg g−1) for zinc (Zn), (58.40 μg g−1 )for vanadium (V), (105.15 μg g−1) for Chrome (Cr), and (48.93 μg g−1)
for Niquel (Ni). Caloplaca species (92.42, μg g−1) for lead (Pb), (172.97 μg g−1 ) for Zinc (Zn), (53.51 μg g−1 ) for vanadium
(V), (91.23 μg g−1 ) for copper (Cu), respectively, and Aspicilia sp (612.91μg g−1) for lead (Pb), (72.24 μg g−1 ) for zinc
(Zn), (56.25 μg g−1) for vanadium (V), (18.24 μg g−1) for copper (Cu). The statistical significance of between Co-V, Ni-
Cr, Ni-Co, Sn-Zn, Co-Cr, Zn-Th, Sn-Th and Co-Zn concentrations confirmed anthropogenic sources mainly due to
emissions from vehicular traffic, fossil fuel combustion correlations, solid waste disposal and other local anthropo-
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
genic activities. Pollution indices were additionally calculated by heavy metals concentrations in order to use lichens
in Guanajuato city as bioindicators of air pollution. The concentration of these metals was observed to be in higher
range as maximum values of Pb, Zn, V, and Cu reported from the lichen samples for the suburban and urban zones
in Guanajuato city. The accumulations of Ni and Cr from both zones are similar in concentration. The contamination
factors or the pollution index factor and the pollution load index criteria revealed high levels of Be, Cu, Co, Zn, Pb,
and Th in Caloplaca species and Aspicilia sp., while Xanthoparmelia species show higher values only in Be, Sb and Pb.
The results revealed that the most sensitive lichens were Aspicilia sp., with the highest levels of Pb. The results obtai-
ned reveal important contributions towards understanding of heavy metal deposition patterns and provide baseline
data that can be used for potential identification of areas at risk from atmospheric heavy metals contamination in the
region. The use of saxicolous lichens provide a cost–effective approach for monitoring regional atmospheric heavy
metal contamination and may be effectively used in large scale air pollution monitoring programmer.
Keywords: lichens, heavy metal pollution, indicator, Guanajuato city.
Resumen
En este estudio se analizaron la concentraciones de algunos metales pesados encontrados en la atmosfera de la ciudad de Guana-
juato, empleando para ello, especies de líquenes saxícolas (Xanthoparmelia mexicana (Gyeln.) Hale, Xanthoparmelia tasmanica
(Hook. f. & Taylor) Hale, Caloplaca aff. brouardii ([Link].) Zahlbr, Caloplaca aff. ludificans Arup, and Aspicilia sp) recolecta-
dos en tres zonas (rural, suburbana y urbana) a lo largo de los meses octubre-noviembre de 2012, abril, julio y octubre de 2013,
así como enero de 2014. Las muestras de líquenes se analizaron usando la técnica de espectrometría de masas de plasma acop-
lado inductivamente. La concentración de metales pesados en las muestras de los líquenes de la especie de Xanthoparmelia van
de (96.21 μg g−1) para plomo (Pb), (95.10 μg g−1) para zinc (Zn), (58.40 μg g−1 ) para vanadio (V), (105.15 μg g−1) para Cromo
(Cr) y (48.93 μg g−1) para Niquel (Ni). Caloplaca especie (92.42, μg g−1) para plomo (Pb), (172.97 μg g−1) para Zinc (Zn), (53.51
μg g−1) para vanadio (V), (91.23 μg g−1) para cobre (Cu), respectivamente, y Aspicilia sp (612.91μg g−1) para plomo (Pb), (72.24
μg g−1) para zinc (Zn), (56.25 μg g−1) para vanadio (V), (18.24 μg g−1) para cobre (Cu). La estadística entre las concentraciones
de Co-V, Ni-Cr, Ni-Co, Sn-Zn, Co-Cr, Zn-Th, Sn-Th and Co-Zn confirma una fuente antropogénica, principalmente debida a
las emisiones del tráfico vehicular, combustión y actividades antropogénicas locales. Los índices de contaminación se calcularon,
con la finalidad de utilizar los líquenes estudiados como indicadores de calidad del aire en la ciudad de Guanajuato. Se observaron
valores altos en Pb, Zn, V y Cu. La concentración de Ni y Cr en ambas zonas son similares. Los factores de contaminación o el
factor de índice de la contaminación y los criterios de índice de contaminación, revelaron niveles elevados de Be, Cu, Co, Zn, Pb
y Th en las especies Caloplaca y Aspicilia sp., mientras que las especies Xanthoparmelia muestran los valores más altos solo en
Be, Sb y Pb. Los resultados revelaron que el liquen más sensible fue la Aspicilia sp, con los niveles más altos de Pb. Los resultados
obtenidos revelan importantes contribuciones para la comprensión de los patrones de deposición de metales pesados y proporcio-
nan a su vez datos de referencia que pueden utilizarse para su posible identificación de zonas de riesgo. El uso de líquenes saxí-
colas ofrece un enfoque rentable para el monitoreo de la contaminación por metales pesados en la atmosfera y pueden utilizarse
con eficacia en la vigilancia de la calidad del aire.
Descriptores: líquenes, contaminación metales pesados, indicador, Ciudad de Guanajuato.
112 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
mining areas is particularly scarce in Mexico and the use Materials and methods
of lichens as bioindicators of atmospheric pollution has
received little attention. Mining activities are widely Study area
known as the main pollutant source by releasing high
concentrations of heavy metals to the environment (Co- Guanajuato is a city in central Mexico and the capital of
nesa et al., 2007). Specifically, mining activities in the the state bearing the same name. The origin and the
Guanajuato city have been developed for over 500 years, economic development of the city come from the disco-
constituting the most important economic activity in the very of mines in the surrounding mountains. The city is
region. For the last five centuries, mineral deposits have located at 2,000 meters above sea level and is characte-
been utilized to obtain commercial quantities of iron, rized by a regional semi-arid climate. Temperature va-
lead, zinc, copper, gold and silver. Iron is present in oxi- ries from 6 to 20° C in winter and 21 to 32°C in summer.
des (hematite, and magnetite) and sulfides (pyrite); lead Most of the rainfall occurs between May and Septem-
and zinc occur in galena and sphalerite, while copper is ber, with a total annual rainfall reaching 700 mm/year.
present in chalcopyrite (Puy et al., 2013). No literature on The trees, which generally reach heights of 15 m, are
heavy metal pollution in the Guanajuato city is available distributed mostly at the summit of the hills where mi-
and neither on the use of lichens as bioindicators for at- ning has not yet spread. The site was selected due to the
mospheric trace elements and heavy metal deposition. anthropogenic activities such as mining, stone crus-
The objective of this work was the use of saxicolous li- hers, and intense vehicular flow. Lichen specimens
chens as bioindicators of atmospheric deposition of trace were collected in three zones in Guanajuato city:
elements and heavy metals in the city of Guanajuato.
Samples were collected at three sites: rural, suburban 1. Urban (four sites were sampled: Plaza de Los Ange-
and urban; the first is considered a pristine area located les (14Q265491-2325649), Belen (14Q265388-
in the forest known as ‘La Bufa’ and the remaining two 2325963), Teatro Principal (14Q265949-2325548) and
are influenced by anthropogenic activities (Figure 1). Po- Universidad de Guanajuato (14Q265748-2325687).
llution indices, such as contamination factors (CFs) and 2. Suburban (two sites were sampled: Pastita
pollution load index (PLI) were used to determine the (14Q266352-2324963), and Music School
pollution state and to assess the possible sources of con- (14Q 267277-2324000), and
taminants. This study can be considered as a basis for 3. Rural (one site were sampled: la Bufa-14Q 266565-
future research on air quality monitoring in areas in- 2323525), (Figure 1).
fluenced by strong human impact (mining extractions,
industry and vehicular emissions).
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 113
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
Sampling and sample treatment calcium hypochlorite solution (C), as well as combina-
tions of both (KC).
Lichen specimens were sampled at the same time du-
ring October and November 2012, April, July, and Oc-
Geology
tober 2013 and January 2014 in three zones (Figure 1):
rural (pristine forest ‘La Bufa’), suburban (influenced The studied lichen grew in sandstones of the Losero
by mining and stone crushers), and urban (influenced Formation located in the pristine forest of ‘La Bufa’ (Fi-
by mining and intense vehicular flow). Collections of gure 2). The sandstone Losero Formation shows great
the lichens, together with their substrates, were care- diversity in grain sizes, from coarse to fine sand. The
fully sampled directly from uncontaminated rocks in grains generally are subrounded to angular and show a
the rural zone (named here standard reference mate- poor selection with respect to grain size. In regards to
rial or SRM) and the other samples were taken from composition, these sandstones are immature, with high
rocks influenced by human activities in suburban and contents of rock fragments in size ranging from some
urban zones (termed here as influenced by human ac- micrometers up to several millimeters. Optical micros-
tivities or IHA). A minimum of three samples for each copy analyses show that all the lithofacies contain
species were collected in selected sites with a steel kni- quartz, feldspars, biotite, volcanic lithics, metamorphic
fe and approximately 3-4 g of the lichen thallus were lithics and iron oxides in small quantities (Figure 3).
taken for each species. The specimens were stored in
perforated plastic bags, maintained at 4°C during the
Geochemical analysis of the sandstones of the losero
transport to the laboratory and refrigerated at -4°C be-
formation
fore use for the analysis of ICP-MS. The lichen speci-
mens were reviewed taxonomically using specialized The stone material used for the geochemical analyses is
keys (Nash et al., 2004, 2007), taking into account their crushed to a size smaller than 75u (200 mesh), and the
vegetative and reproductive characteristics (observed concentration of major elements and trace elements
in the Olympus BX41 optical microscope); and chemi- from the powdered samples were determined using
cal (lichen acids present), based on the reactions of the mass spectrometry-Inductively Coupled Plasma (ICP-
reagents: potassium hydroxide 10% (K) and saturated MS) (Table 1).
Figure 2. a, b, c) Forest La Bufa; d, e) lichen species in Losero Formation; f) Sandstone of Losero Formation
114 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
Figure 3. Mineralogical composition of sandstones of the Losero Formation, a) photomicrograph showing quartz grains consisting
mainly of polycrystalline quartz grains, b) quartz grains consisting mainly of monocrystalline grains, c) volcanic fragments, d)
monocrystalline quartz grains and plane parallel laminae, e) contact grains, f) iron oxides
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 115
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
116 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 117
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
Table 2. Chemical compositions of lichen and sandstone samples collected from rural (pristine forest “La Bufa”) in the Guanajuato city.
--- no detected
The major and trace Sandstone of the Losero Xanthoparmelia species Caloplaca species
Aspicilia sp Mean lichen
elements Formation lichen lichen
Na2O (%) --- --- --- --- ---
Al2O (%) 17.31 12.82 11.13 15.13 13.0266
SiO2 (%) 80.95 32.42 54.28 38.13 41.61
K2O (%) 8.45 13.52 5.68 5.03 8.07666
MgO (%) --- --- 1.85 1.9 1.25
CaO (%) 0.6 27.9 5.88 8.2 13.9933
Fe2O3 (%) 0.43 13.32 22.35 17.7 17.79
Be (μg g−1) 3.41 3.1021 3.1064 14.8402 7.0162
Ni( μg g−1) 7.39 48.9306 6.4322 14.0175 23.1267
Cu( μg g−1) 9.12 21.7104 18.0997 17.8997 19.2366
Co( μg g−1) 2.80 7.0012 3.9307 4.641 5.1909
Sn( μg g−1) 5.70 5.7003 27.5609 6.224 13.1617
Sb( μg g−1) 1.19 1.4497 1.6869 1.3181 1.4849
Zn( μg g−1) 37.88 95.1035 172.971 57.0977 108.3907
Pb( μg g−1) 6.14 21.4031 26.4245 20.0676 22.6317
Cr( μg g−1) 17.71 105.1565 29.3347 15.1433 49.8781
V( μg g−1) 43.55 58.4043 53.5166 56.2548 56.0585
Th( μg g−1) 0.23 5.1625 13.1387 5.1917 7.8309
Statistical analysis a wide range of variation for the first mentioned group
with an extreme value for Pb. For the second group,
Mean, standard deviation and correlation however, fluctuations in values maintain within a mi-
nor range of variation. The statistical significance of co-
The statistical analyses performed in heavy metal accu- rrelations between Co–V; Ni–Cr; Ni-Co; Sn-Zn; Co-Cr;
mulations on lichens are presented in Table 4. Mean, Zn-Th, and Sn–Th concentrations (Table 5) confirm
standard deviation, Min, Max, Med for each of the ele- anthropogenic sources mainly due to emissions from
ven metals analyzed in lichen thallus from rural, subur- vehicular traffic, fossil fuel combustion, solid waste dis-
ban and urban zones resulted significant at < 0.05. posal and other local anthropogenic activities. Subur-
According to statistical tests, the selectivity sequence of ban and urban areas show a high correlation with
metals resulted as follows: Suburban (zone 1): Xanthoparmelia-Aspicilia species (0.9276), (0.9393) res-
pectively. The correlation between these two areas is
mean very good (0.9998). While in the rural area the correla-
Pb>Zn>Cu>V>Cr>Ni>Sb>Sn>Th>Co>Be; tion is good with Caloplaca-Aspicilia sp species (Table 4).
standard deviation show the selectivity series:
Pb>Cu>Zn>V>Sb>Th>Sn>Be>Cr>Co>Ni.
Urban (Zone 2): mean
Pb>Zn>Cu>V>Cr>Ni>Sn>Th>Be>Co>Sb;
standard deviation, Pb>Cu>Zn>V>Sn>Th>Co>Be>Cr>
Ni>Sb.
Rural (Zone 3): mean Cluster analysis
Zn>V>Cr>Ni>Pb>Cu>Sn>Th>Be>Co>Sb;
standard deviation, Zn>Cr>Ni>Sn>Be>Th>Pb>V>Cu>C Cluster analysis in R-mode was performed on the li-
o>Sb. chen data set for all the examined elements. Three clus-
ters were revealed and are shown in Figures 4 and 5.
In summary, the group of metals conforming by Pb, Zn, The lichen sample site shows three statistically signifi-
Cr, V and Cu present the highest values for means, whi- cant clusters. Group 1 consists of the rural zone (sam-
le the group of elements constituted of Co, Sn, Be, Th ples 8 and 7) which is not influenced by vehicular
and Sb show the lowest values. Standard deviation has emissions activity. Group 2 consists of suburban and
118 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 119
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
Figure 5. A dendrogram showing clustering of analyzed lichen sample sites from the study area
urban zones (samples 2 and 5), and group 3 consists of ban and urban zones and it is possible that part of the
suburban, urban, and rural zones (samples 9, 4, 6, 3, 1). released Zn in this area may be derived from mining
These results confirmed that Zn, Pb, and Cu derive and vehicular emissions.
from vehicular traffic, while Co, and V have a main
common source in the metallurgical industry emis- Estimating Pollution impacts
sions. The groups of metals in clusters 1 and 2 might
indicate a geological origin, taking into account that The levels of heavy metals in lichens
elevated levels of these heavy metals were found in
pristine area lichen samples (Xanthoparmelia and Calo- In general, high levels of the heavy metals were obser-
placa species). Conversely, cluster 3 suggests that both ved in all studied lichen samples (Tables 2 and 3). These
metals could derive from a common cultural source, lichens have the ability to accumulate Pb, Cu, Co, Zn,
contributing to the high pollutant contents on lichens and V. The Xanthoparmelia species has the ability to ac-
(Caloplaca species and Aspicilia sp.). Hawksworth et al. cumulate Pb, Zn, and Cu, whilst Caloplaca species and
(2005) and Conti et al. (2008) reported that fossil fuel Aspicilia sp., have accumulated Pb, Cu, Co, Zn, and V,
combustion was an important anthropogenic source for respectively. It is well known that the deposition of Pb
Pb and Cu. In the case of Cu, important pollution sour- is mainly incorporated from street dust, by vehicular
ces include industrial emissions and fossil fuel combus- emissions and fuel combustion from circulating trucks
tion processes (Bernasconi et al., 2000; Cuni et al., 2004). and cars (Hawksworth et al., 2005; Conti et al., 2008).
Particularly, Guanajuato is a mining district well known The Zn and Cu in Caloplaca species indicate that lichen
all over the world for strip mining different types of concentrated these elements. Zn and Cu concentrations
metals, especially Ag and Au associated to sulphide are due to industrial emissions of steel or regular wear
veins. Currently, there is no record of the emissions re- of engines of automobiles, abrasion and burning of tires
leased to the environment by the mineral industry, but (Bernasconi et al., 2000). The Co and V concentrations in
it is known that the particular amount of Pb per year is the Aspicilia sp. lichen show that the above element co-
6.4 μg g−1. This activity, combined with an increase in mes from metallurgical industry emissions (Bernasconi
human population and the related heavy vehicular tra- et al. 2000). According to Bennett & Wetmore (1999), the
ffic could act as the primary sources of Pb and Cu emis- corticolous species would accumulate a higher concen-
sions. Although Zn is a metal commonly considered to tration of atmospheric elements because they are more
be emitted from anthropogenic sources (traffic, meta- exposed to a mixed atmosphere than the saxicolous
llurgy, waste incineration), this study shows much hig- species. However, higher concentrations of heavy me-
her concentrations in pristine zones. However, Aspicilia tals were found in saxicolous species (Xanthoparmelia,
sp. accumulates this metal in higher contents in subur- Caloplaca and Aspicilia sp).
120 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
Table 4. Statistical summary of Heavy Metals concentrations (μg g−1) in lichens collected at 7 sampling sites. SD, standard deviation;
Min, minimum; Max, maximum; Med, median; N, number of samples over detection limits. Mean, standard deviations, minimums,
maximums and medians were calculated using zero for under detection limit values
RURAL Caloplaca
Xanthoparmelia species Aspicilia. sp Mean SD Min Max Median N
Zone 3 species
Be 3.1021 3.1064 14.8402 7.0162 5.5323 3.1021 14.8402 3.1064 7
Ni 48.9306 6.4322 14.0175 23.1267 18.5069 6.4322 48.9306 14.0175 7
Cu 21.7104 18.0997 17.8997 19.2366 1.7511 17.8997 21.7104 18.0997 7
Co 7.0012 3.9307 4.641 5.1909 1.3124 3.9307 7.0012 4.641 7
Sn 5.7003 27.5609 6.224 13.1617 10.1839 5.7003 27.5609 6.224 7
Sb 1.4497 1.6869 1.3181 1.4849 0.1526 1.3181 1.6869 1.4497 7
Zn 95.1035 172.971 57.0977 108.3907 48.2290 57.0977 172.971 95.1035 7
Pb 21.4031 26.4245 20.0676 22.6317 2.7367 20.0676 26.4245 21.4031 7
Cr 105.1565 29.3347 15.1433 49.8781 39.5147 15.1433 105.1565 29.3347 7
V 58.4043 53.5166 56.2548 56.0585 2.0002 53.5166 58.4043 56.2548 7
Th 5.1625 13.1387 5.1917 7.8309 3.7531 5.1625 13.1387 5.1917 7
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 121
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
Table 5. Values of correlation coefficient (P<0.01) between sampling sites (suburban, urban and Rural) and the amounts of heavy
metals found in saxicolous lichen
Be Ni Cu Co Sn Sb Zn Pb Cr V Th
Be 1 0.1247 0.0074 0.3664 0.0498 -0.3069 -0.1479 -0.4139 -0.0315 0.5476 0.1754
Ni 1 -0.1530 0.7577 -0.0150 -0.1267 0.1519 -0.3536 0.9701 0.6661 0.0931
Cu 1 0.3784 0.0241 -0.1653 -0.0220 -0.2203 -0.1925 -0.0924 0.3490
Co 1 0.3167 -0.3181 0.3778 -0.5341 0.7196 0.8191 0.5831
Sn 1 -0.1829 0.9463 -0.3429 0.1494 0.5761 0.9170
Sb 1 -0.2410 -0.1977 -0.1030 -0.2803 -0.3193
Zn 1 -0.1874 0.3346 0.5732 0.8745
Pb 1 -0.3611 -0.5375 -0.4067
Cr 1 0.6646 0.2021
V 1 0.6603
Th 1
122 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
PLI The PLI was calculated for the two areas under in-
vestigation using the eleven investigated metals (Be,
We calculated the Pollution loading Index (PLI) using Ni, Cu, Co, Sn, Sb, Zn, Pb, Cr, V and Th). It was obser-
the following equation ved that the highest PLI was found at Suburban area
(0.037), while the lowest was calculated for Urban
PLI= (CF1 x CF2 x CF3 x CF4 x CF5 x ……..CFn) 1/n (0.024), the calculated PLI were found in the following
sequences: Suburban >Urban (Table 6).
where
Table 6. Lichens and their pollution index factors (PIF): contamination factors (CFs) and pollution load indices (PLI)
URBAN Xanthoparmelia Caloplaca Aspicilia
CFs PLI CFs PLI CFs PLI
(IHA) species species sp
Be 4.241 1.3671 O.0000 4.1234 1.3273 0.0196 1.3245 0.0892 0.0053
Ni 4.1254 0.0843 5.6214 0.8739 3.6721 0.2619
Cu 15.7321 0.7246 91.2346 5.0406 17.9876 1.0049
Co 1.324 0.1891 4.6389 1.180 1.8234 0.3928
Sn 1.4672 0.2573 8.3297 0.3022 2.9314 0.4709
Sb 0.7271 0.5015 1.3425 0.7958 0.8976 0.6809
Zn 45.624 0.4797 77.0293 0.4453 71.2394 1.2476
Pb 96.21 4.4957 92.234 3.4904 600.9243 29.945
Cr 8.9774 0.0853 8.6234 0.2939 5.9874 0.3953
V 12.9876 0.2223 27.4325 0.5125 17.5289 0.3115
Th 0.5432 0.1052 7.3214 0.5572 2.431 0.4682
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 123
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
124 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
Puy-Alquiza María Jesús, Miranda-Aviles Raúl, Zanor Gabriela Ana, Salazar-Hernández Ma. Mercedes, Ordaz-Zubia VeliaYolanda
Province of Tierra del Fuego (southern Patagonia, Argentina). Kinalioglu K., Ozbucak T.B., Kutbay H.G., Huseyinova R., Bilgin
Ecoxitology and Environmental Safety, volume 72, 2008: 1082- A., Demirayak A. Biomonitoring of trace elements with li-
1089. chens in Samsun City, Turkey. Ekoloji, volume 19, 2010:50-53.
Farner A.M., Bates J.W., Bell J.N.B. Ecophysiological effects of Leygonie R. Heavy metals in the atmosphere: sources, emissions,
acid rain on brypphytes and lichens. In: Bryophytes and lichens measurement, effects. Pollution Atmospherique, volume 35,
in a changing environment. Bates J.W. and A.M. Farmer eds. 1993: 63-74.
Clarendon Press, Oxford, 1992, pp. 284-313. Mori L., Gomez-Tuena A., Schaaf P., Goldstein S.L.,Perez-Arvizu
Gartner Lee Limited CE, Jones and Associates, Gamberg Consul- O., Solis-Pichardo G. Lithospheric removal as a trigger for
ting, Geological Survey of Canada, and RWDI Air Inc. Anvil flood basalt magmatism in the Trans-Mexican Volcanic Belt.
Range Mine Complex- Terrestrial effects study: Investigation Journal of Petrology, volume 50 (issue 11), 2009: 2157-2186.
into metal concentrations in vegetation, wildlife and soils. Nash T.H. III and Gries C. Lichens as bio-indicators of sulfur dio-
Draft for Discussion, prepared for Deloitte and Touche Inc. on xide. Symbiosis, volume 33, 2002: 1-21.
behalf of the Faro Mine Closure Planning Office, 2006. Nash T.H. III, Ryan B., Diederich P., Gries C., Bungarts F. Lichen
Garty J. Lichens as biomonitors for heavy metal pollution. In: flora of the greater sonoran desert region, vol. II, Lichens Un-
Markert, B. (Ed.), Plants as Biomonitors: Indicators for Heavy Me- limited, Arizona State University, USA, 2004, pp.742.
tals in the Terrestrial Environment, VCH, Weinheim, 1993, pp. Nash T.H. III, Ryan B.D., Diederich P., Gries C., Bungarts F. Lichen
193-263. flora of the greater sonoran desert region, vol. III, Lichens Unli-
Garty J. Biomonitoring atmospheric heavy metals with lichens: mited, Arizona State University, USA, 2007, pp. 567.
theory and application. Critical Reviews in Plant Sciences, volu- Nieboer E., Richardson D.H.S., Tomassini F.D. Mineral uptake
me 20 (issue 4), 2001: 309-371. and release by lichens: an overview. The Bryologist volume 81
Gries C. Lichens as indicators of air pollution. In: Nash III, T.H. (issue 2), 1978: 226-246.
(Ed.), Lichen Biology. Cambridge University Press, Cambridge, Nimis P.L., Scheidegger C., Wolseley P.A. Monitoring with li-
1996, pp. 240-254. chens-monitoring lichens. NATO Science Series IV, Earth and
Harikumar P.S., and Jisha T.S. Distribution pattern of trace metal Environmental Sciences, vol. 7, Kluwer Academic Publishers,
pollutants in the sediments of an urban wetlands in the 2002, pp. 408.
southwest coast of India. International Journal of Engineering, Nyarko B.J.B., Serfor-Armah Y., Akaho E.H.K., Adomako D., Osae
volume 2 (issue 5), 2010: 840-850. S. Determination of heavy metal pollution levels in lichens at
Harikumar P.S., Nasir U.P., Mujeebu Rahma M.P. Distribution of Obuasi gold mining area in Ghana. Journal Application Science
heavy metals in the core sediments of a tropical wetland sys- Technology, volume 9 (issues 1 , 2), 2004: 28-33.
tem. International Journal. Environmental Science Technology, Puy-Alquiza M.J., Miranda-Avilés R., Salazar-Hernández C., Ve-
volumen 6 (issue 2), 2009: 225-232. ga-González M., Cervantes J. Characterization petrophysical
Hawksworth D., Iturriaga T., Crespo A. Liquenes como bioindica- of the losero formation in the historical architecture of the
dores inmediatos de contaminación y cambios medio-am- Guanajuato city, Mexico. Revista Ingeniería Investigación y tec-
bientales en los trópicos. Revista Iberoamericana Micología, nología, volume XIV (issue 2), 2013: 191-202.
volume 22, 2005: 71-82. Tomlinson D.L., Wilson J.G., Harris C.R., Jeffrey D.W. Problems in
Jeran Z., Jac´imovic R., Batic F., Mavsar R. Lichens as integrating the assessments of heavy-metal levels in estuaries and forma-
air pollution monitors. Environmental Pollution, volume 120, tion of a pollution index. Helgol Meeresunters, volume 33, 1980:
2002: 107-113. 566-575.
I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM 125
Study of the Distribution of Heavy Metals in the Atmosphere of the Guanajuato City: Use of Saxicolous Lichen Species as Bioindicators
126 I ngeniería I nvestigación y T ecnología , volumen XVIII (número 1), enero-marzo 2017: 111-126 ISSN 1405-7743 FI-UNAM
La Revista Ingeniería, Investigación y Tecnología
es impresa por el Departamento de Publicaciones
de la Facultad de Ingeniería de la UNAM.
Derechos Reservados ©
Ingeniería, Investigación y Tecnología is an arbitrated journal About the format of the manuscript
for scientific and technological research, edited and printed
by the Engineering Faculty of UNAM. Its purpose is to • Written in Microsoft Word.
spread original and unpublished research and information • Only one column in Times 11, with multiple line spacing
related to all areas of engineering. It is a forum for profes- set at 1.2, and margins provided as normal in the Word
sors, researchers, professionals and students of the world, Layouts.
which includes articles both in English and Spanish. • Figures, photographs and tables should be inserted in
This guide is an extract of the document available in our their corresponding place throughout the article, not at
website for elaborating and submitting manuscripts, and the end.
can be modified without previous notice; therefore it is the • The figures and photos used lose resolution when in-
responsibility of the authors, to be sure they have consulted serted in Word and are accepted only for evaluation pur-
the last version. It is recommended that you open and read poses. The authors must have TIFF, JPG or EPS versions
it thoroughly before you submit your work, and that you with at least 300 dpi quality (this is printing resolution,
verify the date of the last version, which appears at the be- not screen quality), ready to be sent in case the paper is
gining of the document. accepted for publication. Color figures and photographs
may be included, and will appear in color in the digital
About the papers version of the Journal.
• The article must have a maximum of fifteen pages or up
Only articles related to areas of engineering will be consid- to 6,000 words (whichever is reached first) including ta-
ered for evaluation and its possible publication. bles and figures and/or photographs. The Editor reserves
The following types of papers will not be accepted: the right to ask that long articles be shortened. In excep-
tional cases the Editor can agree to publish longer articles.
• General information articles. • The mathematical expressions must be written clearly
• Historic reviews, even if their topic is focused on some and the International System of Units must be used.
engineering problem.
• Studies of the state of the art (or reviews) in engineering
that do not contribute a solution or an innovation alter-
About the contents
native to the problem in question.
The papers must have the following, not necessarily exhaus-
• Papers whose objective is essentially to promote or at- tive, structure.
tack persons, work teams, products and/or services, or
third parties in general. • Title, complete identification of authors (name, place of
work with complete names, e-mail, postal address), ab-
stract and keywords in English and Spanish.
About submitting
• Introduction, methods, results, discussion, conclusion, ac-
The articles will be received only by e-mail at the flowing knowledgments, appendices (if necessary), and references.
address: [Link]@[Link] • About the authors (breif resume included at the end of
the article).
Together with the article, authors are required to send
the following documents: Statement of Agreement (ethics
code), Article Registration Form (MS Excel sheet filled), and It is the responsibility of the authors:
Knowledge Classifcation Code. All three documents are
provided in the website, and constitute the first of many for- • To have the necessary permissions for publication from
mal communications between the authors and the Editor the institution that financed all or part of the work.
about the paper and they allow you to request the consider- • To timely obtain the necessary reproduction permissions
ation of your work for evaluation. Go to our website and see for the material used from other publications, as well as
the Instructions for Authors mentioned above for the com- to include the appropriate references (our website con-
plete description of all aspects it must include. tains a guide for citing).
At the time the paper is received, an e-mail to the main Once more you are reminded that this is an extract of the
author will be sent with the corresponding Acknowledge- complete version of the document Editorial Policy (instruc-
ment of Receipt as an attachment. In case of a delay in re- tions for authors) available in our website that must be thor-
ceiving such e-mail, you should ask for it at the same oughly read to comply with all the requirements. Authors
address: [Link]@[Link] must also read and comply with our Ethics Code.
Política Editorial / Instrucciones para autores
[Link]
Ingeniería, Investigación y Tecnología es una publicación arbitra- Sobre el formato del manuscrito
da, de investigación científica y tecnológica editada e impresa
por la Facultad de Ingeniería de la UNAM. Su propósito es di- • Escrito en MS Word.
fundir información y conocimientos originales e inéditos rela- • Una sola columna, fuente: Times 11, con interlineado
cionados con la ingeniería en todas sus áreas. Es un foro para múltiple de 1.2 y márgenes tipo normal que poporciona
profesores, investigadores, profesionales y estudiantes de todo el Diseño de página de Word.
el mundo, el cual incluye artículos en español y en inglés. • Las figuras, fotografías y tablas deben ir colocadas en el lu-
La presente guía es un extracto del documento disponible gar que les corresponde a lo largo del artículo, no al final.
en nuestro sitio web, para la elaboración y envío de manuscri- • Las figuras y fotografías utilizadas, al insertarlas en Word
tos y puede ser modificado sin previo aviso, por lo que es res- pierden resolución y se admiten sólo para la evaluación.
ponsabilidad de los autores, asegurarse de haber consultado la Los autores deben contar con las versiones TIFF, JPG o
última versión. Se recomienda leerlo completo en nuestro sitio EPS con calidad mínima de 300 dpi (resolución para im-
web antes de enviar su trabajo y verificar la fecha de la última presión, no pixeles para pantalla), para enviarlas en caso
modificación, la cual aparece al principio de ese documento. de que su artículo sea aceptado para publicación. Pueden
incluirse figuras y fotos a color, las cuales aparecerán así
Sobre el tipo de artículos en la versión digital.
• El artículo debe constar máximo de quince páginas o
Serán considerados para evaluación y posible publicación hasta 6,000 palabras (lo que ocurra primero) incluyendo
únicamente los artículos relacionados con cualquiera de las tablas y figuras. El Editor puede pedir que artículos más
ingenierías. largos se recorten. En casos excepcionales el Editor acep-
tará artículos más largos.
No serán considerados los siguientes tipos de trabajos: • Las expresiones matemáticas deben estar escritas clara-
mente y se debe utilizar el Sistema Internacional de
• Artículos de divulgación. Unidades.
• Reseñas históricas aún cuando la reseña trate sobre al-
gún problema de la ingeniería.
Sobre el contenido del manuscrito
• Estudios del estado del arte de la ingeniería que no apor-
ten alguna solución o alternativa de innovación a lo ana- Los artículos deberán tener la siguiente estructura (no nece-
lizado. sariamente exhaustiva).
• Artículos cuya finalidad esencial sea promover o agredir
a particulares, grupos de trabajo, productos y/o servicios
• Título, identificación completa de cada autor (nombre,
lugar de trabajo con nombres completos, correo electró-
o terceros en general.
nico, dirección postal), resumen y palabras clave en es-
pañol e inglés.
Sobre el envío • Introducción, desarrollo, discusión o resultados y análisis,
conclusiones, agradecimientos, apéndices (en su caso) y
Los artículos serán recibidos sólo por correo electrónico en referencias.
la siguiente dirección: [Link]@[Link] • Semblanza de los autores (resumen curricular breve in-
Con el artículo deberá incluir como archivos anexos, los cluido al final del artículo).
siguientes documentos: Declaratoria de conformidad (códi-
go de ética), Formato de registro de artículos (hoja de MS Es responsabilidad de los autores:
Excel debidamente llenada), y Código de clasifcación del
• Contar con los permisos necesarios para su publicación
conocimiento. Los tres documantos están disponibles en el
por parte de la institución que financió la investigación.
sitio web y constituyen la primera comunicación formal de los
• Obtener oportunamente los permisos de reproducción
autores con el Editor acerca del artículo que se envía y les
parcial de material que se haya utilizado de otras publi-
permite solicitar la consideración del trabajo para su evalua-
caciones, así como citar correctamente su procedencia.
ción. Vaya a la página mencionada y consulte las instruccio-
nes para autores y lea la descripción completa de los puntos Se les recuerda una vez más que este es un extracto del docu-
que debe incluir. mento de la Política Editorial (instrucciones para autores) dis-
Al recibir el artículo, enviaremos por correo electrónico ponible en nuestra página de internet y que deben consultar en
al autor principal, el Acuse de Recibo correspondiente. En la versión completa todos los puntos de esas instrucciones para
caso de no recibirlo a la brevedad, deberá solicitarlo por la cumplir los requerimientos. Los autores también deben leer y
misma vía en la dirección: [Link]@[Link] cumplir nuestro Código de Ética.