Espalhamento geométrico em modelos plano-estratificados

A medição de parâmetros físicos de reservatórios se constitui de grande importância para a detecção de hidrocarbonetos. A obtenção destes parâmetros é realizado através de análise de amplitude com a determinação dos coeficientes de reflexão. Para isto, faz-se necessário a aplicação de técnicas especiais de processamento capazes de corrigir efeitos de divergência esférica. Um problema pode ser estabelecido através da seguinte questão: Qual o efeito relativamente mais importante como responsável pela atenuação de amplitudes, o espalhamento geométrico ou a perda por transmissividade? A justificativa desta pergunta reside em que a correção dinâmica teórica aplicada a dados reais visa exclusivamente o espalhamento geométrico. No entanto, a análise física do problema por diferentes direções põe a resposta em condições de dúvida, o que é interessante e contraditório com a prática. Uma resposta embasada mais fisicamente pode dar melhor subsídio a outros trabalhos em andamento. O presente trabalho visa o cálculo da divergência esférica segundo a teoria Newman-Gutenberg e corrigir sismogramas sintéticos calculados pelo método da refletividade. O modelo-teste é crostal para que se possa ter eventos de refração crítica além das reflexões e para, com isto, melhor orientar quanto à janela de aplicação da correção de divergência esférica o que resulta em obter o então denominado “verdadeiras amplitudes”. O meio simulado é formado por camadas plano-horizontais, homogêneas e isotrópicas. O método da refletividade é uma forma de solução da equação de onda para o referido modelo, o que torna possível um entendimento do problema em estudo. Para se chegar aos resultados obtidos foram calculados sismogramas sintéticos através do programa P-SV-SH desenvolvido por Sandmeier (1998), e curvas do espalhamento geométrico em função do tempo para o modelo estudado como descrito por Newman (1973). Demonstramos como uma das conclusões que a partir dos dados do modelo (velocidades, espessuras, densidades e profundidades) uma equação para a correção de espalhamento geométrico visando às “verdadeiras amplitudes” não é de fácil obtenção. O objetivo maior então deveria ser obter um painel da função de divergência esférica para corrigir as verdadeiras amplitudes.

Estudo comparativo entre estereotomografia e da tomografia da onda NIP: aplicação em dados sintéticos e reais

A determinação de um acurado modelo de velocidades é um requisito fundamental para a realização do imageamento sísmico. Métodos novos como a Estereotomografia préempilhamento e a Tomografia da onda NIP são ferramentas poderosas e bastante sugestivas para este propósito. Basicamente, a Estereotomografia pré-empilhamento se baseia no conceito de eventos localmente coerentes interpretados como reflexões primárias e associados com pares de segmentos de raios, que tem um mesmo ponto de reflexão em profundidade. Na Tomografia da onda NIP um evento sísmico é representado por uma onda hipotética NIP, que está relacionada a um ponto de reflexão em profundidade. Os atributos da onda NIP são determinados no decorrer do Empilhamento de Superfície de Reflexão Comum (empilhamento CRS). Este trabalho tem como objetivo, fazer um estudo comparativo de ambos os métodos de determinação do modelo de velocidades em profundidade. Assim, é realizada uma revisão dos fundamentos teóricos de ambos os métodos tomográficos, destacando as principais diferenças e aplicando estas aproximações em um dado sintético e um dado real marinho (linha sísmica 214-2660 da Bacia do Jequitinhonha). Para avaliar os modelos de velocidades encontrados pelas aproximações, foi utilizada a migração pré-empilhamento em profundidade do tipo Kirchhoff e também as famílias de imagem comum (CIG). Os resultados mostraram que ambos os métodos tomográficos fornecem modelos de velocidades representativos. Contudo, constatou-se que a estereotomografia tem melhor desempenho em meios com variações laterais de velocidades, porém, aplicável somente em dados pré-empilhados com alta razão sinal/ruído.

Estudo do método de empilhamento SRC e aplicação em dados sintéticos e reais

