BRPI0609938A2

BRPI0609938A2 - método para estimar um parámetro de fluido de formação de fundo de poço, aparelho para utilização em fundo de poço e meio legìvel por computador

Info

Publication number: BRPI0609938A2
Application number: BRPI0609938-6A
Authority: BR
Inventors: Rocco Difoggio; Bernardo Pohl
Original assignee: Baker Hughes Inc
Priority date: 2005-04-22
Filing date: 2006-04-21
Publication date: 2010-05-11
Also published as: US20060236758A1; CN101223529A; US20060241866A1; US7299136B2; EA200702234A1; EA014302B1; CN101223529B

Abstract

MéTODO PARA ESTIMAR UM PARáMETRO DE FLUIDO DE FORMAçãO DE FUNDO DE POçO, APARELHO PARA UTILIZAçãO EM FUNDO DE POçO E MEIO LEGìVEL POR COMPUTADOR. A presente invenção refere-se a um sistema e a um método para estimar uma característica de um fluido de formação, que inclui a pureza final, em tempo real. O sistema proporciona a realização de uma pluralidade de medições de um parâmetro do fluido de formação extraído ao longo do tempo e ajustar uma curva não-assintática às medições para se estimar as características do fluido de formação.

Description

Relatório Descritivo da Patente de Invenção para "MÉTODO EAPARELHO PARA ESTIMAR UMA CARACTERÍSTICA DE FLUIDO".

ANTECEDENTE DA INVENÇÃO

CAMPO DA INVENÇÃO

A invenção refere-se, em geral, a um método e aparelho para aquantificação da contaminação de fluido como uma indicação de limpeza deamostra em tempo real em um ambiente do furo do poço acabado. De ma-neira específica, a invenção é um método e aparelho para a medição depropriedades físicas do fluido sendo bombeado a partir da formação quecerca um furo de poço acabado através de um cabo de aço ou monitoramen-to enquanto a broca estima a limpeza de amostra ou prediz o tempo em quepode-se obter uma amostra dotada de uma pureza desejada.

SUMÁRIO DA TÉCNICA RELACIONADA

Na exploração do furo do poço acabado, tipicamente, usam-selamas de perfuração, tais como lamas baseadas em óleo, lamas baseadasem material sintético ou lamas baseadas em água. Os filtrados provenientesdessas lamas, em geral, invadem a formação através da parede do poço atéuma extensão, o que significa que esse filtrado deve ser removido a partir daformação, a fim de acessar os fluidos de formação. A amostragem do poçoaberto é um modo eficaz de adquirir fluidos de reservatório representativos.

A aquisição da amostra permite a determinação de informações críticas paraacessar o valor econômico das reservas. Em adição, pode-se designar asestratégias de produção mais favoráveis para lidar com esses fluidos com-plexos. Na amostragem do poço aberto, inicialmente, o fluxo da formaçãocontém filtrado considerável, porém na medida em que o filtrado é drenadoda formação, o fluxo torna-se, de maneira crescente, mais rico no fluido deformação. Ou seja, o fluxo de amostra proveniente da formação contém umaalta porcentagem de fluido de formação, na medida em que o bombeamentocontinua.

