[go: up one dir, main page]

EP0816629B1 - Méthode et système d'estimation en temps reel d'au moins un paramètre lié au déplacement d'un outil de forage - Google Patents

Méthode et système d'estimation en temps reel d'au moins un paramètre lié au déplacement d'un outil de forage Download PDF

Info

Publication number
EP0816629B1
EP0816629B1 EP97401297A EP97401297A EP0816629B1 EP 0816629 B1 EP0816629 B1 EP 0816629B1 EP 97401297 A EP97401297 A EP 97401297A EP 97401297 A EP97401297 A EP 97401297A EP 0816629 B1 EP0816629 B1 EP 0816629B1
Authority
EP
European Patent Office
Prior art keywords
model
string
displacement
reduced
account
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.)
Expired - Lifetime
Application number
EP97401297A
Other languages
German (de)
English (en)
Other versions
EP0816629A1 (fr
Inventor
Didier Pavone
Jean-Pierre Desplans
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.)
IFP Energies Nouvelles IFPEN
Original Assignee
IFP Energies Nouvelles IFPEN
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 IFP Energies Nouvelles IFPEN filed Critical IFP Energies Nouvelles IFPEN
Publication of EP0816629A1 publication Critical patent/EP0816629A1/fr
Application granted granted Critical
Publication of EP0816629B1 publication Critical patent/EP0816629B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B44/00Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B45/00Measuring the drilling time or rate of penetration