O empilhamento sísmico é um método de imageamento utilizado para simulação de seções sísmicas com afastamento nulo (AN) a partir de dados de cobertura múltipla, produzindo um aumento da razão sinal/ruído na seção simulada. Este método é realizado usando o processamento convencional denominado ponto médio comum (PMC), através das chamadas correções sobretempo normal e sobretempo oblíquo (STN/STO). Alternativamente, um novo método de empilhamento sísmico, denominado superfície de reflexão comum (SRC) pode ser usado para simular seções sísmicas em afastamento nulo (AN). O método SRC tem como principais características: 1) a construção de uma aproximação do tempo de trânsito de reflexões primárias na vizinhança de um raio de reflexão normal; 2) esta aproximação de empilhamento tem como parâmetros os ângulos de emergência do raio normal e os raios de curvaturas de duas ondas hipotéticas denominadas onda ponto de incidência normal (onda PIN) e onda normal (onda N). O método pressupõe ainda que a velocidade da onda longitudinal na parte superior do modelo próxima à superfície da Terra seja conhecida. Para determinar adequadamente os parâmetros supracitados, fundamentais para o método de imageamento SRC, são necessárias estratégias de busca que utilizem casos particulares da aproximação de segunda ordem hiperbólica dos tempos de trânsito. As estratégias de busca apresentadas são: SRC pragmático-extendito e SRC global-local. Neste trabalho é apresentada uma revisão extensa dos fundamentos teóricos de método de empilhamento SRC, como também das estratégias de busca dos parâmetros SRC, dos algoritmos de implementação, e da delimitação de aberturas do operador de empilhamento SRC usando conceitos de zonas de Fresnel e o fator de estiramento da correção STN. Como um exemplo de aplicação dos atributos ou parâmetros SRC é apresentado o resultado da correção STN, sem estiramento do sinal sísmico, obtido utilizando as velocidades estimadas dos parâmetros SRC. Também, são apresentados os resultados da aplicação dos três métodos de empilhamento (PMC (STN/STO), SRC pragmático-extendido, SRC global-local) nos dados sintéticos Marmousi e real terrestre da bacia paleozóica Amazônica. Finalmente apresentam-se as discussões sobre os resultados, quanto à resolução da imagem sísmica e custo computacional.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Inversão quadrática de dados de AVO usando algoritmo de Marquardt

Neste trabalho são apresentadas expressões exatas e aproximações quadráticas das conhecidas equações de Knott-Zöeppritz, que calculam as amplitudes dos coeficientes Rpp e Rsp em termos dos contrastes médios relativos, bem como é procedida a inversão dos parâmetros de dados de AVO a partir destas aproximações quadráticas. Nesta inversão é utilizado o algoritmo de Levenberg-Marquardt, e são considerados apenas os eventos refletidos Rpp e convertidos Rsp, não associados. Nos estudos dos parâmetros físicos dos meios contrastes de impedância (δz), módulo de cisalhamento (δμ) e velocidade da onda p (δα), verificou-se quais desses parâmetros podem ser invertidos. Os resultados obtidos mostram que o contraste de impedância (δz) é muito bem resolvido estando ele relacionado com o contraste de velocidade da onda p (δα) ou com o contraste de cisalhamento (δμ), no caso de eventos refletidos considerando modelos de alto, moderado e baixo contrastes. Por outro lado ao se fixar o contraste de impedância e relacionar os outros dois parâmetros em consideração, os resultados mostram que esses são mal resolvidos, ou seja, a região de ambiguidade torna-se muito grande e os parâmetros tornam-se ambíguos e instáveis. No caso do evento convertido e na combinação do evento refletido com o convertido, para os modelos de baixo e moderado contrastes, (δz) é muito bem resolvido, caso que não acontece para modelo de alto contraste. Diante desses resultados verifica-se que no procedimento de inversão quadrática de dados de AVO, fixado (δμ), a recuperação dos dois parâmetros variados é muito boa, no caso do evento refletido, e razoavelmente boa no caso do evento convertido, por esse motivo optou-se pela fixação do módulo de cisalhamento.

Modelagem de dados MT 2D em multicamadas com anisotropia arbitrária

Neste trabalho apresentamos a modelagem direta bidimensional do método Magnetotelúrico (MT) com anisotropia arbitrária na condutividade elétrica para modelos de multicamadas estratificadas horizontalmente, separando a resposta do campo resultante em duas partes, primária e secundária. A parte primária refere-se ao campo 1D do modelo na ausência de heterogeneidades; a secundária, à contribuição no campo resultante gerada pela anomalia. Esta modelagem foi feita via técnica numérica dos Elementos Finitos (EF). Apresentamos também a modelagem direta do caso MT 1D para o mesmo modelo e anisotropia arbitrária, porém, com solução analítica pelo método da matriz de propagação; tal resultado é usado como fonte na formulação 2D. No estudo a respeito do comportamento dos campos 1D em profundidade na presença da anisotropia, verificamos nas componentes horizontais que suas curvas dentro de cada camada se assemelham a de um campo na presença de fontes, pois experimentam um aumento na amplitude; além de verificar o surgimento da componente normal E_z do campo, sem que haja uma densidade de corrente nesta direção. Com arbitrariedade na anisotropia queremos dizer que seus eixos principais podem assumir qualquer orientação em relação ao eixo de referência adotado para o problema; tal orientação é descrita sob três rotações consecutivas, chamadas de strike, dip e slant, que simulam inclinações bastante presentes na geologia dos ambientes terrestres. Verificamos que a anisotropia provoca alterações no comportamento das pseudo-seções de resistividades aparentes, interferindo na interpretação correta dos dados.

