Добірка наукової літератури з теми "Estimation géostatistique"

Les modélisations régionales proposées sur les modules annuels et les débits moyens d'étiage des sous-bassins de la Moselle française devraient contribuer à l'amélioration des connaissances sur le fonctionnement physique actuel des hydrosystèmes. Elles s'inscrivent dans le contexte des directives de la loi Pêche et de la loi sur l'Eau et plus récemment de la Directive Cadre Européenne. Quarante neuf sous-bassins répartis en trois jeux ont permis de caler et valider un modèle régional des modules annuels et un modèle régional d'étiage. Dans le cas du bassin hydrographique de la Moselle française une certaine dépendance existe entre modèles régionaux dont la loi statistique choisie est la loi de Weibull à 2 paramètres. Une pseudo-dépendance est observée entre la loi régionale des modules annuels et la loi régionale des débits moyens d'étiage pour les années moyennes à sèches. Cette propriété va permettre en particulier l'usage d'une procédure simplifiée commune, établie à partir de la connaissance de jaugeages épisodiques d'étiage, pour l'estimation des descripteurs de débit d'un sous-bassin non observé : le module médian /qa et le débit quotidien minimal médian /vcnd=1. Pour le modèle régional d'étiage un deuxième descripteur local est nécessaire. Il s'agit d'un temps caractéristique d'étiage du sous-bassin ∆e (j) permettant de généraliser le modèle à toute durée d. Le concept débit-durée-fréquence QdF appliqué aux étiages exploite la convergence observée des distributions de différentes durées d et est indépendant de la loi fréquentielle choisie. Le caractère opérationnel de ces modélisations régionales dépend essentiellement de la précision d'estimation des descripteurs de débit du sous-bassin étudié /qa, pour les modules annuels et /vcnd=1 pour les étiages. Ces descripteurs de débit ont été estimés selon deux approches : l'approche classique par régression multiple et selon une approche simple de recherche d'un coefficient de tendance k entre jaugeages épisodiques d'étiage concomitants au sous-bassin étudié (pas ou peu d'observations) et au sous-bassin de référence (chronique de débit continue). Pour cela, un choix de cinq jaugeages d'étiage par an sur les douze dernières années en moyenne a été fait. Le descripteur de débit du sous-bassin étudié est ensuite déduit du produit de k par le descripteur de débit du sous-bassin de référence. Pour /vcnd=1, sb. étudié nous observons dans la majorité des cas une nette amélioration de l'estimation obtenue par régression, notamment une forte réduction des écarts les plus importants. Une similitude des classes de superficie entre sous-bassin étudié et sous-bassin de référence n'est pas exigée. La proximité géographique des sous-bassins semble donner de meilleurs résultats. En ce qui concerne le module médian /qasb. étudié, son estimation par régression multiple est assez performante. Parallèlement à cela, le coefficient k de tendance permet, de même que pour, une estimation cohérente de /qasb. étudié. Ce résultat un peu inattendu laisse supposer que la pseudo-dépendance observée entre modèles régionaux a bien une réalité physique. Nous avons insisté sur cette démarche "de régionalisation" nécessitant un faible investissement en mesures de débit des sous-bassins non observés par rapport au réseau national de suivi hydrométrique. Elle se présente à notre avis comme une alternative (ou complémentarité) intéressante aux méthodes de régionalisation à bases géostatistiques : telles que l'identification du voisinage hydrologique homogène du sous-bassin étudié ou encore la prise en compte de l'effet structurant du réseau hydrographique dans la cartographie du descripteur de débit. L'ensemble des connaissances relative à cette recherche est repris dans un Système d'Information Géographique pour répondre éventuellement à la demande.

