FR2797330A1 - Procede de determination de parametres des formations traversees par un forage - Google Patents
Procede de determination de parametres des formations traversees par un forage Download PDFInfo
- Publication number
- FR2797330A1 FR2797330A1 FR9910178A FR9910178A FR2797330A1 FR 2797330 A1 FR2797330 A1 FR 2797330A1 FR 9910178 A FR9910178 A FR 9910178A FR 9910178 A FR9910178 A FR 9910178A FR 2797330 A1 FR2797330 A1 FR 2797330A1
- Authority
- FR
- France
- Prior art keywords
- parameters
- layer
- process according
- observable
- determining
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 34
- 238000005755 formation reaction Methods 0.000 claims abstract description 34
- 238000005259 measurement Methods 0.000 claims description 18
- 238000005553 drilling Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 9
- 238000009826 distribution Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000009545 invasion Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 238000012549 training Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000706 filtrate Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000010802 sludge Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/20—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with propagation of electric current
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Chemical & Material Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Remote Sensing (AREA)
- Mining & Mineral Resources (AREA)
- Pathology (AREA)
- Geophysics (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
L'invention concerne un procédé de détermination de paramètres des formations traversées par un forage, à partir d'une diagraphie de résistivité (20) enregistrée dans ce forage à l'aide d'un outil d'enregistrement et de mesure, ce procédé comprenant les étapes consistant à déterminer les paramètres de formation par une méthode d'inversion paramétrique, de manière à obtenir un modèle (16) des formations, calculer (17) la réponse de l'outil au modèle, comparer (19) à l'aide d'un critère de comparaison la réponse calculée à la diagraphie enregistrée, et effectuer au moins une nouvelle itération si le critère de comparaison n'est pas satisfait. L'étape (11) de détermination des paramètres à partir des données de diagraphie est effectuée par une méthode de quasi-Newton; et cette méthode de quasi-Newton est mise en oeuvre sur des pseudo-paramètres (12) homogènes déterminés à partir des paramètres de formations.
Description
PROCEDE DE DETERMINATION DE PARAMETRES DES FORMATIONS TRAVERSEES PAR UN FORAGE La présente invention concerne un procédé de détermination de paramètres des formations traversées par un forage, et plus particulièrement un tel procédé de détermination de paramètres à partir d'une diagraphie de résistivité enregistrée dans ce forage à l'aide d'un outil d'enregistrement et de mesure, ledit procédé comprenant les étapes consistant à déterminer lesdits paramètres de formation par une méthode d'inversion paramétrique, de manière à obtenir un modèle des formations, calculer la réponse dudit outil audit modèle, comparer à l'aide d'un critère de comparaison ladite réponse calculée à la diagraphie enregistrée, et effectuer au moins une nouvelle itération si ledit critère de comparaison n'est pas satisfait.
De tels procédés sont connus. D'une manière générale, leur mise en oeuvre est précédée d'une première phase d'initialisation des paramètres.
II est bien connu de réaliser dans des forages des diagraphies de résistivité à l'aide d'outils conçus pour mesurer la résistivité des formations traversées par ce forage en créant dans ces formations des courants galvaniques ou des courants de Foucault aux environs de l'outil. Ces outils donnent à chacune des profondeurs auxquelles sont effectuées les mesures, un ensemble de valeurs de résistivité à cette profondeur pour différentes distances de l'axe du forage.
La relation entre les mesures de résistivité ainsi effectuées et les caractéristiques effectives des formations traversées par le forage sont typiquement affectées par trois types d'effets a) effet de trou dû à la présence du forage, généralement remplis d'une boue de forage beaucoup plus conductrice que les formations ; distribution de conductivité). Le procédé d'interprétation consiste alors à déconvoluer la diagraphie observée à l'aide d'un filtre connu de manière à retrouver la distribution de résistivité. Cette étape correspond à une correction des données de départ des effets d'éponte.
Cette étape est précédée ou suivie d'une étape de correction radiale visant à corriger les résultats de l'effet d'invasion.
Le filtre de convolution peut également être estimé en calculant la réponse de l'outil à une formation possédant une distribution de résistivité proche de la réalité, ce qui correspond à une linéarisation locale du filtre.
Les procédés de ce type sont limités par le fait que les corrections des effets d'éponte et d'invasion y sont considérées comme indépendants alors qu'ils ne peuvent généralement, dans la réalité, être séparés.
Un autre type de procédé connu consiste à diviser la formation en cellules.
On cherche alors à obtenir une image de la formation en termes de résistivité. La formation est divisée en cellules généralement rectangulaires, suivant la direction du forage et la perpendiculaire à cette direction, et une résistivité est déterminée pour chacune des cellules. On utilise alors un algorithme de calcul, par exemple, par éléments finis, par différences finies ou utilisant des réseaux de neurones, pour calculer la réponse de l'outil aux formations ainsi modélisées, de manière à déterminer si la distribution de résistivité supposée explique les résistivités apparentes.
Les frontières horizontales des cellules peuvent être déterminées par un algorithme de segmentation et les frontières verticales sont souvent fixées par l'utilisateur. D'une manière générale, l'utilisation de cellules de grande dimension conduit à une faible définition et donc à une mauvaise approximation de la réalité, et l'utilisation de petites cellules peut conduire à une instabilité qui nécessite de fixer des contraintes de douceur dans la variation de la résistivité, ces contraintes déformant la solution.
En effet, les inconnues du problème inverse qui sont uniquement des inconnues électriques sont ici les résistivités de b) effet d'épontes dû à la présence au-dessus et au-dessous de la zone étudiée, de formations généralement hétérogènes, qui peuvent être bien plus conductrices ou bien plus résistives que les formations qui se trouvent en face de l'appareil (il est à noter que l'effet d'éponte peut être considérable même si la sonde est tout entièrement contenue dans un banc homogène entouré d'épontes de résistivités différentes) ; c) effet d'invasion dû à la présence au voisinage du forage de filtrat de boue de forage qui a remplacé partiellement les fluides initialement présents dans les formations.
En première approximation, on représente généralement la zone envahie par une région s'étendant entre le forage et un cylindre de diamètre dl coaxial au forage et ayant une résistivité Rx0 radialement uniforme et à l'extérieur de laquelle se trouve la zone vierge de résistivité R,. Ce modèle comporte trois inconnues : R,, R_ro et d;, raison pour laquelle on enregistre simultanément au moins trois mesures ayant des profondeurs d'investigation radiales différentes afin de déterminer les trois inconnues.
Le problème consiste donc à déterminer pour chaque profondeur un ensemble de paramètres géométriques et électriques, ici par exemple le diamètre d'invasion d; et les résistivités RXo de zones envahies et R, de zones vierges, à partir d'un ensemble de mesures de résistivité effectuées à des profondeurs séparées par un certain pas de mesure, les résistivités mesurées à chaque profondeur étant chacune caractéristique de la résistivité de la formation à une certaine distance de l'axe du forage.
On connaît un certain nombre de procédés susceptibles de résoudre ce problème. Ces procédés sont essentiellement de trois types.
Un premier type de procédé connu utilise des filtres de déconvolution.
Cette approche suppose que le signal observé est le résultat d'une convolution entre la distribution effective des résistivités et un filtre qui représente la réponse de l'outil à une distribution de résistivité (ou de conductivité si la convolution est réalisée avec une chacune des cellules. Ces résistivités présentent un grand nombre de degrés de liberté, pouvant conduire à l'instabilité précitée. Comme indiqué ci-dessus, une manière de résoudre ces problèmes d'instabilité consiste à imposer des contraintes mais ces dernières ont l'inconvénient de fausser la solution.
Le problème est plus facile à résoudre lorsque le nombre de résistivités décrivant une couche est relativement faible. Mais alors, le modèle de la formation est moins précis. Une couche divisée en un faible nombre de zones radiales ne constitue pas une approximation précise du modèle en Rx0, di et R@ Cette insuffisance introduit des erreurs systématiques dans le processus d'inversion.
Un troisième type de procédé connu consiste à réaliser une inversion locale paramétrique.
L'idée générale sur laquelle est fondée cette approche est d'inverser à l'aide d'une méthode de Newton tous les paramètres inconnus décrivant la formation, tels que les positions des limites des couches ainsi que les valeurs d;, Rxa et R,, c'est-à-dire ici des paramètres géométriques et de résistivité. Afin d'obtenir suffisamment de grandeurs observables, on utilise plusieurs mesures réalisées par l'outil à proximité de la zone où est réalisée l'inversion. Un critère, notamment un critère quadratique, d'évaluation de l'erreur entre les grandeurs observables et les grandeurs reconstituées est minimisé.
L'avantage de ce procédé est sa souplesse. Des modèles de formation très variés peuvent être utilisés et l'on dispose d'une grande liberté dans le choix des mesures d'entrée.
En revanche, du fait de la généralité de l'approche, le problème inverse est souvent largement non linéaire. En conséquence, seule une inversion locale peut être réalisée et des gradients précis doivent être estimés, ce qui est extrêmement coûteux en temps de calcul.
Par ailleurs, lorsque ce procédé est utilisé sur des données de diagraphie réelle, des instabilités ont été observées qui résultent du caractère local de l'inversion.
La présente invention vise à pallier ces inconvénients.
A cet effet, l'invention a pour objet un procédé de détermination de paramètres des formations traversées par un forage, à partir d'une diagraphie de résistivité enregistrée dans ce forage à l'aide d'un outil d'enregistrement et de mesure, ledit procédé comprenant les étapes consistant à déterminer lesdits paramètres de formation par une méthode d'inversion paramétrique, de manière à obtenir un modèle des formations, calculer la réponse dudit outil audit modèle, comparer à l'aide d'un critère de comparaison ladite réponse calculée à la diagraphie enregistrée, et effectuer au moins une nouvelle itération si ledit critère de comparaison n'est pas satisfait, procédé caractérisé par le fait que l'étape de détermination desdits paramètres à partir des données de diagraphie est effectuée par une méthode de quasi- Newton ; et ladite méthode de quasi-Newton est mise en #uvre sur des pseudo-paramètres homogènes déterminés à partir desdits paramètres de formations.
Le procédé selon l'invention suppose que les limites entre couches sont connues.
Par ailleurs, on utilisera tout procédé connu dans l'art antérieur pour reconstituer la réponse de l'outil à une description donnée de la formation (modèle direct).
Pour chaque couche, un modèle de couche peut être choisi comme étant le plus réaliste. Un modèle de couche est une partition de la couche en différentes zones radiales, possédant chacune une résistivité constante. Pour chaque couche, une partie des paramètres peut être choisie comme étant les inconnues du problème inverse. Aussi bien la résistivité que la position des interfaces radiales sont des inconnues possibles.
A chaque couche, est attaché un certain nombre d"'observables". Un "observable" est une mesure de l'outil que l'on considère comme fonction des inconnues de la couche. Par exemple, si l'outil donne deux mesures à chaque profondeur, les deux mesures obtenues à proximité du milieu de la couche peuvent être choisies comme fonction de deux paramètres du modèle de couche.
Les observables peuvent résulter de la même mesure mais prises à différentes profondeurs ou différentes positions de l'outil. La règle du choix des observables dépend de l'outil et du modèle à inverser. La généralité du procédé selon l'invention permet l'inversion d'un large choix d'observables.
L'invention part par conséquent des procédés d'inversion paramétrique locale. Toutefois, on réalise maintenant une optimisation non locale, c'est à dire en même temps sur l'ensemble des couches.
On s'affranchit de l'effet de trou en supposant connue la résistivité de la boue ainsi que le diamètre du forage. Par ailleurs, les positions des frontières horizontales entre couches sont connues, par exemple par le passage à zéro de la dérivée seconde de la conductivité.
On suppose également connues les valeurs des résistivités à certaines profondeurs radiales, par exemple les résistivités RLLS et RLLd ou d'autres distributions radiales de résistivité.
Les inconnues sont à chaque profondeur Rxo, R@ et di.
Si l'on souhaite effectivement déterminer Rx0, il est alors nécessaire de connaître la microrésistivité dans l'hypothèse où les valeurs mesurées sont RLt,, et RLLj.
Quoiqu'il en soit, il est nécessaire, en théorie, que le nombre d'inconnues soit inférieur au nombre d'observables. En pratique, le procédé selon l'invention impose que le nombre d'inconnues soit égal au nombre d'observables. Dans le cas où le nombre d'observables est supérieur au nombre d'inconnues, on peut soit ne pas tenir compte d'un certain nombre d'observables, soit prendre pour inconnue des combinaisons d'observables.
L'invention consiste donc pour l'essentiel tout d'abord à appliquer une méthode de quasi-Newton au problème de l'inversion paramétrique, et d'autre part, à améliorer la convergence, dans une méthode de quasi-Newton généralement affectée par le mélange de paramètres géométriques tels que di, et de paramètres électriques tels que Rx0, et R@, en prenant pour inconnues des combinaisons de ces paramètres géométriques et électriques consistant par exemple toutes en des pseudo-grandeurs électriques et notamment des pseudo-résistivités.
On réalise par conséquent en général deux transformations lors de la mise en #uvre de l'invention.
Une première transformation est généralement réalisée dans l'espace des observables pour réduire le nombre d'observables au nombre d'inconnues. Ceci est nécessaire pour pouvoir appliquer l'algorithme d'inversion qui suppose égaux les nombres d'observables et d'inconnus.
La deuxième transformation est réalisée dans l'espace des inconnues pour obtenir un ensemble d'inconnues homogènes, notamment de nature électrique, à partir d'inconnues de natures inhomogènes, par exemple d'inconnues de nature géométrique et d'inconnues de nature électrique.
On décrira maintenant, à titre d'exemple non limitatif, un mode de mise en #uvre particulier de l'invention, en référence aux dessins schématiques annexés dans lesquels - la figure 1 est un organigramme d'ensemble d'un procédé selon l'invention ; et - la figure 2 est un organigramme d'une itération de l'algorithme de calcul.
Comme montré à la figure 1, la mise en #uvre du procédé selon l'invention commence par une étape 1 de mesure d'un certain nombre de grandeurs physiques de la formation à l'intérieur d'un forage. Le nombre de grandeurs ainsi mesurées doit être, à chaque niveau de mesure, au moins égal au nombre de paramètres dont on souhaite obtenir les valeurs pour chacun de ces niveaux.
A titre d'exemple, dans le cas courant où les paramètres que l'on souhaite déterminer sont les valeurs de résistivité RXo et R, et la distance d'invasion d;, soit trois paramètres, on effectuera les mesures de trois grandeurs au moins à chaque pas.
On pourra par exemple mesurer les grandeurs RLLs et RLLd <I>à</I> l'aide d'un outil connu sous le nom de Dual Laterolog (marque de la société Schlumberger) et mesurer directement Rx0 à l'aide de l'appareil connu sous le nom de MSFL (marque de la société Schlumberger ; "Micro Sphericaly Focused Log"). Dans un autre mode de mise en #uvre de cette étape du procédé, on pourra mesurer les résistivités RLAI, ... RLA5 obtenues à l'aide d'un appareil connu sous le nom de HRLA (marque de la société Schlumberger ; "High Resolution Laterolog Array"), qui correspondent à des résistivités de la formation mesurées à des distances croissantes de l'axe du forage.
L'étape suivante 2 consiste à déterminer l'emplacement des bancs ou couches dont on souhaite déterminer les paramètres. Cette détermination de l'emplacement des bancs peut être effectuée automatiquement à l'aide d'un algorithme de segmentation de type connu, fondé par exemple sur la recherche des points d'inflexion des valeurs de la conductivité apparente.
Au cours de cette étape, on affecte une première valeur d'initialisation aux paramètres.
L'étape suivante 3 consiste à choisir le nombre convenable de grandeurs "observables". Si l'on se place toujours dans le cas où l'on souhaite déterminer Rr0, R@ et di, on prendra comme observable les trois mesures obtenues à l'aide du Dual Laterolog et du MSFL, dans le premier cas.
Dans le deuxième cas où les mesures sont effectuées à l'aide du HRLA, on pourra éliminer une ou deux des valeurs RL,", ... RLA5 etlou obtenir trois observables par combinaison, notamment linéaire, de ces valeurs.
II s'agit au cours de l'étape suivante 4 d'affecter une valeur à chaque observable pour chaque couche déterminée comme indiqué ci-dessus.
On notera en effet que l'on dispose de plus de pas de mesures que de couches puisqu'il faut en particulier au moins trois pas de mesures pour déterminer un point d'inflexion.
II est possible, à cet effet, d'interpoler les mesures effectuées à l'intérieur d'une couche pour obtenir la valeur de chacun des observables au milieu de cette couche. II est également possible d'affecter à chaque couche les valeurs d'observables mesurées au point de mesure le plus proche de milieu de la couche.
L'étape suivante 5 consiste à déterminer les paramètres recherchés.
Mathématiquement, le problème se pose alors de la manière suivante.
Ce problème est le problème inverse du problème direct selon lequel on cherche à déterminer la réponse des outils de mesure, c'est-à-dire les observables du problème de l'invention, à partir des paramètres effectifs de cette formation, c'est-à-dire les inconnues du problème de l'invention.
Soit f(U)=O pour le problème direct et U =#'(O) pour le problème de l'invention, à savoir le problème inverse du problème direct.
On sait résoudre le problème direct, c'est-à-dire évaluer le système d'équations f non linéaires en fonction du modèle de couches et des observables préalablement choisis. A titre d'exemple, dans le cas où les caractéristiques physiques considérées sont des résistivités, ce modèle de couche consiste en la définition de la variation de la résistivité du terrain en fonction de la distance à l'axe du forage, ou encore le modèle Rx0, R, et di.
La solution du problème inverse nécessite l'utilisation d'une méthode itérative c'est-à-dire la détermination d'une suite Ul, U2, ... Ux, qui tend vers la solution Usolution, telle que f(Usolution)-Oobservé Dans les procédés conventionnels, un terme U+ de la suite ci- dessus est obtenu à partir du terme courant U, par la formule suivante U+ = Uc-(#f)-1*[f(Uc)- oobservé où Vf est le jacobien de la fonction f(U).
On voit donc le nombre considérable d'évaluations de la fonction f pour l'estimation de la matrice NxN lorsque l'on sait que le système d'équation peut typiquement comporter 1000 observables et 1000 inconnues.
Selon l'invention, on utilise donc non pas une méthode de Newton selon laquelle, à chaque itération, on estime un nouveau vecteur U et on calcule un nouveau jacobien, mais une méthode dite de quasi-Newton dans laquelle on réalise une estimation à la fois du vecteur U et du jacobien à chaque étape. Si l'estimation courante du jacobien est Bc, l'équation (1) est alors remplacée par U+ = Ur - Bc-1*[f (Uc)- Oobsercé ] (2) L'obtention à chaque itération de la nouvelle estimation du jacobien peut être par exemple obtenue par la méthode de Broyden décrite aux pages 113 à 131 de l'ouvrage de C. T. Kelley "Iterative Method for Linear and Nonlinear Equations" publié par Society for Industrial and Applied Mathematics.
Plus précisément, la matrice B+1 à l'étape suivante est obtenue à partir de la matrice courante Bc#1par la formule suivante:
où lest la matrice identité et où s--U±Uc# II n'est donc plus nécessaire, avec ce procédé, de connaître le jacobien.
où lest la matrice identité et où s--U±Uc# II n'est donc plus nécessaire, avec ce procédé, de connaître le jacobien.
Toutefois, pour obtenir une convergence rapide de la suite des B-', on préconditionne les données initiales de manière que celles-ci soient homogènes. Ceci est obtenu dans le cas présent en utilisant non pas les variables RO, R, et d; qui sont inhomogènes du fait de la présence de d,, mais les résistivités Rxo et R, et la pseudo-résistivité Ra définie par
où dimax et di,,,in sont des paramètres fixés de manière à correspondre à des valeurs maximum et minimum admissibles pour d;.
où dimax et di,,,in sont des paramètres fixés de manière à correspondre à des valeurs maximum et minimum admissibles pour d;.
On voit à la figure 2 l'organigramme complet d'une itération.
On part par exemple d'une matrice identité et en 11 on applique une itération de la méthode de quasi-Newton qui vient d'être décrite pour obtenir les pseudo-résistivités R1, R2 et R3 pour chaque couche, comme montré en 12.
On transforme alors en 13 les pseudo-résistivités R1, R2 et R3 en grandeurs physiques 14 Rxo, Rl et di.
On peut alors procéder en 15 à une mise à jour pour chaque couche du modèle de la formation pour obtenir le modèle de formation 16 mis à jour en Rxo, Rl et di.
On calcule alors en 17 de façon connue, par exemple par une méthode d'éléments finis, la réponse de l'outil à la formation 16 de manière à obtenir une diagraphie simulée 18.
Cette dernière est comparée en 19 à la diagraphie réelle 20. Si l'écart entre la diagraphie simulée et la diagraphie réelle satisfont un critère prédéterminé, l'itération se termine en 21. Dans le cas contraire, une nouvelle itération est entreprise.
Claims (1)
- REVENDICATIONS 1 - Procédé de détermination de paramètres des formations traversées par un forage, à partir d'une diagraphie de résistivité (20) enregistrée dans ce forage à l'aide d'un outil d'enregistrement et de mesure, ledit procédé comprenant les étapes consistant à déterminer lesdits paramètres de formation par une méthode d'inversion paramétrique, de manière à obtenir un modèle des formations, calculer (17) la réponse dudit outil audit modèle, comparer (19) à l'aide d'un critère de comparaison ladite réponse calculée à la diagraphie enregistrée, et effectuer au moins une nouvelle itération si ledit critère de comparaison n'est pas satisfait, procédé caractérisé par le fait que l'étape (11) de détermination desdits paramètres à partir des données de diagraphie est effectuée par une méthode de quasi- Newton ; et ladite méthode de quasi-Newton est mise en #uvre sur des pseudo-paramètres (12) homogènes déterminés à partir desdits paramètres de formations. 2 - Procédé selon la revendication 1, comprenant une étape (2) de détermination des limites entre couches géologiques, préalable à la mise en oeuvre de ladite méthode de quasi-Newton. 3 - Procédé selon la revendication 2, dans lequel la détermination des limites entre couches est effectuée sur la base des points d'inflexion des données de diagraphie. 4 - Procédé selon l'une quelconque des revendications 1 à 3, comprenant une étape de choix d'un modèle de couche pour chaque couche géologique, préalable à la mise en oeuvre de ladite méthode de quasi-Newton. 5 - Procédé selon la revendication 4, dans lequel le modèle de couche est constitué de paramètres de distance à l'axe du forage de manière à délimiter des zones radiales autour de cet axe, et d'un paramètre de résistivité à l'intérieur de chaque zone radiale ainsi délimitée. 6 - Procédé selon l'une quelconque des revendications 1 à 5, comprenant une étape (3) de choix des grandeurs observables. 7 - Procédé selon la revendication 6, dans lequel le choix des grandeurs observables comprend la définition d'une combinaison des données de diagraphie. 8 - Procédé selon l'ensemble des revendications 4 et 6, comprenant une étape d'affectation d'une valeur à chaque observable pour chaque couche géologique. 9 - Procédé selon la revendication 8, dans lequel l'affectation d'une valeur à chaque observable pour chaque couche géologique comprend l'interpolation à l'intérieur de chaque couche des valeurs de l'observable déterminées à l'intérieur de cette couche 10 - Procédé selon la revendication 8, dans lequel l'affectation d'une valeur à chaque observable pour chaque couche géologique comprend l'affectation à cet observable de la valeur qu'il possède au point de mesure le plus proche du milieu de la couche. 11 - Procédé selon l'une quelconque des revendications 1 à 10, dans lequel l'étape (11) de détermination desdits paramètres à partir des données de diagraphie par une méthode de quasi-Newton est effectuée en estimant le jacobien du problème par la méthode de Broyden. 12 - Procédé selon l'une quelconque des revendications 1 à 11, dans lequel la diagraphie utilisée est celle de RLLsetRLLd. 13 - Procédé selon l'une quelconque des revendications 1 à 11, dans lequel la diagraphie utilisée est celle de RLr", <I>...</I> RL,q5. 14 - Application du procédé selon l'une quelconque des revendications 1 à 13 à la détermination des paramètres R,, Rxo et d,.
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9910178A FR2797330B1 (fr) | 1999-08-05 | 1999-08-05 | Procede de determination de parametres des formations traversees par un forage |
AU66984/00A AU6698400A (en) | 1999-08-05 | 2000-07-25 | A method of determining parameters of formations through which a borehole passes |
CA2380382A CA2380382C (fr) | 1999-08-05 | 2000-07-25 | Procede permettant de determiner des parametres de formations traversees d'un trou de forage |
PCT/EP2000/007284 WO2001011391A1 (fr) | 1999-08-05 | 2000-07-25 | Procede permettant de determiner des parametres de formations traversees d'un trou de forage |
GB0200671A GB2369892B (en) | 1999-08-05 | 2000-07-25 | A method of determining parameters of formations through which a borehole passes |
NO20020546A NO20020546L (no) | 1999-08-05 | 2002-02-04 | Fremgangsmåte for å bestemme parametere til formasjoner som et borehull passerer gjennom |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9910178A FR2797330B1 (fr) | 1999-08-05 | 1999-08-05 | Procede de determination de parametres des formations traversees par un forage |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2797330A1 true FR2797330A1 (fr) | 2001-02-09 |
FR2797330B1 FR2797330B1 (fr) | 2002-02-22 |
Family
ID=9548918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR9910178A Expired - Fee Related FR2797330B1 (fr) | 1999-08-05 | 1999-08-05 | Procede de determination de parametres des formations traversees par un forage |
Country Status (6)
Country | Link |
---|---|
AU (1) | AU6698400A (fr) |
CA (1) | CA2380382C (fr) |
FR (1) | FR2797330B1 (fr) |
GB (1) | GB2369892B (fr) |
NO (1) | NO20020546L (fr) |
WO (1) | WO2001011391A1 (fr) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2831917B1 (fr) * | 2001-11-08 | 2004-01-02 | Schlumberger Services Petrol | Procede de determination de la variation de la permeabilite relative a au moins un fluide d'un reservoir contenant des fluides en fonction de la saturation en l'un d'entre eux |
US6611762B1 (en) | 2002-04-04 | 2003-08-26 | Halliburton Energy Services, Inc. | Method for determining parameters of earth formations surrounding a well bore |
US20040019427A1 (en) * | 2002-07-29 | 2004-01-29 | Halliburton Energy Services, Inc. | Method for determining parameters of earth formations surrounding a well bore using neural network inversion |
CN102678102B (zh) * | 2012-03-31 | 2015-02-25 | 中国石油大学(华东) | 基于阵列电成像测井的储层油水识别方法及识别系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4916616A (en) * | 1986-12-08 | 1990-04-10 | Bp Exploration, Inc. | Self-consistent log interpretation method |
US5867806A (en) * | 1996-03-13 | 1999-02-02 | Halliburton Energy Services, Inc. | System and method for performing inversion on LWD resistivity logs with enhanced resolution |
EP0930519A2 (fr) * | 1998-01-15 | 1999-07-21 | Schlumberger Limited | Dispositif et méthode pour évaluer la résistance des formations envahies à angle élevée de pendage apparent |
-
1999
- 1999-08-05 FR FR9910178A patent/FR2797330B1/fr not_active Expired - Fee Related
-
2000
- 2000-07-25 CA CA2380382A patent/CA2380382C/fr not_active Expired - Fee Related
- 2000-07-25 GB GB0200671A patent/GB2369892B/en not_active Expired - Fee Related
- 2000-07-25 AU AU66984/00A patent/AU6698400A/en not_active Abandoned
- 2000-07-25 WO PCT/EP2000/007284 patent/WO2001011391A1/fr active Application Filing
-
2002
- 2002-02-04 NO NO20020546A patent/NO20020546L/no not_active Application Discontinuation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4916616A (en) * | 1986-12-08 | 1990-04-10 | Bp Exploration, Inc. | Self-consistent log interpretation method |
US5867806A (en) * | 1996-03-13 | 1999-02-02 | Halliburton Energy Services, Inc. | System and method for performing inversion on LWD resistivity logs with enhanced resolution |
EP0930519A2 (fr) * | 1998-01-15 | 1999-07-21 | Schlumberger Limited | Dispositif et méthode pour évaluer la résistance des formations envahies à angle élevée de pendage apparent |
Also Published As
Publication number | Publication date |
---|---|
GB0200671D0 (en) | 2002-02-27 |
GB2369892B (en) | 2004-02-25 |
NO20020546D0 (no) | 2002-02-04 |
CA2380382C (fr) | 2010-06-01 |
CA2380382A1 (fr) | 2001-02-15 |
NO20020546L (no) | 2002-02-04 |
GB2369892A (en) | 2002-06-12 |
FR2797330B1 (fr) | 2002-02-22 |
WO2001011391A1 (fr) | 2001-02-15 |
AU6698400A (en) | 2001-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0354112B1 (fr) | Méthode pour obtenir un modèle représentatif d'un milieu hétérogène et notamment du sous-sol | |
EP1669778B1 (fr) | Méthode pour déterminer de l'information spéculaire dans l'imagerie sismique avant sommation | |
FR2910050A1 (fr) | Determination des proprietes d'une formation terrestre a l'aide du tenseur de couplage electromagnetique | |
FR2865810A1 (fr) | Traitement de donnees sismiques | |
EP0131517A2 (fr) | Mécanisme de bras supportant un patin de sonde de diagraphie | |
FR3057903A1 (fr) | Outils a courant de foucault a champ lointain | |
FR2719386A1 (fr) | Procédé et appareil de mesure de la résistivité des formations dans les trous tubes. | |
EP3026430B1 (fr) | Méthode d'imagerie d'un milieu par mesures électriques avec correction d'impédance de contact | |
FR2792731A1 (fr) | Diagraphie de puits par induction electromagnetique et des moyens de mise en oeuvre | |
FR2793032A1 (fr) | Procede et appareil pour determiner la resistivite d'une formation traversee par un puits tube | |
FR2795521A1 (fr) | Procede et dispositif pour determiner la resistivite d'une formation traversee par un puits tube | |
FR2858063A1 (fr) | Procede d'estimation du taux de couverture d'illumination dans le domaine migre | |
CA2464799C (fr) | Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene | |
FR2616920A1 (fr) | Inversion d'un profil sismique vertical en minimisant une fonction du type entropie | |
FR2740169A1 (fr) | Procedes et dispositifs de mesure de caracteristiques d'une formation traversee par un trou de forage | |
FR2797330A1 (fr) | Procede de determination de parametres des formations traversees par un forage | |
FR2666661A1 (fr) | Procede et dispositif de diagraphie pour mesurer la resistivite des formations geologiques. | |
FR2793031A1 (fr) | Procede et appareil pour determiner la resistivite d'une formation traversee par un puits tube | |
EP0944847B1 (fr) | Procede de traitement sismique et notamment procede de prospection sismique 3d mettant en oeuvre une migration des donnees sismiques | |
FR2799552A1 (fr) | Traitement de mesures de formes d'ondes acoustiques provenant d'outils de diagraphie de sondage en matrice | |
EP0410826B1 (fr) | Procédé itératif d'estimation de mouvement, entre une image de référence et une image courante, et dispositif pour la mise en oeuvre de ce procédé | |
FR2853086A1 (fr) | Traitement de donnees sismiques | |
EP0722096B1 (fr) | Dispositif de cartographie d'un milieu par mesure à induction | |
EP1904873B1 (fr) | Procede d'acquisition et de traitement de donnees magnetometriques par des mises a jour locales en temps reel | |
FR2784756A1 (fr) | Methode de detection de ruptures dans des signaux diagraphiques concernant une zone d'un milieu |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |
Effective date: 20160429 |