Sabe-se bem que o fluido que é bombeado a partir de um furode poço acabado submete-se a um processo de limpeza, no qual a purezada amostragem aumenta com o tempo, na medida em que remove-se, gra-dualmente, o filtrado da formação e aparece menos filtrado na amostra. Nopresente, define-se fp como a fração de pureza e fc como a fração de conta-minação, onde fp + fc = 1. Na medida em que o fluido de formação amostra-do muda, o mesmo acontece com as propriedades óptica e física do fluidoamostrado, tal como absorção óptica, fluorescência, índice refrativo, viscosi-dade, densidade, velocidade do som, e módulo à granel. Uma quantidade demedições diferentes são usadas para determinar várias propriedades ópticase físicas de um fluido do interior do poço em tempo real. Ao medir essaspropriedades do fluido, portanto, proporciona a penetração qualitativa emuma pureza da amostragem de fluido, porém não proporciona um valorquantitativo, fp, para a amostra de fluido. Mesmo após bombear o fluido porum longo tempo, a fração da contaminação de fluido não cai, necessaria-mente, a zero. Em muitos casos onde, após o bombeamento por um longoperíodo de tempo, uma propriedade óptica ou física pára de mudar, substan-cialmente, a fração de contaminação (conforme determinada de maneirasubseqüente em um laboratório de superfície pode estar longe de zero e emalguns casos pode ser tão alta quanto 45%. Neste caso, a pureza final podeser da ordem de 55%.

Nos longos momentos de bombeamento, pode-se alcançar umequilíbrio dinâmico, no qual uma amostra de fluido sendo retirada de umazona rosqueada se limpa na mesma taxa que é recontaminada a partir dezonas acima e abaixo da zona rosqueada. Assim, apesar de uma proprieda-de medida no interior do poço (tal como densidade óptica ou "OD") ter para-do, substancialmente, as mudanças, a amostra pode ainda não ter 100% depureza. Esse equilíbrio dinâmico depende de diversos fatores, tais como arazão da permeabilidade da vertical pela horizontal. Portanto, para os propó-sitos desse invento, o ftp é definido com a fração da pureza final, onde a pu-reza final é a pureza alcançada em muitas vezes de longo bombeamento eque é, normalmente, menor que 100%. Assim, pode ser estimar através dasmudanças de monitoramento em OD ou alguma outra propriedade pelo tem-po (ou pelo volume bombeado) é a fração da pureza final, ftp, mas não a fra-ção da pureza do fluido de formação, fp.No momento da extração de fluido a partir de uma formação,deseja-se quantificar o progresso da limpeza, ou seja, o grau de contamina-ção de filtrado ou pureza em tempo real. Sabe-se que há muita contamina-ção de filtrado na amostra (mais que cerca de 5% ou 10%), então não hárazão para coletar a amostra de fluido de formação em um tanque de amos-tra até que o nível de contaminação caia para um nível aceitável. Por outrolado, se através do bombeamento por um longo período, é possível alcan-çar, ligeiramente, apenas o nível mais baixo de contaminação, então podenão haver qualquer necessidade para continuar bombeando. Assim, há umanecessidade de determinar o quanto será preciso bombear para obter umaamostra a partir da formação.

Quando inicia-se o bombeamento, o fluido que é bombeado con-tém uma grande quantidade de contaminação de filtrado de lama e então aporcentagem de filtrado começa a diminuir em uma taxa rápida. Esse pro-cesso de redução de contaminação de filtrado de fluido é referido como pro-cesso de limpeza de amostra. Mais adiante, o fluido bombeado contém me-nos contaminação e a porcentagem de filtrado de fluido cai em uma taxamais lenta. Mullins, et. al. publicaram trabalhos baseados no ajuste de curvade uma aderência da amostra versus tempo para monitorar a limpeza emtempo real, intitulado "Real Time Determination of Filtrate Contamination Du-ring Openhole Wireline Sampling by Optical Spectroscopy," SPWLA, 419 En-contro Anual, Dallas, Texas, Junho de 2000. As patentes U.S 6.274.865 e6.350.986 discutem determinados ajustes de curva.

Neste trabalho, Mullins et al. assume que a taxa da limpeza deamostra, conforme medida através da observação da densidade óptica pro-gride como f/12 onde ré tempo. Essa taxa de limpeza é baseada em experi-ência empírica no Golfo do México e em outros lugares. No entanto, Mullinset al. ainda confirmaram que, para as durações de bombeamento estendido,a taxa de limpeza de amostra para invasão superficial progride como f1/3 eque a taxa de limpeza para invasões mais profundas, com P. Uma suposi-ção de uma taxa de limpeza de amostra de f/12 pode ser rígida e inaplicávelem situações em tempo real. Além do mais, usando o tempo como um pa-râmetro de ajuste, necessariamente, assume uma taxa de bombeamentoconstante. Um outro problema com o monitoramento de limpeza de amostrapelo envelhecimento ao observar a absorção óptica pelo envelhecimento éque partículas de areia e outras partículas podem ocasionar uma dispersãoconsiderável, o que faz com que os valores de absorção mensurados pelotempo "pulem" e pareçam ruidosos. Assim, há uma necessidade para umsistema e método mais flexíveis para obter a estimativa de limpeza de for-mação baseada nas propriedades e características de fluido para o bombe-amento no interior do poço em tempo real.

SUMÁRIO DA INVENÇÃO

A presente invenção proporciona um método e aparelho paraquantificar a limpeza de amostra em tempo real a partir de dados mensura-dos pelo tempo (ou pelo volume) de algumas propriedades ópticas ou físicasde amostras de fluido tomadas a partir da formação que circunda um furo. Ofluido de amostra é extraído a partir da formação que circunda o furo. Namedida em que o fluido continua-se extraindo a partir da formação, a com-posição de fluido de amostra muda, alterando os valores mensurados deuma propriedade óptica ou física para o fluido de amostra.

Em um primeiro aspecto da presente invenção, proporcionam-seum método e aparelho que ajustam os dados de medição de fluido em umacurva não-assintótica. Um exemplo de uma curva não-assintótica é uma cur-va (por exemplo, uma aproximação de séries de força), que proporciona umajuste aperfeiçoado para os dados sob o tempo de bombeamento típico e,que pode, ainda, ser extrapolado várias vezes ao ponto do tempo de bom-beamento, porém que aborda mais ou menos infinidade em tempos infinitos.Um outro exemplo de uma curva não-assintótica é uma equação que possuium componente oscilatório, tal como uma onda senoidal, a qual nunca al-cança um limite fixo. A onda senoidal pode ser ajustada em freqüência, fasee amplitude a fim de proporcionar um ajuste melhorado. Em um terceiro as-pecto da presente invenção, proporcionam-se um método e aparelho querealizam reconhecimento padrão de uma linha reta para um melhor ajustedos dados mensurados em um espaço log-log.Para um melhor desempenho, removem-se primeiro as cavilhasnos dados. Os dados restantes ficam lisos como uma peça através dos cilin-dros de laminação de 100 ou mais pontos adjacentes que usam uma funçãode alisamento. Por exemplo, um ajuste pode ser realizado para a absorçãopelo segmento de tempo de laminação usando uma equação de ajuste não-assintótico, tal como A = b0 + bit + b2f. Então, pelo cálculo, A' = dA/dt = bi +2b2t e A7A = (bi + 2b2t)/( b0 + bit + b2?) pode-se determinar. Então, parauma equação da forma, A(t) = A0 - Ait'p, pode-se fazer uma regressão linearde 1 n(dA/dt) contra 1 n(t) para obter a inclinação e interceptação e a partirdestes calcular, -p = (1 + inclinação) e -A^ = exp(interceptação -1n(1+inclinação)). Dessa maneira, não há uma suposição de um valor de -5/12, de -2/3 (conforme sugeriu Mullins), ou qualquer outro valor fixo para -p. Ao invés disso, pode-se estimar ftp = A/A0 a partir dos valores de melhorajuste para p e A1; e a partir de duas vezes a média de A(t) e Àit"p em umapluralidade de vezes.

Ainda, o método e aparelho da presente invenção pode usaruma equação de ajuste por dados, tal como log(1-ftp) = (-p) log(t) +logíAi/Ao), a qual é a equação de uma linha reta que não possui assíntota (Y= constante), exceto para o caso insignificante de p=0. O método e aparelhorealizam uma série de regressões que usam estimativas diferentes de A0,porém não calculam realmente o próprio A0. Por exemplo, pode-se iniciarcom o valor corrente, A, em um tempo t, como a primeira estimativa de A0,então procede para um valor ligeiramente maior que A+z, então para um va-lor ainda maior que A+2e, e por aí vai. O valor A0 para o qual o ajuste paraos dados mensurados é o mais perto do formato de uma linha reta (baseadono coeficiente de determinação mais alto, ou valor quadrado de R) torna-se,então, a melhor estimativa de um valor A0. Em um terceiro aspecto, propor-cionam-se o método e aparelho da presente invenção que ajustam uma cur-va diferenciável aos dados de medição ou dados de propriedade física deri-vados dos dados de medição. A presente invenção então estima ftp a partirda razão de (dA/dt) para A.

Um outro aspecto da presente invenção, proporcionam-se ummétodo e aparelho que ajustam uma curva assintótica à diferença de duasrespostas, tais como a diferença de duas absorções associadas com com-primentos de ondas diferentes (canais ópticos) ao invés da própria absorção.O uso de uma diferença de absorção remove as saliências da linha baseocasionada pela passagem de partículas de areia ou bolhas.

BREVE DESCRIÇÃO DAS FIGURAS

Outros objetos e vantagens da invenção serão aparentes medi-ante a leitura da descrição detalhada que se segue e mediante a referênciados desenhos anexos, em que:

A figura 1 é um diagrama do Módulo de Caracterização de Fluido;

A figura 2 é uma ilustração de uma modalidade da presente in-venção disposta em um poço que usa uma pluralidade de fontes e sensores;

As figuras 3-10 são diagramas de funções realizadas em modalidades da invenção;

A figura 11 é uma ilustração de uma modalidade da invençãoque usa um transdutor acústico;

A figura 12 é uma ilustração de uma modalidade da invençãoque usa um conjunto piroelétrico; e

A figura 13 é uma ilustração de uma função realizada de acordocom uma outra modalidade da invenção.

BREVE DESCRIÇÃO DA INVENÇÃO

