京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第19695号 / 工博第4150号 / 新制||工||1640(附属図書館) / 32731 / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 小池 克明, 教授 高岡 昌輝, 准教授 村田 澄彦 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

Although carbonates hold more than 60 percent of the world's oil reserves, they, nevertheless, exhibit much lower average recovery factor values than terrigenous sandstone reservoirs. Thus, utilization of advanced enhanced oil recovery (EOR) techniques such as high pressure CO2 injection may normally be required to recover oil in place in carbonate reservoirs. This study addresses how different rock types can influence the seismic monitoring of CO2 sequestration in carbonates. This research utilizes an elastic parameter, defined in a rock physics model of poroelasticity and so-­called as the frame flexibility factor, to successfully quantify the carbonate pore types in core samples available from the Great Bahama Bank (GBB). This study shows that for carbonate samples of a given porosity the lower the frame flexibility factors the higher is the sonic wave velocity. Generally, samples with frame flexibility values of <4 are either rocks with visible moldic pores or intraframe porosity; whereas, samples with frame flexibility values of >4 are rocks with intercrystalline and microporosity. Hence, different carbonate pore geometries can be quantitatively predicted using the elastic parameters capable of characterizing the porous media with a representation of their internal structure on the basis of the flexibility of the frame and pore connectivity. In this research, different fluid substitution scenarios of liquid and gaseous CO2 saturations are demonstrated to characterize the variations in velocity for carbonate-specific pore types. The results suggest that the elastic response of CO2 flooded rocks is mostly governed by pore pressure conditions and carbonate rock types. Ultrasonic P-­wave velocities in the liquid-­phase CO2 flooded samples show a marked decrease in the order of 0.6 to 16 percent. On the contrary, samples flooded with gaseous-­phase CO2 constitute an increase in P-­wave velocities for moldic and intraframe porosities, while establishing a significant decrease for samples with intercrystalline and micro-­porosities. Such velocity variations are explained by the stronger effect of density versus compressibility, accounting for the profound effect of pore geometries on the acoustic properties in carbonates. The theoretical results from this research could be a useful guide for interpreting the response of time-­lapse seismic monitoring of carbonate formations following CO2 injection at depth. In particular, an effective rock-­physics model can aid in better discrimination of the profound effects of different pore geometries on seismic monitoring of CO2 sequestration in carbonates.

[pt] O conhecimento do que acontece no reservatório em produção a partir de variações temporais dos atributos sísmicos devido aos processos dinâmicos vem atingindo um valor crescente na indústria do petróleo, especialmente em arenitos. Este processo possui vários desafios, focados em grande parte a desvendar a superposição dos diferentes efeitos provocados pelas mudanças do reservatório nos dados sísmicos. As propriedades sísmicas são afetadas de maneira complexa por vários fatores, sendo a saturação um dos mais importantes, principalmente em rochas porosas como o arenito. Esta propriedade influencia no módulo elástico da rocha e sua resposta sísmica e, ao mesmo tempo, introduz dispersão da velocidade (variação da velocidade com a freqüência). A transição de fluido efetivo (distribuição homogênea e menores velocidades) para fluido com distribuição heterogênea (e maiores velocidades) estabelece um mecanismo de dispersão presente para freqüências sísmicas in situ, especialmente no arenito. O método mais utilizado para aplicar a técnica de substituição de fluidos se baseia na teoria de Gassmann (1951), que considera o meio poroso estático (estado de isostress), onde o fluido não é afetado pela perturbação da onda. No entanto, pesquisas mostram que as velocidades acústicas em rochas saturadas de fluido dependem da freqüência, do tipo de fluido e sua distribuição no meio poroso, viscosidade e outras propriedades que tornam as ondas dispersivas. Neste trabalho são realizadas simulações de fluxo de reservatórios, transformações de física de rochas, upscaling e modelagem sísmica em cenários de injeção de gás com o objetivo de esclarecer a importância de levar em conta a dispersão da velocidade na análise time-lapse. Para isso, são analisados para cada modelo mapas de saturação, velocidade, impedância e sismogramas sintéticos (seções de contraste) calculados com as teorias de substituição Gassmann (1951) e Mavko E Jizba (1991). Os resultados mostram que a resposta sísmica pode ter um incremento de até 15 por cento quando a dispersão devida ao fluxo local é considerada. Porosidade e tortuosidade são parâmetros essenciais que influenciam de maneira diferente na resposta sísmica. / [en] The evaluation of reservoir dynamics during production through time-lapse interpretation has reached a substantial importance in the petroleum industry, mainly in sandstones. This evaluation presents many challenges, mainly concerned to unmask the overlapping of different effects in seismic data due to reservoir changes. Several factors affect seismic properties and saturation is one of the most important. This property influences the rock bulk modulus and seismic response and also causes a velocity dependence on the frequency. This phenomenon is known as velocity dispersion. Furthermore, the transition from effective homogeneous fluid to heterogeneous saturation represents a dispersion mechanism that appears for seismic frequencies in situ in sandstones. The most commonly method used to perform the fluid substitution technique is based in Gassmann theory (1951). This approach considers a static porous media (isostress condition), where fluid is not affected by wave propagation. However, it is well known that acoustic velocities in fluid saturated rocks depends on frequency, according to fluid type and distribution on porous media, viscosity, and others properties that become waves dispersive. In this work reservoir flow-simulation, rock physics transformations, upscaling and seismic modeling were performed in gas injection scenarios. Synthetic seismograms and some contrast sections were generated using Gassmann (1951) and Mavko & Jizba (1991) substitution theories. The goal is to clarify the relevance of considering velocity dispersion on time-lapse seismic analyzing possible differences in the seismic parameters. Results show that seismic response could increase in 15% when squirt flow dispersion is considered. Porosity and tortuosity are essential parameters to analyze seismic response.

