Reservatórios turbidíticos – simulação da morfologia por estatísticas multiponto e avaliação das propriedades petrofísicas

Esta dissertação apresenta uma metodologia original para simular a morfologia e as propriedades petrofísicas de reservatórios de hidrocarbonetos em sistemas de canais turbidíticos. Estes sistemas são constituídos por complexos, ou seja, conjuntos de canais de arquitetura meandriforme, e são considerados importantes alvos para a indústria petrolífera. A simulação da morfologia divide-se em duas partes, primeiramente é simulada a trajetória do complexo e depois são simuladas as trajetórias dos canais propriamente ditos condicionadas à trajetória do complexo. O algoritmo de simulação utiliza as classes de ângulos azimutais de linhas poligonais de treino como uma variável aleatória. As trajetórias são simuladas também como linhas poligonais, condicionais a estatísticas multiponto das trajetórias de treino e a pontos de controlo. As estatísticas multiponto são organizadas em árvore, que guarda sequências de classes de orientação que ocorrem na trajetória de treino e as respetivas probabilidades de ocorrência. Para avaliar as propriedades petrofísicas, o modelo morfológico das trajetórias é convertido para uma malha de blocos de alta resolução, identificando-se, em cada bloco, a fácies preponderante de acordo com um modelo conceptual de zonamento da secção dos canais. A conversão prioriza os canais mais recentes (do topo) sobre os mais antigos (da base). A cada fácies é associada uma lei de distribuição da porosidade e permeabilidade, assim são geradas imagens destas propriedades petrofísicas por Simulação Sequencial Direta com histogramas locais. Finalmente, o número de blocos da malha é reduzido por upscaling e as simulações são ordenadas para poderem ser utilizadas nos simuladores de fluxo. Para ilustrar a metodologia, utilizaram-se imagens de sísmica 3D de um reservatório turbidítico na Bacia do Baixo Congo para extrair leis de distribuição das dimensões dos canais e trajetórias de treino. Os resultados representam corretamente a arquitetura complexa destes sistemas.

Simulação sequencial de atributos diagnósticos do solo

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Geoestatística, simulação estocástica e sensoriamento remoto na estimativa de produção do café conilon

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modelação estocástica do depósito mineral do Zambujal (Mina de Neves-Corvo): contribuição da densidade como indicador morfológico

Este estudo tem como objetivo construir um modelo estocástico de alta resolução da morfologia e dos teores em metal, do depósito mineral do Zambujal, Mina de Neves-Corvo. O depósito do Zambujal é um corpo vulcanogénico que se localiza no setor Português, na parte mais a Sul da Faixa Piritosa Ibérica, com sulfuretos maciços no topo e fissurais na base. Para construir um modelo estocástico deste depósito, onde os teores em metal evidenciam zonamento, é importante ter em conta a proporção local de sulfuretos na matriz rochosa. Como esta variável não é quantificada em laboratório, propõe-se a utilização da densidade das amostras como um indicador indireto da proporção de minérios na matriz rochosa. Conhecida esta variável, a modelação dos teores em metal pode ser feita para os chamados teores relativos (teores em metal na fração de sulfuretos). As principais etapas da metodologia proposta são: (a) para cada amostra analisada no laboratório, estimação de soluções para a variável proporção de sulfuretos na matriz rochosa ????(????), tendo em conta a paragénese principal do depósito, os teores e a densidade da rocha; (b) construção de um modelo morfológico 3D de baixa resolução com duas regiões, minérios maciços e minérios fissurais, por digitalização de limites em perfis, interpolação de superfícies e conversão para o modelo de blocos do depósito; (c) construção de um modelo morfológico 3D de alta resolução para todo o depósito da variável ????(????) por Simulação Sequencial Direta (SSD), tendo como informação condicionante os histogramas regionais de ????(????) para minérios maciços e minérios fissurais e as soluções obtidas em (a) para ????(????) na localização das sondagens; (d) SSD dos teores relativos em cobre, zinco e prata; (e) discussão dos resultados e quantificação de recursos. Os resultados foram validados por comparação com os equivalentes obtidos numa estimação por krigagem normal dos teores em metal e mostraram ser da mesma ordem de grandeza. O conjunto de imagens simuladas das variáveis ????(????) e teores permite quantificar a incerteza do conhecimento do depósito relativamente à informação disponível.

Desprendimento de vórtices e controle em esteira de cilindros por simulação numérica direta

