Geração de renda e ocupações na cadeia agroindustrial de castanha de caju no Rio Grande do Norte: o caso do sertão de Apodi

This work is part of the research group studying agribusiness chains, by characterizing agroindustrial chain of cashew nuts in Rio Grande do Norte, with special emphasis on Sertão de Apodi Citizenship Territory, showing their different flows and linkages and analyzes generating jobs and income in the production of this good, from the experience of COOPAPI, using as a methodology for analyzing the input-output matrix. And this is justified by the fact that Brazil is in the context as the world's fifth largest producer of cashew nuts and the production of this product focus in the Northeast, especially the Rio Grande do Norte, which is the second largest producer of region, behind only the state of Ceará, which is the largest producer in the country. Moreover, in recent years the state has increased its insertion in international markets this product, particularly the Mossoró-Assu Territory, which constitutes the largest producer in the state, particularly the municipality of Serra do Mel, which has an arrangement different yield, based on the model of integrated rural development, embodied in family farming. Nevertheless, the Sertão de Apodi Citizenship appear as the fourth among the Citizenship territories of Rio Grande do Norte in the ranking of the production of cashew nuts, but the prevailing model of Cooperatives in the Territory, based on the Solidarity Economy, i,e, different from the old cooperative, which is based on the precepts as patronage, with cooperative owners, favor, etc.., encourages scientific research. In this, after a careful analysis of data from field research, it was noticed that among the links in the chain, the production link is what generates the most jobs, while the processing link commands to generate income, but a portion meaningful occupations and income are generated outside the territory, especially out of state RN, denouncing the gap in supply of some specific demands within the chain

Controle de estoques de sal: uma inovação da quantificação e no gerenciamento dos recursos naturais da indústria mineral

The strengthening of the domestic industry in Brazil required the modernization, mechanization and expansion of salt production. Thereafter the production of sea salt started to be made in a process of continuous flow, where the product is constantly stored in yards, with daily movements in and out of salt. Thus far, the major bottleneck found in this production process is the control of production, because due to the large amount produced and variety of losses existing in the various stages of production there are not a regulated and safe way to control inventories with accuracy and speed demanded. In a typical case with a salt marsh company of Rio Grande do Norte state, salt produced is stored in two open courtyards and inventory control of salt made by carrying input / output relationship of salt in each storage yard. This work developed a conceptual model of inventory control, based on topography, adopting surveys into one of the courtyards of the company. There were 25 biweekly survey measurements over a year book to generate digital models representing the stock. For each measurement, results were compared with the values of inventory accounting provided by the salt marsh in order to identify existing losses and mark out the sales department on the actual stock available at each measurement date. Inventories calculated by the model indicated losses of 6,349 tonnes for the period of one year book and 3,279 tonnes for the period between harvests, when compared to the accounting control

O emprego na transição para uma economia verde: cenários para o Brasil com base no consumo de energia pela indústria

Pós-graduação em Economia - FCLAR

Localização de falhas em linhas de transmissão utilizando decomposição harmônica e redes neurais artificiais

Este trabalho propõe a utilização de uma nova metodologia para a localização de falhas em linhas de transmissão (LT). Esta metodologia consiste na utilização da decomposição harmônica da corrente de fuga de uma linha e na aplicação de uma Rede Neural Artificial (RNA) capaz de distinguir padrões da condição normal de funcionamento e padrões de situações de falhas de uma LT. Foi utilizado um modelo Pi capaz de absorver dados reais de tensão e corrente de três fases e de alterar valores de R, L e C segundo modificações ambientais. Neste modelo foram geradas falhas em todas as torres com diferentes valores de capacitância. A saída fornecida pelo modelo é a decomposição da corrente de fuga do trecho considerado. Os dados de entrada e saída do modelo foram utilizados no treinamento da RNA desenvolvida. A aquisição de dados reais de tensão e corrente foi feita através de analisadores de parâmetros de qualidade de energia elétrica instalados nas extremidades de um trecho de LT, Guamá-Utinga, pertencente à Centrais Elétricas do Norte do Brasil ELETRONORTE. O cálculo dos parâmetros construtivos foi feito através do método matricial e melhorado através da utilização do Método de Elementos Finitos (MEF). A RNA foi desenvolvida com o auxílio do software Matlab. Para treinamento da RNA foi utilizado o algoritmo de Retropropagação Resiliente que apresentou um bom desempenho. A RNA foi treinada com dois conjuntos de dados de treinamento para analisar possíveis diferenças entre as saídas fornecidas pelos dois grupos. Nos dois casos apresentou resultados satisfatórios, possibilitando a localização de falhas no trecho considerado.

Evolução agrária do Nordeste paraense por contas sociais ascendentes: a inclusão das redes mercantis de insumos nos processos metodológicos e o cálculo dos níveis de endogenização

A importância deste trabalho está em analisar a dinâmica da evolução econômica do Nordeste Paraense, por suas estruturas fundamentais de produção no agrário, interpretando-se os resultados por meio dos impactos setoriais e os transbordamentos da renda e do produto na economia. Para tanto, aplicou-se o modelo de Contas Sociais Ascendentes de Base Agrária (Matriz de Insumo-Produto) metodologicamente mais abrangente que as primeiras abordagens aplicadas em trabalhos anteriores realizados por Pires (2002) e André (2004) a mesma região, todas derivadas de Costa (2002abcd), inclusive este trabalho. Entre outras diferenciações, o modelo desenvolvido inicialmente estava ligado internamente pelo fluxo de produtos entre as dimensões geográficas do local, ao estadual e nacional, em que o valor da produção é formado no sentido que vai da produção até o consumo final, ou ainda, do rural para o urbano: segue o sentido que vai da produção local até a demanda final nacional. Para este trabalho, a forma metodologicamente adotada, já apresentada por Costa em 2008, no Sudeste Paraense, incorpora, primordialmente, as cadeias backward de insumos e investimentos em que são reconstituídos os trajetos dos agentes produtores até o “primeiro fornecedor”, tanto para as redes de formação de preço, quanto de salários e margem bruta (markup), o que possibilita, uma análise diferenciada dos multiplicadores quanto aos impactos setoriais, implicando nas análises de retenções e transbordamentos da renda e do produto entre as formas de produção e suas ocorrências estruturais geográficas, na verificação das disposições das forças centrípetas e centrífugas, o que permite fazer interpretações quanto às possibilidades de conexões econômicas locais para o desenvolvimento, a partir dos índices de aglomeração.