Under current technical levels and industrial systems, geological carbon sequestration (GCS) is a viable solution to maintain and further reduce carbon dioxide (CO2) concentration and ensure energy security simultaneously. The pre-injection formation characterization and post-injection CO2 monitoring, verification, and accounting (MVA) are two critical and challenging tasks to guarantee the sequestration effect. The tasks can be accomplished using core analyses and well-logging technologies, which complement each other to produce the most accurate and sufficient subsurface information for pore-scale and reservoir-scale studies. In recent years, the unprecedented data sources, increasing computational capability, and the developments of machine learning (ML) and deep learning (DL) algorithms provide novel perspectives for expanding the knowledge from data, which can capture highly complex nonlinear relationships between multivariate inputs and outputs. This work applied ML and DL methods to GCS-related studies at pore and reservoir scales, including digital rock physics (DRP) and the well-logging data interpretation and analysis. DRP provides cost-saving and practical core analysis methods, combining high-resolution imaging techniques, such as the three-dimensional (3D) X-ray computed tomography (CT) scanning, with advanced numerical simulations. Image segmentation is a crucial step of the DRP framework, affecting the accuracy of the following analyses and simulations. We proposed a DL-based workflow for boundary and small target segmentation in digital rock images, which aims to overcome the main challenge in X-ray CT image segmentation, partial volume blurring (PVB). The training data and the model architecture are critical factors affecting the performance of supervised learning models. We employed the entropy-based-masking indicator kriging (IK-EBM) to generate high-quality training data. The performance of IK-EBM on segmentation affected by PVB was compared with some commonly used image segmentation methods on the synthetic data with known ground truth. We then trained and tested the UNet++ model with nested architecture and redesigned skip connections. The evaluation metrics include the pixel-wise (i.e. F1 score, boundary-scaled accuracy, and pixel-by-pixel comparison) and physics-based (porosity, permeability, and CO2 blob curvature distributions) accuracies. We also visualized the feature maps and tested the model generalizations. Contact angle (CA) distribution quantifies the rock surface wettability, which regulates the multiphase behaviors in the porous media. We developed a DL-based CA measurement workflow by integrating an unsupervised learning pipeline for image segmentation and an open-source CA measurement tool. The image segmentation pipeline includes the model training of a CNN-based unsupervised DL model, which is constrained by feature similarity and spatial continuity. In addition, the over-segmentation strategy was adopted for model training, and the post-processing was implemented to cluster the model output to the user-desired target. The performance of the proposed pipeline was evaluated using synthetic data with known ground truth regarding the pixel-wise and physics-based evaluation metrics. The resulting CA measurements with the segmentation results as input data were validated using manual CA measurements. The GCS projects in the Illinois Basin are the first large-scale injection into saline aquifers and employed the latest pulsed neutron tool, the pulsed neutron eXtreme (PNX), to monitor the injected CO2 saturation. The well-logging data provide valuable references for the formation evaluation and CO2 monitoring in GCS in saline aquifers at the reservoir scale. In addition, data-driven models based on supervised ML and DL algorithms provide a novel perspective for well-logging data analysis and interpretation. We applied two commonly used ML and DL algorithms, support vector machine regression (SVR) and artificial neural network (ANN), to the well-logging dataset from GCS projects in the Illinois Basin. The dataset includes the conventional well-logging data for mineralogy and porosity interpretation and PNX data for CO2 saturation estimation. The model performance was evaluated using the root mean square error (RMSE) and R2 score between model-predicted and true values. The results showed that all the ML and DL models achieved excellent accuracies and high efficiency. In addition, we ranked the feature importance of PNX data in the CO2 saturation estimation models using the permutation importance algorithm, and the formation sigma, pressure, and temperature are the three most significant factors in CO2 saturation estimation models. The major challenge for the CO2 storage field projects is the large-scale real-time data processing, including the pore-scale core and reservoir-scale well-logging data. Compared with the traditional data processing methods, ML and DL methods achieved accuracy and efficiency simultaneously. This work developed ML and DL-based workflows and models for X-ray CT image segmentation and well-logging data interpretations based on the available datasets. The performance of data-driven surrogate models has been validated regarding comprehensive evaluation metrics. The findings fill the knowledge gap regarding formation evaluation and fluid behavior simulation across multiple scales, ensuring sequestration security and effect. In addition, the developed ML and DL workflows and models provide efficient and reliable tools for massive GCS-related data processing, which can be widely used in future GCS projects. / Doctor of Philosophy / Geological carbon sequestration (GCS) is the solution to ease the tension between the increasing carbon dioxide (CO2) concentrations in the atmosphere and the high dependence of human society on fossil energy. The sequestration requires the injection formation to have adequate storage capability, injectivity, and impermeable caprock overlain. Also, the injected CO2 plumes should be monitored in real-time to prevent any migration of CO2 to the surface. Therefore, pre-injection formation characterization and post-injection CO2 saturation monitoring are two critical and challenging tasks to guarantee the sequestration effect and security, which can be accomplished using the combination of pore-scale core analyses and reservoir-scale well-logging technologies. This work applied machine learning (ML) and deep learning (DL) methods to GCS-related studies across multiple spatial scales. We developed supervised and unsupervised DL-based workflows to segment the X-ray computed-tomography (CT) image of digital rocks for the pore-scale studies. Image segmentation is a crucial step in the digital rock physics (DRP) framework, and the following analyses and simulations are conducted on the segmented images. We also developed ML and DL models for well-logging data interpretation to analyze the mineralogy and estimate CO2 saturation. Compared with the traditional well-logging analysis methods, which are usually time-consuming and prior knowledge-dependent, the ML and DL methods achieved comparable accuracy and much shorter processing time. The performance of developed workflows and models was validated regarding comprehensive evaluation metrics, achieving excellent accuracies and high efficiency simultaneously. We are at the early stage of CO2 sequestration, and relevant knowledge and tools are inadequate. In addition, the main challenge of CO2 sequestration field projects is the large-scale and real-time data processing for fast decision-making. The findings of this dissertation fill the knowledge gap in GCS-related formation evaluation and fluid behavior simulations across multiple spatial scales. The developed ML and DL workflows provide efficient and reliable tools for massive data processing, which can be widely used in future GCS projects.

