Short-term Hydropower Planning System

The Thesis gives a decision support framework that has significant impact on the economic performance and viability of a hydropower company. The studyaddresses the short-term hydropower planning problem in the Nordic deregulated electricity market. The basics of the Nordic electricity market, trading mechanisms, hydropower system characteristics and production planning are presented in the Thesis. The related modelling theory and optimization methods are covered aswell. The Thesis provides a mixed integer linear programming model applied in asuccessive linearization method for optimal bidding and scheduling decisions inthe hydropower system operation within short-term horizon. A scenario based deterministic approach is exploited for modelling uncertainty in market price and inflow. The Thesis proposes a calibration framework to examine the physical accuracy and economic optimality of the decisions suggested by the model. A calibration example is provided with data from a real hydropower system using a commercial modelling application with the mixed integer linear programming solver CPLEX.

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions

Du principe de continuité de l'exploitation en comptabilité : contribution à la formalisation de la démarche cognitive du réviseur aux comptes dans la validation de ce principe à l'aide d'un système abductif de modélisation de croyance basé sur un processus stochastique

Résumé : J'ai souvent vu des experts être d'avis contraires. Je n'en ai jamais vu aucun avoir tort. Auguste Detoeuf Propos d'O.L. Brenton, confiseur, Editions du Tambourinaire, 1948. En choisissant volontairement une problématique comptable typiquement empirique, ce travail s'est attelé à tenter de démontrer la possibilité de produire des enseignements purement comptables (ie à l'intérieur du schème de représentation de la Comptabilité) en s'interdisant l'emprunt systématique de theories clé-en-main à l'Économie -sauf quant cela s'avère réellement nécessaire et légitime, comme dans l'utilisation du CAPM au chapitre précédent. Encore une fois, rappelons que cette thèse n'est pas un réquisitoire contre l'approche économique en tant que telle, mais un plaidoyer visant à mitiger une telle approche en Comptabilité. En relation avec le positionnement épistémologique effectué au premier chapitre, il a été cherché à mettre en valeur l'apport et la place de la Comptabilité dans l'Économie par le positionnement de la Comptabilité en tant que discipline pourvoyeuse de mesures de représentation de l'activité économique. Il nous paraît clair que si l'activité économique, en tant que sémiosphère comptable directe, dicte les observations comptables, la mesure de ces dernières doit, tant que faire se peut, tenter de s'affranchir de toute dépendance à la discipline économique et aux théories-méthodes qui lui sont liées, en adoptant un mode opératoire orthogonal, rationnel et systématique dans le cadre d'axiomes lui appartenant en propre. Cette prise de position entraîne la définition d'un nouveau cadre épistémologique par rapport à l'approche positive de la Comptabilité. Cette dernière peut se décrire comme l'expression philosophique de l'investissement de la recherche comptable par une réflexion méthodique propre à la recherche économique. Afin d'être au moins partiellement validé, ce nouveau cadre -que nous voyons dérivé du constructivisme -devrait faire montre de sa capacité à traiter de manière satisfaisante une problématique classique de comptabilité empirico-positive. Cette problématique spécifique a été choisie sous la forme de traitement-validation du principe de continuité de l'exploitation. Le principe de continuité de l'exploitation postule (énonciation d'une hypothèse) et établit (vérification de l'hypothèse) que l'entreprise produit ses états financiers dans la perspective d'une poursuite normale de ses activités. Il y a rupture du principe de continuité de l'exploitation (qui devra alors être écartée au profit du principe de liquidation ou de cession) dans le cas de cessation d'activité, totale ou partielle, volontaire ou involontaire, ou la constatation de faits de nature à compromettre la continuité de l'exploitation. Ces faits concernent la situation financière, économique et sociale de l'entreprise et représentent l'ensemble des événements objectifs 33, survenus ou pouvant survenir, susceptibles d'affecter la poursuite de l'activité dans un avenir prévisible. A l'instar de tous les principes comptables, le principe de continuité de l'exploitation procède d'une considération purement théorique. Sa vérification requiert toutefois une analyse concrète, portant réellement et de manière mesurable à conséquence, raison pour laquelle il représente un thème de recherche fort apprécié en comptabilité positive, tant il peut (faussement) se confondre avec les études relatives à la banqueroute et la faillite des entreprises. Dans la pratique, certaines de ces études, basées sur des analyses multivariées discriminantes (VIDA), sont devenues pour l'auditeur de véritables outils de travail de par leur simplicité d'utilisation et d'interprétation. À travers la problématique de ce travail de thèse, il a été tenté de s'acquitter de nombreux objectifs pouvant être regroupés en deux ensembles : celui des objectifs liés à la démarche méthodologique et celui relevant de la mesure-calibration. Ces deux groupes-objectifs ont permis dans une dernière étape la construction d'un modèle se voulant une conséquence logique des choix et hypothèses retenus.

