La planification forestière est un exercice ardu qui demande des efforts humains et matériels très importants. Plusieurs modèles de planification forestière sont proposés dans la littérature. Cependant, ils sont généralement difficiles à implémenter. En outre, leur performance est limitée aux forêts de petite taille ou de taille moyenne. L’objectif général de cette thèse est de proposer des outils simples et efficients qui permettent de réduire la complexité de cet exercice en utilisant la modélisation mathématique. Les modèles proposés dans cette thèse traitent le problème de planification forestière sur un horizon tactique et considèrent principalement deux restrictions spatiales qui sont des restrictions d’adjacence ou de période de régénération et des restrictions de rendement soutenu de récolte. Dans la première partie, nous avons proposé une nouvelle formulation mathématique basée sur la programmation en nombres entiers qui permet de modéliser le problème de planification tactique sous des contraintes spatiales. Nous avons eu recours à une méthode de résolution exacte à l’aide d’un solveur commercial pour résoudre le problème. Nous avons utilisé plusieurs forêts réelles de tailles différentes dont trois sont de très grandes forêts dans la région de la Mauricie au Québec. Pour évaluer la performance de notre modèle proposé, nous avons utilisé deux autres formulations mathématiques qui sont déjà proposées dans la littérature pour résoudre le même problème sous les mêmes contraintes. Les résultats ont montré que notre formulation est plus simple d’un point de vue de la modélisation et de l’implémentation. En ce qui concerne la résolution, les résultats ont montré qu’il n’y pas une formulation qui avait une performance supérieure lorsque la taille des forêts est petite ou moyenne. Cependant, lorsque la taille des forêts devient large, notre formulation avait une performance supérieure. Dans la deuxième partie de la thèse, nous avons proposé deux méthodes pour résoudre le même problème en utilisant les grandes forêts. La première méthode est une heuristique basée sur la relaxation linéaire. Pour évaluer sa performance, nous avons comparé les résultats qu’elle a fournis avec ceux obtenus en utilisant un solveur commercial pendant une durée de résolution prédéfinie. Les résultats ont montré que le solveur est plus performant lorsque le problème est mono-période. Cependant, lorsque le problème est multi-période, contrairement à notre heuristique, le solveur n’était pas capable d’obtenir un résultat pendant le temps de résolution fixé. La deuxième méthode consiste à contrôler la taille des problèmes à résoudre à travers un indicateur qui contrôle la forme des regroupements de peuplements à couper. En plus de sa simplicité, cette méthode permet de réduire significativement la taille des problèmes et d’avoir une planification opérationnelle plus efficiente par l’élimination des regroupements de peuplement ayants des formes non désirables. Dans la troisième partie, nous avons proposé un modèle intégré stochastique bi-niveau qui permet la planification de la coupe tout en considérant les feux de forêts. La formulation mathématique proposée dans la première partie de la thèse représente le premier niveau. Le deuxième niveau est un modèle de propagation de feux. Il considère des points d’ignition du feu, des points critiques tels que des municipalités et des croisements de routes forestières et deux scénarios de différentes directions et vitesses de vent. L’intégration est possible à travers la reformulation primal-dual des conditions d’optimalités du problème du plus court chemin entre les points d’ignition et les points critiques. Nous avons utilisé des forêts hypothétiques pour les expérimentations. Les résultats ont montré que les regroupements de peuplements planifiés sont localisés dans des zones qui permettent de retarder la propagation du feu des points d’ignition vers les points critiques. / Forest planning is a difficult task that needs important human and material efforts. Many forest planning modelsare proposed in the literature. However, their implementation is usually complex. More over, their performanceis limited to small and medium forests. The whole objective of this thesis is to propose simple and efficient tools to reduce the complexity of the forest planning using mathematical modeling. In this thesis, proposed models are applied for tactical forest planning and consider mainly two spatial restrictions that are adjacency or green-up restrictions and sustainable timber yields restrictions.In the first part, we proposed a new mathematical formulation based on integer programming to model the tactical forest planning under spatial restrictions. We solved the model with an exact method using a commercial solver. We used many real forests with different sizes. Three of them are very large forests and are located in the Mauricie region of Quebec. We evaluated the performance of the proposed model using two other formulation sexisting in the literature to model the same problem under the same restrictions. Results showed that there isno formulation that out performed others when small and medium forests are used. However, the proposed formulation out performed others when large forests are used.In the second part, we proposed two heuristic approaches to solve the same problem using large forests. The first one is a linear relaxation heuristic. We compared results obtained with the heuristic with those obtained with a commercial solver when a computational time is predefined. Results showed that the solver is more performant when mono-periods problems were solved. However, when multi-periods problems were solved, unlike ourheuristic, the solver does not find any feasible solution within the fixed computational time. The second heuristicis a size reduction heuristic which uses a shape index to control the shape of the clusters of stands. Besides its simplicity, this method is very beneficial as it reduces significantly problems size and anticipate the operational planning by eliminating non-profitable clusters. In the third part, we proposed a stochastic programming model formulated as an integrated mixed integerprogramming model for scheduling harvesting operations, which directly takes into account a fire spread model. It is a bi-level model. In the first level, we used the tactical forest planning model proposed in the first part. The second level represents the fire spread model that considers fire spread between ignition points and critical points, which could be, for example, municipalities and important road crossings. It considers multiple scenarios, which takes into account wind direction, and wind spread. The integration is possible due to a reformulation of a bi-level optimization problem where the lower level shortest path representation is written explicitly with aprimal-dual reformulation of its optimality conditions. Results based on a set of illustrative instances are presented. Results vary depending on spatial restrictions considered in the model. However, harvesting scheduled are usually allocated to delay fires reaching critical points and tends to allocate harvesting operations around critical points.