[pt] Desde a década de 80, abordagens data-driven têm sido utilizadas para identificação de fluidos e caracterização de reservatórios carbonáticos e siliciclásticos principalmente em relação à análise das amplitudes sísmicas. No entanto, técnicas aplicadas com sucesso para rochas siliciclásticas, como por exemplo: Análise AVO, inversões sísmicas e IDH (Indicadores Diretos de Hidrocarbonetos) revelaram não obter o mesmo êxito para reservatórios carbonáticos heterogêneos. Em contrapartida, diversos artigos mostram que fluxos de caracterização de reservatórios com modelos de física de rochas incorporados têm alcançado grande sucesso para obtenção de propriedades petrofísicas e atributos elásticos de ambas as rochas, utilizando sísmicas e well logs, em uma abordagem model-driven, focada nas características microestruturais do reservatório. Dessa forma, levando em consideração a importância de se utilizar modelos de física de rochas no escopo da caracterização de reservatórios, dois modelos de física de rochas - Xu e Payne e T-Matrix - foram aplicados, comparados e seus parâmetros foram estocasticamente avaliados e otimizados em um arcabouço Bayesiano. Através dessa abordagem, foi possível estimar, de uma forma confiável, os atributos elásticos de um reservatório carbonático (coquinas) levando em consideração diversos tipos de incertezas. Além disso, após a calibração e validação de ambos os modelos de física de rochas para diferentes poços, análises de sensibilidade foram realizadas para compreensão de forma quantitativa do comportamento dos atributos elásticos das coquinas em relação às alterações do conteúdo mineralógico, tipos de poro e fluidos desse reservatório. / [en] Since the 80 s, data-driven approaches have been used for fluids identification and reservoir characterization of siliciclastic and carbonate rocks mainly regarding seismic amplitudes analyses. However, techniques successfully applied for siliciclastic rocks, such as: AVO analysis, seismic inversions and DHI (Direct Hydrocarbon Indicators) ranking revealed not have achieved the same outstanding and reliable results for heterogeneous carbonate rocks. On the other hand, several articles demonstrate that reservoir characterization workflows with rock physics models embedded have been reaching a robust success in order to obtain petrophysical properties and elastic attributes of both rocks, from the seismic and well logs, in a model-driven approach focused on the reservoirs microstructural information. In this way, taking into account the importance of applying rock physics models in the scope of reservoir characterization, two rock physics models - Xu and Payne and T-Matrix - were applied, compared and their parameters were stochastically evaluated and optimized in a Bayesian framework. Through this approach, it was possible to estimate, in a reliable manner, the elastic attributes of a carbonate reservoir (coquinas) taking into consideration different kinds of uncertainties. Furthermore, after the calibration in the well location and validation of both rock physics models for other wells, sensitivity analyses were conducted in order to quantitatively understand how the coquinas elastic attributes behave regarding the variations in the reservoir mineralogical content, pore shapes and fluids.