A figura 1 ilustra uma representação esquemática para um mó-dulo de caracterização de fluido do interior do poço 100 para a obtenção eanálise de dados de medição ópticos. Uma fonte de luz 101 (por exemplolâmpada de tungstênio) emite luz em direção a um fluido 110. A luz pode sercolimada por um dispositivo de lente de colimação 103 disposta entre a fontede luz 101 e um fluido 110. A luz colimada 111 é incidente, em geral perpen-dicular a uma primeira janela de safira 301 adjacente à amostra 110. As ja-nelas de safira 301 e 303 ficam, em geral, perpendiculares aos feixes de luzcolimados e ficam separadas por um espaço ou canal 304 que permite queum fluido 110 flua entre elas. Uma propriedade óptica do fluido, por exemplo,que inclui, mas não se limita ao reflexo, absorção e fluorescência da luz pro-veniente do fluido, é medida pelo tempo através de um sensor óptico, talcomo um espectrômetro 105, porém não se limita a esse. Proporciona-se umprocessador 113 para estimar as propriedades de fluido a partir das medi-ções ópticas. As ferramentas existentes (figura 1) podem ser ajustadas comuma fonte de luz ultravioleta ou infravermelha 112, que pode ser ligadaquando a fonte de luz de tungstênio 101 for desligada. O mesmo espectrô-metro, por exemplo, que compreende filtros de comprimento de ondas úni-cos sob fotodiodos, permite coletar o espectro infravermelho e fluorescênciade óleo cru. O processador 113 inclui memória e executa cálculos usandoequações para estimar características ou propriedades de fluido, tais comocontaminação percentual, a partir das medições ópticas para o fluido, con-forme descrito no presente. Proporciona-se força para diversos componen-tes do módulo 100 através de um suprimento de força.

Conforme mostrado na figura 2, as medições adicionais proveni-entes de fontes e sensores adicionais podem ser adicionadas, que inclui,mas não se limita a um ressonador mecânico flexural, transdutor acústico,conjunto piroelétrico, fonte de luz infravermelha, e sensores para medir índi-ce retroativo. Mais esquemas detalhados do transdutor acústico e do conjun-to piroelétrico são mostrados nas figuras 11 e 12. Essas fontes e sensoresadicionais podem ser proporcionados para medições de parâmetros de fluidoque inclui, porém não limita, viscosidade, densidade, velocidade do som,fluorescência, rarefação, reflexo total, índice refrativo, módulo à granel e re-sistividade. As medições podem ser monitoradas através do tempo para es-timar uma característica do fluido, inclusive pureza final fracional ou conta-minação final fracional, conforme discutido abaixo.

A figura 2 ilustra uma modalidade de um sistema disposto emum poço 12 perfurado a partir da superfície 15 e formado em uma formação16 que pode ser usada para realizar os métodos da invenção. A sonda 14 éproporcionada para a extração do fluido a partir da formação. Os sensoresde medição, tais como o sensor 100 ou outros sensores da presente inven-ção, estão contidos em uma ferramenta de interior de poço 20. A ferramentade interior de poço é disposta a partir de um cabo de aço ou uma coluna deperfuração 18. A ferramenta 20 ainda inclui um controlador que controla aoperação da ferramenta de interior de poço. O controlador inclui memória eprogramas, que inclui algoritmos descritos no presente para executar os mé-todos descritos neste. Durante a operação, a ferramenta 20 é abaixada nofuro do poço e ajustada para obter amostras de fluido a partir da formação16. O sensor 100 mede uma propriedade ou característica desejada do flui-do pelo tempo. O controlador na ferramenta 20 ou um outro controlador nasuperfície, que utiliza os programas, realiza o método descrito no presentedocumento abaixo e proporciona as estimativas desejadas e outros resulta-dos, conforme descrito neste.

Os sistemas anteriores usaram a forma funcional para a limpezaaproximada por In(OD) = C - D/f. Alguns sistemas anteriores calcularam aporcentagem de contaminação ao assumir que, mediante o alcance de ab-sorção óptica assintótica, a amostra atingiu zero em contaminação. Outrossistemas anteriores, no entanto, assumiram que um equilíbrio dinâmico podeser alcançado entre a limpeza de fluido e a incursão de filtrado continuado,que depende da razão da permeabilidade vertical pela horizontal e outrosfatores. Assim, a contaminação pode não cair a zero, porém apenas paraalgum valor mínimo, mesmo depois de períodos muito longos de bombea-mento. Ou seja, o nível de contaminação final representa a contaminaçãomínima, mas não necessariamente contaminação zero.

Na medida em que a composição das propriedades de fluido deformação amostradas muda, o mesmo acontece com as propriedades físicase ópticas do fluido amostrado, tal como absorção óptica, fluorescência, índi-ce refrativo, viscosidade, densidade, velocidade do som, e módulo à granel.Essas propriedades podem ser monitoradas para estimar a fração da purezafinal, a qual é o grau de limpeza de fluido de formação. Pode-se usar medi-ções diferentes com base nas condições reais. Por exemplo, em determina-dos casos, com o monitoramento da limpeza pelo tempo através da obser-vação na absorção óptica pelo tempo (por um trajeto com comprimento de 2mm) pode ser menos desejável porque as partículas de areia e outras parti-cuias podem ocasionar uma dispersão considerável, o que faz com que aabsorção por tempo "pule" muito e pareça bastante ruidosa. Por outro lado,o monitoramento da limpeza pelo tempo através do monitoramento do índicerefrativo (o que é uma técnica baseada em interface) é menos sensível àspartículas na corrente de fluido, porque observa-se apenas uma camada finado fluido que está em contato direto (na interface) com a janela de safira.Similarmente, para os óleos crus, a fluorescência vê apenas uma camadafina de óleo cru próxima ã janela, e dessa maneira, é bastante insensível àspartículas na corrente. No entanto, as medições de qualquer parâmetro ade-quado ou característica do fluido podem ser usadas para o propósito da in-venção.

Como um exemplo de limpeza de fluido de formação, os resulta-dos de simulação e dados de como reais para a densidade óptica podem serajustados em formas tais como, Y = mXp + b ou ln(Y) = mXp + b. Usando aconcentração de contaminação de amostra como Y e o tempo de bombea-mento como X, os resultados de simulação se ajustam bem à essas formas,especialmente a forma logarítmica. Devido à densidade óptica (OD) é umindicador de limpeza, pode-se usar o OD como Y e o tempo de bombeamen-to pode se usado como X. Caso a velocidade de bombeamento mude muitasvezes durante o curso da limpeza, o volume cumulativo bombeado é usadocomo X ao invés do tempo.

Um pequeno valor p indica que o processo de limpeza é lento elevará mais tempo para obter uma amostra de qualidade, enquanto umgrande valor p indica que o processo de limpeza será mais rápido e a chan-ce de obter uma amostra da pureza desejada é alta. O valor de b é usadocom um indicador para a limpeza da melhor qualidade de amostra alcançá-vel (o valor assintótico). Ao comparar o valor OD corrente com o valor b, ob-têm-se a porcentagem de contaminação de amostra corrente. A qualidadede amostra futura é estimada ao usar valores ajustados de m, p, e b,e pode-se tomar uma decisão de continuar ou parar o processo de bombeamento,caso a futura qualidade de amostra estimada for julgada insuficiente.

A força de X, que é -p, (onde p é um número positivo) pode serusado como um indicador da rapidez do processo de limpeza. Por conseqüência, quando o valor p for pequeno, e a contaminação corrente calculada for alta, haverá pouca chance de obter uma amostra de alta qualidade e levará muito tempo, talvez tempo demais para obter a pureza de amostra desejada.

Quando o perfil de invasão for gradual, que significa que uma zona de transição está presente ao invés do limite de invasão bem definido, o valor de p cai para abaixo de 1,0. O valor de p depende da espessura da zona de transição entre a região do filtrado e a região do fluido de formação.

