No description available.

I address the problem of finding a background model for the estimation of resistivities in the earth from controlled-source electromagnetic (CSEM) data by using seismic data and well logs as constraints. Estimation of resistivities is normally done by trial-and-error, in a process called “inversion”, by finding a model of the earth whose responses match the data to within an acceptable error; what comes out of the inversion is what is put into the model by the geophysicist: it does not come out of the data directly. The premise underlying this thesis is that an earth model can be found that satisfies not only the CSEM data but also the seismic data and any well logs. I present a methodology to determine background resistivities from seismic velocities using rock physics, structural constraints, and depth trends. The physical parameters of the seismic wave equation are different from those in the electromagnetic diffusion equation, so there is no direct link between the governing equations. I therefore use a Bayesian framework to incorporate not only the errors in the data and our limited knowledge of the rock parameters, but also the uncertainty of our chosen and calibrated velocity-to-resistivity transform. To test the methodology I use a well log from the North Sea Harding South oil and gas field to calibrate the transform, and apply it to seismic velocities of the nearby Harding Central oil and gas field. I also use short-offset CSEM inversions to estimate the electric anisotropy and to improve the shallow part of the resistivity model, where there is no well control. Three-dimensional modelling of this resistivity model predicts the acquired CSEM data within the estimated uncertainty. This methodology makes it possible to estimate background resistivities from seismic velocities, well logs, and other available geophysical and geological data. Subsequent CSEM surveys can then focus on finding resistive anomalies relative to this background model; these are, potentially, hydrocarbon-bearing formations.

