Dans le cadre de recherches menées par l'ANDRA, qui visent à évaluer la faisabilité de créer un stockage de déchets radioactifs dans la formation argileuse du Callovo Oxfordien, cette thèse s'intéresse à deux étapes indispensables pour obtenir une image 3D du sous sol. La première étape concerne la cartographie de la géométrie de la formation. Celle-ci est déterminée à partir de temps pointés le long de profils sismiques, orientés suivant deux directions principales. Aux intersections, les valeurs enregistrées par deux profils sont différentes. Le faible nombre d'intersections ne permettant pas une étude détaillée du comportement spatial de ces différences, elles sont analysées par le biais du pseudo variogramme croisé. L'analyse variographique permet de choisir un modèle cohérent pour l'ensemble des profils. À partir de ce modèle une estimation du temps corrigée est présentée. La deuxième étape concerne l'étude de l'impédance, variable liée aux propriétés intrinsèques de la roche. Deux types de mesure sont disponibles. D'une part l'impédance sismique résultant de l'inversion géophysique de l'amplitude enregistrée le long de profils de sismique réflexion est largement étendue dans la zone d'intérêt. D'autre part, l'impédance log, résultant du produit de la densité de la roche et de la vitesse des ondes P, n'est connue uniquement dans quelques puits. Le comportement spatial bivariable, le long de la direction verticale, a permis de formuler l'hypothèse suivante : chaque mesure d'impédance est décomposée en une somme d'une variable Z (commune aux deux mesures) et d'un résidu différent pour chaque mesure. Sous certaines hypothèses supplémentaires, il est possible d'estimer la composante commune sur l'ensemble des sections sismiques. En complément à ces études, d'autres points ont été évoqués, par exemple la correspondance entre procédés géophysiques et géostatistiques, ou la prise en compte d'incertitude sur les amplitudes. / In the framework of researches carried out by ANDRA (the French agency for nuclear waste), which aim at evaluating the feasibility of creating a repository of radioactive waste in the Callovo Oxfordian clay formation. This thesis studies two key steps to obtain a 3D image of the formation: its geometry and its internal heterogeneity. The geometry of the sedimentary formation is determined from seismic time picking along profiles. These profiles are oriented along two main directions. At the intersections, the values registered by two profiles are different. The low number of intersections does not allow a precise study of the spatial behavior of these differences; they are analyzed via the pseudo-cross-variogram. The variography analysis allows the choice of a coherent model for all the profiles, from which an estimation of the corrected time is presented. The internal heterogeneity is performed through the study of the impedance, a variable linked to the intrinsic properties of the rock. Two types of measurements are available. On one hand, the seismic impedance resulting from the geophysical inversion of the amplitude recorded along seismic reflection profiles is densely distributed in the area of interest. On the other hand, the impedance log, resulting from the product of rock density and P wave velocity, is known at only few wells. Bivariate spatial behavior along the vertical direction has allowed the formulation of the following hypothesis: each impedance measurement can be decomposed into a sum of a variable Z (common to both measures) and a residual specific to each measure. Under certain supplementary hypotheses, it is possible to estimate the common component along the seismic sections. Some complements to these studies are mentioned, for example the correspondence between geophysical and geostatistical methods or the consideration of uncertainty in the amplitudes.