Quanto mais espessa a zona de transição, mais baixo o valor p. Essa transição gradual possui um efeito similar ao da invasão profunda. Quando a invasão é profunda, então, o fluido limpo proveniente da zona fresca será misturado com o filtrado, na medida em que este escoa em direção à sonda. Portanto, uma invasão profunda terá uma zona de transição espessa, e a limpeza desta zona levará bastante tempo.

O dano de formação pode afetar o processo de limpeza. Pode-se melhorar a limpeza quando a formação próxima ao furo do poço for danificada ou quando a permeabilidade da formação do furo do poço mais próximo for menor que a verdadeira permeabilidade de formação, devido à pequena invasão de partícula.

Uma forma funcional que pode ser melhor ajustada aos dados de simulação, independente da profundidade da invasão ou dano de formação é OD = exp(mt"p+b), de modo que em períodos mais longos, a OD pára de mudar, porque o tempo dependente de tempo cai a zero, conforme o tempo é infinito. Essa forma é equivalente à forma, In(OD) = C - D/tAp, onde C = b, m = -D, e p é um número positivo. O ajuste de curva da limpeza de amostra que usa a forma de ln(Y) = mX"p + B, onde Y = propriedades ópticas ou físicas, tais como absorção ou fluorescencia em alguns comprimentos de onda como a amostra na medida em que efetua-se a limpeza, X = Tempo desde começado o bombeamento da amostra ou, mais corretamente, o volume cumulativo bombeado, quando a taxa de fluxo volumétrico não for constante.Na aplicação da forma funcional In(OD) = A * 1/Tempon + B para resultados de simulação, mostra-se que: a) em um sistema de amostra sem danos de permeabilidade e uma zona de filtragem bem definida (100% da zona de filtrado seguido de 0% de zona), o valor n é 1,1; b) adicionando dano de permeabilidade ao sistema, acelera o processo de limpeza, e o valor de n é 1,3; c) usando uma contaminação de filtrado gradual, em que a contaminação reduz conforme se movimenta para longe do furo do poço (isto é 100, 80, 60.. e 0%), então o valor de n é 0,75; d) adicionando prejuízo ao sistema (c), então o valor de n é 1,0; e) adicionando uma mudança de per-meabilidade devido aos danos de formação, então o valor de n pode variar de 0,25 a 0,5. Ao ajustar os resultados da simulação de limpeza de formação de filtrado e alguns dados de campo (densidade óptica) à forma funcional acima, proporcionam-se as descobertas que se seguem.

Ao usar a concentração de contaminação de amostra como Y e o volume de bombeamento cumulativo como X, os resultados de simulação ajustam a forma muito bem. Devido ao fato de a densidade óptica (OD) ser um indicador de limpeza, os dados de OD podem ser usados como Y e o tempo de bombeamento pode ser usado como X. Caso a velocidade de bombeamento mudar várias vezes durante o curso da limpeza, o volume cumulativo bombeado deve ser usado como X ao invés do tempo. A força de X, que é -p (onde p é um número positivo), pode ser usada como um indicador para o processo de limpeza. Para uma invasão nominal de 15,24 cm (6 polegadas) ou menos, p está perto de 1,0 - 1,1. Quando o perfil de invasão for gradual, isso indica que há uma zona de transição ao invés de um limitede invasão bem definido, e p cai para abaixo de 1,0. O valor de p vai depender da espessura da zona de transição, quanto mais espessa a zona, menor o valor de p. Um efeito similar é encontrado para a invasão profunda. Quando a invasão for profunda, então o fluido limpo proveniente da zona fresca será misturado como filtrado enquanto este flui em direção à sonda. Então, ainvasão profunda terá uma zona de transição espessa, e levará um longo tempo para limpar essa zona.

Quando a formação próxima ao furo do poço está danificada, ouseja, a permeabilidade da formação próxima ao furo do poço é menor que a verdadeira permeabilidade de formação devido à pequena invasão de partícula, a limpeza pode ser melhorada, conforme discutido nos trabalhos na Sociedade dos Engenheiros de Petróleo, trabalhos SPE 39817 e SPE 48958. A anisotropia de formação também ajuda no processo de limpeza (vide trabalhos SPE 39817 e SPE 48958). Quando há um dano ou uma anisotropia de formação, o valor de p aumenta para acima de 1,0 - 1,1 até 1,3 - 1,4. Então, um pequeno valor de p indica que o processo de limpeza será lento e levará mais tempo para obter uma amostra de qualidade da pureza desejada.

As moléculas aromáticas e aromáticas polinucleares são as que fluorescem primariamente. Este é o motivo pelo qual o óleo cru, normalmente, fluoresce muito mais do que o filtrado da lama baseada em óleo (OBM). Para as razões ambientais, designam-se as OBMs sintéticas para serem o mais livres de aromáticos possível, mas podem captar alguma contaminação de aromático da perfuração ou podem ter pequenas quantidades de emulsi-ficadores aromáticos ou materiais de controle de perda de fluido acionados a elas. Ainda, o filtrado da lama baseada em água possui pouca ou nenhuma fluorescência, porque a própria água é não-fluorescente. Ademais, pode-se, deliberadamente, adicionar componentes à lama baseada em água ou baseada em óleo como rastreador fluorescente.

No caso específico onde a propriedade que é ajustada é uma função da absorção óptica, determinadas funções particularmente úteis podem ser selecionadas para a absorção. Tal função é a razão de um pico de óleo ajustado pela linha base para um pico de água ajustado à linha base ou o inverso. Essa função é, em particular, útil no monitoramento da limpeza a partir do filtrado de lama baseada em água para o óleo cru natural. O inverso é particularmente útil no monitoramento da limpeza a partir do filtrado de lama baseado em óleo para águas singenéticas, quando deseja-se colher umaamostra da água.

O pico de óleo ajustado pela linha base é um canal de pico de óleo (próximo de 1740nm) menos um canal adjacente de "referência pelalinha base" de baixa absorção (por exemplo, canais de 1300 ou 1600nm). O pico de água ajustado à linha base é um canal de pico de água (perto de 1420 ou 1935nm) menos um canal adjacente de "referência pela linha base" de baixa absorção (por exemplo, canais de 1300 ou 1600nm). Ao substituir tempo igual infinito dentro do modelo previsto permite a estimativa do valor limitativo da propriedade, P, em tempo infinito. Ao dividir o presente valor da propriedade, P, por seu valor final previsto dá a fração de pureza final.

Em uma primeira modalidade da presente invenção, o método e aparelho desta invenção ajustam os dados de medição de fluido à uma cur-va não-assintótica. Um exemplo de uma curva não-assintótica é uma curva que proporciona um ajuste melhorado para os dados através do tempo de bombeamento típico e, que também pode extrapolada, com sucesso, para diversos tempos do tempo de bombeamento, porém que apresenta mais ou menos infinito nos tempos infinitos, tais como um aproximação de séries deforça. Um outro exemplo de uma curva não-assintótica é uma equação que possui um componente oscilatório, tal como uma onda sinoidal, a qual nunca alcançou um limite fixo. A onda sinoidal pode ser ajustada em freqüência, fase e amplitude para proporcionar um ajuste melhorado aos pulsos na resposta monitorada que estão associados com cada golpe da bomba.