External validation of the ASTRAL score to predict 3- and 12-month functional outcome in the China National Stroke Registry.

BACKGROUND AND PURPOSE: The ASTRAL score was recently introduced as a prognostic tool for acute ischemic stroke. It predicts 3-month outcome reliably in both the derivation and the validation European cohorts. We aimed to validate the ASTRAL score in a Chinese stroke population and moreover to explore its prognostic value to predict 12-month outcome. METHODS: We applied the ASTRAL score to acute ischemic stroke patients admitted to 132 study sites of the China National Stroke Registry. Unfavorable outcome was assessed as a modified Rankin Scale score >2 at 3 and 12 months. Areas under the curve were calculated to quantify the prognostic value. Calibration was assessed by comparing predicted and observed probability of unfavorable outcome using Pearson correlation coefficient. RESULTS: Among 3755 patients, 1473 (39.7%) had 3-month unfavorable outcome. Areas under the curve for 3 and 12 months were 0.82 and 0.81, respectively. There was high correlation between observed and expected probability of unfavorable 3- and 12-month outcome (Pearson correlation coefficient: 0.964 and 0.963, respectively). CONCLUSIONS: ASTRAL score is a reliable tool to predict unfavorable outcome at 3 and 12 months after acute ischemic stroke in the Chinese population. It is a useful tool that can be readily applied in clinical practice to risk-stratify acute stroke patients.

Predicting Non Return to Work after Orthopaedic Trauma: The Wallis Occupational Rehabilitation RisK (WORRK) Model.

BACKGROUND: Workers with persistent disabilities after orthopaedic trauma may need occupational rehabilitation. Despite various risk profiles for non-return-to-work (non-RTW), there is no available predictive model. Moreover, injured workers may have various origins (immigrant workers), which may either affect their return to work or their eligibility for research purposes. The aim of this study was to develop and validate a predictive model that estimates the likelihood of non-RTW after occupational rehabilitation using predictors which do not rely on the worker's background. METHODS: Prospective cohort study (3177 participants, native (51%) and immigrant workers (49%)) with two samples: a) Development sample with patients from 2004 to 2007 with Full and Reduced Models, b) External validation of the Reduced Model with patients from 2008 to March 2010. We collected patients' data and biopsychosocial complexity with an observer rated interview (INTERMED). Non-RTW was assessed two years after discharge from the rehabilitation. Discrimination was assessed by the area under the receiver operating curve (AUC) and calibration was evaluated with a calibration plot. The model was reduced with random forests. RESULTS: At 2 years, the non-RTW status was known for 2462 patients (77.5% of the total sample). The prevalence of non-RTW was 50%. The full model (36 items) and the reduced model (19 items) had acceptable discrimination performance (AUC 0.75, 95% CI 0.72 to 0.78 and 0.74, 95% CI 0.71 to 0.76, respectively) and good calibration. For the validation model, the discrimination performance was acceptable (AUC 0.73; 95% CI 0.70 to 0.77) and calibration was also adequate. CONCLUSIONS: Non-RTW may be predicted with a simple model constructed with variables independent of the patient's education and language fluency. This model is useful for all kinds of trauma in order to adjust for case mix and it is applicable to vulnerable populations like immigrant workers.