Worldwide interest in gas production from shale formations has rapidly increased in recent years, mostly by the successful development of gas shales in North America. The Haynesville Shale is a productive gas shale resource play located in Texas and Louisiana. It produces primarily through enhanced exposure to the reservoir and improved permeability resulting from horizontal drilling and hydraulic fracturing. Accordingly, it is important to estimate the reservoir properties that influence the elastic and geomechanical properties from seismic data. This thesis estimates pore shapes, which affect the transport, elastic, and geomechancial properties, from wellbore seismic velocity in the Haynesville Shale. The approach for this work is to compare computed velocities from an appropriate rock physics model to measured velocities from well log data. In particular, the self-consistent approximation was used to calculate the model-based velocities. The Backus average was used to upscale the high-frequency well log data to the low-frequency seismic scale. Comparisons of calculated velocities from the self-consistent model to upscaled Backus-averaged velocities (at 20 Hz and 50 Hz) with a convergence of 0.5% made it possible to estimate pore aspect ratios as a function of depth. The first of two primary foci of this approach was to estimate pore shapes when a single fluid was emplaced in all the pores. This allowed for understanding pore shapes while minimizing the effects of pore fluids. Secondly, the effects of pore fluid properties were studied by comparing velocities for both patchy and uniform fluid saturation. These correspond to heterogeneous and homogeneous fluid mixing, respectively. Implementation of these fluid mixtures was to model them directly within the self-consistent approximation and by modeling dry-rock velocities, followed by standard Gassmann fluid substitution. P-wave velocities calculated by the self-consistent model for patchy saturation cases had larger values than those from Gassmann fluid substitution, but S-wave velocities were very similar. Pore aspect ratios for variable fluid properties were also calculated by both the self-consistent model and Gassmann fluid substitution. Pore aspect ratios determined for the patchy saturation cases were the smallest, and those for the uniform saturation cases were the largest. Pore aspect ratios calculated by Gassmann fluid substitution were larger because the velocity is inversely related to the aspect ratio in this particular modeling procedure. Estimates of pore aspect ratios for uniform saturation were 0.051 to 0.319 with the average of 0.171 from the velocity modeling using the self-consistent model. For patchy saturation, the aspect ratios were 0.035 to 0.296 with a mean of 0.145. These estimated pore aspect ratios from the patchy saturation case within the self-consistent model are considered the most reasonable set of values I determined. This is because the most likely in-situ fluid distribution is heterogeneous due to the extremely low permeability of the Haynesville Shale. Estimated pore aspect ratios using this modeling help us to understand elastic properties of the Haynesville Shale. In addition, this may help to find zones that correspond to optimal locations for fracturing the shale while considering brittleness and in-situ stress of the formation. / text