Em uma segunda modalidade, o método e aparelho usam reco-nhecimento padrão. Ou seja, o método e aparelho da presente invenção u-sam uma equação tal como log(1-ftp) = (-p) log(t) + log(Ai/A0). O método e aparelho, então, realizam uma série de estimativas diferentes da pureza final ou valor final pois uma propriedade física do fluido é representada por Ao, onde A0 começa com Ai, A + e, A + 2 e, etc. O valor de A0, para o qual o a-juste aos dados é o mais perto do formato de uma linha reta (baseada no valor quadrado de R) torna-se a melhor estimativa de A0. Em uma terceira modalidade, o método e aparelho da presente invenção ajustam uma curva diferenciável aos dados de medição ou dados de propriedade física derivados dos dados de medição. A presente invenção, então, estima A/A0 a partir da razão de (dA/dt) para A. Em uma quarta modalidade, a presente invenção ajusta uma curva assintótica às diferenças de absorção dos canais ópticospróximos (comprimento de onda) ao invés da própria absorção. As diferenças de absorção removem as saliências da linha base ocorridos através da passagem de partículas de areia e bolhas.

Em uma abordagem convencional à contaminação de formação, as equações 1 e 3 são aplicáveis.

Eq. 1 A = A0-A1t"5/12onde A0>0, Ai>0 limA = A0

t-+ °°

Ao invés do tempo, t, o volume, V pode ser usado. Pode-se, ainda, generalizar para o caso onde a melhor força de ajuste, p, é calculada ao invés de assumida.

Eq. 2 A = A0-At t"ponde A0>0, A^O, p>0 limA = A0

t-+ °°

Eq. 3 ftp = A/Ao = fração de absorção final, A0, que é alcançada quando a absorção é A.

Apenas nestes casos onde A0 é a absorção do óleo cru puro a ftp = fração de pureza final também é igual à fp = fração de pureza.

Eq. 4 1 -ftp = [1 - (A/Ao)] = fração longe da absorção final

Para a Eq. 1, a abordagem convencional encontra os melhores A0, Ai ao usar quadros mínimos lineares ajustados aos pontos de dados N, (Ai( tr5/12), onde /' = 1, N. Para a Eq. 2 encontra-se os melhores A0, Ai ao usar quadros mínimos lineares ajustados aos pontos de dados N, (Aif tf13), onde / = 1, N após assumir ou encontrar um melhor valor de ajuste para p, conforme descrito em qualquer outro lugar nesta invenção.

Com referência agora às figuras 3-10, descrevem-se diversas funções realizadas nas modalidades para a invenção. Conforme mostrado na figura 3, em um exemplo da presente invenção, extrai-se o fluido a partir de uma formação 310. Uma propriedade do fluido é mensurada 320, a partir da qual uma estimativa de uma contaminação de fluido é feita a partir de um ajuste da propriedade com uma curva não-assintótica que inclui ajustes executados para obter o declive de dados 330.

Muito embora a discussão abaixo use tempo transcorrido como a variável dependente, entende-se que o volume do fluido bombeado ou al-gum outro parâmetro também pode ser usado. Conforme mostrado na figura 4, em uma modalidade da presente invenção, a presente invenção realiza uma curva não-assintótica em peça ajustada aos dados para determinar valores amenizados e declives de dados nos centros de cada segmento. Realiza-se uma regressão no logaritmo do derivado dos dados pelo tempo versus o logaritmo do tempo para obter uma interceptação e declive de regressão em linha reta. Um valor de ftp de pureza final fracional é estimado a partir da interceptação e declive de regressão em linha reta e a partir das médias de A(t) e Ait"p em uma pluralidade de tempos. Por exemplo, um método e aparelho são proporcionados para ajustar a absorção sobre um segmento do tempo de laminação que usa uma equação não-assintótica tal como as séries de força, A = b0 + bit + b2t2. Então, pelo cálculo, A' = dA/dt = bi + 2b2t e A7A = (bi + 2b2t)/( b0 + bit + b2t2). Para uma equação da forma, A(t) = A0 -Ait"p, pode-se realizar uma regressão em linha reta de ln(dA/dt) versus ln(t) para obter interceptação e declive da linha de melhor ajuste e calcular os valores de melhor ajuste, -p = (1 + declive) e -Ai = exp(interceptação - ln(1 + declive)). Nesta maneira, não há necessidade de assumir um valor de -5/12, ou -2/3, ou de qualquer outro valor fixo para -p. Ao invés disso, pode-se estimar ftp=A/A0 a partir de valores de melhores ajustes para p e Ai, e a partir duas vezes a média de A(t) e Ait"p em uma pluralidade de tempos 410.

Conforme mostrado na figura 5, em uma segunda modalidade da presente invenção, proporcionam-se um método e aparelho que usam uma curva não-assintótica para ajustar os dados 510. Nesta modalidade, o método e aparelho ajustam uma versão modificada da Eq. 1 aos dados, onde a equação modificada não aborda uma assíntota no tempo infinito, tal como os exemplos mostrados nas Equações 8 e 9 abaixo, que usam a forma A = A0 - h(y) onde t —►°° e h(t) não vai a zero.

Eq. 8 A = A0 - A-|Itx, onde x = -n a +m.

Eq. 9 A = A0-A1[fp + k-1seno(wt)].

O temo seno(ojt) pode proporcionar um melhor ajuste aos dados que possuem cavilhas periódicas em respostas que ocorrem, comumente, com todo golpe da bomba, conforme as partículas são agitadas. Logicamen-te, este termo oscilante previne a curva de se estabilizar em um valor fixo, sem importar o tempo, então não se trata de uma curva assintótica. O valor de w pode ser escolhido para coincidir com a freqüência do golpe da bomba. Para a Eq. 9, a presente invenção determina um melhor A0, h-\ ao usar um ' quadro mínimo linear ajustado aos pontos de dados de N, (A„ tf5/12 + k" 1seno(cot)).

Conforme mostrado na figura 6, em uma terceira modalidade, a presente invenção proporciona para um reconhecimento padrão 610. Conforme mostrado na figura 6, a presente invenção realiza um reconhecimento padrão para uma estimativa por tentativa e erro de A0, ao invés de um cálculo direto de A0. Nesta modalidade, o padrão a ser observado é a semelhança mais próxima de uma linha reta, conforme determinado pelo coeficiente de correlação mais alto, R, para um ajuste por quadros mínimos lineares. Um método e aparelho realizam uma série de ajustes por quadros mínimos lineares aos dados de absorção que usam uma série de estimativas diferentes de A0 que começa com A + e, A + 2z, até A + N e, onde A + N £ < 3,5 OD, onde 3,5 é usado como um exemplo para o limite da extensão dinâmica superior da ferramenta. O valor de A0, para o qual o ajuste é o mais próximo de uma linha reta no espaço log-log, torna-se, então,a melhor estimativa de A0. A proximidade do ajuste ao formato de uma linha reta é determinada pela proximidade de R2à unidade, onde R2é o quadrado do coeficiente de correlação que está na faixa de 0 (sem correlação) a 1 (correlação perfeita). Ou seja, para um série de suposições de A0, encontra-se o melhor A0 baseado no melhor fí2em um ajuste por quadros mínimos lineares para os pontos de dados medidos N, (logft], log[1.-(A(ti)/A0)])..

