Physical Properties of Deep Drill Cores : Implications for Meteorite Impact Effects and Crustal Structures

Physical properties provide valuable information about the nature and behavior of rocks and minerals. The changes in rock physical properties generate petrophysical contrasts between various lithologies, for example, between shocked and unshocked rocks in meteorite impact structures or between various lithologies in the crust. These contrasts may cause distinct geophysical anomalies, which are often diagnostic to their primary cause (impact, tectonism, etc). This information is vital to understand the fundamental Earth processes, such as impact cratering and associated crustal deformations. However, most of the present day knowledge of changes in rock physical properties is limited due to a lack of petrophysical data of subsurface samples, especially for meteorite impact structures, since they are often buried under post-impact lithologies or eroded. In order to explore the uppermost crust, deep drillings are required. This dissertation is based on the deep drill core data from three impact structures: (i) the Bosumtwi impact structure (diameter 10.5 km, 1.07 Ma age; Ghana), (ii) the Chesapeake Bay impact structure (85 km, 35 Ma; Virginia, U.S.A.), and (iii) the Chicxulub impact structure (180 km, 65 Ma; Mexico). These drill cores have yielded all basic lithologies associated with impact craters such as post-impact lithologies, impact rocks including suevites and breccias, as well as fractured and unfractured target rocks. The fourth study case of this dissertation deals with the data of the Paleoproterozoic Outokumpu area (Finland), as a non-impact crustal case, where a deep drilling through an economically important ophiolite complex was carried out. The focus in all four cases was to combine results of basic petrophysical studies of relevant rocks of these crustal structures in order to identify and characterize various lithologies by their physical properties and, in this way, to provide new input data for geophysical modellings. Furthermore, the rock magnetic and paleomagnetic properties of three impact structures, combined with basic petrophysics, were used to acquire insight into the impact generated changes in rocks and their magnetic minerals, in order to better understand the influence of impact. The obtained petrophysical data outline the various lithologies and divide rocks into four domains. Based on target lithology the physical properties of the unshocked target rocks are controlled by mineral composition or fabric, particularly porosity in sedimentary rocks, while sediments result from diverse sedimentation and diagenesis processes. The impact rocks, such as breccias and suevites, strongly reflect the impact formation mechanism and are distinguishable from the other lithologies by their density, porosity and magnetic properties. The numerous shock features resulting from melting, brecciation and fracturing of the target rocks, can be seen in the changes of physical properties. These features include an increase in porosity and subsequent decrease in density in impact derived units, either an increase or a decrease in magnetic properties (depending on a specific case), as well as large heterogeneity in physical properties. In few cases a slight gradual downward decrease in porosity, as a shock-induced fracturing, was observed. Coupled with rock magnetic studies, the impact generated changes in magnetic fraction the shock-induced magnetic grain size reduction, hydrothermal- or melting-related magnetic mineral alteration, shock demagnetization and shock- or temperature-related remagnetization can be seen. The Outokumpu drill core shows varying velocities throughout the drill core depending on the microcracking and sample conditions. This is similar to observations by Kern et al., (2009), who also reported the velocity dependence on anisotropy. The physical properties are also used to explain the distinct crustal reflectors as observed in seismic reflection studies in the Outokumpu area. According to the seismic velocity data, the interfaces between the diopside-tremolite skarn layer and either serpentinite, mica schist or black schist are causing the strong seismic reflectivities.

Caracterização de reservatórios não convencionais Shale Gas na Formação Barreirinha Bacia do Amazonas

