Geostatistical methods for definition of geological model of a pre-salt

A presente tese apresenta os resultados da pesquisa do projeto: Métodos Geoestatísticos para a Definição de um Modelo Geológico de um Reservatório do Pré-sal Brasileiro. O objetivo principal do projeto é identificar as propriedades do reservatório e compreender os fatores que afetam a alta heterogeneidade do reservatório como litologia, mineralogia e propriedades do reservatório, conhecer os parâmetros responsáveis ??pelo aumento abrupto da permeabilidade em determinados intervalos, finalmente representando essas características em um modelo geológico, identificando os algoritmos mais adequados para a representação das fácies e propriedades incluindo a representação de fraturas e vugs. Para atingir o objetivo principal, o estudo foi dividido em três partes. A primeira seção trata de compreender as relações entre a composição mineralógica, fácies, as características do reservatório e a produção. Nesta fase, um artigo foi publicado na revista Geoenergy Science and Engineering. “Utilizando Inteligência Artificial Integrada para Caracterizar Mineralogia e Fácies em um Reservatório Carbonatado do Pré-sal, Bacia de Santos, Brasil, usando Núcleos, Registros Wireline e Avaliação Petrofísica Multimineral” Este trabalho nos levou à conclusão de que as melhores zonas dentro do reservatório estão associadas à fácies arbustos-esferulitos, alto volume de dolomitas com presença abundante de porosidades secundárias (vugs). O segundo artigo refere-se à identificação de zonas com permeabilidade e produção extremamente altas, que coincidem com características identificadas no primeiro artigo (altas quantidades de vugs, concentrações altas de dolomita e representadas por fácies como arbustos e esferulites) e podem atender aos requisitos das zonas super K. As abordagens tradicionais de determinação da permeabilidade (amostras de testemunhos, permeabilidade por RMN ou regressões) são incapazes de reproduzir os altos valores de permeabilidade relatados em testes de produção do tipo Drill-stem test DST, ou os altos fluxos de produção distribuidores nos perfis de produção PLTs. Esta pesquisa fornece uma nova metodologia que inclui dados estáticos e contabilização da porosidade das vugs e seus impactos no aumento da permeabilidade, sendo calibrados com dados dinâmicos como testes DSTs, PLT ou Índice de Produtividade (J). além de um método inovador para estimar a probabilidade de encontrar zonas super K. A última parte corresponde à construção do modelo geológico que possa representar as heterogeneidades do reservatório. Para modelagem de fácies avaliaremos dois dos métodos tradicionais Sequential Indicator Simulation SIS e Truncated Gaussinan Simulation TGS, e incorporaremos o uso de Multi-Point Statistics MPS, para corrigir as deficiências apresentadas por métodos baseados em estatísticas de dois pontos (variogramas) como SIS e TGS que são incapazes de reproduzir a continuidade de corpos geológicos, assim como não são capazes de reproduzir corpos geométricos complexos.Pelo contrário, o MPS é capaz de solucionar as fragilidades dos métodos baseados em variogramas, porém é apoiado por imagens de treinamento (TI) que são representações de modelos geológicos conceituais. Os resultados da modelagem de fácies mostraram que o método MPS mostra melhores resultados na modelagem de fácies em reservas complexas como os carbonatos do pré-sal, porem a criação da imagem de treino representa um verdadeiro desafio

A posteriori error estimation for the multiscale hybrid-mixed finite element method

Estimação de erro a posteriori para o método dos elementos finitos misto-híbrido multiescala para problemas de Darcy, MHM-H(div)-E’gama’ , é apresentada. O método MHM-H(div)-E’gama’ adota espaços de elementos finitos de duas escalas: discretizações refinadas são adotadas no interior de subregiões poligonais, mas aproximações do fluxo são restritas sobre as interfaces da malha por um espaço de traços normais coarse. Para estabilidade, aproximações da pressão e do fluxo são compatíveis em termos do divergente, satisfazendo o diagrama de De Rham. A estimação de erro se baseia na reconstrução do potencial, método popular neste tipo de análise no contexto de métodos mistos. Experimentos numéricos são apresentados para ilustrar a eficiência da metodologia proposta

Modeling electrical submersible pump system under multiphase flows using bond graphs

