L’objectif de ce travail est double : d’une part, il s’agit de développer une méthodologie intégrée pour la construction d’un modèle géomécanique ainsi que la représentation des incertitudes associées aux propriétés poro-élastiques des roches constitutives, en exploitant l’ensemble des données disponibles et en s’appuyant de façon cohérente sur les modèles de réservoir statique et dynamique classiquement utilisés par les géologues et les ingénieurs réservoir ; d’autre part, il s’agit d’analyser quel est l’impact des hétérogénéités géologiques, souvent négligées, dans la réponse mécanique du réservoir sollicité par son exploitation, et d’aboutir à des incertitudes sur les champs de contraintes et de déplacements, issues des incertitudes sur ces hétérogénéités et leurs paramètres mécaniques. Pour ce faire, une méthodologie intégrée s’appuyant sur des simulations géostatistiques a été développée. Après la construction du cadre géométrique 3D, le remplissage des propriétés au sein du réservoir suit une démarche de simulations géostatistiques 3D emboîtées, dans laquelle la représentation des hétérogénéités lithologiques conditionne la génération des propriétés poro-élastiques. La démarche consiste ensuite à représenter les incertitudes sur le modèle géomécanique par des ensembles de réalisations géostatistiques dont la réponse mécanique est alors calculée avec un simulateur mécanique aux éléments finis. Les incertitudes sur les champs de contraintes et de déformations sont déduites ensuite des différentes réponses mécaniques obtenues. La démarche a été mise en œuvre sur un réservoir réel, dans un environnement fluvio-deltaïque, produisant en Mer du Nord. Dans ce cadre, il a été démontré que les hétérogénéités du réservoir et leurs incertitudes influencent significativement les calculs des champs de contraintes et de déformations, ainsi que les risques mécaniques de rupture. Des incertitudes sur les quantités mécaniques analysées (premier invariant du tenseur des contraintes et subsidence) ont été aussi estimées / This work has two main objectives. The first one is to develop an integrated methodology allowing to build a 3D geomechanical model and also to image the uncertainties attached to the poro-mechanical properties of the constitutive rocks. This geomechanical model should be based on all related available data and should be consistent with the static and dynamic models, currently built by reservoir geologists and engineers. The second objective is to analyse the impact of geological heterogeneities, which are often neglected, in the mechanical response of the reservoir induced by its exploitation, and furthermore to derive uncertainties on the stress and deformation fields related to the uncertainties on the input properties of the geomechanical model. An integrated methodology based on geostatistical simulations is developed. First, the geometric frame is built; then an approach of embedded stochastic simulations is carried out to infill the different reservoir properties, the lithological description constraining the petrophysical and poro-elastic descriptions. The next step is to generate the mechanical responses of the stochastic realisations, using a finite-element mechanical simulator. The uncertainties on the resulting stress and displacement fields are then deduced from the multiple mechanical responses which are computed. This approach is demonstrated on a real field case, a fluvio-deltaic reservoir in North Sea. It is shown on this example that the reservoir heterogeneities and their uncertainties significantly influence the calculations of stress and strain fields, and also the risks of mechanical failure. Uncertainties on the mechanical quantities under analysis (first invariant of the stress tensor and subsidence) are also derived