Definitions

  • the present invention relates to the field of measurements during drilling, in particular measures concerning the displacement of a drilling tool attached to the end of a drill string.
  • the method according to the invention provides a solution for estimate in particular the instantaneous speed of rotation of the tool at the bottom of the well, said estimates being obtained by means of a calculation program taking into account measurements taken at the top of the drill string, i.e. substantially on the surface of the ground, generally by means of sensors or an instrumented fitting located in the neighborhood of the means for rotating the lining.
  • the information contained in the surface measurements does not alone are not enough to solve the problem posed, that is to say knowing the instantaneous movements of the tool by knowing the instantaneous movements of the surface trim.
  • the surface measurement information must be supplemented by independent information of another kind which takes into account the structure of the drill string and its behavior between the bottom and the surface: this is the role of the knowledge which establishes the theoretical relationships between the bottom and the surface.
  • the methodology of the present invention uses the conjunction of such a model, defined a priori, and surface measurements acquired in real time.
  • the model can take into account essentially the rotational displacements and said reduced model can calculate in real time the instantaneous speed of rotation of the tool of drilling, said parameter measured at the surface can be taken from at least: the speed of packing rotation or torque.
  • the reduced model can be refined by self-adaptive filtering which minimizes the difference between an actual measurement of a parameter related to the displacement of the lining in surface and the corresponding output obtained by said reduced model.
  • Filtering can take into account the torque measured at the surface or the speed of rotation measured at the surface.
  • the invention also relates to a system for estimating the effective displacement. of a drilling tool attached to the end of a drill string and rotated in a well by means of drive located on the surface, in which an installation of calculation includes means for non-linear physical modeling of the drilling process based on general mechanical equations.
  • the parameters of said means modeling are identified by taking into account the parameters of said well and said lining, and the calculation installation includes means for linearizing said model around an operating point, means for reducing said linearized model so to keep only some of the eigen modes of the state matrix of said model, means for calculating, in real time, the displacement of the drilling tool using the means once linearized and reduced and means of measuring at least one parameter linked to the displacement of the lining on the surface.
  • the modeling means may take into account only the torsion, and the parameters can be the speed of rotation and / or the torque.
  • FIG. 1 illustrates a drilling rig on which we will operate the invention.
  • the surface installation comprises a lifting device 1 comprising a lifting tower 2, a winch 3 which allow the displacement of a drilling hook 4.
  • Under the drill hook are suspended drive means 5 for rotating the assembly of the drill string 6 placed in the well 7.
  • These drive means can be of the drive rod or kelly type coupled to a rotation table 8 and the mechanical motorizations, or of the type motorized drive head or "power swivel" suspended directly from the hook and guided longitudinally in the tower.
  • the drill string 6 is conventionally constituted by rods of drilling 10, of part 11 commonly called BHA for "Bottom Hole Assembly” mainly comprising drill rods, a drilling tool 12 in contact with the land being drilled.
  • BHA Bottom Hole Assembly
  • the well 7 is filled with a fluid, called a drilling fluid, which circulates from the surface at the bottom through the inner channel of the drill string and rises to the surface by the annular space between the walls of the well and the drill string.
  • an instrumented fitting can be inserted 13 between the drive means and the top of the lining.
  • This fitting allows measure the torque, the speed of rotation, possibly the longitudinal displacement of the top of the filling.
  • These so-called surface measurements are transmitted by cable or radio towards an electronic installation of recording, processing, display, not shown here.
  • other sensors such as a tachometer on the rotation table to measure the speed of rotation and a measurement of the torque on the motorized device, if the accuracy of the measurements thus obtained is sufficient.
  • Part 11 of the BHA may more specifically include, drill collars, stabilizers, a shock absorber, and a second instrumented connector 14 which will not used only to experimentally control the present invention by allowing the comparison between the displacement of the drilling tool 12 actually measured by the instrumented fitting 14 and the estimated displacement thanks to the implementation of this invention. It is therefore clear that the application of the present invention does not use instrumented fitting placed at the bottom of the well.
  • the driller who conducts a drilling operation with the devices described in the Figure 1 has three main possible actions, which are therefore the command variables possible to control the drilling process. These are, the weight on the tool which is adjusted by the winch which controls the position of the hook, the speed of rotation of the table rotation or equivalent, and the flow of injected drilling fluid.
  • the described model will treat the drill string as a one-dimensional element vertical. Only rotational displacements will be considered, vertical displacements and lateral being neglected.
  • Figure 2 shows the block diagram of the model; here of twist. It's a model classic with finite differences which comprises several tens of meshes represented by blocks 20. Each mesh represents a part of the drill string, composed inter alia of drill rods and drill rods. It is a figured mass-spring-damping triplet by the diagrams referenced 21, 22, 23. Each block is provided with two inputs and outputs represented by the pairs of arrows 24 and 25 which represent the input and output and input and output rotational speeds. This representation shows the way to digitally connect several rods (or meshes) as we connect physically the packing rods.
  • Block 26 represents the boundary conditions at the head of the resulting drill string of the drilling rig.
  • Block 27 represents the friction law of the tool / ground contact determining the boundary conditions at the bottom of the well. This is a type power dissipator resistant couple.
  • the different curves 28 represent a nonlinear law of tool / rock interaction, depending on the instantaneous speed of the tool (input 29) of the resisting torque (outlet 30) and weight on the tool (inlet 31).
  • This model is validated using data recorded on site using the instrumented bottom and surface fittings.
  • the torsion model thus obtained is generally of high order, that is to say of around 50 to 100 to reproduce reality with sufficient finesse.
  • the torsion model retains the same high order. analysis eigen modes of the torsion model allows to quantify the contribution of each fashion on exits of interest. We only keep the relevant modes; that is to say those who have a notable influence on the dynamic behavior represented by said outputs.
  • the reduction method used is the singular disturbance method. It consists in keeping from the state matrix and from the command matrix, the lines and the columns corresponding to the modes to keep. To keep static gains, the modes Rapids are replaced by their static value, which has the consequence of introducing a direct matrix.
  • the method assumes that the fast modes take their equilibrium in a time negligible, that is to say that they are established instantaneously (quasi-static hypothesis).
  • FIG. 3 shows the block diagram of a loop-type estimation system opened.
  • Block 40 shows diagrammatically the means for measuring surface parameters, here the speed of rotation of the lining Vms measured at the table or at the motorized injection head, and the torque Cms also measured at the surface.
  • Block 41 represents the scale model which simulates the physical model of nonlinear torsion by calculating the transfer function between the input (Vms) and the outputs Ces, Cef and Vef representing respectively the torque of estimated area, estimated bottom torque and estimated bottom speed.
  • the transfer function is always an approximation of reality and any mismatch between the model and the actual drilling process can create a divergence between the estimated values and the real values by integration of the differences. Also, in most cases, it is advantageous to carry out a readjustment, or readjustment, using at least one comparison between a series of values from an estimated output and the same actually measured series. In this example, the linear estimator is readjusted from the surface torque.
  • the estimation technique is based on the filtering principles of Luenberger and Kalman ("Automatic linear systems" by P. De Larminat and Y. Thomas-Flammarion science; Paris IV, 1975).
  • the principle of a linear estimator can be illustrated by FIG. 4 where the measurement of the torque Cms and the estimated value Ces are compared in the means 42, the difference between these two values being injected into an adapter 43 in real time.
  • the objective here is to reconstruct the exits as faithfully as possible rather than having an exact model.
  • state registration consists in performing a weighting between the states predicted by the model at time t and the states reconstructed from only measured outputs. This weighting is not a simple average, but it takes into account the degree of precision of the estimates of the states obtained by its two channels independent.
  • the measured torque can be used as input Cms and perform the readjustment by the rotation speed.
  • the state registration technique introduces a slaving of Cms measured on these estimated.
  • the modal study of the knowledge model locates the five first modes at 0.344 Hz, 1.86 Hz, 3.61 Hz, 5.37 Hz, 7.12 Hz.
  • the final choice of reduction order can be made by comparing the results obtained with a model reduced to order 2 and a model reduced to order 4.
  • the reduced estimation model must, preferably meet the technological constraints of real time.
  • Block 50 represents a physical model representing a rotary drilling process, for example illustrated in Figure 2. This model takes into account account for a determined operating situation, in particular by receiving the mechanical characteristics of the drill string used, symbolization referenced 51, well and surface conditions, symbolization referenced 52, and friction laws, symbolized by reference 53.
  • Block 54 represents the main torsion model once linearized and reduced as described above. All these steps gathered under the DF brace run in delayed time compared to the progress of the process rotary drilling, the other steps gathered under the TR brace are executed in time real.
  • Block 55 is directly what has been called the estimator.
  • Means of measure 56 placed on the top of the drill string give the torque measurements and speed of rotation at the top of the rods, i.e. at the surface. These surface measurements are taken into account in the estimator, as described above, to give a estimation of the displacement values of the drilling tool, in particular the speed of rotation of the Vef drilling tool.
  • the present invention is advantageously implemented on a construction site drilling in order to have as precise an estimate as possible of the rotation speed of the drilling tool in real time, using only surface measurements, in particular the rotational speed of conventional means for rotating the gasket drilling, and a surface installation equipped with electronic and computer means. It is very interesting to have an estimate of the background parameters in order to detect, and even to prevent known malfunctions, for example the so-called behavior of "stick-slip” characterized by very sensitive variations in the speed of rotation of the tool at bottom while it is driven through a drill string set rotation from the surface at a substantially constant speed. Tool speed can vary between a practically zero speed and a value of the speed of rotation very higher than the speed applied to the surface. This may have consequences detrimental to the life of the tools, increasing the mechanical fatigue of the train of rods and the frequency of connection breaks.
  • Figure 6 shows a record F of the speed of rotation of a tool drilling, recording made from bottom sensors, for example using means described in document FR / 92-02273.
  • Said means allow measurements of rotational speed and surface torque, synchronized with speed measurements of rotation and torque at the bottom. Area measurements are used as inputs to the estimator according to the invention, and after having carried out all the necessary calculation steps, we obtain the record S corresponding to the speed of the tool at the bottom estimated by the estimator according to the invention.
  • the diagrams F and S are plotted with, in abscissa the time and on the ordinate the speed of rotation of the tool. We can thus compare the measurement actually performed F with the values S obtained by the estimator according to the invention.
  • the present invention is also applicable to the estimate of the instantaneous speed of rotation of elements included in the gasket drilling, elements which can be located at a certain distance from the drilling tool.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Earth Drilling (AREA)

Description