Um exemplo de um declive de uma certa linha seria (-pA^Ao), que para qualquer valor fixo de p, também permite a determinação imediata de A^ Pode-se assumir que um valor fixo de p ou pode-se calcular um valor de melhor ajuste para p a partir do declive da regressão em linha reta de ln(dA/dt) versus ln(t). Nota-se que A0 não é calculado aqui. Apenas R2é calculado para suposições (estimativas) diferentes de Ao. Ou seja, estimativas diferentes de A0 = A + ne, são experimentadas e a que produz o melhor R2éusada. Para estimar A0 para uma solução melhor que z, poder-se-ia usar convergência binaria para testar, iterativãmente, os valores de A0 entre os dois melhores valores de A0 previamente determinados.

Eq. 10 log(1-ftp) = (-p)log(t) + log(A1/Ao)

Para a Eq. 10, para uma série de suposições de A0 diferentes, a presente invenção encontra o melhor R2ao usar um ajuste por quadros mínimos lineares nos pontos de dados N.

Em uma quarta modalidade, o método e aparelho da presente invenção ajustam uma curva diferenciável aos dados mensurados. A presente invenção estima ftp a partir de (dA/dt)/A através do ajuste de uma curva continuamente diferenciável aos dados de absorção (ou dados de absorção regulares). Um ajuste em partes para os diversos segmentos dos dados podem ser realizados. Observar que a curva por ajuste não necessitar abordar um próprio valor final. Este propósito é, simplesmente, para proporcionar uma função de ajuste liso sobre um intervalo de tempo longo o suficiente de pontos de dados, de modo que os valores ajustados tanto de A(t) quanto de dA(t) possam ser calculados a qualquer momento, t, dentro do intervalo e, então, substituídos dentro das equações 14-16.

Conforme mostrados nas equações 14-16, o valor de absorção final pode, agora, ser determinado a partir da razão do atual declive, DA(t)/dt, pelo valor atual, A(t). O ajuste local e as funções de alisamento u-sados para calcular dA(t)/dt e A(t) não necessitam ter os próprios valore finais. Eles podem até tender a somar ou diminuir o infinito, no tempo infinito, conforme ocorreria com um ajuste de série de força ou um grupo de ajustes de séries de força.

Portanto, conforme mostrado nas figuras 7, 8 e 9, sem ao menos calcular A0 ou ajustar uma curva assintótica aos dados de absorção, torna-se impossível determinar a fração da pureza final, ftp, que é alcançada no tempo t, a partir da razão da taxa de mudança da absorção pela absorção.

Por exemplo, tem-se

Eq. 11 A = A0-A1h(t)ondelimh(t) = 0

t-->ooPega-se o primeiro derivado com relação ao tempo Eq. 12 (dA/dt) = -A^dh/dt) portanto At = -(dA/dt)/(dh/dt) Eq. 13 ftp = A/[A + A^t)] = 1/[1 + A^tJ/A] = 1/[1-(dA/dt)(dh/dt)-1h(t)A-1] então'710 Eq. 14 ftp = [1-A"1(dA/ht)g(t)(dh/dt)'.1r1

Exemplo: tem-se h(t) = t"p de modo que dh/dt = -p t"p t"1 Então, 810 Eq. 15 ftp = [1 +A"1(dA/dt)tp1]'1

Para o caso especial onde p=5/12, 910 Eq. 16 f,p = [1+(12/5)A-1(dA/dt)t]-1

Em uma quinta modalidade, conforme mostrado na figura 10, ométodo e aparelho da presente invenção, encontra-se valores finais de dados de diferenças de absorção ao invés de a própria absorção. As diferenças de absorção de canais vizinhos são projetadas para remover desvios de base ocasionados pelas partículas de areia ou bolhas 1010. O método e aparelho da presente invenção realizam um ajuste para as diferenças de absorção ao invés de a própria absorção. As diferenças de absorção são previstas, por exemplo, a diferença entre os canais ópticos, OD16-OD15, que correspondem aos diferentes comprimentos de onda ópticos fora de seus valores finais, ao invés de prever um único canal OD fora de seu valor final.

O dado de diferença de absorção é usado de maneira independente ou em conjunto com a abordagem descrita nas figuras 3-9 a fim de determinar a pureza final fracional, ftp.

Em uma sexta modalidade, conforme mostrado na figura 13, a presente invenção ajusta uma curva não-assintótica, continuamente diferen-ciável 1302 aos dados brutos que podem ser valores mensurados de qualquer parâmetro adequado que esteja relacionado ao fluido ou dados derivados de tais valores mensurados. O ajuste pode ser para o tempo transcorrido ou o ajuste para o volume de fluido bombeado. A presente invenção usa, por exemplo, porém não se limita, ao ajuste de uma curva não-assintóticapara os pontos de dados brutos, tais como A(t) = c^+ c2t1/2 + c3t1/3 + c4t1/4. Ao usar o cálculo, calcula-se de maneira analítica o primeiro derivado como dA/dt = (c2/2)t"1/2 + (c^t"273 + (c4/4)t"3M para estimar o valor "final" do para-metro, tal como absorção, índice refrativo, resistividade, etc. em algum tempo longo (por exemplo, 24 horas) que é muito mais longo que um tempo (2 horas) em que pode-se, normalmente, terminar o bombeamento para alcançar a limpeza de amostra. Na medida em que o tempo progride, tanto (A0-A) quanto t(dA/dt) regridem, onde A é absorção no tempo t. Assumindo que eles regridem em alguma taxa, são proporcionais, o que significa que (A0-A) = mt(dA/dt) onde "m" é uma constante. A presente invenção experimenta diversas suposições para A0 até que encontra uma suposição que produz um ajuste por quadrado mínimo linear melhor entre y = (A0-A) e x = [(dA/dt)].

Dá-se o melhor ajuste através de y = mx+b onde a interceptação, b, é quase zero, que crê-se ser mais sensível do que encontrar o R2 máximo para os ajustes lineares entre duas variáveis que são diretamente proporcionais. A presente invenção seleciona um ponto de dados brutos em algum tempo, t, (de preferência, o tempo mais distante, t) em que os dados reais intersectam (ou ficam próximos) a melhor linha de ajuste. A absorção prevista de um tempo, ligeiramente mais tardio, t+At, usa-se AA = (A0-A)/[1+m(1+t/At)] que é obtido ao substituir dA/dt por AA/ At, ao substituir t por t+At, e substituir A por A+ AA = m t(dA/dt). Aplica-se, recursivamente, a fórmula AA a fim de prever a absorção em t+At e, então usa a nova absorção calculada para computar a absorção em algum tempo ligeiramente mais tardio, t+2At, e assim vai, por todos os tempos futuros.

O declive 1304, m, deste ajuste é positivo, isso significa que uma má seção ou indesejável de dados brutos foi selecionada, a qual se curva ascendente ou descendentemente em direção ao mais ou menos infinito.