Ce travail entre dans le cadre de l’étude des interactions eau-roche dans le cas du stockage du CO2 en milieu géologique. Un intérêt particulier est accordé aux hétérogénéités des paramètres associés aux phénomènes géochimiques. Ces hétérogénéités peuvent s’observer à différentes échelles: celle des grains (les minéraux présentent des défauts de cristallinité et des impuretés), et l’échelle centimétrique/pluri-décamétrique. En particulier, les paramètres thermodynamiques (logK) et de cinétique chimique (dans ce travail nous avons considéré le produit de la constante cinétique k par la surface spécifique S soit kS comme "paramètre de cinétique chimique") sont connus à partir des expériences de laboratoire pour des échantillons de quelques centimètres de dimension, alors que l’on s’intéresse aux réactions minéralogiques à l’échelle des réservoirs.Nous avons évalué les caractéristiques géostatistiques de la variabilité spatiale après réaction à travers des simulations de transport réactif dans lesquelles différents paramètres (logK et kS) sont perturbés avec une première variabilité imposée. Une combinaison de deux approches est ainsi abordée : déterministe et géostatistique. Le code du transport-réactif COORES (IFP-EN et Ecole nationale supérieure des mines de Saint-Etienne) a été utilisé pour les simulations déterministes et le système géochimique étudié concerne la dissolution du diopside avec précipitation de minéraux secondaires comme la calcite et la magnésite.Après analyse par la méthode des plans d’expériences, les résultats montrent qu’une corrélation spatiale élevée combinée avec une grande variance de dispersion des minéraux favorise une réactivité importante des minéraux lorsqu’on perturbe le paramètre de cinétique chimique kS. Par ailleurs une vitesse d’injection élevée accélère le processus de dissolution du minéral étudié. La variabilité spatiale du paramètre thermodynamique n’a cependant pas d’effet significatif sur les résultats, le système se comporte comme dans le cas homogène. Du point de vue de l’homogénéisation du paramètre kS, on retrouve l’influence de l’historique de dissolution. / The present work is based on the study of water-rock interactions in the case of CO2 storage in geological media. Particular attention is devoted to heterogeneities at different observation scales geochemical phenomena. These heterogeneities can be observed at different scales: the grain (mineral crystallinity present defects and impurities), and the centimeter scale / multi- decametric (rocks are heterogeneous at different scales). In particular, the thermodynamic parameters logK and chemical kinetics kS (in this work we considered the product of the rate constant k by the specific surface area S is kS as "chemical kinetics parameter") are known from laboratory experiments to a few centimeters in size, while we are interested in mineralogical reactions across tanks.We propose to evaluate the geostatistical characteristics of the local variability after reaction through simulations of reactive transport on a small scale in which various parameters (logK and kS) are perturbed with a first spatial variability imposed. A combination of both approaches is discussed: deterministic and geostatistical for the study of geochemical problems at different scales. The reactive transport code - COORES (IFP - EN and Ecole nationale supérieure des mines de Saint -Etienne) was used for deterministic simulations and the geochemical system studied concerns the dissolution of diopside with precipitation of secondary minerals such as calcite and magnesite.After analysis by the method of design of experiments, the results show that high spatial correlation variance combined with high dispersion of minerals promotes a high reactivity when minerals chemically disturbing is the kinetic parameter kS. In addition, a high velocity injection accelerates the dissolution of the mineral studied. However, the effect of spatial variability of the thermodynamic parameter, did not significantly affect the results, the system behaves as in the homogeneous case. From the standpoint of homogenizing the parameter kS, include the influence of the history of dissolution.

Dans le cadre de la gestion de la ressource en eau, la connaissance des hétérogénéités des nappes alluviales représente un enjeu primordial. Ce travail de thèse, qui s’inscrit dans le programme du PIREN-Seine, s’est intéressé à la plaine alluviale de la Bassée, à l’amont de Paris, la plus large plaine alluviale de la Seine.Cette plaine alluviale joue le rôle d’une zone tampon contrôlant les flux entre le réseau de surface et un aquifère régional. Connaître les hétérogénéités sédimentaires est primordial pour la compréhension des écoulements et fait l’objet de ce travail.L’architecture à l’échelle kilométrique des dépôts alluviaux et la chronologie des nappes associées est établie. La géométrie du contact nappe-substrat, une vallée incisée, est estimée par krigeage. Une restitution paléo-hydrologique s’appuie sur les paléo-tracés préservés en surface des nappes alluviales. L’architecture fine des dépôts est restituée à partir de données de sondages converties en propriétés pétrophysiques et d’imagerie géoradar.La plaine de la Bassée correspond à deux terrasses étagées (Saalien terminal et Weichsélien) dans lesquelles sont emboîtées deux autres terrasses (Tardiglaciaire et Holocène).L’estimation de la géométrie du contact alluvions-substrat suggère la présence d’une rupture de pente liée à une érosion régressive durant le Weichsélien. Les faciès des différentes nappes sont assez semblables et dominés par des dépôts de forte perméabilité (sables et graviers). De faibles valeurs de perméabilités sont associées aux faciès de colmatage de chenaux et dépôts en masse.Les résultats obtenus sur la plaine de la Bassée méritent d’être confrontés et validés sur des affluents de la Seine (Marne, Yonne, Oise…) afin de proposer un modèle cohérent à l’échelle du bassin. / In the scope of water resource management, the characterization of alluvial plain architecture and heterogeneity represents a challenge. This PhD work which is part of the PIREN-SEINE program focuses on largest alluvial plain of the River Seine, “La Bassée”, located upstream of Paris.This alluvial plain can be considered as a buffer layer that controls the water fluxes between the fluvial system and the regional aquifers. Determining the sedimentary heterogeneities is essential for the comprehension of fluxes, this is the aim of this work.Kilometric scale architecture of the alluvial deposits and the associated alluvial sheets chronology are established. The geometry of the sheet-bedrock contact, an incised valley, is estimated by kriging interpolation. A paleohydrological restitution is based on the paleochannels preserved at the top of the alluvial sheets. The fine scale architecture of deposits is restored from sampling data converted into petrophysics properties and ground-penetrating radar acquisitions.The alluvial plain of La Bassée corresponds to two terraces (terminal Saalian and Weichselian) in which two other terraces are staircased (Late Glacial and Holocene).The estimated sheet-bedrock contact geometry suggests the presence of a knick-point related to a regressive erosion during the Weichselian. Sedimentary facies are very similar from one alluvial sheet to another, and are dominated by high permeability deposits (sand and gravel). Low permeability facies are associated to channel plugs and mass transport deposits.The infilling model of La Bassée region obtained in this study should be compared and validated on the River Seine tributaries (Marne river, Yonne river, Oise river…) to propose a coherent model at the basin scale.

