6533 PDF
6533 PDF
6533 PDF
INGENIERÍA MECÁNICA
Para conocer más información sobre la colección, los tí- congresos nacionales e internacionales, con la puesta en
tulos que la componen y cómo adquirirlos puede visitar marcha de proyectos de innovación educativa en su área.
la web enlace a la página de la colección en
www.lalibreria.upv.es
Y DE MATERIALES
Bombas operando como turbinas (PAT)
Principios de funcionamiento y selección
Editores
Modesto Pérez-Sánchez
Francisco Javier Sánchez-Romero
Helena M. Ramos
P. Amparo López-Jiménez
Colección Manual de Referencia
Los contenidos de esta publicación han sido evaluados mediante el sistema doble ciego,
siguiendo el procedimiento que se recoge en http://bit.ly/Evaluacion_Obras
Para referenciar esta publicación utilice la siguiente cita: Pérez-Sánchez, M.; Sánchez-
Romero, F. J.; Ramos, H. M.; López-Jiménez, P. A. (2020). Bombas operando como
Turbinas (PAT): principios de funcionamiento y selección. Valencia: Editorial Universitat
Politècnica de València
Editores:
Modesto Pérez-Sánchez
Francisco Javier Sánchez-Romero
Helena M. Ramos
P. Amparo López-Jiménez
Si el lector detecta algún error en el libro o bien quiere contactar con los autores,
puede enviar un correo a [email protected]
MODESTO PÉREZ-SÁNCHEZ
Doctor en Ingeniería Hidráulica y Medio Ambiente. Actualmente, es profesor en el
Departamento de Ingeniería Hidráulica y Medio Ambiente de la Universitat
Politècnica de València, donde desarrolla su labor investigadora. Cuenta con
numerosas participaciones en artículos indexados, así como publicaciones en
congresos nacionales e internacionales, tanto en docencia como en investigación,
en el área de ingeniería hidráulica.
FRANCISCO-JAVIER SÁNCHEZ-ROMERO
Doctor ingeniero agrónomo por la Universitat Politècnica de València (UPV).
Profesor en el Departamento de Ingeniería Rural y Agroalimentaria de la UPV,
dentro del área de construcción, donde desarrolla su labor investigadora. Cuenta
con varios artículos indexados, así como participaciones en congresos nacionales e
internacionales. Ha participado en numerosos contratos con la Administración y
entidades privadas en investigación, análisis, proyecto y ejecución de
infraestructuras hidráulicas. Su actividad docente se ha centrado principalmente en
materias relacionadas con la construcción e ingeniería rural.
HELENA M. RAMOS
Profesora en el Instituto Superior Técnico, Universidad de Lisboa, desarrollando su
actividad investigadora en el grupo CERIS dentro del Departamento de Ingeniería
Civil. Ha participado en numerosos proyectos de investigación nacionales e
internacionales y es coautora de gran número de artículos en revistas internacionales,
capítulos de libros y libros. También ha recibido varios reconocimientos científicos
internacionales. Es experta en transitorios, energía hidroeléctrica, nexo agua-energía,
detección de fugas, hidrodinámica y hemodinámica.
P. AMPARO LÓPEZ-JIMÉNEZ
Doctora en Ingeniería Hidráulica y Medio Ambiente. Catedrática de universidad en el
Departamento de Ingeniería Hidráulica y Medio Ambiente de la Universitat Politècnica
de València. Cuenta con numerosas participaciones en artículos indexados, así como
congresos nacionales e internacionales en el ámbito de la investigación en aspectos
de modelación de mecánica y transporte de fluidos y maquinaria hidráulica. En su
faceta docente, tiene diferentes publicaciones indexadas en revistas internacionales
relacionadas con el área de ingeniería hidráulica y mecánica de fluidos, así como
publicaciones en congresos nacionales e internacionales, con la puesta en marcha de
proyectos de innovación educativa en su área.
Resumen
El presente libro tiene por objetivo introducir al lector en el conocimiento de las bombas
operando como turbinas (PAT por su acrónimo en inglés Pump working As Turbines).
En la bibliografía actual, existen multitud de publicaciones referentes a este modo de
operación y la mejora respecto a la sostenibilidad que introducen en los sistemas de
distribución, incrementando la eficiencia energética y permitiendo la generación de
energía renovable, aprovechando los saltos hidráulicos o excesos de presión de los
sistemas. El documento no solo recoge parte de las investigaciones más importantes
desarrolladas respecto a este tema, sino que desarrolla metodologías y casos de estu-
dio que permiten al lector poder seleccionar las máquinas, así como desarrollar estu-
dios previos de recuperación energética en sistemas de distribución. En el primer capí-
tulo hace un análisis de la conveniencia de la utilización de las bombas en modo
turbinas, la motivación de su uso y la historia de su aplicación a partir de la literatura
especializada. El capítulo se centra en los diversos aspectos energéticos relacionado
con su uso en las redes de distribución de agua a presión. En el segundo capítulo se
presenta la tipología posible de la maquinaria susceptible de ser seleccionada así co-
mo la morfología propia de los tipos de máquinas más comunes.
Abordados los dos primeros capítulos donde se profundiza, define y se realiza una
revisión pormenorizada de la literatura, en el tercer capítulo plantea los principios fun-
damentales de funcionamiento de las máquinas hidráulicas, abordando desde la Ecua-
ción de Euler hasta los números adimensionales. Posteriormente estos números serán
utilizados para la determinación de curvas características. En el cuarto capítulo, se
aborda un carácter más técnico, dónde se presentan métodos empíricos que permiten
estimar el punto de operación de una bomba trabajando como turbina y viceversa.
Asimismo, se desarrolla una metodología que permite el análisis energético de una
instalación dados los puntos de operación del sistema (caudal y altura recuperable).
Esta metodología aborda la estimación empírica de puntos de operación en modo
bomba conocidos los puntos de operación como turbina, selección de la bomba a tra-
vés de catálogo comercial, desarrollo de las curvas Q-H a partir de los números adi-
mensionales y estudio de energía recuperable. Finalmente, esta metodología es apli-
cada a un caso de estudio que reproduce la metodología.
III
Bombas funcionando como Turbinas (PAT)
IV
Prólogo
El presente libro tiene por objetivo introducir a toda aquella persona interesada en el
incremento de la sostenibilidad de los sistemas de distribución de agua, concretamente
en la mejora de su eficiencia energética, mediante el empleo de bombas trabajando
como turbinas. Desde el punto de vista de los editores, la lectura y compresión del
documento puede dotar al lector de un texto resumido, descriptivo y detallado que le
permita iniciarse en el uso de las bombas operando en modo turbina en sistemas de
microgeneración. Dentro del libro, se recoge parte de los trabajos científicos publica-
dos en todo el mundo en los últimos años, dando muestra de la importancia que impli-
ca la mejora energética en redes de distribución. Esperamos que el documento ayude
a estudiantes, profesionales y científicos a proponer soluciones que mejoren la soste-
nibilidad de los sistemas de distribución, reduciendo la energía consumida y mejorando
todos los indicadores de sostenibilidad. Finalmente, los editores agradecen profunda-
mente, la participación en este documento, de los profesores Paulo Jose da Costa
Branco y João Filipe Pereira Fernandes, del Instituto Técnico Superior de Lisboa que
han desarrollado la parte eléctrica del contenido. Su aportación muestra la importancia
que existe en estos sistemas de considerar la simbiosis “hidráulica-electricidad”. Final-
mente, deseamos agradecer también a Joan M. Bengochea y Manuel Díaz-Madroñero
que desinteresadamente han contribuido al diseño de la portada.
Los editores
V
Índice
VII
La necesidad de bombas funcionando como
turbinas (PAT) en sistemas de distribución
Resumen
Las redes de distribución de agua tienen por objeto que el agua llegue al punto de con-
sumo en condiciones de presión y caudal (además de calidad) adecuadas. Esto mismo
debe conseguirse de la forma más eficiente, teniendo siempre muy presente energía y po-
tencia en los sistemas de bombeo. Los análisis de los aspectos energéticos de las redes de
distribución hacen intuir que el aprovechamiento de pequeña maquinaria hidráulica pue-
de jugar un papel importante en el conjunto de las mismas. Sin embargo, la experiencia
indica que, en la práctica, generalmente no es así: hay mucha energía que se “pierde” en
las caídas de presión provocadas por la valvulería presente en las redes (y, por otra parte,
completamente necesaria) que no se aprovecha, por regla general, en sistemas de peque-
ñas máquinas hidráulicas de recuperación.
Esta vocación de recuperación energética puede ser realizada (entre otras estrategias)
por bombas funcionando como turbinas (llamadas PAT por sus siglas en inglés: Pumps
working As Turbines) cuando se colocan en las redes para controlar la presión, al menos
en los puntos de funcionamiento en los que puedan adaptarse, en sustitución de ciertas
válvulas. La razón por la cual pueden ser más interesantes las bombas funcionando como
turbinas que las propias turbinas de pequeño tamaño son varias; principalmente de carác-
ter económico y de disponibilidad para el diseñador: suele existir mucha más oferta en
pequeñas máquinas como bombas que como pico o micro turbinas, lo que hace que las
curvas de las máquinas sean más asequibles, los precios menores y la disponibilidad en
el mercado con máquinas producidas en serie haga mucho más fácil el uso de bombas,
aunque sea con este funcionamiento inverso. Ello conllevará retos sobre la selección de
una máquina para funcionar como turbina, infiriendo sus curvas desde las de funciona-
1
Bombas funcionando como turbinas (PAT)
miento como bomba. Esta manera de operar conlleva incertidumbres propias de métodos
semiempíricos que serán presentados en capítulos posteriores.
El presente capítulo describe un análisis de la conveniencia de la utilización de las
bombas en modo turbinas, la motivación de su uso y la historia de su aplicación a partir
de la literatura especializada en los diversos aspectos energéticos relacionados con su uso
en las redes de distribución de agua a presión.
Abstract
The analysis of the energy aspects of water distribution networks suggest that the use of
small hydraulic machinery can play an important role in all of them. However, experi-
ence indicates that this is not the case: there is a lot of energy that is "lost" in the pressure
drops caused by the valves present in the networks (and, on the other hand, completely
necessary). This energy is not traditionally used in small hydropower plants, even when
there is an important amount of recovery possibilities.
This vocation of energy recovery can be carried out using pumps working as tur-
bines, when they are placed in the networks to control the pressure, at least in the operat-
ing points where they can be adapted, in substitution of certain valves. The reason why
pumps operating as turbines can be more interesting than the small turbines themselves
are several. These are mainly related to availability and economy reasons for the design-
er: there is usually much more supply in small machines such as pumps than pico or mi-
cro turbines, which makes the machine curves more affordable, lower prices and availa-
bility in the market with machines produced in series makes pumps much more
affordable, even with this reverse operation. This will entail challenges on the selection
of a machine to operate as a turbine inferring its curves from those of operation as a
pump, which entails uncertainties of semi-empirical methods that will be presented in
later chapters. This chapter makes an analysis of the suitability of the use of pumps in
turbine mode, the motivation for their use as well as the history of their application from
the specialized literature, considering the various energy aspects related to their use in
the pressurized water networks.
2
La necesidad de bombas funcionado como turbinas (PAT) en sistemas de distribución
1. Introducción
En ocasiones, dentro de la gestión de los sistemas de suministro de agua y en la
industria se requiere reducir la presión del fluido. Esta reducción de presión
generalmente se logra a través del uso de una válvula de estrangulamiento, perdiendo
así parte de la energía de presión que lleva el agua en las redes de distribución.
Teniendo en cuenta criterios de optimización energética, se hace necesario estudiar
alternativas que permitan la su recuperación en este tipo de redes.
En ese sentido, una PAT puede recuperar el potencial de la energía hidraúlica (por
poca que sea la que se está perdiendo con la ubicación de una válvula concreta) de
manera eficiente y económica con medios técnicos muy asequibles para los diseñadores.
La utilización de bombas funcionando como turbinas es aconsejable para
aplicaciones en las que se deben reducir las diferencias de presión, o donde se pueden
aprovechar la altura y el caudal de una instalación en situaciones que de otra forma se
perderían. La energía generada se puede utilizar para fines internos (máquinas aisladas)
o para alimentar a la red pública (máquinas conectadas a red). Gracias a los bajos
costos de inversión, este tipo de soluciones tienen una viabilidad económica alta, a
pesar de que los puntos de funcionamiento no siempre tienen los rendimientos óptimos
como bomba que pueden pedirse a la máquina en condiciones de diseño.
El potencial de recuperación de energía en un sistema de distribución de agua está
presente siempre que el agua pierde presión, aunque no siempre hay una máquina
perfectamente diseñada para este punto. En este sentido, las bombas son relativamente
simples y fáciles de mantener y también tienen una eficiencia máxima competitiva en
comparación con las turbinas convencionales, sobre todo cuando se considera también
otros aspectos como su disponibilidad. Quizás el mayor beneficio de las bombas
funcionando como turbinas es que la producción en serie de bombas significa que son
comparativamente, mucho más viables que las turbinas convencionales, a pesar de que
sus rendimientos máximos sean menores.
Por tanto, la labor del diseñador consiste en seleccionar bien la máquina que, a
partir de sus curvas conocidas como bomba, pueda adaptarse al punto de
funcionamiento que se espera de ella cuando funcione como turbina. En ese sentido, se
hace necesario conocer el funcionamiento de la máquina (en régimen permanente y
transitorio), las técnicas de selección, así como el análisis económico derivado de la
optimización de ubicación de la maquinaria que finalmente se seleccione.
3
Bombas funcionando como turbinas (PAT)
4
La necesidad de bombas funcionado como turbinas (PAT) en sistemas de distribución
las necesidades efectivas del sistema (McNabola et al., 2011; 2014). Los sistemas así
diseñados, deberán acoplarse a una red eléctrica síncrona, con generadores de control
de ganancia dinámico, con el fin de producir energía eléctrica sin comprometer las
condiciones de funcionamiento del sistema en términos de presión y caudal en los
puntos de consumo.
La recuperación energética a través de micromaquinaria puede hacerse en diversas
partes de un sistema hidráulico. Puede realizarse en depósitos intermedios de las redes,
en los que en ciertos momentos existan excedentes de presión para el caudal que se
requiera. En este sentido, las investigaciones descritas en Ramos et al., 2009 y Mc
Nabola et al., 2011, atestiguan que podrían recuperarse cientos de miles de euros por
año al instalar este tipo de maquinaria en ciertos casos de estudio presentados en
Portugal e Irlanda.
Puede asimismo, realizarse esta recuperación en el sistema de saneamiento; o de
control de aguas pluviales, tal como el descrito en Saket (2008), en que se utiliza una
pequeña turbina como experiencia piloto en el sistema de saneamiento de la propia
Universidad. Uno de los ejemplos más notables de ello se documenta en Sidney, donde
las instalaciones de un emisario submarino con una diferencia de cota de 60 m
permiten una generación de 12GWh al año (Eco Generation, 2008). Como ejemplo, en
Ramos et al., (2013), se propone la recuperación energética en sistemas con depósitos
de retencion en sistemas pluviales en la zona de Alcântara, en Portugal. En este caso,
los estudios teóricos que indicaban unas recuperaciones potenciales de hasta 17,5
MWh en los escenarios más optimistas en este conjunto.
Sin embargo, es en el control de la presión en la red de distribución a través de la
energía que se disipa en las válvulas de regulación, donde se centra el presente
documento. El control de la presión es necesario para la correcta operación de la red y
contiene un enorme potencial de recuperación energética. Los intentos por controlar la
presión a través de recuperación energética redundarán doblemente en la sostenibilidad
de la gestión final de la red. Las válvulas reductoras de presión son imprescindibles
para la operación de las redes. Son necesarias para mantener los niveles de presión y
caudal de suministro y necesarias hoy en día para el control de fugas. (Fontana et al.,
2012). Por tanto, en este sentido, es en el abastecimiento y fundamentalmente en el
regadío (mundialmente representa casi el 70% en consumo de agua), donde una
aplicación de sistemas de recuperación en su distribución, puede contribuir a una
mejora energética global del conjunto (Figura 1).
5
Bombas funcionando como turbinas (PAT)
6
La necesidad de bombas funcionado como turbinas (PAT) en sistemas de distribución
7
Bombas funcionando como turbinas (PAT)
8
La necesidad de bombas funcionado como turbinas (PAT) en sistemas de distribución
9
Bombas funcionando como turbinas (PAT)
artículo de Ramos y Borga (1999) describe los primeros análisis teóricos sobre las curvas
de las bombas funcionando como turbinas, con un amplio espectro de aplicaciones. El
hecho de utilizar este tipo de maquinaria en sistemas de recuperación, no solo permitirá
desarrollar investigaciones en el ámbito de la mejora de la eficiencia energética, sino que
se profundizará en el estudio de estas máquinas funcionando de forma inversa al de su
diseño, así como determinando las implicaciones con los grupos generadores y conse-
cuencias de eficiencia que se derivan de su conexión y operación.
Desde entonces, el análisis de la bibliografía muestra avances en esta línea, hacia la
recuperación de la energía disipada en las válvulas a través de micromaquinaria hidráu-
lica, particularmente en bombas operadas como turbinas. La utilización de este tipo de
máquinas es muy conveniente y bien referenciada en la bibliografía científica. Si bien
se conoce que el máximo rendimiento obtenido por una bomba funcionando en modo
turbina es menor que las turbinas convencionales, Carravetta et al., (2013) el coste
estimado es 545 €/kW. Sin embargo, una sola máquina no es capaz de adaptarse a las
condiciones de operación de las redes, por lo que las tendencias deberán ir hacia la
utilización de varias máquinas asociadas en paralelo.
Como resultado de este análisis, se enumeran a continuación las estrategias referen-
ciadas en la bibliografía (con algunos ejemplos de artículos publicados en los diversos
ítems) en relación con la recuperación energética en un sentido amplio. Se centran
sobre el control de la presión; y posteriormente, focalizado en el uso de micromaquina-
ria hidráulica para dicha recuperación, con las prestaciones particulares de las PAT.
Reducción de fugas por control de la presión, con la mejora del rendimiento hi-
dráulico de las redes (Colombo y Karney, 2002, Almandoz et al., 2005; Araujo
et al., 2006; Campisano et al., 2009; Christodoulou et al., 2013; Tucciarelli et
al., 1999; Walski et al., 2006; Fontana et al., 2012; Mutikanga et al., 2013;
Coelho et al., 2014.
Instalación de sistemas microhidráulicos de recuperación de energía (Carravetta
et al., 2014; Colombo y Kleiner, (2011); Gonçalves y Ramos, 2012; Lofman et
al., 2012; McNabola et al., 2013; Ramos et al., 2011; Simao et al., 2010; Ra-
mos, 2000; Corcoran, 2013; Pérez-Sánchez et al., 2017a,b,c,d,e,f).
Modelación del funcionamiento de bombas funcionando como turbinas para su
uso óptimo: (Nourbakhsh et al., 1992; Simao et al., 2010; Arriaga, 2010; De-
rakhshan et al., 2010; Frosina et al., 2017, Pérez-Sánchez et al. 2018 a, b).
Utilización de bombas funcionando como turbinas en sistemas de distribución
de agua. Caracterización de su funcionamiento (Ramos et al., 1999; Carravetta
et al., 2012; Carravetta et al., 2013a; Williams, et al., 1998; Rawal, S. et al.,
2007; Yang et al., 2012; Fecarotta et al., 2016.; Ramos y Simao, 2013; Nautiyal
et al., 2011; Lydon et al., 2017).
Desarrollo de sistemas de simulación para conocimiento de la energía poten-
cialmente recuperada en sistemas de distribución de agua (Adhau et al., 2012;
Sitzenfrei et al., 2014; Frijns et al., 2015; Cabrera et al., 2014; Mc Nabola,
2014; Pérez-Sánchez et al., 2016).
10
La necesidad de bombas funcionado como turbinas (PAT) en sistemas de distribución
11
Bombas funcionando como turbinas (PAT)
Conclusiones
Existen varias escalas relacionadas para la gestión de la generación de energía hidráuli-
ca, cuando se considera el agua como un recurso. Sin embargo, incluso hoy en día, la
recuperación de energía es una posibilidad muy atractiva en las redes de distribución,
con bajos costes de inversión para los gestores. El éxito de este nuevo uso depende de
la experiencia adquirida en centrales hidroeléctricas con mayor potencia instalada.
El análisis de este estado del arte descrito en este primer capítulo ha constatado que
la limitación de la utilización de PAT se centra en dos aspectos fundamentalmente: su
bajo rendimiento cuando los sistemas operan con caudales variables y la dificultad de
conjugación de las máquinas con el sistema eléctrico conjunto en la red.
Por tanto, los esfuerzos futuros deberían ser focalizados en conocer el comporta-
miento de las máquinas en condiciones de operación variables, mejorando la adaptabi-
lidad de los sistemas de las máquinas propuestas para optimizar su funcionamiento en
las condiciones variables de caudal en las líneas en que se instalan. Asimismo, deberán
ser conocidos los aspectos eléctricos de la generación de energía para tomar decisiones
sobre autoconsumo o posibles conexiones a la red eléctrica. Finalmente, deberá cono-
cerse el funcionamiento de dos máquinas en paralelo, de forma que su acoplamiento
permita trabajar con mayores rangos de caudal en sistemas más eficientes como con-
junto.
12
Saltos hidráulicos y maquinaria para la
recuperación energética
Pérez-Sánchez, Ma, Sánchez-Romero, FJb, Ramos, HMc, López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
En una red de distribución de agua existe una cierta cantidad de energía potencialmente
recuperable. En este sentido, en el capítulo anterior se ha hecho hincapié en esos aspec-
tos energéticos que impulsan a los diseñadores a determinar el uso de pequeña maquina-
ria hidráulica en lugar de válvulas reguladoras de presión. Sin embargo, todos los puntos
de las redes no son igual de idóneos para la ubicación de una máquina, principalmente
porque las máquinas hidráulicas se caracterizan por tener diferentes puntos nominales de
funcionamiento, con más o menos altura, así como caudal en su punto óptimo. La selec-
ción de la máquina adecuada para cada instalación pasa por el conocimiento de dos con-
ceptos: de un lado se hace imprescindible que el diseñador conozca los diferentes tipos
de máquinas de las que puede disponer para ser seleccionadas; de otro que el diseñador
conozca la morfología de cada tipo de máquinas, a la vez que sus curvas motrices, que
serán determinantes en la selección. El presente capítulo desarrolla ambos aspectos, pre-
sentando la tipología posible de la maquinaria susceptibles de ser seleccionadas y la mor-
fología propia de los tipos de máquinas más comunes.
Abstract
In a water distribution network, there is a potentially recoverable energy. In this sense, in
the previous chapter, emphasis has been placed on those energy aspects that drive designers
to determine the use of small hydraulic machinery instead of pressure regulating valves.
However, all the points of the networks are not equally suitable for the location of a par-
ticular machine, mainly because the hydraulic machines are characterized by having differ-
ent nominal points of operation, with more or less height and flow at their point optimum.
The selection of the right machine involves the knowledge of two concepts: on the one
13
Bombas funcionando como Turbinas (PAT)
hand it is essential that the designer knows the different types of machines that can be used
to be selected; of another that the designer knows the morphology of each type of ma-
chines, at the same time as its model curves, which will be decisive in the selection. This
chapter develops both aspects, presenting the possible typology of the machinery suscepti-
ble to be selected and the morphology of the most common types of machines.
1. Introducción
La maquinaria hidráulica es esencialmente un transformador energético de energía
mecánica en hidráulica yviceversa. En ese sentido, la utilización de pequeña
maquinaria hidráulica en las redes de abastecimiento para transformar en mecánica,
energía de presión, que de otra forma iba a disiparse en válvulas es una eficiente forma
de mejorar el comportamiento energético del sistema conjunto. De esta forma puede
proporcionarse un control de la presión, con lo que implica en la mejora asimismo en
las pérdidas por fugas (mejorando el rendimiento volumétrico de las redes) de forma
paralela.
Sin embargo, para que ello ocurra, los diseñadores tienen que encontrar máquinas
que se ajusten en su punto de funcionamiento (presión recuperada frente a caudal
circulante) a las líneas en las redes en que van a ser instaladas. Ello, en primera
instancia, pasa por conocer la tipología y morfología de las máquinas a disposición del
diseñador.
El presente capítulo se adentra en estos aspectos, con el objetivo de que el lector sea
capaz de conocer los diferentes tipos de máquinas hidráulicas susceptibles de ser
utilizadas como turbina en base a sus puntos de funcionamiento y conocer asimismo la
morfología de la máquina seleccionada.
1.1 ¿Qué es una máquina hidráulica?
Una máquina por definición es un transformador de energía. Por tanto, considerando
una definición totalmente genérica, una máquina absorbe energía y la restituye en
energía de otra clase (e.g., motor eléctrico) o de la misma clase pero transformada (e.g.,
torno) (Mataix, 1990). Dentro de las máquinas se pueden clasificar en: herramientas,
máquinas eléctricas y máquinas de fluido. Las máquinas de fluido son aquellas en las
cuales el elemento intercambiador de energía es un fluido (compresible o
incompresible). Si el fluido es compresible se denominan máquinas térmicas y en este
grupo se encuentran los motores térmicos, turbinas de vapor o turbocompresores como
ejemplo. Si la máquina opera con un fluido incompresible se denomina máquina
hidráulica. En este grupo se encuentra las bombas empleadas en cualquier sistema de
14
Saltos hidráulicos y maquinaria de recuperación energética
distribución, los ventiladores (que aunque operan con fluido compresible, se considera
un fluido incompresible en su análisis) y las turbinas existente en centrales
hidroeléctricas para el aprovechamiento de la energía potencial del agua. Por tanto, las
máquinas de fluido, o bien el fluido proporciona la energía que absorbe la máquina
(e.g., turbina) o bien aquellas en que el fluido es el receptor de energía al que la
máquina restituye la energía mecánica absorbida (e.g., bomba).
1.2 Clasificación de máquinas hidráulicas
En este aspecto, las máquinas hidráulicas pueden clasificarse atendiendo en primer lugar
a su funcionamiento. Por tanto, se habla de máquinas motoras, si las máquinas absorben
energía de un fluido y la restituye en energía de otro tipo (e.g., turbina). Por el contrario,
las máquinas generadoras son aquellas que reciben energía de un elemento exterior y
aportan parte de esa energía al fluido (e.g., bomba). Dentro de cada una de las máquinas
hidráulicas, se clasifican atendiendo al término que varían en la energía del fluido.
Atendiendo a este criterio se definen, según su clasificación general: (Figura 1).
15
Bombas funcionando como Turbinas (PAT)
2. Tipos de instalaciones
La importancia de la energía hidráulica recuperada en sistemas de distribución se basa
en el aprovechamiento de un recurso que actualmente, en la mayoría de los casos, se
está disipando bien mediante válvulas (en flujo a presión) o mediante saltos hidráulicos
(lámina libre) y no contribuye a la mejora del rendimiento energético del sistema. Ac-
tualmente, prácticamente solo son aprovechados los grandes saltos de energía.
Por tanto, es necesario (teniendo en cuenta las instalaciones y equipos disponibles
actualmente), analizar el potencial de recuperación energético en redes de distribución
tanto de lámina libre como de flujo a presión, con el objetivo de reducir la huella ener-
gética del agua en su ciclo “captación-distribución-regeneración” (la regeneración úni-
camente se aplica en el caso del abastecimiento). Además, se destaca la necesidad de
estudiar la capacidad de recuperación en la distribución de redes de riego de flujo a
presión ejecutadas durante la modernización del regadío, dado el gran volumen que se
distribuye (del orden de km3/año), para que contribuya al aumento de la eficiencia
energética de las mismas y por tanto a la sostenibilidad energética y económica de las
explotaciones agrícolas.
2.1. Tipos de instalaciones de generación
Las centrales hidroeléctricas en función de la potencia instalada pueden clasificarse en
micro (Pt<100 kW), mini (100-500 kW) y pequeñas (500-10000 kW). Aunque esta
clasificación es muy importante, ya que es la referencia para la legislación europea,
desde un punto de vista del funcionamiento en los diferentes sistemas de distribución
que se abordarán en los apartados 2.2. y .2.3 de este capítulo, es más significativo clasi-
ficarlas, en función del tipo de central (Ramos, 2000; Castro, 2006):
16
Saltos hidráulicos y maquinaria de recuperación energética
Azud o desviación del cauce
Tubería forzada
Maquinaria
de recuperación
Aliviadero de turbinas
Cauce
Maquinaria
de recuperación
Tubería forzada
Aliviadero de turbinas
17
Bombas funcionando como Turbinas (PAT)
mas inyectados para aprovechar las zonas de alta presión en las zonas bajas, co-
mo es el caso de la Figura 5.
Alturas
Pendiente hidráulica de diseño
H
Hb(c)
Altura Punto
de nudo c
crítico
crítico
cabecera F
A
Z0
E
C
bomba B D
Distancia en tubería desde la cabecera
Distancia a cabecera Di
1400
Turbina Francis
10
M
00
500
M
50
W
0
M
20
300
W
10
0
0
M
50
W
20
W
M
W
M
H (mca)
140
W
Turbina Deriaz
10
5
M
M
Turbina Flujo Cruzado
W
W
Turbina Kaplan
2
M
W
1
50
M
W
10
50
lbo
0
Bu
kW
a
50
20
bin
10
r
0
kW
Tu
0
kW
kW
5
0.5 1 2 5 10 20 50 100 200 500 1000
Q (m 3 /s)
Figura 6. Tipo de turbina en función de caudal, salto y potencia (adaptado de Castro (2006))
Además de las máquinas clásicas utilizadas, en el caso de redes de distribución o con-
ducciones en presión, con potencias inferiores a los 750 kW se pueden utilizar las PAT
que aquí se describen. Estas máquinas presentan la ventaja frente a las máquinas clásicas
de tener un coste €/kW mucho menor. Nourbakhsh et al., (1992), establecieron que el
periodo de retorno de estas instalaciones era inferior a 2 años en intervalos de potencia
entre 5-500 kW. Ya en relación con las PAT, en 1931, Thoma presentó el primer trabajo
de bombas trabajando en condiciones anormales, siendo autores posteriores: Thode
(1984), Isbăşoiu (2008), Derakhshan (2008), presentaron trabajos más detallados en la
descripción, operación y modelo teórico de estas máquinas (Arriaga, 2010).
18
Saltos hidráulicos y maquinaria de recuperación energética
19
Bombas funcionando como Turbinas (PAT)
20
Saltos hidráulicos y maquinaria de recuperación energética
21
Bombas funcionando como Turbinas (PAT)
4.00
3.50
Área regada (millones de ha) 3.00
2.50
2.00
1.50
1.00
0.50
0.00
2004 2005 2006 2007 2008 2009 2010 2011 2012 2013
Año
22
Saltos hidráulicos y maquinaria de recuperación energética
Figura 9. (a) Vista general bomba seccionada. (b) Detalle bloque hidráulico: voluta,
eje y rodete. (c) Vista frontal rodete. (d) Vista longitudinal rodete
En el caso de una bomba operando como turbina, las funciones de cada uno de sus
elementos son:
Brida de impulsión. Es la sección circular por donde se alimenta la PAT y tiene
lugar la entrada de agua.
Voluta. Tiene por función transformar parte de la energía de presión del fluido
en energía cinemática, así como distribuir el caudal uniformemente a lo largo de
la corona directriz.
Corona directriz. Parte fija de la máquina que transforma energía de presión en
energía cinemática.
Rodete. Elemento móvil que, debido a la morfología de sus álabes, transforma la
energía del fluido a través de la variación de momento cinético, transmitiendo
dicho par al eje que gira solidario al mismo.
Eje. Parte móvil rotativa que conecta la parte hidráulica y eléctrica de la máquina. Es
el encargado de acoplar con el motor eléctrico y transmitirle el par de giro al mismo.
Motor eléctrico. En este caso actúa como generador al ser conducido, generando la
energía eléctrica siempre y cuando se den las condiciones de excitación del mismo.
23
Bombas funcionando como Turbinas (PAT)
24
Saltos hidráulicos y maquinaria de recuperación energética
Conclusiones
El presente capítulo hace una definición de máquinas hidráulicas, así como una clasifi-
cación general. Esta clasificación es fundamental para la selección de maquinaria den-
tro de los diferentes sistemas de distribución (lámina libre o presurizado). Del mismo
modo, el capítulo ha mostrado los diferentes tipos de instalaciones existentes de acuer-
do al sistema de instalación del grupo de instalación y la preselección del tipo de ma-
quinaría de acuerdo a las condiciones de caudal y altura. Asimismo, este segundo capí-
tulo ha mostrado la situación actual de las diferentes instalaciones, así como el gran
potencial existente dentro de los sistemas presurizados, máxime si se tiene en cuenta el
gran volumen de agua consumida por el riego. Finalmente, el capítulo ha mostrado las
partes de una PAT y sus diferencias cuando opera como bomba y como turbina. La
descripción general desarrollada y conocida la morfología, permite que en el próximo
capítulo pueda llevarse a cabo la descripción analítica del funcionamiento hidráulico de
las bombas trabajando como turbinas.
25
Ecuaciones fundamentales para la operación de
PAT en las redes de agua a presión
Pérez-Sánchez, Ma, Sánchez-Romero, FJb, Ramos, HMc, López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
Las bombas funcionando como turbinas son, ante todo, máquinas hidráulicas. Así pues, es
necesario conocer sus curvas de funcionamiento para caracterizarlas completamente. En el
presente capítulo se describen estas curvas, a partir de sus principios fundamentales, desde
los triángulos de velocidad para las máquinas operando como bombas y como turbinas. El
funcionamiento particular de una bomba en su forma de turbina, tiene ciertos parámetros ca-
racterísticos que son presentados en el desarrollo del capítulo. Finalmente, se cierra el mismo
con un análisis de energías, fundamental para posteriores aplicaciones de las máquinas en las
redes de abastecimiento o riego.
Palabras clave: Triángulos de velocidad; Altura teórica; Altura real; Parámetros funda-
mentales de selección de PAT
Abstract
Pumps working as turbines are primarily hydraulic machines. Thus, it is necessary to know
their operating curves to characterize them completely. In the present chapter, these curves
are described, starting from their fundamental principles, from the triangles of speed for ma-
chines operating as pumps and as turbines. The particular operation of a pump in its turbine
form has certain characteristic parameters that are presented below in the development of the
chapter. Finally, it closes with an energies analysis, fundamental for subsequent applications
of the machines in the supply or irrigation networks.
27
Bombas funcionando como Turbinas (PAT)
1. Introducción y objetivos
Las características teóricas de funcionamiento de una máquina hidráulica pueden ser
deducidas de forma aproximada del comportamiento del flujo a la entrada y salida de la
misma, sin embargo, la totalidad de los fabricantes de bombas proporcional la curva de
las mismas a partir de ensayos.
El análisis de la curva de una bomba, da una indicación de sus capacidades para
proporcionar altura y caudal; y su rendimiento en función de sus parámetros de diseño,
tales como la velocidad de rotación, la geometría interna, sus pérdidas internas, etc.
El rendimiento cuando esta máquina se instale en el sistema básicamente determina-
rá la compatibilidad de la red con la bomba. De esta forma, las curvas características
determinarán si una bomba específica tiene las capacidades para proporcionar caudal y
presión para operar un sistema concreto de forma adecuada.
Estas curvas, que son fácilmente asequibles por parte del usuario para el comporta-
miento de las turbobombas, no lo son cuando se trata de conocer el funcionamiento de
la misma máquina operando como turbina. La operación como turbina es inversa a la
operación como bomba, lo que en un caso es entrada, se convierte en salida y vicever-
sa, de forma que los triángulos también tienen sentidos opuestos.
Asimismo, el concepto de rendimiento también se invierte; en el caso del funcio-
namiento como bomba el objetivo es producir energía hidráulica mientras que en el
caso del funcionamiento como turbina, se desea producir energía mecánica, con lo que
las prestaciones de una misma máquina no tienen por qué ser necesariamente iguales.
En la actualidad, no todos los fabricantes ensayan sus bombas de este modo inverso,
por lo que, para conocer el funcionamiento de las mismas, se propone aquí la utiliza-
ción de ciertos métodos aproximados basados en números adimensionales que permi-
tan la estimación de estas curvas. Estos métodos constituyen en muchas ocasiones, una
aproximación muy certera al funcionamiento de la bomba como turbina, pues están
basados en cientos de ensayos de máquinas similares y suponen para el diseñador una
herramienta potente para inferir el funcionamiento de las máquinas en este modo inver-
so que se describe en el capítulo.
Así pues, los objetivos del presente capítulo son:
Conocer el principio de funcionamiento de las turbomáquinas a partir de sus
ecuaciones fundamentales: definición de triángulos de velocidades y diagramas
de altura vs. caudal
Definir los parámetros propios y característicos de las bombas funcionando co-
mo turbinas para la determinación de sus curvas fundamentales y ser capaces de
seleccionarlas para una instalación determinada
Definir los conceptos energéticos, de rendimiento y potencias propios de la
bomba funcionando como turbina a través de los esquemas de alturas y poten-
cias propios de la máquina.
28
Selección de PAT desde la curva de bomba
2. Ecuaciones Fundamentales
2.1. Principio de funcionamiento
donde:
es la variación de la propiedad extensiva respecto al tiempo en el sistema. En
este caso de estudio de variación del momento cinético es igual al sumatorio de mo-
mentos exteriores que actúan sobre el sistema ∑ .
, , ; representa la variación de la propiedad específica unitaria
respecto al tiempo dentro del volumen de control. En el caso de encontrarse en régimen
permanente en un sistema indeformable con un fluido incompresible, como es el caso
del estudio de las máquinas hidráulicas, su valor es nulo.
, , , , ∙ ; representa la variación de la propiedad específi-
ca unitaria a través de la superficie de contorno, teniendo en cuenta el intercambio de
propiedad específica.
, es la propiedad específica por unidad de masa en el punto de coordenadas
y en el instante t, situada en el volumen de control o en el . En el caso del análisis de
las máquinas hidráulicas, será como ya se ha citado anteriormente , siendo el
vector posición desde el punto de giro hasta la partícula de estudio y el vector veloci-
dad absoluta de la partícula.
, ; es la densidad de la partícula de fluido en el punto de coordenadas y en el
instante t. En el caso de tratarse de un fluido incompresible esta no varía ni con la posi-
ción ni en el tiempo, por lo que puede tratarse como una magnitud escalar .
, , ; es la velocidad del elemento situado en la superficie de control, en la
posición e instante t. En el caso del estudio de las máquinas hidráulicas, se puede
adoptar la simplificación de propiedades uniformes a la salida y entradas de las diferen-
tes secciones, resultando el producto de , , ∙ igual al caudal circulante en
cada sección de estudio .
29
Bombas funcionando como Turbinas (PAT)
30
Selección de PAT desde la curva de bomba
Teniendo en cuenta que el modulo del vector resultante del producto escalar es
igual a:
Ecuación 5
Por otro lado, la Ecuación 4 puede multiplicarse en ambos miembros por la veloci-
dad angular, quedando la Ecuación 6, en la cual ha sido introducida la relación vecto-
rial de la Ecuación 5:
Ecuación 6
Teniendo en cuenta que el producto entre momento torsor y la velocidad angular es
la potencia aplicada a la máquina (P), en el caso ideal de inexistencia de pérdidas, este
es igual a la Ecuación 7.
= Ecuación 7
siendo: es la constante de la gravedad en / ; altura teórica de Euler aportada
por la bomba en m.
Por lo tanto, la Ecuación 6 puede reescribirse como la Ecuación 8, teniendo en
cuenta la relación de 7.
Ecuación 8
v2
w2
v2m
v1 w1
v1m
v1u v2u
u1
u2
Figura 2. Descomposición de los triángulos de velocidades
En el caso de que la máquina no opere como bomba, sino que funcione como turbi-
na, la altura teórica de Euler ( es igual a la , pero teniendo en cuenta que en
este caso el sentido del flujo se invierte, tal y como muestra la Figura 3
31
Bombas funcionando como Turbinas (PAT)
32
Selección de PAT desde la curva de bomba
Coeficiente de Pfleiderer (μ ) para rodetes de 4 álabes
1
a
0,9
0,8
0,7
μ
0,6
0,5
0,4
0,3
10 14 18 22 26 30 34 38 42 46 50 54 58 62 66 70 74 78 82 86 90
β2 (º)
D1/D2= 0,20 D1/D2= 0,30 D1/D2= 0,40 D1/D2= 0,50
D1/D2= 0,60 D1/D2= 0,70 D1/D2= 0,80 D1/D2= 0,90
Coeficiente de Plfeiderer (μ) para rodetes de 6 álabes
1 b
0,9
0,8
0,7
μ
0,6
0,5
0,4
0,3
10 14 18 22 26 30 34 38 42 46 50 54 58 62 66 70 74 78 82 86 90
β2 (º)
D1/D2= 0,20 D1/D2= 0,30 D1/D2= 0,40 D1/D2= 0,50
D1/D2= 0,60 D1/D2= 0,70 D1/D2= 0,80 D1/D2= 0,90
μ Ecuación 14
Ecuación 15
Ecuación 16
En esta Ecuación 16, es en la que se basan los diferentes métodos empíricos deter-
minación de la curva de operación de la PAT, si los puntos de máxima eficiencia traba-
jando como bomba son conocidos. Del mismo modo, esta ecuación pone de manifiesto
que cualquier máquina operando como turbina para una misma velocidad de giro y
caudal de funcionamiento, la altura recuperada será mayor que la altura proporcionada
cuando opera como bomba.
Al igual que se ha obtenido para los valores de altura, haciendo uso de las leyes de
semejanza puede obtenerse la relación entre caudales, a través de la Ecuación 17.
Ecuación 17
34
Selección de PAT desde la curva de bomba
Figura 6. Pérdidas de carga en el rodete
La altura total ( recuperada por una máquina hidráulica viene definida por la
Ecuación 19:
, ó Ecuación 19
donde , es la altura teórica de Euler en mca, obtenida mediante la ecuación de Eu-
ler.
La Tabla 1 muestra expresiones experimentales para estimar la altura y pérdidas en
una PAT.
Tabla 1 Expresiones para estimar las pérdidas y alturas en una PAT
Ref. Expresión Complemento
Thin et al.,
Ht,E
2008
Ramos et
hf
al., 2000
Thin et al.,
hs kshock=0.2
2008
Subíndice 1 y 2 hace referencia a la entrada y salida del fluido en el rodete; es el ángulo del álabe; es el
número de canales entre álabes; es el ancho del rodete.
35
Bombas funcionando como Turbinas (PAT)
El análisis de pérdidas de carga dentro de una PAT fue abordado por Pérez-Sánchez
et al., 2018 desarrollando un análisis de la altura absorbida por rodete y las pérdidas
que se ocasionaban a la entrada de la máquina, en la voluta y en la salida. La Figura 7
pone de manifiesto el análisis desarrollado mediante fluidodinámica computacional
(CFD) en un modelo calibrado de una PAT cuyo ns=51rpm (m, kW) (Figura 7a). La
Figura 7b muestra los valores para una determinada velocidad de rotación, mientras
que la Figura 7c muestra el porcentaje de la altura recuperada existente en pérdidas y
aprovechada por el rodete en función de la velocidad de giro de la máquina.
Figura 7. (a) Modelo de PAT ns=51 rpm. (b) Alturas en función de localización.
(c) Porcentaje de altura en función de la localización
36
Selección de PAT desde la curva de bomba
Estos parámetros son conocidos a través del ensayo de la máquina hidráulica en la-
boratorio y dependen fundamentalmente de su número específico de revoluciones (ns).
Este valor viene determinado por la Ecuación 27:
/ Ecuación 27
37
Bombas funcionando como Turbinas (PAT)
Dicha ecuación puede reescribirse en la Ecuación 32, teniendo en cuenta las Ecua-
ciones 24, 25 y 26.
tanφ tanφ Ecuación 32
Por tanto, si de una máquina se disponen los parámetros de Suter, la altura y rendi-
miento de la máquina operando como turbina viene establecido por las Ecuaciones 33
y 34.
Ecuación 33
ῃ ῃ ϕ Ecuación 34
La Figura 9 muestra un ejemplo de los valores de los parámetros de Suter para una
máquina con un número específico de revoluciones de 230 rpm.
38
Selección de PAT desde la curva de bomba
Figura 9. Ejemplo de parámetros de Suter para una PAT (ns 230 rpm (m, kW)).
Adaptado de (Ramos y Borga, 2000)
3.2. Coeficientes adimensionales
La aplicación del teorema ∏ de Vaschy-Buckingham a la ecuación fundamental de las
máquinas hidráulicas permite obtener dos parámetros adimensionales que aunque están
referidos a la condición de semejanza restringida (i.e., aquella que únicamente tiene en
cuenta la semejanza geométrica y cinemática). La experimentación con diferentes tipos
de máquinas ha permitido correlacionar las curvas caudal-altura ( ) y eficiencia-
caudal . Estos coeficientes son los denominados coeficiente de caudal ( y coe-
ficiente de altura ( ) (Mataix, 2009) definidos por las Ecuaciones 35 y 36.
Ecuación 35
Ecuación 36
39
Bombas funcionando como Turbinas (PAT)
1.00
nst=10.87 nst=28.1 nst=35.7
0.90 nst=30.1 nst=47.06
nst=18.5 nst=31.9
nst=18.6 nst=38.1
0.80 nst=70
0.70
nst=57.6
0.60
Eficiencia (h )
0.50
0.40
0.30
0.20
0.10
0.00
18
nst=10.87
nst=18.6 nst=30.1
16
nst=18.5 nst=28.1
nst=38.1
14
nst=31.9
nst=41.1
12
Número de Altrua (y )
nst=70
10
nst=57.6
8 nst=47.06
2 Curva tendencia BEP
0
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70
Número de Descarga (ϕ)
4. Potencias y rendimientos
El análisis energético de máquinas hidráulicas operando como motoras es crucial, ya
que es necesario conocer los diferentes términos de alturas, rendimientos y potencias
que intervienen en un sistema de micro generación. Por ello, empleando un esquema de
sistema de generación a pie de presa (Figura 11), se definen las siguientes alturas.
40
Selección de PAT desde la curva de bomba
Altura Bruta ( ); es la altura total disponible entre el punto de origen (en este
caso la presa) y la descarga (en el caso de la figura un río, pero puede ser un
punto presurizado. La determinación de viene definido por un balance de
energía entre la sección B y S, sin considerar pérdidas de carga.
Altura Neta ; es la altura hidráulica disponible para ser aprovechada por la
máquina hidráulica, una vez las perdidas por fricción y singulares han sido de-
traídas de nuestro sistema, puede determinarse mediante el siguiente balance
(Ecuación 37)
Ecuación 37
donde es la energía por unidad de peso en la sección i en mca; es la cota
geométrica en la sección i en m; es la preión en la sección i en mca; es el
término cinético de la sección i en mca.
Altura útil ( ; es la altura realmente aprovechada por la máquina hidráulica y
depende del rendimiento manométrico de la misma. El rendimiento manométri-
co para una turbina viene relacionado por la Ecuación 38
Ecuación 38
41
Bombas funcionando como Turbinas (PAT)
Definidas las alturas, las potencias existentes en un salto hidráulico son (Figura 12):
42
Selección de PAT desde la curva de bomba
Conclusiones
En el presente capítulo se han descrito los fundamentos de la caracterización de las
bombas funcionando como turbinas, concebidas a partir de sus curvas de funciona-
miento. Las curvas, deducidas desde los triángulos de velocidades a la entrada y a la
salida de la máquina y analizados a través de la ecuación de Euler, proporcionan las
curvas teóricas de las máquinas. Sin embargo, para que estas curvas se transformen en
curvas reales, es necesario conocer ciertos parámetros adimensionales estimados a
partir de diversas experimentaciones que permiten conocer el funcionamiento completo
de la máquina, incluidos aspectos energéticos tales como alturas y potencias. El cono-
cimiento de estas curvas características, a falta de ensayos concretos de las bombas
funcionando en modo turbina, permitirán a los diseñadores disponer de información
bien aproximada del comportamiento de las máquinas cuando se instalen en sus siste-
mas correspondientes, permitiendo determinar los puntos de funcionamiento futuros de
las mismas para cada instalación.
43
Selección de máquinas para funcionamiento
como turbina a partir del comportamiento
como bomba
Pérez-Sánchez, Ma,, Sánchez-Romero, FJb, Ramos, HMc, .López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
En epígrafes previos del presente libro, se ha descrito el funcionamiento de las bombas en mo-
do turbinas y su necesidad y aplicabilidad en el conjunto de la red. Sin embargo, es importante
para el lector poder inferir el comportamiento completo de las máquinas cuando funcionen co-
mo turbinas (desconocido a priori) con el de funcionamiento como bombas (proporcionado por
el fabricante y accesible para el público en general). En esa línea, el presente tema describe una
metodología estructurada para la determinación de los puntos concretos de una máquina que
operará como turbina, a partir de aquellos operando como bomba. Esta metodología aplica, tan-
to los puntos de máximo rendimiento o potencia, como el conjunto completo de la curva.
Cuando el lector concluya el aprendizaje de este capítulo, será capaz de obtener (por métodos
empíricos, con un grado de aproximación importante al que luego será el comportamiento
real), el funcionamiento de una máquina cuando se conocen las curvas como bomba y se desea
inferir las de turbina.
Palabras clave: Coeficientes de caudal, altura y potencia; métodos empíricos; curvas
motrices; BEP.
Abstract
The operation of the pumps in turbine mode, and their need and applicability in the networks
as a whole, have been described in previous sections of this book. However, it is important
for the reader to be able to infer the complete behaviour of the machines when they work as
turbines (unknown in advance) if pump curve is only known (provided by the manufacturer
and accessible in general). This chapter describes a structured methodology for the determi-
nation of the specific points of a machine that will operate as a turbine, from those operating
as a pump. This methodology-application of the technology, both the points of maximum ef-
ficiency or power, as well as the complete set of the curve. When readers conclude the learn-
ing of this chapter, they will be able to obtain the operation of a machine in turbine mode,
using empirical methods, which have a large degree of approximation to the real behaviour.
Keywords: Discharge, height and power coefficients. Empirical methods, motor curves, BEP.
45
Bombas funcionando como Turbinas (PAT)
46
Selección de PAT desde la curva de bomba
, H
,
,
, ,
‐Q +Q
Figura 1. Curvas operando como bomba (derecha) y turbina (izquierda)
Ante la falta de catálogos que permitan la selección de bombas operando como tur-
binas, se debe recurrir a estimar el punto de operación de la PAT a partir de la curva
conocida de la máquina hidráulica operando como bomba. Para ello, el usuario debe
servirse de métodos empíricos que estimen el punto de operación de mayor eficiencia
de una máquina operando como turbina a partir de su modo bomba. Este hecho añade
incertidumbre a la selección, como consecuencia que el propio diseño de la máquina
tiene una importante influencia de la eficiencia. El simple hecho que el rodete de la
máquina tenga un número de álabes diferentes, puede causar que la eficiencia sea dife-
rente en dos máquinas similares. Por tanto, desde un principio, el usuario debe asumir
que los métodos que se presentan a lo largo del capítulo tienen un error y que, por tan-
to, deben usarse con cautela a la hora de desarrollar estudios energéticos, debiéndose
siempre recalcular los mismos una vez se conozcan las curvas reales de las máquinas.
El objetivo del presente capítulo está centrado en dos aspectos: en primer lugar, se
pretende que el lector conozca los diferentes métodos existentes para conocer el punto
de máxima eficiencia de una PAT en su funcionamiento como turbina a partir del BEP
de una máquina hidráulica operando como bomba. El segundo objetivo planteado,
consiste en describir detalladamente una metodología de trabajo para, a partir de los
caudales de trabajo, poder seleccionar una bomba que opere como turbina.
47
Bombas funcionando como Turbinas (PAT)
Tabla 1. Métodos empíricos para obtener puntos de operación como turbinas a partir de fun-
cionamiento como bomba (Romero et al., 2017)
Autor
Stephanoff 1 1
1
(Nautiyal et al., 2010) η η
Mc. Claskey
1 1
(Ortiz Flórez & 1
η η
Abella Jiménez 2008)
Alatorre-Frenk 0.85η 0.385 1 0.03
1
(Nautiyal et al., 2010) 2η . 0.205 0.85η 0.385 η
Sharma-Williams 1 1
. . 1
et al., 2010) η η
MICI
(Ortiz Flórez & 0.9-1.0 1.56-1.78 0.75-0.80
Abella Jiménez 2008)
Yang et al. 1.2 1.2
-
(Pugliese et al.,2016) η. η.
(continúa)
48
Selección de PAT desde la curva de bomba
Tabla 1. Métodos empíricos para obtener puntos de operación como turbinas a partir de fun-
cionamiento como bomba (Romero et al., 2017)
Hancock 1 1
-
(Pugliese et al., 2016) η η
Grover
2.379 0.0264n 2.693 0.0229n -
(Nautiyal et al., 2010)
Hergt 1.6 6
1.3 1.3 -
(Nautiyal et al., 2010) n 5 n 3
Ecuación 7
.
.
. para Ecuación 8
En el caso de la expresión para valores mayores a >50 (m, kW) puede con-
ducir a valores de eficiencia erróneos mayores a la unidad.
Los coeficientes expuestos en las ecuaciones anteriores han sido obtenidos a partir
de los puntos de operación de casi un centernar de máquinas, obteniendo coeficientes
de correlación en todos los casos superiores al 96%. La Figura 2 muestra los puntos
experimentales utilizados en todos los ratios así como la predicción del número
específico de la máquina operando como turbina ( , definido por la Ecuación 9.
. Ecuación 9
49
Bombas funcionando como Turbinas (PAT)
Figura 2. Coeficientes para (a), caudal (b), altura (c) y eficiencia (d) en función de
Como en ocasiones, se opera en sentido inverso, a partir de los datos de operación
en modo turbina, se debe seleccionar la bomba, los autores proponen los siguientes
coeficientes en función del número especifico de turbina ( establecido de
acuerdo al punto de operación más probable. Los coeficientes vienen definidos por la
expresiones:
.
Ecuación 10
.
Ecuación 11
.
. para Ecuación 12
En el caso de la expresión para valores mayores a >50 (m, kW) puede con-
ducir a valores de eficiencia erróneos mayores a la unidad.
En este caso, la relación entre y es igual a:
. Ecuación 13
La Figura.3 muestra los coeficientes en función del número específico como turbina
50
Selección de PAT desde la curva de bomba
Figura 3. Coeficientes para (a), caudal (b), altura (c) y eficiencia (d) en función de
2 2
Ecuación 15
51
Bombas funcionando como Turbinas (PAT)
(Derakhshan y Nourbakhsh 2008; Novara y McNabola, 2018; Rossi, et al., 2019; Rossi y
Renzi 2018) proponen curvas adimensionales de ⁄ , ⁄ y ⁄ , que permi-
ten estimar directamente las curvas de funcionamiento en función del punto BEP.
Estas investigaciones relacionan los coeficientes adimensionales con el número
específico de la máquina operando como turbina ( . El uso de estos valores, como
los mostrados en las Figuras 4 y 5, pueden inducir a cierto error, pero pueden ser una
herramienta útil para realizar estudios energéticos preliminares o simulaciones de
operación en redes de distribución que den pie a estudiar, con más detalle, el uso de
PAT en el sistema.
En la actualidad, existe un mayor número de curvas que no han sido mostradas en la
Figura 4 y Figura 5 para poder ser interpretadas. No obstante, el mayor uso de curvas
implica que entre el número específico de turbina obtenido pueda desarrollarse una
interpolación más cercana y de este modo ajustar mejor los valores obtenidos (Pérez-
Sánchez et al., 2017).
52
Selección de PAT desde la curva de bomba
53
Bombas funcionando como Turbinas (PAT)
4. Metodología de selección
Los dos apartados anteriores son importantes para conocer los puntos de operación así
como las curvas completas de operación de las máquinas. No obstante, es importante
desarrollar una metodología de selección de máquinas que vayan a operar como
turbinas a partir de un catálogo de bombas. En esta sección del capítulo, se resume y
modifica la metodología propuesta por Romero et al., (2017) donde a partir de los
valores de caudal y altura potencialmente recuperable, se puede seleccionar la bomba
que operando como turbina permita mejorar la eficiencia energética del sistema. La
metodología se basa en los métodos empíricos para fijar los puntos de operación de la
máquina en su modo turbina así como conocer el valor de . La Figura 4.6 recoge el
diagrama de flujo que define los pasos de la metodología así como la entrada de datos
necesaria. La metodología está establecida en 5 fases:
A) Determinación del punto de operación de la máquina operando como bomba;
para afrontar esta primera fase, debe conocerse el caudal y altura recuperable
por la máquina (o los caudales y alturas si varian a lo largo del tiempo) más
frecuente (Entrada 1, Figura 6). Aplicando, cualquiera de los métodos empíricos
propuestos anteriormente, se obtienen el punto (o puntos) de la máquina
operando en modo bomba ( ; ; ). Se recomienda emplear varios
métodos para comprobar la variabilidad de los puntos obtenidos.
B) Selección de máquina operando en modo bomba; una vez son definidos el punto
de máxima eficiencia de la bomba como resultado del paso A (Salida 1), a partir
de un catálogo de bombas, se lleva a cabo la selección de la máquina.
Seleccionada la bomba, se establecen los valores aportados por el fabricante
como valores de partida que sustituyen los obtenidos en el paso A.
C) Recálculo de los parámetros operando como turbina. Con los valores obtenidos
en B, se determinan los valores , mediante las Ecuaciones 6 a 8. Con
estos valores se determinan los valores empíricos ; ; ,
comparándose con los valores iniciales del caso de estudio (Entrada 1). Si estos
valores son similares, se pasa a la siguiente etapa.
54
Selección de PAT desde la curva de bomba
55
Bombas funcionando como Turbinas (PAT)
concretamente del motor eléctrico que lleve instalado) o directamente las curvas
adimensionales de ⁄ y ⁄ , en función de ⁄ , y a partir del
valor del estimar las curvas de funcionamiento (Figura 5).
F) Determinación de los términos de mejora energética; conocidos los valores de
las diferentes fases (caudal, altura y eficiencia en modo turbina) puede llevarse a
cabo el estudio de la mejora energética. En este caso, el estudio de la energía
recuperable así como la potencia a instalar es fundamental, con el objetivo de
poder llevar a cabo el análisis de viabilidad previo. Para ello, es determinante
que el análisis de la energía recuperable, tenga un valor de eficiencia
conservardor, teniendo en cuenta todos los elementos a instalar (variadores,
condensadores -si el sistema es aislado de la red- motor eléctrico,
transformadores, entre otros). Por tanto, la energía recuperable vendrá definida
por la Ecuación 17, mientras que la potencia por la Ecuación 18:
E Ecuación 17
P γQ H η η Ecuación 18
donde: E es la energía recuperable en kWh/año; es la eficiencia eléctrica del
conjunto del sistema instalado; t es el tiempo de operación en s de la máquina; P
es la potenica de la máquina en kW.
5. Caso de estudio
Esta sección tiene por objetivo desarrollar un ejemplo en el que a partir de unos valores
supuestos de caudal y altura recuperable, se puede seleccionar una máquina que opere
en modo turbina. La solución utiliza la metodología expuesta en el apartado anterior a
través del conocimiento de sus curvas como bomba. Para llevar a cabo el análisis, se
propone el caso siguiente:
Dos depósitos se encuentran conectados mediante una conducción de fundición de
diámetro interior 200 mm, con una rugosidad absoluta de 0.1 mm y una longitud de
6800 m. El depósito superior se encuentra a la cota 90 msnm, el depósito inferior a la
cota 40 msnm. Se conocen que para un caudal de 25 l/s, las pérdidas de carga por
fricción son 21.33 mca y las pérdidas singulares de la conducción, suponiendo un 15%
de las pérdidas por fricción, son 3.20 mca. Por tanto, para regular este caudal, existe
instalada una válvula reguladora que introduce en el sistema 25.47 mca. Los gestores
del sistema pretenden cambiar esta válvula reguladora por una bomba operando como
turbina. Se pide seleccionar una máquina que supla estas condiciones.
56
Selección de PAT desde la curva de bomba
57
Bombas funcionando como Turbinas (PAT)
58
Selección de PAT desde la curva de bomba
59
Bombas funcionando como Turbinas (PAT)
Aplicando los coeficientes calculados según Pérez-Sánchez los puntos BEP en fun-
cionamiento como turbina quedarían como:
214 ⇒ 1.55 17.70 27.44 /
1.74 14.30 24.88
0.91 0.78 0.71
21.55 ,
205 ⇒ 1.52 17 25.92 /
1.71 13 22.27 0.91 0.77 0.70 22.77 ,
En la Tabla 4. se muestran los resultados para las bombas 65-20, 65-26h, 50-26h,
80-16 y 80-20. Se calcula también el error respecto al punto de funcionamiento tanto
en caudal, altura y el total. Los menores errores se obtienen para la bomba 65-26h y
65-20, descartándose el resto de modelos ya que los puntos BEP proporcionados se
alejan del punto de funcionamiento deseado. La elección podría simplificarse a selec-
cionar aquella con un menor error, pero por ejemplo en los casos donde exista varia-
ción del caudal o que las eficiencias sean diferentes, es necesario disponer de las curvas
Q-H y Q-η para la selección final de la bomba.
Tabla 4. Resultados para diferentes bombas de la misma serie
Modelo ⁄ ,
65-20 214 17.70 14.30 0.780 26.23 1.55 1.74 0.91
65-20 205 17.00 13.00 0.770 27.61 1.52 1.71 0.91
65-26h 220 16.00 15.00 0.715 24.06 1.59 1.79 0.89
65-26h 205 15.00 13.00 0.700 25.94 1.55 1.75 0.90
50-26h 220 7.50 15.00 0.550 16.48 1.81 2.03 0.84
80-16 174 27.50 8.75 0.850 47.26 1.31 1.47 0.98
80-20 205 27.00 13.00 0.805 34.80 1.43 1.60 0.94
Modelo ⁄ , % % %
60
Selección de PAT desde la curva de bomba
Figura 10. Curvas Q-H (a) y Q-η (b) de las bombas seleccionadas
Figura 11. Evolución de la eficiencia de las máquinas de inducción con su potencia de salida
(Adaptado de Fernandes et al., 2019)
En la Tabla 5 se muestran los resultados de potencia hidráulica recuperada en fun-
ción de la curva de trabajo seleccionada. Para el modelo 65-26h, la altura recuperada es
61
Bombas funcionando como Turbinas (PAT)
Conclusiones
En el presente capítulo, se han descrito los diferentes métodos empíricos para seleccio-
nar una bomba que opere como turbina, conocidos los puntos de operación. Relaciona-
do con este aspecto, los autores han publicado unas nuevas ecuaciones de regresión
para determinar los coeficientes que permiten determinar los puntos de operación como
bomba a partir de los datos de entrada de máquina operando como turbina. Presentados
los métodos empíricos, se ha abordado la descripción de una metodología modificada a
la presentada por Romero et al., (2017). Esta metodología permite determinar de forma
previa, la selección de máquinas, definición de curvas a partir de los números adimen-
sionales y determinar un primer análisis energético. La metodología propuesta ha sido
aplicada a un caso de estudio en concreto, mostrando las capacidades de la estrategia
de determinación de la curva de la máquina funcionando como turbina, conocido el
modo de funcionamiento como bomba.
62
PAT operando en régimen permanente con
variación de velocidad. Aspectos de
regulación
Pérez-Sánchez, Ma,, Sánchez-Romero, FJb, Ramos, HMc, .López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
Conocidos en profundidad los aspectos de caracterización de las bombas funcionando como
turbinas a través de sus curvas motrices, se hace necesario introducir aspectos de regulación.
El control sobre caudales, alturas y potencias en las máquinas es fundamental para su insta-
lación en redes con regímenes variables. Esto permitirá una adaptación a las posibilidades de
producción de cada punto. El presente capítulo introduce sobre las estrategias de regulación
controlando la curva motriz a través de regulación hidráulica y eléctrica (combinando má-
quinas asociadas en serie y en paralelo con control de velocidad) con regulación con válvu-
las asimismo en serie y paralelo, que permitan la instalación más adecuada a cada punto de
funcionamiento.
Abstract
After an in-depth knowledge of the characterization aspects of the pumps operating as tur-
bines through their motor curves, it is necessary to introduce the regulation aspects. The con-
trol over flows, heights and powers in the machines is fundamental for its installation in
networks with variable regimes, adapting to the production possibilities of each point. This
chapter introduces regulation strategies by controlling the motor curve through hydraulic and
electrical regulation (combining associated machines in series and parallel with speed con-
trol) using serial and parallel valves, in order to control the installation more appropriate to
each operating point.
Keywords: Regulation; Affinity laws; serial and parallel association; Best Efficiency Line;
Best Power Curve.
63
Bombas funcionando como Turbinas (PAT)
1. Introducción y objetivos
En los capítulos anteriores se ha abordado el estudio fundamental de las máquinas hi-
dráulicas, así como la selección de máquinas conocido un punto de operación. No obs-
tante, es conocido que los sistemas de distribución pueden trabajar en condiciones de
caudal y/o altura recuperable variable, por lo que la regulación de los sistemas de recu-
peración se hace igualmente necesaria. Por tanto, desde el punto de vista del gestor de
sistemas de distribución, no solo se debe saber seleccionar una máquina a partir del
punto o zona de operación, sino que debe conocerse la mejor curva de consigna de
funcionamiento para maximizar la energía recuperada. Este capítulo se centra en este
aspecto, y su estudio permitirá acceder a la optimización de sistemas de distribución
que se verá en el siguiente capítulo. Por tanto, los objetivos del presente capítulo son:
Conocer los diferentes escenarios frente a caudal y altura que pueden encon-
trarse un sistema de recuperación mediante PAT.
Definir los diferentes tipos de asociación en serie y paralelo.
Describir la regulación que puede llevarse a cabo en una PAT o grupo de recu-
peración.
Establecer las ecuaciones clásicas para las curvas con variación de velocidad y
modificadas de semejanza
Definir la regulación de sistemas de recuperación sobre la línea de máxima efi-
ciencia y máxima potencia.
64
Aspectos de regulación
Figura 1. Operación PAT a caudal constante (b) sustituyendo válvula regulación (a)
b) Caso 2. Flujo variable a lo largo de tiempo, sustituyendo una válvula regulado-
ra; este segundo caso, la PAT o grupo de PAT van instaladas sustituyendo a la
válvula de regulación (Figura 2, Caso B). A priori, la oscilación de caudal y
altura recuperable a lo largo del tiempo es conocida, como consecuencia de la
instrumentalización que contenga la válvula reguladora (e.g., sensores de
presión, caudalímetro) o conociendo el grado de apertura, diferencia de presión
y curva de coeficiente de caudal de la válvula en función de la apertura.
Figura 2. Operación PAT a caudal variable (b) sustituyendo válvula regulación (a)
c) Caso 3. Flujo variable a lo largo de tiempo sin existir válvula reguladora; en
este tercer caso no se lleva a cabo la sustitución de la válvula regulación, ya que
no existe dentro de la red. Esta situación, se presenta cuando existe una red de
distribución que total o parcialmente, se ve sometida a presiones elevadas como
consecuencia de que, para garantizar la presión en los nudos más desfavorables
del sistema, estas zonas profundas tengan piezométricas por encima de las mí-
nimas necesarias.
Esta situación causa que, a priori, no se conozca la localización exacta de la
PAT o grupo de PAT (Figura 3), ya que en función de su localización y tipo de
máquina dependerá la presión aguas abajo, así como la energía potencialmente
recuperable.
65
Bombas funcionando como Turbinas (PAT)
66
Aspectos de regulación
67
Bombas funcionando como Turbinas (PAT)
60 a 60 b
55 55
50 50
45 45
Altura recuperable (mca)
Altura recuperable (mca)
40 40
35 35
30 30
25 25
20 20
15 15
10 10
5 5
0 0
0 5 10 15 20 25 30 35 40 45 50 0 5 10 15 20 25 30 35 40 45 50
Caudal (l/s) Caudal (l/s)
68
Aspectos de regulación
69
Bombas funcionando como Turbinas (PAT)
Ecuación 1
⁄
donde: es el coeficiente de caudal de la válvula en . ; es el caudal circulante
por la válvula en ⁄ ; densidad del agua (1000 / ); es la aceleración de la
gravedad (9.81 / ); es la pérdida de carga, igual a la altura recuperada por el
grupo de PAT en Pa.
2.2.2 Consignas de operación
El funcionamiento de sistemas de recuperación basado en máquinas en paralelo requie-
re de una serie de consignas de operación, programadas en el sistema y que deben esta-
blecer su funcionamiento en función de la máquina seleccionada, del número de má-
quinas instaladas y del caudal instantáneo en línea. Estas consignas son:
Caudal mínimo ( ). Este parámetro vendrá determinado por la máquina se-
leccionada.
Caudal máximo ( ). Este parámetro dependerá de altura mínima disponible
recuperable ( , así como de la curva de operación elegida. Este parámetro
varía en función de la curva Q-H establecida, provocando que los rangos de
caudales para la puesta en marcha y parada de nuevas máquinas sean diferentes
según la curva de operación establecida. La altura mínima disponible recupera-
ble ( se podrá obtener mediante el análisis hidráulico y energético reali-
zado, o podrá ser establecida por el usuario, siempre y cuando el valor propuesto
sea igual o menor al real existente.
Curva de operación. Se deben establecer las ecuaciones polinómicas de regre-
sión Q-H y Q-, para la determinación en cualquier instante de la potencia gene-
rada. A este parámetro se le denomina “Estrategias de Alturas”, y define para la
máquina seleccionada, el cálculo de la altura recuperada y la eficiencia en fun-
ción del caudal instantáneo turbinado por la máquina, mediante las ecuaciones
de regresión establecidas para estas curvas.
Criterio de puesta en marcha o parada de máquinas. En cualquier instante, el sis-
tema debe decidir cómo distribuye los caudales de la línea entre las máquinas exis-
tentes y el bypass. Esta distribución se puede realizar según diferentes criterios ta-
les como: maximizar el volumen turbinado, maximizar la energía recuperada,
disminuir el número de máquinas en funcionamiento y el tiempo de turbinado. A
este parámetro se le denomina “Estrategias de Caudales”, y establece en el sistema
las normas para la puesta en marcha o parada de las máquinas instaladas.
70
Aspectos de regulación
Figura 8. Tipos de regulación. (a) RH con válvula reductora en serie; (b) RH con
válvula reductora en paralelo; (c) RE mediante variación de velocidad
71
Bombas funcionando como Turbinas (PAT)
72
Aspectos de regulación
En cuanto a los rendimientos, si son máquinas diferentes, cada máquina tendrá una
eficiencia diferente con el caudal circulante total que pase por cada máquina. Si las
máquinas son iguales, la eficiencia será la misma para todas ellas.
Ecuación 9
Si la instalación de las máquinas es en paralelo, al igual que en la asociación en se-
rie, puede darse el caso que las bombas puedan ser iguales o diferentes. En este caso, lo
recomendable es utilizar máquinas idénticas que giren a la misma velocidad (bien sea
la nominal o diferente mediante el uso de variadores de frecuencia (conexión a red) o
en función de la resistencia y condensadores (cuando los sistemas están desconectados
de la red). En este caso, la ecuación que define la asociación de x máquinas iguales es:
Ecuación 10
73
Bombas funcionando como Turbinas (PAT)
Ecuación 13
Ecuación 14
Ecuación 16
Definidas las leyes de semejanza para cambio de velocidad de rotación, las curvas
(altura y eficiencia en función del caudal) quedan según las siguientes ecuaciones:
Ecuación 17
ῃ Ecuación 18
74
Aspectos de regulación
Figura 10. Ejemplo de números adimensionales q (a), h (b), p (c) y e (d) para una
máquina con igual a 51 rpm (m, kW). Adaptado de Pérez-Sánchez et al., (2018a)
Ecuación 20
Ecuación 21
La Figura 11 muestra un ejemplo de dichas curvas para una máquina radial con un
igual a 51 rpm (m, kW), ensayada por Pérez-Sánchez et al., (2018a).
75
Bombas funcionando como Turbinas (PAT)
Figura 11. Ejemplo de números adimensionales q (a), h (b), p (c) y e (d) para una
máquina con igual a 51 rpm (m, kW). Adaptado de Pérez-Sánchez et al., (2018a)
76
Aspectos de regulación
Figura 12. (a) Concepto de BEL y BEH. (b)Ejemplo de BEL y BEH para una máquina
con igual a 51 rpm (m, kW). Adaptado de Pérez-Sánchez et al., (2018a)
En (regimen nominal)
BEP
En min y max
max
BEH
BEH min
BPC
Runaway Curva Runaway
Eficiencia
min
Hmin Zona de operación
de la turbina (VOS)
Figura 13. Esquema curvas y zona de operación para una turbina genérica
77
Bombas funcionando como Turbinas (PAT)
78
Aspectos de regulación
Situación I a Situación II b
H Zona I Zona II Zona III
H Zona I
(no turbinado) (turbinado por 1 PAT) (no turbinado) Zona II Zona III
HTi HTi
Hmax
Altura recuperable
Altura recuperable
HTmin HRi
Hmax
Hmin
HRi
Hmin
Q Q
Qmin Qi Qmax Qmax_en_línea Qmin Qi Qmax Qmax_en_línea
Caudal Caudal
Hmax Hmax
HTmin HTmin
Hmin Hmin
Q Q
Qmin BPC BEH Qmax_en_línea Qmin 2Qmin Qmax 3Qmin 2Qmax 3Qmax Qmax_en_línea
Qmax Qmax Qmax
Caudal Caudal
79
Bombas funcionando como Turbinas (PAT)
Conclusiones
En este capítulo se han presentado los diferentes escenarios de funcionamiento que
pueden darse en un sistema de recuperación, los cuales dependen fundamentalmente,
de la variabilidad de los caudales a lo largo del tiempo, y con ello de las condiciones de
variación de la altura recuperable, que depende a su vez de la energía en cabecera. A
partir de conocer los puntos de operación, se han definido los diferentes tipos de regu-
lación (hidráulica, eléctrica e hidráulica-eléctrica) que pueden desarrollarse. La presen-
tación de los dos últimos tipos de regulación, hace necesario la variabilidad de la velo-
cidad de giro de la máquina y, por tanto, se han abordado las leyes de semejanza, así
como la propuesta de funciones que mejoran la predicción analítica de las curvas. Fi-
nalmente, una vez presentadas las ecuaciones de semejanza, se ha llevado a cabo la
descripción de la regulación sobre las líneas de máxima eficiencia y máxima potencia.
Este capítulo junto con el capítulo anterior son cruciales para poder abordar el siguiente
capítulo. Este llevará cabo la descripción de la metodología de trabajo para la optimi-
zación de sistemas de recuperación, teniendo en cuenta la mejor localización, así como
el desarrollo de un caso de estudio.
80
Localización, selección y funcionamiento de
PAT en sistemas de distribución a presión
Pérez-Sánchez, Ma,, Sánchez-Romero, FJb, Ramos, HMc, .López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
La descripción del comportamiento de las bombas funcionando como turbinas cobra sentido
pleno cuando las máquinas se instalan en las redes, con objeto de sustituir válvulas o como
elementos de control de la presión. Sin embargo, para decidir la ubicación de las mismas su
tipología y el número de máquinas, se hace necesaria la ayuda de herramientas de optimiza-
ción por la complejidad del problema propuesto. El presente capítulo describe una metodo-
logía para realizar esta optimización, basada en el procedimiento del “Recocido Simulado” o
“Simulated Annealing” (proceso en inglés) que se describe en detalle e implementa en una
red real, para mostrar la gran capacidad potencial de la herramienta descrita para optimizar el
comportamiento energético de las redes.
Abstract
The description of the PAT behaviour makes full sense when the machines are installed in
the networks in the replacement of valves, or as pressure control elements. However, in or-
der to decide their location, type and number of machines, the help of optimization tools is
needed, due to the complexity of the proposed problem. This chapter describes a methodolo-
gy to perform this optimization, using Simulated Annealing procedure, which is described in
detail and implementation in a real network, to show the great potential capacity of the tool
described for energy behaviour of networks.
81
Bombas funcionando como Turbinas (PAT)
1. Introducción
A lo largo del presente libro, se ha hecho enfásis que el consumo de energía en los
sistemas de distribución presurizados (i.e., abastecimiento, riego) suponen cerca del 3%
(Nogueira y Perrella, 2014), haciéndose imprescindible mejorar la eficiencia final de
los sistemas. La necesidad de incrementar esta eficiencia ha provocado que en los
últimos años se hayan diseñado y aplicado diferentes estrategias para el aumento de la
eficiencia de este tipo de instalaciones. Estas estrategias dependen del tipo de sistema
de distribución. En sistemas con bombas en cabecera, los métodos aplicados buscan
minimizar la energía inyectada al sistema con diversas técnicas, principalmente
mejorando el funcionamiento y rendimiento de los sistemas de impulsión. En otros
casos, organizando las zonas de distribución en función de su altura u otras variables
(Moreno et al., 2010; Pardo et al., 2013; Jiménez-Bello et al., 2015). En sistemas por
gravedad el incremento de la eficiencia se produce mediante la reducción de presión en
puntos determinados del sistema, uniformizando las presiones existentes y disipando la
energía excedente del sistema, disminuyendo así las fugas y reduciendo la probabilidad
de rotura en las conducciones (Carravetta et al., 2012; Carravetta et al., 2013; Cabrera
et al., 2014b).
Tal como ha sido descrito, la reducción de presión en los sistemas de distribución se
ha realizado tradicionalmente mediante válvulas reductoras de presión (PRV),
colocadas en puntos estratégicos del sistema. La recuperación de esta energía ha sido
estudiada por diferentes autores, proponiendo el reemplazo de las PRV por sistemas
hidroélectricos (Carravetta et al., 2013) de diversos tipos (Chen et al., 2013; Malavasi
et al., 2014; Corcoran et al., 2015; Malavasi et al., 2018). Dentro de estos estudios
destacan el reemplazo de las PRVs mediante bombas trabajando como turbinas (PAT)
(Carravetta et al., 2012; 2013; 2014c).
Estos estudios se han centrado principalmente en sistemas de abastecimiento
urbano. En sistemas de distribución para riego, se han realizado también algunos
estudios (Estrada, 2014; Butera y Balestra, 2015; Pérez-Sánchez et al., 2016). Estos
estudios establecen diferentes escenarios de caudales en las conducciones, y a partir de
estos escenarios determinan las presiones, y por lo tanto el análisis energético del
sistema. La variabilidad de los caudales en sistemas de abastecimiento ha sido
ampliamente estudiada, siendo muy diferente en los sistemas de distribución de agua
para riego, ya que estos sistemas presentan una gran variabilidad temporal debido
principalmente a que dependen de los parámetros agronómicos, siendo compleja su
determinación. Algunos de estos estudios consideran caudales medios (Estrada, 2014)
y otros autores los analizan mediante caudales discretos en determinados períodos de
consumo (Pérez-Sánchez et al., 2016).
La ubicación de estos sistemas de recuperación dentro de la red de distribución ha
sido estudiada por diversos autores, proponiendo diferentes metodologías para su
determinación óptima. Estas metodologías se basan generalmente en la maximización
de la energía recuperada, en la optimización de variables económicas o en la
combinación de diferentes variables de tipo energético y/o económico (Giugni et al,
82
Localización, selección y funcionamiento óptimos de PAT
2014; Samora et al., 2016). La dificultad para encontrar las mejores soluciones es
consecuencia de la naturaleza del problema al resolver, ya que aparecen una gran
cantidad de soluciones posibles, que dependen de los múltiples escenarios de
funcionamiento. Las metodologías que referencia la bibliografía se basan en algoritmos
metaheurísticos (Verdagay et al., 2008) como: ‘Simulated Annealing’(Samora et al.,
2016), el ‘tabu search’ (da Conceição y Ribeiro, 2004) o ‘particle swarm optimization’
(Baghaee 2012; 2016; 2017). Similares metodologías han sido utilizadas para la
ubicación de PRVs (Araujo et al., 2006), dimensionado óptimo de diámetros en las
conducciones (Kuo et al., 2003; Reca et al., 2008; Tospornsampan et al., 2007),
minimización del exceso de energía aportado en sistemas bombeados (Prats et al.,
2011), así como en otro tipo de problemas relacionados con los sistemas de
distribución de agua.
Dentro de estas metodologías destaca por su aplicación a este tipo de problemas el
algoritmo ‘Simulated Annealing’ (Scollen y Hargreaves, 2018). Este método está
inspirado en el método de Montecarlo, basándose en la idea original de Metropolis.
Este método está recomendado para problemas con una gran cantidad posible de
combinaciones. La búsqueda heurística que realiza el algoritmo de las mejores
soluciones tiene como analogía al proceso físico de enfriamiento que ocurre en los
metales.
Una vez determinada la ubicación o ubicaciones de los sistemas de recuperación, la
energía recuperada dependerá de las características de la máquina o máquinas que se
instalen. Las curvas de trabajo Q-H y Q- de las máquinas, determinarán finalmente la
energía recuperada, así como la potencia a instalar. En este sentido, Carravetta et al.,
(2013) propuso diferentes estrategias de regulación en función de la variabilidad de los
caudales y alturas recuperables (resumidas en el capítulo anterior). Diferentes autores,
han simulado la existencia de turbinas en base a curvas adimensionales obtenidas del
estudio de varias máquinas (Derakhshan y Nourbakhsh, 2008; Singh y Nestmann,
2010).
El objetivo del presente capítulo es la descripción del proceso de análisis de la
ubicación óptima (al menos teórica) de la maquinaria hidráulica a instalar en sistemas
de recuperación de energía en redes de distribución de agua, sustituyendo PRV por
sistemas compuestos por PAT, recuperando el exceso de energía, aumentando así la
eficiencia de la red.
2. Materiales y Métodos
En el capítulo que nos ocupa, se va a describir con detalle una propuesta para la utiliza-
ción de una de estas técnicas heurísticas para la búsqueda de la ubicación óptima de los
grupos de recuperación. El proceso de análisis para la determinación de dicha ubica-
ción de sistemas de recuperación en sistemas de distribución de agua, se basa en la
técnica del Simulated Annealing (SA) y queda reflejado en la Figura 1. Este proceso
queda dividido en las siguientes fases:
83
Bombas funcionando como Turbinas (PAT)
Fase I.‐ Determinación escenario de funcionamiento del sistema.
* Cálculo caudales
Características red (Input 1)
circulantes y presiones
en la red.
Hábitos de consumo (Input 2) (EPANET toolkit)
(Output 1/Input 3)
Fase II.‐ Balance energético.
* Cálculo altura teórica
recuperable (H)
Presión mínima de
* Cálculo Energías.
funcionamiento en nudos de la red
(Input 4)
* Cálculo puntos de
funcionamiento (Q,H)
(Output 2/Input 5)
Fase III.‐ Determinación ubicaciones óptimas de
los sistemas de recuperación (Optimización I). * Optimización ubicación
sistemas de recuperación
Parámetros de control algoritmo en función de la variable
Simulated Annealing (Input 6) objetivo seleccionada.
(Output 3/Input 7)
84
Localización, selección y funcionamiento óptimos de PAT
Fase IV. Con las ubicaciones óptimas consideradas, se realiza una segunda optimi-
zación (Optimización II), también basada en el algoritmo heurístico ‘Si-
mulated Annealing’, donde se determinan las máquinas reales a instalar,
así como las consignas de funcionamiento necesarias para optimizar las
funciones objetivo consideradas
2.1. Fase I. Determinación de caudales circulantes
Las características geométricas y topológicas de la red deben ser definidas, determi-
nándose el caudal demandado en cada nudo de consumo. Esta información vendrá
marcada por las características específicas de cada nudo. Por ejemplo, para redes de
riego dependerá de la superficie a abastecer, así como de las características del sistema
de riego instalado en parcela (i.e., marco de plantación, goteros por árbol) (Input 1).
Para el cálculo de los patrones de consumo, diversos autores proponen diferentes
metodologías (Davidson y Bouchart, 2006; Ghiassi et al., 2008; Preis et al., 2009; Sanz
y Pérez, 2014). Para este caso se ha utilizado la metodología propuesta por Derakhshan
y Nourbakhsh (2008). Esta metodología se aplica para un año natural y un intervalo de
tiempo prefijado, determinándose en qué intervalos de tiempo están abiertos o cerrados
cada uno de los nudos de consumo (Input 2). La apertura o cierre de las tomas se esta-
blece en función de:
Consumos históricos obtenidos a partir de los registros de los contadores indi-
viduales para cada punto de consumo. En caso de que no exista esta informa-
ción, puede ser obtenida a partir de los consumos teóricos considerados o por
ejemplo en riego, a partir de las necesidades teóricas de los cultivos. En cual-
quier caso, el conjunto de los volúmenes totales de consumo individuales será
conforme al volumen total consumido anual para todo el sistema.
Curvas de probabilidad de uso para cada una de las tomas de consumo. Esta
información describe los hábitos y costumbres de los usuarios, pudiéndose de-
finir la probabilidad semanal de uso, así como la probabilidad horaria de aper-
tura de la toma e incluso los tiempos y número de aperturas diarias. Esta infor-
mación es obtenida a partir de entrevistas con los usuarios, o directamente de
datos reales recogidos por caudalímetros que se encuentran instalados en las lí-
neas de cabeceras u otras ubicaciones.
La aplicación de esta metodología y el uso EPANET toolkit (Sanz y Pérez, 2014)
permite conocer para cada intervalo de tiempo considerado los caudales para cualquier
línea de la red y presiones de los nudos, lo que permite realizar un balance energético
anual del sistema (Output 1).
2.2. Fase II. Balance energético
El balance energético permite la discretización de los diferentes términos de energía,
obteniendo finalmente la altura recuperable y caudal para cada línea o punto de consu-
mo, en cada instante de tiempo. Conociendo estos pares de datos, es posible el estudio
de las ubicaciones y tipo de máquina a instalar, pero siempre asegurando en cualquier
nudo del sistema la presión mínima de funcionamiento (Input 6), así como el caudal
85
Bombas funcionando como Turbinas (PAT)
requerido en los nudos de consumo. Esta condición debe ser siempre considerada ya
que la recuperación de energía no es el principal objetivo de los sistemas de distribu-
ción de agua (San Bruno et al., 2010). Los términos energéticos para los nudos de con-
sumo (tomas) quedan reflejados en las siguientes ecuaciones:
.
Ecuación 1
. Ecuación 2
. Ecuación 3
. Ecuación 4
. Ecuación 5
donde: es la energía potencial en el nudo de consumo cuando la demanda es nula
para un instante de tiempo considerado en kWh; es la energía disipada por
fricción hasta el nudo de consumo para un instante de tiempo en kWh; es la
energía necesaria en el nudo de consumo para un instante de tiempo para asegurar
el caudal demandado requerido en kWh; es la potencia teórica recuperable en el
nudo de consumo para un instante de tiempo en kW; es la máxima energía
teórica recuperable en el nudo de consumo para un instante de tiempo en kWh;
caudal demandado por el nudo de consumo para un instante de tiempo en m3/s;
es la altura geométrica respecto al plano de referencia del nivel de agua en el punto de
suministro (depósito, balsa, embalses) en m; es la altura geométrica respecto al plano
de referencia del nudo de consumo en m; es el intervalo de tiempo considerado en
s; es la presión de servicio en el nudo de consumo para un instante de tiempo
cuando existe demanda en ese nudo en mca; _ es la presión mínima de servicio
en el nudo de consumo para asegurar el caudal demandado en mca; es la altura
recuperable en el nudo de consumo para un instante de tiempo en mca. Esta altura
se obtendrá como:
_ Ecuación 6
Por lo tanto, para los nudos de consumo el balance final quedará como:
Ecuación 7
donde: es la energía total en el nudo de consumo en kWh/año, obteniéndose co-
mo:
∑ Ecuación 8
86
Localización, selección y funcionamiento óptimos de PAT
∑ Ecuación 10
es la energía total teórica recuperable en el nudo de consumo en kWh/año, obte-
niéndose como:
∑ Ecuación 11
siendo: es el número total de intervalos de tiempo considerados para discretizar un
año.
Para las líneas de distribución, el balance energético es similar, pero se debe consi-
derar el término de energía que no es posible recuperar y que es necesario para mante-
ner la presión mínima en los nudos aguas abajo del sistema. El balance energético para
líneas quedará como:
+ Ecuación 12
donde: es la energía total no recuperada en la línea en kWh/año, obteniéndose
como:
∑ Ecuación 13
Ecuación 14
. _ Ecuación 15
donde: es la energía disponible en la línea para un instante de tiempo en kWh;
es el caudal circulante por la línea para un instante de tiempo , igual a la suma
de los nudos de consumo abiertos en el instante aguas abajo de esa línea en m3/s;
_ es la presión mínima de servicio en la línea para asegurar el caudal demanda-
do aguas abajo en mca; es la presión de servicio en la línea para un instante de
tiempo en mca.
La altura recuperable en la línea para un instante de tiempo en mca, se calculará
como:
_ ; _ Ecuación 16
donde: _ es la presión mínima de servicio en la línea para un instante de tiempo
para asegurar la presión mínima requerida en el nudo de consumo más desfavorable
en mca.
La potencia efectiva instantánea en el eje de la máquina para un instante de tiempo
y energía recuperable en una línea, suponiendo una eficiencia constante de la máqui-
na de , será igual a:
. Ecuación 17
∑ . ∑ Ecuación 18
Para cada línea o nudo de consumo es posible obtener los pares de puntos caudal-
altura recuperable (Output 2). Suponiendo una máquina con 1 (teórica), la expre-
sión anterior nos proporcionará el máximo valor de energía teórica a recuperar para
87
Bombas funcionando como Turbinas (PAT)
cualquier línea o nudo de consumo del sistema. Para el caso de instalación de máquinas
reales, el valor de la eficiencia y el valor de la altura a recuperar vendrán determinados
por las curvas Q- y Q- , específicas de la máquina o máquinas a instalar.
En la Figura2 se esquematizan los términos energéticos considerados.
ELFRik
EFRik
ELTRik
ELTAik
ETRik
Línea de ELNik
ETik
distribución
zL Nudo de
consumo(toma)
ENik
Qik
zi
z=0
88
Localización, selección y funcionamiento óptimos de PAT
89
Bombas funcionando como Turbinas (PAT)
Funciones objetivo
Se estudian dos funciones objetivo diferentes para la búsqueda de las configuraciones
óptimas. La primera función objetivo es la energía total recuperada en el sistema, y
viene definida según:
Ecuación 21
donde: es la energía total recuperada en kWh/año, según la Ecuación 20.
La instalación de un mayor número de sistemas de recuperación, implica aumentar
la energía recuperada. Por ello, es necesario establecer criterios económicos para poder
seleccionar las soluciones más viables. La segunda función objetivo aplica un criterio
económico para definir los sistemas de recuperación óptimos (Castro, 2006; Nicotra et
al., 2018). Esta segunda función objetivo viene definida según:
Ecuación 22
Ecuación 23
Ecuación 24
Ecuación 25
Ecuación 26
Ecuación 27
Ecuación 28
90
Localización, selección y funcionamiento óptimos de PAT
Generación de lista ordenada de
elementos candidatos
Generación de la configuración Asignación de valores de control del
inicial algoritmo (Ti, Tf, a, L0 y N)
Algoritmo Simulated Annealing
n=n+1
Ubicación óptima para n sistemas
de recuperación
Sí
n<N
No
Fase IV
91
Bombas funcionando como Turbinas (PAT)
Como resultado final se obtienen las configuraciones óptimas para los ′ ′ sistemas
de recuperación (Output 3). Estas configuraciones nos indican que líneas y/o nudos de
consumo son los mejores candidatos a que se instalen sistemas de recuperación.
2.4. Fase IV. Determinación, características y funcionamiento de los sistemas
de recuperación (Optimización II)
La energía final recuperada dependerá fundamentalmente del tipo de máquina instala-
da, del número de máquinas y de las consignas de funcionamiento establecidas. Los
caudales y alturas a recuperar determinarán dentro de las máquinas disponibles, cuáles
de ellas se adaptan a las situaciones de recuperación, y finalmente cuáles de ellas ma-
ximizan las funciones objetivo establecidas. Para ello, es necesario el conocimiento de
las curvas Q-H y Q-, así como otros datos relativos al funcionamiento de las
máquinas.
Debido a la variabilidad de los caudales para las redes de distribución, y con el
objetivo de incrementar el rango de operación del sistema y la garantía de
funcionamiento, es frecuente la instalación de varias máquinas en paralelo, aumentando
de esta manera el rango de caudales y volúmenes turbinados (Pérez-Sánchez et al.,
2018b). El número de máquinas a instalar dependerá de sus curvas de funcionamiento,
así como de su velocidad de rotación.
La instalación de varias máquinas y las distintas posibilidades que ofrecen sus
diferentes modos de regulación, obliga al establecimiento de diferentes consignas de
funcionamiento. En este sentido, es necesario establecer los caudales mínimos y
máximos de funcionamiento para cada una de las máquinas, así como cuál de los
modos de regulación posibles es el óptimo para maximizar las funciones objetivo. La
condición necesaria de funcionamiento de máquinas en paralelo será que el caudal
turbinado y altura recuperada por cada una de las máquinas sea exactamente igual
(Pérez-Sánchez et al., 2017a; 2018b).
La disponibilidad de una extensa gama de máquinas, implica que exista un amplio
abanico de soluciones, pero dentro de estas soluciones existirá una solución que
cumpla con las condiciones impuestas, así como que maximize las funciones objetivo
establecidas. Debido a la gran cantidad de combinaciones posibles es necesario
emplear un método heurístico de búsqueda, que permita trabajar con cualquier
configuración en un sistemas de recuperación (tipo de máquina, número de máquinas y
consignas de funcionamiento) y permita analizar cualquier combinación entre sistemas
de recuperación. Debido a la naturaleza del problema se ha seleccionado de nuevo el
algoritmo ‘Simulated Annealing’, que está indicado para este tipo de problemas, siendo
tratado con anterioridad. El objetivo, es desarrollar una combinación de grupos de
máquinas que maximicen el rendimiento global del sistema de recuperación atendiendo
a las funciones objetivo. Este rendimiento busca que ratio energía recuperable vs
energía teórica recuperable en la red sea máximo, teniendo en cuenta la eficiencia de
las máquinas. A continuación se describen con más detalle algunos aspectos ya
comentados.
92
Localización, selección y funcionamiento óptimos de PAT
Ecuación 31
Ecuación 32
93
Bombas funcionando como Turbinas (PAT)
94
Localización, selección y funcionamiento óptimos de PAT
95
Bombas funcionando como Turbinas (PAT)
Fase IV Definición del problema a resolver:
• Número de sistemas de recuperación y su ubicación (N).
• Posibles máquinas a instalar.
• Número de máquinas a instalar.
• Estrategias de alturas.
• Estrategias de caudales.
• Definición de la función objetivo ( 1 o 2)
Orden de cálculo
Generación de la configuración inicial
Curvas y características de
Lista ordenada de
funcionamiento de
máquinas
máquinas reales
Elección de la máquina
n=n+1
Mejora forma turbinado
Sí n<N
No
Asignación de valores de
Algoritmo Simulated
control del algoritmo
Annealing
(Ti, Tf, a y L0)
Configuración óptima para los N
sistemas de recuperación
FIN
96
Localización, selección y funcionamiento óptimos de PAT
3. Caso de estudio
3.1. Descripción
A continuación, se aplica la metodología descrita a una red de distribución para riego.
La red está ubicada en el término municipal de Vallada (Valencia, España). Este siste-
ma abastece 137 hectáreas, principalmente cítricos, árboles frutales y una pequeña
parte de olivar (Figura 7). La red es abastecida mediante un depósito de hormigón pre-
fabricado de 4000 m3 de capacidad. Las alturas de los nudos de consumo varían entre
340 y 250 metros sobre el nivel del mar. El depósito se encuentra a una altura de
360 metros sobre el nivel del mar, lo que asegura una presión mínima en los nudos de
la red para su correcto funcionamiento. Para los nudos de consumo la presión mínima
de funcionamiento es de 30 mca.
El material de las líneas de distribución es fibrocemento (diámetros entre 250 y
500 mm) y polivinilo de cloruro (diámetros entre 63 y 200 mm). Las tomas de consu-
mo se conectan a la red de distribución mediante tuberías de polietileno (diámetros
entre 32 y 140 mm).
Todas las tomas de consumo tienen instalados caudalímetros para la medición de
los volúmenes consumidos. Se dispone de los registros mensuales de volúmenes desde
el año 2001 y hasta el año 2017.
La demanda base de los nudos de consumo se determina a partir de las característi-
cas de cada una de las parcelas. Se tiene en cuenta el tipo de cultivo y las características
de la instalación de riego (i.e., tipo de gotero y/o aspersor, distancia). El número de
sectores de riego viene determinado por la superficie de las parcelas.
97
98
299.57
70 Nombre línea en red de distribución
0 100 200 300 400 500 10
55
Altura de la metros Nudo con toma de consumo
balsa: 360 303.44
Nudo sin toma de consumo
16
347.82
10
323.61
02
197
1074 254.05 Altura del nudo
10
349.73 1014
10
03
356.88 1001 101
1013 338.89 5
345.99 269.54 Dirección del flujo
04
34
306.15 270.49
10
1
268.76
279.61 31
182
347.1 271.31
05 36 10
27
10 339.63 9 Nudo
364
30 279.5 11 270.01 ea
344.47 279.39 11 252.33 final
11
271.28 Lín
24
3
16
22
11 271.47
341.57 279.52
11
24
251.8 Nudo
103
64
02
29
17
0
21
7
278.57 24281
273.84 254.24 17
6
inicial
279.75
1006
298.66 282.56 28 270.69
28 279.51 178
24
337.81 355 196 Nombre línea=Nombre nudo final
4
286.78 10 285.8 29
112
279.51277.58 254.05
5
10
17
23 25
179
18 303.37 278.31
100
77
265.26
2189
304.47 283.5
7
10
17
18
1104 1025
0
5
17
181
303.43
17
24
289.75 1026 287.61 284.9 286.19 1110 Líneas totales 257
353
35
2
17
2 333.13 1108 10 265.78 222
8
0
280.47 276.22
235
376 1107 1071
174
329.13
10
14 290.1 289.16 286.1 236 281.88 270.07 1094259.13 254.68 Líneas distribución 125
329.06 289.18 283.24
18
309
5 11
35
168
326.88 293.94
6
169
1103 70
10
302.04 270.95 10
324.24 289.99 284.5 264.77
85
9 291.23 Líneas a nudos de consumo 132
30
301.5 19
2
6 324.99 296.21 110
100
7
10
23
8
271.32
21
296.34
10
1 10
19
7 322.72 09 Nudos sin toma de consumo 125
69
322.23
10
1082 10 324.86 322.69
319.42 22 295.7
10
298.13
11
8
8
30
20
16
295.73
1
29
10
14
9 321.01 321.69 11 111913 21 Nudos con toma de consumo 132
2 297
106
101 319.38 10 292.15
20
35
1079
29
319.83 321.3 272.51 268.35 264.99
20
9
167
12 319.61 295.52 295.46 290.13 300
109
295.15
15
2
10
318.92 318.3 2 270.81 267.8 Longitud Diámetro Longitud Diámetro Longitud Diámetro
3 7
33
297
10 6 Línea Línea Línea
10
Longitud Diámetro 294.29 291.63 269.34 266.03 (m) (mm) (m) (mm) (m) (mm)
10
291.14
29
93
Línea
6
(m) (mm) 4 12 1033 57.505 90 1065 134.571 200 1097 18.89 53.6
103 349 289.37 11 266.48
1037
1034 288.057 350 1066 218.961 76.6 1098 94.49 53.6
Bombas funcionando como Turbinas (PAT)
10
35
10
1003 75.974 500 1036 326.444 57 1068 165.077 200 1100 45.029 136.2
65
277.24 64
277.89 34 10
3
1004 87.775 350 281.21 4 261.1 1037 37.125 57 1069 61.155 53.6 1101 14.723 53.6
275.29 345
111
1038
303
1005 145.767 300 1038 95.165 300 1070 190.187 144.6 1102 417.012 350
10
271.43 260.11
63
34
2
339 278.8 78 258.09
281.84 10 2
29
1006 148.435 250 285.77 1096 269.8 1039 124.414 350 1071 135.401 126.6 1103 156.813 126.6
2
275.24 274 106
307
279.05
103
1007 126.981 126.6 274.57 272.28
9
104 1040 127.414 57 1072 254.502 350 1104 56.465 99.4
343
2
106
4
106
30
1
277.1
30
1008 115.621 126.6
6
270.83 269.61
1041 45.05 350 1073 30.015 350 1107 83.903 126.6
277.32 1040
1041
60
30
1009 102.902 67.8 340 1095
5
273.49 1042 238.573 99.4 1074 30 500 1108 86.357 126.6
10
273.38 277.96 293 256.25
1010 39.45 57 274.99
32
258.14
2
1043 81.152 350 1075 32.441 93.6 1109 139.474 126.6
294
1043
5 267.27
30
1011 46.206 113 32 273.27 256.91
8
324 266.74 1044 1044 153.495 119.2 1076 37.673 119.2 1110 46.497 126.6
271.49 47 341 263.1
10
310
1053
10 267.88 268.73 23
328
1012 146.142 81.4 11
11
46
3
1045 51.831 119.2 1077 35.412 119.2 1111 41.185 99.4
327
1087 1052 1056
10
323
329
1077 313 320
54
261.34
1086
1013 132.547 400 99
261.21 1046 112.325 119.2 1078 37.175 300 1112 159.876 200
326
269.57 1045266.09 10 264.43
1100
311
1014 345.47 350 319
1058
259.18
10
267.63 0 1048 049 264.1 4 318
1047 70.229 53.6 1079 54.243 113 1113 62.632 76.6
57
259.03 31
109
1015 303.658 350 33 334 1 336 258.19
10
0
266.24 263.69
315
89
263.69 1048 23.984 136.2 1080 27.133 119.2 1114 97.126 126.6
1080
1016 311.706 250 1
1084
338 335
262.46 32 258.83255.54 1049 100.321 136.2 1081 27.221 126.6 1116 222.174 99.4
1017 54.539 350 312 386 391 254.36
262.671075 331 262.12 259.84 257.57
105
01
389
1050 47.984 93.6 1082 84.923 67.8 1117 18.272 53.6
317
254.16
5
11
262.49 337
1050
1018 131.123 350 332 388
262.24
33
1051 254.33 1051 120.881 93.6 1083 167.094 53.6 1118 41.518 93.6
3
254.6 8
390
393
1019 120.256 350
394
1083 109
387
316
261.73 1097 257.81 1052 162.548 153.2 1084 88.287 93.6 1119 31.019 76.6
1118
1020 80.96 350 399 254.11
260.55 255.84
10
396
395
397
11 1053 55.68 350 1085 15.319 126.6 1120 18.124 40.8
1021 144.903 250 27
256.33 256.1 1054 62.235 350 1086 25.471 136.2 1122 89.136 57
1022 50.375 200 398
92 254.63 1055 91.913 53.6 1087 82.125 119.2 1123 25.97 136.2
1023 274.188 162.8 255.69 3
1024 60.683 126.6 1056 228.898 300 1088 14.081 53.6 1124 7.652 93.6
1025 53.351 99.4 1057 18.337 136.2 1089 39.642 93.6 1125 92.803 99.4
1026 73.585 126.6 1058 48.171 53.6 1090 169.181 300 1126 34.368 93.6
1027 61.165 57 1059 96.91 119.2 1091 24.132 53.6 1127 131.416 63.8
412.5
400
296.0
300 248.5
185.2 190.1
200
134.4 133.5
101.8 99.3
100
4. Resultados
4.1. Caudales y presiones en la red (Fase I)
En la Figura 9, se observan los resultados de caudales para la línea de cabecera (línea
1074) a lo largo de todo el año.
El caudal queda discretizado mediante un intervalo horario (8760 registros). El
caudal para sistemas de riego presenta una gran variabilidad a lo largo del año, ya que
queda condicionado a las particularidades agronómicas de los cultivos. Para la línea
1074 el caudal máximo fue de 88.17 l/s (el 29/08 de 09:00-10:00). Para el mes de
agosto, el caudal medio mensual es máximo y con un valor de 34.50 l/s, mientras que
el menor valor se encuentra en el mes de enero con 8.56 l/s. El caudal medio a lo largo
del año presenta un valor de 19.77 l/s.
99
Bombas funcionando como Turbinas (PAT)
100 500000
90 450000
Volumen acumulado (m3)
80 400000
70 350000
Caudal (l/s)
60 300000
50 250000
40 200000
30 150000
20 100000
10 50000
0 0
dic
ene
jul
jun
mar
may
feb
nov
oct
ago
abr
sep
Caudal línea (l/s) Caudal medio mensual (l/s)
Caudal medio anual (l/s) Volumen acumulado(m3)
Figura 9. Caudales para la línea de cabecera a lo largo de todo el año (línea 1074)
En la Figura 10 se observan los resultados a lo largo del día 02/08 para la toma de
consumo 340. Para esta toma se aprecia como el riego se produce durante 3 horas,
entre las 16:00 y 19:00 horas, produciéndose la menor presión en la última hora de
riego, y con un valor de 80.22 mca, siendo la mayor presión durante el día de
82.68 mca.
6.0 83.0
5.44
82.5
5.0
82.0
4.0 81.5
Caudal (l/s)
Presión ( mca)
81.0
3.0
80.5
2.0 80.0
79.5
1.0
79.0
0.0 78.5
10:00
11:00
12:00
13:00
14:00
15:00
16:00
17:00
18:00
19:00
20:00
21:00
22:00
23:00
0:00
1:00
2:00
3:00
4:00
5:00
6:00
7:00
8:00
9:00
Caudal (l/s) Presión (mca)
100
Localización, selección y funcionamiento óptimos de PAT
Tabla 1. Resultados de las funciones (ETR), PSR y (ETR/PSR) para las tomas
Toma PSR Toma PSR Toma PSR
355 2.94 5.61 0.53 317 0.42 21.99 0.02 394 0.13 23.70 0.01
312 2.77 14.17 0.20 316 0.39 63.92 0.01 244 0.12 135.19 0.01
222 2.10 18.52 0.12 391 0.39 18.79 0.03 398 0.10 60.59 0.01
167 1.95 11.52 0.17 176 0.39 29.77 0.02 334 0.10 24.77 0.01
311 1.85 20.15 0.10 321 0.38 25.16 0.02 328 0.10 53.51 0.01
333 1.29 13.59 0.10 10 0.38 44.79 0.01 22 0.10 20.44 0.01
340 1.26 16.92 0.08 326 0.36 15.05 0.03 175 0.10 44.73 0.01
1 1.25 4.11 0.31 305 0.33 19.23 0.02 393 0.10 33.03 0.01
113 1.21 21.61 0.06 314 0.32 23.21 0.02 350 0.10 20.74 0.01
323 1.07 19.47 0.06 21 0.32 14.97 0.03 293 0.10 111.03 0.01
364 0.95 16.78 0.06 336 0.30 74.95 0.01 331 0.10 47.87 0.01
300 0.91 26.95 0.04 235 0.30 65.36 0.01 5 0.09 9.79 0.01
236 0.87 21.82 0.04 349 0.27 10.94 0.03 397 0.09 37.69 0.01
292 0.85 25.99 0.04 27 0.27 30.69 0.01 6 0.08 29.76 0.01
35 0.81 17.57 0.05 13 0.26 12.76 0.03 9 0.08 25.21 0.01
246 0.75 20.24 0.04 20 0.25 56.57 0.01 322 0.08 116.39 0.01
18 0.75 18.69 0.04 389 0.24 10.55 0.03 15 0.08 25.67 0.01
332 0.73 7.59 0.10 399 0.23 43.58 0.01 181 0.07 56.66 0.01
319 0.72 37.08 0.02 179 0.22 20.29 0.02 395 0.06 52.69 0.01
169 0.69 10.67 0.07 306 0.21 32.22 0.01 296 0.06 54.95 0.01
170 0.67 11.96 0.06 356 0.21 53.24 0.01 16 0.06 15.53 0.01
324 0.62 27.48 0.03 329 0.20 61.17 0.01 14 0.06 13.00 0.01
2 0.62 9.20 0.07 343 0.19 27.17 0.01 172 0.06 165.38 0.01
386 0.61 7.41 0.09 55 0.18 46.33 0.01 12 0.04 52.21 0.01
23 0.57 21.20 0.03 180 0.18 55.00 0.01 339 0.04 33.01 0.01
294 0.57 56.16 0.02 337 0.17 17.48 0.01 299 0.03 137.85 0.01
345 0.55 34.69 0.02 313 0.17 40.37 0.01 25 0.03 40.35 0.01
173 0.53 23.24 0.03 196 0.17 63.52 0.01 19 0.02 192.28 0.01
310 0.53 22.77 0.03 245 0.16 60.89 0.01 335 0.02 138.65 0.01
318 0.53 19.05 0.03 320 0.16 65.97 0.01 325 0.02 107.93 0.01
308 0.51 35.21 0.02 174 0.16 32.63 0.01 24 0.02 743.75 0.01
392 0.5 41.53 0.02 327 0.15 60.58 0.01 3 0.01 26.47 0.01
298 0.47 55.56 0.01 338 0.15 13.67 0.02 11 0.01 74.05 0.01
315 0.45 23.14 0.02 34 0.14 104.7 0.01
342 0.43 27.01 0.02 390 0.14 37.13 0.01
101
Bombas funcionando como Turbinas (PAT)
1020 22.00 6.77 3.26 1071 3.21 15.84 0.21 2164 0.90 25.47 0.04
1019 21.05 6.70 3.21 1070 3.03 16.75 0.19 1058 0.76 21.97 0.04
1102 20.55 6.61 3.12 1084 3.02 17.55 0.18 1040 0.71 21.04 0.04
1043 19.76 9.15 2.16 1075 2.86 18.51 0.16 1066 0.66 48.19 0.02
1015 19.52 6.96 2.81 1050 2.60 20.33 0.13 1027 0.61 27.56 0.03
1038 19.31 8.50 2.28 1051 2.60 20.33 0.13 1127 0.60 41.54 0.02
1018 18.92 6.64 2.85 1076 2.60 17.81 0.15 1083 0.57 22.99 0.03
1039 18.41 8.43 2.19 1108 2.42 17.52 0.14 1069 0.47 55.56 0.01
1017 18.00 6.65 2.71 1107 2.36 17.90 0.14 1104 0.46 29.20 0.02
1041 17.81 8.49 2.10 1044 2.20 15.94 0.14 1091 0.45 28.86 0.02
1053 17.71 10.21 1.74 1087 2.20 15.94 0.14 1120 0.41 15.38 0.03
1034 16.81 7.75 2.18 1077 1.95 13.15 0.15 1097 0.40 29.82 0.02
1035 16.81 7.75 2.18 1080 1.89 12.25 0.16 1098 0.37 19.90 0.02
1032 15.76 7.95 1.99 1045 1.85 13.31 0.14 1093 0.36 12.13 0.03
1054 13.01 9.90 1.32 1109 1.66 18.12 0.1 1036 0.36 12.16 0.03
1056 9.35 11.22 0.84 1026 1.66 18.11 0.1 1079 0.33 11.48 0.03
1021 6.46 10.03 0.65 1055 1.64 10.72 0.16 1008 0.27 10.39 0.03
1064 5.60 15.12 0.37 1114 1.60 18.75 0.09 1030 0.27 30.69 0.01
1060 5.20 16.27 0.32 1046 1.53 16.14 0.10 1122 0.27 30.69 0.01
1062 5.20 16.27 0.32 1113 1.48 32.84 0.05 1025 0.25 55.28 0.01
1063 5.20 16.27 0.32 1065 1.48 32.89 0.05 1101 0.24 10.55 0.03
1014 5.17 6.57 0.79 1047 1.44 14.97 0.10 1119 0.22 12.03 0.02
1022 5.12 12.79 0.41 1059 1.41 14.65 0.10 1016 0.18 46.34 0.01
1052 5.11 17.65 0.29 2189 1.39 25.14 0.06 1125 0.16 63.55 0.01
1099 5.11 17.65 0.29 1110 1.29 23.19 0.06 1011 0.15 21.32 0.01
1090 4.76 12.10 0.4 1029 1.22 21.20 0.06 1031 0.14 104.7 0.01
1086 4.73 14.10 0.34 1123 1.21 21.61 0.06 1111 0.14 104.7 0.01
1112 4.32 16.55 0.27 1116 1.09 26.09 0.05 1037 0.10 20.74 0.01
1067 4.27 16.75 0.26 1095 1.08 16.95 0.07 1009 0.06 15.72 0.01
1103 3.47 13.42 0.26 1033 1.05 24.24 0.05 1082 0.06 15.72 0.01
1100 3.44 15.86 0.22 1092 1.05 24.24 0.05 1012 0.06 15.53 0.01
1023 3.43 13.60 0.26 1024 1.04 23.44 0.05 1085 0.05 6.29 0.01
1068 3.42 19.01 0.18 1124 1.00 23.16 0.05 1096 0.04 33.01 0.01
1089 3.40 18.19 0.19 1118 0.96 36.22 0.03 1010 0.03 32.70 0.01
1094 3.33 15.23 0.22 1078 0.94 31.03 0.04 1028 0.02 743.7 0.01
1057 3.33 16.45 0.21 1088 0.93 14.84 0.07 1081 0.01 6.60 0.01
1049 3.33 18.61 0.18 1126 0.92 29.52 0.04
102
Localización, selección y funcionamiento óptimos de PAT
Para la toma de consumo 355 se obtiene la mayor energía teórica recuperable con
un valor de 2.94 MWh/año, siendo también esta toma la que presenta el mayor valor
para el ratio entre la energía teórica recuperable y el PSR (0.53 MWh/año2). La energía
total recuperable para todas las tomas es de 44.9 MWh/año siendo necesario la
colocación de 103 sistemas de recuperación a lo largo de toda la red (uno por cada
toma de consumo).
La línea 1020 es la que presenta una mayor energía teórica recuperable con un valor
de 22.00 MWh/año y un valor máximo de 3.26 MWh/año2 para el ratio entre la energía
teórica recuperable y el PSR.
Figura 11. (a) Puntos Q-H de recuperación para la toma 355. (b) Potencia instantánea y
energía acumulada a lo largo del año para la toma 355
En la Figura 11a se observan los pares de puntos Q-H para la toma 355. La Figura
11b muestra la potencia teórica instantánea (suponiendo =1) y energía acumulada a lo
largo de un año para la misma toma. Para esta toma, el caudal se mantiene constante e
igual 7.42 l/s, mientras que la altura recuperable se distribuye entre un mínimo de
26.96 m y 28.84 m, siendo la media de 28.58 m. La potencia instantánea se mantiene
prácticamente constante con un valor de 2.08 kW.
La Figura 12a muestra los pares de puntos Q-H para la línea 1020, mientras que la
Figura 12b describe la potencia teórica instantánea (suponiendo =1) y energía teórica
103
Bombas funcionando como Turbinas (PAT)
acumulada a lo largo de un año. Para esta línea, el caudal presenta una gran
variabilidad, siendo el máximo de 63.47 l/s. La altura recuperable se distribuye entre un
mínimo de 27.51 m y 33.99 m, siendo la media de 33.70 m. La potencia instantánea
presenta un máximo de 19.02 kW y una media de 2.5 kW.
Figura 12. (a) Puntos Q-H de recuperación para la línea 1020. (b) Potencia instantánea y
energía acumulada a lo largo del año para la línea 1020
104
Localización, selección y funcionamiento óptimos de PAT
45000 5000
40000 4500
35000 4000
3500
30000
2 (Kw∙h/año2)
1 (Kw∙h/año)
3000
25000
2500
20000
2000
15000
1500
10000 1000
5000 500
0 0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140
Número de sistemas de recuperación instalados (N)
Energía teórica recuperable E/PSR
105
Bombas funcionando como Turbinas (PAT)
Tabla 3. Resultados de ubicación de los sistemas de recuperación para las líneas de distribución
según las funciones (ETR), PSR y (ETR/PSR)
Optimización función objetivo (Energía teórica recuperable)
Ti = 100; Tf =0.01; a =0.9; L = 10
PSR (
N Líneas
(kWh/año) (año) (kWh/año2)
1 1020 21993.20 6.76 3252.13
2 1043 1102 30640.51 7.59 4037.05
3 1020 1043 1102 32707.58 7.41 4413.87
4 1020 1043 1090 1102 33819.62 7.82 4323.34
5 1020 1043 1086 1090 1102 34862.68 8.18 4263.70
6 1020 1043 1052 1086 1090 1102 35755.43 8.66 4128.16
7 1020 1043 1052 1086 1090 1094 1102 36433.04 9.04 4030.58
8 1020 1043 1086 1090 1094 1099 1102 1108 37104.49 9.16 4052.65
9 1022 1038 1043 1071 1086 1090 1099 1102 1116 37582.87 9.62 3908.24
10 1020 1038 1043 1086 1090 1094 1099 1102 1103 1116 38213.31 9.25 4132.86
Optimización función objetivo (Ratio entre Energía teórica recuperable y PSR)
Ti= 100; Tf=0.01; a=0.9; L= 10
PSR (
N Líneas
(kWh/año) (año) (kWh/año2)
1 1020 21993.20 6.76 3252.13
2 1043 1102 30640.51 7.59 4037.05
3 1020 1043 1102 32707.58 7.41 4413.87
4 1020 1038 1043 1102 33364.87 7.36 4531.66
5 1019 1038 1043 1102 1103 33950.51 7.48 4538.22
6 1020 1038 1043 1055 1102 1103 33499.61 7.54 4444.95
7 1020 1038 1043 1055 1079 1102 1103 34729.84 7.60 4570.26
8 1019 1038 1053 1055 1077 1102 1103 1119 34623.73 7.72 4485.55
9 1014 1019 1038 1039 1043 1079 1102 1103 1108 34573.43 7.56 4574.42
10 1011 1019 1038 1039 1043 1055 1079 1102 1103 1108 34557.58 7.70 4489.52
106
Localización, selección y funcionamiento óptimos de PAT
han caracterizado sus curvas Q-H y Q-, así como la zona de operación mediante los
caudales máximos y mínimos de la máquina. Para todas las máquinas se han
considerado un rango de velocidades de αmin=0.5 y αmax=1.5. En la Figura 14 se
observan los puntos de funcionamiento BEP para cada máquina, con indicación del
número específico e identificativo de la máquina en la base de datos utilizada.
225 a
11.64(11)
200
175
6.2(1)
150
HBEP_t(m)
125
100
10.3(6)
22.71(33)
75 17.38(20)
17.95(22) 41.83(64)
15.63(15) 29.71(47)
9.3(5) 19.41(26) 21.33(30)
50 8.82(4) 18.6(25) 24.47(34)
20.51(28) 30.55(50)
27.06(42) 35.01(53) 38.21(60)
34.86(52) 31.9(51)
25 35.13(54)
35.7(55)
47.5(65) 70(75)
67.56(73)
38.17(59)
52(67) 62.63(71)
47.75(66) 65.35(72)
0
0 50 100 150 200 250 300
QBEP_t(l/s)
50 b
21.4(31)
45 12.68(13)
15.58(14)
40 15.96(17)
8.43(3) 17.09(19) 17.8(21)
35 12.6(12)
11.21(9)
8.35(2) 52.6(68)
30 10.87(8) 40.85(62)
HBEP_t(m)
16.01(18)
18.5(24)
28.2(44)
25
11.57(10) 21.65(32) 28.1(43)
10.84(7) 26.91(41)
20 29.82(48) 35.91(56)
15.71(16) 25.47(36)
15 18.21(23) 30.1(49) 41.1(63)
25.24(35) 26.03(38) 28.38(45)
26.37(40) 36.84(57)
25.75(37)
10 21.09(29) 38.1(58)
38.49(61) 26.28(39)
57.6(70)
28.55(46) 68.47(74)
5 19.63(27)
54.69(69)
0
0 5 10 15 20 25 30 35 40 45 50
QBEP_t(l/s)
Figura 14. (a) Punto BEP para máquinas reales utilizadas (b) Detalle para la zona
<50 y <50. Se indica el número específico como turbina y entre parénte-
sis el identificativo de la máquina
107
Bombas funcionando como Turbinas (PAT)
El número de máquinas fijado en este caso de estudio es de tres, para todas las
líneas de las configuraciones analizadas. Para el diseño del sistema de recuperación se
ha seguido el esquema de la Figura 15.
108
Localización, selección y funcionamiento óptimos de PAT
determina la(s) máquina(s) a instalar en cada línea, así como la estrategia de alturas
óptima (curva de funcionamiento) determinando el mejor diseño en los sistemas de
recuperación para optimizar la función objetivo establecida, adjuntadas en la Tabla 4.
Tabla 4. Resultados para las configuraciones óptimas respecto de Ψ (fase III), optimizadas en
función de Ψ (energía real recuperable)
PSR
N ηt
Línea 1020
nst 16.01(18)
D 180
1 Ns 2000 9521.82 4.29 2221.29 0.43
ER 9521.82
PSR 4.29
ER/PSR 2221.29
C 1
1043 1102
Línea
nst 25.24(35) 10.84(7)
D 180 250
Ns 1485 1000
2 10612.54 4.84 2192.15 0.35
ER 1629.97 8982.56
PSR 7.73 4.32
ER/PSR 210.73 2081.21
Curva α=1 α=1
Línea 1020 1043 1102
nst 38.49(61) 38.49(61) 10.84(7)
D 142 142 250
Ns 1520 1520 1000
3 11781.48 4.97 2368.31 0.36
ER 2148.93 649.99 8992.57
PSR 6.27 9.78 4.31
ER/PSR 342.59 66.42 2081.21
Curva BEH BEH α=1
Línea 1020 1043 1090 1102
nst 38.49(61) 38.49(61) 16.01(18) 10.84(7)
D 142 142 180 250
Ns 1520 1520 2000 1000
4 12287.44 5.47 2248.2 0.36
ER 2148.93 649.98 505.96 8982.56
PSR 6.27 9.79 16.89 4.32
ER/PSR 342.59 66.43 29.95 2081.21
Curva BEH BEH BEH α=1
Línea 1020 1043 1086 1090 1102
nst 38.49(61) 38.49(61) 10.87(8) 16.01(18) 10.84(7)
D 142 142 101.4 180 250
Ns 1520 1520 2500 2000 1000
5 13057.31 5.97 2186.66 0.37
ER 2148.93 649.98 769.87 505.99 8982.56
PSR 6.27 9.78 14.05 16.89 4.32
ER/PSR 342.59 66.43 54.81 29.95 2081.21
Curva BEH BEH α=1 BEH α=1
109
Bombas funcionando como Turbinas (PAT)
Tabla 5. Resultados para las configuraciones óptimas respecto de (fase III), optimizadas en
función de (ratio entre energía real recuperable y PSR)
N PSR ηt
Línea 1020
nst 11.57(10)
D 176
1 Ns 1520 7632.69 2.98 2557.58 0.35
ER 7632.69
PSR 2.98
ER/PSR 2557.58
CLínea 1
1043 1102
nst 11.57(10) 11.57(10)
D 176 176
Ns 1520 1520
2 8297.45 2.58 3210.86 0.27
ER 1331.19 6966.26
PSR 3.27 2.45
ER/PSR 406.81 2840.25
Curva BEH α=1
Línea 1020 1043 1102
nst 19.63(27) 19.63(27) 11.57(10)
D 139 139 176
Ns 1020 1020 1520
3 8877.14 2.47 3587.23 0.27
ER 1363.43 547.44 6966.26
PSR 2.35 3.06 2.45
ER/PSR 579.96 178.78 2840.25
Curva α=1 BEH α=1
Línea 1020 1038 1043 1102
nst 19.63(27) 19.63(27) -------- 11.57(10)
D 139 139 -------- 176
Ns 1020 1020 -------- 1520
4 8938.07 2.46 3636.65 0.27
ER 1363.43 608.38 0 6966.26
PSR 2.35 2.76 0 2.45
ER/PSR 579.96 220.79 0 2840.25
Curva α=1 BEH ------- α=1
Línea 1019 1038 1043 1102 1103
nst 19.63(27) 19.63(27) -------- 11.57(10) 10.87(8)
D 139 139 -------- 176 101.4
Ns 1020 1020 -------- 1520 2500
5 9159.02 2.53 3615.24 0.27
ER 1363.43 608.38 0 6966.26 220.94
PSR 2.35 2.76 0 2.45 5.59
ER/PSR 579.95 220.79 0 2840.25 39.49
Curva α=1 BEH ------- α=1 BEH
110
Localización, selección y funcionamiento óptimos de PAT
111
Bombas funcionando como Turbinas (PAT)
Tabla 6. Características principales y valores para las máquinas en paralelo para la solución
analizada. Valores para la línea
Línea 1102 1020 1043
nst 11.57(10) 19.63(27) 19.63(27)
ns (rpm) 1520 1020 1020
D (mm) 176 139 139
Qmin (l/s) 2.10 1.59 1.59
Qmax (l/s) 8.81 5.31 5.31
QBEP (l/s) 5.65 3.64 3.64
HBEP (l/s) 21.58 4.50 4.50
ηBEP 0.614 0.611 0.611
PBEP(kW) 0.735 0.098 0.098
Línea 1102 1020 1043
n 3 3 3
ER(kW·h/año) 6966.26 1363.43 547.45
Vt (m3) 203999.41(0.68) 159324.99(0.63) 96535.60(0.68)
Qmin (l/s) 2.10 1.59..5 1.59
Qmax (l/s) 16.91 12.75 9.92
HRmin (mca) 4.53 1.08 1.06
HRmax (mca) 21.46 5.75 3.74
Pmax (kW) 2.19 0.41 0.21
Rango 0.341 0.308 0.289
caudales
Curva
α=1 α=1 BEH
funcionamiento
1 2 3 1 2 3 1 2 3
% ER (kWh/año) 44.53 32.18 23.29 41.94 33.13 24.93 44.70 33.21 54.97
% Vt (m3) 45.76 31.72 22.52 43.09 32.68 24.22 46.18 32.54 21.28
Horas 4945 3266 2289 4860 3506 2559 4007 2695 1745
Qmin (l/s) 2.12 4.40 4.87 1.59 3.17 3.53 1.59 2.65 2.90
Qmax (l/s) 5.64 5.64 5.64 4.25 4.25 4.25 3.31 3.31 3.31
HRmin (mca) 4.53 15.01 17.20 1.08 3.75 4.31 1.06 2.41 2.88
HRmax (mca) 21.46 21.46 21.46 5.75 5.75 5.75 3.74 3.74 3.74
Pmax(kW) 0.73 0.73 0.73 0.14 0.14 0.14 0.07 0.07 0.07
Línea 1102 1020 1043
V (m3) 300427.46 251219.42 141189.20
Qmin (l/s) 0.002 0.002 0.002
Qmax (l/s) 78.21 63.37 46.81
Htmin (mca) 22.61 6.05 3.89
Htmax (mca) 25.52 32.99 52.75
112
Localización, selección y funcionamiento óptimos de PAT
60
Curvas Q‐H para ns=11.57 rpm ‐m,Kw‐ , D=176 mm y N=1520 rpm. a
Número de máquinas instaladas=3
αmax=1.5
Hmax=55.87 m c.a.
50
Runaway
40
BEH
α=1
Altura (m c.a.)
30
20 Curvas Q‐H
BPC QBEP=5.65 l/s
HBEP=21.58 m c.a. 0.741 2.207 10.33
PBEP=0.735 kW . 0.691 1.24 3.50
. 0.837 3.15 18.66
10 αmin=0.5 0.685
Hmin=3.95 m c.a. 0.725
2.205
0
0 1 2 3 4 5 6 7 8 9 10 11 12
Qmin=2.10 l/s Caudal (l/s) Qmax=8.81 l/s
0.90
0.80
0.70 QBEP=5.65 l/s
ηBEP=0.614
0.60
BEH
Rendimiento
0.50
0.30
Curvas Q‐η
η 0.032 0.367 0.424
0.20
η . 0.077 0.526 0.366
η . 0.016 0.228 0.333
0.10 η 0.012 0.136Q 0.230
η 0.018 0.221 0.114
0.00
0 1 2 3 4 5 6 7 8 9 10 11 12
Qmin=2.10 l/s Caudal (l/s) Qmax=8.81 l/s
Figura 16. Curvas Q-H (a) y Q-η (b) para la máquina nst=11.57 (m, kW), instalada
en la línea 1102 (tres máquinas en paralelo)
113
Bombas funcionando como Turbinas (PAT)
12
Curvas Q‐H para ns=19.63 rpm ‐m,Kw‐ , D=139 mm y N=1020 rpm. a
11
Hmax=10.16 m c.a.
10 αmax=1.5
Runaway
9
8
BEH
7
Altura (m c.a.)
6 α=1
BPC
5
4 Curvas Q‐H
0.422 1.281 5.60
3 0.393 0.719 1.20
.
. 0.476 1.83 6.43
2 0.342
αmin=0.5
Hmin=1.05 m c.a.
0.367
1
0.507
0
0 1 2 3 4 5 6 7 8
Qmin=1.59 l/s Caudal (l/s) Qmax=5.31 l/s
1.00
Curvas Q‐η para ns=19.63 rpm ‐m,Kw‐ , D=139 mm y N=1020 rpm. b
0.90
0.80
0.70 QBEP=3.64 l/s
ηBEP=0.611
0.60
BEH
Rendimiento
0.30
Curvas Q‐η
η 0.111 0.808 0.861
0.20
η . 0.264 1.16 0.743
η . 0.054 0.502 0.677
0.10 η 0.036 0.246Q 0.174
η 0.044 0.342 0.116
0.00
0 1 2 3 4 5 6 7 8
Qmin=1.59 l/s Caudal (l/s) Qmax=5.31 l/s
Figura 17. Curvas Q-H (a) y Q-η (b) para la máquina nst=19.63 (m, kW), instalada
en la línea 1020 y 1043 (tres máquinas en paralelo)
114
Localización, selección y funcionamiento óptimos de PAT
55 0.225
50
0.200
45
0.175
40
0.150
Potencia (Kw)
35
Altura (m c.a)
30 0.125
25 0.100
20
0.075
15
0.050
10
5 0.025
0 0.000
0 5 10 15 20 25 30 35 40 45 50
Caudal línea (l/s)
Altura recuperada Q‐H con sistemas de recuperación aguas arriba
Q‐H sin sistemas de recuperación aguas arriba Potencia
Figura 18. Altura recuperada y potencia para el sistema de recuperación en la línea 1043
La Figura 19 corresponde al sistema de recuperación en la línea 1102. Se aprecia la
variación de eficiencia y relación de velocidades de rotación (α) en las máquinas. Para
el coeficiente α se observa como, a partir de un caudal de 2.10 l/s (caudal mínimo de
turbina), empieza a aumentar hasta un caudal de 3.20 l/s. Este caudal corresponde con
el caudal mínimo para la curva α 1. A partir de este caudal las máquinas en
funcionamiento trabajan siempre a α constante. La variación de velocidad en los tramos
de caudales inferiores a los mínimos permite el aprovechamiento de un mayor rango de
caudales, y por lo tanto mayores volúmenes turbinados. Para ello, es necesario la
variación del valor de α en función del caudal por línea. Respecto a los rendimientos, a
partir de una caudal por línea de 16.95 l/s las máquinas trabajan a su máximo
rendimiento ya que el caudal máximo por turbina coincide con el caudal en BEP.
115
Bombas funcionando como Turbinas (PAT)
24 1.6
22
αmax 1.5
1.4
20
1.3
18 1.2
1.1
16
Altura recuperada (m c.a.)
Rendimiento ‐
14
0.9
12 0.8
0.7
10
ηmax 0.6
8
αmin 0.5
6 0.4
0.3
4
0.2
2
0.1
0 0
0 10 20 30 40 50 60 70 80 90
Caudal línea (l/s)
Altura recuperada Rendimiento(η) α
70
35
60
Caudal (l/s)
30
50
25
20 40
15 30
10 20
5 10
0 0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70
Caudal línea (l/s)
Caudal No turbinado Caudal turbinado % Caudal No turbinado % Caudal turbinado
116
Localización, selección y funcionamiento óptimos de PAT
Tabla 7. Funcionamiento del sistema de recuperación para la línea 1020 en función del
valor del caudal por línea
Valor de caudal Descripción del Posición de las válvulas Nº de máquinas
instantáneo (l/s) funcionamiento de control funcionando
VC1, VC2, VC3 y
Todo el caudal es derivado VC4 cerradas
<1.59 aguas abajo a consumo y/o a VC5 abierta 0
otro sistema de recuperación VR1 no es necesario
que actúe.
VC1 y VC4 abiertas
Todo el caudal instantáneo es VC2, VC3 y VC5
1.59 4.25 1
turbinado mediante la PAT1. cerrada
VR1 cerrada
Como se observa para valores inferiores a 1.59 l/s todo el caudal debe ser derivado
y no turbinado. A partir de este valor la puesta o parada de las máquinas se realiza con
el criterio de obtener la máxima potencia instantánea. Este criterio finalmente se
traduce en la parada o puesta en marcha de las diferentes máquinas en función del valor
del caudal por línea. Cuando las máquinas llegan a su máximo caudal de turbinado
(4.25 l/s), se mantiene este caudal pero obliga al sistema a derivar caudal sin turbinar.
En el paso de funcionamiento de la PAT1 a la PAT2, este caudal derivado llega a ser
hasta un 32.35 % del total. A partir de 12.75 l/s el caudal no turbinado empieza a
117
Bombas funcionando como Turbinas (PAT)
aumentar y llega hasta un 79.87% del total (para el caudal máximo por línea). La
válvula reductora de presión, trabaja siempre reduciendo 5.75 m para diferentes
caudales (línea 1020), lo que implica que el valor de varíe con el tiempo. Los
valores de son mostrados en la Figura 21.
300
250
KV(m3/h/bar0.5)
200
150
100
50
0
1
245
489
733
977
1221
1465
1709
1953
2197
2441
2685
2929
3173
3417
3661
3905
4149
4393
4637
4881
5125
5369
5613
5857
6101
6345
6589
6833
7077
7321
7565
7809
8053
8297
8541
Horas
1400 100.00
1196 90.00
1200
80.00
Número de horas
1000 70.00
% Acumulado
800 60.00
50.00
600 40.00
Número de horas
400 % Acumulado 30.00
298
251
216196 20.00
200 143125121
107 83
65 53 45 37
29 24 18 15 13 5
10.00
0 5 1 2 0 1
0 0.00
100
110
120
130
140
150
160
170
180
190
200
210
220
230
240
250
10
20
30
40
50
60
70
80
90
0
Kv(m3/h/bar0.5 )
118
Localización, selección y funcionamiento óptimos de PAT
119
Bombas funcionando como Turbinas (PAT)
Conclusiones
El presente capítulo aborda el desarrollo de una metodología de optimización que
permite, conocida la demanda a lo largo del tiempo de los nudos, la determinación de
los caudales circulantes y a partir de estos, la determinación de las alturas teóricas
recuperables. Conocidos los parametros de caudal y alturas recuperables, la
metodología optimiza la localización de grupos de sistemas de recuperación instalados
en serie que permitan maximizar la energía recuperable. La optimización incluye la
selección de máquinas y optimización en la regulación de acuerdo a la estrategia
establecida y que fue desarrollada en el capítulo anterior. Finalmente, se ha aplicado a
un caso de estudio, concretamente a una red de riego, mostrando el gran potencial de
esta herramienta para los gestores de sistemas de distribución para la mejora de
eficiencia energética en las redes presurizadas.
120
Análisis de viabilidad económica a través del
uso de sistemas con PAT
Pérez-Sánchez, Ma,, Sánchez-Romero, FJb, Ramos, HMc, .López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
Un análisis completo de la instalación de sistemas de bombas funcionando como turbinas
implica conocimiento de aspectos económicos. Es evidente que los beneficios de la recupe-
ración energética son energéticos, ambientales y económicos. Sin embargo, no es sencillo
cuantificar los mismos para proponer la viabilidad del sistema conjunto. El presente capítulo
aborda estos aspectos, haciendo una propuesta sobre cómo estimar esta viabilidad de forma
sistemática y poniéndolo de manifiesto en un caso de estudio concreto, mostrando las capa-
cidades del método para hacer cuantificaciones y comparaciones de diversas configuraciones
y metodologías de estimación económicas.
Abstract
A complete analysis of PAT facility implies the knowledge of the economic aspects. It is
clear that the benefits of energy recovery are energy, environmental and economic. Howev-
er, it is not easy to quantify them to propose the feasibility of the joint system. This chapter
addresses these aspects, making a proposal on how to estimate this feasibility systematically
and highlighting it in a specific case study, showing the capabilities of the method to make
quantifications and comparisons of various configurations and methodologies of economic
data.
121
Bombas funcionando como Turbinas (PAT)
1. Introducción y objetivos
A lo largo del presente libro se ha hecho énfasis en los aspectos técnicos y energéticos
de las bombas funcionando como turbinas, sin embargo, es fundamental en la decisión
final de la instalación el incluir también aspectos económicos y principalmente en la
relación entre costes y beneficios que se producen a lo largo de la vida útil de este tipo
de instalaciones. El cálculo de diferentes índices económicos debe permitir a los gesto-
res decidir sobre las diferentes soluciones y su viabilidad, proporcionando índices eva-
luadores de tipo económico que engloben tanto los aspectos técnicos y energéticos
(e.g., energía recuperada, potencia a instalar, número de máquinas), como los aspectos
económicos involucrados (e.g., coste de la instalación, beneficios).
En el presente capítulo se describen los parámetros económicos básicos a tener en
cuenta para el análisis de este tipo de instalaciones, así como los diferentes métodos de
evaluación económica. Finalmente se presenta un ejemplo práctico.
2. Consideraciones iniciales
La instalación de los sistemas de recuperación genera unos costes que pueden ser clasi-
ficados en dos clases: costes de inversión inicial (IC por su acrónimo en inglés Invest-
ment Costs) y costes anuales (AC por su acrónimo en inglés Annual Costs) (i.e., explo-
tación, mantenimiento y repuestos). De igual forma dentro de los costes anuales
pueden aparecer los costes de sustitución de las máquinas (RC por su acrónimo Repla-
ce Costs). Los costes de sustitución surgen cuando se pretende explotar la red de distri-
bución de riego durante un tiempo superior al tiempo de vida de las bombas (Portela,
1988; 2000). Asimismo, los ingresos que proporciona un sistema de recuperación se
clasifican en dos clases: ingresos anuales (AI por su acrónimo en inglés Annual Inco-
me) e ingresos por el valor residual de la inversión (RI por su acrónimo en inglés de
Residual Income). A continuación, se describen con más detalle cada uno de las clases
de costes e ingresos comentados.
2.1. Costes de inversión inicial (IC)
Estos costes contemplan la inversión inicial para la ejecución, instalación y funciona-
miento de los sistemas de recuperación, siendo igual a:
. Ecuación 1
_ _
123
Bombas funcionando como Turbinas (PAT)
siendo:
í : Ingresos por venta de energía o por autoconsumo en €.
: Beneficio en € que se produce por la disminución de fugas en la red de
distribución.
: Beneficio que se produce por la disminución en emisiones de CO2 en
comparación con otros sistemas de generación de energía.
Estos beneficios de “sostenibilidad” pueden ser importantes, sobre todo el relacio-
nado con la reducción de fugas. No obstante, el análisis de dichos valores no es abor-
dado en este documento, aunque por su importancia pueden ser considerados por el
gestor responsable a la hora de establecer los beneficios de la instalación de sistemas de
generación renovables. La reducción de fugas está fuertemente ligada a la mejora de la
eficiencia energética de los sistemas y, por tanto, su reducción puede conducir a reducir
el consumo “kW/m3” de nuestra instalación, ya no por la generación sino por la reduc-
ción de las fugas.
Finalmente, los ingresos por valor residual (RI) tienen en cuenta los posibles ingre-
sos por venta de los diferentes elementos una vez ha terminado su función en la insta-
lación.
124
Viabilidad económica en redes de agua con PAT
siendo los ingresos anuales netos en € definidos por la Ecuación 3 y los costes
anuales en € definidos por la Ecuación 2. El no considera los flujos de caja y, por
lo tanto, no es capaz de estimar la eficiencia de la inversión a lo largo de la vida útil.
Además, no tiene en cuenta directamente el coste de oportunidad del capital. Este mé-
todo es un indicador rápido y simple de decisión. En el capítulo anterior se utiliza este
indicador para la selección de soluciones factibles.
Retorno de la inversión (ROI)
El método de retorno de la inversión se calcula como el porcentaje del promedio de
los beneficios totales obtenidos menos los costes iniciales de inversión, respecto a los
costes iniciales de inversión.
Ecuación 6
donde: es el beneficio total obtenido en € e es la inversión inicial definida en la
Ecuación 1. Este índice permite una estimación rápida de las ganancias netas del pro-
yecto y proporciona un indicador válido de comparación entre diferentes proyectos.
Para el cálculo de se deben considerar los rendimientos de toda la vida útil del
proyecto (a diferencia del método PSR que solo considera el período que lleva recupe-
rar la inversión original). Este método ignora el valor temporal del dinero al igual que
el .
3.2. Métodos dinámicos
Estos métodos son más complejos que los métodos estáticos, ya que tienen en cuenta el
valor del dinero a lo largo del tiempo, calculando los flujos de caja a partir de los bene-
ficios y costes a lo largo de la vida útil de la instalación.
Valor actualizado neto (VAN)
El VAN representa la suma acumulada de todos los ingresos menos todos los costes,
esperados durante el tiempo de vida del proyecto, con una tasa de descuento. La tasa de
descuento o tasa de actualización, k, permite asignar un valor a los flujos de caja que
ocurren en diferentes momentos. Si n representa el periodo de n años (vida del proyec-
to), el valor de una unidad monetaria en el presente cambiará en el año n por (1+k)n
unidades monetarias. El valor de una unidad monetaria en el año n cambiará en el pre-
sente por 1/(1+k)n unidades monetarias (Portela, 2000). El VAN viene dado por la ex-
presión:
∑ Ecuación 7
125
Bombas funcionando como Turbinas (PAT)
El RB/C es una comparación entre los flujos de ingresos y los flujos de gastos. Los
proyectos con un RB/C menor que 1 se descartan.
El análisis mediante índices financieros es necesario llevarlo a cabo para que, una
vez se ha aplicado la metodología de optimización descrita en el capítulo anterior del
126
Viabilidad económica en redes de agua con PAT
4. Caso de estudio
Lo descrito hasta ahora, se aplica en un caso de estudio concreto. En este caso de estu-
dio se adapta el análisis desarrollado por Montero, (2018). El caso de estudio está cen-
trado en la red de riego ubicada en un municipio de la provincia de Valencia (España),
particularmente en el "sistema de agua de Canyoles" (Figura 1a). La red presurizada
suministra a 383 ha, donde los cultivos principales son cítricos y el agua es almacenada
en una balsa, abastecida desde un sondeo. La longitud total de las conducciones es de
44 km, donde los diámetros están entre 32 y 500 mm con diferentes materiales. Final-
mente, la red cuenta con 404 puntos de riego, abasteciendo a los diferentes cultivos
(Figura 1b).
127
Bombas funcionando como Turbinas (PAT)
128
Viabilidad económica en redes de agua con PAT
Inversión inicial
Con el objetivo de simplificar el estudio económico se diseña una instalación tipo
para todos los sistemas de recuperación. Los elementos que conforman este sistema
tipo son:
Tuberías 300 mm con una longitud total 5 m.
Tuberías 140 mm con una longitud total 9 metros.
1 válvula de mariposa motorizada 200 mm.
7 válvula de mariposa motorizada 150 mm.
1 válvula reguladora motorizada ( variable en función de Q y H).
6 carretes de desmontaje 140 mm.
3 TE 300/300/300.
5 TE 300/140.
5 TE 300/300/140.
3 PAT (no se disponen de PAT suplentes).
La Figura 2 muestra el esquema de los elementos necesarios.
129
Bombas funcionando como Turbinas (PAT)
Cada grupo de PAT irá instalado en una caseta subterránea o arqueta. Del mismo
modo que en la instalación, independientemente del tamaño de las PAT se ha diseñado
una caseta tipo. La arqueta consistirá en un paralelepípedo enterrado de 9x6 metros de
base y de una altura de 2.5 metros, muros y fondo de espesor de 30 cm, con una losa de
15 cm que servirá de tapa movible. En la Figura 9.3 se muestra un esquema de la ar-
queta de alojamiento de las PAT
El coste de la obra civil necesaria para la ejecución de la arqueta se ha estimado en
7936.40 €. Se estima que el coste de los elementos necesarios para la utilización de la
energía es de 1500 € y 500 € el importe total de impuestos. Para los gastos en ingenie-
ría se estima del 3% del total, descontando impuestos. Por lo tanto, la inversión inicial
y sus componentes quedará según la Tabla 2, sin considerar los costes de legalización
que dependen de la normativa.
Tabla 2. Desglose para la inversión inicial
Concepto Cuantía (€)
_ _ 12777.94
Según Ecuación 10
0.10·
7936.40
. 1500
500
0.03· (
Ingresos anuales
Los ingresos anuales vendrán dados por la venta de energía o el propio autoconsu-
mo realizado por la entidad. Se supone que la energía generada se mantiene constante a
lo largo de todo el período de estudio, calculándose los ingresos como:
í Ecuación 11
siendo energía generada en un año en kWh/año y precio de venta de energía en
€/(kWh) En el caso de autoconsumo, el precio de venta será el precio que se paga por
la energía que se deja de consumir. Los precios de energía considerados para cada es-
cenario son para el escenario neutro de P 0.0842 €/kWh, para el escenario optimista
de P 0.10104 €/kWh) y para el escenario pesimista de P 0.07016 €/kWh.
No se consideran otro tipo de ingresos derivados de subvenciones o ayudas a la ins-
talación de sistemas de generación renovables.
Costes anuales
Los costes de mantenimiento y operación se obtendrán como
ó
130
Viabilidad económica en redes de agua con PAT
131
Bombas funcionando como Turbinas (PAT)
132
Viabilidad económica en redes de agua con PAT
Según la Ecuación 7. el valor del VAN será igual a 1900.79 €, que coincide con el
valor final del flujo de caja según Tabla.2. Los años de retorno de la inversión (paso de
flujo de caja negativos a positivos) será igual a 27.474 años. El valor del TIR será de
0.052 muy cercano a la tasa de descuento. El valor para la ratio beneficio-costes queda-
rá como:
37542.03 670.81
⁄ 1.052
6465.08 29846.98
TIR
0.08
0.07
0.06
0.05
0.04
TIR
0.03
0.02
0.01 Escenario Neutro
Escenario Optimista
0.00
Escenario Pesimista
‐0.01
20 25 30
Años vida útil
133
Bombas funcionando como Turbinas (PAT)
134
Viabilidad económica en redes de agua con PAT
135
Bombas funcionando como Turbinas (PAT)
Figura 7. VAN para el SR2 de la solución óptima con dos sistemas de recuperación
TIR
0.02
0.01
0.00
‐0.01
‐0.02
TIR
‐0.03
‐0.04
‐0.05
Escenario Neutro
‐0.06
Escenario Optimista
‐0.07 Escenario Pesimista
‐0.08
20 25 30
Años vida útil
Figura 8. TIR para el SR2 de la solución óptima con dos sistemas de recuperación
En el caso del TIR, los valores obtenidos son negativos, excepto en el escenario op-
timista y 30 años, donde se obtiene un valor de 0.01 (menor que 0.025). Indicando por
tanto que la inversión en ningún caso es rentable (Figura 8.).
136
Viabilidad económica en redes de agua con PAT
Tasa descuento=0.025 Tasa descuento=0.05
25000 0
20000 Escenario Neutro Escenario Neutro
Escenario Optimista ‐5000 Escenario Optimista
15000
Escenario Pesimista Escenario Pesimista
10000
‐10000
5000
VAN
VAN
0 ‐15000
‐5000
‐10000 ‐20000
‐15000
‐25000
‐20000
‐25000 ‐30000
20 25 30 20 25 30
Años vida útil Años vida útil
Tasa descuento=0.075 Tasa descuento=0.1
0 0
Escenario Neutro Escenario Neutro
‐5000 ‐5000
Escenario Optimista Escenario Optimista
‐10000 Escenario Pesimista ‐10000 Escenario Pesimista
‐15000 ‐15000
VAN
VAN
‐20000 ‐20000
‐25000 ‐25000
‐30000 ‐30000
‐35000 ‐35000
‐40000 ‐40000
20 25 30 20 25 30
Años vida útil Años vida útil
Figura 9. VAN para el global de la solución óptima con dos sistemas de recuperación
TIR
0.06
0.05
0.04
0.03
TIR
0.02
0.01
0.00 Escenario Neutro
Escenario Optimista
‐0.01
Escenario Pesimista
‐0.02
20 25 30
Años vida útil
Figura 10. TIR para el global de la solución óptima con dos sistemas de recuperación
En los años de retorno no se obtienen nunca valores por debajo de la vida útil de la
instalación, y para el índice RB/C, los valores obtenidos son siempre menores que 1. Si
calculamos globalmente los resultados para esta solución óptima los valores para el
VAN y el TIR, quedarán según las Figuras 9 y 10 respectivamente.
Para el escenario neutro, el VAN solo sale positivo y con un valor de 6110.05 € para
30 años y una tasa de descuento de 0.025. Este valor es inferior al obtenido para el caso
de solución óptima con un solo sistema de recuperación. En el escenario optimista solo
se obtienen valores positivos para tasas de descuento de 0.025, siendo siempre estos
valores inferiores a los obtenidos para la solución óptima con un solo sistema de recu-
peración.
137
Bombas funcionando como Turbinas (PAT)
Figura 11. TIR para la solución óptima de 3 sistemas de recuperación (a) SR2 (b) SR3
En la Figura 12 se muestran los valores del VAN globales para la solución óptima de
3 sistemas de recuperación. Solo se obtiene un valor positivo (8047.85 €) para un esce-
nario optimista, 30 años y una tasa de descuento de 0.025. Este valor es inferior a los
obtenidos para la solución óptima de uno y dos sistemas de recuperación.
Para los casos óptimos de soluciones para cuatro y cinco sistemas de recuperación,
los resultados obtenidos desde el punto de vista económico son peores que para tres
sistemas de recuperación, por lo que no adjuntan los resultados.
138
Viabilidad económica en redes de agua con PAT
Conclusiones
El presente capítulo pone de manifiesto la necesidad de analizar la viabilidad de los
sistemas de recuperación. Para ello se ha desarrollado una metodología de trabajo para
llevar a cabo dicho análisis, teniendo en cuenta las expresiones necesarias para ello.
Esta metodología ha sido aplicada a un caso de estudio, analizándose la viabilidad
económica de la instalación de sistemas de recuperación mediante el cálculo de índices
financieros. Estos índices permiten estudiar durante la vida útil de la instalación el flujo
de caja y por tanto es posible determinar la rentabilidad de los sistemas. De los casos
estudiados, la solución óptima con un solo sistema de recuperación es la más rentable
desde el punto de vista económico. En el resto de casos, la adicción de nuevos sistemas
de recuperación (con menor energía recuperada), disminuye los valores del VAN y del
TIR, y por lo tanto disminuyen la rentabilidad global de las instalaciones. Un estudio
más detallado de la obra civil e instalaciones complementarias para los diferentes sis-
temas de recuperación podría llevar a resultados más adecuados desde el punto de vista
económico.
139
Aspectos eléctricos de la PAT
Costa Branco, P.J, Fernandes, JFP
Departamento de Engenharia Eletrotécnica e de Computadores. Instituto Superior Técnico.
Resumen
Hoy en día se sabe que se puede lograr una reducción significativa en el coste de capital de
los equipos eléctricos si el motor de inducción de jaula de ardilla utilizado en los sistemas de
bombeo se "transforma" en un generador de inducción. El motor de inducción "listo para
usar" y la electrónica de potencia de motor de estado sólido estándar son relativamente eco-
nómicos en comparación con los generadores síncronos de imanes permanentes personaliza-
dos o los generadores de inducción de rotor bobinado. En este contexto, el capítulo se centra-
rá principalmente en los desarrollos mini-hidroeléctricos fuera de la red para plantas
independientes, similares a los que se encuentran en comunidades aisladas y sitios de fun-
cionamiento.
Saber cómo garantizar la autoexcitación de un generador de inducción (Self Excite In-
duction Generator, SEIG) y cómo esto afecta la eficiencia de la PAT (Pump as Turbine,
PAT), se basa en la suposición de que el ingeniero hidráulico tiene una comprensión del mi-
niproceso completo de desarrollo hidroeléctrico. En consecuencia, este capítulo también cu-
bre los fundamentos de la ingeniería eléctrica y el diseño eléctrico general de los sistemas
minihidroeléctricos aislados. Si bien este no es un tratamiento exhaustivo del tema, es una
indicación del nivel de comprensión requerido para el diseño eléctrico PAT-SEIG
Abstract
Today it is known that a significant reduction in the capital cost of electrical equipment is
achievable if the squirrel-cage induction motor used in pump systems is “transformed” in an
induction generator. The "off-the-shelf" induction motor and standard solid-state motor
power electronics are relatively inexpensive in comparison to custom permanent-magnets
synchronous generators or wound-rotor induction generators. In this context, the chapter will
141
Bombas funcionando como Turbinas (PAT)
mainly focus on off-grid mini-hydro developments for stand-alone plants, similar to those
found in isolated communities and running sites.
Knowing how to guaranty the self-excitation of an induction generator (Self Excite In-
duction Generator, SEIG) and how this affects PAT efficiency (Pump as Turbine, PAT), is
based on the assumption the hydraulic engineer has an understanding of the complete mini-
hydro development process. Consequently, this chapter also covers the electric engineering
fundamentals and general electrical design concerning isolated mini-hydro systems. While
this is not an exhaustive treatment of the subject matter, it is an indication of the level of un-
derstanding required for the electrical PAT-SEIG design.
1. Introducción y objetivos
El hecho de considerar la PAT un elemento de recuperación de energía hidráulica, no
debe impedir que la parte eléctrica esté muy presente para poder alcanzar la eficiencia
máxima. Precisamente, esa simbiosis “hidraulica-electricidad” es lo que conduce a que
el diseño de estas instalaciones sea correcto y por tanto, de utilidad para los usuarios.
No tiene sentido intentar en el máximo punto de eficiencia hidráulica de la máquina, si
después, la eficiencia del generador es baja, bien por una mala selección o sencillmante
porque la potencia del generador hace que su eficiencia no supere el 60%. Por tanto, el
gestor debe buscar la máxima eficiencia global y en este caso, eso conlleva maximizar
la eficiencia hidráulica y eléctrica de la máquina. Entre los objetivos de este capítulo se
encuentra el introducir los conceptos básicos relacionados con la generación de la
energía y la explicación del funcionamiento de un generador. Asi mismo, se
presentarán los esquemas de funcionamiento tanto conectados a red como aislados de
la misma.
capítulo. En este contexto, nos centraremos solo en circuitos eléctricos de corriente alter-
na (AC) lineales, ya que pueden ilustrar y explicar la mayoría de los aspectos clave de la
electricidad en AC necesaria para las unidades de energía hidráulica.
2.1. Tensión alterna y señal de corriente (AC)
En general, la variación del tiempo de la onda de tensión de alterna, , y de la co-
rriente pueden ser definidas, respectivamente, mediante las siguientes expresiones
√2 cos , Ecuación 1
√2 cos , , Ecuación 2
donde [V] e [A] son, respectivamente, los valores de voltaje efectivo (rms) y
corriente; √2 y √2 son los valores de amplitud de voltaje y corriente; [rad/s]
es la frecuencia angular eléctrica dada por 2 , siendo [Hz] la frecuencia
eléctrica de las señales de voltaje y corriente; , es la señal de fase; al tiempo
, con , y , siendo los ángulos de fase del voltaje y corriente para 0.
Si cualquiera de los ángulos de fase iniciales es nulo, , 0o , 0, la fun-
ción se convierte en una onda cosenoidal que pasa por su valor máximo en t = 0. Esto
puede considerarse como una forma de onda de "referencia". Otras formas de onda de
voltaje y corriente con diferentes ángulos de fase iniciales se desplazan a lo largo del
eje de tiempo con respecto a la forma de onda de referencia.
La Figura 1 muestra curvas sinusoidales genéricas con diferentes ángulos iniciales
( , 0, en negro); , /2, en rojo; , /2, en azul). La curva negra es
la referencia, con su punto máximo en t = 0, y puede describirse por 0 . La
curva azul está relacionada de antemano con la referencia, es decir, conduce con un ángulo
de fase inicial de /2, descrita como cos /2 . Su máximo tiene lugar en
(cos /2 1 , por tanto ocurre cuando /2 2 . Por lo tanto, el punto
máximo más cercano a t=0 es /2, como se puede ver en la Figura 8.1. Por el
contrario, la curva roja está retrasada en relación con la referenciada, con un ángulo de
fase inicial de /2. En la figura también se puede ver que el período de las curvas es
de 20 ms, que corresponde a una frecuencia eléctrica de 50 Hz.
143
Bombas funcionando como Turbinas (PAT)
Figura 2. (a) Señal de voltaje continua, (b) Resistencia con intensidad de corriente
Si el voltaje de AC es aplicado a una resistencia (R), representada en la Figura 2(b),
la intensidad de corriente ( cos , ) circulará como se muestra en la Figura
2(a), mediante la ecuación:
,
→ cos , . Ecuación 3
0,
144
Aspectos eléctricos de la PAT
cia eléctrica, es decir, ninguna energía se disipa en calor. Por lo tanto, la bobina "ideal"
solo almacena energía magnética.
cos , , Ecuación 4
2
0,
Por tanto, la expresión que determina la corriente que circula por el inductor será
igual a
cos , Ecuación 5
145
Bombas funcionando como Turbinas (PAT)
→ sin , sin ,
cos , , Ecuación 6
2
0,
Por tanto, la expresión que determina la corriente que circula por el condensador se-
rá igual a:
cos , Ecuación 7
La Figura 4(a) muestra que mientras la corriente retrasa el voltaje en
90º ( , ), en un condensador la adelanta 90º. De forma intuitiva, la corriente
eléctrica “aparecerá” primero para establecer cargas eléctricas en las placas del con-
densador, para que posteriormente, aparezca voltaje entre las placas.
2.2. Potencia en régimen sinusoidal y permanente
Potencia instantánea
Si se define un circuito elemental (resistencia, inductor o condensador) donde circu-
la una corriente eléctrica, √2 , , y una señal de voltaje
√2 , . La potencia instantánea es dada por
. Ecuación 8
Operando considerando las ecuaciones e anteriormente presentadas:
2 , , , , . Ecuación 9
La ecuación anterior muestra que la potencia instantánea está formada por dos par-
tes. La primera de ellas corresponde con la componente oscilante, que tiene una fre-
146
Aspectos eléctricos de la PAT
(o,u o,i )
147
Bombas funcionando como Turbinas (PAT)
Figura 6. Ejemplos de máquinas inductoras instaladas en bombas (a) y (b) Grundfos (2004)
La Figura 6(a) muestra una máquina fabricada en China. Primero se debe buscar la
potencia nominal, a veces denominada potencia útil del motor si se conecta directamente
a la fuente de alimentación. En este caso, la potencia nominal del motor indicada en la
placa de identificación es de 5,5 kW. Este es el valor de potencia neta que el motor tendrá
disponible en su eje cuando opere en sus condiciones nominales (i.e., corriente nominal
del estator, velocidad nominal, voltaje nominal fase-fase, frecuencia nominal)
En la misma placa de identificación, presta atención en la frecuencia nominal del
motor. Indica dos valores de frecuencia eléctrica (50 Hz y 60 Hz), lo que significa que
la bomba puede funcionar en países donde hay una red eléctrica de 50 Hz (e.g., Euro-
pa) o en 60 Hz (e.g., Estados Unidos, América Latina). Los datos de la placa de identi-
148
Aspectos eléctricos de la PAT
ficación se dividen en dos conjuntos de valores nominales que se muestran en dos co-
lumnas, una relativa a 50Hz y otra para 60Hz.
La Tabla 1 muestra los valores nominales del motor para 50 y 60 Hz =s valores de
50Hz y su interpretación se enumeran a continuación. Para 60 Hz es similar y no se
mostrará. Por lo tanto, para 50Hz se tiene:
Rango de valores de tensión línea a línea (VOLT. 380-415) cuando el motor
está conectado en triángulo () a 50 Hz;
Rango de valores de tensión línea a línea (VOLT. 660-690) cuando el motor
está conectado en estrella () a 50 Hz;
La línea RMS y la intensidad nominal (AMP. 11) cuando el motor está conecta-
do en triángulo () a 50 Hz; o (AMP. 6.4) siendo la corriente de fase (en este ca-
so igual a la corriente de línea) cuando se conecta en estrella (Y);
El rango de la velocidad de carga nominal para 50 Hz (R.P.M. 2900-2920) y pa-
ra 60 Hz (R.P.M. 3470-3525);
El rango de factor de potencia nominal a 50 Hz (COS 0.89-0.86) y para 60 Hz
(COS 0.92-0.88), and;
La eficiencia de la máquina (EFF), que es válida para 50Hz o 60Hz, dado por
dos valores: uno para carga completa (100% de carga nominal) igual a 85.7%, y
un segundo valor para 75% de carga completa igual a 86%.
Tabla 1. Ejemplo de interpretación para placas identificativas
50 Hz 60 Hz
VOLT. 380-415/660-690 VOLT. 380-480/660-690
AMP. 11 / 6.4 AMP. 10.5-8.6/6.1-5.0
R.P.M. 2900-2920 R.P.M. 3470-3525
COS 0.89-0.86 COS 0.92-0.88
EFF. (100%FL) 85.7% EFF. (75%FL) 86%
149
Bombas funcionando como Turbinas (PAT)
150
Aspectos eléctricos de la PAT
151
Bombas funcionando como Turbinas (PAT)
El grosor del rotor de estampación y los núcleos del estator generalmente varía de
0.3 a 0.5 mm, estando las ranuras perforadas en ambos. Las perforaciones están reali-
zadas en el rotor, en el lado exterior de los estampados y en el estator, en el lado in-
terno de los estampados. Mientras que, en el rotor, el aluminio se filtra a través de las
ranuras para formar la jaula del rotor, las ranuras del estator llevarán devanados trifási-
cos. Generalmente, estos devanados se fabrican con alambres de cobre esmaltados que
se muestran en la Figura 9(c).
Finalmente, la Figura 9(d) muestra los bobinados de cobre que se colocan en cada
ranura del núcleo del estator. Al ser una máquina eléctrica trifásica, tendrá tres juegos
de bobinado, cada uno para una fase, todos distribuidos alrededor de la periferia exte-
rior del núcleo del estator, como se muestra en la Figura 9(e).
3.2. Principio de operación de un motor de inducción
Para comprender el funcionamiento de un generador de inducción autónomo, es decir,
no conectado a la red, es muy útil tomar el funcionamiento normal del motor como
punto de partida. El principio de funcionamiento de un motor de inducción puede des-
cribirse mediante un pequeño conjunto de pasos secuenciales. Estos se describen a
continuación:
1. La Figura 10 muestra que cuando los bobinados del estator trifásico están conec-
tados a una fuente de alimentación trifásica y equilibrada (tensiones U, V, W),
aparece un campo magnético giratorio en el espacio de aire entre el estator y el
rotor. En la Figura 10, el campo magnético se especifica con un color verde, que
gira en la dirección del reloj con velocidad angular ω [rot/s]. De manera simi-
lar, se indica que el rotor giraba con velocidad angular ω [rot/s]. La unidad
"rot / s" significa una rotación mecánica por segundo. Cuando una máquina de
inducción está conectada a un suministro de CA, la corriente de magnetización
fluye desde el suministro y crea un campo magnético giratorio en la máquina.
en la Figura 11. El campo giratorio corta las barras del rotor en cortocircuito, in-
duciendo corrientes en ellas. Ahora, estas corrientes interactúan con el campo
del estator giratorio y se produce un par. Este par arrastra al rotor con el campo,
pero a una velocidad ligeramente inferior.
Sin un par externo aplicado, el deslizamiento del motor de inducción será muy
pequeño, menor que s = 0.01 (o 1%). Para una máquina de 1 kW, el desliza-
miento a plena carga será entre 0,02 y 0,05 (2%-5%). Las máquinas más grandes
tienen deslizamientos más pequeños.
4. La interacción entre las corrientes eléctricas en las barras de aluminio y el cam-
po magnético deslizante produce una serie de fuerzas tangenciales similares a lo
largo del perímetro del rotor. La Figura 11 muestra estas fuerzas usando flechas
rojas distribuidas por las barras del rotor. Debido a la simetría, se produce un par
electromagnético resultante en el rotor en la misma dirección de rotación del
campo magnético. Cuando se aplica una carga al motor, la diferencia de veloci-
dad aumenta a medida que se produce un mayor par.
153
Bombas funcionando como Turbinas (PAT)
154
Aspectos eléctricos de la PAT
155
Bombas funcionando como Turbinas (PAT)
Figura 13. La P activa y la potencia reactiva Q fluyen en los terminales del estator de la máqui-
na de inducción cuando funcionan como (a) motor y (b) generador. La potencia activa cambia
su dirección del motor al generador, mientras que la potencia reactiva permanece
La Figura 14 muestra dos diagramas esquemáticos típicos para un generador de in-
ducción autoexcitado trifásico. La Figura 14 (a) muestra un generador de inducción
conectado en estrella con un banco de condensadores ya conectado en estrella, mientras
que en la Figura 14 (b) un generador de inducción conectado en triángulo tiene un ban-
co de condensadores conectado en triángulo. Un grupo de varios condensadores de la
misma clasificación que están conectados en serie o en paralelo entre sí para almacenar
energía eléctrica, forma este banco. El banco resultante se usa para suministrar la po-
tencia reactiva al generador de inducción.
156
Aspectos eléctricos de la PAT
Figura 14. Dos diagramas esquemáticos típicos para generadores de inducción autoexcitados
trifásicos. (a) banco de condensadores conectado en estrella, (b) conectado en triángulo
Para comenzar el proceso de autoexcitación del generador de inducción, el banco de
condensadores debe cargarse eléctricamente. Con una batería, se puede conectar prime-
ro en paralelo con todos los condensadores. Después de que la corriente eléctrica en
cada condensador se vuelve casi nula (esto indica que todos los condensadores están
cargados), ahora pueden desconectarse de la batería y conectarse a los devanados del
estator del generador.
Para los casos en que la máquina de inducción se ha utilizado en situaciones anterio-
res (como motor o incluso como generador de inducción), generalmente el hierro en el
rotor todavía está magnetizado, lo que generalmente se denomina "magnetización rema-
nente". En este caso, no es necesario cargar previamente el banco de condensadores.
Con la turbina impulsando la máquina de inducción magnetizada, aparece una señal
de voltaje inicial en los terminales del estator (generalmente voltaje muy pequeño). Dado
que el banco de condensadores está conectado, un conjunto inicial de corrientes comien-
za a cargarlas en el primer medio ciclo y se descarga en el otro medio ciclo, proporcio-
157
Bombas funcionando como Turbinas (PAT)
158
Aspectos eléctricos de la PAT
En el punto de operación P, en la Figura 8.15 (a), se verifica que el voltaje sin carga
se establece por el valor de [F] en el banco de condensadores conectado. Por lo
tanto, su valor viene dado por la Ecuación 10 a la frecuencia nominal del generador
de inducción (el generador también gira a su velocidad nominal) y para una determina-
da corriente de condensador .
Ecuación 15
0
159
Bombas funcionando como Turbinas (PAT)
160
Aspectos eléctricos de la PAT
Figura 16. (a) conexiones atornilladas para un estator conectado en Y y (b) el diagrama esque-
mático correspondiente. (c) conexiones atornilladas para un estator conectado en Y y (d) el
diagrama esquemático correspondiente. (Grunfos, 2004)
Usando solo la información en la placa de identificación del motor, la cantidad de
capacitancia de excitación que debe conectarse en una estrella se calcula de la siguiente
manera:
De la placa de identificación, la corriente de carga completa (nominal) es 1.6 A y el
factor de potencia es 0.74.
1) La potencia aparente total S a plena carga es:
√3 √3 400 1.6 1108 VA Ecuación 18
2) Por lo tanto, la potencia de entrada transmitida al generador de inducción por la
PAT será igual a:
√3 √3 400 1.6 0.74 820 W Ecuación 19
.
Se observa que la potencia a plena carga (nominal) en la placa de identificación,
550 W, debe interpretarse como la potencia nominal producida por el generador. Por lo
tanto, para una potencia de entrada de 820 W suministrada por el PAT, el generador
produce 550 W, lo que muestra una eficiencia del 67%.
3) La potencia reactiva se puede obtener utilizando la relación entre potencia aparen-
te, activa y reactiva dada por:
Ecuación 20
√ 1108 550 ≅ 962 VAR (Voltio Ampere reactiva) Ecuación 21
Por tanto, la potencia reactiva por fase será:
3 → ≅ 321 (Voltio Ampere Reactiva) Ecuación 22
161
Bombas funcionando como Turbinas (PAT)
Ecuación 24
≅ 60 F Ecuación 25
162
Aspectos eléctricos de la PAT
163
Bombas funcionando como Turbinas (PAT)
Figura 19. (a) Eficiencia del SEIG, (b) Eficiencia de la PAT y (c) Eficiencia del Sistema
PAT-SEIG, en función de la carga y velocidad. Fernandes et al., (2019)
164
Aspectos eléctricos de la PAT
165
Bombas funcionando como Turbinas (PAT)
En este ejemplo, se supone que el generador de inducción que funciona con la PAT
tiene dos pares de polos p = 2. La velocidad sincrónica será entonces de 1500 rpm.
Si la PAT gira por debajo de la velocidad síncrona ( <1500 rpm), la máquina
de inducción funcionará como un motor. La unidad PAT-IG comienza a extraer energía
de la red;
Si la PAT gira por encima de la velocidad síncrona ( > 1500 rpm), entra en su
modo de generación. La unidad PAT-IG comienza a suministrar energía a la red.
5.2. Ejemplo de generación energía en la red cuando una unidad PAT-IG está
conectada a una red de 50 Hz/230 V
La Tabla 4 enumera para 60 Hz el par más común de topologías de polos existentes en
el mercado de máquinas de inducción y sus respectivos valores de velocidad síncrona.
Cuando la PAT-IG está conectada a una red de 60 Hz y tiene, por ejemplo, también dos
pares de polos, la nueva velocidad síncrona del generador será de 1800 rpm.
Si la PAT gira por debajo de la velocidad síncrona ( <1800 rpm), la máquina
de inducción funcionará como motor. La unidad PAT-IG comienza a extraer energía de
la red;
Si la PAT gira por encima de la velocidad síncrona ( > 1800 rpm), entra en su
modo de generación. La unidad PAT-IG comienza a suministrar energía a la red.
Tabla 4. Velocidad síncrona en función del número de pares de polos para 60Hz
Pares de polos 60 Hz – Velocidad síncrona rpm 60 Hz – Velocidad síncrona rot/s
1 3600 60.0
2 1800 30.0
3 1200 20.0
4 900 15.0
166
Aspectos eléctricos de la PAT
167
Bombas funcionando como Turbinas (PAT)
Figura 20. (a) Generador de inducción trifásico de conexión en triángulo. (b) generador de
inducción trifásico está conectado a un transformador trifásico antes del interruptor de circuito
168
Aspectos eléctricos de la PAT
169
Bombas funcionando como Turbinas (PAT)
170
Aspectos eléctricos de la PAT
Conclusiones
El presente capítulo ha desarrollado una descripción de los principios básicos eléctricos
que cualquier ingeniero hidráulico que considere la instalación de microsistemas de
generación hidráulica debe conocer. Del mismo modo, se han desarrollado los esque-
mas necesarios para la conexión de los elementos, desarrollando ejemplos de aplica-
ción.
171
Introducción al funcionamiento de sistemas
de agua a presión en régimen transitorio
operando con PAT
Pérez-Sánchez, Ma,, Sánchez-Romero, FJb, Ramos, HMc, .López Jiménez, PAa
a
Departamento de Ingeniería Hidráulica y Medio Ambiente. Universitat Politècnica de València.
b
Departamento de Ingeniería Rural y Agroalimentaria. Universitat Politècnica de València.
c
CERIS, Instituto Superior Técnico, Universidade de Lisboa.
Resumen
El conocimiento de las condiciones en régimen transitorio en redes de agua presurizadas que
están equipadas con bombas que funcionan como turbinas (PAT) es de sumo interés, puesto
que prevé un funcionamiento muy poco tratado por la bibliografía. El presente capítulo
caracteriza la variación de la velocidad de rotación en las máquinas, ligado con situaciones
en régimen no permanente ocurridas en redes de agua en el diseño de los sistemas equipados
con PAT, enfatizando los eventos transitorios que pueden ocurrir durante una operación
normal. Se analizan las decisiones a tomar por los operadores del proyecto hacia una
operación estable y eficiente asociada al comportamiento dinámico normal de un cierre de
válvula de control de flujo o por efecto de sobrevelocidad inducida. Los conceptos básicos
del modelado matemático, la caracterización del comportamiento de las válvulas de control,
los efectos de amortiguación en la propagación de las olas y las condiciones fuera de control
de las PAT están actualmente asociados a un diseño inadecuado. La evaluación correcta de
las reglas básicas de funcionamiento depende del sistema y del tipo de componentes, así
como del nivel de seguridad requerido durante cada operación. Se particulariza en este
capítulo en el análisis del funcionamiento no permanente.
Abstract
The knowledge of transient conditions in water pipe systems which are equipped with PAT
is of the utmost importance. This research characterizes the rotational speed variation, water
hammer phenomenon in the design of PAT systems, emphasizing the transient events that
can occur during a normal operation. This is based on project concerns towards a stable and
efficient operation associated to the normal dynamic behaviour of a flow control valve
closure or by induced overspeed effect. Basic concepts of mathematical modelling,
characterization of control valves behaviour, damping effects in the wave propagation and
173
Bombas funcionando como Turbinas (PAT)
runaway conditions of PAT are currently associated to an inadequate design. The correct
evaluation of basic operating rules depends upon the system and components type, as well as
the safety level required during each operation.
Keywords: Variable rotational speed; runaway conditions; unsteady flow conditions; water
hammer.
1. Introducción
La operación de grupos generadores con rotación variable ha ganado importancia en los
últimos años, impulsado además por el desarrollo de plantas eólicas, donde la gran
variabilidad de la velocidad del viento exige el ajuste de la turbina, evitando velocidades
excesivamente altas (Harrison et al., 2000). En las centrales hidroeléctricas, el uso de la
velocidad variable alcanzó una gran importancia con el desarrollo de plantas reversibles,
donde la misma máquina ejerce la función de bomba, durante el período de baja demanda
energética, y turbina, cuando la demanda es alta, y normalmente la tarifa elevada,
viabilizando económicamente este tipo de central (Deane et al, 2010). Pannatier et al.,
(2010) destacan el aumento de la eficiencia operacional, además de una mayor
estabilidad eléctrica e hidráulica del sistema.
A pesar de las ventajas observadas del uso de la velocidad rotación variable, todavía
se aplica poco en sistemas de PAT. En el caso de bombas comerciales operando como
turbina, Joshi et al., (2005) propusieron su uso para atender comunidades aisladas,
destacando la mayor capacidad de producción de energía. Carravetta et al.,(2013)
desarrollaron estrategias de operación en redes de abastecimiento de agua como forma
de mejorar el control de presión. Los resultados mostraron una buena respuesta en
casos donde no hay una gran variación de la demanda.
El funcionamiento de PAT en régimen de velocidad nominal se encuentra muy
referenciado y ha sido descrita y discutida a lo largo de los capítulos anteriores del
libro. La revisión de las tecnologías disponibles ha sido desarrollada por Senior et al.,
(2010); Elbatran et al., (2015). Nourbakhsh y Jahangiri (1992) abordaron el análisis del
rendimiento y el modelado en PAT; Arriaga (2010); Simão y Ramos (2010); Razan et
al., (2012); Fecarotta et al., (2016) entre otros, mientras que el análisis de estas
máquinas en la red de distribución de agua fue estudiado por Ramos et al., (2009);
Caxaria et al., (2011); Carravetta et al., (2014); Butera y Balestra (2015), diseñando
una estrategia innovadora para maximizar la energía recuperada cuando los flujos
varían a lo largo del día (Carravetta et al.,2013, 2014; Fecarotta et al.,2014, 2016).
Estas estrategias se han aplicado para determinar y maximizar la energía recuperada
teóricamente en los sistemas de agua potable y riego (Samora et al.,2016; Pérez-
Sánchez et al.,2016).
174
Introducción al régimen transitorio en PAT
2. Efectos dinámicos
Los transitorios hidráulicos son los regímenes causados por perturbaciones bruscas de
caudal y presión durante un cambio de un estado permanente a otro. Estos deben
considerarse desde las primeras fases de diseño del esquema de funcionamiento de las
máquinas en modo de generación de energía hidroeléctrica para encontrar las
especificaciones técnicas correspondientes a un diseño más económico y seguro. Los
cambios de caudal son irremediables y por tanto, cualquier máquina motora debe
arrancar, sufrir cambios de carga y/o apagarse. Además, en la instalación pueden surgir
eventos impredecibles, como errores humanos, fallos en los equipos o peligros
ambientales, también pueden causar regímenes inestables graves.
En los sistemas presurizados, todos los componentes deben soportar las
sobrepresiones transitorias máximas y las presiones mínimas debidas a los cambios de
caudal. La solución óptima final debe considerar los transitorios hidráulicos con todos
los efectos dinámicos colaterales así como la interacción con los equipos
hidroeléctricos:
a. Las condiciones máximas y mínimas de presión serán el factor principal en la
definición de las condiciones de seguridad operacional. Estas condiciones
deberan: evitar la rotura de la tubería, la separación de la columna de agua o la
entrada de aire en el sistema y la superposición de las paredes del canal o de la
zona de acceso.
175
Bombas funcionando como Turbinas (PAT)
176
Introducción al régimen transitorio en PAT
177
Bombas funcionando como Turbinas (PAT)
178
Introducción al régimen transitorio en PAT
Ecuación 1
/ ѱ
179
Bombas funcionando como Turbinas (PAT)
3. Regímenes transitorios
3.1. Modelado hidráulico básico
Se puede analizar el flujo en régimen no permanente en sistemas de tuberías
presurizadas con una longitud mayor que el diámetro, considerando modelos
unidimensionales (1D), a través de las ecuaciones de conservación de masa y momento
que se derivan del teorema de transporte de Reynolds (White, 2008). Estos principios
están definidos por las Ecuaciones 5 y6 (Wylie et al., 1993):
0 Ecuación 5
0 Ecuación 6
180
Introducción al régimen transitorio en PAT
; ; | | Ecuación 8
181
Bombas funcionando como Turbinas (PAT)
182
Introducción al régimen transitorio en PAT
Figura 2. Tipo de cierre en función del tipo de válvula del coeficiente de flujo
(Adapatada de Iglesias-Rey et al., 2004)
Los resultados descritos en las Figuras 3 y 4 muestran que para maniobras rápidas
ocurre un pico adicional de presión. Asimismo, cuando el cierre o apertura es parcial, ,
la variación de la presión no es significativa. Sin embargo, para un arranque o parada
en el caso de estar funcionando con PAT, es conveniente un funcionamiento a través de
maniobras lentas. La Figura 3 muestra que el cierre lento tanto parcial como total, no
supone ningún riesgo de sobrepresión, mientras que las maniobras rápidas pueden
conducir a sobrepresiones del 300% respecto a la altura en regimen permanente. En el
caso de maniobras de apertura, Figura 4, los resultados son similares, siendo más
acusada la depresión generada y por tanto, la aparición de riesgo de cavitación.
183
Bombas funcionando como Turbinas (PAT)
184
Introducción al régimen transitorio en PAT
Figura 5. Variación H/H0 (ascendente y descendente) y Q/Q0 en una válvula de bola para flujo
turbulento (izquierda; Re = 100000) y flujo laminar (derecha; Re = 1000). Comparación entre
cierre efectivo y cierre total (adaptada Pérez-Sánchez (2017))
3.4. Condiciones de embalamiento
La velocidad de embalamiento es aquella velocidad de rotación que alcanza una
máquina motora cuando el par resistente se anula, generalmente por desconexión del
generador.En las turbinas de reacción con baja velocidad específica, el caudal
disminuye con la sobrevelocidad transitoria. A la inversa, para turbinas con alta
velocidad específica, la descarga transitoria tiende a aumentar (Ramos, 1995, 2000).
El flujo que atraviesa el rodete se caracteriza por tres tipos de velocidades:
velocidad absoluta del agua (V) con la dirección impuesta por el conjunto de aletas
guía, velocidad relativa (W) a través del corredor y velocidad tangencial C, del
impulsor (Figura 6).
185
Bombas funcionando como Turbinas (PAT)
186
Introducción al régimen transitorio en PAT
Figura 7. Comportamiento experimental de una PAT (ns 51 rpm (m,kW)) cuando se produce la
decarga de la máquina y para diferentes caudales iniciales
Se observa que las presiones aguas abajo no se alteran significativamente. Sin
embargo, hay una pequeña reducción del flujo del sistema, lo que explica el aumento de
la presión de agua. De acuerdo con Chapallaz (1992) y Ramos (1995 y 2000), esta es una
característica de máquinas radiales, que producen mayores pérdidas con el aumento de la
velocidad de rotación debido a las fuerzas centrífugas que actúan en el fluido.. La Figura
8 muestra la reducción del caudal en el sistema en función del aumento de la velocidad
de giro. De esta forma, además de los problemas mecánicos observados en la PAT
causados por el aumento de la vibración, existe también el riesgo de desabastecimiento
de regiones más alejadas en una red real, ya que hay un gran aumento en la pérdida de
carga localizada, reduciendo el caudal disponible para los consumidores.
187
Bombas funcionando como Turbinas (PAT)
Figura 8. Reducción del caudal en función de la velocidad de rotación para una PAT con
ns = 51 (rpm) (m, kW); Qo=3,36 l/s y H=4 m. Caso experimental desarrollado en IST
Según Ramos (1995), la relación entre la descarga de flujo en condiciones fuera de
control ( ) y la descarga para las condiciones iniciales de apertura total ( ). En
este caso el ratio tiende a aumentar linealmente con el aumento de la velocidad
específica, en particular para máquinas de mayor inercia (Figura 9).
188
Introducción al régimen transitorio en PAT
189
Bombas funcionando como Turbinas (PAT)
Conclusiones
En este capítulo final, se ha abordado la introducción al regimen transitorio cuando se
opera con máquinas de microgeneración, que el usuario debería considerar a la hora de
llevar a cabo el análisis. El tipo de análisis estará influenciado por la etapa de diseño y
la complejidad de cada sistema. Por lo tanto, según las características de cada sistema
hidráulico, para las maniobras más predecibles, los diseñadores podrán definir las
reglas de explotación de acuerdo con los niveles de seguridad esperados. De hecho, es
necesario especificar reglas operativas convenientes para controlar las presiones
transitorias (tanto máxima como mínima). Estas especificaciones dependerán
principalmente de los siguientes factores:
Las características del sistema de tuberías a proteger; de hecho, estas
características basadas en la pérdida de carga y la inercia de la columna de agua
pueden modificar de manera adversa el comportamiento del sistema y el mismo
tiempo de cierre de la válvula puede inducir si un cambio de flujo lento o rápido;
Las características intrínsecas de la válvula: una válvula de mariposa (por
ejemplo, para alturas medias) y una válvula esférica (por ejemplo, para cabezas
altas) tienen diferentes efectos en la respuesta de flujo dinámico para la misma
ley de cierre;
190
Introducción al régimen transitorio en PAT
A consecuencia que las PAT no tienen una aleta guía, el control de flujo se
realiza a través de válvulas donde las leyes de cierre y apertura son cruciales en
las condiciones del sistema de seguridad, como el tipo de actuador de válvula;
En función de las características de la bomba como máquina de turbina (radial o
axial), diferentes comportamientos dinámicos se asociarán a:
o la pequeña inercia de las masas en rotación induce un rápido efecto de
sobrevelocidad en condiciones de fuga impuesta por un rechazo de carga
total.
o los efectos de sobrevelocidad provocan variaciones de flujo (es decir,
reducción de flujo en máquinas de bajo ns y aumento en las máquinas de alto
ns) y variaciones de presión que pueden propagar los purgas hacia arriba de
una máquina radial y los purgas hacia abajo de la misma.
Sin duda, el análisis complejo del comportamiento en régimen no permanente de las
máquinas, particularmente en su asociación en serie y en paralelo, incluyendo
condiciones de velocidad variable, entra de lleno en los desarrollos futuros de
la presente tecnología.
191
Bibliografía
R Abbasi, T.; Abbasi, S.A., 2011. Small hydro and the environmental implications of
its extensive utilization. Renew. Sustain. Energy Rev. 2011, 15, 2134-2143.
ABB, 2014. IEC 60034-30-1 standard on efficiency classes for low voltage AC motors.
Technical note. Available from the internet: <https://library.e.abb.com/public/
1018a82e36b29462c1257d41002b3470/TM025%20EN%2008-2014%20IEC60034-
30-1_lowres.pdf > 2 p
Abreu, J., Guarga, R. & Izquierdo, J., 1995. Transitorios y oscilaciones en sistemas
hidráulicos a presión J. Abreu, R. Guarga, & J. Izquierdo, eds., Valencia: U.D. Me-
cánica de Fluidos. Universitat Politècnica de València.
Adhau, S.; Moharil, R.; Adhau, P., 2012. Mini-hydro power generation on existing
irrigation projects: Case study of Indian sites. Renewable and Sustainable Energy
Reviews, vol. 16, pp. 4785-4795.
Agarwal, T. Review of pump as turbine (PAT) for micro-hydropower. Int. J. Emerg.
Technol. Adv. Eng. 2012, 2, 163-169.
Alatorre-Frenk, A., 1994. Cost minimization in micro-hydro systems using pumps as
turbines. phD Thesis. University of Warwick. Reino Unido.
Almandoz, J.; Cabrera, E.; Arregui, F.; Cabrera, E.; Cobacho, R., 2005 Leakage As-
sessment through Water Distribution Network Simulation. J Water Res Pl-ASCE
131(6):458-466.
Almeida, A.B.; Koelle, E., 1992. Fluid transients in pipe networks. Computational
Mechanics Publications. Elsevier Applied Science, Amsterdam, Netherlands.
Alonso-Tristán, C.; González-Peña, D.; Díez-Mediavilla, M., Rodríguez-Amigo, M.,
García-Calderón, T., 2011. Small hydropower plants in Spain: A case study.
Renewable and Sustainable Energy Reviews, vol. 15, nº 6, pp. 2729-2735.
193
Bombas funcionando como Turbinas (PAT)
Andrade, M. A.; Rojano, F.; Romero-Gómez, P.; & Choi, C. Y., 2010. Integrated Wa-
ter Quality Modeling of Water Distribution Systems. In Proc. ASCE 12th Annual
Conference on Water Distribution Systems Analysis, Tucson, Ariz.
Araujo L.; Ramos, H.; Coelho, S., 2006. Pressure control for leakage minimisation in wa-
ter distribution systems management. Water Resources Management 20(1):133-149.
Arias, O.; Fidias, G., 1999. El Proyecto de Investigación: Guía para su elaboración.
Episteme. Caracas.
Arriaga, M., 2010. Pump as turbine – A pico-hydro alternative in Lao People’s
Democratic Republic. Renewable Energy, vol. 35, pp. 1109-1115.
Baghaee, H. R., Mirsalim, M., Gharehpetian, G. B., Kaviani, A. K., 2012.
Security/cost-based optimal allocation of multi-type FACTS devices using multi-
objective particle swarm optimization. SIMULATION, 88(8), 999-1010.
https://doi.org/10.1177/0037549712438715
Baghaee, H. R., Mirsalim, M., Gharehpetian, G. B., Talebi, H. A., 2016.
Reliability/cost-based multi-objective Pareto optimal design of stand-alone
wind/PV/FC generation microgrid system. Energy, 115, 1022–1041.
https://doi.org/10.1016/j.energy.2016.09.007
Baghaee, H. R., Mirsalim, M., Gharehpetian, G. B., 2017. Multi-objective optimal
power management and sizing of a reliable wind/PV microgrid with hydrogen
energy storage using MOPSO. Journal of Intelligent & Fuzzy Systems, 32(3),
1753-1773. https://doi.org/10.3233/JIFS-152372
Barbarelli S., Amelio M., Florio G., 2016. Predictive model estimating the perfor-
mances of centrifugal pumps used as turbines. Energy; 107: 103-121.
https://doi.org/10.1016/j.energy.2016.03.122
Berndtsson, J. C.; Jinno, K., 2008. Sustainability of urban water system: examples from
Fukuoka, Japan. Water Policy 10, 501-513.
Binama M, Su W-T, Li X-B, Li F-C, Wei X-Z, An S. Investigation on pump as turbine
(PAT) technical aspects for micro hydropower schemes: A state-of-the-art review.
Renew Sustain Energy Rev. 2017; 79: 148-179. doi:10.1016/J.RSER.2017.04.071
Bogdanović-Jovanović, J.B.; Milenkovi, D.R.; Svrkota, D.M.; Bogdanovi, B.; Spasi,
I.T. (2014). Pumps used as Turbines. Power ha, Energy Efficiency, CFD Analysis.
Thermal Sc.: Vol 18, No. 3, pp. 1029-1040.
Boyd, G. R.; Wang, H.; Britton, M. D.; Howie, D. C.; Wood, D. J.; Funk, J, E.; Fried-
man, M. J. (2004a). Intrusion within a simulated water distribution system due to
hydraulic transients I: Description of test rig and chemical tracer method. Journal
of Envir. Engin. 130(7), ASCE, pp. 774-777.
194
Bibliografía
Boyd, G. R.; Wang, H.; Britton, M. D.; Howie, D. C.; Wood, D. J.; Funk, J, E.; Fried-
man, M. J. (2004b). Intrusion within a simulated water distribution system due to
hydraulic transients II: Volumetric method and comparison results. Journal of En-
vironmental Engineering 130(7), ASCE, pp. 778-783.
Butera, I.; Balestra, R. (2015). Estimation of the hydropower potential of irrigation
networks. Renewable and Sustainable Energy Reviews, vol. 48, pp. 140-151.
Cabrera, E.; Cobacho, R.; Soriano, J. (2014). Towards an Energy Labelling of Pressurized
Water Networks. Procedia Eng. 2014; 70: 209-217. doi:10.1016/j.proeng.2014.02.024.
Campisano, A.; Creaco, E.; Modica, C. (2009). RTC of valves for leakage reduction in
water supply networks. J Water Res Pl-ASCE 136(1):138-141.
Capelo, Bernardo; Pérez-Sánchez, Modesto; Fernandes, Joao F. P.; Ramos, Helena M.;
López Jiménez, Petra Amparo; Costa Branco, P. J.. (2017) Electrical behaviour of
the pump working as turbine in off grid operation. Applied Energy (28)302-311.
10.1016/j.apenergy.2017.10.039
Carravetta, A.; Del Giudice, G.; Fecarotta, O.; Ramos, H. (2012). Energy production in
water distribution networks: A PAT design strategy. Water Resources Management
26(13):3947-3959.
Carravetta, A.; Del Giudice, G.; Fecarotta, O.; Ramos, H. (2013). PAT design strategy
for energy recovery in water distribution networks by electrical regulation. Ener-
gies 6(1):411-424.
Carravetta A, Conte MC, Fecarotta O, Ramos HM (2014a) Evaluation of PAT perfor-
mances by modified affinity law. Procedia Eng 89: 581-587.
doi:10.1016/j.proeng.2014.11.481
Carravetta, A.; Fecarotta, O.; Del Giudice, G.; Ramos, H. (2014b) Energy recovery in
water systems by PATs: A comparisons among the different installation schemes.
Procedia Eng 70: 275-284.
Carravetta A, Fecarotta O, Martino R, Antipodi L (2014c) PAT efficiency variation with
design parameters. Procedia Eng 70: 285-291. doi:10.1016/j.proeng.2014.02.032
Cassa, A.M.; Van Zyl, J. E. (2014) Predicting the leakage exponents of elastically
deforming cracks in pipes. In: 12th International Conference on Computing and
Control for the Water Industry, CCWI2013. Procedia Engineering 70, 302-310.
Castro, P.; Neves, M. (2003). Chlorine decay in water distribution systems case study–
Lousada network. Electronic Journal of Environmental, Agricultural and Food
Chemistry, 2, 261-266.
Castro, A. (2006). Minicentrales Hidroeléctricas; Instituto para la Diversificación y
Ahorro de la Energía: Madrid, España.
195
Bombas funcionando como Turbinas (PAT)
Caxaria, G.; Mesquita, D.; Ramos, H.M. Small Scale Hydropower: Generator Analysis
and Optimization for Water Supply Systems; World Renewable Energy Congress:
Linköping, Sweden, 2011; pp. 1386-1393.
Chapallaz, J., Eichenberger, P., Fischer, G. “Manual on Pumps Used as Tubines”;
MHPG Series; Vol. 11; Friedr," Vieweg&SohnVerlagsgesellschaftmbH, 1992.
Chen, J., Yang, H. X., Liu, C. P., Lau, C. H., Lo, M. (2013). A novel vertical axis water
turbine for power generation from water pipelines. Energy, 54, 184-193.
https://doi.org/10.1016/j.energy.2013.01.064
Choi, C. Y.; Shen, J. Y.; Austin, R. G. (2008). Development of a comprehensive solute
mixing model (AZRED) for double-tee, cross, and wye junctions. Water Distribution
Systems Analysis 2008; pp. 1-10.
Christodoulou, S.E.; Gagatsis, A.; Xanthos, S.; Kranioti, S.; Agathokleous, A.; Fragi-
adakis, M. (2013). Entropy-based sensor placement optimization for waterloss de-
tection in water distribution networks. Water Resources Management 27(13):
4443-4468.
Coelho, B.; Andrade‐Campos, A. (2014). Efficiency achievement in water supply sys-
tems—A review. Renew. Sustain. Energy Rev. 30, 59-84.
Colombo, A.F.; Karney, B.W. (2002). Energy and Costs of Leaky Pipes: Toward Com-
prehensive Picture. J. Water Resources Plannng and Management 128, 441-450.
Colombo, A.; Kleiner, Y. (2011). Energy recovery in water distribution systems using
microturbines. Probabilistic Methodologies. Water and Wastewater Engineering.
Toronto, Canada.
Corcoran, L.; Coughlan, P.; McNabola, A.; (2013) Energy recovery potential using
micro hydropower in water supply networks in the UK and Ireland. Water Science
and technology. Water Supply. 13 (2) 552-560; doi: 10.2166/ws.2013.050
Comisión Europea (2011). Commission Staff Working Paper. Executive Summary of
the Impact Assessment. Brussels, 30.11.2011. Consultado en Internet en agosto de
2017 en:
http://ec.europa.eu/research/horizon2020/pdf/proposals/horizon_2020_impact_asses
sment_report_executive_summary.pdf
Crespo Chacón M, Rodríguez Díaz J, García Morillo J, et al. Pump-as-Turbine Selec-
tion Methodology for Energy Recovery in Irrigation Networks: Minimising the
Payback Period. Water. 2019; 11(1):149; doi:10.3390/w11010149
da Conceição Cunha, M., & Ribeiro, L. (2004). Tabu search algorithms for water net-
work optimization. European Journal of Operational Research, 157(3),
746-758. https://doi.org/10.1016/S0377-2217(03)00242-X
196
Bibliografía
Davidson, J. W., & Bouchart, F. J.-C. (2006). Adjusting Nodal Demands in SCADA Con-
strained Real-Time Water Distribution Network Models. Journal of Hydraulic Engi-
neering, 132(1), 102-110. https://doi.org/10.1061/(ASCE)0733-9429(2006)132:1(102)
Deane, J. P., Gallachóir, B. Ó., McKeogh, E. J. (2010). Techno-economic review of
existing and new pumped hydro energy storage plant. Renewable and Sustainable
Energy Reviews, 14(4), 1293-1302.
Derakhshan, S., Nourbakhsh, A. (2008). Experimental study of characteristic curves of
centrifugal pumps working as turbines in different specific speeds. Experimental
Thermal and Fluid Sc. 32, (800-807).
Diario Oficial de las Comunidades europeas (2000). Directiva 2000/60/CE del Parla-
mento Europeo y del Consejo de 23 de octubre de 2000 por la que se establece un
marco comunitario de actuación en el ámbito de la política de aguas.
Consultado en agosto de 2017 en: http://www.mapama.gob.es/es/agua/temas/
planificacion-hidrologica/directiva2000_60marcoaguas_tcm7-28986.pdf
Dougherty, D. E., & Marryott, R. A. (1991). Optimal Groundwater Management: 1.
Simulated Annealing. Water Resources Research, 27(10), 2493-2508.
https://doi.org/10.1029/91WR01468
EcoGeneration (2008). Powering Sydney Waters carbon neutral future. News article.
July–August 2008.
Elbatran, A.; Yaakob, O.; Ahmed, Y.M.; Shabara, H. (2015). Operation, performance
and economic analysis of low head micro-hydropower turbines for rural and remote
areas: A review. Renewable and Sustainable Energy Reviews, vol. 43, pp. 40-50.
Estrada Tarragó, F. (2014). Micro-hydro solutions in Alqueva Multipurpose Project
(AMP) towards water-energy-environmental efficiency improvements. Universitat
Politècnica de Catalunya.
FAO. (2017) Aquastat. Consultado en septiembre de 2017 en
http://www.fao.org/nr/water/aquastat/data/query/results.html
Fecarotta, O.; Ariçò, C.; Carravetta, A.; Martino. R.; Ramos, H. (2015). Hydropower
Potential in Water Distribution Networks: Pressure Control by PATs. Water Re-
sources Management. February 2015, Volume 29, Issue 3, pp 699-714.
Fecarotta, O.; Carravett.a A.; Ramos, H.; Martino, R. (2016) An improved affinity
model to enhance variable operating strategy for pumps used as turbines. Journal
of Hydraulic Research 2016; 1686.1-10.
Fecarotta O, Ramos HM, Derakhshan S, Del Giudice G, Carravetta A.(2018) Fine
Tuning a PAT Hydropower Plant in a Water Supply Network to Improve System Ef-
fectiveness. J Water Resour Plan Manag.; 144(8):04018038.
doi:10.1061/(ASCE)WR.1943-5452.0000961
197
Bombas funcionando como Turbinas (PAT)
Fernandes, J F.P.; Pérez-Sánchez, M; Ferreira, F.; Ramos, HM.; López Jiménez, PA;
Costa Branco, P.J. (2019) Optimal energy efficiency of isolated PAT systems by
SEIG excitation tuning. Energy Conversion and Management (183) 391-405.
10.1016/j.enconman.2019.01.016
Fernandes J F.P.; Pérez-Sánchez, M; Ramos, HM.; Ferreira da Silva, F; López Jimé-
nez, PA; P.J. Costa Branco (2019). La Influencia de la Máquina Eléctrica en la Efi-
ciencia de las Bombas Trabajando Como Turbinas. EN XVI Seminario Iberoameri-
cano de Redes de Agua y Drenaje (SEREA 2019). (438-441). Lisboa, Portugal:
Dídia Covas, Laura Monteiro, Joana Carneiro (Eds).
Fontana, N.; Giugni, M.; & Portolano, D.; (2012). Losses reduction and energy pro-
duction in water-distribution networks. Journal of Water Resources Planning and
Management 138(3), 237–244.
Fouquet, D. (2012) Prices for renewable energies in Europe: Report 2011-2012. Tech-
nical Report, European Renewable Energies Federation.
Frijns, J.; Cabrera, E.; Carriço, N.; Covas, D.; Monteiro A. Ramos, H.; Bolognesi, A.;
Bragalli, C.; Baki, S.; Makropoulos, C. (2015) Management tools for hydro energy
interventions in water supply systems. Water practice and technology, 10 (2) 214-228.
Frosina, E.; Buono, D.; Senatore, A. (2017). A Performance Prediction Method for
Pumps as Turbines (PAT) Using a Computational Fluid Dynamics (CFD) Modeling
Approach. Energies 2017, 10, 103.
Gaius-Obaseki, T. (2010). Hydropower opportunities in the water industry. Interna-
tional Journal of Environmental Sciences 1(3), 392-402.
Gallagher, J.; Styles, D.; McNabola, A.; Williams, A.P. Life cycle environmental bal-
ance and greenhouse gas mitigation potential of micro-hydropower energy recovery
in the water industry. J. Clean. Prod. 2015, 99, 152-159.
Gao J, Yao, F.; , Xu, Y.; Sun, G.; Zheng, C.; Qi, S.; Cui, F. (2017).Analysis Model of
Physical Leakage Flow Based on Blind Source Separation Theory. Procedia Engi-
neering, Volume 186, 2017, Pages 269-277, ISSN 1877-7058.
Geldreich, E.E.; LeChevallier, M. (1999) Microbiological quality control in distribution
systems. In: Water Quality and Treatment: A Handbook of Community Water Sup-
plies. Fifth Edition . Letterman, R.D. (Ed.) McGraw-Hill, New York. pp.18.1-18.49.
Generalitat Valenciana (2016). RIS 3. (Research and Innovation Smart Specialisation
Strategies) Estrategia de Especialización Inteligente para la Investigación e Inno-
vación en la Comunitat Valenciana. Consultado en septiembre de 2017 en:
http://www.redidi.es/sites/default/files/biblioteca-documentos/ris3cv.pdf
198
Bibliografía
Ghiassi, M., Zimbra, D. K., & Saidane, H. (2008). Urban Water Demand Forecasting
with a Dynamic Artificial Neural Network Model. Journal of Water Resources
Planning and Management, 134(2), 138-146. https://doi.org/10.1061/(ASCE)0733-
9496(2008)134:2(138)
Giugni, M., Fontana, N., & Ranucci, A. (2014). Optimal Location of PRVs and Turbines
in Water Distribution Systems. Journal of Water Resources Planning and Manage-
ment, 140(9), 06014004. https://doi.org/10.1061/(ASCE)WR.1943-5452.0000418
Gonçalves, F. V.; Ramos, H. M. (2012). Hybrid energy system evaluation in water
supply systems: artificial neural network approach and methodology. Water Sup-
ply: Research and Technology – AQUA 61(2), 59-72.
Grover K (1980) Conversion of pumps to turbines. In: GSA Inter corp. Katonah, New
York
Grundfos (2003). Grundfos Motor Book. Available at:
http://www.grundfos.com/content/dam/Global%20Site/Industries%20%26%20solut
ions/waterutility/pdf/motorbook.pdf[4] Electric Machinery, A.E. Fitzgerald, C.
Kingsley, S.D. Umans, New York: McGraw-Hill, 2003 (sixth edition).
Grundfos (2004) Industry, Grundfos. "Pump Handbook." (2004).
Gubaidullin, (2019) “Magnetic approach to improve the efficiency of water treatment
plants”, Askar Gubaidullin, SKF Magnetic Mechatronics, Vernon, France.
Available from the internet: <http://evolution.skf.com/magnetic-approach-to-
improve-the-efficiency-of-watertreatment-plants/>
Hall, C.; Harakan, J.; Hallock, J.; Cleveland, C.; Jefferson, M. (2003). Hydrocarbons
and the evolution of human culture. NATURE, nº 426, pp. 318-320.
Harrison, R., Hau, E., & Snel, H. (2000). Large wind turbines: design and econom-
ics (Vol. 1). Chichester: Wiley.
Hernandez-Leal, L.; Temmink, H.; Zeeman, G.; Buisman, C. (2010). Bioflocculation of
grey water for improved energy recovery within decentralized sanitation concepts.
Bioresource Technology 101(23), 9065–9070.
Horizonte 2020 (2017). Es Horizonte 2020. COST. Consultado en agosto de 2017 en:
http://www.eshorizonte2020.es/mas-europa/otros-programas/cost
Iglesias-Rey, P., Izquierdo, J., Fuertes, V., & Martínez-Solano, F., 2004. Modelación
de transitorios hidráulicos mediante ordenador. Grupo Mult., Universitat Politècnica
de València.
199
Bombas funcionando como Turbinas (PAT)
200
Bibliografía
Kumar, S.; Forand, S.; Babcock, G.; Richter, W.; Hart, T.; Hwang, S. A. (2014). Total
trihalomethanes in public drinking water supply and birth outcomes: a cross-
sectional study.
Kumar, D., Katoch, S. 2015. Small hydropower development in western Himalayas:
Strategy for faster implementation. Renewable Energy, vol. 77, pp. 571-578.
Kuo, S.-F., Liu, C.-W., & Chen, S.-K. (2003). Comparative study of optimization
techniques for irrigation project planning. Journal of the American Water Resources
Association, 39(1), 59–73. https://doi.org/10.1111/j.1752-1688.2003.tb01561.x
Latchoomun, L.; Ah King, R.; Busawon, K.; Mawooa, D.; Kaully, R.; (2015). Labora-
tory Investigation of the Leakage Characteristics of Unburied HDPE Pipes, Proce-
dia Engin., Volume 119, 2015, 91-100.
Lambert A. (2001) What do we know about pressure: Leakage relationships in distri-
bution systems?. IWA Conference System Approach to Leakage Control and Water
Distribution System Management, Brno, República Checa.
Le Chevallier M. W.; Gullick R. W.; Karim, M. R.; Friedman M.; Funk J. E.(2003).
The Potential for health risks from intrusion of contaminants into the distribution
systems from pressure transients, Journal of Water and Health 01.1. IWA Publish-
ing 2003. pp 3-14.
Li, C.; Wang, D.; Xu, X.; Wang, Z. (2017).Formation of known and unknown disinfec-
tion by-products from natural organic matter fractions during chlorination, chloram-
ination, and ozonation. Science of The Total Environment, Vol. 587. pp 177-184.
Liu, H.; Yuan, Y.; Zhao, M.; Zheng, X.; Lu, J.; & Zhao, H. (2011). Study of Mixing at
Cross Junction in Water Distribution Systems Based on Computational Fluid Dy-
namics. In ICPTT 2011: Sustainable Solutions For Water, Sewer, Gas, And Oil
Pipelines (pp. 552-561).
Lydon, T.; Coughlan, P. J. Hydraul. Eng., McNabola, A. (2017). Pump-As-Turbine:
Characterization as an Energy Recovery Device for the Water Distribution Net-
work, 143(8): 04017020.
Lofman, D., Petersen, M. & Bower, A. (2002). Water, energy and environment nexus:
The California experience. Water Resources Development 18(1), 73–85.
López Jiménez, P.A.; Vela, A.F.; Iglesias, P. (1995) Análisis de seguridad en abaste-
cimientos y su aplicación a la detección de fugas. Ingeniería del agua, 3 (2), 7-22.
López Jiménez, P.A.; Martínez Solano, F.J.; López Patiño, G.; Lara Ledesma, B.
(2007) Modelación computacional del fenómeno de una fuga en tubería de abaste-
cimiento. Ingeniería Hidráulica en México, 2 (22), 43–54.
201
Bombas funcionando como Turbinas (PAT)
López-Jiménez, P.A.; Mora Rodriguez, J.J.; Martínez Solano, F.J.; Pérez García, R.
(2008) Modelación hidrodinámica del fenómeno de intrusión en tubería de abaste-
cimiento. Ingeniería Hidráulica en México, 4 (XXIII), 103-110.
Maestu, J.; Villar, A. 2007.Precios Y Costes de Los Servicios Del Agua En España,
Madrid, Spain, Consultado en septiembre de 2016 en:
http://hispagua.cedex.es/sites/default/files/especiales/Tarifas_agua/precios_costes_s
erviciosagua.pdf
MAGRAMA (2014). El riego localizado alcanza el 48,23% de la superficie regada en
España. Consultado en Internet en septiembre de 2017 en:
http://www.mapama.gob.es/fr/prensa/noticias/el-riego-localizado-alcanza-el-4823--
de-la-superficie-regada-en-espa%C3%B1a-/tcm12-312671-16
Malavasi, S., Ferrarese, G., & Rossi, M. M. A. (2014). A control valve for energy harvest-
ing. Procedia Engineering, 89, 588–594. https://doi.org/10.1016/j.proeng.2014.11.482
Malavasi, S., Rossi, M. M. A., & Ferrarese, G. (2018). GreenValve: hydrodynamics
and applications of the control valve for energy harvesting. Urban Water Journal,
15(3), 200–209. https://doi.org/10.1080/1573062X.2016.1175483
Mataix C (2009) Turbomáquinas Hidráulicas. Universidad Pontificia Comillas, Madrid.
May J. (1994). Pressure dependent leakage. World Water and Environmental Engi-
neering, p 10, 1994.
McNabola, A., Williams, A. P. & Coughlan, P. (2013). Energy recovery in water sup-
ply networks: an assessment of the potential of micro hydropower. Water and Envi-
ronment Journal 27, 435–436.
McNabola, A.; Coughlan, P.; Corcoran, L,; Power C.; Williams,A.P.; Gallagher,H.J.;
Styles, D. (2014). Energy recovery in the water industry using micro-hydropower:
an opportunity to improve sustainability. Water Policy 16 (2014) 168–183.
McNabola, A.; Couglan, P ; Williams, A.P. (2011) The Technical & Economic Feasi-
bility of Energy Recovery in Water Supply Networks RE&PQJ, Vol.1, No.9, May .
Moreno, M., Córcoles, J., Tarjuelo, J., Ortega, J. (2010). Energy efficiency of pressurised
irrigation networks managed on-demand and under a rotation schedule. Biosystems
Engineering, 107(4), 349–363. https://doi.org/10.1016/j.biosystemseng.2010.09.009
Morabito A, Steimes J, Bontems O, Al Zohbi G, Hendrick P. 2017. Set-up of a pump
as turbine use in micro-pumped hydro energy storage: A case of study in Froyennes
Belgium. J. Phys.: Conf. Ser. vol. 813.
Morabito A, Hendrick P. 2019 Pump as turbine applied to micro energy storage and
smart water grids: A case study. Appl Energy. 241:567-579.
doi:10.1016/J.APENERGY.2019.03.018
202
Bibliografía
203
Bombas funcionando como Turbinas (PAT)
Pasten, P., Santamarina, J.C. (2012). Energy and quality of life. Energy Policy, 49,
pp. 468-476.
Pardo, M. A., Manzano, J., Cabrera, E., García-Serra, J. (2013). Energy audit of
irrigation networks. Biosystems Engineering, 115(1), 89–101.
https://doi.org/10.1016/j.biosystemseng.2013.02.005
Pérez-Sánchez, M. (2017) Methodology for Energy Efficiency Analysis in Pressurized
Irrigation Networks, Practical Application. 2017. Available online:
https://riunet.upv.es/bitstream/handle/10251/84012/RESUMEN. pdf?sequence=3
Pérez-Sánchez, M.; López-Jiménez, P.A.; Ramos, H. M. (2018a). Modified Affinity
Laws in Hydraulic Machines towards the Best Efficiency Line. Water Resources
Management., doi. 10.1007/s11269-017-1841-0
Pérez-Sánchez, M.; Sánchez Romero, F.J.; López-Jiménez, P.A.; Ramos, H. (2018b).
PATs selection towards sustainability in irrigation networks: simulated annealing
as a water management tool. Renewable Energy. 10.1016/j.renene.2017.09.060.
Pérez-Sánchez, M; Ferreira, A.R.; López Jiménez, PAo; Ramos, HM.. (2018c) Design
strategy to maximize recovery energy towards smart water grids: case study. Urban
Water Journal, 4 (15), 329 - 337. 10.1080/1573062X.2018.1459747
Pérez-Sánchez, M; Sánchez-Romero, FJ; López Jiménez, PA. (2018d). A Case Study of
Energy recovery systems location in irrigation networks defined by simulated an-
nealing techniques. Simulated Annealing. Introduction, Applications and Theory
(199 - 220). Nova Science Publishers. 978-1-53613-674-6
Pérez-Sánchez, M. (2017). Methodology for energy efficiency improvement analysis in
pressurized irrigation networks. Practical application. López Jiménez, P.A (codir.).
Tesis doctoral. UPV. Mayo 2017.
Pérez Sánchez, Modesto; Sánchez Romero, F.J.; Ramos, H-M.; López Jiménez, P.A.
(2017a). Optimization Strategy for improving the energy efficiency in irrigation sys-
tems by micro hydropower: Practical Application. Water 2017, 9, 799;
doi:10.3390/w9100799.
Pérez-Sánchez, M.; Sánchez-Romero, F.J.; Ramos, H. M.; López Jiménez, P. A.
(2017b) Energy Recovery in Existing Water Networks: Towards Greater Sustaina-
bility. Water, (9), 97. doi:10.3390/w9020097
Pérez-Sánchez, M.; Sánchez-Romero, F.J.; López-Jiménez, P. A. (2017c). Huella
energética del agua en función de los patrones de consumo en redes de distribu-
ción. Ingeniería del agua, v. 21, n. 3, p. 197-212, jul. 2017.
204
Bibliografía
205
Bombas funcionando como Turbinas (PAT)
Ramos, H.M.; Borga, A. (2000). Pumps yielding power. Dam Engineering. 10(4), 197–217
Ramos, H.M.; Borga, A. (1999). Pumps as turbines: an unconventional solution to ener-
gy production. Urban Water, 1(3):261-263. doi:10.1016/S1462-0758(00)00016-9
Ramos, H.M.; Mello,P. (2010). Clean power in water supply systems as a sustainable
solution: from planning to practical implementation. Water Science & Technology:
Water Supply, vol. 10, nº 1, pp. 39-49.
Ramos, H.M.; Simão, M.; Borga, A.(2013) Experiments and CFD Analyses for a New
Reaction Microhydro Propeller with Five Blades. J Energy Eng. 2013;139(2):109-117.
Ramos, H.M.; Kenov, K. N.; Vieira, F. (2011). Environmentally friendly hybrid solu-
tions to improve the energy and hydraulic efficiency in water supply systems. Ener-
gy for Sust. Management 15(4), 436–442.
Ramos, H.M.; Almeida, A.B. Dynamic effects in micro-hydro modelling. Int. Water
Power Dam Constr. 2003, 55, 22–25.
Ramos, H.M.; Amaral, M.P.; Covas, D. (2014). Pumped-storage solution towards
energy efficiency and sustainability: Portugal contribution and real case studies.
Journal of Water Resource and Protection, 6(12), 2014.
Ramos, H. M., Kenov, K. N. & Vieira, F. (2011). Environmentally friendly hybrid
solutions to improve the energy and hydraulic efficiency in water supply systems.
Energy for Sust. Man. 15(4), 436–442.
Ramos, J.; Ramos, H. M. (2009). Sustainable application of renewable sources in wa-
ter pumping systems: Optimized energy system configuration. Energy Journal 37(2),
633–643.
Ramos, H.M.; López Jiménez, P. Amparo; Teyssier, C. (2013) Optimization of Reten-
tion Ponds to Improve the Drainage System Elasticity for Water-Energy Nexus.
Water Resources Man.t (27)2889 - 2901.
Rawal, S., Kshirsagar, J. T., (2007). Numerical Simulation on a Pump Operating in a
Turbine Mode. Proceedings of the Twenty-Third Int. Pump Users Symposium.
USA, Texas AM University.
Razan, J.I.; Islam, R.S.; Hasan, R.; Hasan, S.; Islam, F. A (2012). Comprehensive
Study of Micro-Hydropower Plant and Its Potential in Bangladesh. Renew. Energy
2012, 2012, 635396.
Reca, J., Martínez, J., Gil, C., & Baños, R. (2008). Application of Several Meta-
Heuristic Techniques to the Optimization of Real Looped Water Distribution Net-
works. Water Resources Management, 22(10), 1367–1379.
https://doi.org/10.1007/s11269-007-9230-8
206
Bibliografía
Renault, D.; Facon, T.; Wahaj, R. (2007). Modernizing Irrigation Management: The
MASSCOTE Approach-Mapping System and Services for Canal Operation
Techniques. Paper, FAO Irrigation and Drainage, vol. 63, Roma, p. 177.
Rodríguez, M.,Rodriguez, G.; Serodes, J.; Sadiq, R. (2007). Subproductos de la desin-
fección del agua potable: formación, aspectos sanitarios y reglamentación. Inter-
ciencia 32.11 (2007): 749-756.
Romero, Laura; Pérez-Sánchez, Modesto; López Jiménez, P. Amparo. (2017) Envi-
ronmental and energy problematic in the Mediterranean irrigation regions frame-
work. International Journal of Energy and Environment (IJEE), 1 (8), 51 - 62.
Romero L.; Pérez-Sánchez, M.; López Jiménez, P.A. (2017) Improvement of sustaina-
bility indicators when traditional water management changes: a case study in Ali-
cante (Spain). AIMS Environmental Science, 3 (4), 502 - 522.
Romero-Gómez, P., Lansey, K. E., & Choi, C. Y. (2011). Impact of an incomplete
solute mixing model on sensor network design. Journal of Hydroinformatics, 13 (4),
642-651.
Rossman, L. A. (2000). EPANET 2: Users manual. Environmental Protection Agency
(EPA), EEUU.
Rygaard, M., Binning, P. J. & Albrechtsen, H.-J. (2011). Increasing urban water self-
sufficiency: New era, new challenges. Journal of Environmental Management 92(1),
185–194. Ramos, J. S. & Ramos, H. M. (2010). Multi-criterion optimization of en-
ergy management in drinking systems. Water Science & Technology: Water Supply
10(2), 129–144.
Saket, R.K. (2008) Design, development and reliability evaluation of micro hydro
power generation system based on municipal waste water. IEEE Elect. Power En-
erg. Conf. Energ. Innovat.
Samora, I.; Franca, M.; Schleiss, A.; Ramos, H.M. (2016) Simulated Annealing in Op-
timization of Energy Production in a Water Supply Network. Water Resour. Manag.
2016, 30, 1533–1547.
Samora, I.; Manso, P.; Franca, M.J.; Schleiss, A.J.; Ramos, H.M (2016b). Opportunity
and Economic Feasibility of Inline Microhydropower Units in Water Supply Net-
works. J. Water Resour. Plan. Manag. 2016, 142, 04016052.
Samora, I.; Manso, P.; Franca, M.; Schleiss, A.; Ramos, H. (2016c) Energy Recovery
Using Micro-Hydropower Technology in Water Supply Systems: The Case Study
of the City of Fribourg. Water, 8, 344–360.
San Bruno, G., Choulot, A., Denis, V. (2010). Energy recovery in existing infrastruc-
tures with small hydropower plants. Multipurpose Schemes–overview and Exam-
ples. ESHA, SHAPES & Mhylab.
207
Bombas funcionando como Turbinas (PAT)
Sánchez Zafra, A. (2008). Efectos de los trihalometanos sobre la salud, Hig. Sanid.
Ambient. Pp. 8:280-290.
Sanz, G., Pérez, R. (2014). Demand Pattern Calibration in Water Distribution Networks.
Procedia Engineering, 70, 1495–1504. https://doi.org/10.1016/j.proeng.2014.02.164
Scollen, A., Hargreaves, T. (2018). Simulated Annealing: Introduction, Applications
and Theory. New York: Nova Science Publichsers. 978-1-53613-674-6
Seoane, P.; Allué, R.; Postigo, M.J.; Cordón, M.A. (2013). Boletín Mensual de Estadís-
tica; Instituto Nacional de Estadística: Madrid, Spain, 2013.
Simao, M.; Ramos, H.M. (2010). Hydrodynamic and performance of low power tur-
bines: Conception, modelling and experimental tests. Int. J. Energy Environ. 2010,
1, 431–444.
Simao, M; Pérez-Sánchez, Modesto; Carraveta, Armando; López Jiménez, Petra Am-
paro; Ramos, Helena M. (2018) Velocities in a Centrifugal PAT Operation: Expe-
riments and CFD Analyses. Fluids, 1 (3), 1 - 21. 10.3390/fluids3010003
Simao, M; Pérez-Sánchez, M; Carraveta, A; Ramos, HM. (2019) Flow Conditions for
PATs Operating in Parallel: Experimental and Numerical Analyses. Energies, 91
(12), 1-19. 10.3390/en12050901
Simpson AR, Marchi A (2013). Evaluating the Approximation of the Affinity Laws and
Improving the Efficiency Estimate for Variable Speed Pumps. J Hydraul Eng
139:1314–1317. doi:10.1061/(ASCE)HY.1943-7900.0000776
Singh, P. (2005). Optimization of the internal hydraulic and of system design in pumps
as turbines with field implementation and evaluation.
Singh, P., Nestmann, F. (2010). An optimization routine on a prediction and selection
model for the turbine operation of centrifugal pumps. Experimental Thermal and Flu-
id Science, 34(2), 152–164. https://doi.org/10.1016/j.expthermflusci.2009.10.004
Sitzenfrei, R., Leon, J.V. (2014). Long-time simulation of water distribution systems
for the design of small hydropower systems. Renewable Energy, vol. 72,
pp. 182-187.
Song, I.; Romero-Gómez, P.; & Choi, C. Y. (2009). Experimental verification of in-
complete solute mixing in a pressurized pipe network with multiple cross junctions.
Journal of Hydraulic Engineering, 135(11), 1005-1011.
Thin, K.C.; Khaing, M.; Aye, K.M. Design and performance analysis of centrifugal
pump. World Acad. Sci. Eng. Technol. China 2008, 46, 422–429
Tilmant, A., Goor, Q, Pinte, D. (2009). Agricultural-to-hydropower water transfers:
sharing water and benefits in hydropower-irrigation systems. Hydrology and Earth
System Sciences, vol. 13, pp. 1091-1101.
208
Bibliografía
Tospornsampan, J., Kita, I., Ishii, M., & Kitamura, Y. (2007). Split-Pipe Design of
Water Distribution Network Using Simulated Annealing. World Academy of Sci-
ence, Engineering and Technology, International Journal of Environmental, Chemi-
cal, Ecological, Geological and Geophysical Engineering, 1(4), 28–38.
Tucciarelli T, Criminisi A, Termini D (1999) Leak Analysis in Pipeline Systems by
Means of Optimal Valve Regulation. J Hydraul Eng-ASCE 125(3):277–285.
UPV (2017a). Reglamento para la evaluación de la actividad de investigación, desarrollo
e innovación en la Universitat Politècnica de València. Consultado en agosto de 2017
en: www.upv.es/entidades/VIIT/menu_urlc.html?/entidades/VIIT/info/U0569617.pdf
UPV (2017b). Valoraciones anuales de la actividad Investigadora por estructuras y
Grupos de Investigación. Universitat Politècnica de València. Consultado en octu-
bre de 2017 en: http://www.upv.es/entidades/VIIT/info/216020normalc.html
Van Bloemen Waanders, B.; Hammond, G.; Shadid, J.; Collis, S.; & Murray, R.
(2005). A comparison of Navier-Stokes and network models to predict chemical
transport in municipal water distribution systems. In Proc., World Water and Envi-
ronmental Resources Congress. Anchorage.
Van Zyl, J.E. (2014) Theoretical Modeling of Pressure and Leakage in Water Distribu-
tion Systems, Procedia Engineering, Volume 89, pp. 273-277.
Verdegay, J. L., Yager, R. R., Bonissone, P. P. (2008). On heuristics as a fundamental
constituent of soft computing. Fuzzy Sets and Systems, 159(7), 846–855.
https://doi.org/10.1016/j.fss.2007.08.014
Vicente, S.; Bludszuweit, H. (2012). Flexible design of a pico-hydropower system for
Laos communities. Renewable Energy, vol. 44, pp. 406-413.
Vidal, R.; Gallardo, A.; Martinez, F.; Lopez, P. (1996) Modelización de subproductos
de la desinfección en las redes de distribución de agua potable. Implicaciones
hidráulicas. Tecno Bomb., 1 (1), 32-38.
Vieira, F.; Ramos, H. (2009a). Optimization of operational planning for wind/hydro
hybrid water supply systems. Renewable Energy 34(3), 928–936.
Vieira, F.; Ramos, H. M. (2009b). Optimization of the energy management in water
supply systems. Water Science & Technology: Water Supply 9(1), 59–65.
Vilanova, M.R.N.; Balestieri, J.A.P.; Filho, P.M. (2013) Energy and hydraulic efficien-
cy in conventional water supply systems. Energy Policy 2013, 30, 701–714
Walski, T.; Bezts, W.; Posluzny, E.; Weir, M.; Withman, B. (2006) Modeling leakage
reduction through pressure control. J Am Water Works Ass: 148–155
209
Bombas funcionando como Turbinas (PAT)
210