Selecionar um ponto de dados brutos em algum tempo t, (de preferência, o tempo mais tardio t), em que o dado real intersecte (ou chegue perto) a melhor linha de ajuste. A presente invenção, então, calcula a absorção em algum tempo ligeiramente mais tardio, A(t+At) = A(t)+ AA, em termos de t, A(t), A0, e m ao usar AA = (A0-A)/[1+m(1+t/At)]. A presente invenção, então, aplica de maneira recursiva, a fórmula AA a fim de gerar previsões futuras, A(t), dos dados brutos. Para os dados que surgem e se estabilizam com o passar do tempo, a fração da pureza final em qualquer tempo futuro, t, é dada atra-vés de A(t)/A0. Para a queda e estabilização de dados com o passar do tempo, a fração da pureza terminal em qualquer tempo futuro, t, é dada através de [As-A(t)]/[AS-A0]. Agora, As começa a absorção na margem esquerda (o tempo mais cedo) da janela de dados selecionados pelo usuário.

Em uma sexta modalidade, conforme mostrado na figura 13, não se assume que "m" é um número negativo, que é o motivo de a fórmula de recursão ser escrita como [1 + m(1+ t/At)] ao invés de ser escrita como [1 -m(1+ t/At)]. Na sexta modalidade, usa-se a proximidade do zero da intercep-tação "b". Usa-se a interceptação "b" porque é muito mais sensível que a proximidade do R2à unidade para encontrar a melhor linha de ajuste quando sabe-se que a interceptação desta linha deveria ser zero. Ainda, usa-se uma fórmula recursiva para predizer a absorção em tempos futuros. Para os dados que decaem e se estabilizam através do tempo, usa-se a absorção na margem esquerda da janela selecionada pelo usuário, como a absorção quese inicia ao invés de usar zero como a absorção iniciante, conforme faz-se quando os dados surgem e se estabilizam.

Com referência agora à figura 11, ilustra-se um esquema mais detalhado do conjunto piroelétrico para determinar um espectro infravermelho mediano. Em uma modalidade, a presente invenção proporciona uma fonte de luz 402, tal como uma fonte de luz infravermelha que pode ser uma fonte de luz em estado constante ou uma fonte de luz modulada ou pulsada. No caso de uma fonte de luz em estado constante, proporciona-se um mo-dulador de luz. O modulador pode se qualquer dispositivo adequado que varie a intensidade da fonte de luz, que inclui - mas não se limita - um circuito pulsador eletrônico, bem-conhecido na técnica, que varia a intensidade da fonte de luz ou um picador eletromecânico 404 que interrompe o trajeto da fonte de luz para o fluido do interior do poço. O modulador é proporcionado para modular a intensidade da luz proveniente da fonte de luz que infringe no fluido e no fotodetector. Um refletor ou colimador 403 podem ser proporcionados para focar e/ou concentrar luz proveniente da fonte de luz 402. Uma câmara ou conduto 406 é proporcionado para a apresentação de um fluido do interior do poço para exposição do fluido do interior do poço à luzproveniente da fonte de luz. Uma janela óptica 408 é proporcionada, através da qual o fluido do interior do poço 407 fica exposto à luz. Para os propósitos da presente pedido, o termo "fluido" inclui líquidos, gases e sólidos que podem precipitar a partir de um fluido ou um gás.

A presente invenção inclui, adicionalmente, um detector tal como um detector piroelétrico 412. O detector piroelétrico 412 pode, ainda, compreender um conjunto detector piroelétrico. Um espectrômetro 414 e um processador 422 são proporcionados para analisar sinais a partir do detector piroelétrico a fim de determinar uma propriedade do fluido do interior do poço 407. Um filtro variável linear de infravermelho 416 é proporcionado e interposto entre a radiação de luz 440 do fluido do interior do poço e o detector piroelétrico 412. Proporciona-se um amplificador de alto ganho 420 para amplificar o sinal proveniente do detector piroelétrico 412, quando desejado. O espectrômetro 414 inclui um processador 422 com memória. O processador 422 inclui programas que implementam técnicas de modelagem macias para aplicar uma equação quimiométrica, rede neural ou outros programas logiciais (soft) de modelagem para as medições da luz infravermelha detectada pelo detector piroelétrico para estimar outras propriedades físicas e químicas do fluido do interior do poço a partir do sinal de detector piroelétrico. A saída do espectrômetro com relação ao sinal piroelétrico é, ainda, introduzida no programa logicial de modelagem, rede neural ou equação quimiométrica para estimar as propriedades do fluido do interior do poço.

Voltando-se agora à figura 12, ilustra-se um esquema mais detalhado do transdutor acústico a fim de determinar a velocidade do som nofluido. A presente invenção proporciona um transdutor 701, uma linha de fluxo de amostra 703 ou um trajeto de fluxo de amostra 705 que contém uma amostra de fluido para medir a densidade de fluido e velocidade de som do fluido 708 dentro do tubo ou trajeto de fluxo de amostra ou tanque de amostra 711. A espessura 707 da parede da linha de fluxo 706 é conhecida. Umprocessador 702 e eletrônicos pulsantes 704 são proporcionados para enviar um pulso acústico a partir do pulsador 701a através da parede 706 dentro do fluido 705 no trajeto do fluxo 705 ou a partir do pulso 701 b através da parede706 da espessura 707b para a câmara de amostra 711. O transdutor 701 recebe ecos a partir do pulso acústico, que é monitorado através do processador. A presente invenção compreende, adicionalmente, uma compensação de parede, que é um espaçador acústico interposto entre o transdutor e a parede que é feita do mesmo material que a parede. Esse espaçador, simplesmente, aumenta a distância circular da jornada e o tempo de jornada correspondente para as reverberações de eco de pulso dentro da compensação combinada somado ao material da parede próxima. Isto serve para prolongar o tempo entre os pulsos de eco que decaem sucessivamente e, então, serve para aperfeiçoar a separação do pulso, para evitar a sobreposição de pulsos e para melhorar a quantificação de energia em cada pulso.

O processador diminui a densidade do fluido na linha de fluxo de amostra. A presente invenção captura uma amostra de fluido em uma linha de fluxo a partir da formação ou do furo. A presente invenção, então, envia um pulso acústico para dentro da amostra de fluido na linha de fluxo ou tanque de amostra. O processador da presente invenção, então monitora o retorno do eco dentro da parede da linha de fluxo ou tanque de amostra e integra a energia de cada pulso de eco acústico. O processador determina o declive da redução dos pulsos de eco acústicos integrados que oscilam dentro da parede da linha de fluxo. A presente invenção determina, então, o coeficiente de reflexo para a parede interna/interface de fluido. A presente invenção determina a velocidade do som no fluido. A presente invenção determina a densidade do fluido na linha, conforme descrito acima. A presente invenção determina a viscosidade do fluido na linha de fluxo, conforme supradescrito.