Economia dos produtos florestais não-madeireiros no estado do Amapá: sustentabilidade e desenvolvimento endógeno

O extrativismo vegetal, baseado na exploração sustentável dos produtos florestais não madeireiros é uma das alternativas mais consistentes ecologicamente, no que tange a Conservação da biodiversidade e cobertura natural da Floresta Amazônica. Mas, há problemas sistêmicos que não têm deixado este segmento desenvolver-se de forma satisfatória. O vazamento da economia regional dos produtos extrativistas da Amazônia é uma variável que contribui para o baixo nível de desenvolvimento social e econômico da região. É imperativo que sejam realizadas análises econômicas sobre as possibilidades desse novo modelo A solução provisória para esta área foi estabelecida em 1841, com a definição do espaço entre os rios Oiapoque e Araguari como “Área do Contestado”. Ficou definido então, que a gestão da referida área seria feita de forma compartilhada entre os dois Países. A questão da disputa entre Brasil e França pelas terras do Amapá, só se resolveu definitivamente com a intermediação do presidente suíço Walter Hauser, que em 1o de dezembro de 1900, através do protocolo conhecido como “Laudo Suíço”, concedeu a referida área ao Brasil. de desenvolvimento, apoiado nas bases de um capitalismo contemporâneo, que tem procurado aprimorar os mecanismos econômicos de inclusão social de variáveis ambientais como elementos endógenos do sistema produtivo na Amazônia. O Amapá é o estado brasileiro mais preservado (97% de sua cobertura florestal original) e possui 72% do território como áreas protegidas. Este trabalho pretende identificar o nível de contribuição que os produtos não-madeireiros, extraídos nas florestas por populações tradicionais têm na economia do Amapá. Como método central de análise, foi utilizado o modelo de matriz insumo-produto desenvolvido por Wassily Leontief. Este procedimento analítico tem relevante destaque como instrumento prático de análise e planejamento econômico. Além das análises de programação do crescimento econômico setorial, é adequado para estimar, mediante os efeitos multiplicadores, os impactos do crescimento econômico na produção trabalho e renda setorial de toda a economia. O objeto central do trabalho são as análises estruturais dos arranjos produtivos locais dos produtos florestais não-madeireiros do estado do Amapá, em níveis regional e local. O método desenvolvido pelo Francisco de Assis Costa (Contas Sociais Ascendentes Alfa - CS_) foi usado na construção das matrízes, tendo como base os procedimentos analíticos ascendentes de agregação progressiva dos dados locais O setor extrativista dos produtos florestais não -madeireiros do estado do Amapá, com um PIB de R$ 204 milhões no ano de 2009 tem uma participação pequena, com somente 3,07% das riquezas produzidas no estado. Todavia, quando é observado que 86,77% do valor total do PIB (6,65 bilhões de reais) é decorrente de atividades do Setor Terciário e que todo o Setor Produtivo (Primário e Secundário), tem participação de apenas 13,24%, conclui-se que o valor da produção dos produtos florestais não-madeireiros é bem significativo no estado do Amapá.

Políticas de contenção de desmatamento, produção e mercado de terras na Amazônia: um ensaio sobre a economia local do sudeste paraense usando contas sociais alfa (CSα)

A discussão atual sobre a emissão de carbono associada ao uso agropecuário da terra em prejuízo de florestas se ressente de uma visão sistêmica, no que se refere aos fluxos econômicos propriamente, e suas interações, no que tange ao ambiente institucional que os garante. Dado que os esquemas de compensação implicam a entrada e saída de recursos em contextos econômicos amplos e sistêmicos, fundamental é discutir qual o resultado final desses fluxos sobre as condições gerais de reprodução das economias locais. As questões básicas são: a) como tais políticas poderão, a partir dos setores rurais, afetar a demanda final efetiva e, por essa via, o valor da produção e as variáveis de valor adicionado de toda a economia e b) como as variações na economia afetam as formas de uso da base natural e, portanto, o desmatamento. No que se refere às instituições, o artigo dá especial ênfase às que definem o mercado de terras, porque nele encontra o cerne de questões vitais para o que se discute. O artigo utiliza um modelo ascendente de geração de matrizes de insumo-produto para a economia local do Sudeste Paraense, incorpora nela um balanço de carbono dos setores da produção rural, encontra os multiplicadores e simula quatro situações de politica de contenção de desmatamento e redução das emissões de gases poluentes. A conclusão principal do artigo é que se faz necessário pensar políticas de contenção de desmatamento ligadas indissociavelmente a políticas de produção – a serem operadas por mecanismos que façam convergir as decisões dos agentes com perspectivas macro de desenvolvimento: local, endógeno e sustentável.

Fragmentação geopolítica: uma análise estrutural das propostas de separatismo no Brasil

Este trabalho tem por objetivo dar luz ao debate sobre a criação dos novos entes federativos, trazendo elementos objetivos sobre a caracterização socioeconômica dos novos espaços geopolíticos e seus respectivos papéis no novo contexto espacial do sistema inter-regional brasileiro. O artigo traz resultados inéditos para a discussão, buscando, inicialmente, identificar padrões hierárquicos e de dependência espacial e produtiva neste novo contexto federativo. Tais resultados subsidiama analise subsequente dos impactos da nova configuração das transferências constitucionais, que identifica pormeio de simulações comummodelo inter-regional de insumo-produto especialmente calibrado para as 33 regiões consideradas no estudo, não apenas os potenciais ganhadores e perdedores, mas também os mecanismos de interação espacial subjacentes a estes processos.

Advantages and Applications of Transforming Empirical Model Input Space with Dimensional Models

Model-based calibration of steady-state engine operation is commonly performed with highly parameterized empirical models that are accurate but not very robust, particularly when predicting highly nonlinear responses such as diesel smoke emissions. To address this problem, and to boost the accuracy of more robust non-parametric methods to the same level, GT-Power was used to transform the empirical model input space into multiple input spaces that simplified the input-output relationship and improved the accuracy and robustness of smoke predictions made by three commonly used empirical modeling methods: Multivariate Regression, Neural Networks and the k-Nearest Neighbor method. The availability of multiple input spaces allowed the development of two committee techniques: a 'Simple Committee' technique that used averaged predictions from a set of 10 pre-selected input spaces chosen by the training data and the "Minimum Variance Committee" technique where the input spaces for each prediction were chosen on the basis of disagreement between the three modeling methods. This latter technique equalized the performance of the three modeling methods. The successively increasing improvements resulting from the use of a single best transformed input space (Best Combination Technique), Simple Committee Technique and Minimum Variance Committee Technique were verified with hypothesis testing. The transformed input spaces were also shown to improve outlier detection and to improve k-Nearest Neighbor performance when predicting dynamic emissions with steady-state training data. An unexpected finding was that the benefits of input space transformation were unaffected by changes in the hardware or the calibration of the underlying GT-Power model.