Ce mémoire décrit plusieurs modèles qui cherchent à prédire le vieillissement du cerveau chez des individus âgés ayant une cognition intacte. Nous effectuons d’abord des tests statistiques afin de comparer différents modèles. À partir de ceux-ci, nous présentons des estimés du moment critique à partir duquel certaines structures cérébrales ont décru morphométriquement de façon notable en guise d’indice de fragilité. Nous mesurons également l’asymétrie hémisphérique de la décroissance de chaque région. Nous proposons différentes façons d’ordonnancer les régions du cerveau selon différentes quantités mesurant le taux vieillissement. Finalement, nous établissons un outil qui permet d’évaluer l’âge virtuel correspondant aux mesures cérébrales pour le cerveau d’un individu, selon les estimés de normalité. / In this memoir, we propose different models in order to predict the effect of aging on brain atrophy, computed from neuroimagery data. We begin with different statistical tests in order to chose between models and assert if different regions have significant asymmetry in their development. We then compute critical moments after which we predict that the measures of a region have significantly decreased. We use these critical moments to order the regions, as we suppose it’s a good indicator of fragility in aging. We finish by proposing a method of predicting the neurological age of an individual from it’s neuroimagery data.

Titre de l'écran-titre (visionné le 26 avril 2023) / Ce mémoire vise à étudier le calcul stochastique et certains schémas et applications numériques dans des problèmes réels. Comme nous le savons, au cours des dernières décennies, il y a eu des demandes importantes d'outils et de stratégies de calcul stochastique dans divers domaines. L'un des domaines qui en a le plus besoin est la finance mathématique, où de nombreuses quantités financières, telles que le prix des actifs, sont des variables aléatoires ; Ainsi, les mathématiques algébriques, qui traitent des variables déterministes et non aléatoires, en finance, ont été moins employées à la place dans les équations différentielles stochastiques en tant que branche des mathématiques stochastiques ont reçu une attention significative dans le domaine de la finance. Tout au long de ce travail, en recueillant quelques faits et informations auxiliaires dans la théorie des probabilités, nous étudions le mouvement brownien, qui peut être observé dans des systèmes perturbés dans divers domaines, dont l'électronique, la physique, la biologie, la gestion et l'économie. Comme l'un des principaux objectifs de ce projet , nous introduisons l'intégration par rapport au mouvement brownien et aux processus associés (processus Itô), ce qui nous permettra de définir les équations différentielles stochastiques (EDSs). / This dissertation aims to study stochastic calculus and some numerical schemes and applications in real-life problems. As we know, during the last few decades, there have been significant demands for implements and strategies of stochastic calculus in various fields. One of the fields that has the most need is mathematical finance, where many financial quantities, such as the price of as sets, are random variables; thus, algebraic mathematics, which deals with deterministic and non-random variables, in finance has been employed less and instead in stochastic differential equations as a branch of stochastic mathematics have received significant attention in the field of finance. Throughout this work, by collecting some auxiliary facts and information in the probability theory, we study the Brownian motion, which can be observed in perturbed systems in various fields, including electronics, physics, biology, management, and economics. As one of the main goals of this project, we introduce the integration with respect to the Brownian motion and the related processes (Itô processes), which will enable us to define the stochastic differential equations (SDEs).

