NO336562B1 - Analyse av data fra målinger med nukleær magnetisk resonans - Google Patents
Analyse av data fra målinger med nukleær magnetisk resonansInfo
- Publication number
- NO336562B1 NO336562B1 NO20030223A NO20030223A NO336562B1 NO 336562 B1 NO336562 B1 NO 336562B1 NO 20030223 A NO20030223 A NO 20030223A NO 20030223 A NO20030223 A NO 20030223A NO 336562 B1 NO336562 B1 NO 336562B1
- Authority
- NO
- Norway
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- spins
- sequence
- data
- Prior art date
Links
- 238000005481 NMR spectroscopy Methods 0.000 title claims description 39
- 238000005259 measurement Methods 0.000 title claims description 16
- 238000004458 analytical method Methods 0.000 title description 3
- 230000005291 magnetic effect Effects 0.000 claims description 68
- 238000000034 method Methods 0.000 claims description 33
- 238000002592 echocardiography Methods 0.000 claims description 22
- 239000012530 fluid Substances 0.000 claims description 17
- 239000011435 rock Substances 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 7
- 238000007906 compression Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 239000002023 wood Substances 0.000 claims 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 25
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 23
- 230000005415 magnetization Effects 0.000 description 20
- 230000006870 function Effects 0.000 description 15
- 238000002075 inversion recovery Methods 0.000 description 15
- 230000003068 static effect Effects 0.000 description 13
- 229920006395 saturated elastomer Polymers 0.000 description 10
- 238000011084 recovery Methods 0.000 description 9
- 239000011159 matrix material Substances 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000001351 cycling effect Effects 0.000 description 5
- 239000011148 porous material Substances 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000005755 formation reaction Methods 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 150000001875 compounds Chemical class 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000004181 pedogenesis Methods 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 238000011005 laboratory method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 229910052729 chemical element Inorganic materials 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 150000002170 ethers Chemical class 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000005298 paramagnetic effect Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000000293 three-dimensional nuclear magnetic resonance spectroscopy Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/46—NMR spectroscopy
- G01R33/4625—Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Geochemistry & Mineralogy (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Signal Processing (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Remote Sensing (AREA)
- Geophysics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Description
Denne patentsøknaden krever prioritet fra US foreløpig patentsøknad nr.
60/220 053, innlevert 21. juli 2000, og fra US patentsøknad nr. 09/723 803 innlevert 28. november 2000, som begge ved denne referanse herved inkorporeres i sin helhet.
Denne oppfinnelsen vedrører nukleær magnetisk resonans (NMR)-målinger, og mer bestemt, analyser av NMR-data.
NMR har vært en vanlig laboratorieteknikk i over førti år, og har blitt et viktig verktøy ved formasjonsevaluering. Generell bakgrunn for NMR-brønnlogging kan finnes feks. i US patent nr. 5 023 551 tilhørende Kleinberg et al., som er overdratt til den samme fullmektigen som den foreliggende oppfinnelse, og som ved denne referanse herved inkorporeres i sin helhet.
NMR bygger på det faktum at kjernene i mange kjemiske elementer har impulsmoment («spin») og et magnetisk moment. I et utvendig påført statisk magnetfelt, innretter kjernenes spin seg langs retningen til det statiske felt. Denne likevektssituasjonen kan forstyrres av en puls av et oscillerende magnetfelt (eksempelvis en RF-puls) som vipper spinnene bort fra det statiske feltets retning. Vinkelen som spinnene vippes er gitt av 9=yBitp/2, hvor y er det gyromagnetiske forhold, Bi er den lineære polariserte oscillerende feltstyrke, og tp er pulsens varighet. Vipping av pulser 90 og 180° er det mest vanlige.
Etter vipping skjer to ting samtidig. For det første presiserer spinnene rundt retningen av det statiske felt ved Larmor-frekvensen, gitt ved (o0=yB0, hvor B0er styrken av det statiske feltet og y er det gyromagnetiske forholdet. For hydrogenkjerner er y/27t=4258 Hz/Gauss, slik at, feks., i et statisk felt på 235 Gauss, vil hydrogenspinnene presisere ved en frekvens på 1 MHz. For det andre returnerer spinnene til likevektsretningen i henhold til en svekkelsestid, Ti, som er kjent som spinn-gitter relaksasjonstiden. Fordi spinn-gitter relaksasjon skjer langs likevektsretningen, benevnes Ti også den longitudinale relaksasjonstidskonstant.
De molekylære kjerners spinn er også forbundet med en annen relaksasjonstid, T2, benevnt spinn-spinnrelaksasjonstiden. Ved enden av 90° vippepulsen, peker alle spinnene i en felles retning perpendikulært, eller transversalt på det statiske felt, og alle presiserer ved Larmor-frekvensen. På grunn av små fluktuasjoner ved statiske felt forårsaket av andre spinn eller paramagnetiske urenheter, presiserer imidlertid spinnene med litt forskjellige frekvenser, og den transversale magnetisering avfases med en tidskonstant T2, som også benevnes den transversale
relaksasj onstidskonstant.
En standard teknikk til måling av T2, både i laboratoriet og ved brønnlogging, bruker en RF-pulssekvens som er kjent som CPMG (Carr-Purcell-Meiboom-Gill)-sekvens. Som velkjent, etter en ventetid som kommer foran hver pulssekvens, vipper en 90° puls spinnene inn i det transversale plan og bevirker at spinnene begynner og presessere. Deretter påføres en 180° puls som holder spinnene i måleplanet, men som forårsaker at spinnene, som avfases i det transversale plan, reverserer sin retning og refokuserer. Ved gjentatt reversering av spinnene ved bruk av en serie av 180° pulser fremkommer en serie av «spinnekkoer». Toget av ekkoer måles og prosesseres for å bestemme den irreversible avfasingstidskonstant T2. Ved brønnloggeanvendelser har de detekterte spinnekkoene blitt brukt til å ekstrahere
oljefeltparametere så som porøsitet, porestørrelsesfordeling, og oljeviskositet.
I teorien kan andre NMR-målinger utført i laboratoriet anvendes ved brønnlogging for å ekstrahere ytterligere informasjon om oljefeltet, men i praksis gjør karakteren av brønnlogging og borehullets omgivelse implementering av enkelte laboratorie-NMR-målinger vanskelige. F.eks. har inversjonsgjenvinning en vanlig laboratorieteknikk for måling av Ti. Ved en inversjonsgjenvinningsmåling påføres en 180° puls på et system av spinn som er innrettet langs det statiske magnetfelt for å reversere retningen av spinnene. Systemet av spinn som perturberes på denne måten begynner å svekkes mot sin likevektsretning i henhold til Ti. For å måle den netto magnetiseringen, påføres en 90° puls for å rotere spinnene inn i det transversale plan og således fremkalle et målbart signal. Signalet vil begynne å svekkes når spinnene avfases i det transversale plan, men den initiale amplitude av signalet avhenger av «gjenvinningstiden» mellom 180° pulsen og 90° pulsen. Ved å gjenta dette eksperimentet for forskjellige gjenvinningstider og plotting av den initiale amplitude av signalet mot gjenvinningstiden, kan Ti bestemmes. Selv om denne teknikken har blitt brukt med hell i laboratoriet i flere år, er inversjonsgjenvinning meget tidkrevende, og personer med ordinær fagkunnskap innen området forstår at inversjonsgjenvinning kan være uegnet ved brønnloggeanvendelser.
Det fortsetter følgelig å være et generelt behov for forbedrede NMR-målinger, og, særlig innen olje- og gassutvinningsindustrien, forbedrede NMR-metoder som kan brukes for å ekstrahere informasjon om bergartprøver og brukes ved brønnloggeanvendelser.
Oppfinnelsen tilveiebringer en fremgangsmåte til ekstrahering av informasjon om et system av nukleære spinn, så som i et fluid som kan befinne seg i en bergart eller innenfor et parti av en jordartformasjon som omgir et borehull, som heretter brukt inkluderer uttrykket «bergart» jordarter, jordartformasjoner og et parti av en jordartformasjon), eller andre porøse omgivelser. Fremgangsmåten involverer utførelse av minst én nukleær magnetisk resonansmåling på et system av nukleære spinn og innsamling av nukleære magnetiske resonansdata fra hver av målingene. De nukleære magnetiske resonansdata uttrykkes ved bruk av en kjerne som er separerbar langs minst to dimensjoner, komprimeres langs hver dimensjon av kjernen, og analyseres deretter for å ekstrahere informasjon om systemene av spinn. Foreliggende oppfinnelse angår en metode og en logganordning for ekstrahering av informasjon om et system av nukleære spinn i samsvar med patentkravene log 8.
Ytterligere detaljer og trekk ved oppfinnelsen vil fremgå klarere av den følgende detaljerte beskrivelse.
Oppfinnelsen vil nedenfor bli beskrevet i nærmere detalj sammen med de følgende figurer, hvor: Fig. 1 skjematisk representerer en magnetisk feltpulssekvens som kan brukes i samsvar med oppfinnelsen; Fig. 2 illustrerer en inversjonsgjenvinning-CPMG-sekvens; Fig. 3 illustrerer en utførelse av en drevet likevekts-refokuseringssekvens ifølge oppfinnelsen; Fig. 4A og 4B viser utførelser av en diffusjonsredigerings-CPMG-sekvens ifølge oppfinnelsen; Fig. 5 viser skjematisk en utførelse av en nukleær magnetisk resonansmåling ifølge oppfinnelsen; Fig. 6 viser skjematisk en alternativ utførelse av en nukleær magnetisk resonansmåling ifølge oppfinnelsen; Fig. 7 inneholder en todimensjonal avbildning av Ti mot T2for en vannmettet bergartprøve generert i samsvar med oppfinnelsen; Fig. 8 inneholder en todimensjonal avbildning av Ti mot T2for en oljemettet bergartprøve generert i samsvar med oppfinnelsen; Fig. 9 inneholder en todimensjonal avbildning av <Ti/T2> mot T2for en bergartprøve generert i samsvar med oppfinnelsen; og Fig. 10 illustrerer en rutine som kan brukes ved implementering av en utførelse av en fremgangsmåte ifølge oppfinnelsen. Fig. 1 viser skjematisk en sekvens av magnetfeltpulser, så som RF-pulser, som kan påføres på et system av nukleære spinn, så som i et fluid i en bergart, i samsvar med oppfinnelsen. Magnetfeltpulssekvensen 100 inkluderer to deler. En første del 110 frembringer et system av nukleære spinn i en første tilstand. En annen del 120 genererer fra systemet av spinn en serie av magnetiske resonanssignaler. Som vist på fig. 1 kan magnetfeltpulssekvensen 100 påføres gjentatte ganger, og magnetresonanssignalene detekteres og kombineres for oppbygging av signal/støy-forhold.
Ved generering av seriene av magnetiske resonanssignaler eksiterer den første og annen del fortrinnsvis spinn ved tilnærmet det samme frekvensområdet. F.eks. har en standard CPMG-sekvens med en serie på 180° pulser en tilbøyelighet til å refokusere spinn som befinner seg på resonansfrekvensen, nemlig de spinn som har en frekvens som er hovedsakelig lik Larmor-frekvensen, og det observerte CPMG-signal inkluderer primært slike spinn som befinner seg på resonansfrekvensen. En standard CPMG-sekvens vil følgelig fortrinnsvis pares med en første del som også primært eksiterer spinn på resonansfrekvensen. På en første del som primært påvirker spinn som befinner seg utenfor resonansfrekvensen, så som en konstant RF-irradioasjon (omtalt nedenfor) er den annen del fortrinnsvis konstruert til å eksitere hovedsakelig de samme spinnene som befinner seg utenfor resonansfrekvensen. En CPMG-lignende sekvens, hvor sammensatte pulser som er konstruert til å refokusere trinnene utenfor resonansfrekvensen erstatter 180° pulsen, kan brukes i slike magnetiske pulssekvenser. Fagpersoner innen området vil være i stand til å konstruere andre typer refokuseringssekvenser for å refokusere spinn som har forskjellige frekvenser utenfor resonansfrekvensen.
F.eks. viser fig. 2 en magnetfeltpulssekvens 200 som inkluderer en inversjonsgjenvinningssekvens 210 som er fulgt av en CPMG-sekvens 220. Inversjonsgjenvinningssekvensen inkluderer en 180° puls 212 som reverserer retningen av systemet av spinn, fulgt av en gjenvinningstid, Ti. Etter gjenvinningstiden, som typisk er mindre enn Ti, påføres CPMG-sekvensen 220, med en 90° puls 222 som roterer spinnene inn i det transversale plan, og en serie av 180° pulser (hvorav to er vist med 224) som genererer en serie av spinnekkoer (hvorav to er vist med 226) viss amplitude avhenger av Ti og viss svekking avhenger av den transversale relaksasjonstidskonstant, T2:
hvor Ti, T2og Ti er som definert ovenfor, Mo er en netto magnetisering av spinnene ved termisk likevekt, t2er tiden fra start av CPMG-sekvensen, og f(Ti, T2) er den todimensjonale tetthetsfunksjon for Ti og T2.
Et annet eksempel på en magnetisk feltpulssekvens som kan brukes i samsvar med oppfinnelsen inkluderer en første del som er konstruert til å frembringe systemet av spinn i en drevet likevektstilstand. Generelt vil enhver magnetfeltpulssekvens som gjentatte ganger roterer systemets netto magnetisering mellom den longitudinale og transversale retning frembringe en betydelig drevet likevektsmagnetisering. Så lenge de gjentagende magnetfeltpulsenheter er korte sammenlignet med Ti og T2, kan det vises at den resulterende likevektsmagnetisering avhenger av en funksjon av forholdet, T1/T2:
Etter at spinnene er frembragt i en drevet likevekt, påføres en annen del som er konstruert til å refokusere systmets spinn. Den annen del genererer en serie av magnetiske resonanssignaler som avhenger både av Ti/T2-forholdet og T2:
Ti, T2ogT2er som tidligere definert, og /(T2, — T) er den todimensjonale
^2
tetthetsfunksjon for forholdet T1/T2og T2.
Fig. 3 viser et eksempel på en slik drevet likevekts-refokuseringssekvens. Den
magnetiske feltpulssekvens 300 inkluderer en drevet likevekt fouriertransformasjon (DEFT) sekvens 310 fulgt av en CPMG-sekvens 320. Den viste DEFT-sekvens 310 kan betraktes som en repeterende blokk 315 av magnetfeltpulser generert i henhold til:
hvor ti og t2(«2ti) er tidsavstander mellom magnetiske pulser. Hver repeterende blokk er atskilt med en tid, 81, og hver blokk har en total bredde på 82. Fig. 3 viser sekvensen med to repeterende blokker (i=2); fagpersoner innen området vil imidlertid forstå at sekvensen kan ha ethvert antall repeterende blokker. Ved mange anvendelser kan det være ønskelig å påføre DEFT-sekvensen over en tidsperiode som er lengre enn Ti. Et magnetisk resonansspinnekko 317 fremtrer midtveis mellom de to 180y-graders pulser i hver blokk. Over tid frembringer DEFT-sekvensen systemet av spinn i en drevet likevekt hvor magnetiseringen, innenn grensen 8i«Tiog 82«T2, er gitt som:
Det asymptotiske uttrykk i ligning (5) gir en god approksimasjon for den DEFT-drevne likevektsmagnetisering for alle parametere i området tilnærmet T2/82> 10. Etter at systemet av spinn er frembrakt i en drevet likevekt, påføres CPMG-sekvensen 320, med en 90° puls 322 fulgt av en serie av 180° pulser (hvorav to er vist med 324) for å generere en serie av magnetiske resonansspinnekkoer (hvorav to er vist med 326) hvis signal kan skrives som ligning (3) ovenfor, hvor funksjonen g — er gitt av:
som er mest følsom for endringer i T1/T2i området hvor (5i/52) » <Ti/T2.
Et annet eksempel på en drevet likevekts-refokuseringsmagnetfeltpulssekvens bruker et konstant RF-felt, som typisk påføres i en tid som er større enn Ti, for å frembringe systemet av spinn i en drevet likevekt som avhenger av Ti/T2-forholdet. Det konstante RF-felt har typisk en tilstrekkelig høy frekvens,©1, slik at coiT2»l, i hvilket tilfelle likevektsmagnetiseringen langs y- henholdsvis z-aksene er gitt av:
Ti, T2, (Oi og Mo er som tidligere definert, og Ac&oer forskjellen mellom frekvensen av det påførte RF-felt og spinnenes Larmor-frekvens. Fra ligning (7) og (8) kan det ses at på resonansen, dvs. der hvor frekvensen av det påførte RF-felt er lik Larmor-frekvensen (Aco0=0) blir likevektsmagnetiseringen borte. Ved Aa>0=±a>iblir
M
likevektsmagnetiseringen langs både y- og z-aksen - , som er ekvivalent til
\ + (TJT2)
DEFT-likevektsmagnetiseringen for 81= 82(se ligning (5) ovenfor). For å måle T2svekkelsen av denne likevektsmagnetiseringen følger en CPMG-lignende sekvens som består av serier av sammensatte pulser som er konstruert til å refokusere spinn utenfor resonansfrekvensen fortrinnsvis den konstante RF-irradiasjon. F.eks. kan spinn som har en magnetisering rundt Aa>0=±a>ieffektivt refokuseres ved bruk av en
CPMG-lignende serie av sammensatte pulser av 127x-127.x grader, hvilket
genererer en serie av spinnekkoer ifølge ligning (3) over, hvor funksjonen g(Ti/Ti)
er gitt av
En annen type magnetfeltpulssekvens som kan brukes i samsvar med oppfinnelsen inkluderer en første del som frembringer et system av spinn med en netto magnetisering som avhenger av en diffusjonskoeffisient. En CPMG-sekvens følger typisk den første del og genererer en serie av magnetiske resonansspinnekkoer som generelt svekkes i henhold til:
hvor D er diffusjonskoeffisienten, Ti er en tid som er forbundet med den første del, y er det gyromagnetiske forhold, g er gradientens feltstyrke, og t2er tiden fra starten av CPMG-sekvensen. F.eks., i henhold til en utførelse vist på fig. 4A, inkluderer den første del 410 to spinnekkoer som er atskilt av en forlenget ekkoavstand, tE, fulgt av en standard CPMG-sekvens på 180° pulser 420, og ti<3>= t£3. Ifølge en annen utførelse vist på fig. 4B, inkluderer den første del 430 en stimulert ekkosekvens som erkarakterisert veden tidsavstand 8 og et stimulert ekko som opptrer ved tiden td, idet den annen del 440 inkluderer en standard CPMG-sekvens av 180° pulser, ogTi<3>= 6td8<2>. Diffusjonsredigering av magnetfeltpulssekvenser er beskrevet i nærmere detalj i US patentsøknad nr. 09/723 803, som herved inkorporeres ved referanse.
Fig. 5 illustrerer skjematisk en utførelse av en nukleær magnetisk resonansmåling ifølge oppfinnelsen, hvilken involverer frembringelse av et system av spinn i fluidet med forskjellige initiale tilstander, og deretter måling av svekkelsen av spinnene fra de forskjellige initiale tilstander ved bruk av den samme refokuseringssekvens. En første sekvens 500 som inkluderer en første del 510 fulgt av en annen del 520 påføres på fluidet for å generere en første serie av magnetiske resonansspinnekkoer. En annen sekvens 550 som inkluderer en tredje part 530 fulgt av den annen del 520 påføres for å generere en annen serie av magnetiske resonansspinnekkoer. Den første og annen serie av magnetiske resonansspinnekkoer detekteres og analyseres for å ekstrahere informasjon om fluidet i bergarten. Ytterligere sekvenser (ikke vist) kan påføres på systemet av spinn for å bringe systemet i en flerhet av forskjellige initiale tilstander hvorfra ytterligere serier av magnetiske resonansspinnekkoer kan genereres, detekteres og analyseres.
Som vist på fig. 5 kan den første sekvens 500 påføres gjentatte ganger før den annen sekvens 550 påføres på fluidet. Den annen sekvens 550, kan som vist på fig. 5, også påføres gjentatte ganger. Magnetfeltpulsene i hver sekvens kan være fasesyklisert for å bidra til å redusere støy på grunn av feiltilpassede pulser. Personer med ordinær kunnskap innen fagområdet vil forstå at fasesykli sering kan utføres ved å variere fasen til magnetfeltpulsene innenfor en gitt sekvens, eller ved å variere fasen til magnetfeltpulsene fra sekvens til sekvens.
Den første sekvens 500 omfatter magnetfeltpulssekvensene beskrevet ovenfor, og inkluderer således sekvenser som har en første del 510 som frembringer spinnene initialt i en første tilstand med en magnetisering som avhenger av en første parameter, så som Ti, T1/T2, D, osv. Den annen sekvens 550 frembringer spinnene initialt i en annen tilstand som er forskjellig fra den første tilstand, og kan inkludere de ovenfor beskrevne magnetfeltpulssekvenser for å påføre spinnene en magnetisering som har en forskjellig avhengighet på den første parameter. Alternativt kan den tredje del 530 inkludere en ventetid som er lengre enn en longitudinal relaksasjonstidskonstant, Ti, under hvilken det ikke påføres noen magnetfeltpulser på fluidet. Ventetiden tillater at spinnene når en termisk likevekt med det statiske magnetfelt, Bo.
Den samme refokuseringssekvens 520 brukes i både den første og den annen sekvens for å generere magnetiske resonanssignaler fra de forskjellige initiale tilstander. Der hvor den tredje del inkluderer en ventetid og den annen del inkluderer en CPMG eller en CPMG-lignende sekvens, kan det magnetiske resonanssignal som er generert av den annen sekvens generelt uttrykkes som et standard CPMG-uttrykk, som avhenger av én parameter, T2:
Den annen del 520 blir fortrinnsvis, men ikke nødvendigvis, påført i hovedsakelig den samme varighet i både den første og den annen sekvens for å generere hovedsakelig det samme antallet spinnekkoer i hver av den første og annen serie.
Som vist på fig. 5 kan den første sekvens 500 påføres gjentatte ganger, og deretter kan den annen sekvens 550 gjentatte ganger påføres på fluidet. Alternativt, som vist på fig. 6, kan den første og annen sekvens flettes inn i hverandre for å muliggjøre innsamling av de to magnetiske resonanssignaler under én enkelt gjennomkjøring. Den nukleære magnetiske resonansmåling vist på fig. 6 representerer en alternativ implementering for frembringelse av et system av spinn i forskjellige tilstander, og deretter måling av svekkelsen av spinnene fra de forskjellige tilstander ved bruk av den samme refokuseringssekvens, som beskrevet ovenfor med henvisning til fig. 5. En første magnetfeltpulssekvens 600 inkluderer en første del 610 for å frembringe et system av spinn i en første tilstand, fulgt av en annen del 620 for å generere en første serie av magnetiske spinnekkoer, og følges av en annen magnetisk feltpulssekvens 650 som inkluderer en tredje del 630 for å bringe systemet i en annen tilstand, fulgt av den annen del 620, for å generere en annen serie av magnetiske spinnekkoer. Den første sekvens 600 og den annen sekvens 650 kan alterneres, som vist, og påføres gjentatte ganger, med passende fasesyklisering, for å bygge opp signal/støy-forhold for den første og annen serie av magnetiske spinnekkoer. Ytterligere sekvenser kan flettes inn i den første og annen sekvens for å bringe systemet av spinn i en flerhet av forskjellige initiale tilstander hvorfra ytterligere serier av magnetiske resonansspinnekkoer kan genereres, detekteres og analyseres.
De magnetiske resonansspinnekkoer generert fra magnetfeltpulssekvenser 500, 550 eller 600, 650 detekteres og analyseres for å ekstrahere informasjon om fluidet. F.eks. kan den første sekvens være en inversjonsgjenvinning-CPMG-sekvens som gir spinnene en ikke-likevektmagnetisering som avhenger av Ti og genererer en første serie av magnetiske resonansspinnekkoer som avhenger av Ti og T2(se ligning (1) ovenfor), og den annen sekvens kan være en inversjonsgjenvinning-CPMG-sekvens med en annen gjenvinningstid som forsyner spinnene med en annen ikke-likevektmagnetisering. En slik inversjonsgjenvinning-CPMG-måling kan gjentas et flertall av ganger, med forskjellige gjenvinningstider, og flertallet av serier av magnetiske resonansspinnekkoer kan analyseres for å ekstrahere en todimensjonal Ti-T2fordeling, f(Ti, T2).
Fig. 7 og 8 inneholder en todimensjonal avbildning av Ti og T2for en bergartprøve av en vannmettet henholdsvis en oljemettet sandstein, generert fra et slikt flertall av inversjonsgjenvinning-CPMG-eksperimenter hvor gjenvinningstiden ble variert mellom ca. 100 us og 10 s. Intensiteten av signalet ved hvert punkt på avbildningen tilveiebringer en indikasjon på antallet spinn, og således fluidets volum, hvilket bidrar til signalet ved dette punktet. Ti har muligens nyttig informasjon. I et fluid i bulkform er Ti typisk lik T2, eller Ti/T2= 1. En porøs omgivelse, så som det som finnes i en bergart eller en jordartformasjon, kan introdusere ytterligere relaksasjonseffekter som påvirker de longitudinale og transversale relaksasjonstider forskjellig. F.eks., i en fuktig fase, så som i en vannmettet bergart, relakserer spinnene overveiende på grunn av kollisjoner med kornoverflater, og den longitudinale og transversale relaksasjonstid påvirkes forskjellig, typisk slik at Ti>T2, hvilket fremgår klart av fig. 7 (linjen 701 indikerer Ti=T2). I en ikke-fuktig fase, så som en oljefase, er imidlertid overflaterelaksasjon neglisjerbar, og Ti forblir hovedsakelig lik T2, dvs. at Ti/T2er tilnærmet lik 1. Fig. 8 viser at for lengre T2, som sampler de store oljefylte porene, er T2tilnærmet lik Ti (linjen 801 indikerer at Ti=T2). Ved kortere T2, som sampler de mindre prøvene, vil imidlertid et lag av vann høyst sannsynlig effektivt fylle porene og introdusere
overflaterelaksasjonseffekter som gjør T2kortere i forhold til Ti. Ti inneholder
således informasjon som kan gjøre det mulig å skjelne fuktige og ikke-fuktige faser innenfor en bergartformasjon, eksempelvis vann mot olje.
Enkelte utførelser av oppfinnelsen kan generere et stort antall datapunkter. F.eks. kan flertallet av inversjonsgjenvinning-CPMG-eksperimenter beskrevet ovenfor generere en million eller flere datapunkter. For å analysere NMR-data som kan genereres ved hjelp av fremgangsmåtene ifølge oppfinnelsen i sann tid og uten store krav til datamaskinminne, har det blitt utviklet en inversjonsmetode som trekker fordel av tensorproduktstrukturen til kjernen (se ligning (1), (3) og (9) ovenfor). F.eks. kan inversjonsgjenvinning-CPMG-data uttrykkes med en kjerne som er separerbar i ki = l-2e" t1/T1 og k2=e" tO /TO (se ligning (1)); for de drevne likevekt-refokuseringsdata er kjernen separerbar i ki = g(Ti/T2) og k2= e"<x2/T2>(se ligning 3)); og for diffusjonsredigering-CPMG-dataene er kjernen separerbar i ki = e-Dl,2g2Tl3/6 og k2= e<x2/T2>(se ligning (9)).
Ved implementering av den inverse fremgangsmåte blir hvert datasett diskretisert og skrevet i matrisenotasjon:
hvor R er antallet datasett, og, for de diskretiserte parametere x(<l>), i = 1,...,NXog y(<l>), i = l,...,Ny, Mr er en matrise som inneholder det r-te datasett; matrisene Ki og K2inneholder inngangsverdier som korresponderer med kontinuerlige kjerner ki og k2; matrise Fr er den diskretiserte versjon av en kontinuerlig sannsynlighetstetthetsfunksjon fr(x,y); og matrisen Er representerer støyen som er forbundet med dataene. For å løse med henblikk på Fr, anses ligning (11) som et optimaliseringsproblem som utsettes for en ikke-negativitetsbundethet:
hvor ||-||<2>betegner Frobenius-normen for en matrise. Det første uttrykk i ligning (12) måler differansen mellom dataene og tilpasningen, og det annet uttrykk som er summert over alle datasettene, r = 1,...,R, er et mål for glattheten av den diskretiserte tetthetsfunksjon, Fr. Parameteren, a, er en glattings- eller regulariseringsparameter som styrer glattheten av den estimerte tetthetsfunksjon Fr langs hver dimensjon av dataene, inkludert langs r = 1,...,R.
Inversjonsmetoden involverer først komprimering av dataene Mr langs hver dimensjon av kjernen ved bruk feks. av singulær verdidekomponering av Ki og K2eller én av flere metoder til å beregne gjennomsnitt, som er kjent innen faget, så som vindussummer, osv. Tensorproduktstrukturen til kjernen tillater effektiv komprimering av dataene, som har en tilbøyelighet til å være sterkt overtallige. Datakomprimering reduserer typisk størrelsen av dataene med minst ca. tusen ganger uten å modifisere problemets struktur, og muliggjør en effektiv estimering av den diskretiserte tetthetsfunksjon Fr. I rommet med lavere dimensjon for det komprimerte datasett M, kan tetthetsfunksjonen Fr estimeres ved bruk av én av flere kjente fremgangsmåter, så som minste kvadraters tilpasning, Butler-Reed-Dawson (BRD), osv. F.eks., ved anvendelse av BRD-metoden, blir det bundede optimaliseringsproblem for estimering av Fr (se ligning (12) ovenfor) konvertert til et ikke-bundet optimaliseringsproblem for estimering av en vektor cr, r = 1,...,R, hvis dimensjon er lik størrelsen av de komprimerte data. Vektoren cr kan uttrykkes som:
hvor K er tensorproduktet av K1og K2, som representerer de komprimerte kjerner langs den første henholdsvis den annen dimensjon; og m og F er vektorer som er fremkommet ved leksikografisk ordning av matrisene M henholdsvis Fr. Så snart cr er estimert, kan vektoren Fr estimeres som Fr = max(0, K, cr), hvorfra matrisen Fr kan bestemmes. Optimaliseringen gjøres iterativt inntil den estimerte løsning for Fr konvergerer innenfor en ønsket grense.
Som tidligere nevnt regulerer glattingsparameteren a glattheten til den estimerte tetthetsfunksjon Fr langs hver dimensjon av dataene, inkludert langs r = 1,...,R. En optimal verdi av glattingsparameteren a kan bestemmes ved bruk av én av et antall forskjellige strategier kjent innen faget, så som adhoc, BRD, S-kurve kryss-validering, Morozovs uoverenstemmelsesprinsipp, minimum avgrensning, prediktiv midlere kvadratestimering, osv. F.eks. er anvendelsen av en BRD-metode for å optimalisere a for 2- til 3-dimensjonale data en utvidelse av metoden beskrevet i J. P. Butler et al., «Estimating Solutions of the First Kind Integral Equations with Nonnegative Constraints and Optimal Smoothing», SIAM J. Numerical Analysis, vol. 18, nr. 3, sider 381-397, juni 1981.
Denne inversjonsmetoden kan anvendes på alle to- til tre-dimensjonale NMR-data som generelt kan uttrykkes som:
hvor Me(xi, t2) representerer det r-te datasett; ki og k2representerer kjernen atskilt langs en første dimensjon henholdsvis en annen dimensjon; Ti og t2er en første henholdsvis en annen tid, forbundet med den nukleære magnetiske resonansmåling;
x og y er parametere som er relatert til systemet av spinn, /r(x,y) er en felles sannsynlighetstetthetsfunksjon mellom x og y over de R datasett (/r(x,y) kan også anses som en tredimensjonal sannsynlighetstetthetsfunksjon av x, y, og en uavhengig, eksperimentelt styrt variabel, r); og Er(Ti, t2) representerer støy som er forbundet med det r-te nukleære magnetiske resonansdatasett. I tilfelle hvor R = 1, så er dataene todimensjonale, og den ovenfor beskrevne fremgangsmåte korresponderer med en todimensjonal inversjon. Denne inversjonsteknikken har anvendbarhet utover de NMR-målinger som her er beskrevet. Anvendt på 2- og 2,5-dimensjonale datasett, tilveiebringer denne inversjonsteknikken en beregningsmessig effektiv og robust metode til å ekstrahere en todimensjonal avbildning av hvilke som helst av to parametere som kan transformeres ut av NMR-dataene, så som Ti, T2, Ti/T2-forholdet, D, viskositet, metning, osv. Ytterligere detaljer ved inversjonsmetoden er beskrevet i Appendix C i US foreløpig søknad 60/220 053, som ved denne referanse herved inkorporeres i sin helhet.
Inversjonsgjenvinning-CPMG-målingene i forbindelse med den ovenfor beskrevne inversjonsmetode er i stand til å tilveiebringe en full avbildning av den todimensjonale tetthetsfunksjon, /(Ti, T2), med et detaljeringsnivå som hittil ikke har vært praktisk oppnåelig, som vist på fig. 7 og 8. Andre magnetfeltpulssekvenser og andre inversjonsmetoder kan brukes ved en NMR-brønnloggeanvendelse for å tilveiebringe noe av den samme informasjon som inversjonsgjenvinning-CPMG-sekvensen.
F.eks. forsyner en drevet likevekts-refokuseringssekvens som beskrevet ovenfor systemet av spinn med en likevektsmagnetisering som avhenger av Ti/T2, og genererer deretter en første serie av magnetiske resonanssignaler som avhenger av Ti/T2og T2(se ligning (3) ovenfor). Analyse av den første serie av magnetiske resonanssignaler sammen med en annen serie av genererte magnetiske resonanssignaler, f.eks. ved påføring av en ventetid fulgt av refokuseringssekvensen, muliggjør bestemmelse av et gjennomsnittelig Ti/T2-forhold, dvs. <Ti/T2>. Et gjennomsnittelig Ti/T2-forhold, hvor gjennomsnittet er beregnet over alle Ti-verdiene, eller hvor gjennomsnittet til og med er beregnet over alle spinnene, tilveiebringer en viss indikasjon på om hvorvidt en prøve inneholder en fuktig fase eller en ikke-fuktig fase. Hovedinformasjonen som befinner seg i den fullstendige todimensjonale Ti/T2tetthetsfunksjon kan således ekstraheres i et eksperiment med en enkelt gjennomkjøring som vil være egnet til bruk ved brønnloggeanvendelser.
Fig. 9 illustrerer en todimensjonal avbildning av <Ti/T2> mot T2utledet fra to sandsteinsprøver, hvorav den ene er mettet med vann og den andre er delvis mettet med vann og olje. Disse avbildningene, som kan ha blitt lagt oppå hverandre på et plott, ble ekstrahert fra data generert ved en DEFT-CPMG-sekvens (som gitt av uttrykk (4) ovenfor og som illustrert på fig. 3) med 5i = 82= 400 us. Ved å ta et forhold mellom DEFT-CPMG-dataene (se ligning (3) sammen med ligning (6) ovenfor) og dataene generert fra en standard CPMG (se ligning (10) ovenfor),
fremkom den forventede verdi for g — , hvorfra gjennomsnittet av <Ti/T2> over g r j
Ti for hver T2 ble utledet. Selv om det ikke er så detaljert som avbildningene over Ti-T2vist på fig. 7 og 8, kan en vannfase og en oljefase klart skjelnes i avbildningen over <Ti/T2> mot T2på fig. 9 for den vann-olje-mettede prøven 901. F.eks., ved å se på T2-komponenter rundt 200 ms og høyere for 901, er <Ti/T2> tilnærmet 1, hvilket indikerer tilstedeværelse av olje. Ved kortere verdier av T2, er
<Ti/T2> for 901 klart større enn 1, og faller sammen med avbildningen for <Ti/T2> mot T2 for den vannmettede prøven 902, hvilket viser at de mindre porene i den vann-olje-mettede prøven som korresponderer til de kortere T2-komponentene er fylt med vann. Signalets intensitet ved hvert punkt gir en indikasjon på antallet spinn og således volumet av fluid som bidrar til signalet ved dette punktet.
Den DEFT-CPMG-sekvensen som ble brukt til å generere avbildningene over
<Ti/T2> mot T2som er vist på fig. 9 ble implementert med magnetfeltpulser ca. 1 kHz utenfor resonansfrekvensen. Det har blitt funnet at DEFT-CPMG-signalet i homogene felt er følsomt for fase- og amplitudefeiltilpasninger i pulsene, slik at puls-ufullkommenheter påvirker det observerte signalet. Forskjellige fasesykliserings-fremgangsmåter kan brukes for å redusere effekten av ufullkomne pulser. Signalets følsomhet kan også reduseres ved bruk av magnetfeltpulser som er forskjøvet fra resonansen med Aco«±(l/8i), og ved å bruke to 180° pulser i 52-perioden istedenfor én. En utførelse av disse modifikasjonene på DEFT-CPMG-sekvensen resulterte i godt samsvar mellom de signalene som ble målt i laboratoriet og ligning (3) og (6) ovenfor.
I inhomogene felt så som de som typisk finnes ved de fleste brønnloggeanvendelser, er enkelte spinn imidlertid nødvendigvis utenfor resonansfrekvensen, og, i den raskt pulserende grense (dvs. ti, t2 » Ti, T2), oppviser DEFT-CPMG-signalet ikke den samme følsomhet overfor pulsufullkommenheter. For å tilpasse de observerte DEFT-CPMG-signalene til ligning (3) sammen med ligning (6) ovenfor, tas imidlertid all opptreden utenfor resonansfrekvensen og detaljene ved datainnsamlingen (så som signalbåndbredde, osv.) fortrinnsvis i betraktning vedberegning av g — .
Det skal bemerkes at likevekten som er drevet av en konstant RF-irradiasjon, som blir generert med en enkelt puls, ikke vil være så følsom for pulsfeiltilpasninger og puls-ufullkommenheter som den DEFT-drevne likevekt, og således være velegnet til bruk sammen med innhomogene felt, så som det som finnes i de fleste brønnloggeanvendelser.
Oppfinnelsen kan implementeres i brønnlogging ved bruk av enhver brønnloggeanordning for nukleær magnetisk resonans (NMR) som er kjent innen faget. Utførelser av oppfinnelsen kan implementeres med NMR-brønnlogge-innretninger uten behov for maskinvaremodifikasjoner. Fig. 10 viser et flytskjema over en rutine som kan brukes ved programmering av en prosessor for implementering av visse utførelser av oppfinnelsen. Rutinen kan lagres på eller anordnes på en datamaskin eller et maskinlesbart medium, så som et leseminne (ROM); et direkteminne (RAM); magnetplate eller -bånd; en CD-ROM eller et annet optisk lagringsmedium; elektriske, optiske, akustiske eller andre former for overførte signaler; og lignende. Prosessoren kan være en nedihullsprosessor, en opphullsprosessor, eller en kombinasjon av dette. Prosessoren kan også inkludere en fjerntliggende prosessor som kan brukes til å implementere noe av datainversjonen og tolkningsdelene av rutinen.
Før påbegynnelse av den programmerte rutine, og som vist ved 110, påføres et statisk magnetfelt på et undersøkelsesområde innenfor en jordartformasjon som omgir et borehull. Teoretisk kan jordens magnetfelt brukes som det statiske magnetfelt. For de fleste praktiske formål er imidlertid dette ikke foretrukket. Det statiske magnetfelt og feltgradienten (når dette er ønskelig) påføres typisk ved bruk av et loggeverktøy som har en permanentmagnet eller en oppstilling av permanentmagneter.
Den programmerte rutine begynner ved blokk 1015, som representerer initialisering av parametere for magnetfelt-pulssekvensen. Sekven sp aram eterne kan feks. inkludere i=l, en gjenvinningstid, en ventetid, Ti, t2, osv., avhengig av pulssekvensen som skal påføres. I enkelte utførelser kan også fasesyklisering av instruksjoner innføres i den programmerte rutine ved 1015. Generering av en magnetfeltpulssekvens i undersøkelsesområdet er representert med blokk 1020. Magnetiske resonanssignaler fra undersøkelsesområdet detekteres i blokk 1022 og lagres.
Enkelte utførelser av oppfinnelsen involverer gjentagende påføring av magnetfeltpulssekvensen, eller involverer påføring av et flertall av magnetfeltpulssekvenser. En parameter R kan brukes (settes kanskje ved blokk 1015) for å angi et samlet antall magnetfeltpulssekvenser som skal genereres og påføres. Beslutningsblokk 1025 representerer forespørsel av om hvorvidt en teller, r, er lik R. Hvis ikke inkrementeres r, som representert ved blokk 1026, og, om nødvendig, varieres sekvensparametere, så som gjenvinningstid eller fasesykliseringsinstruksjoner, som representert i blokk 1027, før rutinen returneres til blokk 1020 hvor den neste sekvensen genereres i undersøkelsesområdet. Hvis r=R ved forespørselen i beslutningsblokken 1025, er datainnsamlingen for målingen fullført, og rutinen fortsetter til blokk 1030 hvor datainversjonen begynner. Avhengig av den bestemte NMR-måling som blir utført, kan flere enn én teller og flere enn én iterative sløyfer være nødvendig.
Datainversjonen begynner i blokk 1030 ved å sette inversjonsparameterne, så som rutenettpresisjonen som skal brukes ved evaluering av F, tetthetsfordelingen; typen kjerne; typen optimalisering, støy, osv. I enkelte utførelser blir kjernen og dataene først komprimert, og det første mål ved inversjonen er å ekstrahere tetthetsfordelingen, F, fra dataene, som vist i blokkene 1035-1045. Blokk 1035 representerer komprimering av kjernen langs hver dimensjon i kjernen. I tilfellet hvor kjernen har to dimensjoner, kan 2D-kjernen komprimeres, f.eks. ved først å evaluere hver lD-kjerne med en singulær verdi dekomposisjon, og deretter beregne den komprimerte 2D-kjernen. Kjernen blir typisk komprimert én gang, og anvendt på hvert datasett, slik at når dataene komprimeres i blokk 1040, behøver den komprimerte kjernen ikke å beregnes på ny for hvert datasett. Datasettene Mr, r=l,...,R, blir komprimert langs hver dimensjon av kjernen i blokk 1040. Tetthetsfordelingen Fr for r=l,...,R bestemmes ved 1045.1 enkelte utførelser involverer bestemmelse av Fr å betrakte problemet i et ubundet optimaliseringsnettverk i et rom med en lavere dimensjon (som har en dimensjon som er lik dimensjonen for den komprimerte kjernen) og iterativt å optimalisere Fr inntil den konvergerer innenfor en ønsket grense. Når Fr har blitt bestemt, ekstraheres informasjonen om området av interesse, så som Ti, T1/T2, avbildningen av T1-T2, osv. i blokk 1050. I enkelte utførelser kan informasjon om området av interesse, så som avbildningen av <Ti/T2>-T2ekstraheres uten at det er nødvendig å fullstendig invertere for /, som tidligere beskrevet.
Oppfinnelsen gjør at informasjon om et system av nukleære spinn, så som i et fluid i en bergart, kan ekstraheres, enten i en laboratorieomgivelse eller ved brønnloggeanvendelser. Enkelte utførelser kan brukes til å ekstrahere todimensjonale avbildninger av parametere av interesse, så som Ti og T2, Ti/T2og T2, D og T2, osv., hvilket kan vise seg å være nyttig for å skjelne mellom fluider og/eller mellom forskjellige poreomgivelser i en bergart. Selv om oppfinnelsen her har blitt beskrevet med henvisning til visse eksempler og utførelser, vil det være klart at forskjellige modifikasjoner og endringer kan gjøres med de ovenfor beskrevne utførelser uten å avvike fra oppfinnelsens ramme og idé slik dette er angitt i kravene.
Claims (8)
1. En fremgangsmåte for ekstrahering av informasjon om et system av nukleære spinn, omfattende; a) utførelse av et flertall av nukleære magnetiske resonansmålinger på systemet av nukleære spinn ved gjentakelse av en sekvens av magnetiske feltpulser anvendt på systemet av nukleære spinn, der hver sekvens omfatter: - en første delkarakterisert veden første parameter t1; der nevnte første del er utviklet for å klargjøre systemet av nukleære spinn i en tilstand som avhenger av en første systemparameter x, og at nevnte første del er fulgt av - en andre delkarakterisert veden andre parameter t2, der nevnte andre del er utviklet for å generere en serie av spinn-ekkoer som avhenger av nevnte tilstand og på en andre systemparameter y, der hvert av nevnte spinn-ekkoer erkarakterisert veden verdi på nevnte andre parameter t2;b) innsamling av nukleære magnetiske resonansdata fra hver av de flertall av nukleære magnetiske resonansmålinger, idet de nukleære magnetiske resonansdata er uttrykt som MT( xu x2) =llkl( x, xl) k2( y, X2) fi( x, y) dxdy+ Er( xux2), hvor
MT( xi, x2) representerer de nukleære magnetiske resonansdata for et flertall av verdiene til nevnte første parameter t^
ki og k2er kjerner;
fr( x, y) er en to-dimensjonal sannsynlighetstetthetsfunksjon av x ogy; og £V(ti,t2) representerer støy forbundet med de nukleære magnetiske resonansdata; c) utrykking av de nukleære magnetiske resonansdata som M^ KjF^+ Er, hvor matrisene Ktog K2inneholder inndata som tilsvarer henholdsvis ktog k2, og Fr representerer en diskretisert versjon av fr( x, y) ; d) komprimering av de nukleære magnetiske resonansdata og komprimering av matrisene Kj og K2; e) invertering av de komprimerte nukleære magnetiske resonansdata for å løse to-dimensjonal sannsynlighetstetthetsfunksjon Fr.
2. Fremgangsmåten ifølge krav 1, idet komprimeringen av de nukleære magnetiske resonansdata omfatter omforming av dataene ved bruk av en singulærverdi-dekomposisjon av kjernene.
3. Fremgangsmåte ifølge krav 1, idet inverteringen av de komprimerte nukleære magnetiske resonansdata benytter en Butler-Reed-Dawson optimaliseringsmetode.
4. Fremgangsmåte ifølge krav 1, idet nevnte tilstand omfatter en ikke-likevektstilstand eller en termisk likevektstilstand, avhengig av den første parameter Tl.
5. Fremgangsmåte ifølge krav 1, idet system-parameterene x og y avhenger av et fluid i en bergart.
6. Fremgangsmåte ifølge krav 5, idet system-parameterene x og y omfatter en longitudinal relaksasjonstid som er forbundet med fluidet.
7. Fremgangsmåte ifølge krav 5, idet system-parameterene x og y omfatter en transversal relaksasjonstid som er forbundet med fluidet.
8. En loggeanordning omfattende: et loggeverktøy som er bevegelig gjennom et borehull, og en prosessor som er forbundet med loggeverktøyet,
idet prosessoren er programmert med instruksjoner som, når de utføres av prosessoren, forårsaker at loggeverktøyet utfører fremgangsmåten ifølge ethvert av kravene 1 til 4.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US22005300P | 2000-07-21 | 2000-07-21 | |
US09/723,803 US6570382B1 (en) | 1999-12-10 | 2000-11-28 | Nuclear magnetic resonance method and logging apparatus |
PCT/US2001/023012 WO2002008795A2 (en) | 2000-07-21 | 2001-07-20 | Nuclear magnetic resonance measurements and methods of analyzing nuclear magnetic resonance data |
Publications (3)
Publication Number | Publication Date |
---|---|
NO20030223D0 NO20030223D0 (no) | 2003-01-16 |
NO20030223L NO20030223L (no) | 2003-03-20 |
NO336562B1 true NO336562B1 (no) | 2015-09-28 |
Family
ID=26914524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NO20030223A NO336562B1 (no) | 2000-07-21 | 2003-01-16 | Analyse av data fra målinger med nukleær magnetisk resonans |
Country Status (5)
Country | Link |
---|---|
US (1) | US6462542B1 (no) |
EP (1) | EP1301808B1 (no) |
AU (1) | AU2001280678A1 (no) |
NO (1) | NO336562B1 (no) |
WO (1) | WO2002008795A2 (no) |
Families Citing this family (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6891369B2 (en) * | 1998-08-13 | 2005-05-10 | Schlumberger Technology Corporation | Nuclear magnetic resonance method and logging apparatus for fluid analysis |
US6518757B1 (en) * | 2002-03-08 | 2003-02-11 | Schlumberger Technology Corporation | Use of CPMG sequences with phase cycled refocusing pulses in inside-out NMR for phase encoded imaging and to eliminate coherent ringing within one scan |
US6883702B2 (en) * | 2002-03-21 | 2005-04-26 | Schlumberger Technology Corporation | Method and apparatus for NMR measurement of wettability |
US6838875B2 (en) * | 2002-05-10 | 2005-01-04 | Schlumberger Technology Corporation | Processing NMR data in the presence of coherent ringing |
WO2004021028A1 (en) * | 2002-08-27 | 2004-03-11 | Kennedy Krieger Institute | Microvascular blood volume magnetic resonance imaging |
FR2844355B1 (fr) * | 2002-09-11 | 2005-03-25 | Inst Francais Du Petrole | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
US6960913B2 (en) * | 2003-01-14 | 2005-11-01 | Schlumberger Technology Corporation | Multi-measurement NMR analysis based on maximum entropy |
US6894493B2 (en) * | 2003-01-16 | 2005-05-17 | Schlumberger Technology Corporation | Method and apparatus for NMR measurement of magnetic materials |
US6937014B2 (en) * | 2003-03-24 | 2005-08-30 | Chevron U.S.A. Inc. | Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins |
US6952096B2 (en) * | 2003-09-05 | 2005-10-04 | Schlumberger Technology Corporation | Method and apparatus for determining speed and properties of flowing fluids using NMR measurements |
US7053611B2 (en) * | 2004-06-04 | 2006-05-30 | Schlumberger Technology Corporation | Method and apparatus for using pulsed field gradient NMR measurements to determine fluid properties in a fluid sampling well logging tool |
US8238194B2 (en) * | 2004-09-23 | 2012-08-07 | Schlumberger Technology Corporation | Methods and systems for compressing sonic log data |
US20060062081A1 (en) * | 2004-09-23 | 2006-03-23 | Schlumberger Technology Corporation | Methods and systems for compressing sonic log data |
CA2587896A1 (en) * | 2004-11-22 | 2006-06-01 | Baker Hughes Incorporated | Fast t1 measurement of an earth formation by using driven equilibrium |
US7298142B2 (en) * | 2005-06-27 | 2007-11-20 | Baker Hughes Incorporated | Method and apparatus for reservoir fluid characterization in nuclear magnetic resonance logging |
US7500388B2 (en) * | 2005-12-15 | 2009-03-10 | Schlumberger Technology Corporation | Method and apparatus for in-situ side-wall core sample analysis |
US7603237B2 (en) * | 2006-07-26 | 2009-10-13 | Schlumberger Technology Corporation | Method for analyzing data having shared and distinct properties |
US20080147763A1 (en) * | 2006-12-18 | 2008-06-19 | David Levin | Method and apparatus for using state space differential geometry to perform nonlinear blind source separation |
US8717022B2 (en) * | 2008-07-08 | 2014-05-06 | University Of New Brunswick | Magnetic field gradient monitor apparatus and method |
US9052409B2 (en) * | 2008-07-11 | 2015-06-09 | Schlumberger Technology Corporation | Monte Carlo method for laplace inversion of NMR data |
US8427145B2 (en) | 2010-03-24 | 2013-04-23 | Schlumberger Technology Corporation | System and method for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis |
US20120074934A1 (en) | 2010-09-29 | 2012-03-29 | Derrick Green | Nmr measurements and methods of analyzing nmr data |
WO2012103397A2 (en) * | 2011-01-28 | 2012-08-02 | Schlumberger Technology B.V. | Two dimensional nmr of diffusion and relaxation for material characterization |
US9658358B2 (en) | 2011-12-27 | 2017-05-23 | Schlumberger Technology Corporation | Refocusing pulses and excitation pulses for NMR logging |
US9222902B2 (en) | 2012-01-09 | 2015-12-29 | Schlumberger Technology Corporation | Estimations of nuclear magnetic resonance measurement distributions |
US9201158B2 (en) | 2012-01-24 | 2015-12-01 | Schlumberger Technology Corporation | Estimating and displaying molecular size information of a substance |
WO2014026185A1 (en) * | 2012-08-10 | 2014-02-13 | Maui Imaging, Inc. | Calibration of multiple aperture ultrasound probes |
RU2517762C2 (ru) * | 2012-09-12 | 2014-05-27 | Общество с ограниченной ответственностью "МИКС" | Импульсная последовательность для измерения параметров самодиффузии методом ядерного магнитного резонанса |
CN103116148B (zh) * | 2013-01-30 | 2015-04-01 | 上海理工大学 | 一种核磁共振二维谱反演的方法 |
US9678185B2 (en) | 2013-03-15 | 2017-06-13 | Pepsico, Inc. | Method and apparatus for measuring physico-chemical properties using a nuclear magnetic resonance spectrometer |
EP3052970A1 (en) | 2013-12-19 | 2016-08-10 | Halliburton Energy Services, Inc. | Pore size classification in subterranean formations based on nuclear magnetic resonance (nmr) relaxation distributions |
US10031255B2 (en) | 2014-03-24 | 2018-07-24 | Schlumberger Technology Corporation | Multi-dimensional nuclear magnetic resonance methods for characterizing fluids |
CN103984033B (zh) * | 2014-06-09 | 2017-01-11 | 桂林电子科技大学 | 一种地面核磁共振二维反演方法 |
US9823205B2 (en) | 2014-11-17 | 2017-11-21 | Schlumberger Technology Corporation | Methods and systems for determining surface relaxivity of a medium using nuclear magnetic resonance |
US9851315B2 (en) | 2014-12-11 | 2017-12-26 | Chevron U.S.A. Inc. | Methods for quantitative characterization of asphaltenes in solutions using two-dimensional low-field NMR measurement |
US11080672B2 (en) | 2014-12-12 | 2021-08-03 | Ecoatm, Llc | Systems and methods for recycling consumer electronic devices |
DE102015201985B4 (de) * | 2015-02-05 | 2018-11-15 | Siemens Healthcare Gmbh | Verwendung derselben MR-Signale zur Bestimmung verschiedener Parameter |
MX2018000120A (es) * | 2015-07-31 | 2018-03-22 | Halliburton Energy Services Inc | Aparato y metodo para procesar e interpretar datos de adquisicion de registros de resonancia magnetica nuclear (rmn). |
CN105550442B (zh) * | 2015-12-14 | 2019-05-31 | 中国科学院电子学研究所 | 基于瞬变电磁矩变换的数据处理及三维正演方法 |
SE538834C2 (sv) | 2015-12-29 | 2016-12-20 | Cr Dev Ab | Method of extracting information about a sample by nuclear magnetic resonance measurements |
US10634746B2 (en) | 2016-03-29 | 2020-04-28 | Chevron U.S.A. Inc. | NMR measured pore fluid phase behavior measurements |
US10605950B2 (en) | 2016-04-20 | 2020-03-31 | Halliburton Energy Services, Inc. | Downhole formation evaluation using nuclear magnetic resonance inversion matrix reduction |
CA3125967A1 (en) | 2019-01-08 | 2020-07-16 | Schlumberger Canada Limited | Methods and systems for characterizing a porous rock sample employing combined capillary pressure and nmr measurements |
US11187766B1 (en) | 2020-07-09 | 2021-11-30 | Saudi Arabian Oil Company | Methods and systems for determining fluid content in formation samples using low field nuclear magnetic resonance |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5055788A (en) | 1986-08-27 | 1991-10-08 | Schlumberger Technology Corporation | Borehole measurement of NMR characteristics of earth formations |
US5291137A (en) * | 1992-11-02 | 1994-03-01 | Schlumberger Technology Corporation | Processing method and apparatus for processing spin echo in-phase and quadrature amplitudes from a pulsed nuclear magnetism tool and producing new output data to be recorded on an output record |
US5363041A (en) | 1992-12-31 | 1994-11-08 | Schlumberger Technology Corporation | Determining bound and unbound fluid volumes using nuclear magnetic resonance pulse sequences |
US5680043A (en) | 1995-03-23 | 1997-10-21 | Schlumberger Technology Corporation | Nuclear magnetic resonance technique for determining gas effect with borehole logging tools |
US5977768A (en) * | 1997-06-23 | 1999-11-02 | Schlumberger Technology Corporation | Nuclear magnetic resonance logging with azimuthal resolution |
US6111409A (en) * | 1998-03-02 | 2000-08-29 | Western Atlas International, Inc. | Nuclear magnetic reasonance fluid characterization apparatus and method for using with electric wireline formation testing instruments |
US6184681B1 (en) * | 1998-03-03 | 2001-02-06 | Schlumberger Technology Corporation | Apparatus and method for computing a distribution of spin-spin relaxation times |
US6255818B1 (en) * | 1998-08-18 | 2001-07-03 | Schlumberger Technology Corporation | Method and apparatus for performing magnetic resonance measurements |
-
2001
- 2001-07-20 WO PCT/US2001/023012 patent/WO2002008795A2/en active Application Filing
- 2001-07-20 AU AU2001280678A patent/AU2001280678A1/en not_active Abandoned
- 2001-07-20 US US09/909,664 patent/US6462542B1/en not_active Expired - Lifetime
- 2001-07-20 EP EP01959088A patent/EP1301808B1/en not_active Expired - Lifetime
-
2003
- 2003-01-16 NO NO20030223A patent/NO336562B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
EP1301808B1 (en) | 2009-11-18 |
NO20030223D0 (no) | 2003-01-16 |
US20020067164A1 (en) | 2002-06-06 |
WO2002008795A3 (en) | 2002-05-23 |
EP1301808A2 (en) | 2003-04-16 |
NO20030223L (no) | 2003-03-20 |
US6462542B1 (en) | 2002-10-08 |
WO2002008795A2 (en) | 2002-01-31 |
AU2001280678A1 (en) | 2002-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NO336562B1 (no) | Analyse av data fra målinger med nukleær magnetisk resonans | |
NO336429B1 (no) | Kjernefysiske magnetressonansmetoder for å trekke ut informasjon om et fluid i en bergart | |
US5680043A (en) | Nuclear magnetic resonance technique for determining gas effect with borehole logging tools | |
US6859033B2 (en) | Method for magnetic resonance fluid characterization | |
NO20110236A1 (no) | Monte Carlo-metoden for Laplace-inversjonen av NMR-data | |
NO327685B1 (no) | Kjernemagnetisk resonansmalinger i bronnlegging som benytter en optimalisert refaseregulerende pulssekvens | |
NO336430B1 (no) | Bestemmelse av T1 relaksasjonstid fra multippel ventetid NMR logger frembragt i den samme eller forskjellige passeringer | |
NO315065B1 (no) | Nukle¶r magnetisk resonansloggeapparat for borehull og fremgangsmåte ved logging for å bestemme et volum av hydrokarboner uavhengig av endiffusjonskoeffisient | |
NO322672B1 (no) | Fremgangsmate for a bestemme mengde av bundet fluidum i en grunnformasjon | |
NO312266B1 (no) | Fremgangsmåte for bestemmelse av permeabilitet i poröst medium, spesielt et fluidummettet medium | |
NO337897B1 (no) | Fremgangsmåte og apparat for å benytte NMR-målinger med pulset feltgradient for å bestemme fluidegenskaper i et fluidprøvetakende brønnloggeverktøy | |
NO321432B1 (no) | Fremgangsmate for NMR-malinger av sensitive volumer ved en formasjon ved bruk av et flertall CPMG-ekkotog | |
NO335581B1 (no) | Tilegnelse av NMR-data med flere tidsintervall mellom ekkoene | |
NO335578B1 (no) | Bestemmelse av fuktningsevne i et oljereservoar ved bruk av NMR-borehullsmålinger | |
NO342538B1 (no) | Metodologi for tolkning og analyse av NMR-avbildninger | |
NO336377B1 (no) | Fremgangsmåte og apparatur for måling av fuktningsegenskapene ved kjernemagnetisk resonans | |
WO2013148516A1 (en) | Nuclear magnetic resonance testing for organics and fluids in source and reservoir rock | |
US8532929B2 (en) | Method and apparatus to incorporate internal gradient and restricted diffusion in NMR inversion | |
CN107748126A (zh) | 一种获取岩石孔隙尺寸和孔隙表面弛豫率的核磁共振方法 | |
NO326723B1 (no) | Teknikk for å oppnå estimater med høy oppløsning for petrofysiske egenskaper | |
GB2386190A (en) | Nmr logging based on ssfp | |
NO20141005A1 (no) | Ikke-uniform ekkotog-desimering | |
Mitchell | Can sodium NMR provide more than a tracer for brine in petrophysics? | |
WO2018000044A1 (en) | Downhole diffusion coefficient measurement | |
US10408964B2 (en) | Multiple detector NMR tool system and inversion methods for NMR logging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MK1K | Patent expired |