X-parameters based analytical design of non-linear microwave circuits Application to oscillator design

Resumen El diseño clásico de circuitos de microondas se basa fundamentalmente en el uso de los parámetros s, debido a su capacidad para caracterizar de forma exitosa el comportamiento de cualquier circuito lineal. La relación existente entre los parámetros s con los sistemas de medida actuales y con las herramientas de simulación lineal han facilitado su éxito y su uso extensivo tanto en el diseño como en la caracterización de circuitos y subsistemas de microondas. Sin embargo, a pesar de la gran aceptación de los parámetros s en la comunidad de microondas, el principal inconveniente de esta formulación reside en su limitación para predecir el comportamiento de sistemas no lineales reales. En la actualidad, uno de los principales retos de los diseñadores de microondas es el desarrollo de un contexto análogo que permita integrar tanto el modelado no lineal, como los sistemas de medidas de gran señal y los entornos de simulación no lineal, con el objetivo de extender las capacidades de los parámetros s a regímenes de operación en gran señal y por tanto, obtener una infraestructura que permita tanto la caracterización como el diseño de circuitos no lineales de forma fiable y eficiente. De acuerdo a esta filosofía, en los últimos años se han desarrollado diferentes propuestas como los parámetros X, de Agilent Technologies, o el modelo de Cardiff que tratan de proporcionar esta plataforma común en el ámbito de gran señal. Dentro de este contexto, uno de los objetivos de la presente Tesis es el análisis de la viabilidad del uso de los parámetros X en el diseño y simulación de osciladores para transceptores de microondas. Otro aspecto relevante en el análisis y diseño de circuitos lineales de microondas es la disposición de métodos analíticos sencillos, basados en los parámetros s del transistor, que permitan la obtención directa y rápida de las impedancias de carga y fuente necesarias para cumplir las especificaciones de diseño requeridas en cuanto a ganancia, potencia de salida, eficiencia o adaptación de entrada y salida, así como la determinación analítica de parámetros de diseño clave como el factor de estabilidad o los contornos de ganancia de potencia. Por lo tanto, el desarrollo de una formulación de diseño analítico, basada en los parámetros X y similar a la existente en pequeña señal, permitiría su uso en aplicaciones no lineales y supone un nuevo reto que se va a afrontar en este trabajo. Por tanto, el principal objetivo de la presente Tesis consistiría en la elaboración de una metodología analítica basada en el uso de los parámetros X para el diseño de circuitos no lineales que jugaría un papel similar al que juegan los parámetros s en el diseño de circuitos lineales de microondas. Dichos métodos de diseño analíticos permitirían una mejora significativa en los actuales procedimientos de diseño disponibles en gran señal, así como una reducción considerable en el tiempo de diseño, lo que permitiría la obtención de técnicas mucho más eficientes. Abstract In linear world, classical microwave circuit design relies on the s-parameters due to its capability to successfully characterize the behavior of any linear circuit. Thus the direct use of s-parameters in measurement systems and in linear simulation analysis tools, has facilitated its extensive use and success in the design and characterization of microwave circuits and subsystems. Nevertheless, despite the great success of s-parameters in the microwave community, the main drawback of this formulation is its limitation in the behavior prediction of real non-linear systems. Nowadays, the challenge of microwave designers is the development of an analogue framework that allows to integrate non-linear modeling, large-signal measurement hardware and non-linear simulation environment in order to extend s-parameters capabilities to non-linear regimen and thus, provide the infrastructure for non-linear design and test in a reliable and efficient way. Recently, different attempts with the aim to provide this common platform have been introduced, as the Cardiff approach and the Agilent X-parameters. Hence, this Thesis aims to demonstrate the X-parameter capability to provide this non-linear design and test framework in CAD-based oscillator context. Furthermore, the classical analysis and design of linear microwave transistorbased circuits is based on the development of simple analytical approaches, involving the transistor s-parameters, that are able to quickly provide an analytical solution for the input/output transistor loading conditions as well as analytically determine fundamental parameters as the stability factor, the power gain contours or the input/ output match. Hence, the development of similar analytical design tools that are able to extend s-parameters capabilities in small-signal design to non-linear ap- v plications means a new challenge that is going to be faced in the present work. Therefore, the development of an analytical design framework, based on loadindependent X-parameters, constitutes the core of this Thesis. These analytical nonlinear design approaches would enable to significantly improve current large-signal design processes as well as dramatically decrease the required design time and thus, obtain more efficient approaches.

Análisis, modelado e igualación no lineal de arrays paramétricos acústicos