Bombas centrífugas submersas (ESPs) são amplamente utilizadas em setores industriais que necessitam de altas vazões e ganho de pressão. Na indústria de petróleo e gás, as ESPs frequentemente operam sob condições de escoamento multifásico, tais como emulsões de óleo-água. O comportamento não-Newtoniano das emulsões pode induzir um comportamento dinâmico e instável do sistema. Contudo, a literatura existente foca no comportamento em estado estacionário destes sistemas. O objetivo principal deste trabalho é desenvolver e validar um modelo dinâmico para ESPs utilizando bond-graphs, conduzir análises de identificabilidade e estimar parâmetros do sistema com um conjunto limitado de sensores. A metodologia adotada engloba a coleta de dados experimentais em condições estacionárias e dinâmicas, modelagem por meio de bond-graphs e o emprego de Redes Neurais Informadas por Física (PINNs). O modelo da ESP foi formulado utilizando uma biblioteca desenvolvida, neste trabalho, em Julia. Foram conduzidas análises de identificabilidade, tanto estrutural quanto prática do modelo. As PINNs foram empregadas para abordar o problema inverso e foram avaliadas em dados simulados e experimentais. A validação do modelo foi realizada através da comparação dos dados experimentais com simulações numéricas. A pesquisa resultou em um modelo dinâmico para ESPs baseado em bond-graph que incorpora os subsistemas mecânicos e hidráulicos, resultando em um conjunto de equações diferenciais ordinárias (ODEs). Em condições estacionárias, o modelo demonstrou um elevado coeficiente de determinação e um erro relativamente baixo. Para cenários dinâmicos, ajustes finos nos parâmetros resultaram em melhoria na precisão do modelo, especialmente em relação às pressões, embora algumas discrepâncias nos picos de pressão tenham sido observadas. A análise de identificabilidade estrutural local identificou 12 parâmetros que são determináveis, mas a identificabilidade prática foi alcançada apenas com oito. As PINNs mostraram eficácia na estimativa de parâmetros e estados, sobretudo em cenários de baixa fração de água, mas apresentaram limitações em altas frações de água e em presença de ruído nos dados. O modelo dinâmico, em conjunto com as PINNs, demonstra potencial para aplicações em controle, monitoramento, detecção de falhas e otimização. Embora o modelo capture a dinâmica do sistema, limitações surgem devido às hipóteses feitas quanto à viscosidade da emulsão e ao acoplamento bomba-tubo. Para trabalhos futuros é sugerido aprimorar na modelagem da propagação de ondas de pressão no sistema. Além disso, a adoção de algoritmos mais robustos para PINNs, que poderia melhorar a estimativa dos parâmetros desconhecidos. Este trabalho estabelece uma base para a modelagem de sistemas ESP mais complexos, com potencial de aplicação além do setor de petróleo e gás

Generation of 3D geological model of the pre-salt with karst structures

A sequência carbonática do pré-sal é composta por fácies de camadas finas com heterogeneidades verticais e laterais, principalmente devido à grande variedade de alterações pós-deposicionais que resultam em comportamentos petrofísicos secundários na formação. À medida que a indústria continua a explorar essas jogadas complicadas, é necessário um melhor entendimento da distribuição de características cársticas, incluindo fraturas e como elas afetam a distribuição de propriedades petrofísicas que afetam o desenvolvimento de ativos petrolíferos e projetos de sequestro de carbono. Dados sísmicos, perfis de imagens de poços e perfis de poços convencionais foram usados para investigar a distribuição e características de feições paleocarsticas nos carbonatos aptianos da Formação Barra Velha em uma área piloto da Bacia de Santos, Brasil. Vários atributos sísmicos foram usados para melhorar os detalhes dos dados sísmicos e destacar os principais parâmetros sísmicos, incluindo pontos brilhantes, deformação e geometria de estratos, continuidade de eventos sísmicos e padrões de falha. Também discutimos a ocorrência de permeabilidades secundárias de carstificação e mostramos como elas afetam as propriedades petrofísicas da Formação Barra Velha usando interpretação sísmica integrada, avaliação de formação de wireline e análise de imagem de poço a partir do qual redes de fraturas discretas foram desenvolvidas para explicar o efeito da fraturas e carstes no modelo geológico. Este trabalho fornece novas informações sobre a relação entre as zonas carstificadas e os domínios estruturalmente danificados na sucessão do pré-sal da Bacia de Santos, Brasil. Micro e Mega vugs ocorrem amplamente nos carbonatos da Formação Barra Velha, e as estruturas cársticas parecem ser controladas por fraturas e falhas. As feições cársticas da Formação Barra Velha consistem em vugs que variam de cerca de 3mm a 3cm, fraturas com inclinação, orientação e abertura variáveis, e cavernas/sumidouros alargados preenchidos com materiais brechados em alguns casos. As curvas de proporção vertical 1D foram adotadas como um melhor método para capturar as fácies heterogêneas em grande escala, fornecendo restrições para o modelo de fácies verticalmente com controle limitado na distribuição lateral, porém o número e a distribuição dos poços podem reduzir positivamente as áreas de incertezas na área de estudo . O modelo petrofísico foi gerado usando variogramas avançados com viés de tendência de fácies, enquanto a análise de dados geoestatísticos ajudou a garantir a consistência do modelo com os dados disponíveis. Fraturas e falhas controlam significativamente a distribuição do excesso de porosidades e permeabilidades como resultado da carstificação ao longo desses condutos. A rede de fraturas discretas consistindo predominantemente de fraturas NE-SW e NW-SE ajudou no desenvolvimento do modelo de porosidade-permeabilidade dupla que captura as anisotropias de fluxo matricial e não matricial. Este trabalho apresenta a influência de carstificações e fraturas nas propriedades petrofísicas dos montículos carbonáticos lacustres da Formação Barra Velha e também contribui para o entendimento da distribuição dessas carstificações nessas unidades, úteis em futuras utilidades subsuperficiais e podem influenciar na tomada de decisões quanto aos futuros projetos de exploração/desenvolvimento de hidrocarbonetos e sequestro de carbono na margem continental SE brasileira. As limitações deste estudo estão relacionadas com as dificuldades associadas à integração de múltiplos conjuntos de dados geológicos e geofísicos com várias escalas, no entanto, as imagens de poços que estão em escala sub-sísmica apresentam com maior confiança a ocorrência das feições cársticas

