Brevet FR3074340A1 — publié le 20190531

17-IPM-101756-U1-FR 1 ANALYSE DE LA PROVENANCE DE SÉDIMENTS GÉOLOGIQUES ET SYSTÈME D'AFFICHAGE Domaine technique La présente divulgation concerne d'une manière générale le traitement numérique de données géochronologiques obtenues à partir de matériaux géologiques. Plus précisément, la présente divulgation concerne le traitement d'informations provenant de sédiments au moyen de fonctions générées automatiquement et l'affichage interactif de ces informations sur une carte paléogéographique. Contexte Une analyse de la source au bassin est actuellement une évaluation effectuée au cours de plusieurs semaines ou plusieurs mois par des géoscientifiques qui cherchent à comprendre les risques associés aux réservoirs clastiques d'hydrocarbure. L'analyse est habituellement réalisée par la cartographie de l'origine des sédiments et de leurs réseaux de drainage historiques tout au long du temps géologique. L'analyse implique habituellement l'utilisation de nombreux ensembles de données disparates. Ces dernières années, le nombre d'études publiées portant sur la géochronologie des zircons détritiques a augmenté de manière exponentielle. Les ensembles de données de géochronologie détritique issus de ces études peuvent être des aides importantes pour une analyse de la source au bassin et la caractérisation des réservoirs clastiques. Cependant, comme la quantité de données disponibles augmente, les techniques traditionnelles d'analyse de données utilisées dans les études régionales deviennent intenables. Brève description des figures La figure 1 est un schéma de principe représentant un dispositif de calcul pour effectuer une analyse de provenance géologique incluant une analyse de la source au bassin selon certains aspects. La figure 2 est un graphique illustrant un stade de regroupement agglomératif selon certains aspects. La figure 3 est un graphique illustrant un autre stade de regroupement agglomératif selon certains aspects. La figure 4 est un graphique illustrant des regroupements agglomératifs selon certains aspects. 2017-IPM-101756-U1-FR 2 La figure 5 est un schéma de procédé illustrant un procédé d'analyse de la source au bassin et d'affichage de provenance selon certains aspects. La figure 6 représente un affichage graphique de cartes produites selon certains aspects. 5 La figure 7 représente un affichage graphique incluant une carte et un histogramme produits selon certains aspects. La figure 8 est un graphique illustrant l'application d'une espérance-maximisation à des modèles de mélange gaussien selon certains aspects. La figure 9A et la figure 9B représentent un schéma de procédé illustrant un 10 procédé d'analyse et d'affichage de provenance selon certains aspects. Le schéma de procédé est divisé en deux figures pour clarifier l'illustration. La figure 10 est un schéma de procédé illustrant un procédé d'analyse du bassin au bassin et d'affichage selon certains aspects. La figure 11 représente un affichage graphique incluant un dendrogramme et 15 des distributions d'âges selon certains aspects. Description détaillée Certains aspects de la présente divulgation concernent un système et des procédés qui améliorent, et rendent plus efficace, l'analyse et la comparaison de multiples 20 populations statistiques de données géochronologiques obtenues à partir de matériaux géologiques. Les procédés utilisent des méthodes statistiques et des techniques d'apprentissage par machine pour identifier et caractériser rapidement des sous-populations de données (par exemple des âges de grain unique de zircon détritique) parmi une population d'échantillon plus large ou un groupe de populations. Le système compare des populations distinctes sur la base de 25 ces sous-populations. Les termes « population »et « sous-population » sont utilisés dans le présent document dans le sens statistique. Une population est un vaste groupe de zircons qui peuvent être généralisés à partir d'un nombre limité d'échantillons par rapport à une certaine propriété. Une sous-population est une partie d'une population où les zircons partagent une propriété supplémentaire. Des exemples de propriétés comprennent une plage d'âges, un 30 emplacement géographique, et un temps de dépôt. Le système et les procédés décrits dans le présent document fournissent des études de provenance automatisées, intégrées et rapides, qui peuvent être utilisées dans l'exploration de ressources naturelles pour évaluer la nature et la qualité des roches sédimentaires accueillant les réservoirs d'hydrocarbure. 2017-IPM-101756-U1-FR 3 Les études fournies peuvent comprendre, par exemple, des études de la source au bassin ou du bassin au bassin. Les résultats peuvent être affichés de diverses manières. Par exemple, les résultats peuvent être affichés en utilisant une carte géographique du jour présent, en utilisant une carte paléographique basée sur un modèle géodynamique de tectonique des 5 plaques pour montrer les points de données et leur emplacement pour le laps de temps géologique approprié, ou simultanément sur les deux types de cartes. Comme autre exemple, les résultats peuvent être affichés simultanément sur une carte et avec un histogramme montrant les composantes de pic d'âge. Comme exemple supplémentaire, des comparaisons de distribution d'âges et des dendrogrammes pour une analyse de provenance peuvent être affichés. 10 L'analyse et l'affichage d'informations de la source au bassin selon certains aspects comprennent le regroupement de données géochronologiques associées à un échantillon de sédiment en des sous-populations optimisées au sein d'une population de référence et de populations cibles, et la production de fonctions gaussiennes pour la population de référence et les populations cibles en utilisant les sous-populations comme contraintes a priori. Les fonctions 15 gaussiennes décrivent une distribution de zircons. Les sous-populations au sein de la population de référence et des populations cibles sont comparées sur la base d'au moins un attribut statistique provenant des fonctions gaussiennes pour identifier des zones de provenance de sédiment, et les zones de provenance de sédiment sont affichées de manière géospatiale sur une carte paléographique sous la forme d'un âge de dépôt de l'échantillon de sédiment. 20 La figure 1 représente un exemple d'un dispositif de calcul 100 selon un exemple. Le dispositif de calcul 100 peut inclure un dispositif de traitement 102, un bus 104, une interface de communication 106, un dispositif de mémoire non transitoire ou non volatile 108, un dispositif d'entrée utilisateur 124, et un dispositif d'affichage 126. Dans certains exemples, certains ou l'ensemble des composants montrés dans la figure 1 peuvent être intégrés 25 dans une unique structure, telle qu'un unique boîtier. Dans d'autres exemples, certains ou l'ensemble des composants montrés dans la figure 1 peuvent être distribués (par exemple, dans des boîtiers séparés) et en communication électrique les uns avec les autres. Le dispositif de traitement 102 peut exécuter une ou plusieurs opérations pour fournir une analyse de la source au bassin. Le dispositif de traitement 102 peut exécuter des 30 instructions 111 stockées dans le dispositif de mémoire 108 qui sont exécutables par le dispositif de traitement 102 pour effectuer les opérations. Le dispositif de traitement 102 peut inclure un dispositif de traitement ou de multiples dispositifs de traitement. Des exemples non limitatifs du dispositif de traitement 102 incluent un circuit intégré prédiffusé programmable (« FPGA »), un circuit intégré à application spécifique (« ASIC »), un dispositif de micro-traitement, etc. 2017-IPM-101756-U1-FR 4 Le dispositif de traitement 102 peut être couplé de manière communicative au dispositif de mémoire 108 via le bus 104. Le dispositif de mémoire non volatile 108 peut inclure n'importe quel type de dispositif de mémoire qui conserve des informations lorsqu'il n'est pas alimenté en énergie. Des exemples non limitatifs du dispositif de mémoire 108 incluent une 5 mémoire morte effaçable et programmable électriquement (« EEPROM »), une mémoire flash, ou n'importe quel autre type de mémoire non volatile. Dans certains exemples, au moins une partie du dispositif de mémoire 108 peut inclure un support non transitoire à partir duquel le dispositif de traitement 102 peut lire des instructions. Un support non transitoire lisible par ordinateur peut inclure des dispositifs de stockage électroniques, optiques, magnétiques, ou 10 autres capables de fournir au dispositif de traitement 102 des instructions lisibles par ordinateur ou un autre code de programme. Des exemples non limitatifs d'un support lisible par ordinateur incluent (mais sans s'y limiter) un ou plusieurs disques magnétiques, une ou plusieurs puces mémoires, une mémoire morte (ROM), une mémoire vive (« RAM »), un ASIC, un dispositif de traitement configuré, un dispositif de stockage optique, ou un quelconque autre support à partir 15 duquel un dispositif de traitement informatique peut lire des instructions. Les instructions peuvent inclure des instructions de traitement spécifiques au dispositif générées par un programme de compilation ou un programme d'interprétation à partir d'un code écrit dans n'importe quel langage de programmation informatique adapté, y compris, par exemple, C, C++, C#, etc. 20 Dans certains exemples, le dispositif de mémoire 108 peut inclure des ensembles de données de référence 110, rassemblés et groupés à partir d'une base de données géochronologiques ou de multiples bases de données géochronologiques. Ces données peuvent représenter des données de zircons détritiques provenant d'un échantillon de réservoir ou des âges de minéraux provenant d'une région géographique spécifique. Dans certains exemples, les 25 instructions de programme informatique 111 définissent des sous-populations au sein de chaque groupe de données étudié, caractérisent les données par des attributs statistiques (par exemple, une moyenne, une variance, une proportion de population totale, une distribution asymétrique, un coefficient d'aplatissement, etc.), comparent et corrèlent les sous-populations, et effectuent les autres aspects du procédé. Le dispositif de mémoire 108 peut également inclure des 30 ensembles de données cibles 112, rassemblés et groupés à partir d'une ou de plusieurs bases de données géochronologiques. Le dispositif de mémoire 108 peut inclure des données statistiques de population 113, qui peuvent inclure des attributs statistiques de sous-populations provenant de données géochronologiques. Le dispositif de mémoire 108 peut comprendre une carte paléogéographique 114 stockée et une carte géographique stockée 122. 2017-IPM-101756-U1-FR 5 Dans certains exemples, le dispositif de calcul 100 inclut une interface de communication 106. L'interface de communication 106 peut représenter un ou plusieurs composants qui facilitent une connexion réseau ou autrement facilitent une communication entre des dispositifs électroniques. Des exemples incluent, mais sans s'y limiter, des interfaces 5 câblées telles que des interfaces Ethernet, USB, IEEE 1394, et/ou sans fil telles que des interfaces IEEE 802.11, Bluetooth, de communication en champ proche (CCP), des interfaces RFID, ou des interfaces radio pour un accès à des réseaux téléphoniques cellulaires (par ex., un émetteur-récepteur/une antenne pour un accès à un réseau AMRC, GSM, UMTS, ou un autre réseau de communication mobile). 10 Dans certains exemples, le dispositif de calcul 100 inclut un dispositif d'entrée d'utilisateur 124. Le dispositif d'entrée d'utilisateur 124 peut représenter un ou plusieurs composants utilisés pour entrer des données. Des exemples du dispositif d'entrée d'utilisateur 124 peuvent inclure un clavier, une souris, un pavé tactile, un bouton, ou un affichage à écran tactile, etc. 15 Dans certains exemples, le dispositif de calcul 100 inclut un dispositif d'affichage 126, qui peut afficher les cartes décrites dans le présent document et les autres informations utilisées dans le procédé décrit dans le présent document. Le dispositif d'affichage 126 peut représenter un ou plusieurs composants utilisés pour délivrer des données en sortie. Des exemples du dispositif d'affichage 126 peuvent inclure un affichage à cristaux liquides 20 (LCD), une télévision, un moniteur d'ordinateur, un affichage à écran tactile, etc. Dans certains exemples, le dispositif d'entrée d'utilisateur 124 et le dispositif d'affichage 126 peuvent être un unique dispositif, tel qu'un affichage à écran tactile. Les exemples de techniques décrites ci-dessous utilisent de grands ensembles de données globalement cohérents de données géochronologiques, à la fois pour des zones 25 d'arrière-pays et de bassin, qui sont continuellement sourcés, actualisés et conservés pour permettre des prédictions régionales améliorées. Les grands ensembles de données sont disponibles sous un format standardisé qui permet une comparaison directe, des calculs statistiques et des résultats rapides. Les âges des zircons détritiques pour des échantillons de sédiment publiés dans le domaine public sont extraits et entrés dans une base de données 30 géochronologique relationnelle. Des données d'âge à la fois pour des régions d'arrière-pays et de bassin sont incluses. Des champs contenant l'âge et une erreur analytique de chaque cristal de zircon sont utilisés, comme l'est l'âge de dépôt de l'unité sédimentaire d'où le cristal de zircon a été extrait. Une fois que ces informations sont capturées et standardisées, elles peuvent être importées à des fins d'analyse statistique. 2017-IPM-101756-U1-FR 6 Des regroupements agglomératifs, par exemple, peuvent être utilisés pour regrouper des sous-populations initiales dans les données d'âge de zircon. Tout d'abord, un k petit (nombre de regroupements) est utilisé pour identifier des sous-populations vastes. Les graphiques de la figure 2, de la figure 3 et de la figure 4 illustrent cette technique. La figure 2 5 représente un graphique 200 où des regroupements de sous-population sont présentés en fonction de l'âge pour k=16. Sur un affichage réel, le nombre de sous-populations apparentes pourrait être indiqué avec des couleurs différentes. À des fins d'illustration, la figure 2 nous permet d'observer, d'une manière générale, différents regroupements par la variation de l'ombrage ou du remplissage. Comme le montre l'ombrage unitaire des grands regroupements 202 en dessous 10 d'environ 200 Ma d'âge, davantage de regroupements sont nécessaires pour les sous-populations plus jeunes. Pour les regroupements de sous-populations présentés sur le graphique 300 de la figure 3, k est augmenté et les regroupements sont évalués en utilisant le score de Shapiro. La solution avec le score de Shapiro le plus grand est la solution la plus optimisée et est utilisée pour fixer une valeur pour k. Dans le cas de la figure 3, k=29. Ce paramètre donne une division 15 de regroupements improbable 302 et des sur-interprétations 304, 306 et 308 du nombre de regroupements. Des sous-regroupements supplémentaires sont définis sur la base d'un ensemble de paramètres définis par l'utilisateur. Ces paramètres sont des règles géologiques avec des valeurs qui spécifient la plage d'âges maximale des événements géologiques permis 20 dans un regroupement particulier. Tous les regroupements avec des plages d'âges supérieures à ces valeurs de règle géologique sont décomposés par le procédé de sous-regroupements. Le procédé de sous-regroupement implique l'utilisation d'un regroupement agglomératif et d'un regroupement en k-moyens pour optimiser les sous-regroupements. Une partie finale du procédé de regroupement agglomératif prend en compte l'erreur analytique de chacun des âges des 25 zircons, joint les regroupements qui se trouvent dans l'erreur analytique chaque autre et divise les regroupements qui se trouvent hors de l'erreur analytique de chaque autre. La figure 4 représente un graphique 400 avec des sous-regroupements optimisés en utilisant cette approche de regroupement à deux étapes. Un nombre décroissant de regroupements est attribué à mesure que l'âge augmente de gauche à droite, comme indiqué par la flèche 402. 30 La figure 5 est un schéma de procédé montrant un procédé 500 pour mettre en oeuvre une analyse de la source au bassin et afficher des résultats incluant des informations d'analyse de la source au bassin pour un échantillon de sédiment. Dans cet exemple, le procédé 500 est réalisé par le dispositif de calcul 100 de la figure 1. Au bloc 502, on a accès à des archives numériques, par exemple par l'intermédiaire d'une interface de communication 106 de 2017-IPM-101756-U1-FR 7 la figure 1. Les archives numériques comprennent des données qui peuvent représenter des données de zircons détritiques provenant d'un échantillon de sédiment de réservoir ou des âges de minéraux provenant d'une région géographique spécifique. Les données sélectionnées peuvent également se limiter à une plage d'âges géologiques spécifique. Au bloc 504, les 5 données géochronologiques cibles et les données géochronologiques de référence sont regroupées dans des groupes de population distincts représentant une population de référence et des populations cibles. Les données sont stockées sous la forme d'un ensemble de données de référence 110 et d'ensembles de données cibles 112. Au bloc 506, des sous-populations au sein de la population de référence et des populations cibles sont identifiées et définies par le dispositif 10 de traitement 102. Dans certains aspects, un regroupement agglomératif de données d'âge de zircon, comme précédemment décrit, est utilisé pour définir les sous-populations. Au bloc 508, les sous-populations sont analysées de manière statistique par le dispositif de traitement 102 et caractérisées pour produire un attribut statistiques ou des attributs statistiques sur la base de la distribution des zircons dans l'échantillon. Ces attributs statistiques sont stockés sous la forme 15 des données de populations statistiques 113 représentées sur la figure 1. Les attributs statistiques peuvent comprendre, par exemple, une moyenne, une variance, une proportion de population totale, une distribution asymétrique et un coefficient d'aplatissement. Dans certains aspects, l'analyse statistique comprend la définition de modèles de mélange gaussien qui reproduisent la population de référence et les populations cibles en utilisant les sous-populations comme 20 contraintes a priori. Les modèles de mélanges gaussien décrivent une distribution de zircons dans l'échantillon de sédiment. Les détails d'un exemple d'analyse statistique seront davantage discutés ci-dessous par rapport aux figures 8, 9A et 9B. En continuant avec la figure 5, au bloc 510 le dispositif de traitement 102 compare les sous-populations en utilisant l'attribut ou les attributs statistiques pour déterminer 25 des corrélations entre la population de référence et les populations cibles. De telles corrélations peuvent indiquer une provenance de sédiment, puisqu'elles déterminent la probabilité que chaque sous-population provenant d'un échantillon de référence soit présente dans les échantillons cibles. Au bloc 512, les corrélations sont affichées sur le dispositif d'affichage 126 de manière géospatiale sous la forme de zones de provenance de sédiment montrant les liens 30 entre les régions de source et de bassin. De nombreux types et options d'affichage peuvent être sélectionnés et contrôlés par un utilisateur grâce au dispositif d'entrée utilisateur 124. Certaines de ces options seront discutées avec de plus amples détails par rapport aux figures restantes. Par exemple, un utilisateur peut choisir au bloc 514 que le dispositif de calcul 100 affiche sélectivement un histogramme de distribution d'âge géologique par rapport à une fréquence 2017-IPM-101756-U1-FR 8 d'apparition. Comme autre exemple, les corrélations peuvent être affichées simultanément sur une carte géographique du jour présent et sur une carte paléogéographique à un temps approprié. La carte paléographique peut être basée sur un modèle géodynamique de tectonique des plaques. Un utilisateur peut également choisir au bloc 516 que le dispositif de calcul 100 effectue une 5 analyse du bassin au bassin. Si une telle analyse est déclenchée, le dispositif de traitement 102 peut, au bloc 518, calculer et identifier des valeurs de dissimilarité représentant les différences entre l'échantillon de sédiment et des échantillons supplémentaires qui peuvent indiquer des changements de provenance. Dans certains aspects, un dendrogramme indiquant les valeurs de dissimilarité peut être affiché. Cette caractéristique et le procédé du bloc 518 seront discutés 10 avec de plus amples détails par rapport aux figures 10 et 11. La figure 6 représente un affichage graphique 600 de cartes produites selon certains aspects. Les cartes sont présentées à une échelle continentale (centaines à milliers de km) sur le dispositif d'affichage 126. La zone cartographique 602 comprend une carte géographique du jour présent avec des liens entre des régions de source et des régions de bassin, 15 et la zone cartographique 604 comprend une carte paléogéographique à un temps approprié. La zone cartographique 602 peut être stockée sous la forme de la carte géographique 122 de la figure 1. La zone cartographique 602 présente un ensemble de données géochronologiques de référence 606, basé sur, par exemple, un échantillon de sédiment d'un intervalle de réservoir contenant des informations géochronologiques de zircon détritique ou d'un groupe 20 d'échantillons géochronologiques détritiques provenant de la même formation ou du même intervalle de temps. La zone cartographique 602 présente également des ensembles de données géochronologiques cibles 608, 610 et 612, basés sur, par exemple, des données géochronologiques de roche dure, groupés par régions distinctes, relatifs à la formation ou à une modification de cette région d'arrière-pays. Des sous-populations, d'âge géologique, au sein de 25 chaque groupement sont calculées et peuvent être affichées sous la forme d'histogrammes d'âge par rapport à la fréquence, qui sont schématiquement illustrés sur la figure 6 pour chaque ensemble de données. Les sous-populations sont ensuite statistiquement comparées et un score de probabilité est calculé, par exemple, pour indiquer des zones de source probables qui ont été érodées pour produire les sédiments étudiés, par exemple, la région correspondant à l'ensemble 30 de données 606. Chaque ensemble de données de région est présenté avec une ligne indiquant la région correspondante sur la carte 614. Les flèches indiquent des liens montrant des voies potentielles pour relier les régions de source aux régions de bassin, par exemple, des systèmes de rivières. 2017-IPM-101756-U1-FR 9 En se référant encore à la figure 6, les résultats sont également visualisés dans un contexte paléographique dans une zone cartographique 604 pour évaluer la configuration continentale de tectonique des plaques au moment de la formation de la région correspondant à l'ensemble de données 606, par exemple il y a 200 millions d'années. Dans ce cas, la carte 5 paléographique, stockée sous la forme de la carte paléographique 114 sur la figure 1, est représentée en deux parties, la partie 620 et la partie 622. L'affichage des régions d'intérêt à un temps de référence fournit un aperçu critique supplémentaire et peut être utilisé pour déterminer que certaines zones cibles n'ont pas de connexion avec la zone de référence, par exemple la zone correspondant à l'ensemble de données 612 ne peut pas fournir de matériau à la zone 10 correspondant à l'ensemble de données 606 dans cette fenêtre temporelle géologique particulière, comme le montre la partie cartographique 622, et doit être éliminée pour les études de la source au bassin. Les parties cartographiques paléographiques dans la zone cartographique 604 montrent que deux liens peuvent être écartés. La figure 7 représente un affichage graphique interactif 700 incluant une carte 15 702 et un histogramme 704 produits selon certains aspects et affichés sur le dispositif d'affichage 126. Dans cet exemple, les résultats géospatiaux sont reconstruits sur un modèle de plaque géodynamique à l'âge du dépôt d'un échantillon de sédiment détritique et affichés sur la carte 702. Cet affichage fournit à l'utilisateur une visualisation de la paléogéographie au moment du dépôt et permet à l'utilisateur d'écarter les zones de source de sédiment qui peuvent ne pas avoir 20 été reliées au bassin au moment du dépôt. Un tel affichage peut mettre en évidence des zones de source de sédiment précédemment inconnues qui ne sont plus du tout connectées géographiquement au bassin au jour présent. L'utilisateur peut afficher sélectivement un histogramme montrant des composantes de pic d'âge pour un échantillon spécifique. Le terme « sélectivement » dans ce contexte peut avoir une ou deux significations, que la sélection de 25 l'utilisateur détermine si histogramme doit être affiché ou non et également que la sélection de l'utilisateur peut déterminer quel histogramme d'échantillon est affiché. Dans un cas ou les deux cas, un affichage interactif de données concernant des échantillons géologiques visualisés par un utilisateur a lieu, le dispositif de calcul répondant aux sélections de l'utilisateur. Dans cet exemple, l'utilisateur peut cliquer, toucher, ou « passer une souris » sur un point d'échantillon 30 particulier sur la carte avec le dispositif d'entrée utilisateur 124 et l'histogramme de distribution d'âges pour cet échantillon particulier apparaît. Dans le cas particulier de la figure 7, l'histogramme 704 correspond à l'échantillon 708. Comme autre exemple de l'interactivité de l'affichage 700, un utilisateur peut cliquer, toucher, ou « passer une souris » sur une icône ou un bouton pour amener une palette de 2017-IPM-101756-U1-FR 10 couleur à être affichée. Un affichage visuel réel comme celui représenté sur la figure 7 pourrait inclure une couleur pour mettre en évidence des échantillons et pour mettre en évidence des parties de l'histogramme. Une couleur pourrait également être utilisée pour la superposition de la courbe de distribution illustrée sur l'histogramme 704. Pour la clarté de l'illustration sur la 5 figure 7, un ombrage ou remplissage varié est utilisé à des fins d'illustration. La figure 8 est un graphique illustrant l'application d'une espérance-maximisation à des modèles de mélange gaussien selon certains aspects. Ces fonctions gaussiennes définies par l'utilisation de modèles de mélange gaussien peuvent être utilisées pour l'analyse statistique présentée au bloc 508 de la figure 5. Le graphique 800 présente un certain 10 nombre de sous-populations, qui sont utilisées comme entrée pour la modélisation de mélange gaussien afin de définir des fonctions gaussiennes pour chaque sous-population. Les fonctions gaussiennes peuvent être utilisées pour décrire la distribution statistique de zircons dans chaque échantillon de sédiment. Dans certains aspects, les sous-populations sont utilisées comme contraintes a priori sur un algorithme d'espérance-maximisation (EM) pour générer des modèles 15 de mélange pour chacune des sous-populations. Sur la figure 8, un algorithme d'espérance- maximisation est exécuté pendant trois étapes appliquées aux modèles de mélange gaussien sur des sous-populations de zircon prédéfinies. La fonction gaussien initiale 802 est superposée sur l'histogramme des sous-populations. Grâce à l'algorithme EM-génération de modèles de mélange gaussien, le système utilise un apprentissage par machine pour établir ce qui est 20 nécessaire pour connaître l'échantillon de l'instant afin de produire automatiquement une ou plusieurs fonctions gaussiennes à évaluer afin d'analyser l'échantillon. Des fonctions issues d'une maximisation EM avec diverses probabilités de différence sont présentées sous la forme des lignes 804. Un affichage visuel réel comme celui représenté sur la figure 8 pourrait inclure une couleur pour mettre en évidence des parties de l'histogramme et pour distinguer les 25 fonctions de maximisation EM. Pour la clarté de l'illustration sur la figure 7, un ombrage ou remplissage varié est fourni pour la partie d'histogramme. Les figures 9A et 9B représentent un schéma de procédé illustrant un procédé 900 d'analyse et d'affichage de provenance selon certains aspects. Le schéma de procédé est divisé en deux figures pour clarifier l'illustration. Le schéma de procédé est divisé en outre en 30 trois stades de traitement et ces stades sont indiqués avec des légendes appropriées. Chaque stade est divisé en outre en sections. Le procédé 900 commence au bloc 902 de la figure 9A, où des âges de zircon détritique pour des échantillons de sédiment fournis dans des sources publiquement disponibles sont extraits et importés. Un regroupement agglomératif, comme précédemment décrit par rapport aux figures 1 à 3 et au bloc 506 de la figure 4 a lieu dans le 2017-IPM-101756-U1-FR 11 stade 1 du procédé 900. Des liens basés sur la continuité de Ward i sont utilisés pour chaque regroupement c, nombre j(c), où if max( 6-j): 30 REVENDICATIONS1. Dispositif de calcul comprenant : un dispositif de traitement ; et un dispositif de mémoire non transitoire comprenant des instructions qui sont exécutables par le dispositif de traitement pour amener le dispositif de traitement à effectuer des opérations comprenant : le regroupement de données géochronologiques associées à un échantillon de sédiment en des sous-populations optimisées au sein d'une population de référence et de populations cibles ; la production d'une pluralité de fonctions gaussiennes pour la population de référence et les populations cibles en utilisant les sous-populations comme contraintes a priori, dans lequel la pluralité de fonctions gaussiennes décrit une distribution de zircons dans l'échantillon de sédiment ; la comparaison des sous-populations au sein de la population de référence et des populations cibles sur la base d'au moins un attribut statistique provenant de la pluralité de fonctions gaussiennes pour identifier des zones de provenance de sédiment ; et l'affichage géospatial des zones de provenance de sédiment sur une carte paléographique sous la forme d'un âge de dépôt de l'échantillon de sédiment. 2. Dispositif de calcul selon la revendication 1 dans lequel la carte paléographique est basée au moins en partie sur un modèle géodynamique de tectonique des plaques. 3. Dispositif de calcul selon la revendication 1 dans lequel les opérations comprennent en outre : la réception d'une entrée comprenant une valeur de règle géologique ; et la définition de sous-regroupements au sein de la sous-population en utilisant la valeur de règle géologique. 2017-IPM-101756-U1-FR 19 4. Dispositif de calcul selon l'une quelconque des revendications 1 à 3 dans lequel les opérations comprennent en outre la génération de modèles de mélange gaussien pour les sous-populations dans lequel la pluralité de fonctions gaussiennes est produite à partir des modèles de mélange gaussien. 5 5. Dispositif de calcul selon l'une quelconque des revendications 1 à 3 dans lequel les opérations comprennent en outre le fait d'afficher simultanément un histogramme de distributions de sous-populations de zircon avec des sources associées. 10 6. Dispositif de calcul selon la revendication 1 à 3 dans lequel les opérations comprennent en outre : le calcul de valeurs de dissimilarité représentant les différences entre l'échantillon de sédiments et des échantillons de sédiment supplémentaires ; et l'affichage d'un dendrogramme indiquant les valeurs de 15 dissimilarité. 7. Dispositif de calcul selon la revendication 6 dans lequel les opérations comprennent en outre l'identification d'au moins l'une parmi une relation spatiale ou une relation temporelle entre l'échantillon de sédiments et au moins l'un des échantillons de 20 sédiment supplémentaires. 8. Procédé comprenant : le regroupement, au moyen d'un dispositif de traitement, de données géochronologiques associées à un échantillon de sédiment en des sous-populations 25 optimisées au sein d'une population de référence et de populations cibles ; la production, au moyen du dispositif de traitement exécutant des modèles de mélange gaussien, d'une pluralité de fonctions gaussiennes pour la population de référence et les populations cibles en utilisant les sous-populations comme contraintes a priori, dans lequel la pluralité de fonctions gaussiennes décrit une 30 distribution de zircons dans l'échantillon de sédiment ; la comparaison, au moyen du dispositif de traitement, des sous-populations au sein de la population de référence et des populations cibles sur la base d'au moins un attribut statistique provenant de la pluralité de fonctions gaussiennes pour identifier des zones de provenance de sédiment ; et 2017-IPM-101756-U1-FR 20 l'affichage géospatial, sur un dispositif d'affichage, des zones de provenance de sédiment sur une carte paléographique sous la forme d'un âge de dépôt de l'échantillon de sédiment. 5 9. Procédé selon la revendication 8 dans lequel la carte paléographique est basée au moins en partie sur un modèle géodynamique de tectonique des plaques. 10. Procédé selon la revendication 8 comprenant en outre : la réception d'une entrée comprenant une valeur de règle 10 géologique ; et la définition de sous-regroupements au sein de la sous-population en utilisant la valeur de règle géologique. 11. Procédé selon l'une quelconque des revendications 8 à 10 comprenant en 15 outre le fait d'afficher simultanément un histogramme de distributions de sous-populations de zircon avec des sources associées. 12. Procédé selon l'une quelconque des revendications 8 à 10 comprenant en outre : 20 le calcul de valeurs de dissimilarité représentant les différences entre l'échantillon de sédiments et des échantillons de sédiment supplémentaires ; et l'affichage d'un dendrogramme indiquant les valeurs de dissimilarité. 25 13. Procédé selon la revendication 12 comprenant en outre l'identification d'au moins l'une parmi une relation spatiale ou une relation temporelle entre l'échantillon de sédiments et au moins l'un des échantillons de sédiment supplémentaires. 14. Procédé selon l'une quelconque des revendications 8 à 10 comprenant en 30 outre la génération de modèles de mélange gaussien pour les sous-populations dans lequel la pluralité de fonctions gaussiennes est produite à partir des modèles de mélange gaussien. 15. Support non transitoire lisible par ordinateur qui comprend des instructions qui sont exécutables par un dispositif de traitement pour amener le dispositif de 2017-IPM-101756-U1-FR 21 traitement à afficher des informations d'analyse de provenance de sédiment par la mise en oeuvre du procédé selon l'une quelconque des revendications 8 à 10. 5