La presente Tesis analiza y desarrolla metodología específica que permite la caracterización de sistemas de transmisión acústicos basados en el fenómeno del array paramétrico. Este tipo de estructuras es considerado como uno de los sistemas más representativos de la acústica no lineal con amplias posibilidades tecnológicas. Los arrays paramétricos aprovechan la no linealidad del medio aéreo para obtener en recepción señales en el margen sónico a partir de señales ultrasónicas en emisión. Por desgracia, este procedimiento implica que la señal transmitida y la recibida guardan una relación compleja, que incluye una fuerte ecualización así como una distorsión apreciable por el oyente. Este hecho reduce claramente la posibilidad de obtener sistemas acústicos de gran fidelidad. Hasta ahora, los esfuerzos tecnológicos dirigidos al diseño de sistemas comerciales han tratado de paliar esta falta de fidelidad mediante técnicas de preprocesado fuertemente dependientes de los modelos físicos teóricos. Estos están basados en la ecuación de propagación de onda no lineal. En esta Tesis se propone un nuevo enfoque: la obtención de una representación completa del sistema mediante series de Volterra que permita inferir un sistema de compensación computacionalmente ligero y fiable. La dificultad que entraña la correcta extracción de esta representación obliga a desarrollar una metodología completa de identificación adaptada a este tipo de estructuras. Así, a la hora de aplicar métodos de identificación se hace indispensable la determinación de ciertas características iniciales que favorezcan la parametrización del sistema. En esta Tesis se propone una metodología propia que extrae estas condiciones iniciales. Con estos datos, nos encontramos en disposición de plantear un sistema completo de identificación no lineal basado en señales pseudoaleatorias, que aumenta la fiabilidad de la descripción del sistema, posibilitando tanto la inferencia de la estructura basada en bloques subyacente, como el diseño de mecanismos de compensación adecuados. A su vez, en este escenario concreto en el que intervienen procesos de modulación, factores como el punto de trabajo o las características físicas del transductor, hacen inviables los algoritmos de caracterización habituales. Incluyendo el método de identificación propuesto. Con el fin de eliminar esta problemática se propone una serie de nuevos algoritmos de corrección que permiten la aplicación de la caracterización. Las capacidades de estos nuevos algoritmos se pondrán a prueba sobre un prototipo físico, diseñado a tal efecto. Para ello, se propondrán la metodología y los mecanismos de instrumentación necesarios para llevar a cabo el diseño, la identificación del sistema y su posible corrección, todo ello mediante técnicas de procesado digital previas al sistema de transducción. Los algoritmos se evaluarán en términos de error de modelado a partir de la señal de salida del sistema real frente a la salida sintetizada a partir del modelo estimado. Esta estrategia asegura la posibilidad de aplicar técnicas de compensación ya que éstas son sensibles a errores de estima en módulo y fase. La calidad del sistema final se evaluará en términos de fase, coloración y distorsión no lineal mediante un test propuesto a lo largo de este discurso, como paso previo a una futura evaluación subjetiva. ABSTRACT This Thesis presents a specific methodology for the characterization of acoustic transmission systems based on the parametric array phenomenon. These structures are well-known representatives of the nonlinear acoustics field and display large technological opportunities. Parametric arrays exploit the nonlinear behavior of air to obtain sonic signals at the receptors’side, which were generated within the ultrasonic range. The underlying physical process redunds in a complex relationship between the transmitted and received signals. This includes both a strong equalization and an appreciable distortion for a human listener. High fidelity, acoustic equipment based on this phenomenon is therefore difficult to design. Until recently, efforts devoted to this enterprise have focused in fidelity enhancement based on physically-informed, pre-processing schemes. These derive directly from the nonlinear form of the wave equation. However, online limited enhancement has been achieved. In this Thesis we propose a novel approach: the evaluation of a complete representation of the system through its projection onto the Volterra series, which allows the posterior inference of a computationally light and reliable compensation scheme. The main difficulty in the derivation of such representation strives from the need of a complete identification methodology, suitable for this particular type of structures. As an example, whenever identification techniques are involved, we require preliminary estimates on certain parameters that contribute to the correct parameterization of the system. In this Thesis we propose a methodology to derive such initial values from simple measures. Once these information is made available, a complete identification scheme is required for nonlinear systems based on pseudorandom signals. These contribute to the robustness and fidelity of the resulting model, and facilitate both the inference of the underlying structure, which we subdivide into a simple block-oriented construction, and the design of the corresponding compensation structure. In a scenario such as this where frequency modulations occur, one must control exogenous factors such as devices’ operation point and the physical properties of the transducer. These may conflict with the principia behind the standard identification procedures, as it is the case. With this idea in mind, the Thesis includes a series of novel correction algorithms that facilitate the application of the characterization results onto the system compensation. The proposed algorithms are tested on a prototype that was designed and built for this purpose. The methodology and instrumentation required for its design, the identification of the overall acoustic system and its correction are all based on signal processing techniques, focusing on the system front-end, i.e. prior to transduction. Results are evaluated in terms of input-output modelling error, considering a synthetic construction of the system. This criterion ensures that compensation techniques may actually be introduced, since these are highly sensible to estimation errors both on the envelope and the phase of the signals involved. Finally, the quality of the overall system will be evaluated in terms of phase, spectral color and nonlinear distortion; by means of a test protocol specifically devised for this Thesis, as a prior step for a future, subjective quality evaluation.

Contribución al estudio de las metodologías de cálculo realista con incertidumbre (BEPU), dentro del análisis determinista de seguridad de plantas nucleares