Ce gisement d'uranium en terrain sedimentaire, exploite par la cogema, est situe dans le sud de la france. Nous presentons d'abord les modeles utilises et leurs applications a une zone exploitee en mine a ciel ouvert. Pour proceder a l'estimation des ressources in situ, nous utilisons, a une teneur de coupure fixee, les variables utiles minerai et metal, le minerai representant la proportion mineralisee et le metal la richesse portee par cette proportion mineralisee. Ces variables sont manifestement liees entre elles. Le modele de coregionalisation des variables utiles permet d'effectuer une estimation coherente car prenant en compte, durant l'estimation, du lien entre ces deux variables. En considerant maintenant differentes teneurs de coupure, et en generalisant la demarche precedente aux variables utiles minerai a une teneur de coupure basse et minerai a une teneur de coupure haute, nous aboutissons au modele a residus d'indicatrices. Le principal interet de ce modele isofactoriel est de permettre un krigeage disjonctif, c'est-a-dire une estimation globale et locale des reserves avec changement de support. Nous voyons ensuite une couche uranifere exploitee par tailles en mine souterraine. L'exploitation d'une taille se fait de maniere selective, avec definition d'un premier tir. A la teneur de coupure d'exploitation et en tenant compte des contraintes minieres, des blocs sont regroupes et la zone de premier tir est definie. Dans cette etude, nous avons pu assimiler les valeurs krigees des blocs aux valeurs vraies, et voir le biais qu'entraine une estimation par polygone d'influence. A l'aide des sondages de reconnaissance, le modele a residus d'indicatrices autokrigeables permet d'effectuer la prevision de cette couche exploitee

Comment analyser la distribution spatiale d'espèces marines en relation avec quelques paramètres de l'environnement ? A l'aide de quatre exemples clefs qui permettent de faire le point sur les difficultés de l'analyse structurale en halieutique (densités de poissons, température de surface, courants marins), on discute des choix de modélisation. Il apparaît ainsi que la géostatistique transitive est une méthodologie adaptée quand l'abondance totale est le fait de quelques valeurs extrêmement fortes et qu'il est difficile de délimiter un domaine au sein duquel les variations de densités peuvent être supposées homogènes. De là, on choisit un individu au hasard dans la population et on s'intéresse à la probabilité qu'il se trouve en un certain point de l'espace, à une certaine température ou bien dans une partie de la population ou la densité vaut une certaine valeur (statistiques individuelles). Une description spatiale, basée sur le covariogramme, est obtenue par tirage au sort de deux individus. Des procédures d'estimations du covariogramme adaptées à des échantillonnages irréguliers sont proposées ainsi que les utilisations possibles d'un modèle structural, c'est-à-dire d'un covariogramme modélisé. En particulier, on indique en l'illustrant à partir de densités de merlu du Golfe de Gascogne (France), que le covariogramme contient des éléments structuraux plus fins que ne le laisse envisager sa définition, et on procède à une cartographie par krigeage transitif. L'approche par individu est utilisée pour la description et la modélisation de la répartition des oeufs de maquereau pendant la saison de ponte de 1983 à 1995. L'interprétation dans un contexte spatial des régressions expérimentales (en cloche) étant délicate, on préfère tester l'existence d'un preferendum en température de surface en vérifiant que les individus sont distribues de façon à minimiser leur inertie en température (inertiogrammes). De ce point de vue, les maquereaux apparaissent contraints dans leur régionalisation par des "fronts virtuels d'inertie"

Le but de cette thèse est d'exploiter au mieux les résultats du traitement des données sismiques en profondeur pour réaliser une simulation géostatistique de l'impédance acoustique. La base de la prospection géophysique est la sismique réflexion, une méthode qui fournit les coefficients de réflexion des interfaces géologiques. Ces coefficients fournissent des informations sur la structure du sous-sol et permettent d'estimer les propriétés pétro-acoustiques de chacune des couches géologiques. Nous utilisons une méthode d'imagerie (la migration) qui permet d'obtenir des données sismiques migrées en profondeur, dont les variations de réflectivité selon l'angle sont conservées. L'analyse de ces variations est alors effectuée sur les panneaux Iso-X le long de trajectoires définies par la variation de la réflectivité avec l'angle d'incidence. Ces trajectoires permettent d'obtenir une première interprétation structurale (le squelette) à laquelle sont associés des attributs tels l'intercept et le gradient. Un pointé sur ces deux squelettes complémentaires permet de délimiter des couches dans lesquelles la simulation sera effectuée. Après avoir interprété sur les sections sismiques et d'autres attributs des éventuelles anomalies en hydrocarbures, nous proposons de contraindre localement la simulation géostatistique pour reproduire ces anomalies. Pour cela, nous nous servons du lien qui existe entre l'impédance acoustique et l'intercept. Nous obtenons finalement une section en impédance acoustique en profondeur, valable sur la globalité de la section étudiée, et qui tient compte des anomalies d'hydrocarbures importantes pour l'éventuelle exploration et production du réservoir.

