BRPI0907021B1 - método de avaliação de uma formação terrestre, aparelho para avaliar uma formação terrestre e produto de meio legível por computador - Google Patents

método de avaliação de uma formação terrestre, aparelho para avaliar uma formação terrestre e produto de meio legível por computador Download PDF

Info

Publication number
BRPI0907021B1
BRPI0907021B1 BRPI0907021A BRPI0907021A BRPI0907021B1 BR PI0907021 B1 BRPI0907021 B1 BR PI0907021B1 BR PI0907021 A BRPI0907021 A BR PI0907021A BR PI0907021 A BRPI0907021 A BR PI0907021A BR PI0907021 B1 BRPI0907021 B1 BR PI0907021B1
Authority
BR
Brazil
Prior art keywords
nmr
profiles
sector
distribution
profile
Prior art date
Application number
BRPI0907021A
Other languages
English (en)
Inventor
Romero Pedro
Original Assignee
Baker Hughes 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
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=40898562&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=BRPI0907021(B1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by Baker Hughes Inc filed Critical Baker Hughes Inc
Publication of BRPI0907021A2 publication Critical patent/BRPI0907021A2/pt
Publication of BRPI0907021B1 publication Critical patent/BRPI0907021B1/pt

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/32Electric 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/084Detection of potentially hazardous samples, e.g. toxic samples, explosives, drugs, firearms, weapons
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/3808Magnet assemblies for single-sided MR wherein the magnet assembly is located on one side of a subject only; Magnet assemblies for inside-out MR, e.g. for MR in a borehole or in a blood vessel, or magnet assemblies for fringe-field MR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/383Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using permanent magnets
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/50NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Health & Medical Sciences (AREA)
  • Geophysics (AREA)
  • Remote Sensing (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Toxicology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Vascular Medicine (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Molecular Biology (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

método de avaliação de uma formação terrestre, aparelho para avaliar uma formação terrestre e produto de meio legível por computador sequências de pulso são aplicadas em um fluido em uma formação terrestre em um campo magnético estático e sinais de eco por rotação de rmn são obtidos. os sinais são invertidos para produzir distribuições t2 em uma pluralidade de profundidades. as similaridades entre os perfis dos setores de t2 com resistividade e/ou perfis de raios gama são usadas para identificar e subtrair contribuições para o sinal de rmn do petróleo tendo gradientes internos. a partir dos sinais recebidos, características de relaxamento e difusão do fluido são determinadas. a determinação considera os gradientes do campo interno.

Description

Relatório Descritivo da Patente de Invenção para MÉTODO DE AVALIAÇÃO DE UMA FORMAÇÃO TERRESTRE, APARELHO PARA AVALIAR UMA FORMAÇÃO TERRESTRE E PRODUTO DE MEIO LEGÍVEL POR COMPUTADOR.
CAMPO DA INVENÇÃO [001] A presente invenção refere-se a aparelho e técnicas para fazer medições de ressonância magnética nuclear (RMN) em poços e a métodos para determinar as características magnéticas das formações atravessadas por um poço. Especificamente, a invenção referese ao projeto de uma técnica precisa de medição de RMN, que é planejada na definição das propriedades de difusão do fluido na matriz do poro na presença dos gradientes do campo magnético interno. ANTECEDENTES DA INVENÇÃO [002] Uma variedade de técnicas foi usada na determinação da presença e na estimativa das quantidades de hidrocarbonetos (óleo e gás) nas formações terrestres. Esses métodos são projetados para determinar parâmetros de interesse, incluindo entre outras coisas, porosidade, conteúdo do fluido e permeabilidade da formação de rocha que circunda o poço perfurado para recuperar os hidrocarbonetos. Tipicamente, as ferramentas projetadas para fornecer a informação desejada são usadas para fazer o perfil do poço. Muita da perfilagem é feita depois que os poços foram perfurados. Mais recentemente, poços foram perfilados enquanto na perfuração dos poços. Isso é citado como medição sem interromper a perfuração (MWD) ou perfilagem sem interromper a perfuração (LWD). As medições também foram feitas quando executando a manobra de saída de uma coluna de perfuração de um poço: isso é chamado medição sem interromper a manobra (MWT).
[003] Uma técnica em desenvolvimento usa ferramentas de perfilagem de ressonância magnética nuclear (RMN) e métodos para de
Petição 870190094158, de 20/09/2019, pág. 6/57
2/24 terminar, entre outras coisas, a porosidade, saturação do hidrocarboneto e permeabilidade das formações de rocha. Ferramentas de perfilagem de RMN excitam os núcleos dos fluidos nas formações geológicas na proximidade do poço, de modo que certos parâmetros podem ser estimados, tais como densidade de rotação, tempo de relaxamento longitudinal (geralmente citado na técnica como Ti) e tempo de relaxamento transversal (geralmente citado como T2) das formações geológicas. A partir de tais medições, parâmetros de formação são determinados, tais como porosidade, permeabilidade e saturação do hidrocarboneto, o que proporciona informação valiosa sobre a composição das formações geológicas e a quantidade de hidrocarbonetos extraíveis.
[004] Uma ferramenta de RMN típica gera um campo magnético estático Bo na proximidade do poço e um campo oscilante Bi na direção perpendicular a B0. Esse campo oscilante é geralmente aplicado na forma de pulsos de curta duração. A finalidade do campo Bo é polarizar os momentos magnéticos dos núcleos paralelos ao campo estático e a finalidade do campo B1 é girar os momentos magnéticos por um ângulo controlado pela largura tp e a amplitude B1 do pulso oscilante. Para perfilagem de RMN, a sequência mais comum é a sequência Carr-Purcell-Meiboom-Gill (CPMG) que pode ser expressa como
TW - 90 - (τ - 180 - τ - echo)n (1) onde TW é um tempo de espera, 90 é um pulso de inclinação de 90°, 180 é um pulso de refocalizar de 180° e 2τ = TE é o espaçamento entre ecos.
[005] Depois de ser inclinado por 90°, o momento ma gnético gira ao redor do campo estático em uma frequência particular conhecida como a frequência Larmor ω, dada por ω = γΒο, onde Bo é a intensidade do campo do campo magnético estático e γ é a razão giromagnética. Ao mesmo tempo, os momentos magnéticos retornam para a dire
Petição 870190094158, de 20/09/2019, pág. 7/57
3/24 ção de equilíbrio (isto é, alinhados com o campo estático) de acordo com um tempo de queda conhecido como o tempo de relaxamento da treliça de rotação ou Ti. Não homogeneidades do campo Bo resultam na defasagem dos momentos magnéticos e para remediar isso, um pulso de 180° é incluído na sequência para refocalizar os momentos magnéticos. Essa refocagem fornece uma sequência de n sinais de eco. Essas sequências de eco são então processadas para fornecer informação sobre os tempos de relaxamento. A Patente US No. 6.466.013 para Hawkes e outros e a Patente US No. 6.429.654 para Itskovich e outros, ambas tendo o mesmo cessionário que a presente invenção, ensinam o uso de sequências de pulso CPMG modificadas com exigências de força reduzidas nas quais o ângulo do pulso de refocagem pode ser menor do que 180°. Essas sequência s de pulso modificadas podem ser citadas como sequências de pulso de refasagem ótima (ORPS).
[006] Também está associado com a rotação dos núcleos moleculares um segundo tempo de relaxamento, T2, chamado o tempo transversal ou de relaxamento de rotação-rotação. No fim de um pulso de inclinação de 90°, todas as rotações são apontad as em uma direção comum perpendicular, ou transversal, ao campo estático e elas todas giram na frequência Larmor. Entretanto, por causa das pequenas flutuações no campo estático induzidas por outras rotações ou impurezas paramagnéticas, as rotações precess em frequências ligeiramente diferentes e a magnetização transversal defasam com um tempo de relaxamento T2.
[007] A interpretação dos dados de perfil ou testemunho de RMN é frequentemente iniciada invertendo a queda de eco de CPMG (ou ORPS) no domínio de tempo em uma distribuição de domínio de parâmetro T2. Em geral, o T2 dos fluidos nas rochas porosas depende da distribuição do tamanho do poro e do tipo e número dos fluidos que
Petição 870190094158, de 20/09/2019, pág. 8/57
4/24 saturam o sistema de poro. Por causa da natureza heterogênea dos meios porosos, quedas de T2 exibem um comportamento multiexponencial. A equação básica que descreve o relaxamento transversal da magnetização nos meios porosos saturados de fluido é
1^1 tmx
M(í)= \p{T2}e^dT2 (2) ^2 mio onde M é a magnetização e os efeitos da difusão na presença de um gradiente de campo magnético não foram considerados. A equação (2) é baseada na suposição que os efeitos da difusão podem ser ignorados. Em um campo magnético gradiente, a difusão faz com que os átomos se movam das suas posições originais para novas, o que também faz com que esses átomos adquiram mudanças de fases diferentes comparado com os átomos que não se moveram. Isso contribui para uma taxa de relaxamento mais rápida.
[008] O efeito dos gradientes de campo é dado por uma equação da forma ^2 1.bulk '^‘isurfixe ^diffusion onde os primeiros dois termos no lado direito são relacionados com o relaxamento de volume e relaxamento de superfície enquanto o terceiro termo é relacionado com o gradiente de campo G por uma equação da forma onde TE é o espaçamento entre ecos, C é uma constante e D é a difusibilidade do fluido.
[009] A patente US N° 6.512.371 para Prammer e outros descreve um sistema de perfilagem de poço e método para detectar a presença e estimar a quantidade de hidrocarbonetos gasosos e líquidos na zona próxima do poço. O sistema usa uma ferramenta de perfilaPetição 870190094158, de 20/09/2019, pág. 9/57
5/24 gem de RMN de múltiplas frequências baseada em gradiente para extrair os componentes do sinal característicos para cada tipo de hidrocarboneto. As medições em frequências diferentes são intercaladas para obter, em uma única passagem de perfilagem, múltiplos fluxos de dados correspondendo com tempos de recuperação diferentes e/ou difusibilidade para o mesmo local na formação.
[0010] Uma das dificuldades principais na definição dos parâmetros de autodifusão do fluido na matriz de poro está relacionada com o fato que fluidos diferentes tendo os mesmos tempos de relaxamento e coeficientes de difusão diferentes não podem ser efetivamente separados. Devido à limitação prática da razão de sinal em relação ao ruído, nenhuma das técnicas de inversão existentes permite uma reconstrução efetiva e estável de ambos os espectros de relaxamento e difusão.
[0011] Outra dificuldade na reconstrução dos espectros de relaxamento e difusão é causada pelos gradientes magnéticos internos. Tipicamente, os valores dos gradientes internos são desconhecidos. Dessa forma, os parâmetros de difusão não podem ser corretamente definidos se os gradientes internos não são considerados em ambos os esquemas de medição e interpretação. Esses foram tratados, por exemplo, na Patente US No. 6.597.171 para Hurlimann e outros e na Patente US No. 5.698.979 para Taicher e outros tendo o mesmo cessionário que o presente pedido e cujos conteúdos são incorporados aqui por referência, e a Patente US No. 7.049.815 para Itskovich e outros, tendo o mesmo cessionário do presente pedido e cujos conteúdos são incorporados aqui por referência. Um aspecto comum de todos os métodos da técnica anterior discutidos acima é que eles não consideram medições por outros perfis no processamento dos dados de RMN e as medições em uma profundidade são processadas substancialmente de maneira independente das medições e outras profunPetição 870190094158, de 20/09/2019, pág. 10/57
6/24 didades.
SUMÁRIO DA INVENÇÃO [0012] Uma modalidade da invenção é um método de avaliação de uma formação terrestre.
[0013] O método inclui transportar uma ferramenta de perfilagem para dentro de um poço, obter sinais de ressonância magnética nuclear (RMN) em uma pluralidade de profundidades no poço e processar os sinais de RMN para obter uma distribuição do tempo de relaxamento em cada uma da pluralidade de profundidades, cada uma das distribuições compreendendo uma pluralidade de setores, produzir, a partir da pluralidade de distribuições do tempo de relaxamento, uma pluralidade de perfis de setor de RMN dos valores do tempo de relaxamento correspondendo com cada um da pluralidade de setores, determinar uma medida da similaridade de um primeiro perfil indicativo do conteúdo de hidrocarboneto na formação com cada um da pluralidade de perfis de setor de RMN e usar a medida de similaridade para cada um da pluralidade de perfis de setor de RMN para identificar um subconjunto da pluralidade de perfis de setor de RMN indicativo de hidrocarboneto. Os sinais de RMN podem ser sinais de eco por rotação e o tempo de relaxamento pode ser um tempo de relaxamento transversal T2. A medida de similaridade pode ser um coeficiente de correlação Pearson e/ou um valor p caracterizando uma probabilidade de um coeficiente de correlação observado. O primeiro perfil pode ser um perfil de raios gama e/ou um perfil de resistividade. A determinação da medida de similaridade pode incluir uma correlação cruzada do primeiro perfil com pelo menos um da pluralidade de perfis de setor de RMN e/ou a determinação da entropia mútua entre o primeiro perfil e pelo menos um da pluralidade de perfis de setor de RMN. O método pode incluir obter sinais adicionais de eco por rotação em cada uma da pluralidade de profundidades enquanto aplicando um gradiente de campo magné
Petição 870190094158, de 20/09/2019, pág. 11/57
7/24 tico externo, processar os sinais de eco por rotação e os sinais adicionais de eco por rotação para obter uma distribuição de um coeficiente de difusão (D) em cada uma da pluralidade de profundidades, cada uma da distribuição dos coeficientes de difusão compreendendo uma pluralidade de setores de difusão, produzir, a partir da pluralidade de distribuições do coeficiente de difusão, uma pluralidade de perfis de setor de coeficiente de difusão de valores de setor de coeficiente de difusão correspondendo com cada um da pluralidade dos setores de difusão, determinar uma medida da similaridade do primeiro perfil com cada um da pluralidade de perfis de setor do coeficiente de difusão e usar a medida de similaridade para cada um da pluralidade de perfis de setor de coeficiente de difusão para identificar um subconjunto dos perfis indicativo do hidrocarboneto. O método pode ainda incluir usar o subconjunto identificado dos perfis e uma distribuição T2 em pelo menos uma da pluralidade de profundidades para obter uma distribuição T2 modificada responsiva primariamente à água na formação. O método pode ainda incluir estimar da distribuição modificada, a saturação da água, a água retida na argila, a água retida irredutível, a água móvel e/ou a permeabilidade. O método pode ainda incluir transportar a ferramenta de perfilagem para dentro do poço em um dispositivo de transporte selecionado de um cabo de aço e um tubular de perfuração. [0014] Outra modalidade da invenção é um aparelho para avaliar uma formação terrestre. O aparelho inclui uma ferramenta de perfilagem configurada para ser transportada para dentro de um poço, a ferramenta de perfilagem ainda configurada para obter sinais de ressonância magnética nuclear (RMN) em uma pluralidade de profundidades no poço. O aparelho também inclui pelo menos um processador configurado para processar os sinais de RMN para obter uma distribuição do tempo de relaxamento em cada uma da pluralidade de profundidades, cada uma das distribuições compreendendo uma pluralidade
Petição 870190094158, de 20/09/2019, pág. 12/57
8/24 de setores e usar a pluralidade de distribuições do tempo de relaxamento para produzir uma pluralidade de perfis de setor de RMN dos valores do tempo de relaxamento correspondendo com cada um da pluralidade de setores, determinar uma medida da similaridade de um primeiro perfil indicativo do conteúdo de hidrocarboneto na formação com cada um da pluralidade de perfis de setor de RMN e usar a medida de similaridades para cada um da pluralidade de perfis de setor de RMN para identificar um subconjunto dos perfis de setor de RMN indicativo de hidrocarboneto. Os sinais de RMN podem ser sinais de eco por rotação e o tempo de relaxamento pode ser um tempo de relaxamento transversal (T2). O processador pode ser configurado para usar uma medida de similaridade, que é um coeficiente de correlação Pearson e/ou um valor p caracterizando uma probabilidade de um coeficiente de correlação observado. O aparelho pode ainda incluir uma ferramenta de perfilagem de resistividade e/ou uma ferramenta de perfilagem de raios gama e o processador pode ser configurado para usar como o primeiro perfil, um perfil de raios gama e/ou um perfil de resistividade. O processador pode ser configurado para determinar a medida de similaridade usando uma correlação cruzada do primeiro perfil com pelo menos um da pluralidade de perfis de setor de RMN e/ou a determinação da entropia mútua entre o primeiro perfil e pelo menos um da pluralidade de perfis de setor de RMN. A ferramenta de perfilagem pode ainda ser configurada para obter sinais adicionais de eco por rotação em cada uma da pluralidade de profundidades enquanto aplicando um gradiente de campo magnético externo e o processador pode ser ainda configurado para processar os sinais de eco por rotação e os sinais adicionais de eco por rotação para obter uma distribuição de um coeficiente de difusão em cada uma da pluralidade de profundidades, cada uma da distribuição dos coeficientes de difusão compreendendo uma pluralidade de setores de difusão; produzir, a
Petição 870190094158, de 20/09/2019, pág. 13/57
9/24 partir da pluralidade de distribuições do coeficiente de difusão, uma pluralidade de perfis de setor de coeficiente de difusão de valores de setor de coeficiente de difusão correspondendo com cada um da pluralidade dos setores de difusão; determinar uma medida da similaridade do primeiro perfil com cada um da pluralidade de perfis de setor do coeficiente de difusão e usar a medida de similaridade para cada um da pluralidade de perfis de setor de coeficiente de difusão identificar um subconjunto dos perfis indicativo de hidrocarboneto. O processador pode ser ainda configurado para usar o subconjunto identificado dos perfis e uma distribuição T2 em pelo menos uma da pluralidade de profundidades para obter uma distribuição T2 modificada responsiva primariamente à água na formação. O processador pode ser ainda configurado para estimar da distribuição modificada da saturação da água, água retida na argila, água retida irredutível, água móvel e/ou permeabilidade. O aparelho pode ainda incluir um dispositivo de transporte configurado para transportar a ferramenta de perfilagem para dentro do poço, o dispositivo de transporte sendo um cabo de aço ou um tubular de perfuração.
[0015] Outra modalidade da invenção é um meio legível por computador para uso com um aparelho para avaliar uma formação terrestre. O meio pode ser para uso com um aparelho para avaliar uma formação terrestre na qual o aparelho inclui uma ferramenta de perfilagem configurada para ser transportada para dentro de um poço, a ferramenta de perfilagem ainda configurada para obter sinais de eco por rotação em uma pluralidade de profundidades no poço. O meio inclui instruções que possibilitam que pelo menos um processador processe os sinais de eco por rotação obtidos em uma pluralidade de profundidades no poço na formação terrestre por meio de uma ferramenta de perfilagem para obter uma distribuição do tempo de relaxamento transversal (T2) em cada uma da pluralidade de profundidades, cada
Petição 870190094158, de 20/09/2019, pág. 14/57
10/24 uma das distribuições compreendendo uma pluralidade de setores; usar a pluralidade de distribuições do tempo de relaxamento para produzir uma pluralidade de perfis de setor de RMN dos valores do tempo de relaxamento correspondendo com cada um da pluralidade de setores; determinar uma medida da similaridade de um primeiro perfil indicativo de pelo menos um de: (i) conteúdo de hidrocarboneto ou (ii) conteúdo de água na formação com cada um da pluralidade de perfis de setor de RMN e usar a medida de similaridade para cada um da pluralidade de perfis de setor de RMN para identificar um subconjunto dos perfis de setor de RMN indicativo do conteúdo de pelo menos um de: (i) hidrocarboneto e (ii) água. O meio legível por computador pode ser uma ROM, uma EPROM, uma EAROM, uma memória flash e/ou um disco ótico.
BREVE DESCRIÇÃO DOS DESENHOS [0016] A presente invenção é entendida melhor com referência às figuras seguintes nas quais numerais semelhantes se referem a elementos semelhantes, e nas quais:
[0017] A figura 1 representa diagramaticamente uma ferramenta de perfilagem de RMN em um poço, [0018] A figura 2 (técnica anterior) mostra uma configuração exemplar de ímãs, antena e blindagem adequados para uso com a presente invenção, [0019] A figura 3 mostra um fluxograma ilustrando alguns dos aspectos importantes da presente invenção, [0020] A figura 4 (técnica anterior) mostra uma exibição exemplar de resultados obtidos usando um método da técnica anterior para processamento dos dados de RMN, [0021] A figura 5 mostra uma marcação exemplar de distribuições T2 obtidas pelo processamento dos dados na figura 4 em uma pluralidade de profundidades,
Petição 870190094158, de 20/09/2019, pág. 15/57
11/24 [0022] A figura 6 mostra marcações de resistividade, raios gama, água retida na argila, volume retido irredutível e água móvel através de uma porção do intervalo de profundidade mostrado nas figuras 4 e 5, [0023] A figura 7 mostra marcações de resistividade, raios gama e perfis de setor T2 selecionados sobre o intervalo de profundidade da figura 6, [0024] As figuras 8(a) - (d) ilustram correlações cruzadas exemplares entre perfis medidos em passagens diferentes, [0025] A figura 9(a) mostra uma distribuição T2 exemplar sobre um intervalo de profundidade que inclui petróleo e a figura 9(b) mostra medidas diferentes de similaridade da distribuição T2 com perfis de raios gama e de resistividade, [0026] A figura 10 mostra (a) uma distribuição 2-D dos dados no plano de T2 - difusibilidade, (b) a distribuição em 1-D com relação ao eixo geométrico de T2, (c) a distribuição em 1-D com relação ao eixo geométrico de difusibilidade e (d) uma porção da interface de exibição gráfica, [0027] A figura 11 mostra uma distribuição 2-D dos dados no plano T2 - difusibilidade para o intervalo de profundidade que inclui petróleo, [0028] As figuras 12(a) - (b) mostram o coeficiente de correlação Pearson e o valor p para os dados na figura 11, [0029] A figura 13 mostra a distribuição T2 para os dados da figura 11 codificados em cor de acordo com uma medida da similaridade, [0030] A figura 14 mostra uma distribuição 2-D dos dados no plano T2 - difusibilidade para o intervalo de profundidade que está transportando água, [0031] As figuras 15(a) - (b) mostram o coeficiente de correlação Pearson e o valor p para os dados na figura 14, [0032] A figura 16 mostra as fácies petrográficas e a distribuição T2 correspondente em 100% e saturação de água irredutível,
Petição 870190094158, de 20/09/2019, pág. 16/57
12/24 [0033] A figura 17 mostra resultados de simulação para 100% de saturação de água, [0034] A figura 18 mostra resultados de simulação para 100% de saturação de petróleo e [0035] A figura 19 é um mapa de mobilidade para um reservatório. DESCRIÇÃO DETALHADA DA INVENÇÃO [0036] A figura 1 representa um poço 10 perfurado em um modo típico em uma formação geológica de subsuperfície 12 que será investigada para reservatórios de produção potencial de hidrocarboneto. Uma ferramenta de perfilagem de RMN 14 foi abaixada no poço 10 por meio de um cabo 16 e equipamento de superfície apropriado (representado diagramaticamente por um carretel 18) e está sendo elevada através da formação 12 que compreende uma pluralidade de camadas 12a a 12g de composição diferente, para perfilagem de uma ou mais das características da formação. A ferramenta de perfilagem de RMN pode ser fornecida com molas em arco 22 para manter a ferramenta em uma posição excêntrica dentro do poço com um lado da ferramenta em proximidade com a parede do poço. Os ímãs permanentes 23 produzem o campo magnético estático. Os sinais gerados pela ferramenta 14 são passados para a superfície através do cabo 16 e do cabo 16 através de outra linha 19 para o equipamento de superfície apropriado 20 para processamento, gravação, exibição e/ou para transmissão para outro local para processamento, gravação e/ou exibição. Alternativamente, o processador pode ficar localizado em uma posição adequada (não mostrada) no interior do poço, por exemplo, na ferramenta de perfilagem 14.
[0037] A figura 2 (técnica anterior) ilustra esquematicamente uma modalidade exemplar de um aparelho adequado para uso com o método da presente invenção. Isso é discutido em detalhes na Patente U.S. No. 6.348.792 de Beard e outros, tendo o mesmo cessionário que
Petição 870190094158, de 20/09/2019, pág. 17/57
13/24 a presente invenção e cujos conteúdos são completamente incorporados aqui por referência. A vista de corte da ferramenta na figura 2 ilustra um ímã principal 217, um segundo ímã 218 e uma antena de transceptor compreendendo fios 219 e material de núcleo 210. As setas representam a polarização (por exemplo, do pólo sul para o pólo norte) do ímã principal 217 e do ímã secundário 218. Um aspecto notável da disposição mostrada na figura 2 é que a polarização dos ímãs produzindo o campo estático é para o lado da ferramenta, ao invés de para a frente da ferramenta (o lado direito da figura 2).
[0038] O segundo ímã 218 é posicionado para ampliar a forma do campo magnético estático adicionando um segundo dipólo magnético em proximidade com o dipólo de RF definido pelos fios 219 e o núcleo magnético doce 210. Esse posicionamento move o centro do dipólo estático efetivo mais próximo do dipólo de RF, dessa forma aumentando a extensão azimutal da região de exame. O segundo ímã 218 também reduz o efeito de derivação do núcleo magnético de alta permeabilidade 210 no ímã principal 217. Na ausência do segundo ímã, o campo CC seria efetivamente colocado em curto pelo núcleo 210. Dessa maneira, o segundo ímã, além de agir como um ímã de modelagem para modelar o campo estático para a frente da ferramenta (o lado do ímã principal) também age como um ímã de neutralização com relação ao campo estático no núcleo 210. Esses versados na técnica reconhecerão que a função de neutralização e uma modelagem limitada poderiam ser realizadas simplesmente por terem uma abertura no núcleo; entretanto, desde que algum tipo de modelagem de campo é necessário no lado frontal da ferramenta, em uma modalidade da invenção, o segundo ímã serve tanto para modelagem do campo quanto para a neutralização. Se o campo estático no núcleo 210 fica próximo de zero, então a anelação magnetostritiva do núcleo é substancialmente eliminada.
Petição 870190094158, de 20/09/2019, pág. 18/57
14/24 [0039] Dentro da região de investigação, o gradiente do campo estático é substancialmente uniforme e a intensidade do campo estático se situa dentro de limites predeterminados para produzir uma frequência Larmor substancialmente uniforme. Aqueles versados na técnica reconhecerão que a combinação da modelagem do campo e neutralização poderia ser realizada por outras configurações de ímã do que essas mostradas na figura 2.
[0040] Os sinais de eco de rotação da RMN são obtidos usando o aparelho da figura 2. Opcionalmente, medidas adicionais podem ser feitas usando um campo de gradiente externo como discutido em Reiderman. As medições feitas com o campo de gradiente possibilitam a determinação da difusibilidade. A maneira na qual essas medições são usadas é discutida a seguir.
[0041] É mostrado na figura 3 um fluxograma esboçando algumas das etapas da presente invenção. Medições de RMN são obtidas 303. Além disso, perfis de raios gama indicativos do conteúdo de folhelho da formação e perfis de resistividade são obtidos 301. As medidas dos perfis de resistividade e de raios gama podem ser feitas simultaneamente com a aquisição dos dados da RMN ou, como é mais comum, os perfis de resistividade e/ou raios gama podem ser obtidos em uma execução de perfilagem diferente do que a dos perfis de RMN.
[0042] As medições de RMN são invertidas usando os métodos da técnica anterior para obter uma distribuição T2 em cada uma de uma pluralidade de profundidades 305. A distribuição T2 é caracterizada por seus valores em uma pluralidade de setores. A figura 5 mostra uma marcação exemplar da distribuição T2 em uma pluralidade de profundidades. No exemplo mostrado, aproximadamente 30 setores foram usados com um tempo total maior do que 1024 ms.
[0043] Afastando-se brevemente do assunto principal, a figura 4 mostra uma exibição exemplar do processamento da técnica anterior.
Petição 870190094158, de 20/09/2019, pág. 19/57
15/24
A trilha 1 401 mostra a distribuição T2 de um poço sobre uma faixa de profundidade de aproximadamente 182,9 m (600 pés). A trilha 2 403 mostra o volume retido irredutível 411, a água retida móvel 415 e a água retida na argila 417. A trilha 3 405 mostra a permeabilidade determinada e a trilha 4 407 mostra o desajuste na determinação da permeabilidade e porosidades. A porosidade total inclui a água retida na argila (CBW), água retida capilar (também conhecida como volume irredutível ou BVI), água móvel e hidrocarbonetos. Ver US 6.972.564 para Chen.
[0044] Um aspecto importante da presente invenção é que os dados da RMN são classificados para produzir um perfil de setor da RMN para cada um dos setores da distribuição 307. O método da presente invenção é baseado na identificação e uso de similaridades entre os perfis de setor da RMN e os perfis de resistividade e de raios gama. A figura 6 mostra uma exibição do perfil da resistividade 601 através de um intervalo de 2,44 m (8 pés), o perfil de raios gama 607, um perfil da água retida na argila (da figura 4), uma marcação de BVI 605 e uma marcação de água retida móvel 609. A figura 7 mostra as marcações de raio gama e de resistividade e, além disso, mostra um perfil de setores da RMN 2-8 através do mesmo intervalo de profundidade 701, 703, 705, 707, 709 e 711.
[0045] O presente método é baseado em dois princípios. Primeiro, um aumento no perfil de resistividade é geralmente devido a um aumento no petróleo na formação (desde que o raio gama indica uma pequena quantidade de folhelho). Dessa maneira, nas profundidades nas quais a resistividade é alta, o setor que caracteriza o petróleo na formação tenderá a ter um valor maior. Dessa maneira, com referência novamente à figura 3, um perfil de profundidade é obtido para cada setor da distribuição T2 307. Ainda com referência à figura 3, uma medida da similaridade é determinada entre os perfis de setor T2 e o perfil
Petição 870190094158, de 20/09/2019, pág. 20/57
16/24 de resistividade e o perfil de raios gama 309. A medida da similaridade é discutida abaixo.
[0046] O teste de similaridade pode ser executado usando qualquer correlação linear ou não-linear, dependendo do caso. Entretanto, resultados de testes sugerem o uso da correlação cruzada ou somente da matriz de correlação se não existe mudança de profundidade entre os perfis. Uma vez que o auge do petróleo tenha sido encontrado, a correlação com a viscosidade pode ser estabelecida usando modelos conhecidos. Além do mais, removendo os efeitos do petróleo exemplificado pelo auge do petróleo 311 da distribuição T2 original, a distribuição da água T2 pode ser obtida e, portanto, a saturação da água, Sw. Isso é discutido mais abaixo depois da discussão da figura 14. No caso especial de petróleo pesado, e provavelmente em casos mais gerais, porém dependendo da viscosidade do petróleo e grau de saturação da água, os perfis do setor T2 com correlação negativa com os perfis de resistividade e seus gradientes são mais prováveis de representar a água móvel, que será drenada pelo petróleo quando a saturação do petróleo aumenta. Por esse meio, é possível traçar de volta a distribuição T2 da saturação da água usando um ajuste de curva Gaussiana ou uma abordagem como essa na Patente US 7.363.161 para Georgi e outros que também gera uma escala para correlacionar o setor T2 com o tamanho do grão. Quando um modelo para a distribuição T2 da água é obtido, um modelo para a curva de pressão capilar e qualidade da rocha (P. Romero, SPWLA 2004) e para a permeabilidade relativa (Corey-Burdine) pode ser estabelecido 313. Além do mais, uma marcação cruzada da distribuição T2 de hidrocarboneto vs. a distribuição T2 para os 100% de saturação da água pode ser construída, que indica as zonas de índices de produtividade diferentes. Isso é discutido abaixo com referência às figuras 16-19.
[0047] Esses versados na técnica reconheceríam que perfis feitos
Petição 870190094158, de 20/09/2019, pág. 21/57
17/24 em momentos diferentes no mesmo poço podem não ser perfeitamente registrados. Isso é ilustrado na figura 8. É mostrada na figura 8(a) a correlação cruzada entre os perfis de resistividade e de raios gama. Como esperado, existe uma correlação negativa entre os dois e uma mudança de tempo de zero porque os perfis foram medidos na mesma passagem de perfilagem. A figura 8(b) mostra a correlação cruzada entre os perfis de resistividade e o da água retida na argila. Desde que o último é derivado das medições de RMN feitas em uma passagem de perfilagem separada, existe uma mudança de profundidade de aproximadamente 2,44 m (10 pés) entre os perfis. Uma mudança similar é notada na figura 8(c) do volume retido irredutível e na figura 8(d) para a água móvel. Uma tal etapa de correlação cruzada, embora não ilustrada na figura 3, pode ser necessária antes da determinação da similaridade 309. Como uma alternativa ao uso de um pico da técnica de correlação cruzada para mudanças de tempo diferentes, a entropia mútua entre duas investigações alteradas pode ser usada. A entropia é definida como
H(x) = -E{xlog(P(x))} onde H(.) é a entropia, E{.} é o valor esperado e P(.) é a função da densidade de probabilidade.
[0048] Serão discutidas a seguir as medidas de similaridade que são usadas na presente invenção entre os perfis diferentes. Essas são para finalidades exemplares somente e outras medidas de similaridade poderiam ser usadas com o presente método.
[0049] O primeiro conceito é esse da covariância. Essa é uma medida de quanto as variações das variáveis são correlacionadas. Considere um vetor aleatório
[0050] Para qualquer par de componentes, a covariância é definida pela expectativa cov(X(,yt) = Σλ = E{(X, -A)(Xt - μ.)} (6),
Petição 870190094158, de 20/09/2019, pág. 22/57
18/24
Onde μ, e μκ são os valores médios de Xi e Xk. A covariância é então fornecida pela matriz de covariância:
[0051] O segundo conceito é esse da correlação. A correlação é uma dimensão menor do que a medida da dependência linear das variáveis aleatórias. O coeficiente de correlação Pearson r é a melhor estimativa da correlação de variáveis normalmente distribuídas:
cov(x;,^>
onde gí e Ok são os desvios padrões de Xi e Xk. A matriz de correlação é:
A.» Pu Pu · · P\,n
Pu Pu Pu · P\.»
Pu Pu Pu ·· Pu [0052] Correspondendo com a matriz de correlação existe uma matriz de correlação/valor p.
Pu Pu · · · Pl.n
Pu
Pu Pu
Pu ·· Ρι.„ Pu ··· Pu (10).
[0053] Em termos simples, um valor p é a probabilidade de obter uma descoberta como um resultado do acaso somente.
[0054] Com referência agora à figura 9, é mostrada na figura 9(a) uma distribuição T2 em média exemplar 911 derivada sobre um intervalo de profundidade das medições feitas em um poço. A marcação mostra a T2 em média como uma função do número de setor. 913 é o
Petição 870190094158, de 20/09/2019, pág. 23/57
19/24 coeficiente de correlação Pearson entre os perfis de setor de RMN e o perfil de raios gama, enquanto 915 é o coeficiente de correlação Pearson entre os perfis de setor de RMN e o perfil de resistividade. O valor p da resistividade está representado por 919, enquanto 917 é o valor p dos raios gama. Deve ser enfatizado que os valores para 913, 915, 917 e 919 para qualquer setor particular são obtidos usando o perfil de setor T2 de RMN para esse setor particular. Pode ser observado que exceto para os maiores números de setor, os valores p são próximos de zero, significando que existe uma probabilidade muito pequena que as correlações observadas sejam um resultado do acaso.
[0055] Um processamento similar pode ser executado com as medições feitas com um gradiente de campo para fornecer marcações como uma função da difusibilidade. A figura 10(c) mostra a mesma distribuição T2 911 da figura 9(a). A figura 10(a) é uma marcação 2-D no plano T2 - difusibilidade e a figura 10(b) é uma marcação do coeficiente de correlação Pearson como uma função da distribuição T2 e difusibilidade. A linha 1001 mostra a posição esperada para uma rocha úmida com água enquanto 1003 mostra a tendência esperada para uma rocha saturada com petróleo. O intervalo de profundidade correspondendo com essa marcação era um intervalo saturado de água. Um ponto a observar é que os exemplos fornecidos mais tarde usando distribuições T2 1-D podem também ser feitos com distribuições de difusibilidade 1-D e distribuições 2-D no plano T2 - difusibilidade.
[0056] A figura 11 mostra uma marcação 2-D da correlação Pearson de outro exemplo de dados de campo nos quais existe quantidade considerável de petróleo presente como indicado pelos altos valores perto da linha de petróleo 1003. 1201 na figura 12(a) mostra a distribuição T2 RMN em média para esse intervalo, 1203 mostra o coeficiente de correlação Pearson entre o perfil de raios gamas e o perfil de setor T2 em média, 1205 mostra o coeficiente de correlação Pearson
Petição 870190094158, de 20/09/2019, pág. 24/57
20/24 entre o perfil de resistividade e o perfil de setor T2. 1207 e 1209 mostram os valores p correspondendo com 1203 e 1205, respectivamente. É mostrada na figura 12(b) a difusibilidade em média 1251, as correlações Pearson dos perfis do setor de difusibilidade com os raios gama (1253) e resistividade (1255) e os valores p correspondentes 1257, 1259.
[0057] A figura 13 mostra a distribuição T2 em média que foi codificada em cor para mostrar a correlação do perfil de setor T2 para um setor particular com o perfil de resistividade. As grandes correlações nos setores 8-10 são indicativas de uma quantidade significativa de petróleo no intervalo.
[0058] A figura 14 mostra a correlação Pearson como uma função de T2 e difusibilidade para um intervalo úmido com água. Observe os altos valores perto da linha de água 1001. 1501 na figura 15(a) mostra a distribuição T2 RMN em média para esse intervalo, 1403 mostra o coeficiente de correlação Pearson entre o perfil de raios gama e o perfil de setor T2 em média, 1505 mostra o coeficiente de correlação Pearson entre o perfil de resistividade e o perfil de setor T2. Os valores p não são marcados já que eles são próximos de zero. É mostrada na figura 15(b) a difusibilidade em média 1551, as correlações Pearson dos perfis do setor de difusibilidade com raios gama (1553) e resistividade (1555).
[0059] Com referência novamente agora à figura 11, em uma modalidade da invenção, uma filtragem 2-D é aplicada nos dados para atenuar porções dos dados perto da linha de petróleo 1003. Depois que essa filtragem é feita, a projeção dos dados filtrados sobre o eixo geométrico T2 fornecerá uma estimativa da distribuição T2 da água na formação. Alternativamente, a filtragem pode ser feita para atenuar os dados que não estão perto da linha de petróleo 1003, projetando os dados filtrados sobre o eixo geométrico T2 e subtraindo a distribuição
Petição 870190094158, de 20/09/2019, pág. 25/57
21/24 projetada da distribuição T2 original para fornecer uma estimativa da distribuição T2 da água.
[0060] A figura 16 mostra uma marcação de diversas amostras de rocha da permeabilidade em mD (ordenada) contra a porosidade (ordenada). Também são indicadas na marcação curvas de tamanhos de entrada de poro constante: para 0,1 pm (1601), 0,5 pm (1603), 2,0 pm (1605), 10 pm (1607) e 40 pm (1609). Também são mostradas no painel lateral da figura 16 as marcações da distribuição T2 para os quatro grupos de rochas de tamanhos de grão diferentes. Em cada uma das marcações no painel lateral, a curva dominante é para 100% de saturação de água enquanto a curva colorida clara é para saturação de água irredutível. As curvas na marcação são os limites das fácies (famílias ou grupos) encontradas no conjunto de amostras. Essas fácies são caracterizadas pelo tamanho da entrada do poro e têm distribuições T1 ou T2 características. No caso em que distribuições T1 ou T2 reais não são conhecidas, uma pessoa pode pegar uma distribuição T1 ou T2 características para as fácies (marcação de permeabilidade vs. porosidade) de uma tabela de consulta de fácies-RMN.
[0061] Para tamanhos de grão muito bem sortidos (não necessariamente o caso para os dados da figura 16), a distribuição T2 das fácies petrográficas (com tamanho de grão como parâmetro) pode ser simulada como mostrado na figura 17. Nesse caso da simulação usando os métodos de Georgi é basicamente unimodal, onde a posição do valor médio logaritmo de T2 pode ser correlacionada com a permeabilidade, seguindo, por exemplo, o modelo permanente SDR. A distribuição T2 para três tamanhos de grão diferentes para 100% de saturação de água é mostrada por 1701, 1703, 1705. Os números do setor correspondem com valores T2 variando de maneira logarítmica de 0,25 ms a 16381,37 ms. Rochas reais podem incluir uma combinação de distribuições unimodais, causando uma distribuição multimodal.
Petição 870190094158, de 20/09/2019, pág. 26/57
22/24 [0062] Na figura 18, a simulação de tipos de petróleo diferentes na ausência de qualquer gradiente de campo (isso significa que é igual ao T2 intrínseco - sem efeito de difusão. Quando aumentando a gravidade API do petróleo de aproximadamente 12 (petróleo pesado) até 40 API (leve) a figura 18 mostra que, em geral, quanto mais leve o petróleo e, portanto menor a sua viscosidade, a distribuição T2 ou T1 se desloca para a direita nos eixos geométricos x (MS ou # do setor). A distribuição Ti ou T2 do petróleo pode ser obtida usando mapas de difusibilidade-T2 (mapas DT2) como discutido acima. O gás, que tem viscosidade muito baixa, se situaria na extremidade direita da escala.
[0063] Como a localização da distribuição ao longo dos eixos geométricos de T2 para a amostra 100% saturada de água, não somente corresponde com o tamanho do grão e, portanto fácies, mas também com a permeabilidade (equação SDR Ksdr= a’çmi (T2lm) com a: constante; m, n: parâmetros de expoentes, φ: porosidade e T2lm: média de perfil T2 da distribuição) e a localização da distribuição para amostras de petróleo ao longo dos mesmos eixos geométricos de T2 corresponde com a gravidade API do petróleo - e muito melhor com a viscosidade -, é possível gerar um mapa de mobilidade (permeabilidade sobre viscosidade). Até algumas variáveis geométricas, a mobilidade mostra a proporcionalidade entre taxa de fluxo e pressão diferencial (lei de Darcy): Q=(k/p)*Ap:
onde Q é a taxa de fluxo.
K: permeabilidade.
μ: viscosidade.
Ap: pressão diferencial.
[0064] Como, em geral, o inverso da viscosidade do petróleo é proporcional a T20J, a mobilidade (k/μ) é, portanto, proporcional a k*T2petróleo. Sendo k proporcional a T2água, foi obtido o mapa de mobilidade multiplicando T2água*T2petróleo. O mapa é mostrado na figura 19.
Petição 870190094158, de 20/09/2019, pág. 27/57
23/24 [0065] A abscissa dessa marcação de mobilidade é o indicador de viscosidade derivado de Ti ou T2 como discutido acima. A ordenada é a distribuição Ti ou T2 da água em um meio poroso. 1901, 1903, 1905 e 1907 se referem às distribuições para a microporosidade, mesoporosidade, macroporosidade e megaporosidade, respectivamente. Conforme a figura 16, 1911, 1913, 1915 e 1917 se referem a petróleo extrapesado, petróleo pesado, petróleo médio e petróleo leve, respectivamente.
[0066] Em geral, quando para um dado reservatório ou camada de transporte de petróleo, a marcação no mapa cai no canto direito superior, a mobilidade do fluido é alta 1921 (o melhor caso em termos de produção); entretanto, se a marcação cai no canto esquerdo inferior, 1923, a situação oposta ocorre. Nas situações nas quais uma classificação restrita da qualidade da rocha (fácies) e viscosidade do petróleo (ou gravidade API) em termos das janelas de Ti ou T2, não se firma em cada situação, os mapas têm uma validade universal já que ele representa a lei de Darcy e pode ser usado para mostrar as tendências na avaliação de um reservatório. Quando as marcações têm tendência a ficar acima da diagonal com inclinação positiva, pode ser entendido que a qualidade da rocha excede a qualidade do fluido.
[0067] O método da presente revelação é descrito acima com referência a uma ferramenta de perfilagem de RMN transportada por cabo de aço. O método pode também ser usado em ferramentas de perfilagem transportadas em tubulação espiralada em poços quase horizontais. O método pode também ser usado em sensores de RMN transportados em um tubular de perfuração, tal como uma coluna de perfuração ou tubulação espiralada para aplicações de medição sem interromper a perfuração (MWD). Como é prática padrão na perfilagem de poço, os resultados do processamento são gravados em um meio adequado. Está implícito no processamento dos dados o uso de um
Petição 870190094158, de 20/09/2019, pág. 28/57
24/24 programa de computador implementado em um meio legível por máquina adequado que possibilita que o processador execute o controle e o processamento. O meio legível por máquina pode incluir ROMs, EPROMs. EAROMs, memórias flash e discos óticos.
[0068] Embora a descrição precedente seja direcionada para as modalidades específicas da invenção, várias modificações serão evidentes para aqueles versados na técnica. É planejado que todas as variações dentro do escopo das reivindicações anexas sejam abrangidas pela descrição precedente.

Claims (20)

  1. REIVINDICAÇÕES
    1. Método de avaliação de uma formação terrestre, o método compreendendo:
    (a) transportar uma ferramenta de perfilagem para dentro de um poço (10), (b) obter sinais de ressonância magnética nuclear (RMN) em uma pluralidade de profundidades no poço (10), caracterizado por (c) processar os sinais de RMN para obter uma distribuição do tempo de relaxamento em cada uma da pluralidade de profundidades, cada uma das distribuições compreendendo uma pluralidade de setores, (d) produzir, a partir da pluralidade de distribuições do tempo de relaxamento, uma pluralidade de perfis de setor de RMN dos valores do tempo de relaxamento para cada um da pluralidade de setores, (e) determinar uma medida da similaridade de um primeiro perfil indicativo do conteúdo de hidrocarboneto na formação com cada um da pluralidade de perfis de setor de RMN e (f) usar a medida de similaridade para cada um da pluralidade de perfis de setor de RMN para identificar um subconjunto da pluralidade de perfis de setor de RMN indicativo de hidrocarboneto.
  2. 2. Método, de acordo com a reivindicação 1, caracterizado pelo fato de que os sinais de RMN compreendem sinais de eco por rotação e o tempo de relaxamento compreende um tempo de relaxamento transversal T2.
  3. 3. Método, de acordo com a reivindicação 1, caracterizado pelo fato de que a medida de similaridade é pelo menos um de: (i) um coeficiente de correlação Pearson, (ii) um valor p caracterizando uma probabilidade de um coeficiente de correlação observado e (iii) uma
    Petição 870190094158, de 20/09/2019, pág. 30/57
    2/6 informação mútua.
  4. 4. Método, de acordo com a reivindicação 1, caracterizado pelo fato de que o primeiro perfil é pelo menos um de: (i) um perfil de raios gama e (ii) um perfil de resistividade.
  5. 5. Método, de acordo com a reivindicação 1, caracterizado pelo fato de que determinar a medida de similaridade ainda compreende executar pelo menos um de: (i) uma correlação cruzada do primeiro perfil com pelo menos um da pluralidade de perfis de setor de RMN e (ii) a determinação de uma informação mútua entre o primeiro perfil e pelo menos um da pluralidade de perfis de setor de RMN.
  6. 6. Método, de acordo com a reivindicação 2, caracterizado pelo fato de que ainda compreende:
    (i) obter sinais adicionais de eco por rotação em cada uma da pluralidade de profundidades enquanto aplicando um gradiente de campo magnético externo, (ii) processar os sinais de eco por rotação e os sinais adicionais de eco por rotação para obter uma distribuição de um coeficiente de difusão (D) em cada uma da pluralidade de profundidades, cada uma da distribuição dos coeficientes de difusão compreendendo uma pluralidade de setores de difusão, (iii) produzir, a partir da pluralidade de distribuições do coeficiente de difusão, uma pluralidade de perfis de setor de coeficiente de difusão de valores de setor de coeficiente de difusão correspondendo com cada um da pluralidade dos setores de difusão, (iv) determinar uma medida da similaridade do primeiro perfil com cada um da pluralidade de perfis de setor do coeficiente de difusão e (v) usar a medida de similaridade para cada uma da pluralidade de perfis de setor de coeficiente de difusão para identificar um subconjunto dos perfis de setor do coeficiente de difusão indicativo do
    Petição 870190094158, de 20/09/2019, pág. 31/57
    3/6 hidrocarboneto.
  7. 7. Método, de acordo com a reivindicação 2, caracterizado pelo fato de que ainda compreende usar o subconjunto identificado dos perfis e uma distribuição T2 em pelo menos uma da pluralidade de profundidades para obter uma distribuição T2 modificada responsiva primariamente a pelo menos um de: (i) água e (ii) petróleo na formação.
  8. 8. Método, de acordo com a reivindicação 7, caracterizado pelo fato de que ainda compreende estimar da distribuição modificada pelo menos um de: (i) saturação da água, (ii) água retida na argila, (iii) água retida irredutível, (iv) água móvel, (v) permeabilidade e (vi) hidrocarboneto.
  9. 9. Método, de acordo com a reivindicação 1, caracterizado pelo fato de que ainda compreende transportar a ferramenta de perfilagem para dentro do poço (10) em um dispositivo de transporte selecionado de (i) um cabo de aço (16) e (ii) um tubular de perfuração.
  10. 10. Aparelho para avaliar uma formação terrestre, o aparelho compreendendo:
    uma ferramenta de perfilagem configurada para ser transportada para dentro de um poço (10), a ferramenta de perfilagem ainda configurada para obter sinais de ressonância magnética nuclear (RMN) em uma pluralidade de profundidades no poço (10) e caracterizado por pelo menos um processador configurado para:
    (A) processar os sinais de RMN para obter uma distribuição do tempo de relaxamento em cada uma da pluralidade de profundidades, cada uma das distribuições compreendendo uma pluralidade de setores, (B) usar a pluralidade de distribuições do tempo de relaxamento para produzir uma pluralidade de perfis de setor de RMN dos
    Petição 870190094158, de 20/09/2019, pág. 32/57
    4/6 valores do tempo de relaxamento correspondendo com cada um da pluralidade de setores, (C) determinar uma medida da similaridade de um primeiro perfil indicativo do conteúdo de hidrocarboneto na formação com cada um da pluralidade de perfis de setor de RMN e (D) usar a medida de similaridade para cada um da pluralidade de perfis de setor de RMN para identificar um subconjunto dos perfis de setor de RMN indicativo de hidrocarboneto.
  11. 11. Aparelho, de acordo com a reivindicação 10, caracterizado pelo fato de que os sinais de RMN compreendem sinais de eco por rotação e o tempo de relaxamento compreende um tempo de relaxamento transversal (T2).
  12. 12. Aparelho, de acordo com a reivindicação 10, caracterizado pelo fato de que pelo menos um processador é ainda configurado para usar, como a medida de similaridade, pelo menos um de: (i) um coeficiente de correlação Pearson, (ii) um valor p caracterizando uma probabilidade de um coeficiente de correlação observado e (iii) uma informação mútua.
  13. 13. Aparelho, de acordo com a reivindicação 10, caracterizado pelo fato de que ainda compreende pelo menos uma de uma ferramenta de perfilagem de raios gama e uma ferramenta de perfilagem de resistividade e sendo que o primeiro perfil é pelo menos um de: (i) um perfil de raios gama e (ii) um perfil de resistividade.
  14. 14. Aparelho, de acordo com a reivindicação 10, caracterizado pelo fato de que pelo menos um processador é configurado para determinar a medida de similaridade por ainda executar pelo menos um de: (i) uma correlação cruzada do primeiro perfil com pelo menos um da pluralidade de perfis de setor de RMN e (ii) a determinação de uma informação mútua entre o primeiro perfil e pelo menos um da pluralidade de perfis de setor de RMN.
    Petição 870190094158, de 20/09/2019, pág. 33/57
    5/6
  15. 15. Aparelho, de acordo com a reivindicação 11, caracterizado pelo fato de que a ferramenta de perfilagem é ainda configurada para obter sinais adicionais de eco por rotação em cada uma da pluralidade de profundidades enquanto aplicando um gradiente de campo magnético externo e sendo que o processador é ainda configurado para:
    (i) processar os sinais de eco por rotação e os sinais adicionais de eco por rotação para obter uma distribuição de um coeficiente de difusão em cada uma da pluralidade de profundidades, cada uma da distribuição dos coeficientes de difusão compreendendo uma pluralidade de setores de difusão, (ii) produzir, a partir da pluralidade de distribuições do coeficiente de difusão, uma pluralidade de perfis de setor de coeficiente de difusão de valores de setor de coeficiente de difusão correspondendo com cada um da pluralidade dos setores de difusão, (iii) determinar uma medida da similaridade do primeiro perfil com cada um da pluralidade de perfis de setor do coeficiente de difusão e (iv) usar a medida de similaridade para cada uma da pluralidade de perfis de setor de coeficiente de difusão para identificar um subconjunto dos perfis de setor do coeficiente de difusão indicativo de pelo menos um de: (i) um hidrocarboneto e (ii) água.
  16. 16. Aparelho, de acordo com a reivindicação 11, caracterizado pelo fato de que o processador é ainda configurado para usar o subconjunto identificado dos perfis e uma distribuição T2 em pelo menos uma da pluralidade de profundidades para obter uma distribuição T2 modificada responsiva primariamente à água na formação.
  17. 17. Aparelho, de acordo com a reivindicação 16, caracterizado pelo fato de que o processador é ainda configurado para estimar da distribuição modificada pelo menos um de: (i) saturação da
    Petição 870190094158, de 20/09/2019, pág. 34/57
    6/6 água, (ii) água retida na argila, (iii) água retida irredutível, (iv) água móvel, (v) permeabilidade e (vi) hidrocarboneto.
  18. 18. Aparelho, de acordo com a reivindicação 10, caracterizado pelo fato de que ainda compreende um dispositivo de transporte configurado para transportar a ferramenta de perfilagem para dentro do poço (10), o dispositivo de transporte selecionado de (i) um cabo de aço (16) e (ii) um tubular de perfuração.
  19. 19. Produto de meio legível por computador caracterizado por possuir instruções armazenadas em si que, quando lidas por pelo menos um processador, fazem com que o pelo menos um processador:
    processe os sinais de eco por rotação obtidos em uma pluralidade de profundidades no poço (10) na formação terrestre para estimar uma distribuição de um tempo de relaxamento transversal (T2) em cada uma da pluralidade de profundidades, cada uma das distribuições compreendendo uma pluralidade de setores, use a pluralidade de distribuições do tempo de relaxamento para produzir uma pluralidade de perfis de setor de RMN dos valores do tempo de relaxamento correspondendo com cada um da pluralidade de setores, determinar uma medida da similaridade de um primeiro perfil indicativo do conteúdo de hidrocarboneto na formação com cada um da pluralidade de perfis de setor de RMN, e usar a medida de similaridade para cada um da pluralidade de perfis de setor de RMN para identificar um subconjunto dos perfis de setor de RMN indicativos de hidrocarboneto.
  20. 20. Produto, de acordo com a reivindicação 19, caracterizado pelo fato de que ainda compreende pelo menos um de: (i) uma ROM, (ii) uma EPROM, (iii) uma EAROM, (iv) uma memória flash e (v) um disco ótico.
BRPI0907021A 2008-01-30 2009-01-30 método de avaliação de uma formação terrestre, aparelho para avaliar uma formação terrestre e produto de meio legível por computador BRPI0907021B1 (pt)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/022,755 US7804297B2 (en) 2008-01-30 2008-01-30 Methodology for interpretation and analysis of NMR distributions
PCT/US2009/032606 WO2009097507A2 (en) 2008-01-30 2009-01-30 Methodology for interpretation and analysis of nmr distributions

Publications (2)

Publication Number Publication Date
BRPI0907021A2 BRPI0907021A2 (pt) 2015-07-07
BRPI0907021B1 true BRPI0907021B1 (pt) 2020-02-04

Family

ID=40898562

Family Applications (1)

Application Number Title Priority Date Filing Date
BRPI0907021A BRPI0907021B1 (pt) 2008-01-30 2009-01-30 método de avaliação de uma formação terrestre, aparelho para avaliar uma formação terrestre e produto de meio legível por computador

Country Status (5)

Country Link
US (1) US7804297B2 (pt)
BR (1) BRPI0907021B1 (pt)
GB (1) GB2468824B (pt)
NO (1) NO342538B1 (pt)
WO (1) WO2009097507A2 (pt)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008010810A1 (en) * 2006-07-21 2008-01-24 Halliburton Energy Services, Inc. Fluid saturation estimation
US8330460B2 (en) 2008-01-30 2012-12-11 Baker Hughes Incorporated Method and apparatus for determining multiscale similarity between NMR measurements and a reference well log
US7917294B2 (en) * 2008-03-26 2011-03-29 Baker Hughes Incorporated Determination of irreducible water cut-off using two dimensional nuclear magnetic resonance data
US7893692B2 (en) * 2008-11-03 2011-02-22 Schlumberger Technology Corporation Method for estimating the formation productivity from nuclear magnetic resonance measurements
US8970217B1 (en) 2010-04-14 2015-03-03 Hypres, Inc. System and method for noise reduction in magnetic resonance imaging
US8781745B2 (en) 2010-08-03 2014-07-15 Baker Hughes Incorporated NMR-DNA fingerprint
WO2013066549A1 (en) 2011-10-31 2013-05-10 Baker Hughes Incorporated Hydrocarbon determination in unconventional shale
US20130257424A1 (en) * 2012-03-27 2013-10-03 Schlumberger Technology Corporation Magnetic resonance rock analysis
US9678185B2 (en) 2013-03-15 2017-06-13 Pepsico, Inc. Method and apparatus for measuring physico-chemical properties using a nuclear magnetic resonance spectrometer
WO2015021088A1 (en) * 2013-08-06 2015-02-12 Schlumberger Canada Limited Methods for determining a saturation-height function in oil and gas reservoirs
MX2016006493A (es) 2013-12-19 2017-01-18 Halliburton Energy Services Inc Clasificacion de tamaño de poro en formaciones subterraneas segun distribuciones de relajacion de resonancia magnetica nuclear (rmn).
US10711602B2 (en) * 2015-07-22 2020-07-14 Halliburton Energy Services, Inc. Electromagnetic monitoring with formation-matched resonant induction sensors
US10156531B2 (en) * 2015-10-27 2018-12-18 Saudi Arabian Oil Company Measuring connectivity between different pore types in porous media
US10739489B2 (en) * 2016-01-15 2020-08-11 Baker Hughes, A Ge Company, Llc Low gradient magnetic resonance logging for measurement of light hydrocarbon reservoirs
US10209391B2 (en) 2016-08-23 2019-02-19 Baker Hughes, A Ge Company, Llc Simultaneous inversion of NMR multiple echo trains and conventional logs
US10557962B2 (en) 2016-09-16 2020-02-11 Saudi Arabian Oil Company Method for measurement of hydrocarbon content of tight gas reservoirs
CN110019119B (zh) * 2017-10-12 2021-09-21 中石化石油工程技术服务有限公司 一种储集层含油气性录井解释方法及装置
US11294093B2 (en) * 2018-03-01 2022-04-05 Schlumberger Technology Corporation Nuclear magnetic resonance data acquisition system
US11754746B2 (en) * 2020-02-21 2023-09-12 Saudi Arabian Oil Company Systems and methods for creating 4D guided history matched models
US20230068555A1 (en) * 2021-08-31 2023-03-02 Halliburton Energy Services, Inc. Correction of distorted gradient distributions in nuclear magnetic resonance logging
US12578500B2 (en) 2022-03-10 2026-03-17 Saudi Arabian Oil Company Hydrocarbon reservoir saturation logging

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6392409B1 (en) * 2000-01-14 2002-05-21 Baker Hughes Incorporated Determination of T1 relaxation time from multiple wait time NMR logs acquired in the same or different logging passes
US5680043A (en) * 1995-03-23 1997-10-21 Schlumberger Technology Corporation Nuclear magnetic resonance technique for determining gas effect with borehole logging tools
US6956371B2 (en) * 1995-10-12 2005-10-18 Halliburton Energy Services, Inc. Method and apparatus for detecting diffusion sensitive phases with estimation of residual error in NMR logs
US6512371B2 (en) 1995-10-12 2003-01-28 Halliburton Energy Services, Inc. System and method for determining oil, water and gas saturations for low-field gradient NMR logging tools
US5698979A (en) 1996-02-23 1997-12-16 Western Atlas International, Inc. Method for NMR diffusion measurement
US6147489A (en) * 1997-04-09 2000-11-14 Schlumberger Technology Corporation Method and apparatus for measuring total nuclear magnetic resonance porosity
US6429654B1 (en) 1998-09-11 2002-08-06 Baker Hughes Incorporated Nuclear magnetic resonance pulse sequence for improving signal-to-noise ratio
NO324200B1 (no) 1999-04-19 2007-09-10 Baker Hughes Inc Kjernemagnetiske resonansmalinger i bronnlogging som benytter en optimalisert refaseregulerende pulssekvens
US6331775B1 (en) * 1999-09-15 2001-12-18 Baker Hughes Incorporated Gas zone evaluation by combining dual wait time NMR data with density data
US6337568B1 (en) * 1999-10-25 2002-01-08 Tarek A. Tutunji System and method for enhanced vertical resolution magnetic resonance imaging logs
US6255819B1 (en) * 1999-10-25 2001-07-03 Halliburton Energy Services, Inc. System and method for geologically-enhanced magnetic resonance imaging logs
US6600315B1 (en) * 2000-03-03 2003-07-29 Schlumberger Technology Corporation Method for improving resolution of nuclear magnetic resonance measurements by combining low resolution high accuracy measurements with high resolution low accuracy measurements
EP1301810B1 (en) 2000-07-21 2008-09-10 Services Petroliers Schlumberger Nuclear magnetic resonance methods for extracting information about a fluid in a rock
US6348792B1 (en) 2000-07-27 2002-02-19 Baker Hughes Incorporated Side-looking NMR probe for oil well logging
US6528995B1 (en) * 2001-09-10 2003-03-04 Schlumberger Technology Corporation Methods and apparatus for measuring flow velocity in a wellbore using NMR and applications using same
US6972564B2 (en) 2001-11-06 2005-12-06 Baker Hughes Incorporated Objective oriented methods for NMR log acquisitions for estimating earth formation and fluid properties
US6781371B2 (en) * 2002-09-06 2004-08-24 Schlumberger Technology Corporation High vertical resolution antennas for NMR logging
US6808028B2 (en) * 2002-12-03 2004-10-26 Schlumberger Technology Corporation Method and apparatus utilizing NMR measurements to gather information on a property of the earth formation surrounding a wellbore
WO2004090557A2 (en) * 2003-04-01 2004-10-21 Halliburton Energy Services, Inc. Abnormal pressure determination using nuclear magnetic resonance logging
US7049815B2 (en) 2004-01-27 2006-05-23 Baker Hughes Incorporated Method and apparatus for multi-frequency NMR diffusion measurements in the presence of internal magnetic field gradients
US7309983B2 (en) * 2004-04-30 2007-12-18 Schlumberger Technology Corporation Method for determining characteristics of earth formations
US7196516B2 (en) * 2004-08-16 2007-03-27 Baker Hughes Incorporated Correction of NMR artifacts due to constant-velocity axial motion and spin-lattice relaxation
US7511487B2 (en) * 2007-02-27 2009-03-31 Schlumberger Technology Corporation Logging method for determining characteristic of fluid in a downhole measurement region
US7502691B2 (en) * 2007-07-31 2009-03-10 Baker Hughes Incorporated Method and computer program product for determining a degree of similarity between well log data

Also Published As

Publication number Publication date
US7804297B2 (en) 2010-09-28
GB201012802D0 (en) 2010-09-15
GB2468824A (en) 2010-09-22
WO2009097507A3 (en) 2009-11-26
NO342538B1 (no) 2018-06-11
WO2009097507A2 (en) 2009-08-06
BRPI0907021A2 (pt) 2015-07-07
US20090189604A1 (en) 2009-07-30
GB2468824B (en) 2012-02-22
NO20101090L (no) 2010-08-27

Similar Documents

Publication Publication Date Title
BRPI0907021B1 (pt) método de avaliação de uma formação terrestre, aparelho para avaliar uma formação terrestre e produto de meio legível por computador
US8278922B2 (en) Continuous wettability logging based on NMR measurements
Gao et al. Determination of movable fluid percentage and movable fluid porosity in ultra-low permeability sandstone using nuclear magnetic resonance (NMR) technique
US8330460B2 (en) Method and apparatus for determining multiscale similarity between NMR measurements and a reference well log
AU653632B2 (en) Permeability measure from NMR relaxation measurements for fluids in porous media
Xiao et al. Effect of hydrocarbon on evaluating formation pore structure using nuclear magnetic resonance (NMR) logging
US20130057277A1 (en) Method for Determining Rock Formation Fluid Interaction Properties Using Nuclear Magnetic Resonance Well Logging Measurements
BRPI0410179B1 (pt) método para registro de uma formação de terra e aparelho para utilização em uma formação de terra
Mitchell et al. Mapping oil saturation distribution in a limestone plug with low-field magnetic resonance
BRPI0809121A2 (pt) Determinação de perfil radial de saturação de gás a partir de dados de rmn de múltiplas frequências
US8781745B2 (en) NMR-DNA fingerprint
EP3403078B1 (en) Low gradient nuclear magnetic resonance logging for measurement of light hydrocarbon reservoirs
Zielinski et al. Restricted diffusion effects in saturation estimates from 2D diffusion-relaxation NMR maps
US8532929B2 (en) Method and apparatus to incorporate internal gradient and restricted diffusion in NMR inversion
Hürlimann et al. NMR well logging
AlGhamdi et al. Correlations between NMR-relaxation response and relative permeability from tomographic reservoir-rock images
Tandon et al. A new NMR-based method for wettability assessment in mixed-wet rocks
Sun et al. Two-dimensional NMR logging and field test results
Coman et al. Improved NMR Logging Approach to Simultaneously Determine Porosity, T 2 and T 1
Mitchell et al. Quantitative remaining oil interpretation using magnetic resonance: from the laboratory to the pilot
Nabawy New approaches in reservoir characterization utilizing conventional and special core analyses: a comprehensive review
Liang et al. A potential NMR-based wettability index using free induction decay for rocks
Meridji et al. Fluid identification in complex clastic reservoirs using 2D NMR maps: a case study from Saudi Arabia
Althaus et al. Determine Oil and Water Saturations in Preserved Source Rocks from 2D T 1-T 2 NMR
US6894493B2 (en) Method and apparatus for NMR measurement of magnetic materials

Legal Events

Date Code Title Description
B06T Formal requirements before examination [chapter 6.20 patent gazette]
B06F Objections, documents and/or translations needed after an examination request according [chapter 6.6 patent gazette]
B06A Patent application procedure suspended [chapter 6.1 patent gazette]
B09A Decision: intention to grant [chapter 9.1 patent gazette]
B16A Patent or certificate of addition of invention granted [chapter 16.1 patent gazette]

Free format text: PRAZO DE VALIDADE: 10 (DEZ) ANOS CONTADOS A PARTIR DE 04/02/2020, OBSERVADAS AS CONDICOES LEGAIS.