El análisis determinista de seguridad (DSA) es el procedimiento que sirve para diseñar sistemas, estructuras y componentes relacionados con la seguridad en las plantas nucleares. El DSA se basa en simulaciones computacionales de una serie de hipotéticos accidentes representativos de la instalación, llamados escenarios base de diseño (DBS). Los organismos reguladores señalan una serie de magnitudes de seguridad que deben calcularse en las simulaciones, y establecen unos criterios reguladores de aceptación (CRA), que son restricciones que deben cumplir los valores de esas magnitudes. Las metodologías para realizar los DSA pueden ser de 2 tipos: conservadoras o realistas. Las metodologías conservadoras utilizan modelos predictivos e hipótesis marcadamente pesimistas, y, por ello, relativamente simples. No necesitan incluir un análisis de incertidumbre de sus resultados. Las metodologías realistas se basan en hipótesis y modelos predictivos realistas, generalmente mecanicistas, y se suplementan con un análisis de incertidumbre de sus principales resultados. Se les denomina también metodologías BEPU (“Best Estimate Plus Uncertainty”). En ellas, la incertidumbre se representa, básicamente, de manera probabilista. Para metodologías conservadores, los CRA son, simplemente, restricciones sobre valores calculados de las magnitudes de seguridad, que deben quedar confinados en una “región de aceptación” de su recorrido. Para metodologías BEPU, el CRA no puede ser tan sencillo, porque las magnitudes de seguridad son ahora variables inciertas. En la tesis se desarrolla la manera de introducción de la incertidumbre en los CRA. Básicamente, se mantiene el confinamiento a la misma región de aceptación, establecida por el regulador. Pero no se exige el cumplimiento estricto sino un alto nivel de certidumbre. En el formalismo adoptado, se entiende por ello un “alto nivel de probabilidad”, y ésta corresponde a la incertidumbre de cálculo de las magnitudes de seguridad. Tal incertidumbre puede considerarse como originada en los inputs al modelo de cálculo, y propagada a través de dicho modelo. Los inputs inciertos incluyen las condiciones iniciales y de frontera al cálculo, y los parámetros empíricos de modelo, que se utilizan para incorporar la incertidumbre debida a la imperfección del modelo. Se exige, por tanto, el cumplimiento del CRA con una probabilidad no menor a un valor P0 cercano a 1 y definido por el regulador (nivel de probabilidad o cobertura). Sin embargo, la de cálculo de la magnitud no es la única incertidumbre existente. Aunque un modelo (sus ecuaciones básicas) se conozca a la perfección, la aplicación input-output que produce se conoce de manera imperfecta (salvo que el modelo sea muy simple). La incertidumbre debida la ignorancia sobre la acción del modelo se denomina epistémica; también se puede decir que es incertidumbre respecto a la propagación. La consecuencia es que la probabilidad de cumplimiento del CRA no se puede conocer a la perfección; es una magnitud incierta. Y así se justifica otro término usado aquí para esta incertidumbre epistémica: metaincertidumbre. Los CRA deben incorporar los dos tipos de incertidumbre: la de cálculo de la magnitud de seguridad (aquí llamada aleatoria) y la de cálculo de la probabilidad (llamada epistémica o metaincertidumbre). Ambas incertidumbres pueden introducirse de dos maneras: separadas o combinadas. En ambos casos, el CRA se convierte en un criterio probabilista. Si se separan incertidumbres, se utiliza una probabilidad de segundo orden; si se combinan, se utiliza una probabilidad única. Si se emplea la probabilidad de segundo orden, es necesario que el regulador imponga un segundo nivel de cumplimiento, referido a la incertidumbre epistémica. Se denomina nivel regulador de confianza, y debe ser un número cercano a 1. Al par formado por los dos niveles reguladores (de probabilidad y de confianza) se le llama nivel regulador de tolerancia. En la Tesis se razona que la mejor manera de construir el CRA BEPU es separando las incertidumbres, por dos motivos. Primero, los expertos defienden el tratamiento por separado de incertidumbre aleatoria y epistémica. Segundo, el CRA separado es (salvo en casos excepcionales) más conservador que el CRA combinado. El CRA BEPU no es otra cosa que una hipótesis sobre una distribución de probabilidad, y su comprobación se realiza de forma estadística. En la tesis, los métodos estadísticos para comprobar el CRA BEPU en 3 categorías, según estén basados en construcción de regiones de tolerancia, en estimaciones de cuantiles o en estimaciones de probabilidades (ya sea de cumplimiento, ya sea de excedencia de límites reguladores). Según denominación propuesta recientemente, las dos primeras categorías corresponden a los métodos Q, y la tercera, a los métodos P. El propósito de la clasificación no es hacer un inventario de los distintos métodos en cada categoría, que son muy numerosos y variados, sino de relacionar las distintas categorías y citar los métodos más utilizados y los mejor considerados desde el punto de vista regulador. Se hace mención especial del método más utilizado hasta el momento: el método no paramétrico de Wilks, junto con su extensión, hecha por Wald, al caso multidimensional. Se decribe su método P homólogo, el intervalo de Clopper-Pearson, típicamente ignorado en el ámbito BEPU. En este contexto, se menciona el problema del coste computacional del análisis de incertidumbre. Los métodos de Wilks, Wald y Clopper-Pearson requieren que la muestra aleatortia utilizada tenga un tamaño mínimo, tanto mayor cuanto mayor el nivel de tolerancia exigido. El tamaño de muestra es un indicador del coste computacional, porque cada elemento muestral es un valor de la magnitud de seguridad, que requiere un cálculo con modelos predictivos. Se hace especial énfasis en el coste computacional cuando la magnitud de seguridad es multidimensional; es decir, cuando el CRA es un criterio múltiple. Se demuestra que, cuando las distintas componentes de la magnitud se obtienen de un mismo cálculo, el carácter multidimensional no introduce ningún coste computacional adicional. Se prueba así la falsedad de una creencia habitual en el ámbito BEPU: que el problema multidimensional sólo es atacable desde la extensión de Wald, que tiene un coste de computación creciente con la dimensión del problema. En el caso (que se da a veces) en que cada componente de la magnitud se calcula independientemente de los demás, la influencia de la dimensión en el coste no se puede evitar. Las primeras metodologías BEPU hacían la propagación de incertidumbres a través de un modelo sustitutivo (metamodelo o emulador) del modelo predictivo o código. El objetivo del metamodelo no es su capacidad predictiva, muy inferior a la del modelo original, sino reemplazar a éste exclusivamente en la propagación de incertidumbres. Para ello, el metamodelo se debe construir con los parámetros de input que más contribuyan a la incertidumbre del resultado, y eso requiere un análisis de importancia o de sensibilidad previo. Por su simplicidad, el modelo sustitutivo apenas supone coste computacional, y puede estudiarse exhaustivamente, por ejemplo mediante muestras aleatorias. En consecuencia, la incertidumbre epistémica o metaincertidumbre desaparece, y el criterio BEPU para metamodelos se convierte en una probabilidad simple. En un resumen rápido, el regulador aceptará con más facilidad los métodos estadísticos que menos hipótesis necesiten; los exactos más que los aproximados; los no paramétricos más que los paramétricos, y los frecuentistas más que los bayesianos. El criterio BEPU se basa en una probabilidad de segundo orden. La probabilidad de que las magnitudes de seguridad estén en la región de aceptación no sólo puede asimilarse a una probabilidad de éxito o un grado de cumplimiento del CRA. También tiene una interpretación métrica: representa una distancia (dentro del recorrido de las magnitudes) desde la magnitud calculada hasta los límites reguladores de aceptación. Esta interpretación da pie a una definición que propone esta tesis: la de margen de seguridad probabilista. Dada una magnitud de seguridad escalar con un límite superior de aceptación, se define el margen de seguridad (MS) entre dos valores A y B de la misma como la probabilidad de que A sea menor que B, obtenida a partir de las incertidumbres de A y B. La definición probabilista de MS tiene varias ventajas: es adimensional, puede combinarse de acuerdo con las leyes de la probabilidad y es fácilmente generalizable a varias dimensiones. Además, no cumple la propiedad simétrica. El término margen de seguridad puede aplicarse a distintas situaciones: distancia de una magnitud calculada a un límite regulador (margen de licencia); distancia del valor real de la magnitud a su valor calculado (margen analítico); distancia desde un límite regulador hasta el valor umbral de daño a una barrera (margen de barrera). Esta idea de representar distancias (en el recorrido de magnitudes de seguridad) mediante probabilidades puede aplicarse al estudio del conservadurismo. El margen analítico puede interpretarse como el grado de conservadurismo (GC) de la metodología de cálculo. Utilizando la probabilidad, se puede cuantificar el conservadurismo de límites de tolerancia de una magnitud, y se pueden establecer indicadores de conservadurismo que sirvan para comparar diferentes métodos de construcción de límites y regiones de tolerancia. Un tópico que nunca se abordado de manera rigurosa es el de la validación de metodologías BEPU. Como cualquier otro instrumento de cálculo, una metodología, antes de poder aplicarse a análisis de licencia, tiene que validarse, mediante la comparación entre sus predicciones y valores reales de las magnitudes de seguridad. Tal comparación sólo puede hacerse en escenarios de accidente para los que existan valores medidos de las magnitudes de seguridad, y eso ocurre, básicamente en instalaciones experimentales. El objetivo último del establecimiento de los CRA consiste en verificar que se cumplen para los valores reales de las magnitudes de seguridad, y no sólo para sus valores calculados. En la tesis se demuestra que una condición suficiente para este objetivo último es la conjunción del cumplimiento de 2 criterios: el CRA BEPU de licencia y un criterio análogo, pero aplicado a validación. Y el criterio de validación debe demostrarse en escenarios experimentales y extrapolarse a plantas nucleares. El criterio de licencia exige un valor mínimo (P0) del margen probabilista de licencia; el criterio de validación exige un valor mínimo del margen analítico (el GC). Esos niveles mínimos son básicamente complementarios; cuanto mayor uno, menor el otro. La práctica reguladora actual impone un valor alto al margen de licencia, y eso supone que el GC exigido es pequeño. Adoptar valores menores para P0 supone menor exigencia sobre el cumplimiento del CRA, y, en cambio, más exigencia sobre el GC de la metodología. Y es importante destacar que cuanto mayor sea el valor mínimo del margen (de licencia o analítico) mayor es el coste computacional para demostrarlo. Así que los esfuerzos computacionales también son complementarios: si uno de los niveles es alto (lo que aumenta la exigencia en el cumplimiento del criterio) aumenta el coste computacional. Si se adopta un valor medio de P0, el GC exigido también es medio, con lo que la metodología no tiene que ser muy conservadora, y el coste computacional total (licencia más validación) puede optimizarse. ABSTRACT Deterministic Safety Analysis (DSA) is the procedure used in the design of safety-related systems, structures and components of nuclear power plants (NPPs). DSA is based on computational simulations of a set of hypothetical accidents of the plant, named Design Basis Scenarios (DBS). Nuclear regulatory authorities require the calculation of a set of safety magnitudes, and define the regulatory acceptance criteria (RAC) that must be fulfilled by them. Methodologies for performing DSA van be categorized as conservative or realistic. Conservative methodologies make use of pessimistic model and assumptions, and are relatively simple. They do not need an uncertainty analysis of their results. Realistic methodologies are based on realistic (usually mechanistic) predictive models and assumptions, and need to be supplemented with uncertainty analyses of their results. They are also termed BEPU (“Best Estimate Plus Uncertainty”) methodologies, and are typically based on a probabilistic representation of the uncertainty. For conservative methodologies, the RAC are simply the restriction of calculated values of safety magnitudes to “acceptance regions” defined on their range. For BEPU methodologies, the RAC cannot be so simple, because the safety magnitudes are now uncertain. In the present Thesis, the inclusion of uncertainty in RAC is studied. Basically, the restriction to the acceptance region must be fulfilled “with a high certainty level”. Specifically, a high probability of fulfillment is required. The calculation uncertainty of the magnitudes is considered as propagated from inputs through the predictive model. Uncertain inputs include model empirical parameters, which store the uncertainty due to the model imperfection. The fulfillment of the RAC is required with a probability not less than a value P0 close to 1 and defined by the regulator (probability or coverage level). Calculation uncertainty is not the only one involved. Even if a model (i.e. the basic equations) is perfectly known, the input-output mapping produced by the model is imperfectly known (unless the model is very simple). This ignorance is called epistemic uncertainty, and it is associated to the process of propagation). In fact, it is propagated to the probability of fulfilling the RAC. Another term used on the Thesis for this epistemic uncertainty is metauncertainty. The RAC must include the two types of uncertainty: one for the calculation of the magnitude (aleatory uncertainty); the other one, for the calculation of the probability (epistemic uncertainty). The two uncertainties can be taken into account in a separate fashion, or can be combined. In any case the RAC becomes a probabilistic criterion. If uncertainties are separated, a second-order probability is used; of both are combined, a single probability is used. On the first case, the regulator must define a level of fulfillment for the epistemic uncertainty, termed regulatory confidence level, as a value close to 1. The pair of regulatory levels (probability and confidence) is termed the regulatory tolerance level. The Thesis concludes that the adequate way of setting the BEPU RAC is by separating the uncertainties. There are two reasons to do so: experts recommend the separation of aleatory and epistemic uncertainty; and the separated RAC is in general more conservative than the joint RAC. The BEPU RAC is a hypothesis on a probability distribution, and must be statistically tested. The Thesis classifies the statistical methods to verify the RAC fulfillment in 3 categories: methods based on tolerance regions, in quantile estimators and on probability (of success or failure) estimators. The former two have been termed Q-methods, whereas those in the third category are termed P-methods. The purpose of our categorization is not to make an exhaustive survey of the very numerous existing methods. Rather, the goal is to relate the three categories and examine the most used methods from a regulatory standpoint. Special mention deserves the most used method, due to Wilks, and its extension to multidimensional variables (due to Wald). The counterpart P-method of Wilks’ is Clopper-Pearson interval, typically ignored in the BEPU realm. The problem of the computational cost of an uncertainty analysis is tackled. Wilks’, Wald’s and Clopper-Pearson methods require a minimum sample size, which is a growing function of the tolerance level. The sample size is an indicator of the computational cost, because each element of the sample must be calculated with the predictive models (codes). When the RAC is a multiple criteria, the safety magnitude becomes multidimensional. When all its components are output of the same calculation, the multidimensional character does not introduce additional computational cost. In this way, an extended idea in the BEPU realm, stating that the multi-D problem can only be tackled with the Wald extension, is proven to be false. When the components of the magnitude are independently calculated, the influence of the problem dimension on the cost cannot be avoided. The former BEPU methodologies performed the uncertainty propagation through a surrogate model of the code, also termed emulator or metamodel. The goal of a metamodel is not the predictive capability, clearly worse to the original code, but the capacity to propagate uncertainties with a lower computational cost. The emulator must contain the input parameters contributing the most to the output uncertainty, and this requires a previous importance analysis. The surrogate model is practically inexpensive to run, so that it can be exhaustively analyzed through Monte Carlo. Therefore, the epistemic uncertainty due to sampling will be reduced to almost zero, and the BEPU RAC for metamodels includes a simple probability. The regulatory authority will tend to accept the use of statistical methods which need a minimum of assumptions: exact, nonparametric and frequentist methods rather than approximate, parametric and bayesian methods, respectively. The BEPU RAC is based on a second-order probability. The probability of the safety magnitudes being inside the acceptance region is a success probability and can be interpreted as a fulfillment degree if the RAC. Furthermore, it has a metric interpretation, as a distance (in the range of magnitudes) from calculated values of the magnitudes to acceptance regulatory limits. A probabilistic definition of safety margin (SM) is proposed in the thesis. The same from a value A to other value B of a safety magnitude is defined as the probability that A is less severe than B, obtained from the uncertainties if A and B. The probabilistic definition of SM has several advantages: it is nondimensional, ranges in the interval (0,1) and can be easily generalized to multiple dimensions. Furthermore, probabilistic SM are combined according to the probability laws. And a basic property: probabilistic SM are not symmetric. There are several types of SM: distance from a calculated value to a regulatory limit (licensing margin); or from the real value to the calculated value of a magnitude (analytical margin); or from the regulatory limit to the damage threshold (barrier margin). These representations of distances (in the magnitudes’ range) as probabilities can be applied to the quantification of conservativeness. Analytical margins can be interpreted as the degree of conservativeness (DG) of the computational methodology. Conservativeness indicators are established in the Thesis, useful in the comparison of different methods of constructing tolerance limits and regions. There is a topic which has not been rigorously tackled to the date: the validation of BEPU methodologies. Before being applied in licensing, methodologies must be validated, on the basis of comparisons of their predictions ad real values of the safety magnitudes. Real data are obtained, basically, in experimental facilities. The ultimate goal of establishing RAC is to verify that real values (aside from calculated values) fulfill them. In the Thesis it is proved that a sufficient condition for this goal is the conjunction of 2 criteria: the BEPU RAC and an analogous criterion for validation. And this las criterion must be proved in experimental scenarios and extrapolated to NPPs. The licensing RAC requires a minimum value (P0) of the probabilistic licensing margin; the validation criterion requires a minimum value of the analytical margin (i.e., of the DG). These minimum values are basically complementary; the higher one of them, the lower the other one. The regulatory practice sets a high value on the licensing margin, so that the required DG is low. The possible adoption of lower values for P0 would imply weaker exigence on the RCA fulfillment and, on the other hand, higher exigence on the conservativeness of the methodology. It is important to highlight that a higher minimum value of the licensing or analytical margin requires a higher computational cost. Therefore, the computational efforts are also complementary. If medium levels are adopted, the required DG is also medium, and the methodology does not need to be very conservative. The total computational effort (licensing plus validation) could be optimized.