Actuellement, la restauration des oeuvres d’art, notamment des statues fracturées, repose sur des techniqueséprouvées, mais empiriques. Les statues endommagées comportent souvent des parties brisées. Leur restaurationconsiste la plupart du temps à les rassembler. Ainsi apparait la nécessité de prendre en compte lespropriétés mécaniques des interfaces entre les différentes parties brisées, ce qui permet de limiter l’ampleur desréparations et ainsi, de mieux conserver l’intégrité de l’oeuvre. Par ailleurs, les techniques numériques d’acquisition3D font leur entrée au service de la conservation du patrimoine. Cette thèse propose une méthodologiecapable d’utiliser des données issues d’une acquisition 3D pour simuler les opérations de restauration et leurseffets sur la structure de l’oeuvre. Les processus de restauration peuvent ainsi être testés et optimisés.Un scanner laser est utilisé pour l’acquisition de la géométrie des oeuvres, ce qui nous permet de reconstruireun modèle 3D pour la simulation numérique. Les calculs sont menés dans le cadre de la mécanique des milieuxcontinus déformables avec FLAC3D. Pour vérifier tous les points clés garantissant la stabilité mécanique, lecomportement des éléments de renforts et celui des interfaces entre les blocs ont été considérés. À partirdes résultats de ces études, une critique des stratégies de restauration mises en oeuvre ou envisageables estproposée.De plus, plusieurs méthodes de caractérisation visant à retrouver la provenance du matériau et/ou estimerles propriétés mécaniques de l’oeuvre sont proposées : caractérisations physico-chimiques et minéralogiques,essais non destructifs et destructifs. Une campagne expérimentale visant à caractériser le comportement desfractures en contact frottant avec acquisition de l’état de surface a été réalisée. Une analyse des corrélationsentre les propriétés mécaniques et morphologique des interfaces est ensuite élaborée. Enfin nous proposonsdes modèles prédictifs construits par régressions linéaires multiples et multivariées. Cette étude permet desimuler le comportement d’une oeuvre fracturée. / Currently, the artworks restoration, including broken statues, are based on proven but empirical techniques. Damaged statues often have broken parts. Restoration mainly consists in joining these parts. Thus appears the necessity of taking into account the interfaces mechanical properties between the different broken parts, reducing thereby the extent of repairs and thus better maintain the integrity of the work. Beside this, digital 3D acquisition techniques become more and more popular among the cultural heritage community. This thesis proposes a methodology able to use 3D acquisition datas to simulate the restoration operations and their effects on the structure. The restoration strategy can thus be tested and optimized. A laser scanner is used to acquire the statue geometry, allowing to reconstruct a 3D model for numerical simulation. Computations are carried out within the framework of deformable continuum mechanics using the FLAC3D software. In order to check all the key points ensuring mechanical stability, the behavior of reinforcing elements and those of the interfaces between blocks were considered. From the results of the present studies, a critical analysis of implemented or possible restoration strategy is proposed. In addition, several characterization methods to find the material origin and/or estimate the mechanical properties of the work are proposed : physico-chemical and mineralogical characterizations, non-destructive and destructive tests. An experimental campaign to characterize the behavior of fractures in contact with surface acquisition was performed. An analysis of correlations between mechanical and morphological interfaces properties is then developed. Finally, we propose a predictive modelling based on multiple linear and multivariate regressions. This study allow the simulation of fractured artworks behaviour.