Titre de l'écran-titre (visionné le 8 août 2023) / Les anodes de carbone sont des composants essentiels du procédé Hall-Héroult utilisé dans l'industrie de l'aluminium primaire. La mise en forme des anodes par compactage ou vibrocompactage a un impact significatif sur leur qualité et leurs différentes propriétés. La présence de gradients de densité dans les anodes peut causer leur surconsommation dans les cuves d'électrolyse et réduire leur performance. Ce projet vise à étudier le comportement thermomécanique de la pâte d'anode lors de sa mise en forme, à modéliser ce comportement et à examiner l'effet de la géométrie du moule sur la distribution de la densité dans l'anode. L'effet de la température de la pâte d'anode sur son comportement mécanique lors de la compaction a été étudié expérimentalement. Dans ce contexte, des essais de compaction monotones et des essais de friction ont été menés à des températures situées entre 130 °C et 170 °C. Les résultats ont révélé les effets importants de la température sur le comportement non linéaire de la pâte et ont également permis de déterminer les coefficients de friction statique et dynamique de l'interface pâte/acier. Une loi de comportement thermoviscoplastique non linéaire a été développée afin de modéliser le comportement de la pâte d'anode lors de sa compaction. Cette loi est basée sur un cadre thermodynamique, la théorie des grandes déformations et le concept de configuration naturelle. Les paramètres de la loi évoluent durant la compaction en fonction de la densité et de la température. Ils ont été identifiés à l'aide d'une procédure d'identification inverse. La loi de comportement a été implémentée dans le logiciel de simulation par éléments finis Abaqus à travers une sous-routine VUMAT, destinée à l'analyse dynamique explicite. Les simulations par éléments finis portant sur la mise en forme par compaction de géométries complexes ont été en mesure de prédire avec une bonne précision les profils de densité mesurés par tomodensitométrie. L'effet de la température ainsi que de différentes caractéristiques géométriques, telles que des arêtes chanfreinées et arrondies, des tourillons et des rainures de différentes dimensions, ont été étudié à l'aide de simulations par éléments finis. Dans cette optique, de petites anodes ont été compactées par éléments finis afin de déterminer si les paramètres susmentionnés peuvent améliorer l'uniformité de la densité. Par la suite, des simulations par éléments finis ont été réalisées sur une anode industrielle et différentes modifications de la géométrie du moule ont été proposées. Les résultats de ces simulations suggèrent que la réduction de la taille des rainures et la modification de leur géométrie pourraient améliorer de façon significative l'uniformité de la densité des anodes industrielles. / Carbon anodes are essential components of the Hall-Héroult process used in the primary aluminium industry. The anode forming process by compaction or vibrocompaction has a significant impact on their quality and properties. The presence of density gradients in the anode can cause carbon overconsumption in the electrolysis cells and negatively affect their performance. This project aims to study the thermomechanical behavior of the anode paste during compaction, to model this behavior and to examine the effect of the mold geometry on the density distribution in the anode. The effect of the anode paste temperature on its mechanical behavior during compaction was studied experimentally. In this context, monotonic compaction tests and friction tests were conducted at temperatures between 130 °C and 170 °C. The results revealed the significant effect of temperature on the nonlinear behavior of the paste and also allowed the determination of the static and kinetic friction coefficients of the paste/steel interface. A nonlinear thermo-viscoplastic constitutive law was developed to model the behavior of the anode paste during compaction and is based on a thermodynamic framework, the finite strain theory and the concept of natural configuration. The constitutive law's material parameters evolve during compaction as a function of density and temperature and were identified using an inverse identification procedure. The constitutive law was implemented in the finite element simulation software Abaqus through a VUMAT subroutine, intended for explicit dynamic analysis. The finite element simulations of the compaction of complex geometries were able to predict with good accuracy the density profiles measured by X-ray computed tomography. The effect of temperature as well as of different geometrical features, such as chamfered and rounded edges, stub holes and slots of different dimensions, was investigated using finite element analysis. To this end, small anodes were compacted through finite element simulations to determine if the aforementioned parameters can improve the density uniformity. Additional finite element simulations were performed on an industrial anode and different modifications to the mold geometry were proposed. The results of these simulations show that reducing the size of the slots and modifying their geometry could improve the density uniformity of industrial anodes.