O desprendimento de vórtices é responsável por vibrações prejudiciais a estruturas e pelo aumento de seus coeficientes de arraste. Portanto o correto entendimento e o controle do desprendimento de vórtices é de grande interesse de modo a diminuir custos de manutenção de estruturas estáticas e custos com combustíveis de estruturas móveis, por exemplo. O caráter tridimensional do desprendimento de vórtices a baixos números de Reynolds, entre 100 e 300, é aqui estudado através de Simulações Numéricas Diretas (DNS) do escoamento ao redor de cilindros circulares em duas e três dimensões. O controle do desprendimento de vórtices foi estudado nesta dissertação, através do uso de placas de separação do escoamento em simulações numéricas bidimensionais para números de Reynolds 100, 160 e 300, mostrando que para estes números de Reynolds é possível reduzir as freqüências de desprendimento, e em alguns casos até suprimi-las. A simulação tridimensional para Re = 300 mostrou um escoamento aparentemente bidimensional que afeta sensivelmente a freqüência de desprendimento quando em comparação o caso 2D. Para a realização das simulação numéricas foi utilizado um esquema compacto de diferenças finitas em associação com um método de representação virtual de obstáculos imersos em escoamentos, que permite a representação de geometrias complexas sobre uma malha cartesiana.

Transição à turbulência na camada de mistura estavelmente estratificada utilizando simulação numérica direta e simulação de grandes escalas

A transição à turbulência em uma camada de mistura estavelmente estratificada é de grande interesse para uma variedade de problemas geofísicos e de engenharia. Esta transição é controlada pela competição entre o cisalhamento do escoamento de base e as forças de empuxo, devido à estratificação em densidade do ambiente. Os efeitos do empuxo atuam no escoamento reduzindo a taxa de crescimento das perturbações e retardando a transição à turbulência, enquanto o cisalhamento fornece energia cinética ao escoamento. O presente trabalho investiga a natureza da transição à turbulência em uma camada de mistura temporal estavelmente estratificada através de Simulação Numérica Direta (DNS) e Simulação de Grandes Escalas (LES). O propósito da investigação é analisar o efeito da estratificação estável no desenvolvimento da instabilidade de Kelvin-Helmholtz (K-H) e na formação dos vórtices longitudinais, que se formam após a saturação dos turbilhões primários de K-H. Além deste propósito, é examinado, utilizando de DNS e LES, o desenvolvimento das instabilidades secundárias de K-H na camada baroclínica. Os testes numéricos tridimensionais são realizados com diferentes tipos de condições iniciais para a flutuação de velocidade transversal, enquanto uma condição forçada é usada para as outras duas componentes de flutuação de velocidade. Em particular, o efeito do comprimento transversal do domínio de cálculo é testado empregando diferentes comprimentos, enquanto são usadas as mesmas dimensões para a direção longitudinal e vertical. As simulações bidimensionais mostram que o aumento da estratificação inibe o processo de emparelhamento, reduz a troca de energia entre os turbilhões de K-H e o escoamento, atenua a instabilidade de K-H e diminui o fluxo vertical de massa. A instabilidade secundária do tipo K-H é identtificada na camada baroclínica para Re ¸ 500 quando há o processo de emparelhamento dos vórtices simulados. Na simulação a Re = 500, a instabilidade secundária de K-H aparece tanto para Ri = 0.07 (fraca estratificação) como para Ri = 0.167 (forte estratificação). Os resultados tridimensionais demonstram que os vórtices longitudinais são claramente formados na camada a Ri = 0. Por outro lado, nos casos estratificados os vórtices são enfraquecidos, devido ao gradiente longitudinal de densidade, que diminui a vorticidade nos turbilhões de K-H enquanto aumenta na região entre eles.

Simulação sequencial na interpolação dos dados de entrada ou saída do modelo de lixiviação do software Araquá

Pós-graduação em Agronomia (Irrigação e Drenagem) - FCA

Métodos geoestatísticos na caracterização espacial de óxidos de ferro em diferentes pedoformas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Rede de modelos termodinâmicos locais

A simulação computacional de processos é uma ferramenta valiosa em diversas etapas da operação de uma planta química, tais como o projeto, a análise operacional, o estudo da estratégia de controle e a otimização, entre outras. Um programa simulador requer, freqüentemente, o cálculo das propriedades termofísicas que fazem parte das equações correspondentes aos balanços de massa, energia e quantidade de movimento. Algumas dessas quantidades, típicas de operações baseadas no equilíbrio de fases (tais como a destilação fracionada e a separação flash), são as fugacidades e entalpias das fases líquida e vapor em equilíbrio. Contudo, o uso de modelos e correlações cada vez mais sofisticadas para representar acuradamente as propriedades termofísicas dos sistemas reais fez com que grande parte do tempo gasto numa simulação de processos seja devida à avaliação destas propriedades. Muitas estratégias têm sido propostas na literaturas a fim de se reduzir a carga computacional envolvida na simulação, estática e dinâmica, dos problemas de Engenharia Química. Esta diminuição do tempo de processamento pode muitas vezes ser determinante na aplicação ou não da simulação de processos na prática industrial, como, por exemplo, no controle automático. Esta dissertação aborda uma das alternativas para a redução do tempo computacional gasto na simulação de processos: a aproximação das funções que descrevem as propriedades termofísicas através de funções mais simples, porém de fácil avaliação. Estas funções, a fim de se garantir uma aproximação adequada, devem ser corrigidas ou atualizadas de alguma forma, visto que se tratam de modelos estritamente válidos em uma pequena região do espaço das variáveis, sendo, por isto, chamados de modelos locais. Partindo-se do estado atual desta técnica, é proposta nesta dissertação uma nova metodologia para se efetuar esta aproximação, através da representação global da propriedade mediante múltiplas funções simples, constituindo, desse modo, uma rede de modelos locais. Algumas possibilidades para a geração efetiva destas redes são também discutidas, e o resultado da metodologia é testado em alguns problemas típicos de separação por equilíbrio de fases.