O presente trabalho de mestrado está inserido nos estudos de refração profunda do Projeto Temático 'Estudos Geofísicos e Modelo Tectônico dos Setores Central e Sudeste da Província Tocantins, Brasil Central'. Nesses estudos foram levantadas três linhas de refração de aproximadamente 300 km de extensão, duas no setor Central da Província Tocantins e uma no setor Sudeste, que é o objeto de estudo deste trabalho. Foram utilizados 111 sismógrafos digitais SGR pertencentes ao programa PASSCAL, instrumentos auxiliares do USGS, e 13 sismógrafos digitais e instrumentos auxiliares do IAG/USP. A linha sísmica teve aproximadamente 300 km de extensão com pontos de registro separados a cada 2,5 km, distribuídos ao longo de estradas principais e secundárias. A cada 50 km, aproximadamente, foi realizada uma explosão, nas explosões dos extremos da linha foram utilizados 1000 kg de explosivo e para a explosão central uma carga de 500 kg. Para a determinação das coordenadas geográficas dos pontos de tiro e de registro, foi utilizado o método diferencial com medidas de GPS. O principal objetivo deste trabalho foi obter como produto final um modelo de velocidades sísmicas contendo as características físicas das principais descontinuidades na crosta terrestre e no manto superior. Para análise e processamento dos dados foram utilizados os pacotes SAC, SU, SEIS. Para a modelagem foram utilizados a teoria do raio e a elaboração de sismogramas sintéticos, do pacote SEIS. Para a elaboração do modelo final foram utilizados os dados das explosões dos pontos extremos e central, tendo em vista que devido a problemas técnicos não foram registrados os sinais das outras 4 explosões. Além disso, as explosões registradas não apresentaram sinais claros em toda a extensão da linha. Devido a tudo isso e considerando as unidades geológicas presentes na região de estudo são sugeridos três modelos de velocidades sísmicas. O primeiro modelo refere-se ao tiro direto (EX31) localizado no extremo sudoeste da linha, sobre a Bacia do Paraná. Para este modelo obteve-se para superfície (0 km) a velocidade inicial de 2 km/s (coberturas); para a profundidade de 0,086 km a velocidade inicial é de 5,15 km/s (basalto); para a profundidade de 0,350 km obteve-se a velocidade inicial de 4,6 km/s (arenito - camada de baixa velocidade); para profundidade de 0,650 km a velocidade inicial é de 5,75 km/s e para profundidade de 4 km obteve-se a velocidade inicial de 6,07 km/s. O segundo modelo refere-se ao tiro reverso (EX34) localizado no centro da linha sobre granitóides do Grupo Araxá. Para este modelo obteve -se para superfície (0 km) a velocidade inicial de 2 km/s; para a profundidade de 0,06 km a velocidade inicial é de 5,69 km/s e para a profundidade de 0,860 km obteve-se a velocidade inicial de 6,25 km/s. Finalmente, o terceiro modelo refere-se ao tiro direto para toda a extensão da linha (300 km). Este modelo foi definido a partir de fases secundárias lidas nos registros e modelos anteriores propostos na literatura. Da superfície até os 4 km iniciais de profundidade este modelo é igual ao primeiro, para uma profundidade de 20 km obteve-se a velocidade inicial de 6,70 km/s e para uma profundidade de 40 km a velocidade é de 8,00 km/s (descontinuidade de MOHO). / This work to fulfil the degree of Master of Sciences is inserted among the deep seismic refraction studies of the Thematic Project 'Geophysical Studies and Tectonic Model of the Tocantins Province Central and Southeast Sectors, Central Brazil'. Three refraction lines, of around 300 km long each, were deployed, two of them in the Central sector and the other in the SE sector, that is subject of the present work. The equipment used in this experiment was composed by 111 SGR digital seismographs belonging to the PASSCAL Program. Complemented with auxiliary instruments from USGS and 13 seismographs belonging to IAG/USP. The space among the recording points was 2.5 km, which were located along main and secondary roads. Every 50 km was fired an explosion with 1000 kg of emulsion in each extreme and 500 kg in the central point. The geographical co-ordinates were determined by using the GPS differential method. The main objective of this work is to obtain as a final product a seismic velocity model with the physical characteristics of the main discontinuities in the crust and upper mantle. The packages SAC, SU and SEIS were used to perform the data analysis and processing. To carry on the modelling were used the ray theory and the synthetic seismograms construction, belonging to the SEIS package Data from the extreme and middle points of the seismic line were used to elaborate the final model, considering that due to technical problems signals from the other four explosions were not recorded. Apart from that, the recorded explosions did not present clear signals all along the extension of the line. Due to these facts, and considering also the geological units present in the studied region, are suggested three seismic velocity models. The first model is referred to the direct shot (EX31), which is localised in the Southwest extreme of the line on the Parana Basin province. In this model we obtained the P wave velocity (VP) of 2 km/sec at the surface, corresponding to the unconsolidated sediments and soil on the top of that basin. At a depth of 86 m we found VP of 5,15 km/sec and at a depth of 350 m the velocity VP of 4,6 km/sec, corresponding to the basalt and sand layers of the Parana Basin. Underlying them, at 650 m of depth we found the basement with VP of 5,75 km/sec and finally at a depth of 4 km there is a layer with VP of 6,07 km/sec, corresponding to a typical upper crust P wave velocity. The second model corresponds to the reverse shot (EX34) that is localised in the middle point of the line on the granitoides of the Araxa Group. For this model we obtained VP of 2 km/sec for the superficial layers, then at a depth of 60 m was obtained V P of 5,69 km/sec and for a depth of 860 m the value of V P is 6,25 km/sec. Finally, the third model belongs to the whole line section (300 km) from the direct shot (EX31). This model was obtained by using the arrivals of secondary phases and the results of models proposed in other works. From the surface down to 4 km of depth this model is similar to the first one. At 20 km of depth there is a layer with VP of 6,70 km/sec, corresponding to the lower crust, with Moho at a depth of 40 km with VP of 8,00 km/sec.