High-performance liquid chromatography of the renal blood flow marker p-aminohippuric acid (PAH) and its metabolite N-acetyl PAH improves PAH clearance measurements.

PAH (N-(4-aminobenzoyl)glycin) clearance measurements have been used for 50 years in clinical research for the determination of renal plasma flow. The quantitation of PAH in plasma or urine is generally performed by colorimetric method after diazotation reaction but the measurements must be corrected for the unspecific residual response observed in blank plasma. We have developed a HPLC method to specifically determine PAH and its metabolite NAc-PAH using a gradient elution ion-pair reversed-phase chromatography with UV detection at 273 and 265 nm, respectively. The separations were performed at room temperature on a ChromCart (125 mmx4 mm I.D.) Nucleosil 100-5 microm C18AB cartridge column, using a gradient elution of MeOH-buffer pH 3.9 1:99-->15:85 over 15 min. The pH 3.9 buffered aqueous solution consisted in a mixture of 375 ml sodium citrate-citric acid solution (21.01 g citric acid and 8.0 g NaOH per liter), added up with 2.7 ml H3PO4 85%, 1.0 g of sodium heptanesulfonate and completed ad 1000 ml with ultrapure water. The N-acetyltransferase activity does not seem to notably affect PAH clearances, although NAc-PAH represents 10.2+/-2.7% of PAH excreted unchanged in 12 healthy subjects. The performance of the HPLC and the colorimetric method have been compared using urine and plasma samples collected from healthy volunteers. Good correlations (r=0.94 and 0.97, for plasma and urine, respectively) are found between the results obtained with both techniques. However, the colorimetric method gives higher concentrations of PAH in urine and lower concentrations in plasma than those determined by HPLC. Hence, both renal (ClR) and systemic (Cls) clearances are systematically higher (35.1 and 17.8%, respectively) with the colorimetric method. The fraction of PAH excreted by the kidney ClR/ClS calculated from HPLC data (n=143) is, as expected, always <1 (mean=0.73+/-0.11), whereas the colorimetric method gives a mean extraction ratio of 0.87+/-0.13 implying some unphysiological values (>1). In conclusion, HPLC not only enables the simultaneous quantitation of PAH and NAc-PAH, but may also provide more accurate and precise PAH clearance measurements.

Project for Active and Robot Vision Course

Tässä työssä raportoidaan harjoitustyön kehittäminen ja toteuttaminen Aktiivisen- ja robottinäön kurssille. Harjoitustyössä suunnitellaan ja toteutetaan järjestelmä joka liikuttaa kappaleita robottikäsivarrella kolmiuloitteisessa avaruudessa. Kappaleidenpaikkojen määrittämiseen järjestelmä käyttää digitaalisia kuvia. Tässä työssä esiteltävässä harjoitustyötoteutuksessa käytettiin raja-arvoistusta HSV-väriavaruudessa kappaleiden segmentointiin kuvasta niiden värien perusteella. Segmentoinnin tuloksena saatavaa binäärikuvaa suodatettiin mediaanisuotimella kuvan häiriöiden poistamiseksi. Kappaleen paikkabinäärikuvassa määritettiin nimeämällä yhtenäisiä pikseliryhmiä yhtenäisen alueen nimeämismenetelmällä. Kappaleen paikaksi määritettiin suurimman nimetyn pikseliryhmän paikka. Kappaleiden paikat kuvassa yhdistettiin kolmiuloitteisiin koordinaatteihin kalibroidun kameran avulla. Järjestelmä liikutti kappaleita niiden arvioitujen kolmiuloitteisten paikkojen perusteella.