La présente invention concerne le domaine des mesures en cours de forage, en particulier des mesures concernant le déplacement d'un outil de forage fixé à l'extrémité d'un train de tiges de forage. La méthode selon l'invention propose une solution pour estimer notamment la vitesse de rotation instantanée de l'outil au fond du puits, lesdites estimations étant obtenues par le moyen d'un programme de calcul prenant en compte des mesures effectuées au sommet du train de tiges, c'est-à-dire sensiblement à la surface du sol, généralement par le moyen de capteurs ou d'un raccord instrumenté situés dans le voisinage des moyens d'entraínement en rotation de la garniture.
On connaít des techniques de mesure pour l'acquisition d'informations liées au comportement dynamique de la garniture de forage, qui utilisent un ensemble de capteurs de fond reliés à la surface par un conducteur électrique. Dans le document FR-2 688 026, il est utilisé deux ensembles de capteurs de mesure reliés par un câble du type logging, l'un étant situé au fond du puits, l'autre au sommet de la garniture de forage. Cependant, la présence d'un câble le long de la garniture de forage est gênante pour les opérations de forage proprement dites.
On connaít le document EP-A-709546 qui décrit une méthode prédictive de mesures de fond de puits consécutives à des mesures de surface, à partir d'une équation prédictive obtenue à la suite de collectes de mesures de fond et de surface. Cette méthode impose une collecte de mesures de fond et deux des moyens de collecte de fond de puits.
On connaít par les documents FR 2645205 ou FR 2666845 des dispositifs de surface placés au sommet de la garniture qui déterminent certains dysfonctionnements de forage en fonction de mesures de surface, mais sans prendre en compte, de manière physique, le comportement dynamique de la garniture et de l'outil de forage dans le puits.
Entre le fond d'un puits et la surface du sol, il existe un train de tiges le long duquel ont lieu des phénomènes dissipatifs d'énergie (frottement sur la paroi, amortissement de torsion,...), des phénomènes conservatifs de flexibilité, notamment en torsion. Il y a ainsi une distorsion entre les mesures des déplacements de fond et de surface qui dépend principalement des caractéristiques intrinsèques de la garniture (longueur, raideur, géométrie), des caractéristiques de frottement à l'interface tiges/paroi et de phénomènes aléatoires.
C'est pourquoi, les informations contenues dans les mesures de surface ne suffisent pas à elles seules à résoudre le problème posé, c'est-à-dire connaítre les déplacements instantanés de l'outil en connaissant les déplacements instantanés de la garniture en surface. Il faut compléter les informations de mesures de surface par des informations indépendantes, d'une autre nature, qui prennent en compte la structure du train de tiges et son comportement entre le fond et la surface: c'est le rôle du modèle de connaissance qui établit les relations théoriques entre le fond et la surface.
La méthodologie de la présente invention utilise la conjonction d'un tel modèle, défini a priori, et de mesures de surface acquises en temps réel.
Ainsi, la présente invention concerne une méthode d'estimation du déplacement effectif d'un outil de forage fixé à l'extrémité d'une garniture de forage et entraíné en rotation dans un puits par des moyens d'entraínement situés en surface, dans laquelle on utilise un modèle physique non linéaire du processus de forage fondé sur des équations générales de la mécanique. Dans la méthode, on effectue les étapes suivantes:
  • on identifie les paramètres dudit modèle en prenant en compte les paramètres dudit puits et de ladite garniture,
  • on linéarise ledit modèle autour d'un point de fonctionnement,
  • on réduit ledit modèle linéarisé en ne conservant que certains des modes propres de la matrice d'état dudit modèle,
  • on calcule, en temps réel, le déplacement de l'outil de forage à l'aide du modèle réduit et d'au moins un paramètre mesuré en surface.
Le modèle peut prendre en compte essentiellement les déplacements en rotation et ledit modèle réduit peut calculer en temps réel la vitesse de rotation instantanée de l'outil de forage, ledit paramètre mesuré en surface peut être pris parmi au moins: la vitesse de rotation de la garniture ou le couple.
Le modèle réduit peut être affiné par un filtrage auto adaptatif qui minimise la différence entre une mesure réelle d'un paramètre lié au déplacement de la garniture en surface et la sortie correspondante obtenue par ledit modèle réduit.
Le filtrage peut prendre en compte le couple mesuré en surface ou la vitesse de rotation mesurée en surface.
L'invention concerne également un système d'estimation du déplacement effectif d'un outil de forage fixé à l'extrémité d'une garniture de forage et entraíné en rotation dans un puits par des moyens d'entraínement situés en surface, dans lequel une installation de calcul comporte des moyens de modélisation physique non linéaire du processus de forage fondé sur des équations générales de la mécanique. Les paramètres desdits moyens de modélisation sont identifiés en prenant en compte les paramètres dudit puits et de ladite garniture, et l'installation de calcul comporte des moyens de linéarisation dudit modèle autour d'un point de fonctionnement, des moyens de réduction dudit modèle linéarisé afin de ne conserver que certains des modes propres de la matrice d'état dudit modèle, des moyens de calcul, en temps réel, du déplacement de l'outil de forage à l'aide des moyens de modélisation une fois linéarisés et réduits et des moyens de mesure d'au moins un paramètre lié au déplacement de la garniture en surface.
Les moyens de modélisation peuvent ne prendre en compte que la torsion, et les paramètres peuvent être la vitesse de rotation et/ou le couple.
La présente invention sera mieux comprise et ses avantages apparaítront clairement à la lecture de la description d'un exemple, nullement limitatif, illustré par les figures ci-après annexées, parmi lesquelles:
  • la figure 1 représente schématiquement les moyens mis en oeuvre pour une opération de forage,
  • la figure 2 représente un exemple de diagramme d'un modèle physique en torsion,
  • la figure 3 représente un diagramme d'un estimateur en boucle ouverte,
  • la figure 4 représente un diagramme d'un estimateur avec recalage,
  • la figure 5 représente schématiquement la méthodologie de la constitution de l'estimateur selon l'invention.
  • la figure 6 illustre les résultats obtenus avec l'estimateur.