Os folhelhos pretos devonianos da Formação Barreirinha caracterizamse pela alta radioatividade na porção basal, grande extensão areal, espessura e profundidade de soterramento variável que vão de exposição na superfície até mais de 3000 m. Eles são as principais rochas geradoras do sistema petrolífero convencional da Bacia do Amazonas, e recentemente foram consideradas como promissores plays de gás não convencional. Folhelhos são geralmente caracterizados por uma matriz fechada, que faz com que sejam relativamente impermeáveis em relação ao fluxo de gás, a menos que ocorram fraturas, e dependendo das suas características geológicas e geoquímicas podem funcionar com um Sistema Petrolífero autossuficiente, atuando tanto como rocha fonte, quanto como reservatório de gás (reservatório Shale Gas). Assim, o gás natural termogênico ou biogênico gerado pode ser armazenado em folhelhos ricos em matéria orgânica na forma livre, adsorvida, ou em estado dissolvido. Em contraste com os sistemas petrolíferos convencionais, reservatórios Shale Gas, possuem mecanismos de aprisionamento e armazenamento únicos, sendo necessária a utilização de técnicas de avaliações específicas. No entanto, folhelhos prolíficos geralmente podem ser reconhecidos a partir de alguns parâmetros básicos: arquitetura geológica e sedimentar, propriedades geoquímicas e petrofísicas e composição mineralógica. Tendo em vista a carência de pesquisas de caráter descritivo, com cunho exploratório dos folhelhos geradores da Formação Barreirinha, esta dissertação tem como objetivo introduzir uma metodologia de identificação de intervalos de folhelho gerador com potencial para reservatório Shale Gas. Começando com uma investigação regional sobre o contexto geológico e sedimentar, seguido de uma avaliação abrangente enfocando as características geoquímicas, petrofísicas e litofácies dos folhelhos a partir da integração de parâmetros obtidos de perfis geofísicos de poço, análises geoquímicas e aplicação dos conceitos de Estratigrafia de Sequencia.

鄂尔多斯盆地上古生界低渗透砂岩有效储层识别与预测

Sulige Gasfield, with a basically proven reserve as high as one trillion cubic meters, is one giant gas field discovered in China. The major gas -bearing layers are Upper Paleozoic strata with fluvial-lacustrine sedimentary facies. Generally, gas reservoirs in this field are characteristic by "five low" properties, namely low porosity, low permeability, low formation pressure, low productivity and low gas abundance. Reservoirs in this field also feature in a large distribution area, thin single sandbody thickness, poor reservoir physical properties, thin effective reservoir thickness, sharp horizontal and/or vertical changes in reservoir properties as well as poor connectivity between different reservoirs. Although outstanding achievements have been acquired in this field, there are still several problems in the evaluation and development of the reservoirs, such as: the relation between seismic attributes and reservoir property parameters is not exclusive, which yields more than one solution in using seismic attributes to predict reservoir parameters; the wave impedance distribution ranges of sandstone and mudstone are overlapped, means it is impossible to distinguish them through the application of post-stack impedance inversion; studies on seismic petrophysics, reservoir geophysical properties, wave reflection models and AVO features have a poor foundation, makes it difficult to recognize the specific differences between tight sandstone and gas-bearing sandstone and their distribution laws. These are the main reasons causing the low well drilling success rate and poor economic returns, which usually result in ineffective development and utilization of the field. Therefore, it is of great importance to perform studies on identification and prediction of effective reservoirs in low permeable sandstone strata. Taking the 2D and 3D multiwave-multicomponent seismic exploration block in Su6-Su5 area of Sulige field as a study area and He 8 member as target bed, analysis of the target bed sedimentary characteristics and logging data properties are performed, while criteria to identify effective reservoirs are determined. Then, techniques and technologies such as pre-stack seismic information (AVO, elastic impedance, wave-let absorption attenuation) and Gamma inversion, reservoir litological and geophysical properties prediction are used to increase the precision in identifying and predicting effective reservoirs; while P-wave and S-wave impedance, ratio of P/S wave velocities, rock elastic parameters and elastic impedance are used to perform sandstone gas-bearing property identification and gas reservoir thickness prediction. Innovative achievements are summarized as follows: 1. The study of this thesis is the first time that multiwave-multicomponent seismic data are used to identify and predict non-marine classic reservoirs in China. Through the application of multiwave-multicomponents seismic data and integration of both pre-stack and post-stack seismic data, a set of workflows and methods to perform high-precision prediction of effective reservoirs in low permeable sandstone is established systematically. 2. Four key techniques to perform effective reservoir prediction including AVO analysis, pre-stack elastic wave impedance inversion, elastic parameters inversion, and absorption attenuation analysis are developed, utilizing pre-stack seismic data to the utmost and increasing the correct rate for effective reservoir prediction to 83% from the former 67% with routine methods. 3. This thesis summarizes techniques and technologies used in the identification reservoir gas-bearing properties using multiwave-multicomponent seismic data. And for the first time, quantitative analysis on reservoir fluids such as oil, gas, and/or water are carried out, and characteristic lithology prediction techniques through the integration of pre-stack and post-stack seismic prediction techniques, common seismic inversion and rock elastic parameters inversion, as well as P-wave inversion and converted wave inversion is put forward, further increasing the correct rate of effective reservoir prediction in this area to 90%. 4. Ten seismic attribute parameters are selected in the 3D multi-wave area to perform a comprehensive evaluation on effective reservoirs using weighted-factor method. The results show that the first class effective reservoir covers an area of 10.08% of the study area, while the second and the third class reservoirs take 43.8% and 46% respectively, sharply increasing the success rate for appraisal and development wells.