Valorización energética de biomasas en el marco de la política energética española. Incentivos económico-financieros y políticos, aportación de valor añadido y prospectiva estratégica de desarrollo

El modelo económico actual basado en el consumo y en la búsqueda permanente de una mayor calidad de vida, unido a una población mundial en aumento, contribuye a incrementar la demanda de servicios energéticos para cubrir las necesidades de energía de las personas y las industrias. Desde finales del siglo XIX la energía se ha generado fundamentalmente a partir de combustibles fósiles (carbón, petróleo y gas), convertidos en el suministro energético predominante mundialmente. Las emisiones de gases de efecto invernadero que genera la prestación de servicios energéticos han contribuido considerablemente al aumento histórico de las concentraciones de esos gases en la atmósfera, hasta el punto de que el consumo de combustibles fósiles es responsable de la mayoría de las emisiones antropogénicas (IPCC, 2012). Existen diversas opciones para disminuir las emisiones de gases de efecto invernadero del sector energético y con ello contribuir a mitigar el cambio climático, entre otras sería viable aumentar la eficiencia energética y sustituir combustibles de origen fósil por combustibles de origen renovable, pudiendo garantizar un suministro de energía sostenible, competitivo y seguro. De todas las energías renovables susceptibles de formar parte de una cartera de opciones de mitigación, esta tesis se centra en la bioenergía generada a partir de la valorización energética de las biomasas agrícolas, forestales, ganaderas o de otro tipo, con fines eléctricos y térmicos. Con objeto de mostrar su capacidad para contribuir a mitigar el cambio climático y su potencial contribución al desarrollo socioeconómico, a la generación de energía distribuida y a reducir determinados efectos negativos sobre el medio ambiente, se ha analizado minuciosamente el sector español de la biomasa en su conjunto. Desde el recurso biomásico que existe en España, las formas de extraerlo y procesarlo, las tecnologías de valorización energética, sus usos energéticos principales y la capacidad de implementación del sector en España. Asimismo se ha examinado el contexto energético tanto internacional y europeo como nacional, y se han analizado pormenorizadamente los instrumentos de soporte que han contribuido de manera directa e indirecta al desarrollo del sector en España. Además, la tesis integra el análisis de los resultados obtenidos mediante dos metodologías diferentes con fines también distintos. Por un lado se han obtenido los resultados medioambientales y socioeconómicos de los análisis de ciclo de vida input-output generados a partir de las cinco tecnologías biomásicas más ampliamente utilizadas en España. Y por otro lado, en base a los objetivos energéticos y medioambientales establecidos, se han obtenido distintas proyecciones de la implementación del sector a medio plazo, en forma de escenarios energéticos con horizonte 2035, mediante el modelo TIMES-Spain. La tesis ofrece también una serie de conclusiones y recomendaciones que podrían resultar pertinentes para los agentes que constituyen la cadena de valor del propio sector e interesados, así como para la formulación de políticas y mecanismos de apoyo para los agentes decisores, tanto del ámbito de la Administración General del Estado como autonómico y regional, sobre las características y ventajas de determinadas formas de valorización, sobre los efectos sociales y medioambientales que induce su uso, y sobre la capacidad de sector para contribuir a determinadas políticas más allá de las puramente energéticas. En todo caso, esta tesis doctoral aspira a contribuir a la toma de decisiones idóneas tanto a los agentes del sector como a responsables públicos, con objeto de adoptar medidas orientadas a fomentar modificaciones del sistema energético que incrementen la proporción de energía renovable y, de esta forma, contribuir a mitigar la amenaza que supone el cambio climático no solo en la actualidad, sino especialmente en los próximos años para las generaciones venideras. ABSTRACT The current economic model based on both, consumption and the constant search for greater quality of life, coupled with a growing world population, contribute to increase the demand for energy services in order to meet the energy needs of people and industries. Since the late nineteenth century, energy has been basically generated from fossil fuels (coal, oil and gas), which converted fossil fuels into the predominant World energy supply source. Emissions of greenhouse gases generated by the provision of energy services have contributed significantly to the historical increase in the concentrations of these gases in the atmosphere, to the extent that the consumption of fossil fuels is responsible for most of the anthropogenic emissions (IPCC, 2012). There are several options to reduce emissions of greenhouse gases in the energy sector and, thereby, to contribute to mitigate climate change. Among others, would be feasible to increase energy efficiency and progressively replacing fossil fuels by renewable fuels, which are able to ensure a sustainable, competitive and secure energy supply. Of all the renewable energies likely to form part of a portfolio of mitigation options, this thesis focuses on bioenergy generated from agricultural, forestry, farming or other kind of biomass, with electrical and thermal purposes. In order to show their ability to contribute to mitigate climate change and its potential contribution to socio-economic development, distributed energy generation and to reduce certain negative effects on the environment, the Spanish biomass sector as a whole has been dissected. From the types of biomass resources that exist in Spain, ways of extracting and processing them, energy production technologies, its main energy uses and the implementation capacity of the sector in Spain. It has also examined the international, European and national energy context, and has thoroughly analyzed the support instruments that have contributed directly and indirectly to the development of the sector in Spain, so far. Furthermore, the thesis integrates the analysis of results obtained using two different methodologies also with different purposes. On the one hand, the environmental and socio-economic results of the analysis of input-output cycle life generated from the five biomass technologies most widely used in Spain, have been obtained. On the other hand, different projections of the implementation of the sector in the medium term, as energy scenarios with horizon 2035, have been obtained by the model TIMES-Spain, based on several energy and environmental objectives. The thesis also offers a series of conclusions and recommendations that could be relevant to the agents that constitute the value chain of the biomass sector itself and other stakeholders. As well as policy and support mechanisms for decision-makers, from both: The Central and Regional Governments, on the characteristics and advantages of certain forms of valorization, on the social and environmental effects that induce their use, and the ability of the biomass sector to contribute to certain policies beyond the purely energy ones. In any case, this thesis aims to contribute to decision making, suitable for both: Industry players and to public officials. In order to adopt measures to promote significant changes in the energy system that increase the proportion of renewable energy and, consequently, that contribute to mitigate the threat of climate change; not only today but in the coming years, especially for future generations.