Le but de cette thèse est d’exploiter au mieux les résultats du traitement des données sismiques en profondeur pour réaliser une simulation géostatistique de l’impédance acoustique. La base de la prospection géophysique est la sismique réflexion, une méthode qui fournit les coefficients de réflexion des interfaces géologiques. Ces coefficients fournissent des informations sur la structure du sous-sol et permettent d’estimer les propriétés pétro-acoustiques de chacune des couches géologiques. Nous utilisons une méthode d’imagerie (la migration) qui permet d’obtenir des données sismiques migrées en profondeur, dont les variations de réflectivité selon l’angle sont conservées. L’analyse de ces variations est alors effectuée sur les panneaux Iso-X le long de trajectoires définies par la variation de la réflectivité avec l’angle d’incidence. Ces trajectoires permettent d’obtenir une première interprétation structurale (le squelette) à laquelle sont associés des attributs tels l’intercept et le gradient. Un pointé sur ces deux squelettes complémentaires permet de délimiter des couches dans lesquelles la simulation sera effectuée. Après avoir interprété sur les sections sismiques et d’autres attributs des éventuelles anomalies en hydrocarbures, nous proposons de contraindre localement la simulation géostatistique pour reproduire ces anomalies. Pour cela, nous nous servons du lien qui existe entre l’impédance acoustique et l’intercept. Nous obtenons finalement une section en impédance acoustique en profondeur, valable sur la globalité de la section étudiée, et qui tient compte des anomalies d’hydrocarbures importantes pour l’éventuelle exploration et production du réservoir
One of the main goal of geosciences in the petroleum industry is to build accurate and predictive reservoir models, expressed in terms of properties such as geological facies, porosity or acoustic impedance. The seismic reflection technique, at the base of geophysical prospecting, provides the reflection coefficients associated to the geological interfaces. These coefficients give information about the subsurface structures and allow us to estimate the layers petroacoustic properties. We use a preserved amplitude prestack depth migration which gives the depth migrated seismic data, where reflectivity variations versus angle are conserved. The variation analysis is done with the common image gathers sections along trajectories defined by the reflectivity variations with the angle. These trajectories provide a geology depth structural interpretation (the skeleton) with two major AVA attributes: intercept and gradient. A horizon picking phase using both skeletons defines the intervals where the simulation will be achieved. Then, after the interpretation of the hydrocarbons anomalies over the seismic and AVA attributes sections, we propose to locally constrain the geostatistical simulation in order to reproduce these anomalies. To handle it, we use the relationship between the acoustic impedance and the intercept. This method finally provides a depth acoustic impedance section, valid over the whole studied area, which contains the hydrocarbon anomalies for a possible oil and gas exploration & production

Jusqu'a recemment, les equations aux derivees partielles et la geostatistique sont des domaines qu'on a peu l'habitude d'associer. Pourtant, parmi les sciences de la terre, les phenomenes regis par des equations aux derivees partielles sont nombreux: on peut citer les ecoulements de fluides (air en meteorologie, eau en hydrogeologie) et les transferts de chaleur (geothermie). Ce type de problemes a ete, depuis peu, l'objet de nombreuses etudes sous l'angle probabiliste. Cependant, souvent dans ces travaux, seules les solutions stationnaires sont considerees et peu d'auteurs ont fait appel aux methodes d'estimation telles que le krigeage ou le cokrigeage pour resoudre leurs problemes. Nous nous proposons d'etudier les deux equations suivantes dans le contexte de la theorie des fonctions aleatoires intrinseques d'ordre k(fai-k): z=y (1) z=y/x (2). Designe l'operateur laplacien et x une direction de l'espace. L'equation (1) sera examinee sous l'angle du "probleme direct" (estimation de z a partir des donnees de y) tandis que l'equation (2) sera abordee sous l'angle du "probleme inverse" (estimation de y a partir des donnees de z). La premiere etape est consacree a la determination des couples de covariances de z et y compatibles avec l'equation etudiee. En interpretant l'equation (2) en termes de distributions, on pourra s'inspirer de l'analyse structurale de l'equation (1) et utiliser des modeles non differentiables pour y. Dans le cas du "probleme direct", on construira des estimateurs par cokrigeage de z qui: 1) tiennent compte des conditions aux limites du type dirichlet, neumann et surfaces libres et plus generalement des conditions du type lineaire; 2) verifient l'equation (1) en tout point ou y est connue. Contrairement aux methodes traditionnelles, le cokrigeage peut prendre en compte des conditions aux limites ne situant pas a la frontiere et n'exige pas la connaissance de la total