On Effort Estimation in Software Projects

Software engineering is criticized as not being engineering or 'well-developed' science at all. Software engineers seem not to know exactly how long their projects will last, what they will cost, and will the software work properly after release. Measurements have to be taken in software projects to improve this situation. It is of limited use to only collect metrics afterwards. The values of the relevant metrics have to be predicted, too. The predictions (i.e. estimates) form the basis for proper project management. One of the most painful problems in software projects is effort estimation. It has a clear and central effect on other project attributes like cost and schedule, and to product attributes like size and quality. Effort estimation can be used for several purposes. In this thesis only the effort estimation in software projects for project management purposes is discussed. There is a short introduction to the measurement issues, and some metrics relevantin estimation context are presented. Effort estimation methods are covered quite broadly. The main new contribution in this thesis is the new estimation model that has been created. It takes use of the basic concepts of Function Point Analysis, but avoids the problems and pitfalls found in the method. It is relativelyeasy to use and learn. Effort estimation accuracy has significantly improved after taking this model into use. A major innovation related to the new estimationmodel is the identified need for hierarchical software size measurement. The author of this thesis has developed a three level solution for the estimation model. All currently used size metrics are static in nature, but this new proposed metric is dynamic. It takes use of the increased understanding of the nature of the work as specification and design work proceeds. It thus 'grows up' along with software projects. The effort estimation model development is not possible without gathering and analyzing history data. However, there are many problems with data in software engineering. A major roadblock is the amount and quality of data available. This thesis shows some useful techniques that have been successful in gathering and analyzing the data needed. An estimation process is needed to ensure that methods are used in a proper way, estimates are stored, reported and analyzed properly, and they are used for project management activities. A higher mechanism called measurement framework is also introduced shortly. The purpose of the framework is to define and maintain a measurement or estimationprocess. Without a proper framework, the estimation capability of an organization declines. It requires effort even to maintain an achieved level of estimationaccuracy. Estimation results in several successive releases are analyzed. It isclearly seen that the new estimation model works and the estimation improvementactions have been successful. The calibration of the hierarchical model is a critical activity. An example is shown to shed more light on the calibration and the model itself. There are also remarks about the sensitivity of the model. Finally, an example of usage is shown.

Stratification of the severity of critically ill patients with classification trees