Identificação de sistemas para o estudo de controle motor.

Qualquer tarefa motora ativa se dá pela ativação de uma população de unidades motoras. Porém, devido a diversas dificuldades, tanto técnicas quanto éticas, não é possível medir a entrada sináptica dos motoneurônios em humanos. Por essas razões, o uso de modelos computacionais realistas de um núcleo de motoneurônios e as suas respectivas fibras musculares tem um importante papel no estudo do controle humano dos músculos. Entretanto, tais modelos são complexos e uma análise matemática é difícil. Neste texto é apresentada uma abordagem baseada em identificação de sistemas de um modelo realista de um núcleo de unidades motoras, com o objetivo de obter um modelo mais simples capaz de representar a transdução das entradas do núcleo de unidades motoras na força do músculo associado ao núcleo. A identificação de sistemas foi baseada em um algoritmo de mínimos quadrados ortogonal para achar um modelo NARMAX, sendo que a entrada considerada foi a condutância sináptica excitatória dendrítica total dos motoneurônios e a saída foi a força dos músculos produzida pelo núcleo de unidades motoras. O modelo identificado reproduziu o comportamento médio da saída do modelo computacional realista, mesmo para pares de sinal de entrada-saída não usados durante o processo de identificação do modelo, como sinais de força muscular modulados senoidalmente. Funções de resposta em frequência generalizada do núcleo de motoneurônios foram obtidas do modelo NARMAX, e levaram a que se inferisse que oscilações corticais na banda-beta (20 Hz) podem influenciar no controle da geração de força pela medula espinhal, comportamento do núcleo de motoneurônios até então desconhecido.