Modelo de mudança organizacional: estudo de caso no setor de compras

O presente trabalho apresenta um modelo de mudança organizacional para melhorar o processo empresarial de compra de uma construtora, proveniente do sistema de comunicação entre o escritório central e suas obras, a fim de identificar as oportunidades e deficiências existentes no sistema, seja de ordem tecnológica e/ou gerencial. O mapeamento das informações permite detectar os percursos desnecessários que impossibilitam o sistema de comunicação funcionar de forma ágil e otimizada, bem como perceber as oportunidades de melhorias contínuas ao processo. O método de pesquisa se baseou no aperfeiçoamento de processos empresariais, que consiste na modelagem do processo citado, por meio da coleta de dados, análise de documentos, desenho de fluxos de informações existentes no modelo e discussão dos resultados. Para validação da aplicação do método, fez-se o estudo de caso único em uma construtora do subsetor edificações, o qual gerou fluxogramas do processo de compras atual e proposto como modelo de mudança organizacional.

Mapeamento dos parâmetros florísticos e estruturais de floresta de mangue com dados LIDAR e SRTM

Este estudo apresenta a estimativa dos parâmetros florísticos e estruturais (determinação da espécie, altura, diâmetro a altura do Peito - DAP e biomassa) do mangue a partir de informações da superfície adquiridas remotamente com os sensores Laser Detection and Range (LIDAR), Shuttle Radar Topography Mission (SRTM) e ortofotos na Ilha dos Guarás, conjunto de arquipélagos localizado a 30 km da desembocadura do rio amazonas. Para esse trabalho foram utilizadas informações do SRTM, LIDAR e fotografias aéreas processadas e ortorretificadas durante dois sobrevôos realizados entre o mês de julho e agosto de 2011. Com a ortofoto foi feito o mapa do reconhecimento de unidades geobotânicas que delimitou apenas a classe mangue. Em seguida, foi realizada a correção da altura elipsoidal para a altura ortométrica, onde a nuvem de pontos foi interpolada pelo método vizinho mais próximo, gerando Modelo Digital de Elevação (MDE) LIDAR (full points) com RMSE de 0,88 cm e por meio de uma linguagem macro foi estatisticamente separadas as informações do último pulso da superfície, conhecido também por ground points. Em seguida, os dados foram interpolados pelo método de krigeagem que gerou o valor de Modelo Digital de Superfície (MDS), o qual foi subtraído do MDE. Com base no Modelo Digital de Vegetação (MDV) foram definidos os sítios de coleta e selecionadas as árvores ascendentes, intermediárias e emergentes, porte no qual foi medido o DAP e altura. No total foram coletadas 212 amostras individuais de mangue e para assegurar o nível de acurácia do conjunto coletado, foi realizado o cálculo de RMSE entre as alturas do LIDAR e Campo, que resultou em RMSE= 1,10 m. Os modelos escolhidos para calibração LIDAR e altura de campo foi do tipo linear, com R² = 91% e RMSE= 0,98 cm e para calibração da DAP e altura de campo foi escolhido o modelo Logarítmico R² = 74,1%. Nos resultados da calibração do SRTM o modelo logarítmico também foi o mais adequado para a relação entre altura média e SRTM com R² = 91% e RMSE de 2,2 m e DAP Médio e SRTM, com R² = 88% e RMSE 2,2 cm. A partir de um inventário foi realizada a estimativa da biomassa por espécie por meio das equações alométricas de Fromard e posteriormente os resultados foram espacializados em forma de mapas com alto nível de detalhamento oriundo das informações LIDAR e SRTM corrigido e ortofotos.

Proveniência das rochas metassedimentares do Cinturão Araguaia com base em datações Pb-Pb em zircão e idades-modelo Sm-Nd

Este trabalho apresenta dados geocronológicos ²⁰⁷Pb/²⁰⁶Pb de grãos detríticos de zircão obtidos pelo método de evaporação de chumbo e idades-modelo Sm-Nd (T_DM) de rochas metassedimentares do Cinturão Araguaia, e discute as possíveis áreas-fonte dessas rochas, buscando investigar a história evolutiva deste cinturão no contexto da amalgamação do Gondwana. As datações em grãos detríticos de zircão de quartzitos da Formação Morro do Campo apontaram idades arqueanas (3,0-2,65 Ga) para o domínio norte (região de Xambioá) e, para o domínio sul (região de Paraíso do Tocantins), revelaram idades meso-neoproterozoicas (1,25-0,85 Ga) e, secundariamente, paleoproterozoicas (1,85-1,70 Ga), sugerindo a existência de áreas fontes distintas para os dois domínios. As idades-modelo Sm-Nd (T_DM) obtidas em metapelitos (ardósias, filitos, micaxistos) dos grupos Estrondo e Tocantins apresentaram distribuição bimodal com maior frequência de idades entre 2,1 e 1,4 Ga, com moda entre 1,7 e 1,6 Ga, e outras menos frequentes entre 2,7 e 2,4 Ga, sugerindo mistura de fontes de idade Paleoproterozoica (ou até Arqueana) com fontes mais jovens, provavelmente do Meso-Neoproterozoico. Os principais candidatos a fonte para as rochas do Cinturão Araguaia seriam os segmentos crustais situados a sudeste (Cráton São Francisco, Maciço de Goiás e Arco Magmático de Goiás). Toda a sucessão de rochas sedimentares da bacia oceânica Araguaia e rochas magmáticas associadas a estes segmentos foram transportados, posteriormente, em direção à margem oriental do Cráton Amazônico, durante a tectônica principal de estruturação do Cinturão Araguaia, resultante da amalgamação do supercontinente Gondwana.