Cette thèse a pour but de contribuer à la qualification du formulaire de calcul neutronique NARVAL, dédié aux coeurs de réacteurs à plaques. En particulier, l’objectif est de développer des méthodes innovantes permettant d’utiliser les données expérimentales inédites du programme HIPPOCAMPE pour évaluer la précision du profil de puissance calculé. La complexité provient de la localisation de l’instrumentation (chambres à fission placées entre les assemblages) et des hétérogénéités caractéristiques de ce type de coeurs (géométrie à plaques, poisons consommables et de contrôle solides). Pour aborder ce problème deux voies ont été mises en oeuvre : la première voie consiste à « combiner puis extrapoler » les écarts C/E observés afin de déterminer les incertitudes associées aux facteurs de puissance. Nous avons utilisé, pour ce faire, la méthode « P/A », traditionnellement employée dans les REP électrogènes mais jamais appliquée aux coeurs à plaques à ce jour. La deuxième voie passe, en revanche, par la reconstruction d’une nappe de puissance à utiliser comme référence (comparaison calcul/« expérience-reconstruite ») : nous avons focalisé notre travail sur des techniques géostatistiques. Après avoir constaté que les deux méthodes conduisent à des résultats satisfaisants (erreur comparable à l’incertitude expérimentale cible) nous avons continué notre recherche, en explorant les possibles développements et en introduisant en particulier une nouvelle méthode hybride (associant les techniques géostatistiques à la méthode P/A) qui permet d'améliorer ultérieurement la qualification du profil de puissance (écart-type des écarts C/E cohérent avec la constatation expérimentale). / The aim of this doctoral thesis work is to contribute to the experimental validation of a neutron physic code, called NARVAL, devoted to the analysis of slab cores reactors. The primary objective is to develop some innovative methods in order to validate the computed power map starting from the original experimental data, provided by the HIPPOCAMPE campaign. The particular position of the instrumentation (fission chambers located between the assemblies) and the strong heterogeneities, characterising this specific core design (slab geometry, burnable and control neutron absorbers in solid state) represent the main challenge of this work. Two different approaches are investigated : the first one consists in “combining and extrapolating” the observed calculated/experimental results in order to evaluate the uncertainty of power coefficients. Among different solutions, the “P/A” method is chosen : it is usually employed to perform conventional PWR plant analysis and has never been applied before to slab cores. The latter aims to reconstruct a power map that could be used as a direct reference for code validation : in this case the geostatistical techniques are selected. These methods provide satisfactory results as estimated errors are in good agreement with the experimental uncertainty target. Nevertheless, in this work a new hybrid method, applying the geostatistical technics to the P/A scheme, is investigated and developed. The good agreement between the experimental and the estimated validations of the computed power map attests the noteworthy performance of this new method.

La problématique du calage d'historique en ingénierie de réservoir, c'est-à-dire le calage des modèles géostatistiques aux données de production, est un problème inverse mal posé. Dans un cadre bayésien, sa résolution suppose l'inférence de la distribution de probabilité du modèle géostatistique conditionné aux données dynamiques, rendant compte à la fois de l'a priori géologique, exprimé dans le modèle géostatistique, et de l'adéquation aux données de production. Typiquement, la résolution de ce problème passe par la génération d'un ensemble de réalisations calées aux données, échantillon représentatif de cette distribution. Les modèles géostatistiques sont en général discrétisés sur des grilles de plusieurs centaines de milliers, voire des millions de blocs ; les approches classiques tentent de résoudre le problème inverse en considérant l'ensemble des blocs comme paramètres du modèle. Sa dimension est alors considérable et les méthodes d'échantillonnages deviennent impraticables sur un cas réel. Il convient alors de choisir une paramétrisation susceptible de réduire la dimension du problème. Dans la première partie de cette thèse, nous présentons une méthode de paramétrisation optimale des modèles géostatistiques basés sur les champs aléatoires gaussiens, à partir de leur décomposition de Karhunen-Loève (KL). Nous en décrivons les fondements théo- riques, puis, sur des applications aux modèles de champs aléatoires gaussiens courants en géostatistique, selon des critères d'abord statistiques puis hydrodynamiques, nous quantifions la réduction de la dimension du problème offerte par cette paramétrisation. Dans la seconde partie, nous présentons les principes des méthodes de Monte-Carlo par Chaînes de Markov (MCMC) et les défauts des méthodes classiques pour la résolution du problème inverse dans le cadre bayésien. Nous développons alors l'approche par chaînes de Markov en interaction dont nous exposons les avantages. Enfin, les résultats obtenus par l'emploi conjoint de ces deux méthodes sont présentés dans deux articles. Une approche différente, passant par l'emploi de méthodes d'analyse de sensibilité, est également décrite dans un troisième article.