Characterization and modeling of fractures using borehole image logs and seismic attributes for the pre-salt section of the Santos basin, Brazil

A Bacia offshore de Santos é uma bacia prolífica com grande potencial para exploração de hidrocarbonetos, em parte devido à sua evolução tectônica e estratigrafia associada a carbonatos lacustres descobertos nas águas profundas brasileiras embaixo do sal. A alta permeabilidade nesta bacia é provavelmente controlada por altas intensidades de fratura e dissolução de carbonato, portanto, prever o comportamento de fraturas naturais nesses depósitos ajudará a entender seus efeitos no fluxo de fluido e na produtividade do reservatório. Este estudo se concentra na construção de um modelo de rede de fratura discreta com base na compreensão da intensidade, orientação e distribuição espacial das fraturas, isso foi realizado por meio do uso de ferramentas como registros de imagem de poço de 4 poços e sísmica 3-D. A partir dessas ferramentas foi possível identificar os principais horizontes (topo e base do play do pré-sal), falhas regionais, interpretação das fraturas e do estado de tensão in situ no play do pré-sal da Bacia de Santos. Neste trabalho, descobrimos que, de acordo com os breakouts e as fraturas induzidas por perfuração, a tensão horizontal máxima (SHmax) tem um azimute NE-SW, que comparado com a orientação e regime de estruturas de grande escala sísmica no reservatório, é possível destacar que a área de estudo é influenciada por um padrão normal de falha e configuração de tensão extensional. Além disso, de acordo com as fraturas naturais, identificamos quatro sets, os sets 1 e 2 seguem uma direção paralela a SHmax NE-SW e os sets 3 e 4 seguem uma direção paralela a Shmin. Em relação à caracterização das redes de fraturas, distribuímos os atributos de fratura no grid geológica de acordo com parâmetros como densidade de fraturas, abertura de fratura, comprimento de fratura e orientação de fratura, esses parâmetros permitiram quantificar a porosidade e permeabilidade da fratura, destacando assim regiões com as melhores respostas de permeabilidade de fratura

Carbonate reservoir characterization with integrated rock, special well logs, and machine learning methods : a study case in a Pre-Salt field from the Santos Basin

