EP2742365B1 - Kompensiertes cross-well-tomografieverfahren und systeme - Google Patents

Kompensiertes cross-well-tomografieverfahren und systeme Download PDF

Info

Publication number
EP2742365B1
EP2742365B1 EP11873738.6A EP11873738A EP2742365B1 EP 2742365 B1 EP2742365 B1 EP 2742365B1 EP 11873738 A EP11873738 A EP 11873738A EP 2742365 B1 EP2742365 B1 EP 2742365B1
Authority
EP
European Patent Office
Prior art keywords
receivers
transmitters
data
compensated value
ikjm
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.)
Not-in-force
Application number
EP11873738.6A
Other languages
English (en)
French (fr)
Other versions
EP2742365A1 (de
EP2742365A4 (de
Inventor
Burkay Donderici
Baris GUNER
Michael S. Bittar
Luis E. San Martin
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.)
Halliburton Energy Services Inc
Original Assignee
Halliburton Energy Services Inc
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 Halliburton Energy Services Inc filed Critical Halliburton Energy Services Inc
Publication of EP2742365A1 publication Critical patent/EP2742365A1/de
Publication of EP2742365A4 publication Critical patent/EP2742365A4/de
Application granted granted Critical
Publication of EP2742365B1 publication Critical patent/EP2742365B1/de
Not-in-force legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/42Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling

Definitions

  • crosswell tomography was performed using seismic transmitters and receivers, but more recently the focus has been on the use of electromagnetic (EM) transmitters and receivers.
  • EM electromagnetic
  • noise and inaccuracies in the survey system will negatively impact image quality.
  • One important cause of this degradation in crosswell tomography is mismatches between assumed and actual values for sensor gains and phases. These mismatches can result from a number of causes including differences in sensor types, calibration errors, variation in borehole configurations (e.g., mud composition and borehole radius), and environmental parameters such as temperature and pressure. Phase mismatches can be caused or exacerbated by latencies in electronic systems and synchronization errors for the widely-distributed components of a typical crosswell tomography system. Traditional methods crosswell tomography do not appear to adequately address such inaccuracies.
  • a method for improving image quality of crosswell tomography comprising: obtaining data from at least two receivers in response to transmissions from at least two transmitters, wherein the transmitters and receivers are distributed among two or more boreholes; deriving at least one compensated value that combines the data from the at least two receivers' responses to the at least two transmitters as a ratio that cancels errors between actual and assumed characteristics of the at least two receivers and the at least two transmitters; and performing an inversion based at least in part on said compensated value to obtain a sub surface image, wherein said compensated value compensates for inaccuracies including sensor gains and phase variations, wherein said compensated value is expressible by V ij ′ V km ′ V im ′ V kj ′ , wherein i and k are receiver indices and j and m are transmitter indices, and wherein said inversion employs a forward model that relates subsurface formation properties to said at least one
  • a system for improving image quality of crosswell tomography comprising: a memory unit that stores inversion software; and a processor coupled to said memory to execute the software, wherein said software configures said processor to: obtain data from at least two receivers in response to transmissions from at least two transmitters, wherein one of said receivers or transmitters are in a borehole; derive at least one compensated value based on responses of said receivers to each of said transmitters that are combined as a ratio that cancels errors between actual and assumed characteristics of said receivers and said transmitters; and perform an inversion based at least in part on said compensated value to obtain a subsurface image, wherein said compensated value compensates for inaccuracies including sensor gains and phase variations, wherein said compensated value is expressible by V ij ′ V km ′ V im ′ V kj ′ , wherein i and k are receiver indices and j and m are transmitter indices, wherein the in
  • data is obtained from at least two receivers in response to each of at least two transmitters, the transmitters and receivers being located in boreholes near the region of interest.
  • the data is combined to form at least one compensated value that compensates for any transmitter and receiver gain or phase inaccuracies.
  • An inversion is then performed to identify the subsurface structure that best matches these compensated values.
  • At least some system embodiments implement this method using software that then outputs one or more images derived from the subsurface structure.
  • the disclosed systems and methods employ at least one borehole in or near the subsurface region of interest. Such boreholes may be drilled as illustrated in Fig. 1 .
  • a drilling platform 2 is equipped with a derrick 4 that supports a hoist 6 for raising and lowering a drill string 8.
  • the hoist 6 suspends a top drive 10 that is used to rotate the drill string 8 and to lower the drill string through the well head 12. Sections of the drill string 8 are connected by threaded connectors 7.
  • Connected to the lower end of the drill string 8 is a drill bit 14 that rotates to create and extend a borehole 16 along a desired path 18.
  • a pump 20 circulates drilling fluid through a supply pipe 20 to top drive 10, downhole through the interior of drill string 8, through orifices in drill bit 14, back to the surface via the annulus around drill string 8, and into a retention pit 24.
  • the drilling fluid transports cuttings from the borehole into the pit 24 and aids in maintaining the integrity of the borehole 20.
  • Logging instruments 26 may be positioned on the drill string 8.
  • a telemetry module 28 may provide communication between the logging instruments 26 and a surface interface 34 using any one of various available telemetry technologies.
  • acoustic sensors 30 and 32 detect mud pulse telemetry carrying data from the logging instruments.
  • a data processing system 38 is coupled to the surface interface 34 by a wired or wireless communications link 36 to control and acquire measurements from the logging instruments 26.
  • the illustrated data processing system 38 is a general purpose computer with one or more input devices 42 and one or more output devices 44 for interacting with a user.
  • Software on information storage media 40 (or in the computer's internal memory) configures the computer's internal processor(s) to acquire crosswell tomography measurements and derive a model of the subsurface structure for display to the user.
  • Logging instruments 26 preferably include acoustic transducers for sending and/or receiving acoustic signals, or electromagnetic antennas for sending and/or receiving electromagnetic ("EM") signals. Though illustrated on a drill string, the logging instruments can additionally be conveyed on a wire line or a coiled tubing string. Data processing system 38 acquires measurements using multiple transmitters and receivers in various positions in and around the region of interest.
  • EM electromagnetic
  • Fig. 2A shows an illustrative configuration in which various transmitters ("T") and receivers ("R”) are positioned in each of the branches of a multilateral well, including a substantially vertical branch penetrating the region of interest and five substantially horizontal branches extending outward in various directions along the top of the reservoir.
  • Fig. 2B shows a different illustrative configuration in which multiple transceiver (transmitter and receiver) nodes are positioned throughout a region of interests using multiple vertical wells, some of which may include lateral branches having additional transceiver nodes.
  • Some or all of the boreholes may be open (i.e., uncased). Similarly, some or all of the boreholes may be cased and the sensors positioned inside the casing. Where the casing causes too much attenuation, the sensors may be configured to send/receive low frequency magnetic waves.
  • Some systems may employ additional sensors positioned at the surface.
  • the number and position of the nodes depends on the region's volume, properties (resistivity or acoustic impedance), and the desired resolution of the model. Generally, however, a greater number of nodes offers better model accuracy, and a three dimensional arrangement of nodes may be preferred to a two dimensional arrangement, which in turn may be preferred to a one dimensional arrangement. Where some of the nodes are spaced too far from the region of interest to be useful, they may be omitted to avoid introducing unnecessary noise into the inversion process.
  • Fig. 3 is a function block diagram of a compensated crosswell tomography system.
  • Each of the nodes in Figs. 2A-2B includes at least one sensor 302 for transmitting and/or receiving acoustic or EM signals.
  • the sensors are antennas or magnetometers.
  • Suitable antenna configurations include a coil antenna (see, e.g., Fig. 4A ) with an optional bucking coil (see, e.g., Fig. 4B ), a toroidal antenna, an electrical gap antenna, a spaced-electrode antenna, and a horn antenna.
  • the suitable sensors include piezoelectric transducers, magnetic coil speakers, bender bars, and conventional seismic sources.
  • Some of the sensors are coupled to a transmitter 304 to transmit signals into the formation, some are coupled to a receiver 306 to receive signals from the formation, some are switchable between transmit and receive mode. In the latter case, a switch 308 may be provided between the sensor 302 and its associated transmitter 304 and receiver 306. Some embodiments may even provide multiple transmitters and receivers for at given sensor, e.g., for different signal frequencies. Conversely, the system may include multiple sensors 302 coupled to a given transmitter 304 or receiver 306, e.g., to enable the transmitter or receiver to send or receive signals in different directions or polarizations. Tilted antennas are contemplated for making measurements of orthogonal or linearly independent EM signal components.
  • Each node includes a telemetry unit 310 that controls the associated transmitters 304 and receivers 306 (and switches 308) in accordance with commands from a system control center 314.
  • a network 312 couples the control center 314 to the various nodes via conventional wired or wireless communication links to the well head of each well in the system, and an interface unit at each well head completes the communications chain via an intrawell telemetry channel such as a wireline, an electromagnetic telemetry link, an acoustic telemetry link, or a mud pulse telemetry technique. Any suitable configuration can be used for network 312.
  • control center 314 iterates through each of the available transmitters 304, causing them to fire in turn while all of the other nodes are set to detect and capture any resulting receive signal.
  • the telemetry units 310 capture the receiver signals for each transmitter firing and communicate them to the control center 314. Such measurements may be repeated and combined to improve signal to noise ratio.
  • the data processing unit 38 acquires the set of receiver responses to each transmitter firing and processes them as outlined below to produce a model of the subsurface structure. A user can view and interact with the model via a visualization unit 44.
  • Fig. 5 shows a high-level block diagram of an inversion process.
  • Raw data (the amplitude and phase of each receiver's responses to each transmitter) is processed in block 502 to obtain compensated data.
  • an iterative inversion process adapts a structural model of the region of interest until the compensated data are matched by predicted data.
  • a forward model block 506 generates the predicted data based on the current structural model, and the inversion block responsively adapts the structural model until the error between the predicted and compensated data values falls below a threshold.
  • the adapted structural model is then output by block 504 for visualization and/or analysis to determine the amount and distribution of fluids in the reservoir.
  • the compensated data is formed by constructing ratios of the raw data. These ratios are designed to compensate for any systematic errors in the transmit/receive characteristics of the various sensors. Because of this compensation, the results of the inversion process are much improved relative to those of an inversion process performed directly on the raw data. Moreover, this approach relaxes system design constraints, enabling the use of a variety of sensor types, and operation in diverse borehole and environmental circumstances. The synchronization requirements are relaxed, and the system can more easily accommodate miscalibration or parameter drift of the sensors.
  • the system includes a total of T transmit antennas and a total of R receive antennas, there are a total of T !/(2( T -2)! transmitter pairs and a total of R !/(2( R -2)!) receiver pairs, yielding up to T !R!/(4(T-2)!( R -2)! compensated values.
  • the measurements may be made at multiple frequencies, further increasing the number of measurements to be matched in the inversion process.
  • Fig. 6 is a more detailed block diagram of an inversion process.
  • the raw data is denoted as v t,r,f where t , r and f in the subscript represent a particular transmitter, receiver and frequency, respectively.
  • t , r and f in the subscript represent a particular transmitter, receiver and frequency, respectively.
  • p t1,t2,r1,r2,f the compensated data
  • some embodiments normalize the raw data with respect to measurements predicted by an assumed (perhaps constant resistivity) background model. Such normalization may be useful for reducing the dominance of strong signal channels in the inversion process.
  • a parallel track in Fig. 6 begins with an initial model of the subsurface structure in the region of interest, represented in the figure by the conductivity values ⁇ i , j , k 1 , where the superscript indicates an iteration number and i, j, k represent spatial indices in the x, y, and z directions for cells in the conductivity grid.
  • the conductivity grid needs to be chosen large enough to cover the region of interest and other features that might affect the measurements.
  • a forward model 604 operates on the model to predict the signal measurements for every transmitter-receiver combination and frequency using an initial guess for the formation conductivity.
  • the predicted data due to this conductivity profile (denoted as v F t,r,f ) are processed by a block 606 that performs the same operations as block 602.
  • the output of block 606 is denoted p t 1 , t 2 , r 1 , r 2 F .
  • Block 608 determines the difference between the outputs of blocks 602 and 608, producing p DIF t1,t2,r1,r2,f .
  • Block 610 evaluates the magnitude of this difference to determine whether the model has converged.
  • convergence is determined by comparing the magnitude of the difference to the magnitude of compensated data p t1,t2,r1,r2,f , e.g., testing whether ⁇ P DIF ⁇ ⁇ c ⁇ P ⁇ , where P DIF and P are vectors whose elements are producing p DIF t1,t2,r1,r2,f . and p t1,t2,r1,r2,f , respectively and c is a constant. If convergence is satisfied, the block 611 outputs the conductivity model.
  • block 612 derives a set of changes for the conductivity model based on p DIF
  • block 614 applies the changes to the conductivity model to arrive at an updated model ⁇ n + 1 i,j,k for iteration step n.
  • Fig. 7 shows a top view of a simulation model for comparing the performance of the disclosed method to a tomography method that fails to compensate the raw data.
  • the simulation model is a formation bed with nine vertical boreholes each having a transceiver positioned 70 meters above the bottom of the formation bed. The boreholes are evenly spaced apart by 500 meters in X and Y directions.
  • the formation bed is assumed to have two resistive anomalies in an otherwise homogeneous medium of 1 ⁇ -m resistivity. Both anomalies have a thickness of 65 m in z-direction. If the middle transceiver is taken as the origin, the first anomaly is the volume defined between 155m and 285 m in x-direction, and -130 to 130 m. in y-direction.
  • This anomaly has a resistivity of 25 ⁇ -m.
  • the second anomaly is defined between -290 and -150 m in x-direction and 120 and 330 m in y-direction. Resistivity of this anomaly is taken as 20 ⁇ -m.
  • a multi-frequency system with z-directed coil transmitters was considered with operating frequencies of 1, 2, 5 and 10 Hz.
  • Synthetic data was created with a 3-D finite difference time domain (FDTD) code according to the conductivity profile described above. A 5% multiplicative gain noise was applied to the synthetic data.
  • FDTD finite difference time domain
  • Fig. 8 The results from the traditional inversion of this data are shown in Fig. 8 .
  • the true locations of the anomalies are indicated with dashed lines.
  • the inverted image exhibits severe artifacts due to noise. Although several peaks close to the targets are observed, results also show many false positives and it is not possible to reliably determine the shape and boundaries of the reservoir.
  • Fig. 9 The results from inversion using compensated values in accordance with the method disclosed above are shown in Fig. 9 .
  • the conductivity image is unaffected by the presence of sensor gain and phase errors. It can be seen that higher intensity regions match very accurately with reservoir boundaries and any false positives occurring due to sensor gain and phase noises are eliminated.
  • Fig. 10 shows an illustrative flow diagram for compensated crosswell tomography.
  • data is obtained from at least two receivers in response to transmissions from at least two transmitters. At least one of the receivers or transmitters are located inside a borehole. The transmitters and receivers can also be distributed among two or more boreholes and on the surface.
  • data is obtained from additional pairs of receivers in response to additional pairs of transmitters.
  • at least one compensated value is derived based on the responses of the receivers to the transmitters.
  • the compensated value is expressible as: V ij ⁇ V km ⁇ V im ⁇ V kj ⁇ , where i and k are receiver indices and j and m are transmitter indices.
  • An inversion is performed based at least in part on the derived compensated value(s), as shown in block 1008.
  • a sub-surface image of a target area is produced based on the inversion.
  • the image can provide information about at least one formation property, such as acoustic velocity, attenuation, seismic reflectivity, electromagnetic resistivity, or electrical permittivity.
  • the operations represented in Fig. 10 can be implemented as software on a general purpose computer.
  • the processor(s) in such a computer are configured to execute software stored on disk or in memory, the software causing the processors to obtain the data, process it, and display resulting images in accordance with user commands.
  • the compensated measurement formula is not the only one that yields cancellation of the gain and phase errors, as combinations of measurements by three receivers in response to three transmitters, or four receivers in response to four transmitters, etc, will yield similar cancellations.
  • the following claims embrace such alternative forms, equivalents, and modifications.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Claims (11)

  1. Computerimplementiertes
    Verfahren zum Verbessern der Bildqualität von Crosswell-Tomografie, umfassend:
    Erlangen von Daten von wenigstens zwei Empfängern (306) als das Ansprechen auf Übertragungen von wenigstens zwei Sendern (304), wobei die Sender (304) und Empfänger (306) auf zwei oder mehr Bohrlöcher verteilt sind;
    Ableiten wenigstens eines kompensierten Werts (Pikjm), der Daten aus dem Ansprechen der wenigstens zwei Empfänger (306) auf die wenigstens zwei Sender (304) kombiniert, als ein Verhältnis, das Fehler zwischen tatsächlichen und angenommenen Charakteristiken der wenigstens zwei Empfänger (306) und der wenigstens zwei Sender (304) aufhebt; und
    Durchführen einer Inversion wenigstens teilweise auf Grundlage des kompensierten Werts, um ein unterirdisches Bild zu erlangen, wobei der kompensierte Wert (Pikjm) Ungenauigkeiten einschließlich Sensorverstärkungen und Phasenvariationen kompensiert, wobei der kompensierte Wert (Pikjm) durch V ij V km V im V kj ,
    Figure imgb0011
    ausgedrückt werden kann, wobei i und k Empfängerindizes sind und j und m Senderindizes sind, v' die erlangten Daten darstellt, und wobei die Inversion ein Vorwärtsmodell benutzt, das Eigenschaften der unterirdischen Formation mit dem wenigstens einen kompensierten Wert (Pikjm) in Beziehung setzt.
  2. Verfahren nach einem der vorhergehenden Ansprüche, wobei das Erlangen von Daten ferner das Erlangen von Daten von weiteren Empfängern (306) als das Ansprechen auf weitere Paare von Sendern (304) beinhaltet.
  3. Verfahren nach einem der vorhergehenden Ansprüche, wobei die Daten als ein komplexer Wert ausgedrückt werden können, um Verstärkung und Phase darzustellen.
  4. Verfahren nach einem der vorhergehenden Ansprüche, wobei die Empfänger (306) wenigstens eins der Folgenden beinhalten: einen magnetischen Dipol, einen elektrischen Dipol, beabstandete Elektroden.
  5. Verfahren nach einem der vorhergehenden Ansprüche, wobei das Bild Informationen zu wenigstens einer Formationseigenschaft bereitstellt, die ausgewählt wird aus einer Gruppe bestehend aus: akustischer Geschwindigkeit, Dämpfung, seismischem Reflexionsgrad und elektromagnetischem Widerstand.
  6. System zum Verbessern der Bildqualität von Crosswell-Tomografie, umfassend:
    eine Speichereinheit (40), die Inversionssoftware speichert; und
    einen Prozessor, der an den Speicher (40) gekoppelt ist, um die Software auszuführen, wobei die Software den Prozessor zu Folgendem konfiguriert:
    Erlangen von Daten von wenigstens zwei Empfängern (306) als das Ansprechen auf Übertragungen von wenigstens zwei Sendern (304), wobei sich einer der Sender (304) oder Empfänger (306) in einem Bohrloch befindet;
    Ableiten wenigstens eines kompensierten Werts (Pikjm) auf Grundlage des Ansprechens der Empfänger (306) auf die einzelnen Sender (304), die als ein Verhältnis kombiniert werden, das Fehler zwischen tatsächlichen und angenommenen Charakteristiken der Empfänger (306) und der Sender (304) aufhebt; und
    Durchführen einer Inversion wenigstens teilweise auf Grundlage des kompensierten Werts (Pikjm), um ein unterirdisches Bild zu erlangen, wobei der kompensierte Wert Ungenauigkeiten einschließlich Sensorverstärkungen und Phasenvariationen kompensiert, wobei
    der kompensierte Wert durch V ij V km V im V kj ,
    Figure imgb0012
    ausgedrückt werden kann, wobei i und k Empfängerindizes sind und j und m Senderindizes sind, v' die erlangten Daten darstellt, wobei die Inversion ein Vorwärtsmodell benutzt, das Eigenschaften der unterirdischen Formation mit dem wenigstens einen kompensierten Wert (Pikjm) in Beziehung setzt und wobei die Sender (304) und Empfänger (306) auf zwei oder mehr Bohrlöcher verteilt sind.
  7. System nach Anspruch 6, ferner umfassend ein Netz, das elektrisch an den Prozessor gekoppelt ist, wenigstens zwei Sender (304) und wenigstens zwei Empfänger (306), wobei wenigstens einer der Sender oder einer der Empfänger innerhalb des Bohrlochs angeordnet ist, und wobei wahlweise das Netz elektrisch an wenigstens einen Sensor, wenigstens einen Sender und wenigstens einen Empfänger gekoppelt ist.
  8. System nach Anspruch 6 oder 7, wobei die Software den Prozessor ferner zum Erlangen von Daten von weiteren Empfängern als das Ansprechen auf weitere Paare von Sendern konfiguriert.
  9. System nach einem der Ansprüche 6 bis 8, wobei die Daten als ein komplexer Wert ausgedrückt werden können, um Verstärkung und Phase darzustellen.
  10. System nach einem der Ansprüche 6 bis 9, wobei das Bild Informationen zu wenigstens einer Formationseigenschaft bereitstellt, die ausgewählt ist aus einer Gruppe bestehend aus: akustischer Geschwindigkeit, Dämpfung, seismischem Reflexionsgrad und elektromagnetischem Widerstand.
  11. System nach einem der Ansprüche 6 bis 10, wobei die Sender und Empfänger jeweils an eine Antenne gekoppelt sind, wobei die Antenne in der Gruppe bestehend aus Spulen, Solenoiden, Drahtantennen, Elektroden, Gap-Antenne und Toroiden ist.
EP11873738.6A 2011-10-06 2011-10-06 Kompensiertes cross-well-tomografieverfahren und systeme Not-in-force EP2742365B1 (de)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2011/055020 WO2013052049A1 (en) 2011-10-06 2011-10-06 Compensated cross-well tomography methods and systems

Publications (3)

Publication Number Publication Date
EP2742365A1 EP2742365A1 (de) 2014-06-18
EP2742365A4 EP2742365A4 (de) 2015-01-14
EP2742365B1 true EP2742365B1 (de) 2017-03-01

Family

ID=48044027

Family Applications (1)

Application Number Title Priority Date Filing Date
EP11873738.6A Not-in-force EP2742365B1 (de) 2011-10-06 2011-10-06 Kompensiertes cross-well-tomografieverfahren und systeme

Country Status (9)

Country Link
US (1) US9008970B2 (de)
EP (1) EP2742365B1 (de)
CN (1) CN103842852B (de)
AU (1) AU2011378454B2 (de)
BR (1) BR112014005881A2 (de)
CA (1) CA2845627C (de)
IN (1) IN2014DN01793A (de)
RU (1) RU2577418C2 (de)
WO (1) WO2013052049A1 (de)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2845627C (en) 2011-10-06 2016-11-29 Halliburton Energy Services, Inc. Compensated cross-well tomography methods and systems
WO2015050954A1 (en) * 2013-10-03 2015-04-09 Halliburton Energy Services, Inc. Downhole tool with radial array of conformable sensors for downhole detection and imaging
WO2015088563A1 (en) 2013-12-13 2015-06-18 Halliburton Energy Services, Inc. Methods and systems of electromagnetic interferometry for downhole environments
CA2947143C (en) 2014-05-01 2020-03-24 Halliburton Energy Services, Inc. Casing segment having at least one transmission crossover arrangement
EP3134607B1 (de) 2014-05-01 2019-03-13 Halliburton Energy Services, Inc. Multilaterale produktionssteuerungsverfahren und systeme unter einsatz eines gehäusesegments mit mindestens einer übertragungsübergangsanordnung
US10145233B2 (en) 2014-05-01 2018-12-04 Halliburton Energy Services, Inc. Guided drilling methods and systems employing a casing segment with at least one transmission crossover arrangement
BR112016025406B1 (pt) * 2014-05-01 2022-11-01 Halliburton Energy Services, Inc Método de tomografia interpoço e sistema de tomografia interpoço
CN104391333B (zh) * 2014-10-21 2017-04-26 安徽理工大学 多井间地质信息的探测处理系统与方法
US10302796B2 (en) 2014-11-26 2019-05-28 Halliburton Energy Services, Inc. Onshore electromagnetic reservoir monitoring
CN104360396B (zh) * 2014-12-04 2016-10-19 中国海洋石油总公司 一种海上井间tti介质三种初至波走时层析成像方法
WO2017095418A1 (en) * 2015-12-03 2017-06-08 Halliburton Energy Services, Inc. Crosswell tomography using an array of optical fiber transducers
US10502044B2 (en) 2016-08-12 2019-12-10 Halliburton Energy Services, Inc. Multistage processing and inversion of corrosion detection tools
CN110192004B (zh) * 2017-02-06 2023-11-10 哈利伯顿能源服务公司 到矿床边界的距离反演解的像素化
US11723579B2 (en) 2017-09-19 2023-08-15 Neuroenhancement Lab, LLC Method and apparatus for neuroenhancement
US11717686B2 (en) 2017-12-04 2023-08-08 Neuroenhancement Lab, LLC Method and apparatus for neuroenhancement to facilitate learning and performance
EP3731749A4 (de) 2017-12-31 2022-07-27 Neuroenhancement Lab, LLC System und verfahren zur neuroverstärkung zur verbesserung der emotionalen reaktion
US11364361B2 (en) 2018-04-20 2022-06-21 Neuroenhancement Lab, LLC System and method for inducing sleep by transplanting mental states
US11452839B2 (en) 2018-09-14 2022-09-27 Neuroenhancement Lab, LLC System and method of improving sleep
RU2706205C1 (ru) * 2019-04-30 2019-11-14 Общество С Ограниченной Ответственностью "Радионда" Система радиоволновой геоинтроскопии межскважинного пространства
US11817799B2 (en) * 2020-11-20 2023-11-14 Halliburton Energy Services, Inc. Single crystal ultrasonic transducer with charge mode receiver

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2628027A1 (de) * 2010-11-12 2013-08-21 Halliburton Energy Services, Inc. System und verfahren für umweltmessungen

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4446434A (en) 1978-12-20 1984-05-01 Conoco Inc. Hydrocarbon prospecting method with changing of electrode spacing for the indirect detection of hydrocarbon reservoirs
US5355088A (en) * 1991-04-16 1994-10-11 Schlumberger Technology Corporation Method and apparatus for determining parameters of a transition zone of a formation traversed by a wellbore and generating a more accurate output record medium
BR9805949A (pt) * 1997-03-17 1999-08-31 Yamamoto Engineering Corp Dispositivo e método para transmissão de onda acústica no subsolo, para recebimento de onda acústica, e método de exploracão de subsolo usando o dispositivo acima mencionado.
US6703838B2 (en) 1998-04-13 2004-03-09 Schlumberger Technology Corporation Method and apparatus for measuring characteristics of geological formations
US6534986B2 (en) 2000-05-01 2003-03-18 Schlumberger Technology Corporation Permanently emplaced electromagnetic system and method for measuring formation resistivity adjacent to and between wells
US6393363B1 (en) * 2000-06-28 2002-05-21 Schlumberger Technology Corporation Method and apparatus for the measurement of the electrical resistivity of geologic formations employing modeling data
US6538447B2 (en) 2000-12-13 2003-03-25 Halliburton Energy Services, Inc. Compensated multi-mode elctromagnetic wave resistivity tool
US8296113B2 (en) 2001-05-18 2012-10-23 Halliburton Energy Services, Inc. Virtual steering of induction tool attenuation and phase difference measurements
US6777940B2 (en) * 2002-11-08 2004-08-17 Ultima Labs, Inc. Apparatus and method for resistivity well logging
US7046581B2 (en) 2003-12-01 2006-05-16 Shell Oil Company Well-to-well tomography
WO2007053251A1 (en) 2005-11-01 2007-05-10 Exxonmobil Upstream Research Company Method for phase and amplitude correction in controlled source electromagnetic survey data
KR100861084B1 (ko) * 2006-12-08 2008-09-30 한국지질자원연구원 물리탐사 자료의 4차원 역산 방법 및 이를 이용한지하구조의 4차원 영상화 방법.
US8131526B2 (en) * 2007-04-14 2012-03-06 Schlumberger Technology Corporation System and method for evaluating petroleum reservoir using forward modeling
US7565244B2 (en) * 2007-06-27 2009-07-21 Schlumberger Technology Corporation Method and system for removing effects of conductive casings and wellbore and surface heterogeneity in electromagnetic imaging surveys
US8129993B2 (en) * 2007-07-10 2012-03-06 Schlumberger Technology Corporation Determining formation parameters using electromagnetic coupling components
US7852087B2 (en) 2007-08-10 2010-12-14 Schlumberger Technology Corporation Removing effects of near surface geology from surface-to-borehole electromagnetic data
US8401796B2 (en) 2008-09-29 2013-03-19 Schlumberger Technology Corporation Methods and systems for acoustically monitoring formations
US8456166B2 (en) * 2008-12-02 2013-06-04 Schlumberger Technology Corporation Single-well through casing induction logging tool
US8812237B2 (en) 2009-02-05 2014-08-19 Schlumberger Technology Corporation Deep-reading electromagnetic data acquisition method
WO2012005737A1 (en) 2010-07-09 2012-01-12 Halliburton Energy Services, Inc. Imaging and sensing of subterranean reservoirs
CA2845627C (en) 2011-10-06 2016-11-29 Halliburton Energy Services, Inc. Compensated cross-well tomography methods and systems

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2628027A1 (de) * 2010-11-12 2013-08-21 Halliburton Energy Services, Inc. System und verfahren für umweltmessungen

Also Published As

Publication number Publication date
RU2014113422A (ru) 2015-11-20
BR112014005881A2 (pt) 2017-03-28
US9008970B2 (en) 2015-04-14
IN2014DN01793A (de) 2015-05-15
CA2845627C (en) 2016-11-29
AU2011378454B2 (en) 2015-05-21
EP2742365A1 (de) 2014-06-18
CN103842852B (zh) 2017-11-10
US20140244175A1 (en) 2014-08-28
CA2845627A1 (en) 2013-04-11
RU2577418C2 (ru) 2016-03-20
AU2011378454A1 (en) 2014-03-06
WO2013052049A1 (en) 2013-04-11
CN103842852A (zh) 2014-06-04
EP2742365A4 (de) 2015-01-14

Similar Documents

Publication Publication Date Title
EP2742365B1 (de) Kompensiertes cross-well-tomografieverfahren und systeme
US10627541B2 (en) System, method and computer-program product for in-situ calibration of a wellbore resistivity logging tool
AU2003278866B2 (en) Fixed-depth of investigation log for multi-spacing multi-frequency LWD resistivity tools
US10451765B2 (en) Post-well reservoir characterization using image-constrained inversion
US20180203151A1 (en) Measuring petrophysical properties of an earth formation by regularized direct inversion of electromagnetic signals
US10359535B2 (en) Electrode-based tool measurement corrections based on measured leakage currents
US10031254B2 (en) Electrode-based tool measurement corrections based on leakage currents estimated using a predetermined internal impedance model or table
AU2017263252B2 (en) Methods and systems employing look-around and look-ahead inversion of downhole measurements
US10914859B2 (en) Real-time true resistivity estimation for logging-while-drilling tools
NO20160142A1 (en) Surface calibration of a wellbore resistivity logging tool
US10598816B2 (en) Systems and methods for relative dip correction
US10302802B2 (en) Apparatus and methods of skin effect correction
WO2017078740A1 (en) Downhole logging systems and methods employing adjustably-spaced modules
US10459110B2 (en) Flexible conductive shield for downhole electromagnetic noise suppression
AU2016433067A1 (en) Use of gap subs behind a coil antenna in electromagnetic induction tools

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

17P Request for examination filed

Effective date: 20140313

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

RIN1 Information on inventor provided before grant (corrected)

Inventor name: SAN MARTIN, LUIS, E.

Inventor name: BITTAR, MICHAEL, S.

Inventor name: GUNER, BARIS

Inventor name: DONDERICI, BURKAY

A4 Supplementary search report drawn up and despatched

Effective date: 20141215

DAX Request for extension of the european patent (deleted)
RIC1 Information provided on ipc code assigned before grant

Ipc: G01V 1/42 20060101ALN20141209BHEP

Ipc: G01V 3/18 20060101AFI20141209BHEP

17Q First examination report despatched

Effective date: 20160205

REG Reference to a national code

Ref country code: DE

Ref legal event code: R079

Ref document number: 602011035585

Country of ref document: DE

Free format text: PREVIOUS MAIN CLASS: G01V0001000000

Ipc: G01V0003180000

RIC1 Information provided on ipc code assigned before grant

Ipc: G01V 1/42 20060101ALN20160810BHEP

Ipc: G01V 3/18 20060101AFI20160810BHEP

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

RIC1 Information provided on ipc code assigned before grant

Ipc: G01V 1/42 20060101ALN20160907BHEP

Ipc: G01V 3/18 20060101AFI20160907BHEP

RIC1 Information provided on ipc code assigned before grant

Ipc: G01V 3/18 20060101AFI20160912BHEP

Ipc: G01V 1/42 20060101ALN20160912BHEP

INTG Intention to grant announced

Effective date: 20160928

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: AT

Ref legal event code: REF

Ref document number: 872001

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170315

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602011035585

Country of ref document: DE

REG Reference to a national code

Ref country code: NO

Ref legal event code: T2

Effective date: 20170301

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20170301

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 872001

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170301

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

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170602

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

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

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170601

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

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 FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

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 FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

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

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170701

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170703

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602011035585

Country of ref document: DE

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

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

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

26N No opposition filed

Effective date: 20171204

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

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 602011035585

Country of ref document: DE

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

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: IE

Ref legal event code: MM4A

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20180629

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

Ref country code: LU

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

Effective date: 20171006

Ref country code: DE

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

Effective date: 20180501

Ref country code: LI

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

Effective date: 20171031

Ref country code: CH

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

Effective date: 20171031

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20171031

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

Ref country code: FR

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

Effective date: 20171031

Ref country code: BE

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

Effective date: 20171031

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

Ref country code: MT

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

Effective date: 20171006

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

Ref country code: IE

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

Effective date: 20171006

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

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20111006

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

Ref country code: CY

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

Effective date: 20170301

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

Ref country code: NO

Payment date: 20190925

Year of fee payment: 9

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

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

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

Ref country code: GB

Payment date: 20190823

Year of fee payment: 9

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

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

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

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170301

REG Reference to a national code

Ref country code: NO

Ref legal event code: MMEP

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

Effective date: 20201006

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

Ref country code: NO

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

Effective date: 20201031

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: 20201006