L'interférométrie sismique, tout comme l'interférométrie optique, étudie les phénomènes d'interférence entre des couples de signaux sismiques afin de mettre en évidence des différences entre ces signaux (par exemple Curtis et al. 2006). Les traitements utilisés consistent le plus souvent à corréler les enregistrements entre différents capteurs pour remonter aux fonctions de Green, (ou réponse impulsionnelle) entre ces récepteurs (par exemple, Derode et al. (2003), Wapenaar et al. (2004), Larose et al. (2006), Sanchez-Sesma et Campillo (2006)). Ce principe a déjà été appliqué avec succès dans les domaines de la sismologie (Campillo et Paul (2003)), des ultrasons (Weaver et Lobkis (2001)), de l'exploration sismique (Schuster (2001) et Wapenaar et al. (2004)), et même de l'hélio-sismologie (Duvall et al. (1993)) Dans tous ces cas, l'analyse des corrélations a conduit à une description détaillée des milieux de propagation, en l'occurrence l'intérieur de la Terre dans le cas de la sismologie. La sismique passive, par opposition à la sismique active utilisant les sources cohérentes artificielles (explosifs, canons à air, vibrateurs...), exploite les sources cohérentes naturelles (séismes...). Depuis peu la sismique passive exploite également des champs d'ondes aléatoires engendrées à des temps inconnus par une multitude de sources inconnues dans le sous-sol qui sont enregistrées à différentes positions de récepteurs. L'analyse par interférométrie sismique des enregistrements à deux capteurs permet de remonter aux fonctions de Green, ou réponse impulsionnelle, entre ces deux récepteurs (Derode et al. (2003)). Cette thèse, à vocations à la fois expérimentale et applicative, a deux buts : - vérifier au laboratoire, sur des expériences ultrasonores et avec différents types de roches, l'efficacité du monitoring de substitution de fluides par l'analyse interférométrique ultrasonore - appliquer sur le terrain les méthodes d'interférométrie sismique passive à des expériences ponctuelles de surveillance sismique passive d'exploitation du sous-sol. Ce manuscrit présente, après une synthèse bibliographique, la mise au point d'une nouvelle méthode de mesure des constantes élastiques d'un échantillon de roche (sèche ou saturée de fluide) basée sur les principes de l'interférométrie ultrasonore et de la spectroscopie par résonance ultrasonore. La méthode a été testée et validée (reproductibilité, fidélité, fiabilité…) sur un matériau standard de propriétés connues (aluminium). On expose que les effets de substitution fluide sont tout-à-fait mesurables avec la méthode sur divers échantillons de roches sèches puis saturées (en eau ou en éthylène glycol) et les résultats sont en accord raisonnable avec la théorie poroélastique de Biot-Gassmann. En outre, un certain nombre de faiblesses de la méthode ont été mises en évidence, à savoir la limitation à des roches assez homogènes et peu atténuantes. La dernière partie de ce manuscrit met en évidence des variations de vitesse des ondes dans un champ d'hydrocarbures (informations relatives à ce champ confidentielles) concomitantes avec le début de l'injection de vapeur dans celui-ci (pour récupération assistée de l'huile). / Seismic interferometry, like optic interferometry, studies the interferences phenomena between couples of seismic signals, with the aim of pointing at differences between those signals (Curtis et al. 2006 for instance). The data processing consists, generally, in correlating the recordings between different stations to retrieve the Green's function between these sensors (Derode et al. (2003), Wapenaar et al. (2004), Larose et al. (2006), Sanchez-Sesma et Campillo (2006)). This principle has already been successfully applied in the field of seismology (Campillo and Paul (2003)), ultrasound (Weaver et Lobkis (2001)), seismic exploration and even helioseismology (Duvall et al. (1993)). In all these cases, the analysis of the correlations leads to a detailed description of the propagation medium. Contrary to active seismic which uses artificial coherent sources (explosives, air guns…), passive seismic exploits natural coherent sources (seisms…). Since a few years, passive seismics also exploits random wave fields generated at unknown times by many unknown sources in the ground, and recorded at different stations positions. The analysis with cross-correlation of pairs of recordings, from pairs of sensors, leads to the Green's function between the two sensors (Derode et al. (2003)). This thesis has two objectives: -to check, at the lab scale, the effectiveness of monitoring of fluids substitution with noise correlation (ultrasound scale) in rocks -to apply noise correlation methods to passive seismic monitoring of a hydrocarbons field. This thesis presents, after a state of the art, the set-up of a new method to measure elastic constants of a rock sample (dry or fluid-saturated), based on ultrasound interferometry principle and resonant ultrasound spectroscopy. The method has been tested and validated (reproducibility, accuracy, precision…) on a standard material (aluminium). We show that the effects of a fluids substitution are measurable on various rock samples (dry or saturated, with water or with ethylene glycol) with this method. Plus the results are in agreement with Biot-Gassmann's theory. Besides, several weaknesses of the method were pointed, that is to say the method does not work on heterogeneous or attenuating medium. The last part of this thesis exposes speed of waves variations in a hydrocarbons field, when steam is injected simultaneously inside the reservoir (enhanced oil recovery operation).