Ce rapport présente une méthode pour estimer localement, au niveau des récepteurs, les paramètres d'anisotropie à partir de données sismiques de puits enregistrées en trois composantes. On considère des milieux transversalement isotropes d'axe de symétrie contenu dans le plan vertical. Les données sont constituées par les lenteurs apparentes verticales et les directions de polarisation des ondes p et sv mesurées en profondeur dans un puits. Les paramètres d'anisotropie sont inversés à partir de ces données par des techniques d'optimisation

Des argiles kaoliniques sont exploitées depuis une centaine d'années, dans la partie occidentale du bassin des Charentes. Leurs dépôts sont liés aux formations sédimentaires d'un réseau paléo-deltaïque. La genèse des gisements de kaolin du Bassin des Charentes obéit à des règles géologiques simples mais leur géométrie détaillée présente une grande complexité, renforcée par le fait qu'on doit distinguer des argiles de qualités très différentes.
Les gisements récemment mis en exploitation, sont relativement plus riches en alumine que ceux déjà épuisés. Ces nouveaux gisements, plus profonds que les anciens, ont subi des phénomènes postsédimentaires plus intenses. Dans cette étude on compare des données chimiques des différents gisements de Bassin des Charentes. La formation des argiles hyperalumineuses est particulièrement étudiée. On essaie d'améliorer la connaissance des gisements en utilisant les analyses statistiques. Les données disponibles sont les résultats d'analyses chimiques des échantillons carottés des argiles kaoliniques.
Des méthodes géostatistiques sont proposées pour l'étude de ces gisements et sont appliquées à certains d'entre eux. Les problèmes de l'utilisation des outils de la géostatistique pour ce genre de dépôt sont abordés et on essaie de présenter des démarches adaptées.
Ces outils nous permettent de mettre en évidence l'importance d'une connaissance de la variabilité spatiale des éléments chimiques dans ce gisement pour mieux gérer les opérations d'exploration, d'extraction et de classification des dépôts de kaolin.
Des simulations géostatistiques permettent de visualiser la variabilité spatiale des différentes catégories d'argile. Elles ont été mises en oeuvre dans le cadre de la méthode des plurigaussiennes seuillées. Deux cas ont été traités : simulation de faciès argileux fondés principalement sur des critères visuels, et simulation de classes d'argile définies à partir de leur teneur en Al2O3.