O desenvolvimento da produção de petróleo nos reservatórios do Pré-Sal requer suporte de modelos geológicos e petrofísicos consistentes. O principal objetivo desta pesquisa é desenvolver fluxos de trabalho de caracterização de reservatórios, combinando geologia, petrofísica e aprendizado de máquina. O intervalo estudado, pertencendo à Formação Barra Velha, apresenta uma sucessão de fácies in-situ, composta por shrubstones, mudstones laminados e spherulitestones, e fácies retrabalhadas compostas por grainstones e rudstones intraclásticos. A primeira parte do estudo integra descrição de testemunho, petrografia e perfis de imagem acústica para segmentar o intervalo em Petrofácies de Reservatório. Posteriormente avalia-se a qualidade de reservatório com perfis de ressonância magnética nuclear (RMN) e Micro-CT. A melhor qualidade de reservatório pertence a rudstones com baixa cimentação e grainstones com moderada cimentação por dolomita/quartzo, com valores de porosidade entre 15-20% e permeabilidades entre 100 e 500 mD. As melhores fácies in-situ pertencem aos shrubstones com cimentação por dolomita moderada, com porosidade entre 12-17% e permeabilidade em torno de 100 mD. O método integrado mostrou-se eficaz na transferência de informações, em escala de testemunho para perfil do poço, além de permitir ao intérprete fazer inferências fora do intervalo testemunhado, respeitando o contexto geológico. A segunda parte do estudo aborda a estimativa imprecisa da porosidade em intervalos dominados por porosidade secundária, aplicando aprendizado de máquina supervisionado. Integra-se dados de perfis convencionais, RMN e um índice de porosidade vugular extraído de perfis de imagem de acústica, e valida-se o modelo com análises petrofísica básica. Foram testados os algoritmos Random Forest, XGBoost e Support Vector Regression, com métricas de R2 de 0,835; 0,836; 0,800 e raiz do erro quadrático médio (REQM) de 1,75; 1,68; 1,82; respectivamente. Todos os três modelos superaram a estimativa original do perfil RMN de R2 de 0,64 e REQM de 2,78. Limitações metodológicas podem surgir, especialmente relacionadas a artefatos nos perfis de imagem acústicas, que precisam ser mitigados antes do treinamento do modelo. A incorporação de métodos de aprendizado de máquina à caracterização de reservatórios pode aprimorar a precisão na estimativa de parâmetros de reservatório, além do potencial para redução no tempo de análise o custo com amostragem. Em última análise, a integração destes métodos com petrofísica e geologia aprimora o fluxo de trabalho no desenvolvimento de reservatórios complexos

Multi-objective production strategy optimization methodology under uncertainties

A otimização do desenvolvimento de campos de petróleo desempenha um papel fundamental na maximização do potencial dos reservatórios de hidrocarbonetos. A tomada de decisão nesse domínio complexo pode se basear em diversas funções-objetivo. O aumento da complexidade dos sistemas de reservatórios e as incertezas inerentes em fatores geológicos e econômicos têm sugerido o uso da otimização multiobjetivo (MOO) para melhorar os processos decisórios. Os potenciais benefícios da MOO incluem uma tomada de decisão mais robusta, que considera os compromissos entre diferentes objetivos, e que pode acomodar melhor as incertezas, resultando em estratégias de produção mais resilientes e adaptáveis. O objetivo principal deste estudo é projetar um framework robusto multiobjetivo para otimizar estratégias de produção considerando incertezas nos parâmetros do reservatório e nas condições econômicas, com foco específico em campos de petróleo do pré-sal, aplicando o reservatório de referência UNISIM-II-D. Para enfrentar os desafios decorrentes da aplicação da MOO, foram realizados três estudos principais para desenvolver e avaliar frameworks que contemplam esses objetivos complexos. O primeiro estudo apresenta uma revisão abrangente da evolução dos algoritmos de MOO, destacando sua crescente relevância na otimização do desenvolvimento de campos. Ele categoriza os algoritmos em métodos a priori e a posteriori, abordando desafios como o custo computacional elevado das simulações de reservatórios sob incerteza. O segundo estudo investiga o comportamento de diferentes funções-objetivo dentro da otimização do desenvolvimento de campos, com foco principal no valor monetário esperado (VME) como o principal critério. Essa análise busca identificar objetivos conflitantes, formando a base para a última etapa do estudo. Os resultados destacam uma forte correlação entre o VME e a produção cumulativa de óleo (PCO), demonstrando que otimizar o VME pode impulsionar melhorias na PCO e no fator de recuperação. Além disso, a adoção de estratégias que priorizam objetivos ambientais (como a redução da produção cumulativa de água e gás), mesmo que impliquem em uma leve queda no VME, pode proporcionar benefícios substanciais, incluindo redução das emissões de gases de efeito estufa, menores custos de tratamento de água e maior vida útil do reservatório. O estudo mostra que confiar exclusivamente no VME na otimização do desenvolvimento de campos traz riscos inerentes, pois nem todos os modelos representativos (MRs) apresentam aumento no valor presente líquido (VPL). Em particular, no caso mais pessimista, o VPL não demonstra melhoria quando o VME é maximizado por meio de uma otimização de objetivo único. A terceira etapa foi iniciada após a identificação dos objetivos conflitantes na segunda etapa para abordar e equilibrar sistematicamente esses conflitos no planejamento do desenvolvimento do campo. Nesta etapa, a Otimização Multiobjetivo (MOO) foi realizada utilizando o algoritmo genético de classificação não dominada II (NSGA-II) para equilibrar simultaneamente o EMV com o VPL sob incerteza econômica (VPLeco) para o RM pessimista, gerando assim um conjunto diversificado de soluções Pareto-ótimas que ilustram efetivamente as compensações entre diferentes estratégias de produção possíveis, considerando incertezas geológicas e econômicas. Essa abordagem não apenas identifica um conjunto mais amplo de estratégias ótimas, mas também oferece uma compreensão mais refinada dos compromissos, melhorando o desempenho geral ao capturar sinergias entre os objetivos. Com o framework proposto, foi alcançada uma melhoria de 3% no VME e um aumento de 28% no VPLeco do MR pessimista em comparação com a otimização de objetivo único do VME. Esses resultados destacam a força e a robustez do framework em proporcionar decisões informadas e equilibradas que guiem o desenvolvimento de reservatórios rumo a melhores resultados econômicos e geológicos