孤东油田馆陶组河流相储层流动单元模型与剩余油分布规律研究

Since C.L. Hearn presented the concept of flow unit in 1984, its generation mechanisms and controlling factors have been studied in many aspects using different methods by researchers. There are some basic methods to do the research, and there are several concepts and classification standards about flow unit. Based on previous achievements and using methodologies from sedimentary geology, geophysics, seismic stratigraphy, and reservoir engineering, the author systemically studies the factors controlling flow unit, puts forward a series of methods for recognition, classification and evaluation of flow unit. The results obtained in this paper have important significance not only for understanding the flow unit, but also for revealing the distribution of remaining oil. As a case, this paper deals with the reservoir rocks in Guantao Group of Gudong Oilfield. Zhanhua Sag, Jiyang Depression in Bohaiwan Basin. Based on the study of stratigraphic, depositional and structural characteristics, the author establishes reservoir geological models, reveals the geological characteristics of oil-bearing reservoir of fluvial facies, points out the factors controlling flow unit and geological parameters for classification of flow unit. and summarizes methods and technologies for flow unit study when geological, well-logging and mathematical methods are used. It is the first attempt in literatures to evaluate reservoir by well-logging data constrained by geological conditions, then a well-logging evaluation model can be built. This kind of model is more precise than ever for calculating physical parameters in flow unit. In a well bore, there are six methods to recognize a flow unit. Among them, the activity function and intra-layer difference methods are the most effective. Along a section, the composition type of flow unit can be located according amplitude and impedance on seismic section. Slice method and other methods are used to distinguish flow unit. In order to reveal the distribution laws of flow unit in space, the author create a new method, named combination and composition of flow unit. Based on microscopic pore structure research, the classification methods of flow unit are developed. There are three types of flow unit in the reservoir of fluvial facies. They have their own lithology, petrophysics and pore structure character. Using judgement method, standard functions are built to determine the class of flow unit of fluvial facies. Combining reservoir engineering methods, the distribution laws of remaining oil in different types, or in different part of a flow unit are studied. It is evident that the remaining oil is controlled by the type of flow unit. The author reveals the relationship between flow unit and remaining oil distribution, builds the flowing models, predicts the variation of reservoir parameters in space, put forward different methods developing remaining oil in different flow unit. Especially, based on the results obtained in this paper, some suggestions for the adjustment of the developing flow units have been applied in Districts No.4 and No.7, and good results have been yielded. So, the results of this paper can guide oil field development. They are useful and significant for developing the remaining oil and enhancing the oil recovery efficiency.

东营凹陷北部陡坡带砂砾岩体性油气藏形成机理及分布规律研究