Submitted by Jesiel Ferreira Gomes (jesielgomes@ufcg.edu.br) on 2018-04-17T23:15:52Z No. of bitstreams: 1 FELIPE DE OLIVEIRA APOLINÁRIO – DISSERTAÇÃO (PPGEPM) 2016.pdf: 5151929 bytes, checksum: f4706d54cb97c9b01a64299ddb28cd7d (MD5) / Made available in DSpace on 2018-04-17T23:15:52Z (GMT). No. of bitstreams: 1 FELIPE DE OLIVEIRA APOLINÁRIO – DISSERTAÇÃO (PPGEPM) 2016.pdf: 5151929 bytes, checksum: f4706d54cb97c9b01a64299ddb28cd7d (MD5) Previous issue date: 2016-09-06 / Capes / O presente trabalho teve como objetivo analisar a influência da saturação fluida nas propriedades elásticas de rochas carbonáticas, bem como a eficácia dos modelos de substituição de fluidos e de simulação computacional. Foram estudadas 9 amostras de rochas carbonáticas, sendo dois calcários laminados e sete tufas. As medições de velocidade foram realizadas em amostras secas, saturadas com água ou com óleo, sob diferentes pressões efetivas. A simulação de propagação de ondas foi feita no COMSOL Multiphysics 5.1, utilizando o Avizo Fire 8.1 para a criação das amostras digitais. Análise por difração de raios X (DRX) foi realizada para determinar a composição das amostras de rocha. Os resultados obtidos nos ensaios laboratoriais e nas simulações computacionais foram comparados com as estimativas dos modelos de substituição de fluidosGassmann, Biot e Brown & Korringa. Foi observado que a saturação das amostras com agua ou óleo geraram aumentos nas velocidades de propagação de onda P, porém sem apresentar um comportamento padrão. Para o caso das ondas S, a saturação por óleo predominantemente gerou aumentos nas velocidades, com exceção para os casos em que as amostras possuíam porosidade secundária do tipo vugular, devido à pouca influência que o óleo oferece para o módulo de cisalhamento nestes casos. A saturação por água resultou em diminuições das velocidades de propagação de onda S devido ao aumento da densidade total. Também foi constatado que o modelo de Gassmann foi o mais efetivo na estimativa de velocidades de onda P e S, enquanto que o de Biot mostrou-se eficaz apenas para a estimativa de velocidades de ondas S, sendo ineficiente para a estimativa de velocidades de ondas P, com erros de até 20%. A simulação computacional gerou resultados superdimensionados, porém que evidenciam que um aperfeiçoamento da metodologia, tal como o aumento do número de pontos de leitura pode gerar resultados mais próximos dos obtidos laboratorialmente e de maior confiabilidade. / This research aimed to analyze the influence of the saturating fluid in carbonate rocks, as well as verify the effectiveness of the fluid substitution models and computational simulations of wave propagation. To do so, nine carbonate rock samples were analyzed, which two of them were laminated limestones and seven were carbonate tufas. The measurement of velocities were made in dry, water saturated and oil saturated samples, under different effective pressures. The wave propagation simulations were made in COMSOL Multiphysics 5.1 using Avizo Fire 8.1 to generate the digital rock samples. The results obtained in lab procedures and in computer simulations were compared with the estimated velocities of the fluid substitution models of Gassmann, Biot and Brown & Korringa. It was observed that the saturation of the samples with water or oil resulted in an increasing of P-wave velocities, however without a pattern. The saturation with oil resulted in most cases in an increasing of S-wave velocities, the exceptions occurred in samples which had vugular porosity, due to the small influence of the oil in the shear modulus in this cases. T he saturation with water resulted in a decreasing of S-wave velocities due to the increment of the bulk density. It was found that the Gassmann’s model was more effective than the other two models in estimating P-wave and S-wave velocities. Biot’s model generated unsatisfying results to P-wave velocities, with errors up to 20%. However, this model had a good accuracy in estimating S-wave velocities. The computer simulations produced mainly overestimated results, though it was shown that an optimization of methodology, such as and addition in the number of the measure points, could improve the quality of the data, providing more representative results.