L'AMÉLIORATION DES TECHNIQUES D'OBSERVATION DES PRÉCIPITATIONS À L'ÉCHELLE RÉGIONALE EST L'UN DES ENJEUX FORTS DE L'OBSERVATOIRE HYDROMÉTÉOROLOGIQUE MÉDITERRANÉEN CÉVENNES-VIVARAIS. SUITE À UNE PREMIÈRE EXPÉRIMENTATION RÉALISÉE À BOLLÈNE À L'AUTOMNE 2002, LES DONNÉES RECUEILLIES PAR LE RADAR MÉTÉOROLOGIQUE OPÉRATIONNEL DE BOLLÈNE DU RÉSEAU ARAMIS ONT SERVI DE SUPPORT À UNE CHAÎNE D'ALGORITHMES MISE AU POINT POUR L'ESTIMATION QUANTITATIVE DES PRÉCIPITATIONS AU SOL. <br />NOTRE CONTRIBUTION À CES TRAITEMENTS CONCERNE PARTICULIÈREMENT UNE SOURCE D'ERREUR POUR LA MESURE DES PRÉCIPITATIONS PAR RADAR MÉTÉOROLOGIQUE : L'HÉTÉROGÉNÉITÉ VERTICALE DES PRÉCIPITATIONS. NOUS NOUS PLAÇONS DANS LA CONTINUITÉ DE LA MÉTHODE À BASE STATISTIQUE D'IDENTIFICATION DU PROFIL VERTICAL DE RÉFLECTIVITÉ PAR INVERSION PROPOSÉE PAR ANDRIEU ET CREUTIN (1995). NOUS AVONS ÉTENDU CETTE MÉTHODE AU CAS DE SUPPORTS GÉOGRAPHIQUES ÉVOLUTIFS DANS LE TEMPS AFIN DE CONSIDÉRER LE TYPE DE PRÉCIPITATIONS (CONVECTIF, STRATIFORME). <br />L'ÉVALUATION DES TRAITEMENTS, PAR LA CARACTÉRISATION DE L'INCERTITUDE DE L'ESTIMATION DE PLUIE PAR RADAR, EST UN ASPECT IMPORTANT DU TRAVAIL. NOUS AVONS DÉVELOPPÉ UNE TECHNIQUE BASÉE SUR LA GÉOSTATISTIQUE POUR ÉTABLIR UNE PLUIE DE RÉFÉRENCE À PARTIR DE DONNÉES PLUVIOMÉTRIQUES. NOUS AVONS ÉGALEMENT DÉVELOPPÉ UN MODÈLE POUR CARACTÉRISER LA STRUCTURE STATISTIQUE DES ERREURS AINSI QUE LEUR STRUCTURE SPATIO-TEMPORELLE. <br />NOUS AVONS ENFIN DÉVELOPPÉ UNE APPROCHE D'IDENTIFICATION DU PVR À BASE PHYSIQUE. EN DIMINUANT LE NOMBRE DE PARAMÈTRE PAR RAPPORT À LA MÉTHODE INITIALE ET EN INTRODUISANT DES CONTRAINTES PLUS PHYSIQUES, CETTE APPROCHE VISE À AMÉLIORER LA ROBUSTESSE D'IDENTIFICATION DU PVR.