Aplicações sintéticas dos adutos de Morita-Baylis-Hillman em reações mediadas por paládio

Esse trabalho descreveu o uso de adutos de Morita-Baylis-Hillman como substratos para reações mediadas por paládio. Inicialmente, desenvolvemos uma estratégia eficiente para preparar a-benzilcetoesteres a partir de adutos de Morita-Baylis-Hillman. A reação de Mizoroki-Heck foi empregada nessa estratégia e esses compostos foram obtidos em elevados rendimentos e seletividades. Além disso, os paladaciclos de Nájera foram usados com sucesso pela primeira vez como catalisadores nessa reação. Em seguida, os derivados a-(4-hidroxibenzil)-b-cetoesteres obtidos dessa reação foram tratadas com um reagente de iodo hipervalente para fornecer novas espirocicloexadienonas em bons rendimentos (21- 60%). Inúmeros grupos funcionais mostraram ser toleráveis. As atividades biológicas contra bactérias Gram-positivas mostraram que esses compostos são mais ativos que a meticilina. A segunda parte desse trabalho reportou o uso de adutos de MBH em reações de ciclocarbonilação para sintetizar 3-alquenil-ftalídeos. Nesse método, o uso de grupos de proteção no grupo hidroxila não era necessário e as seletividades observadas foram acima de 95:05 E/Z. Quando os derivados de adutos de MBH-nitro o-halogenados foram utilizados, obtivemos indanonas e, neste caso, propomos que uma reação de Heck intramolecular e favorecida. Na sequência, os adutos de MBH foram tratados com o sistema Pd/norborneno para formar compostos policíclicos. Quando os ácidos borônicos foram utilizados, isolamos compostos em baixo rendimento a partir de um processo de arilação tandem. Na última parte dessa tese os adutos de MBH também foram empregados para preparar 4 fragrâncias. Essa estratégia direta e simples originou esses compostos em bom rendimento global (34-41%)

Beyond the Wax Concentration; How Micro/Macro Wax Ratios Change Thermal, Rheological, and Morphological Properties of the Waxy System

Knowledge of the rheological and thermal properties of wax oils is central to understanding their behavior and their impact under the diverse operational conditions to which these oils are subjected in the O&G industry. A better understanding of the compositional effects on the behavior of wax deposits is essential for the improvement of prediction and mitigation techniques for this type of phenomenon. In this study, we evaluated different model systems based on micro/macro wax and how these combinations affect the thermal and rheological properties. As expected, combinations richer in micro wax showed higher wax appearance temperature and gelation temperature (Tgel), but gel strength is inversely proportional to micro wax content, as evidenced by rheological tests. The morphology and aspect ratio of the crystals are parameters that directly impact the rheological properties of the formulations; crystals with a high aspect ratio and needle morphology tend to form stronger gels due to greater interaction between crystals. The rheological and morphological tests suggest that aggregates of micro wax crystals can weaken the three-dimensional network of macro wax crystals and consequently reduce the yield stress. The results suggest that the addition of micro wax can minimize the pressure required for oil field flow restart.