Aplicação de modelos de defeitos na geração de conjuntos de teste completos a partir de Sistemas de Transição com Entrada/Saída

O Teste Baseado em Modelos (TBM) emergiu como uma estratégia promissora para minimizar problemas relacionados à falta de tempo e recursos em teste de software e visa verificar se a implementação sob teste está em conformidade com sua especificação. Casos de teste são gerados automaticamente a partir de modelos comportamentais produzidos durante o ciclo de desenvolvimento de software. Entre as técnicas de modelagem existentes, Sistemas de Transição com Entrada/Saída (do inglês, Input/Output Transition Systems - IOTSs), são modelos amplamente utilizados no TBM por serem mais expressivos do que Máquinas de Estado Finito (MEFs). Apesar dos métodos existentes para geração de testes a partir de IOTSs, o problema da seleção de casos de testes é um tópico difícil e importante. Os métodos existentes para IOTS são não-determinísticos, ao contrário da teoria existente para MEFs, que fornece garantia de cobertura completa com base em um modelo de defeitos. Esta tese investiga a aplicação de modelos de defeitos em métodos determinísticos de geração de testes a partir de IOTSs. Foi proposto um método para geração de conjuntos de teste com base no método W para MEFs. O método gera conjuntos de teste de forma determinística além de satisfazer condições de suficiência de cobertura da especificação e de todos os defeitos do domínio de defeitos definido. Estudos empíricos avaliaram a aplicabilidade e eficácia do método proposto: resultados experimentais para analisar o custo de geração de conjuntos de teste utilizando IOTSs gerados aleatoriamente e um estudo de caso com especificações da indústria mostram a efetividade dos conjuntos gerados em relação ao método tradicional de Tretmans.