Black shales are sedimentary rocks with a high content of organic carbon, which leads to a dark grayish to black color. Due to their potential to contain oil or gas, black shales are of great interest for the support of the worldwide energy supply. An integrated seismic investigation of the Lower Palaeozoic black shales was carried out at the Danish island Bornholm to locate the shallow-lying Alum Shale layer and its surrounding formations and to characterize its potential as a source rock. Therefore, two seismic experiments at a total of three crossing profiles were carried out in October 2010 and in June 2012 in the southern part of the island. Two different active measurements were conducted with either a weight drop source or a minivibrator. Additionally, the ambient noise field was recorded at the study location over a time interval of about one day, and also a laboratory analysis of borehole samples was carried out. The seismic profiles were positioned as close as possible to two scientific boreholes which were used for comparative purposes. The seismic field data was analyzed with traveltime tomography, surface wave inversion and seismic interferometry to obtain the P-wave and S-wave velocity models of the subsurface. The P-wave velocity models which were determined for all three profiles clearly locate the Alum Shale layer between the Komstad Limestone layer on top and the Læså Sandstone Formation at the base of the models. The black shale layer has P-wave velocities around 3 km/s which are lower compared to the adjacent formations. Due to a very good agreement of the sonic log and the vertical velocity profiles of the two seismic lines, which are directly crossing the borehole where the sonic log was conducted, the reliability of the traveltime tomography is proven. A correlation of the seismic velocities with the content of organic carbon is an important task for the characterization of the reservoir properties of a black shale formation. It is not possible without calibration but in combination with a full 2D tomographic image of the subsurface it gives the subsurface distribution of the organic material. The S-wave model obtained with surface wave inversion of the vibroseis data of one of the profiles images the Alum Shale layer also very well with S-wave velocities around 2 km/s. Although individual 1D velocity models for each of the source positions were determined, the subsurface S-wave velocity distribution is very uniform with a good match between the single models. A really new approach described here is the application of seismic interferometry to a really small study area and a quite short time interval. Also new is the selective procedure of only using time windows with the best crosscorrelation signals to achieve the final interferograms. Due to the small scale of the interferometry even P-wave signals can be observed in the final crosscorrelations. In the laboratory measurements the seismic body waves were recorded for different pressure and temperature stages. Therefore, samples of different depths of the Alum Shale were available from one of the scientific boreholes at the study location. The measured velocities have a high variance with changing pressure or temperature. Recordings with wave propagation both parallel and perpendicular to the bedding of the samples reveal a great amount of anisotropy for the P-wave velocity, whereas the S-wave velocity is almost independent of the wave direction. The calculated velocity ratio is also highly anisotropic with very low values for the perpendicular samples and very high values for the parallel ones. Interestingly, the laboratory velocities of the perpendicular samples are comparable to the velocities of the field experiments indicating that the field measurements are sensitive to wave propagation in vertical direction. The velocity ratio is also calculated with the P-wave and S-wave velocity models of the field experiments. Again, the Alum Shale can be clearly separated from the adjacent formations because it shows overall very low vP/vS ratios around 1.4. The very low velocity ratio indicates the content of gas in the black shale formation. With the combination of all the different methods described here, a comprehensive interpretation of the seismic response of the black shale layer can be made and the hydrocarbon source rock potential can be estimated. / Schwarzschiefer sind Sedimentgesteine, die einen hohen Gehalt an organischem Kohlenstoff aufweisen, was zu einer dunkelgrauen bis schwarzen Färbung führt. Da Schwarzschiefer das Potenzial besitzen, Öl oder Gas zu enthalten und somit zur weltweiten Energieversorgung beitragen könnten, sind sie von großem Interesse. Mit Hilfe der Kombination verschiedener seismischer Messverfahren wurden die Schwarzschiefer des Unteren Paläozoikums auf der dänischen Insel Bornholm untersucht um den oberflächennahen Alaunschiefer und dessen Umgebungsgestein dort zu lokalisieren und sein Potenzial als Muttergestein abzuschätzen. Dafür wurden im Oktober 2010 und im Juni 2012 im südlichen Teil der Insel zwei seismische Experimente auf insgesamt drei sich kreuzenden Profilen durchgeführt. Für zwei aktive seismische Messungen wurden ein Fallgewicht und ein Minivibrator als Quellen genutzt. Zusätzlich wurde im Messgebiet noch das Wellenfeld des umgebenden Rauschens über einen Zeitraum von etwa einem Tag aufgezeichnet. Außerdem wurden Labormessungen an Bohrkernen aus dem Alaunschiefer durchgeführt. Die seismischen Messprofile befanden sich so nah wie möglich an zwei wissenschaftlichen Bohrungen, die für Vergleichszwecke genutzt wurden. Um die P- und S-Wellengeschwindigkeitsmodelle des Untergrundes zu erhalten wurden die seismischen Felddaten mittels Laufzeittomographie, Oberflächenwelleninversion und seismischer Interferometrie ausgewertet. Die P-Wellenmodelle, die für alle drei seismischen Profile erstellt wurden, zeigen den Alaunschiefer zwischen dem Komstad Kalkstein, der den Alaunschiefer überdeckt, und der Læså Sandsteinformation, die die Basis der Modelle bildet. Für die Schwarzschieferschicht ergeben sich mit rund 3 km/s deutlich geringere P-Wellengeschwindigkeiten als für die umgebenden Gesteine. Zwei seismische Profile liegen direkt an einer der Bohrungen, für die verschiedene Bohrloch-Logs durchgeführt wurden. Der Vergleich des Sonic-Logs mit den vertikalen Geschwindigkeitsprofilen beider Modelle am Bohrpunkt zeigt eine sehr gute übereinstimmung aller Geschwindigkeiten. Dies ist ein Indiz für die Plausibilität der durchgeführten Laufzeittomographie. Um die Reservoireigenschaften der Schwarzschieferschicht einordnen zu können, wurde versucht, die seismischen Geschwindigkeiten mit dem Gehalt an organischem Material zu korrelieren. Ohne geeignete Kalibrierung ist diese Korrelation schwierig, kann aber mit Hilfe der Tomographieergebnisse ein zweidimensionales Abbild der Verteilung des organischen Materials im Untergrund liefern. Auch das S-Wellengeschwindigkeitsmodell, welches mit der Oberflächenwelleninversion der Vibroseisdaten erstellt wurde, bildet den Alaunschiefer gut ab. Hierbei zeigen sich S-Wellengeschwindigkeiten um 2 km/s. Obwohl jeweils nur 1D-Modelle für jede Quellposition bestimmt wurden, ergibt sich für die gesamte Untergrundstruktur des untersuchten Profils ein einheitliches Bild der Geschwindigkeiten. Einen sehr neuen Ansatz bildet die Anwendung der seismischen Interferometrie auf ein sehr kleines Untersuchungsgebiet und über einen sehr kurzen Zeitraum. Neu ist außerdem, dass für die Bestimmung der endgültigen Interferogramme nur Zeitfenster der Kreuzkorrelationen ausgewählt werden, in denen die Signalqualität hinreichend gut ist. In den berechneten Kreuzkorrelationen sind sogar P-Wellen enthalten, was auf die geringen Abstände der seismischen Rekorder zurück zu führen ist. Bei den Labormessungen wurden die Raumwellen für verschiedene Drücke und Temperaturen aufgezeichnet. Die Messungen der Geschwindigkeiten sowohl parallel als auch senkrecht zur Schichtung der Proben zeigen eine starke Anisotropie für die P-Welle. Dagegen scheint die S-Wellengeschwindigkeit fast unabhängig von der Ausbreitungsrichtung der Wellen zu sein. Auch das Verhältnis der Geschwindigkeiten weist starke Anisotropie auf. Für die Wellenausbreitung senkrecht zur Schichtung zeigen sich sehr niedrige Werte, die Werte für die Messungen parallel zur Schichtung sind dagegen deutlich erhöht. Ein interessanter Aspekt der aus den Labormessungen resultiert ist, dass die Geschwindigkeit der Messungen senkrecht zur Schichtung mit den Geschwindigkeitswerten der Feldmessungen übereinstimmen. Damit scheinen die Feldmessungen besonders die Ausbreitung der Wellen in vertikaler Richtung zu registrieren. Das Geschwindigkeitsverhältnis wurde auch mit den P- und S-Wellenmodellen der Feldexperimente berechnet. Auch hier hebt sich der Alaunschiefer mit deutlich verringerten Werten um 1.4 vom Umgebungsgestein ab. Solch geringe Werte für das Verhältnis der Geschwindigkeiten deutet auf den Gehalt von Gas im Schwarzschiefer. Mit der Kombination der verschiedenen Methoden ist es möglich, die seismische Antwort der Schwarzschieferschicht umfassend zu beschreiben und Schlussfolgerungen darüber zu ziehen, ob die hier untersuchte Schwarzschieferschicht das Potenzial hat als Kohlenwasserstofflagerstätte zu fungieren.

