Brevet FR2981169A1 — publié le 20130412

Dispositif et procédé d'identification de mécanisme source Contexte de l'invention Domaine technique Les modes de réalisation de l'objet décrit ici concernent généralement des procédés et des systèmes et, plus particulièrement, des mécanismes et des techniques d'identification d'un mécanisme source associé à un évènement microsismique souterrain. Description du contexte Les opérations de fracture simulées sont destinées à augmenter la productivité d'un puits de réservoir à 15 hydrocarbures. Ces opérations comprennent l'injection d'un fluide à haute pression dans une couche du sous-sol dans lequel se trouve le réservoir. L'injection du fluide produit des micro-fractures au sein de la couche. Cette technique permet d'accroître la 20 perméabilité du réservoir en favorisant la circulation des hydrocarbures, via les micro-fractures, vers le puits. Les micro-fractures des couches sont une source possible d'évènements microsismiques. D'autres 25 évènements microsismiques existent et doivent être surveillés pour différentes raisons, comme la gestion d'un puits. Ces évènements microsismiques peuvent comprendre les effets sismiques générés dans une formation non seulement par fracture, mais également 30 par déplétion, inondation, traitement, mouvement de faille, effondrement, brèche, compaction, ou tout(e) autre intervention ou effet souterrain(e) similaire. Pour déterminer ces évènements microsismiques, les géologues ou les physiciens enregistrent généralement, 35 à la surface de la Terre, les ondes générées par ces évènements. Les ondes sont enregistrées en fonction du temps par un ou plusieurs récepteur(s). Les capteurs sismiques sont déployés à proximité du réservoir. La méthode conventionnelle pour déterminer les évènements microsismiques dans la zone explorée consiste à déployer quelques capteurs à l'intérieur du/des puits d'observation et à surveiller les évènements de fracture. Il est également possible de déployer des capteurs sismiques au niveau ou près de la surface de la Terre. Cependant, les procédés connus ne permettent pas de traiter rapidement les données, et ne conviennent pas pour la surveillance du sous-sol en temps réel, comme cela est décrit ci-après. Un procédé connu applique un algorithme de semblance aux données enregistrées par un groupe de récepteurs en surface. Cet algorithme teste toutes les positions possibles de la source dans un volume donné au sein de la sous-surface. Pour chaque position possible de la source, l'algorithme calcule les retards de propagation correspondants des ondes directes entre la position source étudiée et chaque récepteur du groupe. Ensuite, les données brutes enregistrées par le groupe de récepteurs (voir figure 1) sont retardées dans le temps sur la base du retard de propagation calculé de chaque récepteur. La figure 1 illustre, pour un exemple synthétique, des traces enregistrées 2 à différents moments (axe Y) le long de différentes positions (axe X). Le résultat de ce processus est illustré sur la figure 2, sur laquelle les traces 2 sont alignées sur une ligne droite. Ainsi, une sommation de ces traces est constructive et détermine un pic 4, comme cela est illustré sur la figure 3. On remarque que l'alignement en ligne droite des 35 traces 2 sur la figure 2 a lieu uniquement lorsque la position sélectionnée pour la source correspond à la position réelle de la source. Pour les autres cas, c'est-à-dire lorsque la position sélectionnée ne correspond pas à la position réelle de la source, une répartition de traces différente est obtenue. Une fonction de semblance est utilisée pour calculer la position sélectionnée qui correspond à la position réelle de la source. Si l'on considère que ri est un échantillon de la dème trace à un moment donné et est corrigé pour une position source donnée, un facteur de numération Semb Nom de la fonction de semblance est défini par : Nr Semb Nom= E r. Nr étant le nombre total de traces enregistrées par le groupe de récepteurs. Un facteur de dénomination de la fonction de semblance est considéré comme étant : Nr Semb _Denom= Eri2 i=1 et la fonction de semblance est définie comme étant le rapport entre le carré du facteur de nomination et du facteur de dénomination : Semb = Semb Nom2 Semb _ Denom La fonction de semblance est calculée de manière 30 répétée pour différentes positions sélectionnées de la source, et la valeur la plus élevée de cette fonction 20 25 - 4 indique la position source réelle. Cependant, cette technique a ses propres limites. Par exemple, cette technique rencontre un phénomène physique appelé « mécanisme source ». Le mécanisme source (ou mécanisme focal) d'un évènement sismique décrit la déformation non élastique de la région source qui génère les ondes sismiques. Les mécanismes source sont dérivés d'une solution du tenseur de moments pour l'évènement sismique, qui est lui-même estimé par une analyse des formes d'ondes sismiques observées. Selon le mécanisme source, les déplacements du sol mesurés par les récepteurs peuvent être soit uniquement positifs (en cas d'explosion dans la sous-surface), soit positifs et négatifs (en cas de tremblement de terre), ou peuvent avoir d'autres répartitions. En d'autres termes, les ondes produites par la source sismique peuvent avoir, dans deux directions différentes, des polarisations différentes. Par exemple, la polarisation des ondelettes peut inversée, comme cela est illustré sur la figure figure 4 illustre un changement de polarisation la trace 6 au niveau de la position 17 le long de être 4. La entre l'axe X et la trace 8 au niveau de la position 18 le long de l'axe X. dans ce cas, la sommation impliquée par la fonction de semblance comprend les ondelettes caractérisées par des polarisations positives et négatives. Ainsi, après avoir retardé les traces, comme cela est illustré sur la figure 5, la sommation est partiellement destructive, si bien que l'amplitude de crête diminue, comme cela est illustré sur la figure 6. Plus spécifiquement, la crête 4 passe d'environ 0,27 (voir figure 3) à 0,12 (voir figure 6). Dans cet exemple, la réduction de la crête 4 représente une perte d'énergie de plus de 7 dB, ce qui est élevé pour des évènements microsismiques. Cette perte devient plus prononcée pour un cas réel qui est affecté par un bruit, comme cela est illustré sur la figure 7. La figure 7 illustre les traces 2 affectées par un bruit 9. Après avoir corrigé le retard temporel, comme cela est illustré sur la figure 8, une crête est déterminée sur la figure 9 comme n'étant pas correctement définie. En d'autres termes, ladite crête est difficile à déterminer car elle ne possède pas de forme correctement définie. Le pire cas se produit lorsque la polarisation est telle qu'il existe un même nombre de traces affectées par la polarisation négative que le nombre de traces affectées par la polarisation positive. Dans ce cas, la valeur de crête théorique après la sommation est sensiblement nulle. Dans une application, chaque signal d'évènement est pondéré par un coefficient lié à la forme du mécanisme source. Ce phénomène et l'inversion de signe d'ondelette résultante restent un problème irrésolu dans le domaine de la détection de source. Différentes techniques ont été développées pour résoudre ce problème. Avec certaines solutions, l'algorithme de semblance est modifié afin que seule la valeur absolue des échantillons soit considérée. Cette solution évite d'additionner les ondelettes à signes opposés. Cependant, étant donné que les échantillons de bruit sont toujours positifs, le niveau de bruit est augmenté, ce qui n'est pas souhaitable. Le rapport signal/bruit (SNR) d'un récepteur du groupe de récepteurs (c'est-à-dire l'énergie maximale comparée à l'énergie sonore) diminue sensiblement avec cette solution. Avec une technique améliorée, on peut considérer que la polarisation possède une cohérence locale. Ainsi, les traces sont réparties en N sous-blocs locaux dans lesquels la polarisation est globalement constante. La somme dans chaque sous-bloc devient globalement constructive. Ensuite, le résultat final est calculé à l'aide de la somme de la valeur absolue de tous les sous-blocs. Dans ce cas, le niveau de bruit est toujours élevé, étant donné que N valeurs positives de bruit sont additionnées, mais il est moins élevé par rapport à la solution précédente. Cette solution améliorée ne convient pas en cas d'enregistrements très bruyants, et crée une perte de SNR de 3 dB (pour un ensemble en 2 dimensions) ou 6 dB (pour un ensemble en 3 dimensions). Une autre solution consiste à calculer la semblance pour chaque position possible des inversions de signe. Cela donne un résultat final précis, mais le calcul de semblance nécessite un algorithme coûteux. Par exemple, pour la surveillance d'un réservoir, la taille du groupe de récepteurs est importante et le nombre de traces est élevé (plusieurs milliers). Ainsi, le coût de calcul devient prohibitif. Par conséquent, il serait souhaitable de prévoir 25 des systèmes et des procédés qui empêchent les problèmes et les inconvénients décrits ci-dessus. Résumé Selon un exemple de mode de réalisation, il existe 30 un procédé d'identification d'un mécanisme source associé à un évènement microsismique qui a lieu dans une sous-surface de la Terre. Le procédé comprend une étape de réception de traces associées à une pluralité de récepteurs situés dans un réseau au-dessus ou au- 35 dessous de la sous-surface ; une étape de division du réseau en un ensemble de sous-blocs, chaque sous-bloc comprenant un ou plusieurs récepteur(s) de la pluralité de récepteurs ; une étape de détermination, pour chaque sous-bloc, d'un élément de sommation correspondant sur la base des traces associées au sous-bloc correspondant ; une étape d'identification d'un nombre prédéterminé de modèles qui sont associés à plusieurs mécanismes source ; une étape de calcul, dans le dispositif de calcul, pour chaque modèle, d'une fonction de semblance dédiée auxdits modèles sur la base des éléments de sommation des sous-blocs ; et une étape d'identification du mécanisme source de l'évènement microsismique sur la base de la fonction de semblance. Un modèle comprend un premier sous-ensemble de sous-blocs à polarité positive et un second sous-ensemble de sous-blocs à polarité négative, le premier sous-ensemble et le second sous-ensemble formant l'ensemble de sous-blocs entier. Selon un autre exemple de procédé, il existe un procédé mis en place dans un dispositif de calcul afin d'identifier un mécanisme source associé à un évènement microsismique qui a lieu dans une sous-surface de la Terre. Le procédé comprend une étape de réception de traces associées à une pluralité de récepteurs situés dans un réseau au-dessus de la sous-surface ; une étape de division du réseau en un ensemble de sous-blocs, chaque sous-bloc comprenant un ou plusieurs récepteur(s) de la pluralité de récepteurs ; une étape de détermination, pour chaque sous-bloc, d'un élément de sommation correspondant sur la base des traces associées au sous-bloc correspondant ; une étape d'identification d'un nombre prédéterminé de modèles qui sont associés à plusieurs mécanismes source ; une étape de multiplication de chaque élément de sommation par un coefficient de polarisation selon un modèle correspondant ; une étape de calcul, dans le dispositif de calcul, pour chaque modèle, d'une fonction de semblance sur la base des éléments de sommation des sous-blocs ; et une étape d'identification du mécanisme 5 source de l'évènement microsismique sur la base de la fonction de semblance. Un modèle comprend un premier sous-ensemble des sous-blocs à polarité positive et un second sous-ensemble des sous-blocs à polarité négative, le premier sous-ensemble et le second sous-10 ensemble formant l'ensemble de sous-blocs entier. Selon un autre exemple de mode de réalisation, il existe un dispositif de calcul configuré pour identifier un mécanisme source associé à un évènement microsismique qui a lieu dans une sous-surface de la 15 Terre. Le dispositif de calcul comprend un processeur configuré pour recevoir des traces associées à une pluralité de récepteurs situés sur un réseau au-dessus de la sous-surface ; diviser le réseau en un ensemble de sous-blocs, chaque sous-bloc comprenant un ou 20 plusieurs récepteur(s) de la pluralité de récepteurs ; déterminer, pour chaque sous-bloc, un élément de sommation correspondant sur la base des traces associées au sous-bloc correspondant ; recevoir un nombre prédéterminé de modèles qui sont associés à 25 plusieurs mécanismes source ; multiplier chaque élément de sommation par un coefficient de polarisation selon un modèle prédéterminé ; calculer, pour chaque modèle, une fonction de semblance sur la base des éléments de sommation des sous-blocs ; et identifier le mécanisme 30 source de l'évènement microsismique sur la base de la fonction de semblance. Un modèle comprend un premier sous-ensemble de sous-blocs à polarité positive et un second sous-ensemble de sous-blocs à polarité négative, le premier sous-ensemble et le second sous-ensemble 35 formant l'ensemble de sous-blocs entier. Brève description des dessins Les dessins joints, qui sont intégrés à et font partie de la demande, illustrent un ou plusieurs mode(s) de réalisation et, avec la description, expliquent lesdits modes de réalisation. Sur les dessins : La figure 1 est un graphique qui illustre des données sismiques brutes collectées à partir de 10 plusieurs récepteurs situés à la surface de la Terre ; La figure 2 est un graphique qui illustre des données sismiques corrigées dans le temps, qui correspondent à la figure 1 ; La figure 3 est un graphique qui illustre une 15 trace associée à un évènement sismique après sommation ; La figure 4 est un graphique qui illustre des données sismiques brutes à changements de polarité, collectées à partir de plusieurs récepteurs situés à la 20 surface de la Terre ; La figure 5 est un graphique qui illustre des données sismiques corrigées dans le temps, qui correspondent à la figure 4 ; La figure 6 est un graphique qui illustre une 25 trace associée à un évènement sismique en cas de changement de polarité ; La figure 7 est un graphique qui illustre des données sismiques brutes ayant un niveau de bruit élevé ; 30 La figure 8 est un graphique qui illustre les données sismiques corrigées dans le temps de la figure 7 ; La figure 9 est un graphique qui illustre une trace associée à un évènement sismique illustré sur la 35 figure 7 ; La figure 10 est un schéma d'un système d'enregistrement de données sismiques brutes associées à un évènement sismique selon un exemple de mode de réalisation ; La figure 11 est un schéma de plusieurs récepteurs situés à la surface de la Terre afin d'enregistrer un évènement sismique ; La figure 12 est un organigramme d'un procédé de détermination d'un mécanisme source associé à un 10 évènement sismique selon un exemple de mode de réalisation ; La figure 13 est un schéma qui illustre les emplacements d'un évènement sismique et des récepteurs selon un exemple de mode de réalisation ; 15 Les figures 14a à 14d sont des schémas de modèles associés à différents mécanismes sources selon un exemple de mode de réalisation ; Les figures 15a à 15c sont des graphiques qui illustrent des crêtes calculées en présence d'un faible 20 bruit afin de déterminer un mécanisme source d'un évènement microsismique ; Les figures 16a à 16c sont des graphiques qui illustrent des crêtes calculées en présence d'un bruit élevé afin de déterminer un mécanisme source d'un 25 évènement microsismique ; La figure 17 est un organigramme d'un procédé de calcul d'un mécanisme source d'un évènement microsismique selon un exemple de mode de réalisation ; et 30 La figure 18 est un schéma d'un dispositif de calcul capable d'exécuter le procédé illustré sur la figure 17. Description détaillée La description suivante des exemples de modes de réalisation se rapporte aux dessins joints. Les mêmes numéros de référence sur les différents dessins identifient des éléments identiques ou similaires. La description détaillée suivante ne limite en aucun cas l'invention. L'étendue de l'invention est en revanche définie par les revendications jointes. Les modes de réalisation suivants sont décrits, pour des raisons de simplicité, pour des évènements microsismiques générés en sous-sol. Cependant, les modes de réalisation décrits ci-après ne se limitent pas à ces évènements, mais peuvent être appliqués à d'autres sources d'évènements sismiques. La référence, tout au long de la demande, à « un mode de réalisation » signifie qu'une caractéristique ou une structure particulière décrite en relation avec un mode de réalisation est incluse dans au moins un mode de réalisation de l'objet de l'invention. Ainsi, l'apparition de l'expression « dans un mode de réalisation » en divers endroits de la demande ne fait pas nécessairement référence au même mode de réalisation. En outre, les caractéristiques et les structures particulières peuvent être combinées de n'importe quelle manière adéquate dans un ou plusieurs mode(s) de réalisation. Selon un exemple de mode de réalisation, il existe un procédé qui localise des évènements microsismiques observés dans une sous-surface de la Terre. Le procédé 30 divise en sous-blocs une surface sur laquelle se trouvent les récepteurs. Pour chaque emplacement source possible, le procédé applique des retards de propagation selon un modèle de vitesse de sous-surface. Ensuite, il additionne les traces retardées associées à 35 chaque bloc. Les sous-blocs de la surface sont associés à des modèles de polarisation d'ondes. Chaque modèle est associé à un ou plusieurs mécanisme(s) source. Une fonction de semblance est calculée en additionnant les sommes pour chaque sous-bloc. Ainsi, différents modèles de polarisation d'ondes peuvent être testés pour un coût de calcul marginal. Le modèle qui génère la crête la plus élevée est sélectionné et indique le mécanisme source correct. On remarque que le procédé détermine le mécanisme source au niveau de détection de l'évènement, alors que les procédés existants détectent tout d'abord l'emplacement de l'évènement, puis déterminent le mécanisme source. Cette nouvelle solution peut détecter et déterminer le mécanisme source même pour des évènements insignifiants qui sont ratés par les procédés existants en raison de la dégradation du SNR. Ainsi, ce nouveau procédé fonctionne dans les environnements bruyants et augmente légèrement le coût de calcul de l'algorithme de semblance. Le nouveau procédé peut offrir un gain de plusieurs dB, comme cela est décrit ci-après. Avant de décrire en détail le nouvel algorithme, une brève présentation du matériel impliqué dans cette identification de mécanisme source s'impose. Un système 10 d'enregistrement et/ou de détermination des positions des évènements microsismiques est illustré sur la figure 10. Le système 10 est déployé au-dessus d'une sous-surface qui comprend des couches géologiques 12, 14 et 16. La couche 16 est une couche qui contient des hydrocarbures. Un puits 18 peut être foré à travers les couches géologiques, jusqu'à la couche 16. Le système 10 comprend un dispositif de fracture 20 et un dispositif de surveillance 50. Le dispositif de fracture 20 comprend une colonne d'injection de fluide 22 qui s'étend dans le puits 18 et une unité de pompage 24 positionnée sur la surface 26 de la sous-surface. La colonne d'injection 22 comprend une extrémité inférieure 28 munie d'ouvertures 30 et une extrémité supérieure 32 reliée à l'unité de pompage 24. L'unité de pompage 24 comprend une pompe 34 configurée pour injecter un fluide de fracture à haute pression 36 dans la couche 16 via la colonne d'injection 22. Le fluide de fracture se compose généralement d'un mélange de liquide (comme de l'eau) et de particules (comme du sable). Le fluide de fracture est évacué par la colonne 22 via les ouvertures 30 créées par des perforations dans le revêtement du puits. Le fluide de fracture pénètre dans la couche de la sous-surface 16, induisant une fracture dans la couche 16, c'est-à-dire l'apparition de craquelures à l'intérieur de la couche 16. Les particules contenues dans le fluide de fracture sont forcées dans les craquelures et y restent en place lorsque la pression du fluide est relâchée, formant ainsi un réseau perméable qui permet aux hydrocarbures de s'écouler dans la couche 16, jusqu'au puits 18. Le dispositif de surveillance 50 peut être placé séparément du dispositif de fracture 24, ou sur le dispositif de fracture 24. Le dispositif de surveillance peut comprendre un réseau 52 de récepteurs (comme des géophones) 54, une unité d'enregistrement 56 et une unité de traitement 58. Le réseau de récepteurs 52 comprend une pluralité de récepteurs 54 disposés à la surface du sol 26 ou à proximité de cette surface. Les récepteurs 54 peuvent être disposés sur les noeuds d'un réseau 60, comme cela est illustré sur la figure 11, ou peuvent être disposés selon n'importe quelle autre configuration. L'unité de traitement 58 peut être configurée pour exécuter n'importe lequel des procédés décrits ci-après. - 14 Les mouvements du sol détectés par les récepteurs 54 sont transformés en signaux (comme des tensions électriques) et sont transmis à l'unité d'enregistrement 56 afin d'y être enregistrés. Les 5 enregistrements obtenus représentent les données sismiques brutes. Par exemple, selon un exemple de mode réalisation, les signaux ri,j (i indique une coordonnée sur l'axe X et j indique une coordonnée sur l'axe Y) sont détectés par les récepteurs 54. Le signal « r » 10 est une fonction de la position du récepteur et du moment t auquel il est enregistré. Selon un exemple de mode de réalisation, des données brutes sont collectées à l'étape 1200 (voir figure 12) à partir des récepteurs 54, et ces données 15 sont fournies à l'unité de traitement 58. À l'étape 1202, les données reçues sont retardées (corrigées dans le temps) sur la base d'une position source sélectionnée 70. En d'autres termes, comme cela est illustré sur la figure 13, si la source est considérée 20 comme se trouvant à la position 70, une trace enregistrée doit alors être corrigée sur la base d'un retard correspondant qui correspond à un trajet 72 de l'onde de pression qui se propage de l'emplacement de la source 70 jusqu'au récepteur. Dans un autre mode 25 de réalisation, des ondes partagées peuvent être considérées, avec ou à la place des ondes de pression. Chaque signal enregistré peut être corrigé selon l'emplacement de la source. Ce fait indique que l'identification du mécanisme source est effectuée 30 pendant que l'emplacement de la source est déterminé. À l'étape 1204, la surface 76 (voir figure 13) qui comprend les récepteurs 54 est divisée en sous-blocs 78. On remarque que cette étape peut être effectuée une seule fois, avant l'étape 1200. Par exemple, si un 35 groupe de récepteurs possède N = 40 x 40 récepteurs, ce groupe peut être divisé en 25 sous-blocs, comme cela est illustré sur la figure 14a. Dans une application, le groupe peut être divisé en 20 à 30 blocs. Les sous-blocs sont numérotés nl à n25, et chaque sous-bloc peut comprendre 40 x 40 / 25 = 64 récepteurs. Dans une autre application, un sous-bloc peut comprendre un nombre de récepteurs différent d'un autre sous-bloc. Dans ce cas particulier (utilisé uniquement pour illustrer le nouveau concept et en aucun cas destiné à limiter les nouveaux modes de réalisation), il est supposé que les sous-blocs 80 à fond blanc (voir les figures 14a à 14d) comprennent des récepteurs qui sont affectés par une polarisation positive, et que les sous-blocs 82 à fond gris comprennent des récepteurs qui sont affectés par une polarisation négative. Ainsi, il est possible d'avoir 225 modèles possibles pour les 25 sous-blocs illustrés sur les figures 14a à 14d. Cependant, la majorité des modèles possibles (voir par exemple celui illustré sur la figure 14d) ne correspond à aucun mécanisme source. Dans la pratique, environ 20 à 30 modèles sont associés à un mécanisme source spécifique. Des exemples de ces modèles sont illustrés sur les figures 14b et 14c. Dans une application, le nombre total de modèles est supérieur à 1 000 000 et le nombre de modèles (nombre prédéterminé) associés à un mécanisme source spécifique est inférieur à 100. Dans une application, il existe une correspondance unique entre le nombre de mécanismes source et le nombre de modèles réels. Dans une autre application, plusieurs modèles peuvent correspondre au même mécanisme source. À l'étape 1206, les données sismiques corrigées dans le temps des récepteurs provenant de chaque sous-bloc sont additionnées. Par exemple, si l'on considère 35 le sous-bloc nl, la somme des traces est donnée par : - 16 n=64 = j i=1 où rj.,i est la trace qui provient du récepteur 5 « i » dans le sous-bloc nl. Pour le sous-bloc nk, la somme est donnée par : n=60 Sk = Erij i=i 10 La semblance classique peut être calculée en additionnant les traces de tous les sous-blocs, comme cela est illustré ci-dessous : N=1600 Semb _Nom = Erg . 15 Cela équivaut à calculer : 25 Semb _Nom =Es, k=1 20 Pour la sommation classique, 1599 sommations sont effectuées. Cependant, pour le nouveau procédé, 25 sous-blocs ayant chacun 64 traces sont calculés, et les 25 résultats sont additionnés, ce qui ajoute 24 opérations supplémentaires. Il y a ainsi 25 x 63 + 24 25 opérations, soit un total de 1599 opérations. Dans une application, les sous-blocs susceptibles de provoquer du bruit lors du calcul de semblance (comme les sous-blocs qui correspondent à une inversion de polarité) peuvent être ignorés. Plus généralement, dans une autre application, chaque élément de la semblance peut être multiplié par un coefficient afin d'améliorer les performances de l'algorithme. Le coefficient peut agir comme une pondération et peut être un coefficient basé sur la puissance sonore, un coefficient d'apodisation, etc. Une fois les éléments de sommation de chaque sous-bloc obtenus, le procédé illustré sur la figure 12 passe à l'étape 1208, à laquelle différents modèles sont sélectionnés. Les modèles peuvent être sélectionnés, par exemple, sur la base d'une précédente détection de mécanisme source. Les modèles peuvent être stockés dans un dispositif de stockage afin d'être utilisés ultérieurement. À l'étape 1210, pour chaque modèle, une fonction de semblance correspondante est calculée, comme cela est décrit ci-après. Par exemple, pour le modèle illustré sur la figure 14b, un nominateur Pattern n Nom est calculé comme suit : Pattern _n _Nom -s,+ s2+ s, + + 86 + s6 - - sa + s9 + s,, + - SI3 SI4 SI7 SI6 SI7 + S20 + 821 + 822 + S23 + S24 e2s, qui représente la sommation de tous les sous-blocs. On remarque que chaque élément de sommation est multiplié par un coefficient de polarisation (plus ou 25 moins un) selon un modèle correspondant. On remarque que PatternnNom peut être reformulé comme suit : Pattern _n _Nom = Semb _Nom- 2s, - 2s, -2s, - 2s,, - 28 1, - 231, -2s' 30 ce qui représente seulement 18 sommations supplémentaires. Ainsi, en déduisant un modèle donné à partir d'un autre, il est possible de réduire le nombre de sommations et, implicitement, le coût de calcul. Étant donné que le dénominateur de la fonction de semblance n'est pas affecté par les signes de polarisation des sous-blocs, n semblances peuvent être calculées, ce qui correspond aux n modèles sélectionnés. La semblance pour un modèle donné « n » est obtenue par : Semb _n= Pattern - n - Nom2 Semb _Denon: Une fois que les fonctions de semblance ont été calculées pour les n modèles considérés, celui qui possède la valeur maximale est sélectionné à l'étape 1212 comme identifiant le mécanisme source pour la position source considérée. Avec moins de 30 modèles possibles, le nouveau procédé réduit de manière significative le coût de calcul qui permet de déterminer le mécanisme source. Les modèles peuvent être prédéterminés par ceux qui analysent les données sismiques, par exemple, sur la base de leur expérience. Dès que les modèles ont été prédéterminés, ils sont stockés dans la mémoire du dispositif d'analyse afin d'être utilisés ultérieurement. En outre, on remarque que les étapes illustrées sur la figure 12 peuvent être répétées lorsque l'emplacement de la source est modifié. Les résultats produits par le nouveau procédé décrit ci-dessus sont désormais comparés avec des algorithmes existants. Trois algorithmes de détection de source sont considérés en présence d'un mécanisme source. Des données synthétiques sont considérées pour un groupe de récepteurs en 2 dimensions. Un premier algorithme utilise la semblance classique, le second est le nouveau mécanisme décrit ci-dessus, et le troisième utilise des sous-blocs avec des valeurs absolues (aucun modèle identifié), comme cela est décrit dans la section « Contexte ». Lorsque le niveau de bruit est faible, les trois algorithmes identifient correctement l'emplacement de la source, comme cela est à 15c. On remarque que l'algorithme de semblance illustré sur les figures 15a la figure 15a correspond à classique, que la figure 15b correspond au nouvel algorithme, et que la figure 15c correspond à l'algorithme qui attribue des valeurs absolues pour chaque sous-bloc. La position maximale P sur chaque graphique fournit les coordonnées X et Z de la source. Les valeurs maximales sont normalisées à l'aide des valeurs maximales de semblance classiques. Les graphiques indiquent tous une erreur de position de la source, qui est inférieure à 20 mètres pour le second et le troisième algorithmes. Cependant, la situation est assez différente lorsque le niveau de bruit est élevé, comme cela est illustré sur les figures 16a à 16c. Dans ce cas, les deux algorithmes classiques illustrés sur les figures 16a et 16c offrent plusieurs crêtes et de fortes erreurs associées aux crêtes. En raison des crêtes multiples, ces algorithmes ne sont pas capables de déterminer la position de la source et, ainsi, implicitement, du mécanisme source. À l'inverse, comme cela est illustré sur la figure 16b, le nouvel algorithme possède une crête claire et son erreur est toujours inférieure à 20 mètres. Ainsi, les exemples ci-dessus indiquent que le nouvel algorithme est capable, même lorsque le niveau de bruit est élevé, de déterminer correctement l'emplacement de la source et du mécanisme source associé, et toutes ces caractéristiques sont réalisées avec une augmentation marginale du coût de calcul. Le nouvel algorithme susmentionné peut être mis en oeuvre comme un procédé illustré sur la figure 17. Le procédé illustré sur la figure 17 comprend une étape 1700 de réception de traces associées à une pluralité de récepteurs situés dans un réseau au-dessus de la sous-surface ; une étape 1702 de division du réseau en un ensemble de sous-blocs, chaque sous-bloc comprenant un ou plusieurs récepteur(s) de la pluralité de récepteurs ; une étape 1704 de détermination, pour chaque sous-bloc, d'un élément de sommation correspondant sur la base des traces associées au sous-bloc correspondant ; une étape 1706 d'identification d'un nombre prédéterminé de modèles qui sont associés à plusieurs mécanismes source ; une étape 1708 de calcul, dans le dispositif de calcul, pour chaque modèle, d'une fonction de semblance sur la base des éléments de sommation des sous-blocs ; et une étape 1710 d'identification du mécanisme source de l'évènement microsismique sur la base de la fonction de semblance. Un modèle comprend un premier sous-ensemble de sous-blocs à polarité positive et un second sous-ensemble de sous-blocs à polarité négative, le premier sous-ensemble et le second sous-ensemble formant l'ensemble de sous-blocs entier. Les procédés décrits ci-dessus peuvent être appliqués dans des dispositifs dédiés (comme des réseaux, des ordinateurs ou des réseaux de cloud computing dédiés, etc.) afin d'être exécutés. Une combinaison de logiciel et de matériel peut être utilisée pour mettre en oeuvre les procédés décrits ci-dessus. Ledit dispositif peut être un dispositif de calcul configuré pour identifier un mécanisme source associé à un évènement microsismique qui a lieu au niveau d'une 35 sous-surface de la Terre. Le dispositif de calcul - 21 comprend un processeur qui est configuré pour recevoir des traces associées à une pluralité de récepteurs situés sur un réseau au-dessus de la sous-surface ; diviser le réseau en un ensemble de sous-blocs, chaque sous-bloc comprenant un ou plusieurs récepteur(s) de la pluralité de récepteurs ; déterminer, pour chaque sous-bloc, un élément de sommation correspondant sur la base des traces associées au sous-bloc correspondant ; recevoir un nombre prédéterminé de modèles qui sont associés à plusieurs mécanismes source ; multiplier chaque élément de sommation par un coefficient de polarisation selon un modèle correspondant ; calculer, pour chaque modèle, une fonction de semblance sur la base des éléments de sommation des sous-blocs ; et identifier le mécanisme source de l'évènement microsismique sur la base de la fonction de semblance. Un modèle comprend un premier sous-ensemble de sous-blocs à polarité positive et un second sous-ensemble de sous-blocs à polarité négative, le premier sous- ensemble et le second sous-ensemble formant l'ensemble de sous-blocs entier. Le nombre prédéterminé de modèles peut être inférieur à 100, et le nombre total de modèles possibles peut être supérieur à 1 000 000. Le processeur peut en outre être configuré pour : recevoir des données sismiques brutes de la part des récepteurs ; corriger, pour chaque position source, un temps de propagation correspondant à un trajet de récepteur source, selon un modèle de vitesse de sous- surface donné afin d'obtenir des données sismiques corrigées ; et additionner les données sismiques corrigées afin d'obtenir les traces. Optionnellement, le processeur peut être configuré pour additionner les traces correspondantes de chaque sous-bloc afin de former l'élément de sommation correspondant pour chaque sous-bloc, l'ensemble de sous-blocs pouvant comprendre entre 20 et 30 sous-blocs. Dans une application, il existe une correspondance unique entre le nombre prédéterminé de modèles et les 5 mécanismes sources. Optionnellement, le processeur est en outre configuré pour comparer les résultats de la fonction de semblance pour le nombre de modèles prédéterminé ; et pour sélectionner le modèle qui possède la valeur la plus élevée pour la fonction de 10 semblance, comme étant le modèle qui correspond au mécanisme source, ou pour considérer un emplacement donné de l'évènement microsismique ; et pour sélectionner uniquement les traces associées à l'emplacement donné. 15 Une machine dédiée qui peut mettre en oeuvre un ou plusieurs des exemples de modes de réalisation décrits ci-dessus va maintenant être décrite en référence à la figure 18 Un exemple de dispositif de calcul 1800 qui 20 convient pour réaliser les activités décrites dans les exemples de modes de réalisation peut comprendre un serveur 1801. Ledit serveur 1801 peut comprendre un processeur central (CPU) 1802 couplé à une mémoire vive (RAM) 1804 et à une mémoire morte (ROM) 1806. La ROM 25 1806 peut également être d'autres types de supports de stockage destinés à stocker des programmes, comme une ROM programmable (PROM), une PROM effaçable (EPROM), etc. Le processeur 1802 peut communiquer avec d'autres composants internes et externes par le biais de 30 circuits d'entrée/de sortie (I/O) 1808 et de bus 1810, afin de fournir des signaux de commande et similaires. Le processeur 1802 exécute une variété de fonctions, comme cela est connu dans l'art, conformément à des instructions de logiciels et/ou de firmware. Le serveur 1801 peut également comprendre un ou plusieurs dispositif(s) de stockage de données, comme des disques durs et des disquettes 1812, des CD-ROM 1814, et d'autres dispositifs capables de lire et/ou de stocker des informations, comme des DVD, etc. Dans un mode de réalisation, un logiciel d'exécution des étapes décrites ci-dessus peut être stocké et distribué sur un CD-ROM 1816, une disquette 1818 ou toute autre forme de support capable de stocker des informations de manière portative. Ces supports de stockage peuvent être insérés dans, et lus par, des dispositifs tels qu'un lecteur de CD-ROM 1814, un lecteur de disque 1812, etc. Le serveur 1801 peut être couplé à un afficheur 1820, qui peut être n'importe quel type d'afficheur ou d'écran de présentation connu, comme un afficheur LCD, un afficheur plasma, des tubes cathodiques (CRT), etc. une interface utilisateur d'entrée 1822 est prévue, et comprend un ou plusieurs mécanisme(s) d'interface utilisateur, comme une souris, un clavier, un microphone, un pavé tactile, un écran tactile, un système de reconnaissance vocale, etc. Le serveur 1801 peut être couplé à d'autres dispositifs de calcul, comme des terminaux terrestres et/ou sans fil, via un réseau. Le serveur peut faire partie d'un plus grand réseau, comme un réseau global (GAN), tel qu'Internet 1828, qui permet une connexion ultime aux nombreux clients terrestres et/ou mobiles. Comme cela sera également apprécié par l'homme du métier, les exemples de modes de réalisation peuvent 30 être intégrés à un dispositif de communication sans fil ou à un réseau informatique en tant que procédé, ou à un programme informatique. Par conséquent, les exemples de modes de réalisation peuvent prendre la forme d'un mode de réalisation totalement matériel ou d'un mode de 35 réalisation qui allie des aspects matériels et logiciels. En outre, les exemples de modes de réalisation peuvent prendre la forme d'un programme informatique stocké sur un support de stockage lisible par ordinateur, qui comprend des instructions lisibles par ordinateur. N'importe quel support lisible par ordinateur peut être utilisé, y compris des disques durs, des CD-ROM, des DVD, des dispositifs de stockage optiques, ou des dispositifs de stockage magnétiques tels que des disquettes ou des bandes magnétiques. D'autres exemples non limitatifs de supports lisibles par ordinateur comprennent des mémoires Flash ou d'autres mémoires connues. Les exemples de modes de réalisation décrits prévoient un dispositif de calcul, un logiciel et un procédé de détermination d'un mécanisme source associé à un évènement microsismique. Il doit être compris que cette description n'est pas destinée à limiter l'invention. A l'inverse, les exemples de modes de réalisation sont destinés à couvrir les variantes, les modifications et les équivalents, qui font partie de l'esprit et de l'étendue de l'invention définie par les revendications jointes. En outre, dans la description détaillée des exemples de modes de réalisation, plusieurs détails spécifiques sont fournis afin de permettre une parfaite compréhension de l'invention. Cependant, l'homme du métier comprendra que différents modes de réalisation peuvent être mis en pratique sans lesdits détails spécifiques. Bien que les caractéristiques et les éléments des 30 présents exemples modes de réalisation soient décrits dans des combinaisons particulières dans les modes de réalisation, chaque caractéristique ou élément peut être utilisé(e) seul(e), sans les autres caractéristiques et éléments des modes de réalisation, ou dans différentes combinaisons, avec ou sans les autres caractéristiques et éléments décrits ici. La présente description écrite utilise des exemples de l'objet décrit afin de permettre à l'homme 5 du métier de le mettre en pratique, y compris de fabriquer et d'utiliser n'importe quel(s) dispositif(s) ou système(s) et d'exécuter n'importe quel procédé inclus. L'étendue brevetable de l'objet de la présente invention est définie par les revendications, et peut 10 comprendre d'autres exemples qui apparaîtront à l'homme du métier. Ces autres exemples sont destinés à faire partie de l'étendue des revendications.