Simulação espacial da perda de solo e nutrientes por erosão em um argissolo sob cultivo de cana-de-açúcar

Pós-graduação em Agronomia (Ciência do Solo) - FCAV

Métodos geoestatísticos na modelagem espacial do diâmetro médio do cristal da goethita

Uma das necessidades da agricultura de precisão é avaliar a qualidade dos mapas dos atributos dos solos. Neste sentido, o presente trabalho objetivou avaliar o desempenho dos métodos geoestatísticos: krigagem ordinária e simulação sequencial gaussiana na predição espacial do diâmetro médio do cristal da goethita com 121 pontos amostrados em uma malha de 1 ha com espaçamentos regulares de 10 em 10 m. Após a análise textural e da concentração dos óxidos de ferro, calcularam-se os valores do diâmetro médio do cristal da goethita os quais foram analisados pela estatística descritiva e geoestatística; em seguida, foram utilizadas a krigagem ordinária e a simulação sequencial gaussiana. Com os resultados avaliou-se qual foi o método mais fiel para reproduzir as estatísticas, a função de densidade de probabilidade acumulada condicional e a estatística epsilon εy da amostra. As estimativas E-Type foram semelhantes à krigagem ordinária devido à minimização da variância. No entanto, a krigagem deixa de apresentar, em locais específicos, o grau de cristalinidade da goethita enquanto o mapa E-Type indicou que a simulação sequencial gaussiana deve ser utilizada ao invés de mapas de krigagem. Os mapas E-type devem ser preferíveis por apresentar melhor desempenho na modelagem.

Simulação geoestatística e classificação numérica para o delineamento de unidade de mapeamento de atributos e erosão do solo

Pós-graduação em Agronomia (Ciência do Solo) - FCAV

Uncertainties in the prediction of spatial variability of soil CO2 emissions and related properties

A emissão de CO2 do solo apresenta alta variabilidade espacial, devido à grande dependência espacial observada nas propriedades do solo que a influenciam. Neste estudo, objetivou-se: caracterizar e relacionar a variabilidade espacial da respiração do solo e propriedades relacionadas; avaliar a acurácia dos resultados fornecidos pelo método da krigagem ordinária e simulação sequencial gaussiana; e avaliar a incerteza na predição da variabilidade espacial da emissão de CO2 do solo e demais propriedades utilizando a simulação sequencial gaussiana. O estudo foi conduzido em uma malha amostral irregular com 141 pontos, instalada sobre a cultura de cana-de-açúcar. Nesses pontos foram avaliados a emissão de CO2 do solo, a temperatura do solo, a porosidade livre de água, o teor de matéria orgânica e a densidade do solo. Todas as variáveis apresentaram estrutura de dependência espacial. A emissão de CO2 do solo mostrou correlações positivas com a matéria orgânica (r = 0,25, p < 0,05) e a porosidade livre de água (r = 0,27, p <0,01) e negativa com a densidade do solo (r = -0,41, p < 0,01). No entanto, quando os valores estimados espacialmente (N=8833) são considerados, a porosidade livre de água passa a ser a principal variável responsável pelas características espaciais da respiração do solo, apresentando correlação de 0,26 (p < 0,01). As simulações individuais propiciaram, para todas as variáveis analisadas, melhor reprodução das funções de distribuição acumuladas e dos variogramas, em comparação à krigagem e estimativa E-type. As maiores incertezas na predição da emissão de CO2 estiveram associadas às regiões da área estudada com maiores valores observados e estimados, produzindo estimativas, ao longo do período estudado, de 0,18 a 1,85 t CO2 ha-1, dependendo dos diferentes cenários simulados. O conhecimento das incertezas gerado por meio dos diferentes cenários de estimativa pode ser incluído em inventários de gases do efeito estufa, resultando em estimativas mais conservadoras do potencial de emissão desses gases.

Incerteza espacial da emissão de 'CO IND.2' do solo em dois sistemas de manejo da palhada da cana-de-açúcar

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Incertezas na estimativa da variabilidade espacial da emissão de CO2 do solo e propriedades edáficas em área de cana crua

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)