Refinamento da representação de raízes no modelo de biosfera SiB2 em área de floresta na Amazônia

O objetivo desta pesquisa foi estudar a sensibilidade dos fluxos de superfície e hidrologia do solo em relação à representação e distribuição de raízes no modelo de biosfera para uma floresta de terra firme na Amazônia. A finalidade foi avaliar o impacto na representatividade dos fluxos de energia considerando a sazonalidade da região amazônica, usando como suporte medidas intensivas realizadas em uma reserva biológica (Reserva Biológica do Cuieiras, em Manaus). Foram realizadas oito simulações com o modelo de biosfera SiB2 (“Simple Biosphere Model” – versão 2) , onde cada simulação representou um cenário diferente de distribuição de raízes em uma profundidade de 4 m de solo, dividido em três camadas: 0,5 m, 1,5 m e 2,0 m. As raízes foram distribuídas privilegiando a concentração de raízes na camada superficial, em seguida, na camada intermediária e, por fim, uma concentração de raízes abaixo de 2,0 m de profundidade. As simulações foram realizadas para o período de 2003 a 2006, enfatizando o ano de 2005 para avaliar o efeito da representação de raízes nos fluxos de energia (calor latente – LE e calor sensível – H) e de dióxido de carbono (CO₂). A partir da análise integrada dos fluxos simulados com dados observacionais medidos no sítio experimental foi possível perceber que uma redução na precipitação no ano de 2005, apesar de ter sido menor na parte central da Amazônia, implicou na diminuição da umidade do solo, mostrando que a floresta passou por um período de estresse hídrico maior do que os outros anos analisados. O modelo representou a energia disponível com valores muito próximos aos observados, variando sazonalmente em concordância com os dados medidos em 2005. No entanto, LE é superestimado durante a estação chuvosa, mas mostra juntamente com o fluxo de CO₂, a redução com a umidade do solo na estação seca, enquanto H é superestimado em até 20 W m^-2 durante todo o período simulado. Estes resultados mostram que, a consideração de raízes rasas é mais apropriada para regiões que possuem estação seca curta, conforme se caracteriza a área de estudo, e raízes profundas devem favorecer a modelagem dos processos para superfície de áreas com estação seca mais pronunciada. Com isto os resultados revelam que há necessidade de obter mais informações de propriedades físicas do solo, apropriadas às condições da região, para que outros refinamentos sejam efetivos na distinção do comportamento de florestas tropicais sob diferentes regimes de disponibilidade de água no solo.

Desenvolvimento de um modelo experimental de hemodiálise em cães

Com o objetivo de desenvolver um modelo de hemodiálise (HD) em cães, foram estudados 18 animais, sem raça definida, machos, clinicamente sadios, com peso corporal variando entre sete e 14 kg. O acesso vascular foi obtido através de implantação do cateter de duplo lúmen em veia jugular externa. As sessões de HD, em número de cinco por animal, com até três horas de duração, foram realizadas em hemodialisadora de sistema proporcional com ultrafiltração (UF) controlada, com solução dialisante padrão e tampão bicarbonato. A UF foi ajustada para HD isovolêmica, utilizou-se perfil de sódio, e para anticoagulação heparina sódica. Os animais foram mantidos anestesiados com cloridrato de levomepromazina e propofol. Foram avaliados dados hematológicos, bioquímicos, hemogasometria, pressão arterial sistêmica e tempo de coagulação ativado. Foi observada diminuição do número global de hemácias, volume globular, hemoglobina e leucócitos. Em relação aos exames bioquímicos, houve manutenção nos níveis de sódio sérico, e quanto à hemogasometria, a manutenção da SO₂. A pressão arterial sistêmica manteve-se constante. Os resultados obtidos no presente trabalho permitiram concluir que foi possível o desenvolvimento do modelo proposto e mostrou que a HD em cães é um método viável e seguro, que poderá contribuir para o tratamento clínico da insuficiência renal nesta espécie.