La figure 1 illustre un appareil de forage sur lequel on mettra en oeuvre l'invention. L'installation de surface comprend un appareil de levage 1 comprenant une tour de levage 2, un treuil 3 qui permettent le déplacement d'un crochet de forage 4. Sous le crochet de forage sont suspendus des moyens d'entraínement 5 en rotation de l'ensemble de la garniture de forage 6 placée dans le puits 7. Ces moyens d'entraínement peuvent être du type tige d'entraínement ou kelly accouplée à une table de rotation 8 et les motorisations mécaniques, ou du type tête d'entraínement motorisée ou "power swivel" suspendue directement au crochet et guidée longitudinalement dans la tour.
La garniture de forage 6 est constituée conventionnellement par des tiges de forage 10, d'une partie 11 appelée couramment BHA pour "Bottom Hole Assembly" comportant principalement des masses-tiges, un outil de forage 12 en contact avec le terrain en cours de forage. Le puits 7 est rempli d'un fluide, dit de forage, qui circule de la surface au fond par le canal intérieur de la garniture de forage et remonte en surface par l'espace annulaire entre les parois du puits et la garniture de forage.
Pour la mise en oeuvre de l'invention, on peut intercaler un raccord instrumenté 13 entre les moyens d'entraínement et le sommet de la garniture. Ce raccord permet de mesurer le couple, la vitesse de rotation, éventuellement le déplacement longitudinal du sommet de la garniture. Ces mesures, dites de surface, sont transmises par câble ou radio vers une installation électronique d'enregistrement, de traitement, d'affichage, non représentée ici. A la place du raccord 13, on pourra utiliser d'autres capteurs tels un tachymètre sur la table de rotation pour mesurer la vitesse de rotation et un appareil de mesure du couple sur l'appareil de motorisation, si la précision des mesures ainsi obtenues est suffisante.
La partie 11 de la BHA peut plus précisément comporter, des masses-tiges, des stabilisateurs, un amortisseur de chocs, et un second raccord instrumenté 14 qui ne sera utilisé que pour contrôler expérimentalement la présente invention en permettant la comparaison entre le déplacement de l'outil de forage 12 effectivement mesuré par le raccord instrumenté 14 et le déplacement estimé grâce à la mise en oeuvre de la présente invention. Il est donc clair que l'application de la présente invention n'utilise pas de raccord instrumenté placé au fond du puits.
Le foreur qui conduit une opération de forage avec les appareils décrits sur la figure 1 a trois principales actions possibles, qui sont donc les variables de commande possibles permettant de contrôler le processus de forage. Ce sont, le poids sur l'outil qui est réglé par le treuil lequel contrôle la position du crochet, la vitesse de rotation de la table de rotation ou équivalent, et le débit de fluide de forage injecté.
Pour illustrer un exemple de la présente invention, on utilisera un modèle du système mécanique composé des éléments technologiques suivants:
  • un ensemble d'entraínement: organe de régulation et motorisation,
  • un ensemble de tiges,
  • un ensemble de masses-tiges,
  • un terrain représentant le contact outil/roche.