A presente invenção foi descrita como método e aparelho que operam em um ambiente do fundo do poço na modalidade preferida, no entanto, a presente invenção pode, ainda, ser incorporada como um conjunto de instruções em um meio legível por computador, que compreende, ROM, RAM, CD ROM, Flash ou qualquer outro meio legível por computador, agora, conhecido ou não pelo fato de fazer com que um computador implemente o método da presente invenção, quando executado. Muito embora uma moda-lidade preferida da invenção tenha sido mostrada através da invenção acima, é, apenas, para o propósito de exemplificaçao e não pretende-se limitar o escopo da invenção, que é definido pelas reivindicações que se seguem.

Claims

1. Método para estimar uma característica de um fluido do interior do poço que compreende:receber o fluido de uma formação pelo tempo; fazer uma pluralidade de medições,relacionadas à uma propriedade do fluido;ajustar uma curva não-assintótica à uma pluralidade de medições; eestimar a característica do fluido a partir da curva ajustada.

2. Método, de acordo com a reivindicação 1, que compreende adicionalmente:tomar um primeiro derivado da curva ajustada; e fazer estimativas para um valor final do parâmetro, A0, até que uma estimativa seja encontrada que produza um ajuste, substancialmente melhor, por quadrado mínimo entre y e x, onde y = valor do parâmetro A no tempo t, menos o valor final, A0, e onde x = t(dA/dt), e t=tempo.

3. Método, de acordo com a reivindicação 2, onde o melhor ajuste é dado através de y=mx + b, onde a interceptação b é a mais perto de zero e m é um declive.

4. Método, de acordo com a reivindicação 3, que compreende adicionalmente:selecionar um ponto de dados que é perto da curva de ajuste; e usar A e tempo, t, para prever um valor futuro do parâmetro A em um tempo mais tardio, t + At.

5. Método, de acordo com a reivindicação 4, onde a presivão de um valor futuro de um parâmetro A compreende adicionalmente:determinar AA = (A0-A)/[1 + t/At)].

6. Método, de acordo com a reivindicação 5, que compreende adicionalmente:determinar, recursivamente, AA para prever o parâmetro A para um tempo futuro.

7. Método, de acordo com a reivindicação 1, onde o ajuste deuma curva não-assintótica compreende a execução de uma curva não-assintótica em partes para obter valores regulares e declive de dados em centros de segmentos de dados, e o método compreende adicionalmente:regredir o logaritmo de um derivado de dados relacionado a um tempo ou volume contra o logaritmo de um de um conjunto que consiste no tempo e volume para obter uma interceptação e declive de regressão em linha reta.

8. Método, de acordo com a reivindicação 7, que compreende adicionalmente:computar uma pureza final fracional, ftp, a partir de uma interceptação e declive de regressão em linha reta e a partir de médias de valores de dados regulares em uma pluralidade de tempos.

9. Método, de acordo com a reivindicação 1, onde a curva não-assintótica é definida através da equação A = A0 - Ailtx, onde x = -n a + m.

10. Método, de acordo com a reivindicação 1, onde a curva não-assintótica é definida através da equação A = A0 - A^f" + k"1seno (cot)].

11. Método, de acordo com a reivindicação 1, onde o ajuste de uma curva não-assintótica compreende ajustar uma série de estimativas por quadros mínimos lineares de um valor final, A0 são ajustados à equação, log fc = log[1-(A/A0)] = (-5Ai/12A0) log t, onde a estimativa que possui o ajuste mais perto de uma linha reta em espaço log) é a melhor estimativa de A0 na contaminação de fluido final.

12. Método, de acordo com a reivindicação 1, onde a característica é uma de uma pureza final fracional, ftp e contaminação final fracional ftc.

13. Método, de acordo com a reivindicação 1, onde a característica é a pureza final fracional que é estimada a partir de ftp da equação ftp = {l-A-1(dA/dt)h(t) (dh/dt)-1]-1.

14. Método, de acordo com a reivindicação 1, onde a característica é pureza final fracional, ftp e é determinada a partir da equação ftp = [1+A"1(dA/dt)tp1]-1

15. Método, de acordo com a reivindicação 1, que compreende adicionalmente estimar a pureza final fracional, ftp e é determinada a partirda equação ftp = [1 + (12/5) A"1(dA/dt)t]"V

16. Método, de acordo com a reivindicação 1, onde a propriedade do fluido compreende pelo menos viscosidade, densidade, velocidade do som, fluorescência, índice refrativo, módulo à granel, resistividade e diferen-ças de propriedade ópticas.

17. Aparelho para estimar uma característica do interior do poço de fluido de formação que compreende:uma sonda em comunicação de fluido com uma formação; um sensor que faz uma pluralidade de medições de uma propriedade do fluido por tempo; eum processador configurado para ajustar as medições da propriedade para uma curva não-assintótica para estimar a característica do fluido.

18. Aparelho, de acordo com a reivindicação 17, onde o processador é configurado, adicionalmente, para tomar um primeiro derivado da curva ajustada e fazer estimativas para um valor final, A0, até que uma estimativa do valor final seja encontrada, que produza um ajuste, substancialmente melhor, por quadrado mínimo entre y e x, onde y = valor de p da propriedade A no tempo t, menos o valor final, A0, e onde x = t(dA/dt), e t=tempo.

19. Aparelho, de acordo com a reivindicação 18, onde o melhor ajuste é dado através de y=mx+b, onde a interceptação, b, é mais perto de zero e m é o declive.

20. Aparelho, de acordo com a reivindicação 19, onde o processador é configurado, adicionalmente, para selecionar um ponto de dados que é próximo à curva ajustada e usa A, e tempo, t, para prever uma absorção futura em um tempo mais adiante, t + At.

21. Aparelho, de acordo com a reivindicação 20, onde para prever o valor futuro, o processador usa a equação AA = (A0-A)/[1 + m(1 + t/At)].

22. Aparelho, de acordo com a reivindicação 21, onde o processador é configurado, adicionalmente, para determinar, de maneira recursiva, AA para prever o valor de A.

23. Meio legível por computador que contém instruções quequando executadas por um computador realiza um método para estimar uma característica de um fluido de formação que compreende:receber o fluido a partir de uma formação;fazer uma pluralidade de medições para uma propriedade do fluido pelo tempo; eajustar uma curva não-assintótica para as medições da propriedade para estimar a característica do fluido.

24. Meio, de acordo com a reivindicação 23, onde o método compreende adicionalmente: tomar um primeiro derivado da curva; efazer estimativas para valores finais, A0, até que uma estimativa seja encontrada, que produza um ajuste, substancialmente melhor, por quadro mínimo entre yex, onde y = valor final A, em tempo t, menos o valor final, A0, e onde x = t(dA/dt), e t=tempo.

25. Meio, de acordo com a reivindicação 24, onde o método para o melhor ajuste é dado através de y=mx+b, onde a interceptação b,é próxima de zero e m é o declive.

26. Meio, de acordo com a reivindicação 25, onde o método compreende adicionalmente:selecionar um ponto de dados que é próximo à curva ajustada e usa A, e tempo t, para prever um valor futuro em um tempo mais adiante, t + At.

27. Método, de acordo com a reivindicação 26, onde na previsão do método, A consiste, adicionalmente, em determinar AA = (A0-A)/[1 + t/At)].

28. Meio, de acordo com a reivindicação 27, onde o método compreende adicionalmente:determinar, recursivamente, AA para prever o parâmetro A para um tempo futuro.