No description available.

La thèse a pour but d’étudier les effets de changements de mouillabilité de roches dans des conditions d’injections d’eau douce en tant que méthode de récupération d’hydrocarbures. Afin d’identifier le ou les mécanismes à l’origine du gain additionnel de récupération nous avons utilisé un microtomographe RX. Nous avons ainsi imagé les états de saturations finales d’un milieu poreux rempli de saumures et d’huiles. Une fois le drainage primaire réalisé nous avons effectué deux phases d’imbibitions : avec une saumure (récupération secondaire) puis une imbibition d’eau douce (récupération tertiaire). L’analyse de la mouillabilité à l’échelle du pore a permis de mettre en évidence l’effet de la température et de la salinité sur la mouillabilité. Nous avons montré que les changements de mouillages des roches n’étaient pas occasionnées par la seule expansion de la couche électrique en revanche des changements de mouillabilité ont été montrés. Ces changements s’expliquant par des transitions de mouillages de second ordre observées non seulement pour des gouttes d’huiles sur de l’eau mais également sur un substrat en verre. Au final, la mouillabilité en milieux poreux doit être mise en évidence à une échelle sous-Micrométrique ce qui est relativement nouveau dans le domaine pétrolier. / The aim of the thesis is to study rock wettability change effects caused by Low Salinity brine injection as tertiary recovery method. To identify the underlying mechanism or mechanisms of additional oil recovery X-Ray imaging technology was applied. We have also imaged the end-Point saturations of filled by brine and water core samples. Once the primary drainage is realized we carried out two phases imbibitions: with high salinity brine (waterflooding) and with low salinity brine (tertiary recovery mode). The wettability analysis at pore scale permitted to put in evidence the thermal and saline effects playing a decisive role in rock wettability. We have showed wettability changes are not caused by only electrical double layer expansion, however wettability changes was shown. These changes are explained by wettability transition of second order and observed not only for oil droplet on brine, but also for oil deposited on glass substrate. Finally, the pore space wettability needs to be evidenced at sub-Micrometric scale that is new for the petroleum domain.