Thèse en cotutelle, doctorat en génie civil : Université Laval, Québec, Canada, Philosophiæ doctor (Ph. D.) et Université Toulouse III - Paul Sabatier, Toulouse, France / Un modèle de comportement mécanique des bétons fibrés est développé et implanté dans le logiciel éléments finis Cast3m. Il vient compléter le modèle de fissuration orthotrope du béton Fluendo3D, en ajoutant la capacité de traiter des matériaux à fibres courtes, cylindriques et rectilignes. Ce modèle, qui s'appuie sur des données d'essais et certains mécanismes modélisés issus de la littérature, permet d'apporter des éléments de compréhension du comportement de ce type de matériau, et notamment concernant le phénomène de multi-fissuration. La caractérisation des phénomènes mis enjeu durant l'extraction des fibres est le point de départ de cette étude. Les effets de l'inclinaison des fibres par rapport à la direction d'extraction sont pris en compte et interviennent dans le comportement du modèle qui présente la capacité d'utiliser des orientations préférentielles de fibres. Le phénomène de multi-fissuration est représenté grâce à une loi d'effet d'échelle de Weibull qui permet de tenir compte de la dispersion des résistances à la traction du béton et d'expliquer le développement de la multi-fissuration. Cette représentation permet d'obtenir des distributions d'ouvertures de fissures dans des macro-éléments finis, et apporte donc une nouvelle précision dans le calcul des ouvertures de fissures dans des structures de grandes dimensions. / A model for the mechanical behavior of fiber-reinforced concrete is developed and implemented in Cast3m finite element software. It completes the orthotropic concrete cracking model Fluendo3D, by adding the ability to treat materials with short, cylindrical and straight fibers. This model, which is based on tests results and assumptions of mechanisms modelled from litterature, provides elements of understanding of the behavior of this type of material, especially concerning the phenomenon of multi-cracking. The characterization of the phenomena involved during the extraction of the fibers is the starting point of this study. The effects of the inclination of the fibers with respect to the direction of extraction are taken into account and intervene in the behavior of the model which presents the capacity to use preferential orientations of fibers. The multi-cracking phenomenon is represented by a Weibull scaling law that allows to take into account the dispersion of the concrete tensile strengths and to explain the development of multi-cracking process. This representation makes it possible to obtain crack opening distributions in macro-finite elements, and thus brings a new precision in the calculation of cracks openings in large structures.