Background: Development of three classification trees (CT) based on the CART (Classification and Regression Trees), CHAID (Chi-Square Automatic Interaction Detection) and C4.5 methodologies for the calculation of probability of hospital mortality; the comparison of the results with the APACHE II, SAPS II and MPM II-24 scores, and with a model based on multiple logistic regression (LR). Methods: Retrospective study of 2864 patients. Random partition (70:30) into a Development Set (DS) n = 1808 and Validation Set (VS) n = 808. Their properties of discrimination are compared with the ROC curve (AUC CI 95%), Percent of correct classification (PCC CI 95%); and the calibration with the Calibration Curve and the Standardized Mortality Ratio (SMR CI 95%). Results: CTs are produced with a different selection of variables and decision rules: CART (5 variables and 8 decision rules), CHAID (7 variables and 15 rules) and C4.5 (6 variables and 10 rules). The common variables were: inotropic therapy, Glasgow, age, (A-a)O2 gradient and antecedent of chronic illness. In VS: all the models achieved acceptable discrimination with AUC above 0.7. CT: CART (0.75(0.71-0.81)), CHAID (0.76(0.72-0.79)) and C4.5 (0.76(0.73-0.80)). PCC: CART (72(69- 75)), CHAID (72(69-75)) and C4.5 (76(73-79)). Calibration (SMR) better in the CT: CART (1.04(0.95-1.31)), CHAID (1.06(0.97-1.15) and C4.5 (1.08(0.98-1.16)). Conclusion: With different methodologies of CTs, trees are generated with different selection of variables and decision rules. The CTs are easy to interpret, and they stratify the risk of hospital mortality. The CTs should be taken into account for the classification of the prognosis of critically ill patients.

Prediction of crude protein and classification of the growth stage of wheat plant samples from NIR spectra

Near-infrared spectroscopy (NIRS) was used to analyse the crude protein content of dried and milled samples of wheat and to discriminate samples according to their stage of growth. A calibration set of 72 samples from three growth stages of wheat (tillering, heading and harvest) and a validation set of 28 samples was collected for this purpose. Principal components analysis (PCA) of the calibration set discriminated groups of samples according to the growth stage of the wheat. Based on these differences, a classification procedure (SIMCA) showed a very accurate classification of the validation set samples : all of them were successfully classified in each group using this procedure when both the residual and the leverage were used in the classification criteria. Looking only at the residuals all the samples were also correctly classified except one of tillering stage that was assigned to both tillering and heading stages. Finally, the determination of the crude protein content of these samples was considered in two ways: building up a global model for all the growth stages, and building up local models for each stage, separately. The best prediction results for crude protein were obtained using a global model for samples in the two first growth stages (tillering and heading), and using a local model for the harvest stage samples.

A priori parameterisation of the CERES soil-crop models and tests against several European data sets

Mechanistic soil-crop models have become indispensable tools to investigate the effect of management practices on the productivity or environmental impacts of arable crops. Ideally these models may claim to be universally applicable because they simulate the major processes governing the fate of inputs such as fertiliser nitrogen or pesticides. However, because they deal with complex systems and uncertain phenomena, site-specific calibration is usually a prerequisite to ensure their predictions are realistic. This statement implies that some experimental knowledge on the system to be simulated should be available prior to any modelling attempt, and raises a tremendous limitation to practical applications of models. Because the demand for more general simulation results is high, modellers have nevertheless taken the bold step of extrapolating a model tested within a limited sample of real conditions to a much larger domain. While methodological questions are often disregarded in this extrapolation process, they are specifically addressed in this paper, and in particular the issue of models a priori parameterisation. We thus implemented and tested a standard procedure to parameterize the soil components of a modified version of the CERES models. The procedure converts routinely-available soil properties into functional characteristics by means of pedo-transfer functions. The resulting predictions of soil water and nitrogen dynamics, as well as crop biomass, nitrogen content and leaf area index were compared to observations from trials conducted in five locations across Europe (southern Italy, northern Spain, northern France and northern Germany). In three cases, the model’s performance was judged acceptable when compared to experimental errors on the measurements, based on a test of the model’s root mean squared error (RMSE). Significant deviations between observations and model outputs were however noted in all sites, and could be ascribed to various model routines. In decreasing importance, these were: water balance, the turnover of soil organic matter, and crop N uptake. A better match to field observations could therefore be achieved by visually adjusting related parameters, such as field-capacity water content or the size of soil microbial biomass. As a result, model predictions fell within the measurement errors in all sites for most variables, and the model’s RMSE was within the range of published values for similar tests. We conclude that the proposed a priori method yields acceptable simulations with only a 50% probability, a figure which may be greatly increased through a posteriori calibration. Modellers should thus exercise caution when extrapolating their models to a large sample of pedo-climatic conditions for which they have only limited information.

Use of the heat dissipation method for sap flow measurement in citrus nursery trees1

Sap flow could be used as physiological parameter to assist irrigation of screen house citrus nursery trees by continuous water consumption estimation. Herein we report a first set of results indicating the potential use of the heat dissipation method for sap flow measurement in containerized citrus nursery trees. 'Valencia' sweet orange [Citrus sinensis (L.) Osbeck] budded on 'Rangpur' lime (Citrus limonia Osbeck) was evaluated for 30 days during summer. Heat dissipation probes and thermocouple sensors were constructed with low-cost and easily available materials in order to improve accessibility of the method. Sap flow showed high correlation to air temperature inside the screen house. However, errors due to natural thermal gradient and plant tissue injuries affected measurement precision. Transpiration estimated by sap flow measurement was four times higher than gravimetric measurement. Improved micro-probes, adequate method calibration, and non-toxic insulating materials should be further investigated.

Kuumalanka-anemometri virtausmittauksissa

Tässä työssä selostetaan kuumalanka-anemometrin käyttö virtausmittauksissa. Kuumalanka-anemometrilla saadaan mitattua virtausnopeuden ja -suunnan lisäksi nopeusheilahteluja. Mittaustaajuus on tyypillisesti useita kymmeniä tuhansia mittauksia sekunnissa ja signaali on jatkuva. Nykytekniikalla pystytään helposti tallentamaan mittauslaitteistolta saatu viesti tietokoneelle ja muuntamaan se nopeudeksi. Hetkellisten nopeuksien avulla voidaan laskea turbulenttisen virtauksen ominaisuuksia, kuten turbulenssin intensiteetti ja spektri. Kuumalanka-anemometrissa lämmitetään sähköisesti ohutta lankaa, joka on mitattavassa virtauksessa. Langan sähköteho on suunnilleen yhtäsuuri kuin langasta konvektiolla siirtyvä lämpöteho. Tällöin on teoreettisesti mahdollista laskea virtausnopeus lämpötehosta lämmönsiirtokorrelaatioilla. Käytännössä laitteisto joudutaan kuitenkin erikseen kalibroimaan, mutta sähkötehon teoreettista riippuvuutta konvektiosta käytetään hyväksi. Kuumalangan lämmitettävä osuus on tyypillisesti halkaisijaltaan 5 µm ja pituudeltaan noin 1 mm. Sitä käytetään pääasiassa kaasuvirtausten mittaamiseen ja valtaosassa mittauksissa virtausaineena on ilma. Kuumalanka voi olla toteutettu kuumakalvotekniikalla, jossa halkaisijaltaan noin 50 - 70 µm paksuinen kuitu on päällystetty ohuella sähköä johtavalla kalvolla. Kuumakalvoanturin ei tarvitse olla muodoltaan sylinterimäinen, se voi olla mm. kartiomainen tai kiilamainen. Erikoispäällystetyllä kuumakalvoanturilla on mahdollista mitata myös nestevirtauksia. Mitattaessa kaasuvirtauksia kuumakalvon etuna on selvästi parempi kestävyys verrattuna kuumalankaan. Nimitystä kuumalanka-anemometri käytetään yleisesti molemmista anturityypeistä Tämän työn alussa käsitellään sylinterin yli tapahtuvaan virtaukseen liittyvää virtausmekaniikkaa ja lämmönsiirtoa. Anemometrin sähköinen osa, laitteisto ja sen kalibrointi käydään läpi. Langan suuntariippuvuuden laskentaan esitetään tarvittavat yhtälöt. Työssä esitellään kolme laitteistolla tehtyä perusmittausta: anturin kohtauskulman muuttaminen, pyörähdyssymmerisen suihkun nopeuskenttä ja tuulitunnelin rajakerros. Lisäksi esitellään yksi käytännöllinen ja vaativampi mittaus, jossa on mitattu nopeusprofiili radiaalikompressorin diffuusorin loppuosassa.

Unification of Emission Measurements in Paper and Pulp Mills

Tässä diplomityössä tavoitteena on selvittää käytännöt ja menetelmät, joilla UPM-Kymmenen sellu-, paperi- ja vaneritehtaat ympäri maailman mittaavat ja laskevat ilmapäästönsä. Tämä tehdään kyselylomakkeella, joka lähetetään tehtaiden ympäristöpäälliköille. Kaikki tärkeimmät seikat ilmapäästöihin liittyen, kuten vaaditut jatkuvatoimiset mittaukset, jaksottaiset mittaukset, raportointikäytäntö, kalibrointi jne. kysytään lomakkeessa. Kyselylomakkeessa painotetaan mittauskäytäntöä sellutehtaissa sekä energiantuotannossa. Saatujen tulosten perusteella annetaan ehdotuksia sekä ohjeita tulevaisuutta varten, jotta mittaustulosten kokoaminen helpottuisi ja vertailukelpoisuus paranisi. Työn kirjallisuusosuudessa selvitetään yleisimmät päästölähteet sekä päästökomponentit paperi –ja selluteollisuudessa. Näiden ei toivottujen yhdisteiden syntymekanismit sekä menetelmät niiden poistamiseksi savukaasuista on myös lyhyesti kuvailtu. Myös erilaiset analysointi- ja näytteenottomenetelmät on kerrottu. Erot tehtaiden ympäristöluvissa käydään läpi, jakaen tehtaat kolmeen maantieteelliseen ryhmään. Lupakäytäntöjen osalta Suomen osuutta on painotettu, sillä UPM-Kymmene on varsin Suomikeskeinen yhtiö tehtaiden lukumääriin ja sijainteihin katsottuna. Viranomaismääräykset sekä päästörajat muutamista tehtaista on esitetty havainnollistaakseen alueellisia eroja.

Tuotantolaitteiden kvalifiointi ja prosessin validointi GMP-tuotantoa varten

GMP-säädösten mukaan aktiivisten lääkeaineiden, kriittisten lääkeaineintermediaattien ja lääkeapuaineiden valmistusprosessit pitää validoida. Validointityöhön kuuluu oleellisesti tuotantolaitteiden kvalifiointi ja prosessin validointi. Käytännössä tuotantolaitteiden kvalifiointi toteutetaan tekemällä laitteille suunnitelmien tarkastus (DQ), asennus- ja käyttöönottotarkastus (IQ), toiminnan testaus (OQ) sekä suorituskykytestit (PQ). Tuotantolaitteiden kvalifiointiin kuuluu myös laitteiden asianmukaisten kalibrointi-, kunnossapito- ja puhdistusohjeiden sekä työ- ja toimintaohjeiden (SOP:ien) laatiminen. Prosessin validoinnissa laaditaan dokumentoidut todisteet siitä, että prosessi toimii vakaasti ja tuotteelle asetetut vaatimukset täyttyvät johdonmukaisesti. GMP-tuotantolaitteiden kvalifiointiin ja lääkevalmistusprosessin validointiin on laadittu erilaisia GMP-säädöksiä noudattavia yleisiä validointiohjeita, kuten PIC/S:n ja FDA:n ohjeet kvalifioinnista ja validoinnista. IVT/SC on laatinut yksiselitteiset validointistandardit validointityön selventämiseksi. Validoinnin tilastolliseen tarkasteluun on käytettävissä GHTF:n laatimat tilastolliset validointimenetelmät. Yleensä tuotantolaitteiden kvalifiointi ja prosessin validointi tehdään ennen lääkevalmisteen kaupallisen tuotannon aloittamista. Kvalifiointi- ja validointityö voidaan tehdä kuitenkin myös tuotannon yhteydessä (konkurrentisti) tai retrospektiivisesti käyttäen hyväksi valmistettujen tuotantoerien prosessitietoja. Tässä työssä laadittiin Kemira Fine Chemicals Oy:n Kokkolan GMP-tuotantolinjan lääkeaineintermediaattiprosessin validoinnin yleissuunnitelma (VMP), joka sisältää sekä tuotantolaitteiden kvalifiointisuunnitelman että prosessin validointisuunnitelman. Suunnitelmissa huomioitiin tuotantolaitteiden aikaisempi käyttö muuhun hienokemikaalituotantoon ja tuotantolinjan muuttaminen GMP-vaatimusten mukaiseksi. Työhön kuului myös tuotantolaitteiden kvalifiointityön tekeminen laaditun suunnitelman mukaisesti.