La question des inégalités environnementales et sociales de santé (IESS) à l’échelle des territoires est rapportée à deux dimensions cumulatives : un niveau d’exposition des populations à un environnement de vie dégradé et un niveau de vulnérabilité face à ces facteurs de risques, capables d’affecter la santé. Ces travaux de recherche concernent les Systèmes d’Information Géographique (SIG) appliqués au domaine de la santé environnementale. Ils ont pour objectifs généraux la caractérisation de la qualité des milieux environnementaux et l’évaluation des IESS. L’atteinte de ces objectifs passe préalablement par la récolte, la genèse et la mise en forme de bases de données environnementales spatialisées. Ces bases issues de la surveillance physico-chimique et de la biosurveillance des milieux sont ensuite cartographiées à l'aide des outils SIG, notamment les méthodes géostatistiques d’interpolation spatiale. En parallèle, les disparités spatiales d'incidence de pathologies sont étudiées grâce à des méthodes de cartographie des maladies (Ratios d'Incidence Standardisés : SIR) et de détection de clusters atypiques d'évènements (statistique de scan) mises en œuvre sur la base de registres de maladies. Finalement des études écologiques spatialisées sont élaborées afin de rapprocher les cartographies environnementales générées aux indicateurs sanitaires et socio-économiques. Ces travaux s’attachent à répondre à la question « les populations présentant un état de santé dégradé vivent-elles dans un environnement dégradé ? » Ce questionnement a été étudié à travers trois axes de recherche. Le premier concerne la caractérisation de l’imprégnation de l’environnement par les éléments traces métalliques (ETM) et l’évaluation des IESS à fine résolution sur trois bassins de vie de la région Nord-Pas de Calais (NPdC). Cette recherche est menée à partir de mesures d’imprégnations biologiques réalisées à la fois chez des lichens épiphytes et chez l’Homme et d’un indicateur de défaveur localisé. Les ETM mesurés ont été considérés de manière individuelle, mais également de manière globale à travers le développement d’un indice intégré multimétallique, afin d’exprimer le niveau d’imprégnation général par plusieurs métaux. Des inégalités environnementales ont été mises en évidence à l’échelle des quartiers du bassin de Dunkerque. Nos résultats supposent également l’implication de l’environnement dans l’imprégnation des populations pour certains ETM.Dans le deuxième axe, nous avons mis en évidence des disparités spatiales d’incidence de l’insuffisance rénale chronique terminale entre les cantons du NPdC et souligné l’implication du statut socio-économique dans l’apparition de ces disparités, contrairement à d’autres facteurs (diabète, maladies cardiovasculaires, pratiques médicales). Seule une partie des disparités de cette pathologie est aujourd’hui expliquée. Il est donc nécessaire de poursuivre sur le rôle des contaminants de l’environnement. Enfin, le troisième axe s’intéresse à l’analyse spatiale et spatio-temporelle de contaminations des eaux souterraines par les ETM, en vue d’identifier de potentiels facteurs de risques environnementaux des maladies inflammatoires chroniques de l’intestin. La réalisation de ces travaux s’appuie sur plusieurs collaborations établies avec le Réseau REIN, le Registre EPIMAD et des équipes de recherche (EA4483 et EA2694 Université Lille 2, TVES EA4477 ULCO).A la suite de ces travaux de thèse, les perspectives de recherche sont de poursuivre le développement d’indicateurs intégrés d’exposition aux contaminations multiples des milieux environnementaux. Les résultats montrent également un manque de renseignement des bases de données environnementales par rapport aux registres de santé. Un travail est donc nécessaire afin de définir le contenu de ces bases, nécessaire pour caractériser la qualité des milieux et faciliter l’étude de l’interaction de l’Homme avec son environnement. / Environmental and social inequalities in health (ESIH) over territories are related to two cumulative dimensions: populations exposed to their living poor quality environment and the vulnerability of these populations to the environmental risk factors, which can affect health. This research deals with the Geographic Information Systems (GIS) applied to the field of environmental health. General purposes are the characterization of environmental media quality and the assessment of ESIH.Achieving these objectives requires a first step of harvest, genesis and formatting spatialized environmental databases. Such data are resulting from physico-chemical monitoring and biomonitoring. They were then mapped using GIS tools, including geostatistical spatial interpolation methods. On the over hand, spatial variability in the incidence of diseases were investigated using disease mapping methods (Standardized Incidence Ratios: SIR) and the detection of atypical clusters of events (scan statistics), which are based on disease registries. Finally, geographical ecological studies are developed to associate the environmental maps generated to health and socio-economic status. Thus, this work aims to answer the question \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\"do people with poor state of health live in a poor quality environment?\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\" This question has been studied through three main researches.The first relates the characterization of trace elements burdens in the environment and the assessment of ESIH on neighborhood scale, over three territories in the Nord-Pas de Calais (NPdC) region. This research is conducted from measurements of biological burdens performed both in epiphytic lichens and humans and from a localized index of deprivation. The measured metals were considered individually, but also holistically by developing an integrated multimetallic index, in order to describe the general status of environmental pollution by metals. Environmental inequalities were observed on neighborhood scale in Dunkerque. Our results assume that trace elements burdens in populations are affected by environmental burdens.In our second research we revealed spatial disparities in the incidence of end stage renal disease (ESRD) on small area in the NPdC. Unlike other factors (diabetes, cardiovascular disease, medical practices), we highlighted the role of socio-economic status in the occurrence of such disparities. Only a part of the ESRD variability is currently explained. It is therefore necessary to focus on the environmental hypothesis.The third research focuses on the spatial and spatio-temporal analysis of groundwater contamination by trace elements, in order to identify potential environmental risk factors in the incidence of chronic inflammatory bowel disease.This work is based on several collaborations with the REIN network, the EPIMAD registry, and several research teams (EA4483 and EA2694 Université Lille 2, TVES EA4477 ULCO).Following this thesis, the research prospects are to pursue the development of integrated indicators to assess population exposure to the multiple environmental media contamination. The results also indicate a lack of information in environmental databases compared to health registries. A work is thus needed to define the content of such databases. These are necessary to characterize the environmental quality and to help the assessment of interaction between the populations and their living environment.

Cette étude est consacrée à la modélisation spatio-temporelle de la pollution atmosphérique urbaine en utilisant un ensemble de méthodes statistiques exploitant les mesures de concentrations de polluants (NO2, O3) fournies par un réseau de surveillance de la qualité de l'air (AIRPARIF). Le principal objectif visé est l'amélioration de la cartographie des champs de concentration de polluants (le domaine d'intérêt étant la région d'Île-de-France) en utilisant, d'une part, des méthodes d'interpolation basées sur la structure spatiale ou spatio-temporelle des observations (krigeage spatial ou spatio-temporel), et d'autre part, des algorithmes, prenant en compte les mesures, pour corriger les sorties d'un modèle déterministe (Filtre de Kalman d'Ensemble). Les résultats obtenus montrent que dans le cas du dioxyde d'azote la cartographie basée uniquement sur l'interpolation spatiale (le krigeage) conduit à des résultats satisfaisants, car la répartition spatiale des stations est bonne. En revanche, pour l'ozone, c'est l'assimilation séquentielle de données appliquée au modèle (CHIMERE) qui permet une meilleure reconstitution de la forme et de la position du panache pendant les épisodes de forte pollution analysés. En complément de la cartographie, un autre but de ce travail est d'effectuer localement la prévision des niveaux d'ozone sur un horizon de 24 heures. L'approche choisie est celle mettant en œuvre des méthodes de type réseaux neuronaux. Les résultats obtenus en appliquant deux types d'architectures neuronales indiquent une précision correcte surtout pour les 8 premières heures de l'horizon de prédiction / This study is devoted to the spatio-temporal modelling of air pollution at a regional scale using a set of statistical methods in order to treat the measurements of pollutant concentrations (NO2, O3) provided by an air quality monitoring network (AIRPARIF). The main objective is the improvement of the pollutant _elds mapping using either interpolation methods based on the spatial or spatio-temporal structure of the data (spatial or spatiotemporal kriging) or some algorithms taking into account the observations, in order to correct the concentrations simulated by a deterministic model (Ensemble Kalman Filter). The results show that nitrogen dioxide mapping based only on spatial interpolation (kriging) gives the best results, while the spatial repartition of the monitoring sites is good. For the ozone mapping it is the sequential data assimilation that leads us to a better reconstruction of the plume's form and position for the analyzed cases. Complementary to the pollutant mapping, another objective was to perform a local prediction of ozone concentrations on a 24-hour horizon; this task was performed using Artificial Neural Networks. The performance indices obtained using two types of neural architectures indicate a fair accuracy especially for the first 8 hours of prediction horizon

Le béluga du Saint-Laurent est une espèce menacée au Canada et protégée par la Loi sur les espèces en péril du Canada. La détermination des fonctions biologiques de ses habitats essentiels est nécessaire afin d’assurer le rétablissement de la population. Parcs Canada a entamé en 2009 un suivi des proies du béluga dans deux de ses aires de fréquentation intensive situées dans le Parc marin du Saguenay–Saint-Laurent : l’embouchure de la rivière Saguenay et la baie Sainte-Marguerite. L’étude de l’abondance et de la distribution des proies est réalisée par sondage hydroacoustique le long de transects à l’aide d’un échosondeur multifréquences. Un protocole d’observations systématiques du béluga est mené simultanément aux sondages hydroacoustiques à partir de sites terrestres. Le premier objectif de cette étude est de développer la méthodologie concernant le traitement, la classification et la cartographie des données hydroacoustiques échantillonnées. L’objectif principal consiste à déterminer si l’abondance et la distribution des proies pélagiques ont une influence sur l’utilisation de ces deux habitats par le béluga. La cartographie de la biomasse relative de poissons a été réalisée pour la couche de surface, la couche en profondeur et pour l’ensemble de la colonne d’eau par krigeage ordinaire pour les deux habitats pour les 29 transects. À la baie Sainte-Marguerite, le nombre de bélugas observés augmente avec la biomasse relative des proies en surface et en profondeur. À l’embouchure de la rivière Saguenay, les résultats n’ont pas été concluants. Les résultats suggèrent que l’alimentation pourrait être l’une des fonctions biologiques de la baie Sainte-Marguerite. / The St. Lawrence beluga is a threatened species in Canada and protected by the Species at Risk Act in Canada. Determining the biological functions of its critical habitat is needed to ensure the recovery of the population. Parks Canada began in 2009 a monitoring of beluga preys in two of its areas of intensive use located in the Saguenay–St. Lawrence Marine Park: the Saguenay River mouth and the Sainte-Marguerite bay. The study of the abundance and distribution of prey is achieved by hydroacoustic survey along transects using a multifrequency echo sounder. A systematic observation protocol of the beluga is conducted simultaneously with hydroacoustic surveys from land-based stations. The first objective of this study is to develop the methodology for the processing, classification and mapping of these hydroacoustic data. The main objective is to determine whether the abundance and distribution of pelagic prey influence the use of these two habitats by beluga. The mapping of the fish relative biomass was conducted for the surface layer, the depth layer and the entire water column by ordinary kriging for the two habitats for the 29 transects. At Sainte-Marguerite bay, the number of whales observed increases with fish relative biomass in the surface and deep layers. At the Saguenay River mouth, the results are not conclusive. The results suggest that feeding may be one of the biological functions of the Sainte-Marguerite bay.

Ces travaux présentent des avancées théoriques pour l'application de l'approche EDPS (Équation aux Dérivées Partielles Stochastique) en Géostatistique. On considère dans cette approche récente que les données régionalisées proviennent de la réalisation d'un Champ Aléatoire satisfaisant une EDPS. Dans le cadre théorique des Champs Aléatoires Généralisés, l'influence d'une EDPS linéaire sur la structure de covariance de ses éventuelles solutions a été étudiée avec une grande généralité. Un critère d'existence et d'unicité des solutions stationnaires pour une classe assez large d'EDPSs linéaires a été obtenu, ainsi que des expressions pour les mesures spectrales associées. Ces résultats permettent de développer des modèles spatio-temporels présentant des propriétés non-triviales grâce à l'analyse d'équations d'évolution présentant un ordre de dérivation temporel fractionnaire. Des paramétrisations adaptées de ces modèles permettent de contrôler leur séparabilité et leur symétrie ainsi que leur régularité spatiale et temporelle séparément. Des résultats concernant des solutions stationnaires pour des EDPSs issues de la physique telles que l'équation de la Chaleur et l'équation d'Onde sont présentés. Puis, une méthode de simulation non-conditionnelle adaptée à ces modèles est étudiée. Cette méthode est basée sur le calcul d'une approximation de la Transformée de Fourier du champ, et elle peut être implémentée de façon efficace grâce à la Transformée de Fourier Rapide. La convergence de cette méthode a été montrée théoriquement dans un sens faible et dans un sens fort. Cette méthode est appliquée à la résolution numérique des EDPSs présentées dans ces travaux. Des illustrations de modèles présentant des propriétés non-triviales et reliés à des équations de la physique sont alors présentées. / This dissertation presents theoretical advances in the application of the Stochastic Partial Differential Equation (SPDE) approach in Geostatistics. This recently developed approach consists in interpreting a regionalised data-set as a realisation of a Random Field satisfying a SPDE. Within the theoretical framework of Generalized Random Fields, the influence of a linear SPDE over the covariance structure of its potential solutions can be studied with a great generality. A criterion of existence and uniqueness of stationary solutions for a wide-class of linear SPDEs has been obtained, together with an expression for the related spectral measures. These results allow to develop spatio-temporal covariance models presenting non-trivial properties through the analysis of evolution equations presenting a fractional temporal derivative order. Suitable parametrizations of such models allow to control their separability, symmetry and separated space-time regularities. Results concerning stationary solutions for physically inspired SPDEs such as the Heat equation and the Wave equation are also presented. A method of non-conditional simulation adapted to these models is then studied. This method is based on the computation of an approximation of the Fourier Transform of the field, and it can be implemented efficiently thanks to the Fast Fourier Transform algorithm. The convergence of this method has been theoretically proven in suitable weak and strong senses. This method is applied to numerically solve the SPDEs studied in this work. Illustrations of models presenting non-trivial properties and related to physically driven equations are then given.

Les grands réaménagements urbains génèrent d’importants volumes de terres excavées, dont la valorisation est désormais possible si elles sont compatibles avec le fond pédo-géochimique du site receveur. Afin de faciliter cette valorisation, cette thèse s’intéresse à des méthodologies de traitement de données acquises à d'autres fins que la détermination du fond pédo-géochimique, afin de fournir un référentiel de qualité géochimique des sols urbains à l'échelle d'un quartier ou d'une ville. Les premiers tests montrent les limites des calculs statistiques, usuellement utilisés avec des données dédiées au fond pédo-géochimique du fait de leur sensibilité à la limite de quantification, à la loi de distribution ou encore au pourcentage d'anomalie. Le découpage du territoire en entités géographiques cohérentes à partir de la nature des sols et sous-sols et des pressions anthropiques (actuelles ou passées) permet la constitution de différents niveaux de fond suivant la qualité géochimique des sols, mais une des difficultés concerne le renseignement de l'ensemble des entités. La classification statistique des échantillons, même spatialisée, n'apporte, actuellement pas de résultats exploitables pour la gestion des terres excavées. Enfin, le filtrage géostatistique via un modèle linéaire de corégionalisation présente une alternative solide pour le filtrage des anomalies et l'estimation du fond à l'échelle d'un quartier. / Large urban redevelopments generate large volumes of excavated soils whose reusability is now possible if they are compatible with the pedo-geochemical background of the receiving site. In order to facilitate this valorization, this thesis focuses on processing methodologies on data acquired for other purposes than the determination of the pedo-geochemical background, in order to provide a geochemical quality reference system for urban soils at the scale of one neighborhood or city. The first tests show the limits of the statistical stresholds, usually used with data dedicated to the pedo-geochemical background because of their sensitivity to the limit of quantification, to the distribution law or to the percentage of anomaly. The segmentation of the territory into coherent geographical entities based on the nature of the soil and subsoil and the anthropogenic pressures (current or past) allows the constitution of different geochemical background levels according to the geochemical quality of the soils, but one of the difficulties concerns the information of all entities. The statistical classification of the samples, even if spatialized, does not currently provide usable results for the management of the excavated earth. Lastly, geostatistical filtering via a linear model of coregionalisation provides a solid alternative for anomaly filtering and background estimation at a neighborhood scale.

La direction des vents nourriciers, obtenue à partir des gradients granulométriques et des morphologies dunaires, semble avoir changé au cours du temps. Celle-ci est orientée O-E au Pléniglaciaire moyen, puis NO-SE au Pléniglaciaire supérieur, tandis qu’au Dryas récent la direction est plutôt SSO-NNE. Dans ce travail, l’hypothèse d’un changement de la saison à laquelle les sables sont remobilisés est privilégiée. Ainsi, par comparaison avec la distribution azimutale actuelle des vents à l’échelle de l’année et en accord avec les reconstitutions pour le LGM, il est suggéré qu’au Pléniglaciaire supérieur la déflation éolienne était décalée vers la saison estivale, c’est-à-dire après la phase de dégel des sols et de fonte de la neige. Au contraire, les orientations établies pour le Pléniglaciaire moyen et le Dyras récent refléteraient une déflation active en période hivernale. / Le Pleniglacial, the north-west during the Upper Pleniglacial, and from south-southwest throughout the Younger Dryas. Here, we connect changing wind directions to seasonal variations in aeolian transport. Comparing annual azimuthal distribution of current wind regimes with wind regime models for LGM suggests that during the Upper Pleniglacial the dominant aeolian transport season shifted towards the summer, i.e. after seasonally frozen ground thawed and snow cover melted. Conversely, Middle Pleniglacial and Younger Dryas wind regimes reflect a general deflation during the winter.

La simulation numérique est un outil essentiel pour mieux comprendre la réalité des phénomènes qui se produisent au sein des milieux discrets, en particulier des milieux rocheux fracturés. Notre travail a pour objectif d'étudier le comportement mécanique et dynamique des milieux rocheux en réalisant des modèles discontinus dans lesquels les mesures de l'orientation et l'espacement interfractura, effectuées sur le terrain, sont intégrées au modèle de façon stochastique. Pour la génération des milieux fracturés numériques des algorithmes basés sur les méthodes statistiques et géostatistiques sont développés. Le comportement mécanique de modèles académiques et de structures réelles en maçonnerie est étudié grâce à la méthode Non Smooth Contact Dynamics (NSCD). Des structures en maçonneries telles que le monument historique de l'amphithéâtre de Nîmes et l'aqueduc d'Arles sont étudiées, à échelle réelle, en deux et trois dimensions pour leurs comportements dynamiques pendant une sollicitation sismique. Les modèles de milieux rocheux fracturés engendrés à partir d'études de cas sont analysés selon deux axes : l'étude de la stabilité de pentes rocheuses et la conception de travaux souterrains par la méthode NSCD en utilisant le code LMGC90. La stabilité de ces modèles est étudiée en deux et trois dimensions. Les analyses effectuées sur ces modèles ont permis de proposer des recommandations pour l'installation d'éventuels soutènements. Elles ont également permis la distinction des familles des fractures jouant un rôle prépondérant dans l'instabilité de ces massifs. La modélisation de certaines méthodes de protection en travaux au rocher telles que le boulonnage et les filets pare-pierres est ensuite abordée dans les perspectives envisagées comme développement de nos travaux.

Actuellement, la distribution spatiale des thons tropicaux pêchés par les senneurs français dans l'Océan Indien est représentée à partir des données de captures et d'effort issues des livres de bord. Or ces données ne comportent qu'une position par jour et le nombre de calées réalisées pendant cette journée : l'ensemble des zones prospectées sans action de pêche, synonyme d'absence d'agrégations de thon n'est pas connu. Ces zones seraient pourtant utiles à la cartographie des abondances. Cette flottille étant équipée de systèmes VMS (Vessel Monitoring Systems), les trajectoires données par les VMS sont utilisées pour déduire du comportement des navires, les zones de présence (pêche) et d'absence (non pêche) d'agrégations de thonidés. L'objectif est donc d'estimer l'activité d'un navire à partir de sa trajectoire. L'application d'un modèle markovien permet de distinguer les activités de route, de recherche, et d'arrêt, et de pêche. Pour la majorité des sorties en mer (exceptées celles où un observateur est à bord), l'activité réelle du navire n'est pas connue (variable latente), ce qui justifie l'approche par chaînes de Markov à états cachés. Le modèle a été développé dans un contexte bayésien incluant des priors sur certains des paramètres (vitesses, angles de changement de cap et matrice de transition). Il a été calibré et validé grâce aux données collectées par des observateurs embarqués sur 10% de la flottille. Le taux d'erreur est de 10% pour des positions VMS enregistrées toutes les heures. Une fois estimée, l'activité des senneurs en chaque point de leur trajectoire a permis (1) d'obtenir des indicateurs spatiaux et temporels de l'effort de la flottille de senneurs français, et (2) d'étudier la structure spatiale des zones de présence de bancs de thons (variogrammes d'indicatrices) et d'en déduire, par cokrigeage d'indicatrices, des cartes de présence/présence potentielle/absence de thon capturable, pouvant être interprétées comme indices d'abondance de thon.

Les enjeux économiques de l'industrie pétrolière imposent un important besoin de décrire la structure géologique des réservoirs d'hydrocarbures et des caractéristiques pétrophysiques des roches qui les constituent. La modélisation des réservoirs consiste à construire des modèles numériques qui représentent les hétérogénéités dans le réservoir à différentes échelles. Le schéma booléen est un modèle très utilisé pour établir un modèle lithologique de l'architecture interne du réservoir, que ce soit à l'échelle d'une unité génétique ou à l'échelle granulométrique à l'intérieur de ce dernier. Ce modèle présente néanmoins un problème important lors de son application à la modélisation de réservoirs présentant une distribution non stationnaire des lithofaciès qui les composent. Les proportions de ces lithofaciès fournissent une information sur leur distribution spatiale et constituent un outil très efficace pour déterminer la présence d'une non stationnarité. Elles permettent de quantifier l'information géologique des dépôts et constituent ainsi une importante contrainte à respecter pour obtenir des simulations géologiquement réalistes. Dans cette thèse, nous avons intégré dans le modèle booléen cette information contenue dans les proportions. Le problème abordé est l'inférence des paramètres du modèle à partir des proportions. Deux aspects ont été traités. D'une part, la définition du problème d'inférence du point de vue théorique, en établissant le rapport entre les variables expérimentales et les paramètres constitutifs du modèle. D'autre part, le développement d'un outil informatique opérationnel permettant d'automatiser le passage entre les données de proportion et les paramètres du modèle. Nous proposons une méthode d'inférence, basée sur un processus de déconvolution, qui permet d'introduire l'information des proportions dans le modèle et ainsi de reproduire dans les simulations une distribution non stationnaire des hétérogénéités de l'unité sédimentaire. Ces simulations respectent ainsi par construction la distribution latérale et verticale des proportions dans l'unité de même que la valeur moyenne globale de la proportion sur tout le volume étudié. Celle-ci est associée à des paramètres pétrophysiques telle que la porosité, et il est important de la retrouver dans les simulations.