FR2859539A1 - Procede de resistivite electromagnetique directionnelle - Google Patents

Procede de resistivite electromagnetique directionnelle Download PDF

Info

Publication number
FR2859539A1
FR2859539A1 FR0411139A FR0411139A FR2859539A1 FR 2859539 A1 FR2859539 A1 FR 2859539A1 FR 0411139 A FR0411139 A FR 0411139A FR 0411139 A FR0411139 A FR 0411139A FR 2859539 A1 FR2859539 A1 FR 2859539A1
Authority
FR
France
Prior art keywords
resistivity
formation
directional
measurements
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR0411139A
Other languages
English (en)
Other versions
FR2859539B1 (fr
Inventor
Qiming Li
Lawrence Chou
Dzevat Omeragic
Libo Yang
Alain Dumont
Lingyun Hu
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Services Petroliers Schlumberger SA
Original Assignee
Services Petroliers Schlumberger SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Services Petroliers Schlumberger SA filed Critical Services Petroliers Schlumberger SA
Publication of FR2859539A1 publication Critical patent/FR2859539A1/fr
Application granted granted Critical
Publication of FR2859539B1 publication Critical patent/FR2859539B1/fr
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/26Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
    • G01V3/28Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/30Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves

Abstract

Une nouvelle technique de traitement de données à la volée est utile pour extraire des signaux à partir de la variation azimutale des mesures directionnelles acquises par un outil de diagraphie dans un sondage. Les signaux de frontière, d'anisotropie et de fracture appropriés sont extraits de la réponse de formation par ajustement de la variation azimutale des tensions mesurées à certaines fonctions sinusoïdales. L'orientation de la stratification est également obtenue en tant que résultat. Les signaux directionnels extraits sont utiles pour obtenir des distances de frontière et prendre des décisions de géoguidage. Deux techniques mettant en oeuvre l'inversion et la génération de tracé croisé peuvent être utilisées, suivant la nature de la frontière. Une interface d'utilisateur graphique (GUI) fait partie d'un système pour faciliter la définition flexible d'objectifs d'inversion, pour améliorer les résultats d'inversion, et pour la visualisation du modèle de formation ainsi que des mesures d'inversion.

Description

<Desc/Clms Page number 1>
PROCEDE DE RESISTIVITE ELECTROMAGNETIQUE DIRECTIONNELLE CONTEXTE DE L'INVENTION Domaine de l'invention
La présente invention concerne généralement le domaine de la diagraphie de sondage. Plus particulièrement, l'invention concerne des techniques améliorées dans lesquelles des instruments équipés de systèmes d'antenne ayant des représentations de dipôle magnétique transversal ou incliné sont utilisés pour des mesures électromagnétiques de formations souterraines et pour positionner des sondages par rapport à des frontières géologiques dans un réservoir.
L'invention a une application générale dans l'art de la diagraphie de sondage, mais est particulièrement utile dans une opération de diagraphie en cours de forage.
Etat de l'art antérieur
Différentes techniques de diagraphie de sondage sont connues dans le domaine de la prospection et la production d'hydrocarbures. Ces techniques utilisent typiquement des instruments ou des outils équipés de sources adaptées pour émettre de l'énergie dans une formation souterraine qui a été pénétrée par un sondage. Dans cette description, "instrument" et "outil" sont utilisés de manière interchangeable pour désigner, par exemple, un instrument (ou outil) électromagnétique, un outil (ou instrument) de travail au câble, ou un outil (ou instrument) de diagraphie en cours de forage. L'énergie émise interagit avec la
<Desc/Clms Page number 2>
formation environnante pour produire des signaux qui sont ensuite détectés et mesurés par un ou plusieurs capteurs. En traitant les données de signal détecté, un profil des propriétés de la formation est obtenu.
Les diagraphies par induction et par propagation électromagnétique (EM) sont des techniques connues. Les instruments de diagraphie sont disposés dans un sondage pour mesurer la conductivité électrique (ou son inverse, la résistivité) de formations rocheuses entourant le sondage. Dans la présente description, toute référence à la conductivité est destinée à inclure son inverse, la résistivité, ou réciproquement.
Un outil de résistivité électromagnétique type comprend une antenne émettrice et une ou plusieurs (typiquement une paire) d'antennes réceptrices disposées à une certaine distance de l'antenne émettrice le long de l'axe de l'outil (voir la figure 1).
Les outils d'induction mesurent la résistivité (ou conductivité) de la formation en mesurant la tension induite dans la ou les antenne(s) réceptrice(s) due au flux magnétique induit par les courants circulant dans l'antenne émettrice (ou de transmission). Un outil de propagation EM fonctionne de manière similaire mais, typiquement, à des fréquences plus élevées que les outils à induction pour des espacements d'antenne comparables (environ 106 Hz pour les outils à propagation comparé à environ 104 Hz pour les outils à induction). Un outil à propagation type peut fonctionner dans une plage de fréquence de 1 kHz à 2 MHz.
<Desc/Clms Page number 3>
Les émetteurs et récepteurs conventionnels sont des antennes formées de bobines constituées d'un ou plusieurs tours de fil conducteur isolé enroulés autour d'un support. Ces antennes peuvent typiquement fonctionner en tant que sources et/ou récepteurs. Il apparaîtra à l'homme du métier que la même antenne peut être utilisée en tant qu'émetteur à un moment et en tant que récepteur à un autre. Il apparaîtra également que les configurations d'émetteur-récepteur présentement décrites sont interchangeables en raison du principe de réciprocité, c'est-à-dire que "l'émetteur" peut être utilisé en tant que "récepteur", et vice-versa.
Les antennes fonctionnent selon le principe qu'une bobine portant un courant (par exemple, une bobine émettrice) génère un champ magnétique. L'énergie électromagnétique provenant de l'antenne émettrice est émise dans la formation environnante et cette émission induit des courants de Foucault circulant dans la formation autour de l'émetteur (voir la figure 2A). Les courants de Foucault induits dans la formation, qui sont des fonctions de la résistivité de la formation, génèrent un champ magnétique qui, à son tour, induit une tension électrique dans les antennes réceptrices.
Si une paire de récepteurs espacés est utilisée, les tensions induites dans les deux antennes réceptrices auront des phases et des amplitudes différentes en raison de la propagation géométrique et l'absorption par la formation environnante. La différence de phase (déphasage #) et de rapport d'amplitude (atténuation, A) des deux récepteurs peuvent être utilisées pour
<Desc/Clms Page number 4>
déduire la résistivité de la formation. Le déphasage (#) et l'atténuation (A) détectés dépendent non seulement de l'espacement entre les deux récepteurs et des distances entre l'émetteur et les récepteurs, mais également de la fréquence des ondes EM générées par l'émetteur.
Dans des instruments de diagraphie à induction et à propagation conventionnels, les antennes émettrices et réceptrices sont montées avec leurs axes le long de l'axe longitudinal de l'instrument. Par conséquent, ces outils sont mis en oeuvre avec des antennes ayant des représentations de dipôle magnétique longitudinal (LMD). Une technique émergente dans le domaine de la diagraphie de sondage est l'utilisation d'instruments comprenant des antennes ayant des bobines inclinées ou transversales, c'est-à-dire que l'axe de la bobine n'est pas parallèle à l'axe longitudinal de l'outil.
Ces instruments sont donc mis en oeuvre avec une antenne à dipôle magnétique transversal ou incliné (TMD). Il apparaîtra à l'homme du métier que différentes manières sont possibles pour incliner ou biaiser une antenne. Des instruments de diagraphie équipés d'antennes TMD sont décrits, par exemple, dans les brevets U.S. n 6 163 155 ; 147 496 ; 5 115 198 ; 4 319 191 ; 5 508 616 ; 5 757 191 ; 5 781 436 ; 6 044 325 ; et 6 147 496.
La figure 2A présente une représentation simplifiée de courants de Foucault et d'énergie électromagnétique (EM) circulant depuis un instrument de diagraphie disposé dans une partie ou un segment de sondage qui pénètre dans une formation souterraine dans
<Desc/Clms Page number 5>
une direction perpendiculaire aux couches de sédimentation. Ce n'est cependant pas une description exacte de tous les différents segments qui constituent un sondage - en particulier lorsque le sondage a été foré de manière directionnelle comme décrit ci-dessous. A savoir, des segments d'un sondage pénètrent souvent dans des couches de formation à un angle autre que 90 degrés, comme décrit sur la figure 2B. Lorsque cela se produit, le plan de formation est dit avoir une relation de pendage. Un angle de pendage relatif, #, est défini comme étant l'angle entre l'axe du sondage (axe d'outil) BA et la normale N par rapport au plan P d'une couche de formation souhaitée.
Les techniques de forage connues dans l'art comprennent le forage de sondages depuis une position géographique choisie à la surface du sol, le long d'une trajectoire choisie. La trajectoire peut s'étendre à d'autres positions géographiques choisies à des profondeurs particulières dans le sondage. Ces techniques sont appelées collectivement techniques de "forage directionnel". Une application du forage directionnel est le forage de sondages fortement déviés (par rapport à la verticale) ou même horizontaux, des sondages dans et le long de formations rocheuses contenant des hydrocarbures relativement minces (appelées "zones productives") sur des distances étendues. Ces sondages fortement déviés sont censés augmenter fortement le drainage d'hydrocarbures depuis la zone productive par rapport aux sondages "conventionnels" qui pénètrent "verticalement" (de manière sensiblement perpendiculaire à l'agencement de
<Desc/Clms Page number 6>
couches de la formation, comme décrit sur la figure 2A) dans la zone productive.
Dans un puits de sondage fortement dévié ou horizontal dans une zone productive, il est important de maintenir la trajectoire du sondage de telle manière qu'il reste dans une position particulière dans la zone productive. Des systèmes de forage directionnel sont connus dans l'art qui utilisent des "moteurs à boues" et des "foreuses pliées", ainsi que d'autres moyens, pour contrôler la trajectoire d'un sondage par rapport à des références géographiques, telles que le nord magnétique, la gravité terrestre (verticale) et la vitesse de rotation de la terre (par rapport à l'espace cosmique). L'agencement de couches des formations, cependant, peut être tel que la zone productive ne soit pas située le long d'une trajectoire prévisible à des positions géographiques distantes de l'emplacement en surface du sondage. Typiquement, l'opérateur de sondage utilise des informations (telles que des diagraphies LWD) obtenues lors du forage de sondage pour maintenir la trajectoire du sondage dans la zone productive et pour vérifier que le sondage est, en fait, foré dans la zone productive.
Des techniques connues dans l'art pour maintenir une trajectoire sont décrites, par exemple, dans Tribe et al., Precise Well Placement using Rotary Steerable Systems and LWD measurement, SOCIETY OF PETROLEUM ENGINEERS, publication 71396,30 septembre 2001. La technique décrite dans cette référence est basée sur des réponses de capteur de conductivité LWD. Si, par exemple, la conductivité de la zone productive est
<Desc/Clms Page number 7>
connue avant la pénétration par le sondage et si les conductivités de zones sus-jacentes et sous-jacentes présentent un contraste significatif par rapport à la zone productive, une mesure de la conductivité de formation effectuée durant le forage peut être utilisée en tant que critère pour "guider" le sondage de manière à rester dans la zone productive. Plus spécifiquement, si la conductivité mesurée s'écarte significativement de la conductivité de la zone productive, c'est une indication que le sondage approche de, ou a même pénétré, l'interface de la formation rocheuse susjacente ou sous-jacente. A titre d'exemple, la conductivité d'un sable saturé en pétrole peut être significativement plus faible que celle d'une argile litée sous-jacente ou sus-jacente. Une indication que la conductivité en position adjacente au sondage augmente peut être interprétée comme une indication que le sondage approche de la couche de formation susjacente ou sous-jacente (de l'argile litée dans cet exemple). La technique de forage directionnel utilisant une mesure de propriété de formation pour guider l'ajustement de trajectoire est généralement appelée "géoguidage".
En plus des mesures EM, des mesures acoustiques et radioactives sont également utilisées en tant que moyen de géoguidage. En utilisant à nouveau l'exemple d'une zone de production de pétrole avec de l'argile litée sous-jacent et sus-jacent, la radioactivité gamma naturelle dans la zone productive est généralement nettement inférieure à l'activité de rayonnement gamma naturelle des formations d'argile litée au-dessus et
<Desc/Clms Page number 8>
au-dessous de la zone productive. En conséquence, une augmentation de l'activité de rayonnement gamma naturelle mesurée par un capteur de rayons gamma LWD indiquera que le sondage dévie du centre de la zone productive et approche ou pénètre l'interface supérieure ou inférieure avec l'argile litée.
Si, comme dans les exemples précédents, la conductivité et la radioactivité naturelle des formations d'argile litée sus-jacente et sous-jacente sont similaires l'une à l'autre, les techniques de géoguidage précédemment décrites indiquent uniquement que le sondage quitte la zone productive, mais n'indiquent pas si le sondage s'écarte de la zone productive par le sommet de la zone ou par la base de la zone. Cela pose un problème à l'opérateur de forage, qui doit corriger la trajectoire du forage pour maintenir la position choisie dans la zone productive.
Des instruments de diagraphie à induction EM sont correctement adaptés pour les applications de géoguidage parce que leur profondeur latérale (radiale) d'investigation dans les formations entourant le sondage est relativement élevée, en particulier par rapport aux instruments nucléaires. L'investigation radiale plus profonde permet aux instruments à induction de "voir" à une distance latérale (ou radiale) significative de l'axe du sondage. Dans les applications de géoguidage, cette plus grande profondeur d'investigation permet la détection de l'approche de frontières de couches de formation à des distances latérales plus élevées du sondage, ce qui laisse à l'opération de forage plus de temps pour
<Desc/Clms Page number 9>
effectuer d'éventuelles corrections de trajectoire. Les instruments de type à propagation conventionnels sont capables de détecter des variations axiales et latérales (radiales) de conductivité des formations entourant l'instrument, mais la réponse de ces instruments ne peut généralement pas résoudre les variations azimutales de la conductivité des formations entourant l'instrument. De plus, de tels instruments sont incapables de détecter une anisotropie dans des puits verticaux.
Deux marchés émergents importants rendent la résolution de ces inconvénients plus urgente. Le premier domaine émergent est le besoin croissant d'un positionnement de puits précis, qui requiert des mesures directionnelles pour prendre des décisions de guidage pour placer le sondage de manière optimale dans le réservoir. Le second est la réponse de résistivité faible dans des formations stratifiées où l'identification et la caractérisation précises de réserves d'hydrocarbures n'est pas possible sans connaître l'anisotropie de résistivité. De nombreux brevets récents décrivent des procédés et appareils pour effectuer des mesures directionnelles et obtenir une anisotropie de résistivité. Pour les applications de diagraphie en cours de forage, le brevet U.S. n 5 508 616 de Sato et al. décrit un outil à bobine de type à induction avec deux bobines inclinés dans des directions différentes non alignées avec l'axe longitudinal de l'outil. Le caractère directionnel de la mesure est illustré par un argument simple en ce que la fonction de sensibilité des deux bobines inclinées
<Desc/Clms Page number 10>
est concentrée vers la région de chevauchement de la zone de sensibilité de chaque bobine. Par rotation de l'outil, Sato et al. revendiquent qu'une image de résistivité azimutale plus profonde de la formation peut être obtenue. Cependant, cette référence de brevet ne donne aucun détail sur la manière d'obtenir la résistivité azimutale et ne décrit pas de techniques de détection/caractérisation de frontière requises pour prendre des décisions de géoguidage quantitatives.
Le brevet U.S. n 6 181 138 de Hagiwara et Song étend les bobines directionnelles fixes simples de Sato et al. à des triples bobines à induction orthogonales toutes colocalisées avec l'émetteur et le récepteur.
Aucune rotation de l'outil n'est indiquée comme étant requise, étant donné que la direction de focalisation peut être ajustée dans une orientation arbitraire par combinaison linéaire des réponses des bobines orthogonales. Il n'est pas clairement indiqué s'il existe une conception d'écran qui permette le passage de toutes les composantes EM requises sans distorsion incontrôlable importante de la forme d'onde pour les applications "en cours de forage".
Le brevet U.S. n 6 297 639 de Clark et al., cédé au cessionnaire de la présente invention, décrit des procédés et un appareil pour effectuer des mesures directionnelles en utilisant différentes conceptions d'écran pour produire une atténuation choisie de l'énergie d'onde EM pour des bobines d'antennes axiales, inclinées et transversales. Cette référence de brevet décrit, entre autres choses, des mesures par induction et propagation directionnelles générales avec
<Desc/Clms Page number 11>
des bobines inclinées et des écrans appropriés, ainsi qu'un procédé pour conduire la compensation de sondage pour celles-ci qui n'est pas trivial. Une combinaison de bobines émettrice/réceptrice - une axiale et une inclinée - est explicitement décrite par Clark et al., ainsi que son application pour la détection de direction de frontière de couche en observant la variation azimutale du signal induit lorsque l'outil tourne. La variation azimutale du couplage peut être utilisée pour guider des puits en cours de forage. Des brevets d'écrans supplémentaires ont été cédés depuis, comprenant le brevet U.S. n 6 351 127 de Rosthal et al. et le brevet U.S. n 6 566 881 de Omeragic et al., tous deux étant cédés au cessionnaire de la présente invention.
Le brevet U.S. n 6 476 609 de Bittar étend un brevet d'anisotropie antérieur décrivant des émetteurs et des récepteurs pouvant avoir un angle d'inclinaison, le brevet U.S. n 6 163 155 également de Bittar, au domaine d'application du géoguidage. La réponse à une stratification d'un appareil à induction et propagation incliné haut/bas est décrite par la différence ou le rapport de signaux dans deux orientations différentes, mais aucun écran n'est mentionné. Les effets de l'anisotropie ou du pendage ne sont pas non plus considérés. Il manque également une description de la manière d'utiliser ces mesures pour dériver une distance précise par rapport à une frontière de couche de formation. Le brevet US 6 476609 suppose implicitement que l'orientation de la stratification est précisément connue de manière à calculer la réponse
<Desc/Clms Page number 12>
haute/basse. Aucune technique, cependant, n'est décrite pour localiser la direction vers le haut ou vers le bas précise avant le calcul des signaux directionnels hautbas.
La publication de demande de brevet U.S. n 2003/0085707 de Minerbo et al., cédée au cessionnaire de la présente invention, décrit des configurations d'outil et des techniques de symétrisation qui simplifient la réponse des mesures directionnelles au point qu'elle devienne pratiquement indépendante de l'anisotropie ou de l'angle de pendage.
Les réponses à la distance de frontière de couche avec différents pendages et anisotropies se chevauchent sensiblement, sauf à proximité de la frontière de couche. Des mesures de type à induction à deux bobines (un émetteur et un récepteur : "TR") et de type à propagation à trois bobines (un émetteur et deux récepteurs : "TRR") peuvent être symétrisées pour obtenir cette simplification. La symétrisation est effectuée entre deux paires TR inclinées de même espacement, mais l'angle incliné d'émetteur et l'angle incliné de récepteur étant échangés. Seuls les cas dans lesquels les moments magnétiques des émetteurs et des récepteurs sont situés dans le même plan sont pris en compte. Cela présente l'inconvénient de ne pas permettre d'obtenir le signal requis pour le géoguidage en permanence durant le coulissage, ce qui est le cas pour le positionnement de puits avec un moteur à boue durant la définition d'angle sur la trajectoire. S'il s'avère que le moment magnétique de l'outil est situé parallèlement à l'agencement de couches durant le
<Desc/Clms Page number 13>
coulissage, le signal directionnel haut/bas généré est égal à zéro indépendamment de la distance à la frontière. Par conséquent, aucune surveillance de la distance n'est possible.
La publication de brevet U.S. n 2003/0200029 de Omeragic et al., également cédée au cessionnaire de la présente invention, décrit des mesures directionnelles de type à propagation pour la détermination d'anisotropie dans des puits presque verticaux avec compensation de sondage. Des techniques d'inversion sont également utilisées pour obtenir la propriété anisotrope de la formation. La publication de brevet U. S. n 2003/0184302 de Omeragic et Esmersoy, cédée au cessionnaire de la présente invention, décrit également des techniques d'investigation avec des mesures directionnelles.
Les publications de brevet U.S. n 2004/0046560A1 et 2004/0046561A1 de Itzkovic et al., décrivent l'utilisation d'antennes quadripolaires, et des mesures de type couplage dipôle-quadripôle transversal et induction avec des caractéristiques directionnelles similaires à une réponse XZ de dipôle croisé conventionnelle. La réalisation pratique sur une massetige métallique et le blindage adéquat de telles antennes ne sont pas clairement décrits. De plus, l'effet sur le sondage de telles mesures et de ses interaction/couplage avec l'effet de frontière peuvent être différents des mesures de type XZ.
Aucune des références de brevet ci-dessus ne décrit l'utilisation de réponses azimutales du signal mesuré ou des techniques pour extraire de telles
<Desc/Clms Page number 14>
réponses. Ces références ne décrivent pas la manière d'utiliser la mesure directionnelle pour obtenir des distances à une frontière pour utilisation en géoguidage. Seule la mesure appelée haut/bas, qui est la différence de signal mesuré entre l'outil dirigé directement vers et depuis la couche de formation, est mentionnée. Les informations de pendage et d'azimut de stratification précises ne sont généralement pas connues avant le forage, et, de plus, elles varient fréquemment dans des situations de positionnement de puits complexes lorsqu'un géoguidage est requis.
L'utilisation d'une direction haut/bas de stratification prédéfinie produit, au mieux, une mesure dégradée et, au pire, peut conduire à des décisions de géoguidage erronées lorsque azimut de stratification change brusquement. En principe, les mesures peuvent être mises en fichier en fond de manière azimutale.
Cette technique présente plusieurs inconvénients comprenant des difficultés à aligner les fichiers supérieur et inférieur précisément avec l'orientation de la stratification de formation et l'incapacité à utiliser (c'est-à-dire, perdre) les données qui ne sont pas dans les fichiers haut et bas. La mémoire importante requise pour enregistrer les données azimutales avec une précision suffisante est également un problème.
Plus important, l'art existant pour le géoguidage utilisant des mesures directionnelles fonctionne uniquement pour le guidage vers le haut et vers le bas.
Il existe de nombreux cas dans lesquels le sondage doit
<Desc/Clms Page number 15>
se déplacer de manière azimutale pour éviter de sortir de la zone productive.
Il existe donc un besoin de procédés et de techniques d'extraction et d'analyse de la dépendance azimutale des mesures de diagraphie directionnelles, utilisant des mesures effectuées à tous les angles azimutaux, pour caractériser la formation rocheuse et pour guider les puits durant le forage avec une précision améliorée.
De plus, il existe un besoin d'obtenir azimut de stratification à partir des mesures directionnelles et de générer des mesures qui peuvent être utilisées pour le positionnement de puits en guidage haut/bas ou azimutal.
De plus, il existe un besoin de procédés d'utilisation de ces mesures directionnelles en temps réel pour obtenir des distances de frontière de couche et pour obtenir des modèles souterrains précis de telle manière que des décisions de géoguidage puissent être prises pour le positionnement de puits.
De plus, il existe un besoin d'un procédé de détection de la présence d'anisotropie de résistivité dans des couches de formation adjacentes à des puits sensiblement verticaux.
De plus, il existe un besoin d'un système efficace qui permette de telles mesures directionnelles, les analyse au fond et transmette les informations appropriées à la surface pour faciliter le géoguidage haut/bas ou azimutal lors du positionnement de puits.
De plus, il serait avantageux qu'un tel système puisse fournir des informations de distance à une frontière
<Desc/Clms Page number 16>
durant les phases de coulissage du forage (c'est-àdire, sans rotation du train de tiges de forage) ainsi que lorsque le système/ outil est en rotation.
DEFINITIONS
Certains termes sont définis tout au long de cette description tels qu'ils sont utilisés la première fois, tandis que certains autres termes utilisés dans cette description sont définis ci-dessous :
La "convergence" désigne l'état dans lequel des valeurs calculées de manière itérative sont proches de valeurs observées ou de limites finies au fur et à mesure que le nombre de cycles d'itération augmente.
Un "tracé croisé" désigne un graphique, ou la génération d'un tel graphique, qui représente la relation entre deux mesures différentes du même sujet ou échantillon.
"Inversion" ou "inverser" désigne la dérivation d'un modèle (également appelé "modèle d'inversion") à partir de données mesurées (par exemple, des données de diagraphie) qui décrit une formation souterraine et est cohérent avec les données mesurées.
"Face de coupe" désigne l'orientation angulaire d'un instrument autour de son axe longitudinal et représente un angle sous-tendu entre une référence choisie sur le logement de l'instrument (par exemple, une masse-tige de forage) et la paroi la plus haute sur le plan gravitationnel du forage ou le nord géographique.
"Symétrie" ou "symétrique", dans le présent contexte, désigne une configuration dans laquelle des
<Desc/Clms Page number 17>
ensembles d'agencements d'émetteur-récepteur sont disposés dans des orientations opposées le long de l'axe longitudinal d'un outil, de telle manière que les ensembles émetteur-récepteur puissent être corrélés à une opération de symétrie standard (par exemple, une translation, un plan de miroir, une inversion et une rotation) par rapport à un point sur l'axe de l'outil ou un plan de symétrie perpendiculaire à l'axe de l'outil.
RESUME DE L'INVENTION
Des aspects de la présente invention concernent la caractérisation de formation grâce à l'utilisation de tracés croisés. Un procédé comprend l'étape de génération de tracé croisé de deux mesures de diagraphie directionnelles acquises par un instrument disposé dans un sondage croisant la formation pour obtenir une distance à au moins une frontière de formation et une résistivité pour au moins une couche de formation. Le tracé croisé est réalisé en utilisant une frontière. La résistivité obtenue est la résistivité de couche d'épaulement et la distance obtenue est la distance la plus proche de la couche d'épaulement.
Un autre procédé comprend l'étape de génération de tracé croisé d'une résistivité et d'une mesure directionnelle déterminées en utilisant un instrument disposé dans un sondage croisant la formation pour obtenir une distance à au moins une frontière de formation et une résistivité pour au moins une couche de formation. Le tracé croisé est réalisé en utilisant
<Desc/Clms Page number 18>
un modèle à une frontière. La résistivité obtenue est la résistivité de couche d'épaulement et la distance obtenue est la distance la plus proche de la couche d'épaulement.
Un procédé supplémentaire comprend l'étape consistant à générer un tracé croisé d'une résistivité et deux mesures directionnelles déterminées en utilisant un instrument disposé dans un sondage croisant la formation pour obtenir une distance à au moins une frontière de formation et une résistivité pour au moins deux couches de formation. Le tracé croisé est réalisé en utilisant un modèle à une frontière. Les résistivités obtenues sont les résistivités de couche et de couche d'épaulement, et la distance obtenue est la distance la plus proche de la couche d'épaulement. La distance à la frontière et la résistivité de couche déterminées peuvent être utilisées pour prendre des décisions de forage.
Un autre aspect de la présente invention concerne l'utilisation d'une technique d'inversion pour interpréter les mesures directionnelles pour des applications de géoguidage. Présentement, l'algorithme d'ajustement itéractif est utile pour des mesures directionnelles en temps réel sélectionnées comme étant utiles en géoguidage. Un modèle d'inversion approprié est choisi pour les mesures directionnelles en temps réel sélectionnées. Une fois qu'il est vérifié que le modèle choisi est cohérent avec d'autres informations, il est utilisé pour prendre des décisions de forage.
<Desc/Clms Page number 19>
L'étape de sélection de modèle comprend l'exécution d'une pluralité de modèles comprenant des types de modèle : isotrope homogène (un seul paramètre : résistivité) ; anisotrope homogène (deux paramètres : Rh et Rv) ; formation isotrope à frontière unique, frontière au-dessus ou au-dessous (trois paramètres : Rcouche, Répaulement et distance à la frontière) ; formation anisotrope à frontière unique, frontière au-dessus ou au-dessous (quatre paramètres : Rcouche~h, Rcouche~v, Répaulement et distance à la frontière ; formation isotrope à deux frontières et trois paramètres : (cinq paramètres : Rcouche, Répaulement~haut, Répaulement~bas et distance à la frontière au-dessus et au-dessous de l'outil) ; formation anisotrope à deux frontières et trois paramètres : (six paramètres : Rcouche~h, Rcouche~v, Répaulement~~haut, Répaulement bas et distance à la frontière au-dessus et au-dessous de l'outil).
L'étape de sélection de modèle comprend, de préférence, la création d'une visualisation des mesures directionnelles sélectionnées.
Dans un mode de réalisation particulier, l'étape de sélection de modèle comprend l'identification de paramètres de formation connus, le choix interactif des modèles avec lesquels inverser les mesures directionnelles, et le choix du modèle le plus simple qui est ajuté aux informations connues.
L'étape de vérification de modèle comprend la comparaison du modèle choisi aux caractéristiques
<Desc/Clms Page number 20>
géologiques connues et aux autres paramètres de formation mesurés, et la mise à jour du modèle choisi si le modèle choisi n'est pas cohérent avec les informations connues.
L'étape de mise à jour comprend l'affinement du modèle choisi sur la base de l'un parmi les tendances, les connaissances antérieures, des informations externes et une combinaison de ceux-ci. Des paramètres d'inversion appropriés sont choisis et des plages pour les paramètres choisis sont définies. Le modèle d'inversion est, de préférence, mis à jour en ajoutant des couches de formation supplémentaires. Certaines des mesures directionnelles en temps réel peuvent être repondérées ou éliminées et les mesures directionnelles en temps réel résultantes réinversées dans le modèle mis à jour.
Les procédés de génération de tracé croisé selon la présente invention comprennent, de préférence, la définition d'un modèle approprié, le choix de mesures directionnelles appropriées, l'entrée des mesures sélectionnées dans le modèle défini pour générer le tracé croisé, et la génération d'une représentation graphique du tracé croisé. Le tracé croisé peut être mis à jour avec des mesures en temps réel supplémentaires.
BREVE DESCRIPTION DES DESSINS
Afin que les caractéristiques et avantages de la présente invention mentionnés ci-dessus puissent être compris en détail, une description plus particulière de l'invention, brièvement résumée ci-dessus, peut être
<Desc/Clms Page number 21>
faite en référence aux modes de réalisation de celle-ci qui sont illustrés dans les dessins annexés. Il doit être noté, cependant, que les dessins annexés illustrent uniquement des modes de réalisation types de cette invention et, par conséquent, ne doivent pas être considérés comme une limitation de sa portée, car l'invention peut admettre d'autres modes de réalisation également efficaces.
La figure 1 représente des schémas d'outils à induction et à propagation de l'art antérieur.
Les figures 2A et 2B sont des vues en élévation représentant des courants de Foucault induits par un outil de diagraphie dans un sondage pénétrant dans une formation sans et avec un pendage relatif, respectivement.
La figure 3 est une représentation en élévation d'un train de tiges de forage rotatif conventionnel dans lequel la présente invention peut être avantageusement utilisée.
La figure 4 est une représentation schématique d'un outil de diagraphie de mesure directionnelle de base ayant des paires d'antennes émettrices et réceptrices symétriques.
La figure 5A est une représentation schématique d'un outil de diagraphie de mesure directionnelle ayant une configuration TRR qui est insensible à l'anisotropie à un angle de pendage quelconque, selon un aspect de la présente invention.
La figure 5B représente des tracés de la réponse de propagation directionnelle pour une formation à
<Desc/Clms Page number 22>
trois couches en utilisant un outil de diagraphie selon la figure 5A.
La figure 6 représente une comparaison graphique d'une technique de mise en fichier et d'une technique d'ajustement à la volée selon un aspect de la présente invention.
La figure 7 présente des tracés représentant la convergence et l'erreur dans un ensemble de coefficients cibles obtenus par une technique d'ajustement selon un aspect de la présente invention, mise en oeuvre par un algorithme entier dans un processeur de signal numérique selon un autre aspect de la présente invention.
La figure 8 représente un tracé croisé représentant une mesure de résistivité conventionnelle et une mesure directionnelle utilisées pour obtenir la résistivité d'une couche de formation et la distance aux frontières de couche, pour une résistivité de couche d'épaulement connue.
La figure 9 représente les résultats d'une technique d'inversion appliquée pour obtenir la résistivité et les localisations de frontière point par point d'une couche de formation.
La figure 10 représente une technique d'inversion à base de tracé croisé appliquée pour interpréter les mesures de résistivité directionnelles.
La figure 11 est un organigramme pour le géoguidage selon un aspect de la présente invention.
La figure 12 est un organigramme pour un tracé croisé utile pour déterminer des distances de frontière de couche et l'application de telles distances en
<Desc/Clms Page number 23>
géoguidage en temps réel, selon des aspects supplémentaires de la présente invention.
La figure 13 représente un schéma d'organigramme représentant la structure d'une GUI pour afficher des modèles de couche de formation basés sur une inversion.
La figure 14 représente une visualisation générée par ordinateur d'une inversion de géoguidage associée à la GUI de la figure 13.
DESCRIPTION DETAILLEE DE L'INVENTION
La figure 3 représente un appareil de forage et un train de tiges de forage conventionnels dans lesquels la présente invention peut être avantageusement utilisée. Une plate-forme terrestre et un ensemble de derrick 10 sont positionnés sur un sondage 11 pénétrant dans une formation souterraine F. Dans le mélange de réaction illustré, le sondage 11 est formé par forage rotatif d'une manière qui est connue. Il apparaîtra à l'homme du métier dont cette description bénéficie, cependant, que la présente invention peut également être appliquée dans des applications de forage directionnel ainsi que de forage rotatif et n'est pas limité à des appareils de forage terrestres.
Un train de tiges de forage 12 est suspendu dans le sondage 11 et comprend un trépan 15 à son extrémité inférieure. Le train de tiges de forage 12 est mis en rotation par une table de rotation 16, entraînée par des moyens non représentés, qui met en prise une tige d'entraînement 17 à l'extrémité supérieure du train de tiges de forage. Le train de tiges de forage 12 est suspendu à un crochet 18, attaché à une moufle mobile
<Desc/Clms Page number 24>
(également non représentée), par l'intermédiaire de la tige entraînement 17 et une tête d'injection rotative 19 qui permet la rotation du train de tiges de forage par rapport au crochet.
Un liquide ou boue de forage 26est stocké dans un bassin de décantation 27 formé au site du puits. Une pompe 29délivre la boue de forage 26à l'intérieur du train de tiges de forage 12 par l'intermédiaire d'un orifice dans la tête d'injection 19, ce qui amène la boue de forage à s'écouler vers le bas dans le train de tiges de forage 12 comme indiqué par la flèche directionnelle 9. La boue de forage sort du train de tiges de forage 12 par des orifices dans le trépan 15, puis circule vers le haut à travers la région entre l'extérieur du train de tiges de forage et la paroi du sondage, appelée espace annulaire, comme indiqué par les flèches de direction 32. De cette manière, la boue de forage lubrifie le trépan 15 et entraîne les découpes de formation à la surface lorsqu'elle retourne au bassin de décantation 27 pour recirculation.
Le train de tiges de forage 12 comprend, de plus, un ensemble de fond de puits, généralement désigné par 34, près du trépan 15 (en d'autres termes, à plusieurs longueurs de masse-tige du trépan). L'ensemble de fond de puits comprend des dispositifs pour mesurer, traiter et stocker des informations, ainsi que communiquer avec la surface. L'ensemble de fond de puits 34 comprend donc, entre autres choses, un appareil de mesure et de communication locale 36pour déterminer et communiquer la résistivité de la formation F entourant le sondage 11. L'appareil de communication 36, également appelé
<Desc/Clms Page number 25>
outil de résistivité, comprend une première paire d'antennes émettrice/ réceptrice T, R, ainsi qu'une seconde paire d'antennes émettrice/ réceptrice T', R'.
La seconde paire d'antennes T', R' est symétrique par rapport à la première paire d'antennes T, R, comme décrit de manière plus détaillée ci-dessous. L'outil de résistivité 36 comprend, de plus un dispositif de commande pour commander l'acquisition de données, comme il est connu dans l'art.
L'ensemble de fond de puits 34 comprend, de plus, des instruments logés dans des masses-tiges 38,39 pour remplir différentes autres fonctions de mesure, telles que la mesure de la radiation naturelle, de la densité (rayons gamma ou neutrons), et la pression interstitielle de la formation F. Au moins certaines masses-tiges sont équipées de stabilisateurs 37, comme il est connu dans l'art.
Un sous-ensemble de communication avec la surface/locale 40 est également inclus dans l'ensemble de fond de puits 34, juste au-dessus de la masse-tige 39. Le sous-ensemble 40 comprend une antenne toroïdale 42 utilisée pour la communication locale avec l'outil de résistivité 36 (bien que d'autres moyens de communication locale puissent être avantageusement utilisés), et un type connu de système de télémétrie acoustique qui communique avec un système similaire (non représenté) à la surface terrestre par l'intermédiaire de signaux transmis dans le liquide ou boue de forage. Par conséquent, le système de télémétrie dans le sous-ensemble 40 comprend un émetteur acoustique qui génère un signal acoustique
<Desc/Clms Page number 26>
dans la boue de forage (également appelé "boue pulsée") qui est représentatif des paramètres de fond mesurés.
Le signal acoustique généré est reçu à la surface par des capteurs représentés par le numéro de référence 31. Les capteurs, par exemple, des capteurs piézoélectriques, convertissent les signaux acoustiques reçus en signaux électroniques. La sortie des capteurs 31 est couplée à un sous-système de réception de surface 90, qui démodule les signaux transmis. La sortie du sous-système de réception 90 est ensuite couplée à un processeur 85 et un enregistreur 45 informatiques. Le processeur 85 peut être utilisé pour déterminer le profil de résistivité de formation (entre autres choses) sur une base en "temps réel" durant la diagraphie ou ultérieurement en accédant aux données enregistrées provenant de l'enregistreur 45. Le processeur informatique est couplé à un moniteur 92 qui utilise une interface d'utilisateur graphique ("GUI") par l'intermédiaire de laquelle les paramètres de fond mesurés et des résultats particuliers dérivés de ceuxci (par exemple, des profils de résistivité) sont représentés graphiquement pour l'utilisateur.
Un système de transmission de surface 95 est également disposé pour recevoir des commandes entrées par l'utilisateur (par exemple, par l'intermédiaire de la GUI dans le moniteur 92) et est fonctionnel pour interrompre sélectivement le fonctionnement de la pompe 29d'une manière qui est détectable par des capteurs 99 dans le sous-ensemble 40. De cette manière, il existe une communication bilatérale entre le sous-ensemble 40 et l'équipement de surface. Un sous-ensemble 40 adapté
<Desc/Clms Page number 27>
est décrit de manière plus détaillée dans les brevets U. S. n 5 235 285 et 5 517 464, tous deux étant cédés au cessionnaire de la présente invention. il apparaîtra à l'homme du métier que d'autres techniques acoustiques, ainsi que d'autres moyens de télémétrie (par exemple, électromécaniques, électromagnétiques) peuvent être utilisés pour la communication avec la surface.
Dépendance azimutale de la mesure directionnelle et du nouvel équivalent de type à propagation
Deux types d'antenne à bobine sont utilisés pour comparer des mesures avec une sensibilité directionnelle. Un type obtient sa sensibilité directionnelle par le fait que l'antenne est décalée, par exemple, par rapport au centre de l'axe longitudinal d'un outil de diagraphie, ou partiellement couverte. Des mesures directionnelles peuvent également être effectuées avec une antenne configurée de telle manière que son moment magnétique ne soit pas aligné avec l'axe longitudinal de l'outil portant l'antenne.
La présente invention concerne le second type d'antenne à sensibilité directionnelle.
La figure 4 représente schématiquement un outil de résistivité de base 36 pour la mesure d'onde électromagnétique (EM) directionnelle. L'outil 36 comprend une antenne émettrice T qui déclenche une onde EM d'une certaine fréquence f et une antenne réceptrice R qui est éloignée d'une certaine distance L. Il est également décrit la paire symétrique (T', R') conformément aux enseignements de la publication de
<Desc/Clms Page number 28>
demande de brevet U.S. n 2003/0085707 ("Minerbo et al.") cédée au cessionnaire de la présente invention.
Par souci de clarté et de simplification, la discussion faite ci-après est limitée à l'antenne émettrice T et l'antenne réceptrice R, bien qu'elle soit généralement applicable à la paire d'antennes symétriques, T' et R'.
Il doit être noté que, bien que les moments inclinés des deux paires de symétrisation soient dans le même plan sur la figure 4, cela n'est pas requis dans la présente invention. Comme il apparaîtra dans la discussion ci-après, des signaux provenant de deux paires qui ont leurs moments dans des plans différents peuvent encore être additionnés les uns aux autres pour obtenir des résultats équivalents si les coefficients ou le déphasage ou l'atténuation directionnels extraits sont utilisés dans l'opération de symétrisation.
En fonctionnement, l'antenne réceptrice R enregistre une tension VRT induite par l'onde EM provenant de la première antenne émettrice T et ses courants secondaires générés dans la formation dans laquelle pénètre le sondage contenant l'outil de diagraphie 36. Les deux antennes T et R sont fixées sur l'outil 36 et, par conséquent, tournent avec l'outil.
Il peut être supposé que les orientations d'antenne peuvent former des angles #T pour l'antenne émettrice T, et #R pour l'antenne réceptrice R. la variation azimutale de la tension de couplage lorsque l'outil est en rotation peut ensuite être exprimée en terme de couplage de composants cartésiens des dipôles magnétiques par :
<Desc/Clms Page number 29>
VRT (#) = [V@cos#T cos#R+1/2(Vxx+V@y)sin#Tsin#R] + [Vxxsin#Tcos#R+Vxxcos#Tsin#R]cos#+[Vjesin#Tcos#R+V@@cos#Tsin#R]sin # + [1/2(Vyz+Vxy)sin#Tsin#R]sin2#+[1/2(V@x-Vy@)sin#Tsin#R]son2# = C@ (#T,#R)+Cic(#T,#R)cos#+Ci@(#T,#R)sin#+C@@(#T,#R)cos2#+C2,(#T,#R)sin2# (1.1) où un ensemble de coefficients complexes Co, C1c, C1s, C2c, C2s, a été défini pour représenter les amplitudes des différents composants de la réponse de formation mesurée. Les coefficients complexes sont donc définis par .
C@ (#T,#R) # [Vncos#Tcos#R+1/2(Vxx+Vyy)sin#Tsin#R] C1c (#T,#R) # [Vxxsin#Tcos#R+Vxxcos#Tsin#R] C@@ (#T,#R) # [Vy@sin#Tcos#R+V@ycos#Tsin#R] C2c (#T,#R) # [1/2(Vxx-Vyy)sin#Tsin#R] C2@ (#T,#R) # [1/2(Vyz+Vxy)sin#Tsin#R] (1.2)
Selon un aspect de la présente invention, il est reconnu que ces coefficients sont des fonctions de la résistivité de formation, la déviation du sondage et l'angle azimut à l'emplacement de l'outil.
Avec une opération de symétrisation, c'est-à-dire (#T # OR), l'éq. (1. 1) est simplifiée en : # (#) # VRT(#,#T,#R)-VRT(#,#R,#T) = 2 [V@@-V2x]sin(#T-#r)cos#+2[Vya-Vxy]sin(#T-#R)sin# #1c(#@ ,#R)cos#+#1s(#T,#R)sin# (1.3)
<Desc/Clms Page number 30>
Toutes les harmoniques de deuxième ordre (C2c, C2s) disparaissent après la soustraction parce qu'elles sont symétriques en terme d'échange d'angles d'inclinaison d'émetteur et de récepteur. Par conséquent, une antisymétrisation simplifie la variation azimutale du signal anti-symétrisé.
A ce stade, le point de référence de l'angle azimutal est arbitraire. Pour La planification de géométrie, si nous choisissons un point de référence à l'angle en tant que direction normale au plan de stratification, alors Vyz = Vzy = 0 par symétrie et # (#) aurait une dépendance pure vis à vis de cos#. Dans une application pratique, l'orientation de la stratification est inconnue. Cependant, quelle que soit la référence, l'orientation de stratification peut être calculée par :
Figure img00300001
Avec la rotation de #bed, # est normal à la stratification et donc, # (#) est exactement [Vxz - Vzx] hormis une constante de multiplication 2sin(OT - OR) -
Une fois que la tension à chacune des bobines de récepteur due à chacune des bobines d'émetteur est déterminée, la mesure totale peut être déterminée en additionnant les tensions dans le cas d'un outil à induction ; ou en calculant le rapport complexe des tensions dans le cas d'un outil à propagation. Par
<Desc/Clms Page number 31>
exemple, pour le dispositif de diagraphie à propagation de la figure 4, la valeur absolue de la tension à chaque récepteur peut être obtenue comme étant la racine carrée de la somme des carrés des parties réelle et imaginaire de la tension complexe (éq. 1. 1) et le rapport des valeurs absolues donne l'atténuation, à partir de laquelle la résistivité déterminée à partir de l'atténuation Rad peut être obtenue (résistivité de formations à une profondeur relativement élevée d'investigation autour des récepteurs). La phase pour chaque récepteur est obtenue à partir de l'arctangente du rapport des parties réelle et imaginaire de la tension complexe et le déphasage est la différence de phase aux deux récepteurs. La résistivité déterminée à partir du déphasage Rps peut ensuite être obtenue (résistivité de formations à une profondeur relativement faible d'investigation autour des récepteurs).
Pour les mesures de type à propagation, la différence de logarithme des tensions (ou le rapport) entre deux mesures est utilisée. Selon les enseignements de Minerbo et al., nous utilisons l'amplitude de la réponse azimutale, c'est-à-dire, la différence de déphasage et d'atténuation de mesure, à un angle et celle à (# + 180), évaluées au maximum de la réponse de tension. Cela conduit approximativement, à partir des équations (1.1-2) :
<Desc/Clms Page number 32>
Figure img00320001
Le maximum de #V# est obtenu à # = 0 si x est choisi comme étant la direction normale à la stratification. Evaluée à l'angle 0 = 0, l'éq. (1.5) donne :
Figure img00320002
Cependant, ce n'est pas encore le type xz-zx pur de réponses qui est souhaité, c'est-à-dire, qui sont insensibles à l'anisotropie et l'angle de pendage de stratification.
La présente invention concerne des mesures directionnelles qui sont insensibles à l'anisotropie de la formation dans une large plage d'angles de pendage et dans une plage de fréquence large. Comme mentionné ci-dessus, des modes de réalisation particuliers de l'invention sont basés sur des configurations ou systèmes d'antennes anti-symétrisées. Après une procédure de symétrisation (#T # OR), comme décrit par Minerbo et al., nous obtenons :
<Desc/Clms Page number 33>
Figure img00330001
Cela est, à nouveau, similaire à la réponse de type à induction, bien que le dénominateur ait encore des composantes qui ne soient pas simplement [xz-zx]. Cela prouve que la procédure de symétrisation pour une mesure de type à propagation peut produire des réponses similaires à celles de type à induction symétrisée, mais pas un type pur. Il est également vrai qu'une mesure de propagation peut être effectuée dans deux orientations arbitraires de la réponse azimutale.
La réponse finale de cette analyse contient des composantes de différents types de couplages, ce qui, à première vue, semble indésirable compte tenu des techniques de diagraphie conventionnelles. Cependant, cette réponse permet d'obtenir un procédé amélioré de conduite de la mesure de propagation qui est simple et beaucoup plus proche du type à induction. La mesure directionnelle de l'invention est réalisée en observant que différents types de couplages se séparent naturellement en différentes dépendances azimutales. Il doit être noté que nous pouvons utiliser Co, C1c, Cis pour effectuer une mesure directionnelle plus pure.
Une orientation de stratification déterminée #bed peut être supposée à l'aide de l'éq. (1.4), de telle manière que tous les angles 0 et #,# soient référencés dans cette direction. Dans ce cas l'éq. (1. 1) est simplifiée en :
<Desc/Clms Page number 34>
Figure img00340001
Une mesure de propagation peut alors être simplement définie par :
Figure img00340002

et
Figure img00340003
Ensuite, avec la symétrisation, la première expression est obtenue sous la forme :
Figure img00340004

qui est seulement légèrement différente de l'éq. (1.7) .
<Desc/Clms Page number 35>
Par conséquent, pour le type à induction et à propagation de mesures à deux bobines (TR), l'analyse des différentes composantes doit être effectuée sur le niveau de tension VRT(). Cela produit les comportements exacts de sin, cos, sin2 et cos2 qui peuvent être extraits par l'intermédiaire d'un nouvel algorithme de traitement décrit dans la section suivante.
Il doit être noté, à nouveau, que dans des cas particuliers où #T = 0 ou OR = 0, les termes sin2 et cos2 disparaissent. La dépendance vis à vis de la tension est simplement basée sur les termes sin et cos#.
Un aspect important des mesures de type à propagation de déphasage et d'atténuation est qu'elles sont naturellement adaptées pour des mesures "en cours de forage", pour lesquelles une caractérisation détaillée de la dérive électronique thermique dans des conditions de fond est difficile à obtenir. Les mesures de déphasage et d'atténuation directionnelles présentement définies présentent l'avantage d'un outil de résistivité à propagation compensée de sondage conventionnel : la caractéristique d'antenne émettrice et réceptrice et la dérive de l'électronique de réception sont complètement éliminées de la mesure.
Il doit également être noté que dans cette invention, le processus de symétrisation est exécuté avec les coefficients extraits qui peuvent être extraits indépendamment de azimut relatif effectif des deux paires TR. Si l'orientation de paire TR est physiquement tournée à un angle supplémentaire o
<Desc/Clms Page number 36>
autour de l'axe d'outil, la réponse est décrite exactement avec la même expression sauf que est remplacé par # # # + o. Cela a des implications importantes pour les configurations de bobine effectives. En utilisant deux paires de mesures TR, mais avec des orientations azimutales différentes, le processus de symétrisation durant la rotation peut être exécuté exactement comme si les deux paires étaient dans le même plan. Cependant, lorsque l'outil coulisse, le signal mesuré dans les deux orientations peut être utilisé pour effectuer les mesures directionnelles requises en supposant simplement que azimut de l'outil et la stratification ne changent pas par rapport à la dernière fois que l'outil était encore en rotation. Par conséquent, nous pouvons utiliser l'éq. 1. 3 pour obtenir les amplitudes Clc et Cls requises pour effectuer les mesures de propagation.
Ces analyses peuvent être étendues directement au type TRR conventionnel de mesures, comme décrit dans Minerbo et al. Il apparaîtra aisément à l'homme du métier que cette procédure produit essentiellement la même réponse que celle décrite ci-dessus, mais avec deux fois le signal lorsque l'espacement entre la paire de récepteurs est beaucoup plus faible comparé à l'espacement TR. Les signaux directionnels des deux récepteurs sont simplement additionnés.
La figure 5A représente une configuration TRR qui est insensible à l'anisotropie à un angle de pendage quelconque, et la figure 5B représente des réponses selon cette configuration. L'antenne émettrice Tl est excitée et le déphasage et l'atténuation des antennes
<Desc/Clms Page number 37>
réceptrices Rll, R12 sont mesurés. Ensuite, l'antenne émettrice T2 est excitée et le déphasage et l'atténuation des antennes réceptrices R21, R22 sont mesurés. La lecture de l'outil correspond aux différences entre ces deux ensembles de mesures. Etant donné que les mesures individuelles sont identiques dans un milieu homogène à un angle quelconque et avec une anisotropie quelconque, la lecture de l'outil est zéro dans un milieu homogène à un pendage quelconque.
Les réponses de mesure dans une formation anisotrope à trois couches sont décrites sur la figure 5B. La lecture de l'outil est zéro loin de la frontière à un angle de pendage quelconque, et il y a peu de sensibilité à l'anisotropie à proximité de la frontière. La séparation des réponses provient du fait que les réponses de propagation ne sont pas symétriques si les emplacements de l'émetteur et du récepteur sont échangés. Une mesure haut/bas contient uniquement les informations directionnelles, même à proximité de la frontière. Il doit être observé que les réponses d'atténuation sont pratiquement chevauchantes pour différents pendages si toutes les antennes sont dans le même milieu, de manière similaire à la mesure d'induction XZ-ZX idéale (décrite ci-dessus). Les mesures de déphasage sont également chevauchantes, bien que les réponses aient une valeur doublée dans la couche conductrice (1 S/m).
<Desc/Clms Page number 38>
Algorithmes de processeur de signal numérique (DSP) pour extraire le signal de la dépendance azimutale
Il apparaît clairement dans l'analyse précédente que le meilleur moyen d'extraire les mesures directionnelles est d'analyser les signaux de tension et extraire les mesures appropriées, puis de les synthétiser par symétrisation pour utilisation finale dans le géoguidage. Le procédé conventionnel pour traiter les données azimutales consiste à les mettre en fichiers dans des petits fichiers de portée azimutale égale (voir la partie gauche de la figure 6), comme cela est effectué avec des images de sondage de densité et de résistivité. Ensuite, les fichiers haut et bas peuvent être identifiés en localisant le maximum et le minimum. Cependant, il existe de nombreux inconvénients dans cette technique :
1. la mise en fichiers réduit la résolution d'angle de la mesure ;
2. la mise en fichiers conduit à la perte de données étant donné qu'aucune des données hors des fichiers haut et bas n'est utilisée pour le calcul de distance à la frontière ;
3. dans des cas de broutage, la mise en fichiers n'est pas uniforme et donc, les fichiers haut et bas peuvent être vides ou avoir des petits échantillons, induisant une erreur dans les mesures haut/bas et produisant éventuellement l'identification erronée des valeurs de crête dans le pire des cas ;
4. la mise en fichiers consomme plus de mémoire.
Un aspect de la présente invention concerne principalement la détermination des coefficients des
<Desc/Clms Page number 39>
termes sin#, cos#, sin2 et cos2, qui définissent un nouvel algorithme à utiliser au lieu de la mise en fichiers conventionnelle. Cette technique de l'invention est appelée ajustement "à la volée" de la réponse azimutale en fonction des extractions des termes sin et con correspondants des mesures directionnelles, effectuées de manière itérative (voir les points représentés graphiquement sur la figure 9). Un tel algorithme d'ajustement est exécuté dans un DSP par l'intermédiaire d'un algorithme entier, par conséquent, il est suffisamment rapide pour être exécuté pour tous les canaux dans le temps d'échantillonnage de 4 ms. L'utilisation précise des informations d'angle azimut, et la randomisation des séquences d'acquisition, rendent l'algorithme robuste pour tolérer une rotation d'outil irrégulière ainsi qu'un broutage dans des conditions de forage rigoureuses. De cette manière, toutes les données sont utilisées pour obtenir le signal haut/bas au lieu d'uniquement les données dans les deux fichiers, ce qui améliore le rapport signal/bruit dans la mesure.
L'utilisation d'angles azimut précis rend également l'orientation de stratification déterminée plus précise.
L'algorithme détaillé peut être décrit comme suit.
Implémentation de virgule flottante : commencer par une valeur initiale de matrice Po et de vecteur Uo, puis passer à l'algorithme décrit ci-dessous (également représenté graphiquement dans la partie droite de la figure 6) avec la mesure y(#1) et la base r = (1 cos#1 sin#1 cos2#1 sin2i)T, où P est une matrice de dimension
<Desc/Clms Page number 40>
M x M et U et r sont des vecteurs de dimension M. M est la dimension de la fonction de base. Après N itérations, U converge vers une valeur qui représente les coefficients de l'expression. Cet algorithme est stable et la convergence est généralement obtenue en 10 à 15 itérations.
L'algorithme détaillé est décrit ci-dessous : initialiser Po et Uo ; pour m = 1 à Nsamples
Figure img00400001

m suivant ; retour (U) ; où :
Nsamples est le nombre totale d'échantillons acquis dans un cycle,
M est la dimension du vecteur de fonction d'approximation (nombre de fonctions d'approximation),
U est le vecteur de coefficients d'ajustement de dimension M, r est le vecteur de valeurs de fonction d'approximation à chaque position de mesure de dimension M et
P est une matrice de dimension M x M.
<Desc/Clms Page number 41>
Algorithmes entiers dans le DSP
Dans de nombreux cas, l'implémentation de virgule flottante est trop coûteuse pour être effectuée avec les unités centrales de fond actuellement disponibles parce qu'il faut ajuster des centaines de canaux et l'acquisition de données pour chaque angle azimut doit être très courte (ms) afin que l'angle soit précis à une vitesse de rotation élevée. Dans ce cas, une implémentation d'entier peut être appliquée, avec une certaine modification, pour améliorer la précision (par exemple, utiliser 32 bits pour la multiplication), effectuer la remise à l'échelle afin d'éviter un dépassement de capacité et pour accélérer la convergence. Les valeurs de la fonction de base peuvent également être prégénérées et stockées en mémoire de manière à être interpolées ultérieurement pour obtenir la valeur de l'angle réel #1. En résumé, les algorithmes entiers sont : - basés sur une représentation de données entières 16 bits avec des variables intermédiaires entières 32 bits ; - plus efficaces mais moins précis ; - adaptatifs à des stratégies pour améliorer la précision et la vitesse de convergence : mise à l'échelle, initialisation et réinitialisation ; - adaptatifs à des stratégies pour adapter l'environnement de DSP : division, arrondi.
Un flux de convergence de la version entière de l'algorithme est décrit sur la figure 7. Il doit être noté que l'ajustement est très précis avec des erreurs généralement inférieures à 1 %. Un point très important
<Desc/Clms Page number 42>
de l'ajustement est que pratiquement toutes les données sont utilisées pour dériver les coefficients (signaux de mesure directionnelle), ce qui améliore significativement le rapport signal/bruit. Par exemple, si 32 fichiers sont utilisés, seulement 1/16 des données sont utilisées si la mise en fichiers haut/bas est mise en oeuvre. En utilisant l'ajustement à la volée, par contre, pratiquement toutes les données sont prises en compte (hormis la partie de convergence initiale).
Etant donné que seulement les signaux appropriés sont extraits dans la technique d'ajustement (décrite ci-dessus), seuls les coefficients utiles doivent être enregistrés. Par conséquent, dans ce cas, il n'est nécessaire d'enregistrer que 5 coefficients, comparés à 32 si on met en fichier toutes les données en utilisant l'exemple de 32 fichiers. Il apparaîtra à l'homme du métier les avantages de la technique de l'invention qui comprennent la précision du signal extrait et une amélioration particulière de la précision de l'angle azimutal.
Dans un autre aspect, la présente invention propose la caractérisation du bruit dans des mesures directionnelles. A savoir, dans des modes de réalisation dans lesquels seulement l'une parmi la première paire d'antennes (TR) a un dipôle magnétique incliné, le bruit des signaux de tension mesurés peut être caractérisé en utilisant les coefficients de seconde harmonique. Dans des modes de réalisation dans lesquels chacune des premières antennes a un dipôle magnétique incliné ou transversal, et les coefficients
<Desc/Clms Page number 43>
d'ajustement comprennent les coefficients de troisième harmonique. Dans d'autres cas supplémentaires, le bruit des signaux de tension mesurés peut être caractérisé en combinant les signaux.
La présente invention est également adaptative à des conditions dans lesquelles l'instrument de diagraphie n'est pas mis en rotation, telles que, par exemple, lors de l'arrêt de la rotation du train de tiges de forage lors d'un forage directionnel avec un ensemble de moteur à boue. Dans de telles conditions, la seconde paire d'antennes (T'R') est, de préférence, symétrique par rapport à la première paire d'antennes.
L'azimut d'une couche de formation souhaitée est déterminé en combinant ces couplages d'antenne (comme décrit ci-dessus) et en déterminant les coefficients constant et de première harmonique à partir des signaux de tension mesurés acquis lorsque l'instrument n'est pas en rotation. Les coefficients déterminés peuvent ensuite être utilisés pour exécuter l'ajustement itératif lorsque l'outil de diagraphie n'est pas à nouveau en rotation.
Distance de frontière pour utilisation en géoguidage
Selon un autre aspect de la présente invention, les coefficients déterminés peuvent être utilisés pour obtenir l'orientation de la stratification de formation. Pour chaque canal de mesures directionnelles, l'orientation de la stratification peut être déterminée à l'aide de l'éq. 1. 4. Il peut être noté, de plus, que la valeur moyenne de nombreux canaux, pondérée par une fonction de l'intensité de
<Desc/Clms Page number 44>
signal relative de chaque canal qui comprend une mesure, peut également être utilisée pour augmenter la précision étant donné que cette orientation doit être identique pour tous les canaux. L'amplitude et la phase du signal de tension VRT(#bed) dans la direction normale supposée par rapport à la couche peuvent ensuite être calculés.
Le déphasage et l'atténuation directionnels peuvent être calculés avec l'éq. 1. 10. Une symétrisation est ensuite effectuée pour générer le signal final nécessaire pour obtenir la distance à la frontière. Il apparaîtra, cependant, que l'ordre de ces étapes peut être modifié avec des réponses similaires ou identiques.
Pour obtenir la distance de frontière en temps réel, deux techniques peuvent être utilisées. Pour des modèles simples (une seule frontière), un tracé croisé des deux mesures directionnelles nous permet d'obtenir la distance à la frontière et la résistivité de formation d'une des deux couches. Un graphique de tracé croisé représentatif est décrit sur la figure 8, en utilisant la réponse d'une atténuation à 84 pouces, 100 kHz en fonction d'une résistivité de déphasage à 28 pouces, 2 MHz (la résistivité d'épaulement étant Rh = 0,8 #m, Rv = 3,6 #m). Dans ce cas, la symétrisation élimine efficacement l'influence d'autres paramètres tels que l'anisotropie et le pendage. Les graphiques de mesure symétrisée sont simples.
L'utilisation de différentes combinaisons de mesures appariées pour obtenir une image cohérente augmente le niveau de confiance de l'interprétation.
<Desc/Clms Page number 45>
La figure 9 illustre l'utilisation d'une inversion à base de tracé croisé pour interpréter les mesures directionnelles. Deux mesures directionnelles à 84 pouces, 400 kHz sont utilisées pour interpréter les mesures en utilisant un modèle à une frontière avec une résistivité de couche fixe de 100 Qm et une résistivité (au-dessous de l'outil) et une distance de couche d'épaulement variables. La réponse n'est pas sensible à la résistivité de couche et le signal est principalement déterminé par la résistivité et la distance de couche d'épaulement. Les points représentés sur l'image d'écran capturée représentent des mesures de résistivité de couche d'épaulement et des lectures de distance. Les valeurs de résistivité et de distance sont lues sur le tracé croisé et transférées à l'écran.
Pour le cas d'une frontière à une couche dans lequel les résistivités des couches des deux côtés de la frontière sont connues, on peut utiliser trois entrées, dont au moins une est une mesure de la résistivité autour de l'outil (par exemple, l'outil 36 sur la figure 3).
Pour des modèles plus complexes mettant en oeuvre plusieurs frontières, un programme d'inversion basé sur la technique décrite du brevet U.S. n 6 594 584 est utilisé. La trajectoire du sondage est projetée sur un modèle initial de la formation. Un outil de diagraphie est disposé dans un segment de la trajectoire de sondage et ses réponses le long du segment sont mesurées. Les réponses prévues de l'outil, selon le modèle, sont également déterminées. Les différences résultant de la comparaison des réponses prévues et
<Desc/Clms Page number 46>
mesurées le long du segment sont ensuite utilisées pour ajuster le modèle, et le cycle de comparaison et d'ajustement est répété jusqu'à ce que les différences diminuent au-dessous d'un seuil choisi. La robustesse du modèle d'inversion est améliorée avec des points de départ et des critères physiques multiples pour distinguer les solutions.
Un algorithme d'inversion est ensuite appliqué.
Cet algorithme accepte n'importe quelle mesure en entrée et détermine ensuite le modèle le plus cohérent pour les données. Des modèles multiples sont appliqués et le meilleur modèle est automatiquement sélectionné.
Une reconstruction de profil de résistivité type est décrite sur la figure 10A, tandis que la reconstruction de frontière (structure) correspondante est décrite sur la figure 10B. Ces résultats indiquent que les solutions sont plus précises à l'approche de couches conductrices qu'à l'approche de couches résistives.
Cela est prévisible parce que les réponses directionnelles sont principalement définies par la couche conductrice et ont une sensibilité nettement supérieure à une variation de résistivité dans les couches conductrices au-dessus et au-dessous, puis à la variation de résistivité dans des couches d'épaulement résistives.
Un autre aspect de la présente invention est l'utilisation des distances de frontière déterminées pour prendre des décisions de géoguidage.
L'organigramme de géoguidage global peut être décrit en référence à la figure 11. Premièrement, des mesures directionnelles en temps réel sont sélectionnées dans
<Desc/Clms Page number 47>
la case 110, en utilisant la technique d'ajustement à la volée comme décrit ci-dessus. Les mesures directionnelles sont ensuite visualisées, dans la case 120 en affichant les réponses de diagraphie brutes, ou en utilisant des graphiques de tracé croisé comme décrit sur les figures 8 et 10 (et présentement décrites par ailleurs).
Le point de décision 130 dirige ensuite l'exécution, selon qu'un modèle approprié ait été identifié ou non. Si l'agencement de formation ou d'autres paramètres sont connus, cette approche permet à l'utilisateur de fixer les paramètres connus. Par exemple, la résistivité de couche d'épaulement supérieure ou la stratification peuvent être connues comme étant stables, et il peut être connu avec une certitude élevée qu'il n'y a pas de frontière audessous. De telles informations permettent, par exemple, la sélection de modèle avec une seule frontière et une résistivité de couche d'épaulement fixe. Mathématiquement, cela signifie que seuls trois paramètres doivent être inversés, la résistivité horizontale et verticale de la couche (Rh et Rv) et la distance à la frontière supérieure. C'est un exemple d'inversion de modèle complète (case 140), qui assure une interprétation cohérente et évite toute confusion dans des cas limites dans lesquels un modèle non physique avec un meilleur ajustement mathématique de mesure peut devoir être sélectionné.
En variante, des inversions par approximation rapides avec des modèles multiples sont exécutées (case 150) ainsi qu'un algorithme de sélection de modèle
<Desc/Clms Page number 48>
automatique (case 160). L'algorithme applique différents modèles, du plus simple (pas de frontière, formation isotrope) au plus complexe (deux distances et formation anisotrope), comprenant les suivants : isotrope homogène (un seul paramètre : résistivité) ; anisotrope homogène (deux paramètres : Rh et Rv) ; formation isotrope à frontière unique, frontière au-dessus ou au-dessous (trois paramètres : Rcouche, Répaulement et distance à la frontière) ; formation anisotrope à frontière unique, frontière au-dessus ou au-dessous (quatre paramètres : Rcouche~h, Rcouche v, Répaulement et distance à la frontière) ; formation isotrope à deux frontières et trois paramètres : (cinq paramètres : Rcouche, Répaulement~haut, Répaulement~bas et distance à la frontière au-dessus et au-dessous de l'outil) ; formation anisotrope à deux frontières et trois paramètres : (six paramètres : Rcouche~h, Rcouche~v, Répaulement~haut, Répaulement~bas et distance à la frontière au-dessus et au-dessous de l'outil).
Une solution est obtenue pour chaque modèle.
L'algorithme de sélection de modèle utilise des contraintes physiques pour les mesures directionnelles et impose des conditions qui sélectionnent le "modèle le plus simple le mieux corrélé aux données" (règle de Bayes). Le critère d'information de Akaike (AIC) conventionnel ou le critère d'information bayesien peuvent être utilisés pour pénaliser la complexité du modèle.
<Desc/Clms Page number 49>
Une validation de modèle est ensuite appliquée, au point de décision 180, pour déterminer si le modèle choisi est cohérent avec les connaissances antérieures sur la formation géologique ou d'autres mesures telles que les rayons gamma, ou d'autres mesures disponibles en temps réel comme décrit par la case 170. Si le modèle est cohérent avec les autres données, il est accepté en entrée (case 190) pour la décision de forage (case 200), et un rapport approprié est généré à la case 210. La mise en oeuvre de la décision de forage conduit à des mesures directionnelles en temps réel supplémentaires (case 100), qui sont à nouveau entrées dans les étapes d'ajustement et de visualisation des cases 110 et 120. Si le modèle d'inversion sélectionné est incohérent avec les autres mesures de la case 170, une mise à jour du modèle est appropriée. Dans ce cas, un sous-programme en mode interactif 220 est appliqué qui est cohérent avec le programme d'inversion complexe du brevet US 6 594 584.
Souvent, l'inversion automatique peut créer une interprétation qui n'est pas cohérente. Elle peut être causée par le bruit de mesure causé par l'électronique ainsi que le "bruit de modèle", c'est-à-dire, le fait que le modèle réel est différent de tous les modèles utilisés dans l'inversion. Par exemple, le sondage, l'invasion, les couches minces, la stratification croisée et la taille d'outil finie ne sont pas inclus dans les modèles et ils peuvent causer une erreur dans le processus d'ajustement.
La flexibilité et la capacité d'un utilisateur à sélectionner de manière interactive le modèle commun
<Desc/Clms Page number 50>
pour un segment de données sont essentielles pour l'interprétation correcte des mesures. Le logiciel définissant le sous-programme en mode interactif (case 220) a les fonctionnalités suivantes : affinement du modèle (case 221) sur la base des tendances, connaissances antérieures, ou une source d'informations externe ; limitation ou fixation de certains paramètres d'inversion (case 222) ; élimination de certaines mesures qui peuvent être plus affectées par l'environnement non inclus dans le modèle ; et retraitement des données (case 223).
La figure 12 illustre un organigramme général pour générer des tracés croisés de mesures directionnelles qui fournissent des distances de frontière, comme présentement décrit par ailleurs. La génération ou la modification d'un tracé croisé est représentée par le sous-programme 20, qui est initié en définissant un modèle approprié (case 30). Les mesures directionnelles appropriées sont ensuite sélectionnées en fonction de la capacité à déterminer/prédire la résistivité d'épaulement (case 40). Si la résistivité de couche d'épaulement est connue de manière fiable comme c'est souvent le cas dans les grands champs dans lesquels de nombreux puits sont forés, l'interprétation est basée sur une résistivité et une mesure directionnelle (case 50), pour déterminer la résistivité de couche réelle (couche d'épaulement corrigée) et la distance à la frontière. Si la résistivité de couche d'épaulement n'est pas connue, il est recommandé d'utiliser le tracé croisé de deux mesures directionnelles (case 60), comme décrit sur la figure 9, pour déterminer la résistivité
<Desc/Clms Page number 51>
de couche d'épaulement et la distance à la frontière.
Les mesures sélectionnées sont entrées dans le modèle défini pour générer des réponses, à la case 70, qui peuvent être visualisées (case 80) comme indiqué par les figures 8 et 9. Dans les cas où aucune résistivité n'est connue, mais les plages de leurs différences sont connues, il est possible de générer des tracés croisés tridimensionnels, en combinant une résistivité et deux mesures directionnelles.
Une fois que le tracé croisé a été créé ou mis à jour (sous-programme 20), il peut être mis à jour en continu avec des mesures en temps réel supplémentaires (case 100) et visualisé (case 120'). Le traitement du tracé croisé de mesures directionnelles donne la distance à une ou plusieurs frontières de formation et la résistivité d'une ou plusieurs couches (case 220), qui peuvent être transmises et visualisées (case 230) pour permettre la décision de forage, à la case 200.
Selon un autre aspect de la présente invention, une interface d'utilisateur graphique (GUI) a été créée pour faciliter l'utilisation du programme et afficher les modèles de stratification inversés pour visualisation par l'utilisateur. Un schéma de principe illustrant la structure de la GUI et différentes manières de traiter les données sont représentés sur la figure 13. A savoir, la GUi permet : - l'application de modèle direct pour une stratification donnée, - la génération de tracés croisés pour l'analyse de sensibilité de mesure,
<Desc/Clms Page number 52>
- l'inversion de données en temps réel caractérisée par : *l'inversion à base de tracé croisé d'une résistivité et une distance (pour une résistivité d'épaulement ou de couche donnée), *l'interprétation automatique, en utilisant une inversion rapide basée sur la superposition de réponses de frontière unique et la sélection de modèle automatique, des contraintes physiques combinées avec le critère d'information de Akaike ou un modèle mixte, *l'interprétation interactive basée sur l'inversion complète de paramètres quelconques (jusqu'à 6) d'un modèle à trois couches (une couche et deux épaulements - deux distances et quatre résistivités), l'application d'un modèle de milieu stratifié complet dans la boucle d'inversion, permettant des affinements de modèle, avec des options pour sélectionner/désélectionner certains paramètres d'inversion et limiter ceux-ci et/ou sélectionner et repondérer les mesures disponibles.
*le basculement entre le mode rapide et interactif pour l'affinement de modèle.
- la visualisation de mesures et d'inversion.
Un exemple d'un écran d'interprétation est décrit sur la figure 14. les mesures qui sont utilisées dans l'inversion et la trajectoire de sondage sont affichées à l'écran. Les résultats inversés, comprenant la distance aux frontières supérieure et inférieure, les résistivités de couche et deux couches d'épaulement sont représentés graphiquement. Les résultats sont
<Desc/Clms Page number 53>
représentés point par point, et basés sur des positions reconstruites des frontières. Il est possible de déduire le pendage stratigraphique. Il doit être noté que, bien que le modèle direct appliqué dans la boucle d'inversion soit un milieu stratifié 1D, l'approche permet la génération d'images de formation plus complexes, comprenant : - des frontières non parallèles, - la détection de défauts, comprenant des défauts sous-sismiques.
Evaluation de l'anisotropie
Avec une configuration TRR dans laquelle l'antenne émettrice et au moins une des antennes réceptrices sont inclinées, des mesures d'anisotropie peuvent être effectuées, conformément aux enseignements de la publication de demande de brevet U.S. n 2003/020029 ("Omeragic et al."). Les signaux EM détectés relatifs à un déphasage ou un rapport d'amplitude sont combinés pour déterminer l'anisotropie. Il doit être noté qu'une nouvelle technique de compensation de sondage est décrite par Omeragic et al. pour une telle mesure.
La partie directionnelle de telles mesures combinées peut être analysée de la même manière que celle décrite ci-dessus. En conséquence, les coefficients de composante peuvent être utilisés pour dériver l'anisotropie par inversion de la technique décrite dans le brevet U.S. n 6 594 584 (également décrite ci-dessus). Cela est vrai pour tous les angles de pendage.
<Desc/Clms Page number 54>
Il apparaîtra à l'homme du métier que dans des sondages fortement déviés, la mesure de propagation directionnelle avec des bobines axiales présente déjà une sensibilité à l'anisotropie acceptable, ce qui élimine le besoin de la technique de compensation de sondage dans ces sondages.
Il apparaîtra à l'homme du métier que cette invention peut être mise en oeuvre en utilisant un ou plusieurs ordinateurs d'usage général adaptés ayant un matériel approprié et programmés pour exécuter les processus de l'invention. La programmation peut être effectuée en utilisant un ou plusieurs dispositifs de stockage lisibles par le processeur informatique et codant un ou plusieurs programmes d'instructions exécutables par l'ordinateur pour effectuer les opérations décrites ci-dessus. Le dispositif de stockage de programme peut prendre la forme, par exemple, d'une ou plusieurs disquettes ; un CD-ROM ou un autre disque optique ; une bande magnétique ; une puce de mémoire morte (ROM) ; et d'autres formes du type connu dans l'art ou développées ultérieurement. Le programme d'instructions peut être en "code objet", c'est-à-dire, sous forme binaire qui est exécutable plus ou moins directement par l'ordinateur ; en "code source" qui requiert une compilation ou une interprétation avant exécution ; ou sous une forme intermédiaire telle qu'un code partiellement compilé.
Les formes précises du dispositif de stockage de programme et du codage d'instructions sont présentement immatérielles. Par conséquent, ces moyens de traitement peuvent être mis en oeuvre dans l'équipement de
<Desc/Clms Page number 55>
surface, dans l'outil ou partagés par les deux comme il est connu dans l'art. Il apparaîtra également que les techniques de l'invention peuvent être utilisées avec un type quelconque de système de diagraphie de forage, par exemple, des outils de travail au câble, des outils LWD/MWD ou des outils LWT.
Il apparaîtra à la lecture de la description cidessus que différents modifications et changements peuvent être effectués dans les modes de réalisation préférés et autres de la présente invention sans s'écarter de sa portée réelle.
Cette description est destinée à des fins d'illustration uniquement et ne doit pas être considérée dans un sens limitatif. La portée de cette invention doit être déterminée uniquement par le contenu des revendications ci-après. Le terme "comprenant" dans les revendications est censé signifier "comprenant au moins" de telle manière que la liste d'éléments mentionnée dans une revendication soit un groupe ouvert. "Un" et les autres termes singuliers sont destinés à inclure les formes plurielles de ceuxci, sauf exclusion spécifique.

Claims (9)

REVENDICATIONS
1. Procédé pour caractériser une formation souterraine, comprenant l'étape consistant à : générer un tracé croisé de deux mesures de diagraphie directionnelles obtenues à l'aide d'un instrument disposé dans un sondage croisant la formation pour obtenir la distance à au moins une frontière de formation et la résistivité pour au moins une couche de formation.
2. Procédé selon la revendication 1, dans lequel : le tracé croisé est généré en utilisant un modèle à une frontière ; la résistivité obtenue est la résistivité de couche d'épaulement ; et la distance obtenue est la distance la plus proche de la couche d'épaulement.
3. Procédé selon la revendication 1, dans lequel l'étape de génération de tracé croisé comprend les étapes consistant à : définir un modèle approprié ; sélectionner des mesures directionnelles appropriées ; entrer les mesures choisies dans le modèle défini pour générer le tracé croisé ; et générer une représentation graphique du tracé croisé.
<Desc/Clms Page number 57>
4. Procédé selon la revendication 1, dans lequel la distance et la résistivité obtenues sont utilisées pour prendre des décisions de forage.
5. Procédé pour caractériser une formation souterraine, comprenant l'étape consistant à : générer un tracé croisé d'une mesure de résistivité et directionnelle déterminée en utilisant un instrument disposé dans un sondage croisant la formation pour obtenir la distance à au moins une frontière de formation et la résistivité pour au moins une couche de formation.
6. Procédé pour caractériser une formation souterraine, comprenant l'étape consistant à : générer un tracé croisé d'une résistivité et deux mesures directionnelles déterminées en utilisant un instrument disposé dans un sondage croisant la formation pour obtenir la distance à au moins une frontière de formation et la résistivité pour au moins deux couches de formation.
7. Procédé selon la revendication 6, comprenant, de plus, les étapes consistant à : choisir un modèle d'inversion approprié pour les mesures directionnelles en temps réel choisies ; vérifier que le modèle choisi est cohérent avec d'autres informations ; et utiliser le modèle vérifié pour prendre des décisions de forage.
<Desc/Clms Page number 58>
8. Procédé selon la revendication 7, dans lequel l'étape de choix de modèle comprend la création d'une visualisation des mesures directionnelles choisies.
9. Procédé selon la revendication 7, dans lequel l'étape de choix de modèle comprend : l'identification de paramètres de formation connus ; le choix interactif des modèles avec lesquels effectuer l'invention des mesures directionnelles choisies ; et choisir le modèle le plus simple qui est ajusté aux informations connues;
FR0411139A 2003-05-22 2004-10-20 Procede de resistivite electromagnetique directionnelle Expired - Fee Related FR2859539B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US47268603P 2003-05-22 2003-05-22
US10/709,212 US7382135B2 (en) 2003-05-22 2004-04-21 Directional electromagnetic wave resistivity apparatus and method

Publications (2)

Publication Number Publication Date
FR2859539A1 true FR2859539A1 (fr) 2005-03-11
FR2859539B1 FR2859539B1 (fr) 2013-01-18

Family

ID=32511770

Family Applications (2)

Application Number Title Priority Date Filing Date
FR0451007A Expired - Fee Related FR2855272B1 (fr) 2003-05-22 2004-05-21 Appareil et procede de resistivite electromagnetique directionnelle
FR0411139A Expired - Fee Related FR2859539B1 (fr) 2003-05-22 2004-10-20 Procede de resistivite electromagnetique directionnelle

Family Applications Before (1)

Application Number Title Priority Date Filing Date
FR0451007A Expired - Fee Related FR2855272B1 (fr) 2003-05-22 2004-05-21 Appareil et procede de resistivite electromagnetique directionnelle

Country Status (9)

Country Link
US (1) US7382135B2 (fr)
CN (2) CN101082276B (fr)
BR (1) BRPI0401805B1 (fr)
DE (1) DE102004024969A1 (fr)
EA (1) EA007587B1 (fr)
FR (2) FR2855272B1 (fr)
GB (1) GB2402489B (fr)
MX (1) MXPA04004692A (fr)
NO (1) NO335727B1 (fr)

Families Citing this family (148)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6163155A (en) * 1999-01-28 2000-12-19 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations
US7659722B2 (en) 1999-01-28 2010-02-09 Halliburton Energy Services, Inc. Method for azimuthal resistivity measurement and bed boundary detection
GB2407170B (en) * 2002-12-31 2005-06-22 Schlumberger Holdings System and method for locating a fracture in an earth formation
US7202670B2 (en) * 2003-08-08 2007-04-10 Schlumberger Technology Corporation Method for characterizing a subsurface formation with a logging instrument disposed in a borehole penetrating the formation
GB2417783B (en) * 2003-08-08 2006-09-20 Schlumberger Holdings Methods and apparatus for characterising a subsurface formation
US7423426B2 (en) * 2004-02-09 2008-09-09 Baker Hughes Incorporated Selective excitation in earth's magnetic field nuclear magnetic resonance well logging tool
US7663363B2 (en) * 2004-02-09 2010-02-16 Baker Hughes Incorporated Method and apparatus for high signal-to-noise ratio NMR well logging
US7239145B2 (en) * 2004-03-29 2007-07-03 Schlumberger Technology Center Subsurface electromagnetic measurements using cross-magnetic dipoles
US7027926B2 (en) * 2004-04-19 2006-04-11 Pathfinder Energy Services, Inc. Enhanced measurement of azimuthal dependence of subterranean parameters
US7848887B2 (en) * 2004-04-21 2010-12-07 Schlumberger Technology Corporation Making directional measurements using a rotating and non-rotating drilling apparatus
US8736270B2 (en) 2004-07-14 2014-05-27 Schlumberger Technology Corporation Look ahead logging system
US7755361B2 (en) * 2004-07-14 2010-07-13 Schlumberger Technology Corporation Apparatus and system for well placement and reservoir characterization
US7786733B2 (en) * 2004-07-14 2010-08-31 Schlumberger Technology Corporation Apparatus and system for well placement and reservoir characterization
US20060017443A1 (en) * 2004-07-23 2006-01-26 Baker Hughes Incorporated Deep reading propagation resistivity tool for determination of distance to a bed boundary with a transition zone
US7630872B2 (en) * 2004-09-16 2009-12-08 Schlumberger Technology Corporation Methods for visualizing distances between wellbore and formation boundaries
CA2584068C (fr) * 2004-10-22 2011-05-10 Baker Hughes Incorporated Mesures magnetiques pendant la rotation
US7471088B2 (en) * 2004-12-13 2008-12-30 Baker Hughes Incorporated Elimination of the anisotropy effect in LWD azimuthal resistivity tool data
US7536261B2 (en) * 2005-04-22 2009-05-19 Schlumberger Technology Corporation Anti-symmetrized electromagnetic measurements
CN100429529C (zh) * 2006-01-10 2008-10-29 中国石油天然气集团公司 利用全波列、偶极横波测井资料确定气层的方法
US8931335B2 (en) * 2006-04-07 2015-01-13 Baker Hughes Incorporation Processing of multi-component induction data in the presence of borehole abnormalities
BRPI0711054A2 (pt) 2006-05-04 2011-08-23 Shell Int Research métodos para analisar uma formação subterránea atravessada por um furo de poço e para produzir um fluido de hidrocarboneto mineral de um formação geológica, e, meio legìvel por computador
AU2006344741B2 (en) * 2006-06-19 2011-07-07 Halliburton Energy Services, Inc. Antenna cutout in a downhole tubular
RU2394270C1 (ru) 2006-07-11 2010-07-10 Хэллибертон Энерджи Сервисиз, Инк. Модульный приборный узел для геоуправления
BRPI0713267B1 (pt) * 2006-07-12 2018-06-19 Halliburton Energy Services, Inc. “método para construir uma antena inclinada sobre um segmento recuperável de uma ferramenta de registro durante perfuração, ferramenta de registro durante perfuração, e, ferramenta de furo abaixo recuperável”
AU2007277410B2 (en) * 2006-07-25 2012-02-16 Exxonmobil Upstream Research Company Method for determining physical properties of structures
US8593147B2 (en) 2006-08-08 2013-11-26 Halliburton Energy Services, Inc. Resistivity logging with reduced dip artifacts
WO2008024153A2 (fr) * 2006-08-24 2008-02-28 Exxonmobil Upstream Research Company Système de traitement de données électromagnétiques
WO2008033184A2 (fr) * 2006-09-13 2008-03-20 Exxonmobil Upstream Research Company Inversion rapide des données d'arpentage par reconnaissance électromagnétique
US8466683B2 (en) * 2006-12-14 2013-06-18 Schlumberger Technology Corporation Determining properties of earth formations using the electromagnetic coupling tensor
CN101460698B (zh) 2006-12-15 2013-01-02 哈里伯顿能源服务公司 具有旋转天线结构的天线耦合元件测量工具
WO2008115229A1 (fr) * 2007-03-16 2008-09-25 Halliburton Energy Services, Inc. Systèmes et procédés d'inversion robuste pour instruments de diagraphie de résistivité sensibles en azimut
US8049508B2 (en) * 2007-03-16 2011-11-01 Baker Hughes Incorporated Method and apparatus for determining formation boundary near the bit for conductive mud
US7751280B2 (en) * 2007-03-27 2010-07-06 Schlumberger Technology Corporation Determining wellbore position within subsurface earth structures and updating models of such structures using azimuthal formation measurements
US8497685B2 (en) 2007-05-22 2013-07-30 Schlumberger Technology Corporation Angular position sensor for a downhole tool
US7657377B2 (en) 2007-05-31 2010-02-02 Cbg Corporation Azimuthal measurement-while-drilling (MWD) tool
US7558675B2 (en) * 2007-07-25 2009-07-07 Smith International, Inc. Probablistic imaging with azimuthally sensitive MWD/LWD sensors
US7912648B2 (en) * 2007-10-02 2011-03-22 Baker Hughes Incorporated Method and apparatus for imaging bed boundaries using azimuthal propagation resistivity measurements
US8906700B2 (en) 2007-11-06 2014-12-09 Ambergen, Inc. Methods and compositions for phototransfer
US8932879B2 (en) * 2007-11-06 2015-01-13 Ambergen, Inc. Methods and compounds for phototransfer
RU2475780C2 (ru) * 2007-11-30 2013-02-20 Эксонмобил Апстрим Рисерч Компани Способ для регистрации трещиноватости коллектора и диагональных пластов, использующий трехосные/многокомпонентные измерения анизотропии удельного сопротивления
US20090157361A1 (en) * 2007-12-12 2009-06-18 Toghi Farid Method of well placement modeling and geosteering
US8547783B2 (en) 2007-12-12 2013-10-01 Exxonmobil Upstream Research Company Method and apparatus for evaluating submarine formations
JP5662804B2 (ja) * 2007-12-18 2015-02-04 シュルンベルジェ ホールディングス リミテッドSchlnmberger Holdings Limited 表面電磁探査を改善するシステム及び方法
CA2680869C (fr) 2008-01-18 2011-07-12 Halliburton Energy Services, Inc. Forage guide electromagnetiquement (em) par rapport a un trou de forage existant
CN101525999B (zh) * 2008-03-06 2013-04-24 中国石油化工股份有限公司 电磁随钻测量系统的适应性分析方法
US8278931B2 (en) * 2008-07-14 2012-10-02 Baker Hughes Incorporated Anisotropy orientation image from resistivity measurements for geosteering and formation evaluation
US7991555B2 (en) * 2008-07-30 2011-08-02 Schlumberger Technology Corporation Electromagnetic directional measurements for non-parallel bed formations
US8427162B2 (en) * 2008-08-25 2013-04-23 Baker Hughes Incorporated Apparatus and method for detection of position of a component in an earth formation
US8278928B2 (en) * 2008-08-25 2012-10-02 Baker Hughes Incorporated Apparatus and method for detection of position of a component in an earth formation
WO2010039357A2 (fr) * 2008-10-01 2010-04-08 Schlumberger Canada Limited Outil de diagraphie à antennes ayant des angles d’inclinaison égaux
AU2008365630B2 (en) 2008-12-16 2012-05-03 Halliburton Energy Services, Inc. Azimuthal at-bit resistivity and geosteering methods and systems
US8089268B2 (en) * 2009-03-24 2012-01-03 Smith International, Inc. Apparatus and method for removing anisotropy effect from directional resistivity measurements
US8195400B2 (en) * 2009-05-08 2012-06-05 Smith International, Inc. Directional resistivity imaging using harmonic representations
US8159227B2 (en) * 2009-05-11 2012-04-17 Smith International Inc. Methods for making directional resistivity measurements
US7990153B2 (en) * 2009-05-11 2011-08-02 Smith International, Inc. Compensated directional resistivity measurements
US8497673B2 (en) * 2009-09-28 2013-07-30 Schlumberger Technology Corporation Directional resistivity antenna shield
US8466682B2 (en) * 2009-09-29 2013-06-18 Schlumberger Technology Corporation Apparatus and method for downhole electromagnetic measurement while drilling
US8433518B2 (en) * 2009-10-05 2013-04-30 Schlumberger Technology Corporation Multilevel workflow method to extract resistivity anisotropy data from 3D induction measurements
US8271199B2 (en) * 2009-12-31 2012-09-18 Smith International, Inc. Binning method for borehole imaging
GB2486759B (en) 2010-01-22 2014-09-03 Halliburton Energy Serv Inc Method and apparatus for resistivity measurements
MY159410A (en) 2010-04-15 2017-01-13 Halliburton Energy Services Inc Processing and geosteering with a rotating tool
BR112012027443A2 (pt) * 2010-04-27 2016-07-19 Halliburton Energy Services Inc aparelho, sistema, método de processamento implementado, e, artigo.
EA026512B1 (ru) * 2010-05-21 2017-04-28 Лэндмарк Грэфикс Корпорейшн Способ геонавигации горизонтальной скважины и компьютерный носитель информации, содержащий команды для геонавигации горизонтальной скважины
US8600115B2 (en) 2010-06-10 2013-12-03 Schlumberger Technology Corporation Borehole image reconstruction using inversion and tool spatial sensitivity functions
US9372276B2 (en) 2010-06-10 2016-06-21 Schlumberger Technology Corporation Combinations of axial and saddle coils to create the equivalent of tilted coils for directional resistivity measurements
US9933541B2 (en) * 2010-06-22 2018-04-03 Schlumberger Technology Corporation Determining resistivity anisotropy and formation structure for vertical wellbore sections
US9115569B2 (en) 2010-06-22 2015-08-25 Halliburton Energy Services, Inc. Real-time casing detection using tilted and crossed antenna measurement
US8749243B2 (en) * 2010-06-22 2014-06-10 Halliburton Energy Services, Inc. Real time determination of casing location and distance with tilted antenna measurement
US8844648B2 (en) 2010-06-22 2014-09-30 Halliburton Energy Services, Inc. System and method for EM ranging in oil-based mud
US8756018B2 (en) 2010-06-23 2014-06-17 Schlumberger Technology Corporation Method for time lapsed reservoir monitoring using azimuthally sensitive resistivity measurements while drilling
US8558548B2 (en) 2010-07-28 2013-10-15 Schlumberger Technology Corporation Determining anisotropic resistivity
US9360582B2 (en) 2010-07-02 2016-06-07 Halliburton Energy Services, Inc. Correcting for magnetic interference in azimuthal tool measurements
EP2598915B1 (fr) 2010-07-27 2018-11-07 Exxonmobil Upstream Research Company Inversion de données géophysiques pour l'obtention de paramètres géologiques ou de la lithologie
EP2606452A4 (fr) 2010-08-16 2017-08-16 Exxonmobil Upstream Research Company Réduction de la dimensionnalité du problème de l'inversion conjointe
US9273517B2 (en) 2010-08-19 2016-03-01 Schlumberger Technology Corporation Downhole closed-loop geosteering methodology
US20120109527A1 (en) * 2010-09-17 2012-05-03 Baker Hughes Incorporated Apparatus and Methods for Drilling Wellbores by Ranging Existing Boreholes Using Induction Devices
CN101982734A (zh) * 2010-10-29 2011-03-02 南昌航空大学 一种用于地下磁导航的计算方法
US8536871B2 (en) 2010-11-02 2013-09-17 Schlumberger Technology Corporation Method of correcting resistivity measurements for toll bending effects
US9658360B2 (en) 2010-12-03 2017-05-23 Schlumberger Technology Corporation High resolution LWD imaging
NO346394B1 (no) * 2010-12-23 2022-07-04 China Petroleum & Chem Corp En anordning og fremgangsmåte for brønnlogging og data prosesseringsanordning
US9103936B2 (en) 2011-01-20 2015-08-11 Halliburton Energy Services, Inc. Method and apparatus for determining and displaying geosteering information
US8626446B2 (en) 2011-04-01 2014-01-07 Schlumberger Technology Corporation Method of directional resistivity logging
EP3410160A1 (fr) 2011-04-18 2018-12-05 Halliburton Energy Services Inc. Procédé de traitement et de détection de fond de puits en temps réel de limite de couche pour application de géoguidage
US9453929B2 (en) 2011-06-02 2016-09-27 Exxonmobil Upstream Research Company Joint inversion with unknown lithology
EP2721478A4 (fr) 2011-06-17 2015-12-02 Exxonmobil Upstream Res Co Gel de domaine dans une inversion conjointe
WO2013012470A1 (fr) 2011-07-21 2013-01-24 Exxonmobil Upstream Research Company Pondération adaptative de types de données géophysiques dans une inversion conjointe
EP2737170B1 (fr) * 2011-07-29 2017-08-16 Landmark Graphics Corporation Procédé et système pour corréler une diagraphie mesurée à une diagraphie prédite
US10539009B2 (en) * 2011-08-10 2020-01-21 Scientific Drilling International, Inc. Short range data transmission in a borehole
MX358888B (es) 2011-08-18 2018-08-28 Halliburton Energy Serv Inc Herramientas y métodos mejorados de detección de revestimiento.
US9239403B2 (en) * 2011-08-29 2016-01-19 Hallibburton Energy Services, Inc. Apparatus and methods of controlling recordation of resistivity-related readings in determining formation resistivity
US10823871B2 (en) * 2011-11-15 2020-11-03 Halliburton Energy Services, Inc. Enhanced resistivity measurement with at-bit resistivity sensor
US10371852B2 (en) * 2011-12-21 2019-08-06 Schlumberger Technology Corporation Formation properties from conductivity tensor
EP2802363B1 (fr) 2012-01-09 2018-10-24 Zimmer, Inc. Dispositif composite combinant métal poreux et stimuli osseux
US9540922B2 (en) 2012-03-29 2017-01-10 Schlumberger Technology Corporation Electromagnetic method for obtaining dip azimuth angle
US9547102B2 (en) 2012-06-25 2017-01-17 Halliburton Energy Services, Inc. Resistivity logging systems and methods employing ratio signal set for inversion
AU2012383577B2 (en) 2012-06-25 2015-07-09 Halliburton Energy Services, Inc. Tilted antenna logging systems and methods yielding robust measurement signals
WO2014025354A1 (fr) * 2012-08-10 2014-02-13 Halliburton Energy Services, Inc. Procédés et systèmes pour la modélisation directe d'image de fond de trou de propriétés de formation
RU2601232C2 (ru) 2012-08-20 2016-10-27 Лэндмарк Графикс Корпорейшн Способы и системы для включения мест выделенных точек псевдоповерхности в модели скорости сейсмических волн
AU2012394434A1 (en) * 2012-11-13 2015-04-30 Landmark Graphics Corporation System, method and computer program product for a rug plot for geosteering applications
CN103015995B (zh) * 2012-12-14 2015-06-24 中国电子科技集团公司第二十二研究所 通过带倾斜线圈的电阻率设备的地层倾角测量方法及装置
RU2628000C2 (ru) * 2012-12-31 2017-08-14 Халлибертон Энерджи Сервисез, Инк. Глубокая азимутальная система с применением многополюсных датчиков
MX351130B (es) 2012-12-31 2017-10-03 Halliburton Energy Services Inc Imagenología de formaciones con antenas multipolares.
US10591638B2 (en) 2013-03-06 2020-03-17 Exxonmobil Upstream Research Company Inversion of geophysical data on computer system having parallel processors
RU2627947C2 (ru) * 2013-03-15 2017-08-14 Хэллибертон Энерджи Сервисиз, Инк. Выявление нетрадиционных пластов
US9389332B2 (en) * 2013-04-01 2016-07-12 Oliden Technology, Llc Method and tool for directional electromagnetic well logging
US10345475B2 (en) 2013-04-02 2019-07-09 Schlumberger Technology Corporation Extended 1D inversion of electromagnetic measurements for formation evaluation
US9846255B2 (en) 2013-04-22 2017-12-19 Exxonmobil Upstream Research Company Reverse semi-airborne electromagnetic prospecting
WO2014201297A2 (fr) 2013-06-12 2014-12-18 Well Resolutions Technology Appareil et procédés permettant d'effectuer des mesures de résistivité azimutales
US9207053B2 (en) * 2013-06-21 2015-12-08 Rosemount Aerospace Inc. Harmonic shuttered seeker
WO2015050985A1 (fr) 2013-10-01 2015-04-09 Schlumberger Canada Limited Procédés, supports lisibles par ordinateur et systèmes permettant d'appliquer un traitement unidimensionnel (1d) dans une formation non 1d
US10190408B2 (en) * 2013-11-22 2019-01-29 Aps Technology, Inc. System, apparatus, and method for drilling
CN103670387A (zh) * 2013-12-12 2014-03-26 贝兹维仪器(苏州)有限公司 一种地层定向电阻率测量方法及装置
CN103711474B (zh) * 2013-12-19 2016-08-17 天津大学 一种正交偶极子声电组合测井仪器
US9765613B2 (en) 2014-03-03 2017-09-19 Aps Technology, Inc. Drilling system and electromagnetic telemetry tool with an electrical connector assembly and associated methods
MX2016010526A (es) 2014-03-11 2016-10-31 Halliburton Energy Services Inc Metodos y sistemas de registro de induccion con componentes multiples utilizando inversion de modelo fusionado.
US9581721B2 (en) * 2014-03-29 2017-02-28 Schlumberger Technology Corporation Method for making downhole electromagnetic logging while drilling measurements
US10215878B2 (en) 2014-03-29 2019-02-26 Schlumberger Technology Corporation Gain compensated directional propagation measurements
US10302805B2 (en) 2014-03-30 2019-05-28 Schlumberger Technology Corporation System and methods for obtaining compensated electromagnetic measurements
WO2015152955A1 (fr) 2014-04-03 2015-10-08 Halliburton Energy Services, Inc. Systèmes de diagraphie d'induction à plusieurs composants et procédés utilisant une inversion de fréquence choisie
US9790784B2 (en) 2014-05-20 2017-10-17 Aps Technology, Inc. Telemetry system, current sensor, and related methods for a drilling system
CN105464650A (zh) * 2014-09-04 2016-04-06 中国石油化工股份有限公司 一种随钻测井解释方法
US10401529B2 (en) 2014-10-17 2019-09-03 Halliburton Energy Services, Inc. Fast-changing dip formation resistivity estimation
US10267945B2 (en) * 2014-10-20 2019-04-23 Schlumberger Technology Corporation Use of transverse antenna measurements for casing and pipe detection
AU2014411434B2 (en) * 2014-11-13 2017-12-07 Halliburton Energy Services, Inc. Resistivity logging tools with tilted ferrite elements for azimuthal sensitivity
US9784880B2 (en) * 2014-11-20 2017-10-10 Schlumberger Technology Corporation Compensated deep propagation measurements with differential rotation
CA2968501C (fr) * 2015-01-07 2022-05-24 Halliburton Energy Services, Inc. Parametrage de modele terrestre fonctionnel pour l'inversion de resistivite
US9976413B2 (en) 2015-02-20 2018-05-22 Aps Technology, Inc. Pressure locking device for downhole tools
CN104732562A (zh) * 2015-03-24 2015-06-24 中国石油集团川庆钻探工程有限公司 一种电成像测井图像相似度比对方法
US11307323B2 (en) * 2015-09-30 2022-04-19 Schlumberger Technology Corporation Methods and systems to analyze bed boundary detection
WO2017132098A1 (fr) * 2016-01-25 2017-08-03 Schlumberger Technology Corporation Forage à l'aide d'informations caractérisant des hétérogénéités latérales sur la base de mesures de résistivité directionnelle profonde
GB2561479B (en) * 2016-02-16 2021-08-11 Halliburton Energy Services Inc Generating an earth model from spatial correlations of equivalent earth models
CN106324689B (zh) * 2016-06-24 2018-05-11 杭州迅美科技有限公司 一种水平井地层环境下电阻率各向异性识别方法
CN107762494A (zh) * 2016-08-19 2018-03-06 中国石油集团长城钻探工程有限公司 随钻偶极声波测井仪器以及声波测井方法
CN110191999B (zh) * 2017-02-06 2021-03-16 哈里伯顿能源服务公司 用多个初始猜测进行的多层地床边界距离(dtbb)反演
US9995840B1 (en) * 2017-04-17 2018-06-12 Nabors Drilling Technologies Usa, Inc. Azimuthal minor averaging in a wellbore
RU2655279C1 (ru) * 2017-06-19 2018-05-24 Публичное акционерное общество "Татнефть" имени В.Д. Шашина Способ определения геомеханических параметров горных пород
US11391859B2 (en) 2018-06-29 2022-07-19 Halliburton Energy Services, Inc. Determining formation properties in a geological formation using an inversion process on a modified response matrix associated with a downhole tool
WO2020013812A1 (fr) * 2018-07-10 2020-01-16 Halliburton Energy Services, Inc. Correction d'azimut relatif pour inversion de résistivité
CN109668938B (zh) * 2018-12-13 2024-01-09 长江水利委员会长江科学院 堤防渗漏通道三维磁测电阻率探测装置及方法
US11340378B2 (en) 2019-01-14 2022-05-24 Halliburton Energy Services, Inc. Azimuthal borehole rendering of radioelement spectral gamma data
WO2020159484A1 (fr) * 2019-01-29 2020-08-06 Halliburton Energy Services, Inc. Correction itérative de trou de forage
CN110296652B (zh) * 2019-07-17 2021-05-25 伟志股份公司 一种城市建筑基坑常规变形监测系统及方法
CA3169795A1 (fr) * 2020-02-28 2021-09-02 Stephen Douglas BUTT Appareil pour georadar et procede
US20210285297A1 (en) * 2020-03-13 2021-09-16 Baker Hughes Oilfield Operations, Llc Automated geosteering based on a distance to oil-water contact
CN111827986B (zh) * 2020-09-03 2023-10-10 北京贝威通能源科技集团有限公司 一种近钻头无线短传系统和方法
CN112963142B (zh) * 2021-03-26 2024-02-02 北京吉星恒大能源科技有限公司 一种双制式无线传输高低电阻率测量近钻头系统
CN113359199B (zh) * 2021-07-06 2022-07-08 北京航空航天大学 一种基于聚焦磁场的井周电阻率测量方法
CN113885085B (zh) * 2021-08-05 2023-06-09 中煤科工集团西安研究院有限公司 一种井下直流轴向偶极动源超前探测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5883515A (en) * 1993-07-21 1999-03-16 Western Atlas International, Inc. Method of determining formation resistivity utilizing combined measurements of inductive and galvanic logging instruments
US5992228A (en) * 1996-10-23 1999-11-30 Dunham; Lanny L. Method for determining resistivity derived porosity and porosity derived resistivity
WO2002071101A1 (fr) * 2001-03-01 2002-09-12 Baker Hughes Incorporated Procede de focalisation base sur un fond non homogene permettant de mesurer l'induction de rangees dans un puits de forage devie

Family Cites Families (47)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4319191A (en) 1980-01-10 1982-03-09 Texaco Inc. Dielectric well logging with radially oriented coils
SU998995A1 (ru) * 1981-09-11 1983-02-23 Всесоюзный научно-исследовательский институт нефтепромысловой геофизики Устройство дл электромагнитного каротажа
US4689569A (en) * 1984-12-17 1987-08-25 Southwest Research Institute Directional antenna system for use in a borehole incorporating antenna dipole elements
US4766384A (en) 1986-06-20 1988-08-23 Schlumberger Technology Corp. Well logging apparatus for determining dip, azimuth, and invaded zone conductivity
US5115198A (en) 1989-09-14 1992-05-19 Halliburton Logging Services, Inc. Pulsed electromagnetic dipmeter method and apparatus employing coils with finite spacing
US4980643A (en) * 1989-09-28 1990-12-25 Halliburton Logging Services, Inc. Induction logging and apparatus utilizing skew signal measurements in dipping beds
AU654346B2 (en) 1991-05-28 1994-11-03 Schlumberger Technology B.V. Slot antenna having two nonparallel elements
US5235285A (en) 1991-10-31 1993-08-10 Schlumberger Technology Corporation Well logging apparatus having toroidal induction antenna for measuring, while drilling, resistivity of earth formations
RU2069878C1 (ru) 1992-02-10 1996-11-27 Кузьмичев Олег Борисович Способ электромагнитного каротажа скважин
JP2534193B2 (ja) 1993-05-31 1996-09-11 石油資源開発株式会社 指向性インダクション検層法および装置
US5517464A (en) 1994-05-04 1996-05-14 Schlumberger Technology Corporation Integrated modulator and turbine-generator for a measurement while drilling tool
US5757191A (en) 1994-12-09 1998-05-26 Halliburton Energy Services, Inc. Virtual induction sonde for steering transmitted and received signals
AU707977B2 (en) 1996-07-01 1999-07-22 Shell Internationale Research Maatschappij B.V. Electrical logging of a laminated earth formation
RU2107313C1 (ru) 1996-07-12 1998-03-20 Дворецкий Петр Иванович Способ геофизических исследований скважин сложной конфигурации, основанный на применении направленных широкополосных электромагнитных импульсов, возбуждаемых щелевой цилиндрической антенной решеткой
US5781436A (en) 1996-07-26 1998-07-14 Western Atlas International, Inc. Method and apparatus for transverse electromagnetic induction well logging
US6044325A (en) 1998-03-17 2000-03-28 Western Atlas International, Inc. Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument
US6068957A (en) 1998-06-29 2000-05-30 Eastman Kodak Company Lubricating layer in photographic elements
US6476609B1 (en) 1999-01-28 2002-11-05 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for geosteering within a desired payzone
US6163155A (en) 1999-01-28 2000-12-19 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations
US6181138B1 (en) 1999-02-22 2001-01-30 Halliburton Energy Services, Inc. Directional resistivity measurements for azimuthal proximity detection of bed boundaries
US6304086B1 (en) 1999-09-07 2001-10-16 Schlumberger Technology Corporation Method and apparatus for evaluating the resistivity of formations with high dip angles or high-contrast thin layers
US6594584B1 (en) 1999-10-21 2003-07-15 Schlumberger Technology Corporation Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument
US6297639B1 (en) 1999-12-01 2001-10-02 Schlumberger Technology Corporation Method and apparatus for directional well logging with a shield having sloped slots
US6351127B1 (en) 1999-12-01 2002-02-26 Schlumberger Technology Corporation Shielding method and apparatus for selective attenuation of an electromagnetic energy field component
US6566881B2 (en) 1999-12-01 2003-05-20 Schlumberger Technology Corporation Shielding method and apparatus using transverse slots
US6727705B2 (en) 2000-03-27 2004-04-27 Schlumberger Technology Corporation Subsurface monitoring and borehole placement using a modified tubular equipped with tilted or transverse magnetic dipoles
US6836218B2 (en) 2000-05-22 2004-12-28 Schlumberger Technology Corporation Modified tubular equipped with a tilted or transverse magnetic dipole for downhole logging
US6393363B1 (en) * 2000-06-28 2002-05-21 Schlumberger Technology Corporation Method and apparatus for the measurement of the electrical resistivity of geologic formations employing modeling data
US6509738B1 (en) * 2000-07-14 2003-01-21 Schlumberger Technology Corporation Electromagnetic induction well logging instrument having azimuthally sensitive response
US6573722B2 (en) 2000-12-15 2003-06-03 Schlumberger Technology Corporation Method and apparatus for cancellation of borehole effects due to a tilted or transverse magnetic dipole
US6541979B2 (en) 2000-12-19 2003-04-01 Schlumberger Technology Corporation Multi-coil electromagnetic focusing methods and apparatus to reduce borehole eccentricity effects
US6636045B2 (en) * 2001-04-03 2003-10-21 Baker Hughes Incorporated Method of determining formation anisotropy in deviated wells using separation of induction mode
US6584408B2 (en) 2001-06-26 2003-06-24 Schlumberger Technology Corporation Subsurface formation parameters from tri-axial measurements
US6969994B2 (en) 2001-09-26 2005-11-29 Schlumberger Technology Corporation Directional electromagnetic measurements insensitive to dip and anisotropy
US6556015B1 (en) 2001-10-11 2003-04-29 Schlumberger Technology Corporation Method and system for determining formation anisotropic resistivity with reduced borehole effects from tilted or transverse magnetic dipoles
US6819110B2 (en) 2002-03-26 2004-11-16 Schlumberger Technology Corporation Electromagnetic resistivity logging instrument with transverse magnetic dipole component antennas providing axially extended response
US6667620B2 (en) 2002-03-29 2003-12-23 Schlumberger Technology Corporation Current-directing shield apparatus for use with transverse magnetic dipole antennas
US6930652B2 (en) 2002-03-29 2005-08-16 Schlumberger Technology Corporation Simplified antenna structures for logging tools
US6690170B2 (en) 2002-03-29 2004-02-10 Schlumberger Technology Corporation Antenna structures for electromagnetic well logging tools
US6998844B2 (en) 2002-04-19 2006-02-14 Schlumberger Technology Corporation Propagation based electromagnetic measurement of anisotropy using transverse or tilted magnetic dipoles
US6794875B2 (en) 2002-05-20 2004-09-21 Halliburton Energy Services, Inc. Induction well logging apparatus and method
US6903553B2 (en) 2002-09-06 2005-06-07 Baker Hughes Incorporated Method and apparatus for a quadrupole transmitter for directionally sensitive induction tool
US6937022B2 (en) 2002-09-06 2005-08-30 Baker Hughes Incorporated Method and apparatus for a quadrupole transmitter for directionally sensitive induction tool
US6819111B2 (en) * 2002-11-22 2004-11-16 Baker Hughes Incorporated Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electro-magnetic antenna combination and orientation with additional rotation and position measurements
US6937021B2 (en) 2002-12-09 2005-08-30 Schlumberger Technology Corporation Method and apparatus for determining the presence and orientation of a fraction in an earth formation
US6924646B2 (en) 2002-12-31 2005-08-02 Schlumberger Technology Corporation System and method for locating a fracture in an earth formation
US7202670B2 (en) 2003-08-08 2007-04-10 Schlumberger Technology Corporation Method for characterizing a subsurface formation with a logging instrument disposed in a borehole penetrating the formation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5883515A (en) * 1993-07-21 1999-03-16 Western Atlas International, Inc. Method of determining formation resistivity utilizing combined measurements of inductive and galvanic logging instruments
US5992228A (en) * 1996-10-23 1999-11-30 Dunham; Lanny L. Method for determining resistivity derived porosity and porosity derived resistivity
WO2002071101A1 (fr) * 2001-03-01 2002-09-12 Baker Hughes Incorporated Procede de focalisation base sur un fond non homogene permettant de mesurer l'induction de rangees dans un puits de forage devie

Also Published As

Publication number Publication date
FR2855272B1 (fr) 2013-02-08
US7382135B2 (en) 2008-06-03
BRPI0401805A (pt) 2005-01-18
EA200400573A1 (ru) 2004-12-30
FR2855272A1 (fr) 2004-11-26
CN1573013B (zh) 2010-09-01
FR2859539B1 (fr) 2013-01-18
NO20042106L (no) 2004-11-23
GB2402489A (en) 2004-12-08
NO335727B1 (no) 2015-02-02
US20050140373A1 (en) 2005-06-30
EA007587B1 (ru) 2006-12-29
CN1573013A (zh) 2005-02-02
CN101082276A (zh) 2007-12-05
DE102004024969A1 (de) 2006-01-05
MXPA04004692A (es) 2004-11-25
BRPI0401805B1 (pt) 2016-07-26
CN101082276B (zh) 2012-03-28
GB0410076D0 (en) 2004-06-09
GB2402489B (en) 2005-11-30

Similar Documents

Publication Publication Date Title
FR2859539A1 (fr) Procede de resistivite electromagnetique directionnelle
Plessix et al. Multiparameter full-waveform inversion: Marine and land examples
EA006075B1 (ru) Электромагнитный способ определения углов падения независимо от типа бурового раствора и околоскважинного пространства
CN103975118B (zh) 钻头应用的前瞻预测
BR112018072718B1 (pt) Método e aparelho para estimativa de um parâmetro invertido de uma formação subsuperficial
WO2013036509A1 (fr) Evaluation de pendage et d&#39;anisotropie de formation en temps réel à l&#39;aide de mesures d&#39;induction multiaxiale
RU2615219C2 (ru) Обработка анизотропии в пологих скважинах
US11604301B2 (en) Methods and systems for automated sonic imaging
Di et al. Estimating subsurface properties using a semisupervised neural network approach
US10598816B2 (en) Systems and methods for relative dip correction
Bittar et al. Advancement and economic benefit of geosteering and well-placement technology
BR112018017328B1 (pt) Método e sistema para caracterização de uma formação subterrânea
CA2846626A1 (fr) Appareil, procedes et systemes d&#39;estimation analytique
US20220351037A1 (en) Method and system for spectroscopic prediction of subsurface properties using machine learning
FR3035723A1 (fr) Estimation de parametre anisotrope a partir de donnees vsp autonomes utilisant une evolution differentielle
US10705241B2 (en) Determining sea water resistivity
Looms et al. Geostatistical inference using crosshole ground-penetrating radar
FR3062673A1 (fr) Pixellisation par solution d inversion de la distance par rapport au lit
BR0116453B1 (pt) Aparelho e método para converter os sinais de saída de uma ferramenta de registro
Bi et al. Geologic-time-based interpolation of borehole data for building high-resolution models: Methods and applications
US20220397699A1 (en) System and method for enhancing petrophysical characterization of porous media
FR3035145A1 (fr)
FR3035146A1 (fr) Appareils, procedes, et systemes de mesure de proprietes de formation
Hilburn et al. Uncertainty in orthorhombic model building: Analysis, mitigation, and validation
CN105672999B (zh) 钻头应用的前瞻预测

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 13

ST Notification of lapse

Effective date: 20180131