The north steep slope zone of Dongying Depression has great potential in oil resource and as the usage of 3-d seismic data in the last decade, the exploration of oil and gas has get into the stage of sandy glavel body lithological oil-gas pool exploration. In this thesis, writer take the north steep slope zone of Dongying Depression as target area and take Sha-III and Sha-[V Menber as purpose stratum, study on sequence stratigraphy, depositional system, reservoir description, emphasesing on analyzing of forming of lithological oil-gas pool, especially the dynamics principle of oil and gas preliminary movement from the source rock to the reservoir form lithological oil-gas pools. The aim of this work is to give some quantitatively explanation for the mechanism of lithological oil-gas pool forming, and set up the theory of pool form with characteristic terrestrial faulted basin. There are main conclusions and views as follow. 1. Applying with principle of sequence stratigrapgy, according to the depositional cycles of Dongying Depression, the sequence stratigraphical partition of Tertiary was finished, stressing on dismembering Sha-III and Sha-IV Menber as 5system tracts. 2. The structure of Dongying Depression especially of the north steep slope zone has accomplished, including the analyzing the structural cortroling to depositional condition of the north steep slope zone of Dongying Depression, discussed relationship between the structure of the north steep slope zone and the pool-forming. 3. The horizontal and vertical exchanges of ancient climates and ancient physiognomy of the all stratum units and studies on characteristic of depositional system distribution have been finished, found that there are five depositional systems in the north steep slope zone of Dongying Depression as fluvial, delta (tan-delta), sub-water fluvial fan lacustrine, gravitive flow, and seven formations of sandy glavel body, and forecasting of all kinds of sandy glavel body has been made. 4. Seismic stratigraphy and log stratigraphy have been made, described and forecasted all kinds of reservoir of objective stratum by means of physical geography method, setup a series means of sandy glavel body description suit to target area. 5. The pool-forming system has been studied, analyzing all the elements in petroleum sub-system of Sha-III and Sha-IV Menber of Dongying Depression with view of source controlling, estimated the petroleum system applying source rock potential index combining with distribution ofreservior. 6.Through studying types of pool, the controlling factors of pool-forming of sandy glavel body were discussed by deposition stages, formation types, structure ect. as a conclusion that the characteristics of pool forming in the north steep slope zone of Dongying Depression are, the controlling factor of the pools is mainly lithology, petrophysics of oil sands vary greatly, with a large heterogeneity, all kind of reservoir with different formation has different pool-forming conditions, and as a result, formed various pools of sandy glavel body along the steep slope with regular combination, distribution and constituted the multiple petroleum accumulative pattern. 7. It's the first time to cauculate and estimate the fluid pressure in source rock of Dongying Depression, set up the stratum fluid pressure in Dongying Depression, and firstly use equivalent charging pressure and reservoir forming index to quantitatively evaluate the pool-forming condition of lithological pool.8. Above all studies, follow up the scent of the exploration combined with practice a lot of explorative targets were found, and got geat economic and social benefit.

Determinação da distribuição de um nível de turbidito na formação Macaé no campo de Namorado através de uso de atributos sísmicos

This study was developed on the basis of petrophysics characterization, 3D seismic interpretation, and rock-wire-log integration. The results were integrated with amplitude attribute through cross plot (rock attribute x amplitude value). 3D seismic data also allowed inteipreting the top of Lagoa Feia, Macaé, and Carapebus formations, and Outeiro Member. The analyses of the maps of structures and attributes, integrated with wire-log data, suggest that Namorado turbidite level is candidate to be mapped with stratigraphic attributes like AverageTroughAmplitude, RMS velocity and Total AbsoluteAmplitude. Lagoa Feia Formation attribute maps (Dip-Azimuth and TotalEnergy) allows to interpret a low fault density at the carbonate coquina level in the Namorado field. This level is also considered one at the best seismic reflector and an important reservoir of the Campos Basin.

Fluxo da água subterrânea em sistema de encosta-rio, município de Paulínia (SP): caracterização hidrogeológica e simulação numérica

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

Determinação automática da porosidade e zoneamento de perfis através da rede neural artificial competitiva

Duas das mais importantes atividades da interpretação de perfis para avaliação de reservatórios de hidrocarbonetos são o zoneamento do perfil (log zonation) e o cálculo da porosidade efetiva das rochas atravessadas pelo poço. O zoneamento é a interpretação visual do perfil para identificação das camadas reservatório e, consequentemente, dos seus limites verticais, ou seja, é a separação formal do perfil em rochas reservatório e rochas selante. Todo procedimento de zoneamento é realizado de forma manual, valendo-se do conhecimento geológico-geofísico e da experiência do intérprete, na avaliação visual dos padrões (características da curva do perfil representativa de um evento geológico) correspondentes a cada tipo litológico específico. O cálculo da porosidade efetiva combina tanto uma atividade visual, na identificação dos pontos representativos de uma particular rocha reservatório no perfil, como a escolha adequada da equação petrofísica que relaciona as propriedades físicas mensuradas da rocha com sua porosidade. A partir do conhecimento da porosidade, será estabelecido o volume eventualmente ocupado por hidrocarboneto. Esta atividade, essencial para a qualificação de reservatórios, requer muito do conhecimento e da experiência do intérprete de perfil para a efetiva avaliação da porosidade efetiva, ou seja, a porosidade da rocha reservatório, isenta do efeito da argila sobre a medida das propriedades físicas da mesma. Uma forma eficiente de automatizar estes procedimentos e auxiliar o geofísico de poço nestas atividades, que particularmente demandam grande dispêndio de tempo, é apresentado nesta dissertação, na forma de um novo perfil, derivado dos perfis tradicionais de porosidade, que apresenta diretamente o zoneamento. Pode-se destacar neste novo perfil as profundidades do topo e da base das rochas reservatório e das rochas selante, escalonado na forma de porosidade efetiva, denominado perfil de porosidade efetiva zoneado. A obtenção do perfil de porosidade efetiva zoneado é baseado no projeto e execução de várias arquiteturas de rede neural artificial, do tipo direta, com treinamento não supervisionado e contendo uma camada de neurônios artificiais, do tipo competitivo. Estas arquiteturas são projetadas de modo a simular o comportamento do intérprete de perfil, quando da utilização do gráfico densidade-neutrônico, para as situações de aplicabilidade do modelo arenito-folhelho. A aplicabilidade e limitações desta metodologia são avaliadas diretamente sobre dados reais, oriundos da bacia do Lago Maracaibo (Venezuela).

Modelo tectono-estrutural dos campos de petróleo Namorado, Cherne e Albacora

The main goal in this research is a tectono-estructural characterization of the Cherne, Albacora and Namorado Fields, located at Campos Basin, in order to investigate the relationship between the geologic evolution and the rock´s physical properties of the reservoir, and how they affect the hydrocarbon accumulation in those fields. Well correlations show that the inferior turbidites have a regional lateral continuity. Basic petrophysics analysis, calculated here, shows that the three fields present porosity values that range from 15 to 20%, shale volume range from 26 to 30% and formation water saturation range from 23 to 45%, based on formation water resistivity dada from Albacora Field. Petrophysics maps feature a trend in Albacora Field that increase the porosity values to SE, and in Cherne and Namorado Field the trend increase towards N. Seismic horizons where interpreted between the first appearance of the Namorado Sandstone and the top of Quissamã Formation. This interval presents normal listric faulting, in Cherne and Namorado Field with NWSE and NE-SW direction, and sedimentation trend to NW-SE, in Albacora Field the faulting presents a NNE-SSW and N-S direction, with a sedimentation trend to NE-SW. Seismic attribute maps present amplitude anomalies close to the producing wells, and on Namorado Field, it indicates a potential hydrocarbon accumulation in the NE region. For each field is indicated laboratory tests for a better characterization of the petrophysical properties, since that they don’t form the same reservoir level, therefore, not influencing the water saturation calculation

Ab initio quantum mechanical investigation of structural and chemical-physical properties of selected minerals for minero-petrological, structural ceramic and biomaterial applications

The purpose of this thesis is the atomic-scale simulation of the crystal-chemical and physical (phonon, energetic) properties of some strategically important minerals for structural ceramics, biomedical and petrological applications. These properties affect the thermodynamic stability and rule the mineral-environment interface phenomena, with important economical, (bio)technological, petrological and environmental implications. The minerals of interest belong to the family of phyllosilicates (talc, pyrophyllite and muscovite) and apatite (OHAp), chosen for their importance in industrial and biomedical applications (structural ceramics) and petrophysics. In this thesis work we have applicated quantum mechanics methods, formulas and knowledge to the resolution of mineralogical problems ("Quantum Mineralogy”). The chosen theoretical approach is the Density Functional Theory (DFT), along with periodic boundary conditions to limit the portion of the mineral in analysis to the crystallographic cell and the hybrid functional B3LYP. The crystalline orbitals were simulated by linear combination of Gaussian functions (GTO). The dispersive forces, which are important for the structural determination of phyllosilicates and not properly con-sidered in pure DFT method, have been included by means of a semi-empirical correction. The phonon and the mechanical properties were also calculated. The equation of state, both in athermal conditions and in a wide temperature range, has been obtained by means of variations in the volume of the cell and quasi-harmonic approximation. Some thermo-chemical properties of the minerals (isochoric and isobaric thermal capacity) were calculated, because of their considerable applicative importance. For the first time three-dimensional charts related to these properties at different pressures and temperatures were provided. The hydroxylapatite has been studied from the standpoint of structural and phonon properties for its biotechnological role. In fact, biological apatite represents the inorganic phase of vertebrate hard tissues. Numerous carbonated (hydroxyl)apatite structures were modelled by QM to cover the broadest spectrum of possible biological structural variations to fulfil bioceramics applications.

Reservoir characterization of the Cardium Formation in the Ferrier Oilfield, west-central Alberta, Canada.

This thesis has the main aim of defining the lithostratigraphy, depositional architecture, post-depositional modifications and reservoir characteristics of the Cardium Formation in the Ferrier Oilfield, and how these characteristics can have great impact over production rates, GOR and produced fluid discrimination. In the Ferrier area, the Cardium Formation is composed by a NE prograding clastic sequence made up of offshore to shoreface deposits sealed by marine shales. The main reservoir is composed by sandstones and conglomerates interpreted to have deposited in a shoreface depositional environment. Lithofacies and net reservoir thickness mapping led to more detailed understanding of the 3D reservoir architecture, and cross-sections shed light on the Cardium depositional architecture and post-deposition sediment erosion in the Ferrier area. Detailed core logging, thin section, SEM and CL analyses were used to study the mineralogy, texture and pore characterization of the Cardium reservoir, and three main compartments have been identified based on production data and reservoir characteristics. Finally, two situations showing odd production behaviour of the Cardium were resolved. This shed light on the effect of structural features and reservoir quality and thickness over hydrocarbon migration pathways. The Ferrier example offers a unique case of fluid discrimination in clastic reservoirs due both to depositional and post-depositional factors, and could be used as analogue for similar situations in the Western Canadian Sedimentary Basin.

Petrophysical investigation on sediment core CRP-3 from the Ross Sea, Antarctica

A suite of petrophysical properties - velocity, resistivity, bulk density, porosity, and matrix density - was measured on 88 core plugs from the CRP-3 drillhole. Core-plug bulk densities were used to recalibrate both whole-core and downhole bulk density logs. Core-plug measurements of matrix density permit conversion of the whole-core and downhole bulk density logs to porosity. Both velocity and formation factor (a normalized measure of resistivity) are strongly correlated with porosity. The velocity/porosity pattern is similar to that for the lower part of CRP-2A and is consistent with the empirical relationship for sandstones. Core-plug and whole-core measurements of P-wave velocity at atmospheric pressure exhibit excellent agreement. Measurements of velocity as a function of pressure indicate a significantly higher velocity sensitivity to pressure than has been observed at CRP-1 and CRP-2A; rebound or presence of microcracks at CRP-3 may be responsible. The percentage difference between velocities at in situ pressures and atmospheric pressures increases downhole from 0% at the seafloor to 9% at the bottom. This pattern can be used to correct whole-core velocity data, measured at atmospheric pressure, to in situ velocities for depth-to-time conversion and associated comparison to the seismic profile across the drillsite