O presente trabalho de mestrado está inserido nos estudos de refração profunda do Projeto Temático 'Estudos Geofísicos e Modelo Tectônico dos Setores Central e Sudeste da Província Tocantins, Brasil Central'. Nesses estudos foram levantadas três linhas de refração de aproximadamente 300 km de extensão, duas no setor Central da Província Tocantins e uma no setor Sudeste, que é o objeto de estudo deste trabalho. Foram utilizados 111 sismógrafos digitais SGR pertencentes ao programa PASSCAL, instrumentos auxiliares do USGS, e 13 sismógrafos digitais e instrumentos auxiliares do IAG/USP. A linha sísmica teve aproximadamente 300 km de extensão com pontos de registro separados a cada 2,5 km, distribuídos ao longo de estradas principais e secundárias. A cada 50 km, aproximadamente, foi realizada uma explosão, nas explosões dos extremos da linha foram utilizados 1000 kg de explosivo e para a explosão central uma carga de 500 kg. Para a determinação das coordenadas geográficas dos pontos de tiro e de registro, foi utilizado o método diferencial com medidas de GPS. O principal objetivo deste trabalho foi obter como produto final um modelo de velocidades sísmicas contendo as características físicas das principais descontinuidades na crosta terrestre e no manto superior. Para análise e processamento dos dados foram utilizados os pacotes SAC, SU, SEIS. Para a modelagem foram utilizados a teoria do raio e a elaboração de sismogramas sintéticos, do pacote SEIS. Para a elaboração do modelo final foram utilizados os dados das explosões dos pontos extremos e central, tendo em vista que devido a problemas técnicos não foram registrados os sinais das outras 4 explosões. Além disso, as explosões registradas não apresentaram sinais claros em toda a extensão da linha. Devido a tudo isso e considerando as unidades geológicas presentes na região de estudo são sugeridos três modelos de velocidades sísmicas. O primeiro modelo refere-se ao tiro direto (EX31) localizado no extremo sudoeste da linha, sobre a Bacia do Paraná. Para este modelo obteve-se para superfície (0 km) a velocidade inicial de 2 km/s (coberturas); para a profundidade de 0,086 km a velocidade inicial é de 5,15 km/s (basalto); para a profundidade de 0,350 km obteve-se a velocidade inicial de 4,6 km/s (arenito - camada de baixa velocidade); para profundidade de 0,650 km a velocidade inicial é de 5,75 km/s e para profundidade de 4 km obteve-se a velocidade inicial de 6,07 km/s. O segundo modelo refere-se ao tiro reverso (EX34) localizado no centro da linha sobre granitóides do Grupo Araxá. Para este modelo obteve -se para superfície (0 km) a velocidade inicial de 2 km/s; para a profundidade de 0,06 km a velocidade inicial é de 5,69 km/s e para a profundidade de 0,860 km obteve-se a velocidade inicial de 6,25 km/s. Finalmente, o terceiro modelo refere-se ao tiro direto para toda a extensão da linha (300 km). Este modelo foi definido a partir de fases secundárias lidas nos registros e modelos anteriores propostos na literatura. Da superfície até os 4 km iniciais de profundidade este modelo é igual ao primeiro, para uma profundidade de 20 km obteve-se a velocidade inicial de 6,70 km/s e para uma profundidade de 40 km a velocidade é de 8,00 km/s (descontinuidade de MOHO). / This work to fulfil the degree of Master of Sciences is inserted among the deep seismic refraction studies of the Thematic Project 'Geophysical Studies and Tectonic Model of the Tocantins Province Central and Southeast Sectors, Central Brazil'. Three refraction lines, of around 300 km long each, were deployed, two of them in the Central sector and the other in the SE sector, that is subject of the present work. The equipment used in this experiment was composed by 111 SGR digital seismographs belonging to the PASSCAL Program. Complemented with auxiliary instruments from USGS and 13 seismographs belonging to IAG/USP. The space among the recording points was 2.5 km, which were located along main and secondary roads. Every 50 km was fired an explosion with 1000 kg of emulsion in each extreme and 500 kg in the central point. The geographical co-ordinates were determined by using the GPS differential method. The main objective of this work is to obtain as a final product a seismic velocity model with the physical characteristics of the main discontinuities in the crust and upper mantle. The packages SAC, SU and SEIS were used to perform the data analysis and processing. To carry on the modelling were used the ray theory and the synthetic seismograms construction, belonging to the SEIS package Data from the extreme and middle points of the seismic line were used to elaborate the final model, considering that due to technical problems signals from the other four explosions were not recorded. Apart from that, the recorded explosions did not present clear signals all along the extension of the line. Due to these facts, and considering also the geological units present in the studied region, are suggested three seismic velocity models. The first model is referred to the direct shot (EX31), which is localised in the Southwest extreme of the line on the Parana Basin province. In this model we obtained the P wave velocity (VP) of 2 km/sec at the surface, corresponding to the unconsolidated sediments and soil on the top of that basin. At a depth of 86 m we found VP of 5,15 km/sec and at a depth of 350 m the velocity VP of 4,6 km/sec, corresponding to the basalt and sand layers of the Parana Basin. Underlying them, at 650 m of depth we found the basement with VP of 5,75 km/sec and finally at a depth of 4 km there is a layer with VP of 6,07 km/sec, corresponding to a typical upper crust P wave velocity. The second model corresponds to the reverse shot (EX34) that is localised in the middle point of the line on the granitoides of the Araxa Group. For this model we obtained VP of 2 km/sec for the superficial layers, then at a depth of 60 m was obtained V P of 5,69 km/sec and for a depth of 860 m the value of V P is 6,25 km/sec. Finally, the third model belongs to the whole line section (300 km) from the direct shot (EX31). This model was obtained by using the arrivals of secondary phases and the results of models proposed in other works. From the surface down to 4 km of depth this model is similar to the first one. At 20 km of depth there is a layer with VP of 6,70 km/sec, corresponding to the lower crust, with Moho at a depth of 40 km with VP of 8,00 km/sec.