L'aménagement de la pêcherie céphalopodière au Maroc est principalement axé sur la régulation des captures du poulpe (Octopus vulgaris), l'espèce prioritaire. Son potentiel exploitable est déterminé sur la base des évaluations directes à travers des campagnes de prospection scientifiques. Il s'agit de (i) fournir un estimateur de bonne qualité – sans biais et précis – de l'abondance globale, accompagné de sa variance d'estimation et des cartes de distribution du stock, (ii) suivre la dynamique spatio-temporelle du poulpe dont la connaissance est indispensable pour la gestion. Pour y parvenir : • Les propriétés spatiales de la densité du poulpe sont décrites et les conditions de stationnarité sont analysées sur une série de données historique. • Des méthodes d'estimation globale – transitives et intrinsèques – sont appliquées en tenant compte des conditions particulières des données, notamment la non-stationnarité concernant l'approche probabiliste intrinsèque et l'échantillonnage non contrôlé concernant l'approche déterministe transitive. • La stratégie d'occupation de l'espace par le poulpe est décrite à l'aide d'indicateurs géostatistiques. Elle se caractérise par un schéma de reproduction étendu sur le plateau continental et un succès du recrutement essentiellement côtier. D'importantes implications en matière de gestion sont discutées. Pour finir, ce travail est replacé dans son contexte plus global de l'aménagement de la pêcherie céphalopodière, et de la nouvelle stratégie de la pêche. De nouvelles perspectives pour l'application des techniques géostatistiques et des statistiques spatiales sont envisagées pour les évaluations de stocks, en particulier du poulpe
The regulation of the cephalopod fishery off Morocco is mainly focused on catch limitation of the target specie, octopus (Octopus vulgaris). Its exploitable biomass is estimated by direct assessments through scientific surveys. These have to (i) provide a good estimator - unbiased and accurate – of the global abundance with its estimation variance and stock distribution maps, (ii) monitor the spatio-temporal dynamic of the octopus, whose knowledge is essential for management. To achieve this: • The spatial features of the octopus densities are described and the conditions of stationarity are analyzed on the basis of historical data. • Methods for global estimation - transitive and intrinsic - are applied taking into account the data characteristics, including non-stationarity on the probabilistic intrinsic approach and the uncontrolled sampling on the deterministic transitive approach. • The spatial pattern of the octopus is described using geostatistical indicators. It is characterized by a widespread distribution of the spawning stock over the continental shelf and a coastal recruitment success. Some important implications for management are discussed. Finally, this work is seen within larger context of the cephalopod fishery management and the new fishing policy. New perspectives for the application of the geostatistical techniques and the spatial statistics are considered for stock assessments, especially of the octopus

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éostatiques, appliquant la théorie des fonctions aléatoires. Après avoir introduit des éléments de géomatique et de géostatique, et précisé la nature spécifique de l'autocorrélation spatiale, nous introduisons les inférences design-basedet model-based. Á l'aide de fonctions aléatoires, nous étudions l'efficacité de l'échantillonnage, optimisons les 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. 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.

Hors des zones urbanisées, la propagation à grande distance du bruit émis par des sources routières ou industrielles est fortement influencée par les propriétés du sol ainsi que par les conditions atmosphériques de la couche de surface. Il est donc admis que la caractérisation d'une situation de pollution sonore doit inclure en particulier l'information micrométéorologique. Or la relation entre conditions de propagation et niveaux sonores mesurés est statistiquement incertaine. L'objectif de cette thèse est d'apporter les premiers éléments de réflexion sur l'estimation temporelle et spatiale du champ acoustique et de ses paramètres d'influence grâce à une analyse objective par les méthodes géostatistiques. La première partie de nos travaux portent sur l'analyse de séries temporelles acoustiques finement échantillonnées. L'autocorrélation est observée grâce au variogramme expérimental dont on tente d'expliquer l'allure par les conditions micrométéorologiques correspondantes. La deuxième partie de nos travaux se focalise sur les propriétés acoustiques du sol. Le paramètre phénoménologique de résistance spécifique au passage de l'air est étudié temporellement et spatialement. Nous obtenons des interpolations spatiales optimisées par Krigeage ainsi que des simulations du champ de ce paramètre ouvrant sur des perspectives méthodologiques pour des modèles numériques évolués de propagation acoustique. En acceptant d'emblée le champ acoustique comme un champ stochastique grâce aux outils géostatistiques d'analyse de la structure spatiale des mesures de terrain, nous obtenons dans la troisième partie de nos travaux, une première méthode de cartographie acoustique
In non-urban zones, the long range sound propagation for road or industrial sources is highly dependent on the boundary conditions (ground characteristics) and on the stability conditions of the surface boundary layer of the atmosphere. Therefore the characterization of the sound propagation conditions are meant to include the micrometeorological information. Besides, the relation between propagation conditions and the acoustic level is statistically uncertain. The purpose of the thesis is to bring new ideas and methods on the problem of time and space estimation of the acoustic field an its related parameters of influence within the geostatistical framework. First, acoustic time series are analysed. The experimental variogram of the series gives an image of the autocorrelation, the shape of which we try to explain with the corresponding micrometeorological conditions. The second part of the work is focused on both time and space evolution of the ground properties of an average field. The ground specifie air flow resistivity is the main parameter of the study. An ordinary Kriging technique is used for mapping. Conditional simulations are also obtained and can be utilized in future studies as input data in complex numerical models. Third, when considering the acoustic field as a stochastic field, it is possible to analyse in situ data. We develop a method of space estimation of the sound level using geostatistical tools combining both a simple physical model and a random model of the fluctuations