Le modèle décrit traitera le train de tiges comme un élément monodimensionnel vertical. Seuls les déplacements en rotation seront considérés, les déplacements verticaux et latéraux étant négligés.
La figure 2 représente le schéma-bloc du modèle; ici de torsion. C'est un modèle classique aux différences finies qui comporte plusieurs dizaines de mailles représentées par les blocs 20. Chaque maille représente une partie du train de tiges, composé entre autres de tiges de forage et de masses-tiges. Il s'agit d'un triplet masse-ressort-amortissement figuré par les schémas référencés 21, 22, 23. Chaque bloc est muni de deux entrées et sorties représentées par les couples de flèches 24 et 25 qui représentent les couples d'entrée et de sortie et les vitesses de rotation d'entrée et de sortie. Cette représentation montre la manière de connecter numériquement plusieurs tiges (ou mailles) comme on connecte physiquement les tiges de la garniture.
Le bloc 26 représente les conditions aux limites en tête du train de tiges résultant de l'appareil de forage.
Le bloc 27 représente la loi de friction du contact outil/terrain déterminant les conditions aux limites au fond du puits. Il s'agit ici d'un dissipateur de puissance de type couple résistant. Les différentes courbes 28 représentent une loi non linéaire de l'interaction outil/roche, en fonction de la vitesse instantanée de l'outil (entrée 29) du couple résistant (sortie 30) et du poids sur l'outil (entrée 31).
Ce modèle est validé en utilisant des données enregistrées sur chantier à l'aide des raccords instrumentés de fond et de surface.
Le fluide de forage et les parois du puits n'interviennent que dans la mesure où ils génèrent un couple résistant de friction. Par expérience, et en utilisant les mesures de fond et de surface, on pourra établir une loi de friction le long des tiges linéaire en fonction de vitesse de rotation.
Le modèle de torsion ainsi obtenu est généralement d'ordre élevé, c'est-à-dire de l'ordre de 50 à 100 pour reproduire la réalité avec une finesse suffisante.
Pour obtenir un modèle rapidement exécutable et robuste au changement de conditions de forage, par exemple le changement de terrains traversés, on procède aux étapes ci-après décrites.
On linéarise le modèle généralement non linéaire. Dans l'exemple ci-dessus décrit, on linéarise le modèle en choisissant un point de fonctionnement (une vitesse de rotation et un poids sur l'outil) représentatif des conditions de forage réelles. On peut vérifier que le comportement du modèle torsionel de connaissance, une fois linéarisé, est correct dans le voisinage du point de fonctionnement.
La linéarisation autour d'un point de fonctionnement consiste à calculer le Jacobien du système d'état non-linéaire. Le système d'état linéaire obtenu est alors de la forme: x = A.x + B.e s = C.x + D.e    avec:
x = X-X0
X 0= valeurs des états au point de fonctionnement
e = E-E0
E 0= valeurs des entrées au point de fonctionnement
s = S-S0
S 0= valeurs des sorties au point de fonctionnement
La mise sous forme modale se fait d'abord par un changement de base : z = P.x    z = P. x
Après résolution on obtient: z = P-1.A.P.z + P-1.B.e    z = Λ.z + BΛ.e s = C.P.z+D.e   s = CΛ.z+D.e
  • P est la matrice des vecteurs propres de la matrice d'état A,
  • Λ est la matrice diagonale des valeurs propres de la matrice d'état A.
  • Après linéarisation, le modèle de torsion conserve le même ordre élevé. L'analyse des modes propres du modèle de torsion permet de quantifier la contribution de chaque mode sur les sorties dignes d'intérêt. On ne conserve alors que les modes pertinents; c'est-à-dire ceux qui ont une influence notable sur le comportement dynamique représenté par lesdites sorties.
    Le modèle réduit doit reproduire les phénomènes dans une certaine bande de fréquences. Les critères de sélection des modes sont donc de deux ordres et reposent sur des concepts d'observabilité:
    • suppression des modes non ou peu observables sur les sorties mesurées,
    • suppression des modes hautes fréquences, n'entrant pas dans la bande de fréquence de la commande ou de l'estimateur.
    La méthode de réduction employée est la méthode des perturbations singulières. Elle consiste à garder de la matrice d'état et de la matrice de commande, les lignes et les colonnes correspondant aux modes à garder. Pour conserver les gains statiques, les modes rapides sont remplacés par leur valeur statique, ce qui a pour conséquence d'introduire une matrice directe.
    La méthode suppose que les modes rapides prennent leur équilibre en un temps négligeable, c'est-à-dire qu'ils s'établissent instantanément (hypothèse quasi-statique).
    La figure 3 montre le bloc diagramme d'un système d'estimation du type boucle ouverte. Le bloc 40 schématise les moyens de mesures de paramètres de surface, ici, la vitesse de rotation de la garniture Vms mesurée à la table ou à la tête d'injection motorisée, et le couple Cms également mesuré en surface. Le bloc 41 représente le modèle réduit qui simule le modèle physique de torsion non linéaire en calculant la fonction de transfert entre l'entrée (Vms) et les sorties Ces, Cef et Vef représentant respectivement le couple de surface estimé, le couple au fond estimé et la vitesse de fond estimée.
    Cependant la fonction de transfert est toujours une approximation de la réalité et toute désadaptation entre le modèle et le processus réel de forage peut créer une divergence entre les valeurs estimées et les valeurs réelles par intégration des écarts. Aussi, dans la plupart des cas, il est avantageux d'effectuer un réajustement, ou recalage, à l'aide d'au moins une comparaison entre une série de valeurs d'une sortie estimée et la même série réellement mesurée. Dans cet exemple, l'estimateur linéaire est recalé à partir du couple de surface.
    Selon le même principe et de façon équivalente, il sera possible d'entrer le couple mesuré en surface Cms et de recaler à partir de la comparaison entre la vitesse de rotation estimée et la vitesse de rotation mesurée en surface.
    La technique d'estimation repose sur les principes de filtrage de Luenberger et de Kalman ("Automatique des systèmes linéaires" par P. De Larminat et Y. Thomas-Flammarion Sciences; Paris IV, 1975). Le principe d'un estimateur linéaire peut être illustré par la figure 4 où la mesure du couple Cms et la valeur estimée Ces sont comparées dans les moyens 42, l'écart entre ces deux valeurs étant injecté dans un adaptateur 43 en temps réel. L'objectif est ici de reconstituer le plus fidèlement possible les sorties plutôt que d'avoir un modèle exact. C'est pourquoi on effectue un recalage d'état. Comme les sorties sont reliées directement aux états, le recalage d'état consiste à effectuer une pondération entre les états prédits par le modèle à l'instant t et les états reconstitués à partir des seules sorties mesurées. Cette pondération n'est pas une simple moyenne, mais elle prend en compte le degré de précision des estimations des états obtenus par ses deux voies indépendantes. Comme précisé plus haut, on peut utiliser comme entrée le couple mesuré Cms et effectuer le recalage par la vitesse de rotation.
    Une fois recalés les états du modèle qui représentent la dynamique du processus de forage, toutes les sorties, qu'elles soient mesurées ou non peuvent être recalculées.
    Cette estimation n'est pas seulement intéressante pour les variables non mesurées comme Cef et Vef, elle s'applique également aux variables mesurées (par exemple Cms ou Vms) qui ont servi au recalage. La valeur estimée Ces est l'équivalent d'une valeur filtrée sur la base d'un modèle: c'est pourquoi on utilise généralement le terme de filtrage (filtrage de Luenberger, filtrage de Kalman...).
    La technique de recalage d'états, telle que décrite précédemment introduit un asservissement de Cms mesuré sur Ces estimé.
    Ce bouclage supprime le risque de divergence mentionné ci-dessus, lorsque le modèle est simulé en boucle ouverte (figure 3).
    Il y a ainsi une désensibilisation des variables estimées vis à vis des imperfections du modèle. Dans ce contexte, on n'a plus besoin d'avoir un modèle parfait : un modèle approché est suffisant.
    En outre, on ne dispose ici que d'une mesure, le couple C, pour effectuer le recalage : il ne paraít pas possible de recaler un grand nombre d'états à partir de cette mesure. C'est pourquoi, le modèle de torsion non linéaire ne convient pas en dépit de sa plus grande précision.
    Il existe donc un compromis à effectuer entre la précision et l'ordre du système. Il faut rechercher le modèle d'ordre minimum qui respecte les tolérances de précision souhaitables, et qui soit également facile à régler et robuste.
    Le choix de l'ordre du modèle réduit dépend des critères qualitatifs suivants :
    • il faut sauvegarder les modes propres de vibration en torsion qui sont prépondérants dans les sorties à ré estimer ;
    • pour des raisons de cohérence et de stabilité numérique, il faut rejeter les modes de fréquences élevées supérieures à fmax = fe/2 où fe est la fréquence d'échantillonnage des entrées et des sorties. On peut préciser ces notions en prenant par exemple le cas d'un train de tiges composée de 800m de tiges 5" (127 mm) et de 200m de masses-tiges 8" (203,2 mm) pour lequel le système d'acquisition a une fréquence d'échantillonnage de 10 Hz.
    L'étude modale du modèle de connaissance, discrétisé en 35 mailles, situe les cinq premiers modes à 0.344 Hz, 1.86 Hz, 3.61 Hz, 5.37 Hz, 7.12 Hz.
    fmax étant égal à la moitié de la fréquence d'échantillonnage, fmax= 5 Hz. Mais compte tenu des possibilités des filtres anti-repliement, la fréquence maximale observable est de 2,8 Hz.
    Il est donc superflu de choisir un modèle réduit d'ordre supérieur à 4 (deux modes) si on veut intégrer le modèle à la cadence d'échantillonnage fe.
    Le choix définitif de l'ordre de réduction peut se faire par comparaison des résultats obtenus avec un modèle réduit à l'ordre 2 et un modèle réduit à l'ordre 4.
    De plus, il ne faut pas oublier que le modèle réduit d'estimation doit, de préférence, satisfaire les contraintes technologiques du temps réel.
    L'estimateur est donc construit selon les étapes suivantes :
    • génération du modèle réduit,
    • implantation des filtres passe-haut,
    • agrégation des filtres passe-haut et du modèle réduit, l'ensemble devient le modèle d'estimation,
    • calcul des gains de recalage,
    • construction de l'estimateur complet.
    La méthodologie pour la construction de l'estimateur selon l'invention peut être illustrée par la figure 5. Le bloc 50 représente un modèle physique représentant un processus de forage rotary, par exemple illustré par la figure 2. Ce modèle prend en compte une situation de fonctionnement déterminée en recevant notamment les caractéristiques mécaniques de la garniture de forage utilisée, symbolisation référencée 51, les conditions de puits et de surface, symbolisation référencée 52, et des lois de friction, symbolisation référencée 53. Le bloc 54 représente le modèle de torsion principal une fois linéarisé et réduit selon la description ci-dessus. Toutes ces étapes rassemblées sous l'accolade DF s'exécutent en temps différé par rapport au déroulement du processus de forage rotary, les autres étapes rassemblées sous l'accolade TR sont exécutées en temps réel.
    Le bloc 55 est directement ce que l'on a appelé l'estimateur. Des moyens de mesure 56 placés sur le sommet de la garniture de forage donnent les mesures de couple et vitesse de rotation au sommet des tiges, c'est-à-dire en surface. Ces mesures de surface sont prises en compte dans l'estimateur, comme décrit plus haut, pour donner une estimation des valeurs de déplacement de l'outil de forage, en particulier la vitesse de rotation de l'outil de forage Vef.
    La présente invention est avantageusement mise en oeuvre sur un chantier de forage afin d'avoir une estimation aussi précise que possible de la vitesse de rotation de l'outil de forage en temps réel, et cela à partir des seules mesures de surface, notamment la vitesse de rotation des moyens conventionnels de mise en rotation de la garniture de forage, et d'une installation de surface équipée de moyens électroniques et informatiques. Il est très intéressant d'avoir une estimation des paramètres de fond de façon à détecter, et même à prévenir des dysfonctionnements connus, par exemple le comportement dit de "stick-slip" caractérisé par des variations très sensibles de la vitesse de rotation de l'outil au fond alors que celui-ci est entraíné par l'intermédiaire d'une garniture de forage mise en rotation à partir de la surface à une vitesse sensiblement constante. La vitesse de l'outil peut varier entre une vitesse pratiquement nulle et une valeur de la vitesse de rotation très supérieure à la vitesse appliquée en surface. Cela peut avoir pour conséquences des effets néfastes sur la durée de vie des outils, sur l'augmentation de la fatigue mécanique du train de tiges et la fréquence des ruptures des connexions.
    La figure 6 montre un enregistrement F de la vitesse de rotation d'un outil de forage, enregistrement effectué à partir de capteurs de fond, par exemple à l'aide des moyens décrits dans le document FR/92-02273. Lesdits moyens permettent des mesures de la vitesse de rotation et du couple en surface, synchronisées avec des mesures de la vitesse de rotation et du couple au fond. Les mesures de surface sont utilisées comme entrées dans l'estimateur selon l'invention, et après avoir exécuté toutes les étapes de calcul nécessaires, on obtient l'enregistrement S correspondant à la vitesse de l'outil au fond estimée par l'estimateur selon l'invention. Sur la figure 6, les diagrammes F et S sont tracés avec, en abscisse le temps et en ordonnée la vitesse de rotation de l'outil. On peut ainsi comparer la mesure réellement effectuée F avec les valeurs S obtenues par l'estimateur selon l'invention.
    Il est clair que le dysfonctionnement de "stick-slip" qui apparaít au temps 80 s est distinctement estimé sur le diagramme S, ce qui permet au foreur d'être averti de l'existence de ce dysfonctionnement grâce à l'estimateur et sans avoir besoin de mesures de fond en temps réel. On peut remarquer sur cette figure, qu'entre 0 et 20 secondes, il y a une mauvaise estimation qui s'explique par le fait de l'initialisation du filtrage auto adaptatif. Cette constatation montre tout l'intérêt de la boucle de recalage.
    Il faut également comprendre que la présente invention est également applicable à l'estimation de la vitesse de rotation instantanée d'éléments inclus dans la garniture de forage, éléments qui peuvent être situés à une certaine distance de l'outil de forage.

    Claims (5)

    1. Méthode d'estimation du déplacement effectif d'un outil de forage (12) fixé à l'extrémité d'une garniture de forage (6) et entraíné en rotation dans un puits (7) par des moyens d'entraínement (8) situés en surface, dans laquelle on utilise un modèle physique (50) non linéaire du processus de forage fondé sur des équations générales de la mécanique, caractérisée en ce que l'on effectue les étapes suivantes:
      on identifie les paramètres dudit modèle en prenant en compte les paramètres dudit puits et de ladite garniture (51, 52, 53),
      on linéarise ledit modèle autour d'un point de fonctionnement,
      on réduit ledit modèle linéarisé en ne conservant que certains des modes propres de la matrice d'état dudit modèle,
      on calcule, en temps réel, le déplacement de l'outil de forage à l'aide du modèle réduit (54) et d'au moins un paramètre mesuré en surface, et en ce que ledit modèle prend en compte essentiellement les déplacements en rotation et en ce que ledit modèle réduit calcule en temps réel la vitesse de rotation instantanée de l'outil de forage, ledit paramètre mesuré en surface est au moins pris parmi: la vitesse de rotation de la garniture ou le couple mesuré en surface.
    2. Méthode selon la revendication 1, dans laquelle le modèle réduit est affiné par un filtrage auto adaptatif qui minimise la différence entre une mesure réelle d'un paramètre lié au déplacement de la garniture en surface et la sortie correspondante obtenue par ledit modèle réduit.
    3. Méthode selon la revendication 2, dans laquelle ledit filtrage prend en compte le couple mesuré en surface.
    4. Méthode selon la revendication 3, dans laquelle ledit filtrage prend en compte la vitesse mesurée en surface.
    5. Système d'estimation du déplacement effectif d'un outil de forage (12) fixé à l'extrémité d'une garniture de forage (6) et entraíné en rotation dans un puits (7) par des moyens d'entraínement (8) situés en surface, dans lequel une installation de calcul (9) comporte des moyens de modélisation physique non linéaire du processus de forage fondé sur des équations générales de la mécanique, en ce que des paramètres desdits moyens de modélisation sont identifiés en prenant en compte les paramètres dudit puits et de ladite garniture, en ce que l'installation de calcul comporte des moyens de linéarisation dudit modèle autour d'un point de fonctionnement, des moyens de réduction dudit modèle linéarisé afin de ne conserver que certains des modes propres de la matrice d'état dudit modèle, des moyens de calcul, en temps réel, du déplacement de l'outil de forage à l'aide des moyens de modélisation une fois linéarisés et réduits et des moyens de mesure d'au moins un paramètre lié au déplacement de la garniture en surface (13) : vitesse de rotation et/ou le couple, et dans lequel les moyens de modélisation ne prennent en compte que la torsion.
    EP97401297A 1996-06-24 1997-06-09 Méthode et système d'estimation en temps reel d'au moins un paramètre lié au déplacement d'un outil de forage Expired - Lifetime EP0816629B1 (fr)

    Applications Claiming Priority (2)

    Application Number Priority Date Filing Date Title
    FR9607915 1996-06-24
    FR9607915A FR2750160B1 (fr) 1996-06-24 1996-06-24 Methode et systeme d'estimation en temps reel d'au moins un parametre lie au deplacement d'un outil de forage

    Publications (2)

    Publication Number Publication Date
    EP0816629A1 EP0816629A1 (fr) 1998-01-07
    EP0816629B1 true EP0816629B1 (fr) 2003-05-14

    Family

    ID=9493416

    Family Applications (1)

    Application Number Title Priority Date Filing Date
    EP97401297A Expired - Lifetime EP0816629B1 (fr) 1996-06-24 1997-06-09 Méthode et système d'estimation en temps reel d'au moins un paramètre lié au déplacement d'un outil de forage

    Country Status (5)

    Country Link
    US (1) US5852235A (fr)
    EP (1) EP0816629B1 (fr)
    CA (1) CA2209056C (fr)
    FR (1) FR2750160B1 (fr)
    NO (1) NO972932L (fr)

    Families Citing this family (21)

    * Cited by examiner, † Cited by third party
    Publication number Priority date Publication date Assignee Title
    EP0870899A1 (fr) * 1997-04-11 1998-10-14 Shell Internationale Researchmaatschappij B.V. Ensemble de forage avec une tendance réduite de stick-slip
    US6327539B1 (en) * 1998-09-09 2001-12-04 Shell Oil Company Method of determining drill string stiffness
    US6347292B1 (en) * 1999-02-17 2002-02-12 Den-Con Electronics, Inc. Oilfield equipment identification method and apparatus
    RU2183266C1 (ru) * 2000-09-27 2002-06-10 Общество с ограниченной ответственностью "ЮганскНИПИнефть" Способ определения работоспособности породоразрушающего инструмента
    CA2724453C (fr) * 2008-06-17 2014-08-12 Exxonmobil Upstream Research Company Procedes et systemes permettant d'attenuer les vibrations de forage
    EA033087B1 (ru) 2008-11-21 2019-08-30 Эксонмобил Апстрим Рисерч Компани Способ и система для моделирования, проектирования и проведения буровых работ, которые учитывают вибрации
    BRPI0917046B1 (pt) 2008-12-02 2020-11-10 National Oilwell Varco, L.P. método para estimar a velocidade rotacional instantânea de uma estrutura inferior do poço
    WO2010063982A1 (fr) * 2008-12-02 2010-06-10 National Oilwell Varco, L.P. Procédé et appareil de réduction d'un glissement saccadé
    WO2011017627A1 (fr) 2009-08-07 2011-02-10 Exxonmobil Upstream Research Company Procédés pour estimer des indices de vibrations de forage de fond de trou à partir d'une mesure de surface
    EA201270258A1 (ru) 2009-08-07 2012-09-28 Эксонмобил Апстрим Рисерч Компани Способы оценки амплитуды вибраций на забое при бурении по результатам измерений на поверхности
    MX2016002003A (es) * 2013-09-17 2016-05-18 Halliburton Energy Services Inc Eliminacion de vibraciones de atascamiento y deslizamiento en un ensamblaje de perforacion.
    CN103941632B (zh) * 2014-04-29 2016-08-24 西北工业大学 一种基于自适应法抑制钻杆粘滑振动的系统及方法
    CA2980277C (fr) * 2015-04-29 2023-07-18 Halliburton Energy Services, Inc. Systemes et procedes pour estimation d'etat sans capteur, estimation de perturbation, et adaptation de modele pour des systemes de forage orientables rotatifs
    WO2017127045A1 (fr) 2016-01-19 2017-07-27 Halliburton Energy Services, Inc. Procédé de minimisation de réponse d'outil pour opérations de diagraphie de fond
    US10851639B2 (en) 2017-03-31 2020-12-01 Exxonmobil Upstream Research Company Method for drilling wellbores utilizing a drill string assembly optimized for stick-slip vibration conditions
    US11536128B2 (en) 2017-03-31 2022-12-27 Exxonmobil Upstream Research Company Method for drilling wellbores utilizing drilling parameters optimized for stick-slip vibration conditions
    US10648321B2 (en) 2017-04-04 2020-05-12 Nabors Drilling Technologies Usa, Inc. Surface control system adaptive downhole weight on bit/torque on bit estimation and utilization
    US10370902B2 (en) 2017-04-04 2019-08-06 Nabors Drilling Technologies Usa, Inc. Downhole steering control apparatus and methods
    CN111655969B (zh) 2017-12-23 2023-08-25 诺伊蒂克技术股份有限公司 使用实时测量和建模来优化管子运行操作的系统和方法
    US11704453B2 (en) * 2019-06-06 2023-07-18 Halliburton Energy Services, Inc. Drill bit design selection and use
    US11748531B2 (en) 2020-10-19 2023-09-05 Halliburton Energy Services, Inc. Mitigation of high frequency coupled vibrations in PDC bits using in-cone depth of cut controllers

    Family Cites Families (12)

    * Cited by examiner, † Cited by third party
    Publication number Priority date Publication date Assignee Title
    US4845628A (en) * 1986-08-18 1989-07-04 Automated Decisions, Inc. Method for optimization of drilling costs
    US4794535A (en) * 1986-08-18 1988-12-27 Automated Decisions, Inc. Method for determining economic drill bit utilization
    FR2614360B1 (fr) * 1987-04-27 1989-06-16 Forex Neptune Procede de mesure de la vitesse d'avancement d'un outil de forage
    FR2645205B1 (fr) * 1989-03-31 1991-06-07 Elf Aquitaine Dispositif de representation auditive et/ou visuelle des phenomenes mecaniques dans un forage et utilisation du dispositif dans un procede de conduite d'un forage
    FR2666845B1 (fr) * 1990-09-14 1997-01-10 Elf Aquitaine Procede de conduite d'un forage.
    GB2264562B (en) * 1992-02-22 1995-03-22 Anadrill Int Sa Determination of drill bit rate of penetration from surface measurements
    FR2688026B1 (fr) * 1992-02-27 1994-04-15 Institut Francais Petrole Systeme et methode d'acquisition de donnees physiques liees a un forage en cours.
    US5305836A (en) * 1992-04-08 1994-04-26 Baroid Technology, Inc. System and method for controlling drill bit usage and well plan
    GB9216740D0 (en) * 1992-08-06 1992-09-23 Schlumberger Services Petrol Determination of drill bit rate of penetration from surface measurements
    GB9218836D0 (en) * 1992-09-05 1992-10-21 Schlumberger Services Petrol Method for determining weight on bit
    NO315670B1 (no) * 1994-10-19 2003-10-06 Anadrill Int Sa Fremgangsmåte og anordning for måling av boretilstander ved kombinasjon avnedihulls- og overflatemålinger
    US5581024A (en) * 1994-10-20 1996-12-03 Baker Hughes Incorporated Downhole depth correlation and computation apparatus and methods for combining multiple borehole measurements

    Also Published As

    Publication number Publication date
    NO972932L (no) 1997-12-29
    FR2750160B1 (fr) 1998-08-07
    US5852235A (en) 1998-12-22
    CA2209056A1 (fr) 1997-12-24
    NO972932D0 (no) 1997-06-23
    FR2750160A1 (fr) 1997-12-26
    CA2209056C (fr) 2006-11-14
    EP0816629A1 (fr) 1998-01-07

    Similar Documents

    Publication Publication Date Title
    EP0816629B1 (fr) Méthode et système d'estimation en temps reel d'au moins un paramètre lié au déplacement d'un outil de forage
    EP1046781B1 (fr) Méthode et système de détection du déplacement longitudinal d'un outil de forage
    EP0816630B1 (fr) Méthode et système d'estimation en temps réel d'au moins un paramètre lié au comportement d'un outil de fond de puits
    FR3027339A1 (fr)
    EP3936669B1 (fr) Procédé, système et produit-programme informatique pour déterminer des propriétés du sol
    FR2869067A1 (fr) Systeme et procede de synthese de champ pour l'optimisation d'un dispositif de forage
    FR3043481A1 (fr)
    FR2570757A1 (fr) Procede et dispositif pour estimer les caracteristiques de formation de la formation exposee au fond d'un trou
    FR3077142A1 (fr) Optimisation d’une longueur de référence à des fins de préservation de signal et traitement d’une longueur de référence à des fins de détection de vibration répartie
    EP0125164A1 (fr) Méthode de détermination des caractéristiques d'une formation souterraine produisant un fluide
    US10190406B2 (en) Formation fracturing potential using surrounding pore pressures
    WO2012038647A1 (fr) Procede d'estimation de parametres elastiques par inversion de mesures sismiques 4d
    FR3039196A1 (fr)
    FR3066224A1 (fr) Procede et systeme pour forer un puits de forage et identifier une defaillance du trepan de forage par deconvolution des donnees de capteur
    FR3084102A1 (fr) Ajustement du fonctionnement d’un outil de puits pour manipuler le taux de pénétration (rop) d’un trépan de forage sur la base de multiples prévisions de rop
    FR2611804A1 (fr) Procede de controle des operations de forage d'un puits
    FR2963111A1 (fr) Procede d'estimation de parametres elastiques par inversion de mesures sismiques 4d
    FR3152900A1 (fr) Procédé, appareil, et système d’estimation de coordonnées d’une pointe de dent de godet, excavatrice, et support de stockage
    FR3035147A1 (fr)
    FR3034191B1 (fr) Determination de pression de formation
    Wiktorski et al. Comparative study of surface and downhole drillstring vibrations measurements on a laboratory-scale drilling rig
    FR3035145A1 (fr)
    EP1341094B1 (fr) Procédé pour la détermination d'un indice de qualité spatiale de données regionalisées
    FR3076626A1 (fr) Prediction de proprietes de roches sismiques en temps avance sur la base d'une analyse sismique 4d
    EP3502751A1 (fr) Procédé de modelisation d'un bassin sédimentaire

    Legal Events

    Date Code Title Description
    PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

    Free format text: ORIGINAL CODE: 0009012

    AK Designated contracting states

    Kind code of ref document: A1

    Designated state(s): GB IT NL

    17P Request for examination filed

    Effective date: 19980707

    AKX Designation fees paid

    Free format text: GB IT NL

    RBV Designated contracting states (corrected)

    Designated state(s): GB IT NL

    REG Reference to a national code

    Ref country code: DE

    Ref legal event code: 8566

    17Q First examination report despatched

    Effective date: 20010313

    GRAH Despatch of communication of intention to grant a patent

    Free format text: ORIGINAL CODE: EPIDOS IGRA

    GRAH Despatch of communication of intention to grant a patent

    Free format text: ORIGINAL CODE: EPIDOS IGRA

    GRAH Despatch of communication of intention to grant a patent

    Free format text: ORIGINAL CODE: EPIDOS IGRA

    GRAA (expected) grant

    Free format text: ORIGINAL CODE: 0009210

    AK Designated contracting states

    Designated state(s): GB IT NL

    REG Reference to a national code

    Ref country code: GB

    Ref legal event code: FG4D

    Free format text: NOT ENGLISH

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: GB

    Payment date: 20030602

    Year of fee payment: 7

    GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)

    Effective date: 20030515

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: NL

    Payment date: 20030627

    Year of fee payment: 7

    PLBE No opposition filed within time limit

    Free format text: ORIGINAL CODE: 0009261

    STAA Information on the status of an ep patent application or granted ep patent

    Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

    26N No opposition filed

    Effective date: 20040217

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: GB

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20040609

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: NL

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20050101

    GBPC Gb: european patent ceased through non-payment of renewal fee

    Effective date: 20040609

    NLV4 Nl: lapsed or anulled due to non-payment of the annual fee

    Effective date: 20050101

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: IT

    Payment date: 20070622

    Year of fee payment: 11

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: IT

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20080609