De nos jours, la conception parasismique des structures en acier repose sur la dissipation d'énergie sismique par l'action de grandes déformations plastiques. Les cadres à diagonales ductiles confinées sont un type de système à contreventements concentriques caractérisé par des membrures qui plastifient autant en traction qu'en compression. Grâce à un mécanisme de restreinte latérale qui renferme un noyau métallique, les diagonales ductiles confinées (DDC) tirent avantage de la ductilité cyclique du matériau. Cependant, les DDC affichent fréquemment une faible rigidité après plastification, causant ainsi d'importantes déformations latérales durant un séisme. De plus, la plastification du noyau est généralement associée à un seul objectif de performance; son comportement sous d'autres intensités sismiques est donc essentiellement inconnu. Les DDC hybrides sont étudiées en tant qu'alternative aux systèmes traditionnels pour pallier leurs limitations. Le concept d'hybridité vise à maîtriser différentes caractéristiques de métaux minutieusement sélectionnés pour obtenir une réponse souhaitable. Cette étude numérique comporte trois objectifs. Les matériaux de noyaux potentiels sont d'abord évalués afin de déterminer les meilleures combinaisons de deux métaux avec des propriétés mécaniques complémentaires. L'analyse de données expérimentales indique que le comportement plastique des DDC hybrides est amélioré par l'emploi d'un acier au carbone 350WT conjointement avec un autre métal qui possède une faible limite élastique et une grande capacité d'écrouissage (p. ex. l'acier inoxydable 304L, l'alliage d'aluminium 5083-O, l'acier au carbone A36, ou l'acier LYP-100). Par la suite, deux nouvelles configurations de DDC hybrides sont conçues afin d'accueillir trois noyaux ductiles connectés en parallèle. La première option proposée comprend un mécanisme de restreinte fait de tubes en acier remplis de béton, tandis que la deuxième option est fabriquée exclusivement à partir de composantes métalliques. Finalement, plusieurs simulations par la méthode des éléments finis sont réalisées sur des modèles numériques pour valider quantitativement l'accroissement de la performance. Comparées aux DDC conventionnelles, les DDC hybrides présentent une réponse d'écrouissage améliorée, une légère augmentation de la rigidité axiale, ainsi qu'une plus grande capacité à dissiper l'énergie. / Contemporary seismic-resistant design of steel structures relies on the dissipation of earthquake energy through significant inelastic deformations. Buckling-restrained braced frames (BRBFs) are a type of concentrically-braced system characterized by braces that yield both in tension and in compression. Thanks to a restraining mechanism that confines a ductile steel core, buckling-restrained braces (BRBs) can take advantage of the cyclic ductility of the steel material. However, BRBs commonly display a low post-yield stiffness, causing substantial interstory drifts and large residual drifts after seismic events. Moreover, yielding of the core is often tied to only a single performance objective, thus making its response at other levels of seismicity largely unpredictable. Hybrid BRB solutions are explored as an alternative to the traditional BRB system to overcome its limitations. The hybrid concept is hinged on harnessing different characteristics from different materials that are carefully combined into one ductile design to achieve a desirable response. This numerical study has three main objectives. Potential core metals are first evaluated to determine the best combinations of two materials with complementary engineering properties. Analysis of experimental data indicates that the post-yield behavior of hybrid BRBs is improved by employing 350WT carbon steel in conjunction with another metal, which possesses a low-yield and high-strain-hardening capacity (e.g., 304L stainless steel, 5083-O aluminum alloy, A36 carbon steel, or LYP-100 low-yield-point steel). Afterwards, two new hybrid BRB systems are designed to accommodate the complex deformation pattern of three core plates connected in parallel. The first proposed option has a restraining mechanism made from concrete-filled steel tubes, while the second hybrid BRB option is fabricated exclusively from metal plate components. Lastly, multiple finite element simulations are carried out on numerical models to quantitatively validate the performance enhancement. Compared to conventional BRBs, hybrid BRBs exhibit an improved strain hardening response, a slight increase in axial stiffness, and a greater energy dissipation capability for an equivalent brace strength.

Titre de l'écran-titre (visionné le 5 septembre 2023) / Au vu de l'urgence climatique actuelle et des différents rapports scientifiques au sujet du dérèglement climatique, il est primordial et même vital de réduire drastiquement la pollution générée par l'Homme. Le dernier rapport du GIEC (Groupe d'experts Intergouvernemental sur l'Évolution du Climat) (2022) indique qu'il faut effectivement réduire de moitié les émissions d'ici 2030 et insiste fortement sur la nécessité d'agir immédiatement afin de préserver la planète. Dans ce sens, l'industrie de la production de béton est responsable de 4 à 8% des rejets totaux humains de dioxyde de carbone et doit donc urgemment évoluer afin de réduire son impact environnemental. L'objectif principal de cette étude est de participer à cette transition industrielle en développant un modèle numérique fiable et exploitable afin d'optimiser la production, diminuer le temps de malaxage et réduire le gaspillage du béton à l'aide des outils technologiques de contrôle qualité déjà mis en place. En effet, développer une simulation numérique permettant de mieux comprendre le comportement et les profils d'écoulement du béton frais à l'intérieur d'une toupie de malaxage est extrêmement prometteur puisqu'elle permet alors d'encore mieux optimiser les temps et les coûts de malaxage. Afin d'espérer exploiter un tel outil numérique complexe, la mise en place de modèles élémentaires d'écoulement du béton frais est essentielle afin de valider, caractériser et calibrer les simulations numériques. Dans ce mémoire, le développement de trois modèles simples d'écoulement est traité et les résultats obtenus ont permis de valider le comportement numérique de l'écoulement du béton frais. Chacun de ces modèles présente des forces et faiblesses et participe à la création d'un environnement de travail numérique qui permet de beaucoup mieux appréhender la rhéologie et le comportement d'écoulement du béton frais. Ce projet de recherche est donc une véritable porte d'entrée vers la modélisation complète de la production du béton frais. / In view of the current climate emergency and the various scientific reports on climate change, it is essential and even vital to drastically reduce man-made pollution. The latest IPCC (Intergovernmental Panel on Climate Change) report (2022) indicates that emissions must be halved by 2030 and strongly emphasizes the need to act immediately to preserve the planet. In this sense, the concrete production industry is responsible for 4-8% of total human carbon dioxide emissions and therefore urgently needs to evolve to reduce its environmental impact. The main objective of this study is to participate in this industrial transition by developing a reliable and exploitable numerical model to optimize the production, reduce mixing time and also reduce concrete waste by using technological quality control tools already available. Indeed, developing a numerical simulation allowing to better understand the behavior and flow profiles of fresh concrete inside a mixing-truck is extremely promising as it allows for further optimization of mixing times and costs. In order to be able to exploit such a complex numerical tool, the implementation of elementary fresh concrete flow models is essential to validate, characterize and calibrate the numerical simulations. In this thesis, the development of three simple flow models is discussed and the results obtained are used to validate the numerical behavior of fresh concrete flow. Each of these models has strengths and weaknesses and contributes to the creation of a numerical working environment that provides a much better understanding of the rheology and flow behavior of fresh concrete. This research project is therefore a real gateway to a full modelling of fresh concrete production.

Ce mémoire présente les travaux de mise en œuvre d'un modèle numérique par éléments finis de la cellule unitaire représentative d'un composite interlock 3D, à partir de tomographies. Ce matériau se compose de torons tissés dans les directions chaine et trame, ainsi qu'à travers l'épaisseur. Ce type d'interlock 3X a été largement étudié ces dix dernières années, sans pour autant que des modèles numériques performants soient réalisés, ce qui rend nécessaire une modélisation numérique de sa cellule unitaire. Il présente pourtant des avantages vis-à-vis des autres types de composites, notamment grâce à l'absence du phénomène de délamination, cause principale de dégradation des composites. Ce travail présente la démarche et les moyens mis en œuvre pour obtenir la cellule unitaire numérique du matériau. Des tomographies de différents états du tissu sec, non compacté et compacté, ont été réalisées afin de générer un modèle éléments finis à partir du logiciel de traitement d'image AVIZO ainsi que du logiciel de génération de textile TEXGEN. Un modèle simplifié d'une tranche de cellule unitaire a été choisi afin de simuler et valider le procédé de compaction numérique avec le logiciel élément finis Abaqus. Les résultats du modèle de compaction ont été comparés aux tomographies du tissu sec compacté. Les résultats de cette simulation simplifiée sont prometteurs. La déformation et les déplacements des torons obtenus sur la cellule unitaire du tissu sec simulé sont très proches de ceux observés sur les tomographies du tissu compacté de la sorte. / This thesis presents the work conducted toward the development of a finite element model of the representative unit cell of a 3D interlock composite, based on tomographies. This material is composed of yarns woven in the warp and weft directions, as well as through the thickness. This kind of 3X interlock has been very little studied to date, which makes numerical modelling of its unit cell necessary. However, it present advantages over other types of composites, in particular the absence of delamination phenomenon, the main cause of composite degradation. This work presents the approach and the means employed to obtain the numerical unit cell of the material. Tomographies of different states of the dry, uncompacted and compacted fabric were carried out in order to generate a finite element model from the image processing software AVIZO and the textile generation software TEXGEN. A simplified model of a unit cell slice was chosen in order to simulate and validate the numerical compaction process with the Abaqus finite element software. The results of the compaction model were compared with tomographies of dry texile. The results of this simplified simulation are promising. The deformation and displacements of the yarns obtained on the unit cell of the simulated dry textile are very close to those observed on tomographies of the compacted textile.

Il a été noté que les approches de planification forestière au niveau stratégique ont traditionnellement ignoré l'interaction avec le reste de la chaîne de valeur [57]. En effet, bien que les décisions prises à ce niveau soient à l'heure actuelle supportées par des modèles d'approvisionnement forestier, elles portent peu d'attention à la valeur et aux coûts des produits forestiers, aux capacités de transformation et à la localisation des usines et des marchés [27]. Afin de faciliter l'intégration de ces éléments dans la planification forestière, ce mémoire propose un formalisme de représentation et de description des unités d'affaires, des produits, des flux de matière, des processus et des marchés. Le modèle proposé dans ce mémoire a comme objectif de supporter la planification stratégique dans l'industrie forestière.

La prévision hydrologique consiste à évaluer quelle sera l'évolution du débit au cours des prochains pas de temps. En utilisant les systèmes actuels de prévisions hydrologiques déterministes, il est impossible d'apprécier simplement l'incertitude associée à ce type de prévision, ce que peut nuire à la prise de décisions. La prévision hydrologique d'ensemble (PHE) cherche à étayer cette incertitude en proposant, à chaque pas de temps, une distribution de probabilité, la prévision probabiliste, en place et lieu d'une estimation unique du débit, la prévision déterministe. La PHE offre de nombreux bénéfices : elle informe l'utilisateur de l'incertitude; elle permet aux autorités qui prennent des décisions de déterminer des critères d'alerte et de mettre en place des scénarios d'urgence; elle fournit les informations nécessaires à la prise de décisions tenant compte du risque. L'objectif principal de cette thèse est l'évaluation de prévisions hydrologiques d'ensemble, en mettant l'accent sur la performance et la fiabilité de celles-ci. Deux techniques pour construire des ensembles sont explorées: a) une première reposant sur des prévisions météorologiques d'ensemble (PME) et b) une seconde exploitant simultanément un ensemble de modèles hydrologiques (multimodèle). En termes généraux, les objectifs de la thèse ont été établis afin d'évaluer : a) les incertitudes associées à la structure du modèle : une étude qui repose sur des simulations journalières issues de dix-sept modèles hydrologiques globaux, pour plus de mille bassins versants français; b) les incertitudes associées à la prévision météorologique : une étude qui exploite la PME du Service Météorologique du Canada et un modèle hydrologique opérationnel semi-distribué, pour un horizon de 3 jours sur douze bassins versants québécois; c) les incertitudes associées à la fois à la structure du modèle et à la prévision météorologique : une étude qui repose à la fois sur la PME issue du ECMWF (European Centre for Medium-Range Weather Forecasts) et seize modèles hydrologiques globaux, pour un horizon de 9 jours sur 29 bassins versants français. Les résultats mets en évidence les avantages des systèmes probabilistes par rapport aux les déterministes. Les prévisions probabilistes sont toutefois souvent affectées par une sous dispersion de leur distribution prédictive. Elles exigent alors un post traitement avant d'être intégrées dans un processus de prise de décision. Plus intéressant encore, les résultats ont également montré le grand potentiel de combiner plusieurs sources d'incertitude, notamment celle associée à la prévision météorologique et celle associée à la structure des modèles hydrologiques. Il nous semble donc prioritaire de continuer à explorer davantage cette approche combinatoire.