Face à la contradiction consistant à traiter les variables régionalisées écologiques sans tenir compte de leurs propriétés spatiales, nous développons des méthodes géomatiques, utilisant des techniques informatiques, et géostatistiques, appliquant la théorie des fonctions aléatoires. Après avoir introduit des éléments de géomatique et de géostatistique, et avoir précisé la nature spécifique de l'autocorrélation spatiale, nous introduisons les inférences design-based et model-based. A l'aide de fonctions aléatoires, nous étudions l'efficacité de l'échantillonnage, optimisons des prédicteurs, et calculons des intervalles de prédiction. Nous proposons une procédure d'optimisation des classes de distances lors du calcul du variogramme. Nous examinons également l'utilisation de l'intégrale du variogramme, et justifions la modélisation du variogramme par ajustement aux moindres carrés pondérés. Nous discutons de la précision du variogramme dans les cadres design-based et model-based, et au sens du jackknife. L'optimisation de l'échantillonnage d'une population finie en vue de l'estimation de la moyenne spatiale ou du variogramme est examinée à l'aide de plusieurs heuristiques d'optimisation combinatoire et de simulations de fonctions aléatoires. Le problème du test de la corrélation ou de l'association entre deux variables régionalisées est étudié, à nouveau en utilisant des simulations de fonctions aléatoires. Nous passons en revue plusieurs méthodes et recommandons les tests qui font explicitement référence à l'autocorrélation spatiale des variables régionalisées. Dans le cadre de la définition de l'association spatiale entre variables régionalisées, nous proposons une méthode hybride utilisant des quadtrees et une distance d'édition entre arborescences récursives. Enfin, nous étudions des mesures de la complexité spatiale, critiquons l'analyse fractale et proposons des méthodes alternatives, notamment une mesure de complexité topologique d'une carte en isolignes.

Le point de départ de ce mémoire est la technique d'Analyse en Composantes Principales de Processus (ACPP), extension à des champs spatialement et/ou temporellement continus de l'Analyse en Composantes Principales (ACP) classique. Elle conduit à rechercher les fonctions propres d'un noyau de covariance, solutions d'une équation intégrale de Fredholm et permet ainsi de tenir compte explicitement de la forme et de la taille du domaine d'étude. Dans le chapitre l, après un rappel de la méthode et la présentation de l'approximation numérique utilisée lorsqu'on travaille avec un échantillon de taille fInie et un processus discrétisé spatialement, on s'attache à montrer les liens de cette technique avec la recherche de modes normaux de systèmes dynamiques puis avec les techniques d'interpolation optimale et la géostatistique. On montre ainsi que l'ACPP peut s'interpréter comme un cas particulier de krigeage s'appuyant sur un modèle de covariance spécifique, non nécessairement stationnaire. Dans le chapitre 2, afIn de quantifIer les erreurs sur l'estimation des fonctions propres, dues à l'échantillonnage statistique et à la discrétisation spatiale, une solution analytique, pour un modèle particulier de covariance et un domaine circulaire, est proposée puis comparée aux solutions numériques. Le chapitre 3 pose ensuite les principes de l'utilisation de l'ACPP pour la simulation de champs bidimensionnels et analyse la sensibilité de la méthode au choix des paramètres en terme de restitution d'une covariance cible. La technique est aussi comparée à d'autres méthodes (Bandes tournantes, méthodes matricielles), puis utilisée pour la simulation de champs non stationnaires. Enfin, le chapitre 4 illustre tous ces développements méthodologiques sur des données de température mensuelle de surface de la mer sur l'Atlantique intertropical entre 1964 et 1987.

L'insolation est étudiée, à l'échelle synoptique (la France) à l'aide d'un ensemble de mesures de la durée quotidienne d'insolation réalisées au sol et par télédétection satellitaire. L'objectif poursuivi est d'analyser une série de données suffi* samment longue pour caractériser de manière statistique le comportement climatologique de l'insolation dans l'espace. Les données sont fournies par 88 héliographes du réseau synoptique français ainsi que par le programme " Cactus " de la Météorologie Nationale, qui utilise les images visibles du satellite Météosat. L'étude porte sur une année de données. La similarité de la structure spatiale des mesures du sol et du satellite est montrée par analyse objective, et par analyse en composantes principales de processus. La fine perception de la variabilité spatiale du satellite est alors utilisée pour optimiser un réseau de mesures au sol. Les performances du programme Cactus sont comparées avec celles de méthodes d'interpolation classiques (basées sur 23 héliographes), pour l'estimation du rapport d'insolation. Elles se sont révélées plus performantes en terme de cofluctuation. Il apparaît cependant une tendance du satellite à exagérer les valeurs extrêmes. L'étalonnage des mesures satellitaires par des mesures sol a pour but de conjuguer les qualités de précision des hélio* graphes, et la définition spatiale du satellite. Une approche géostatistique est proposée. L'effet recherché est observé à partir d'une certaine densité du réseau d'étalonnage