Seismic velocity models for heat zones in Athabasca tar sands

Geophysics ◽  
1990 ◽  
Vol 55 (8) ◽  
pp. 1108-1112 ◽  
Author(s):  
Larry R. Lines ◽  
Ronald Jackson ◽  
James D. Covey

Recent laboratory and field studies indicate that the P-wave velocity in Athabasca tar sands decreases when temperature increases during steam injection. In this paper we derive time variant velocity models from seismic traveltime inversions of both reflection and borehole data. Prior to steam injection, three‐dimensional (3-D) reflector velocity‐depth models are established using image‐ray conversions of traveltimes to depth. The changes in velocity due to steam injection are modeled by inverting traveltime data from seismic monitor surveys after steam injection and comparing these results to velocities computed prior to steam injection. Velocity models are essentially determined by traveltimes from the 3-D seismic reflection survey. The surface‐to‐wellbore data traveltimes show the expected delay caused by steam injection but do not significantly alter the velocity model produced by reflection traveltimes. For seismic monitor surveys, low‐velocity zones show a very good correlation with zones of temperature increase at injector well positions. The results indicate that velocity models obtained from seismic traveltimes may prove useful in detecting steam fronts in tar sands.

2011 ◽  
Vol 2011 ◽  
pp. 1-6 ◽  
Author(s):  
Haiou Li ◽  
Xiwei Xu ◽  
Wentao Ma ◽  
Ronghua Xie ◽  
Jingli Yuan ◽  
...  

Three-dimensional P wave velocity models under the Zipingpu reservoir in Longmenshan fault zone are obtained with a resolution of 2 km in the horizontal direction and 1 km in depth. We used a total of 8589 P wave arrival times from 1014 local earthquakes recorded by both the Zipingpu reservoir network and temporary stations deployed in the area. The 3-D velocity images at shallow depth show the low-velocity regions have strong correlation with the surface trace of the Zipingpu reservoir. According to the extension of those low-velocity regions, the infiltration depth directly from the Zipingpu reservoir itself is limited to 3.5 km depth, while the infiltration depth downwards along the Beichuan-Yingxiu fault in the study area is about 5.5 km depth. Results show the low-velocity region in the east part of the study area is related to the Proterozoic sedimentary rocks. The Guanxian-Anxian fault is well delineated by obvious velocity contrast and may mark the border between the Tibetan Plateau in the west and the Sichuan basin in the east.


2019 ◽  
Vol 23 (3) ◽  
pp. 209-223 ◽  
Author(s):  
Caglar Ozer ◽  
Mehmet Ozyazicioglu

Erzurum and its surroundings are one of the seismically active and hydrothermal areas in the Eastern part of Turkey. This study is the first approach to characterize the crust by seismic features by using the local earthquake tomography method. The earthquake source location and the three dimensional seismic velocity structures are solved simultaneously by an iterative tomographic algorithm, LOTOS-12. Data from a combined permanent network comprising comprises of 59 seismometers which was installed by Ataturk University-Earthquake Research Center and Earthquake Department of the Disaster and Emergency Management Authority  to monitor the seismic activity in the Eastern Anatolia, In this paper, three-dimensional Vp and Vp/Vs characteristics of Erzurum geothermal area were investigated down to 30 km by using 1685 well-located earthquakes with 29.894 arrival times, consisting of 17.298 P- wave and 12.596 S- wave arrivals. We develop new high-resolution depth-cross sections through Erzurum and its surroundings to provide the subsurface geological structure of seismogenic layers and geothermal areas. We applied various size horizontal and vertical checkerboard resolution tests to determine the quality of our inversion process. The basin models are traceable down to 3 km depth, in terms of P-wave velocity models. The higher P-wave velocity areas in surface layers are related to the metamorphic and magmatic compact materials. We report that the low Vp and high Vp/Vs values are observed in Yedisu, Kaynarpinar, Askale, Cimenozu, Kaplica, Ovacik, Yigitler, E part of Icmeler, Koprukoy, Uzunahmet, Budakli, Soylemez, Koprukoy, Gunduzu, Karayazi, Icmesu, E part of Horasan and Kaynak regions indicated geothermal reservoir.


Geology ◽  
2019 ◽  
Vol 48 (3) ◽  
pp. 231-235 ◽  
Author(s):  
B.G. McVey ◽  
E.E.E. Hooft ◽  
B.A. Heath ◽  
D.R. Toomey ◽  
M. Paulatto ◽  
...  

Abstract Despite multidisciplinary evidence for crustal magma accumulation below Santorini volcano, Greece, the structure and melt content of the shallow magmatic system remain poorly constrained. We use three-dimensional (3-D) velocity models from tomographic inversions of active-source seismic P-wave travel times to identify a pronounced low-velocity anomaly (–21%) from 2.8 km to 5 km depth localized below the northern caldera basin. This anomaly is consistent with depth estimates of pre-eruptive storage and a recent inflation episode, supporting the interpretation of a shallow magma body that causes seismic attenuation and ray bending. A suite of synthetic tests shows that the geometry is well recovered while a range of melt contents (4%–13% to fully molten) are allowable. A thin mush region (2%–7% to 3%–10% melt) extends from the main magma body toward the northeast, observed as low velocities confined by tectono-magmatic lineaments. This anomaly terminates northwest of Kolumbo; little to no melt underlies the seamount from 3 to 5 km depth. These structural constraints suggest that crustal extension and edifice loads control the geometry of magma accumulation and emphasize that the shallow crust remains conducive to melt storage shortly after a caldera-forming eruption.


2014 ◽  
Vol 51 (4) ◽  
pp. 358-372 ◽  
Author(s):  
Draga Talinga ◽  
Andrew J. Calvert

Across the Nechako–Chilcotin plateau of British Columbia, the distribution of Cretaceous sedimentary rocks, which are considered prospective for hydrocarbon exploration, is poorly known due to the surface cover of glacial deposits and Tertiary volcanic rocks. To constrain the subsurface distribution of these Cretaceous rocks, in 2008 Geoscience BC acquired seven long, up to 14.4 km, offset vibroseis seismic reflection lines across a north-northwest-trending belt of exhumed sedimentary rocks inferred to be part of the Taylor Creek Group. P-wave velocity models, which are consistent with sonic logs from nearby wells, have been estimated using three-dimensional first-arrival tomography to depths ranging from 1 to 4 km. Igneous basement can be identified on most lines using the 5.5 km/s isovelocity contour, which locates the top of the basement to an accuracy of ∼400 m where its depth is known in exploration wells. There is no general distinction on the basis of seismic velocity between Cretaceous sedimentary and Paleocene–Eocene volcanic–volcaniclastic rocks, both of which appear to be characterized in the tomographic models by velocities of 3.0–5.0 km/s. The geometry of the igneous basement inferred from the velocity models identifies north-trending basins and ridges, which correlate with exposed rocks of the Jurassic Hazelton Group. Identified Cretaceous sedimentary rocks occur beneath less negative Bouguer gravity anomalies, but the original distribution of these rocks has been disrupted by later Tertiary extension that created north-trending basins associated with the most negative gravity anomalies. We suggest that Cretaceous sedimentary rocks, if deposited, could be preserved within these basins if the rocks had not been eroded prior to Tertiary extension.


Author(s):  
Fumiaki Nagashima ◽  
Hiroshi Kawase

Summary P-wave velocity (Vp) is an important parameter for constructing seismic velocity models of the subsurface structures by using microtremors and earthquake ground motions or any other geophysical exploration data. In order to reflect the ground survey information in Japan to the Vp structure, we investigated the relationships among Vs, Vp, and depth by using PS-logging data at all K-NET and KiK-net sites. Vp values are concentrated at around 500 m/s and 1,500 m/s when Vs is lower than 1,000 m/s, where these concentrated areas show two distinctive characteristics of unsaturated and saturated soil, respectively. Many Vp values in the layer shallower than 4 m are around 500 m/s, which suggests the dominance of unsaturated soil, while many Vp values in the layer deeper than 4 m are larger than 1,500 m/s, which suggests the dominance of saturated soil there. We also investigated those relationships for different soil types at K-NET sites. Although each soil type has its own depth range, all soil types show similar relationships among Vs, Vp, and depth. Then, considering the depth profile of Vp, we divided the dataset into two by the depth, which is shallower or deeper than 4 m, and calculated the geometrical mean of Vp and the geometrical standard deviation in every Vs bins of 200 m/s. Finally, we obtained the regression curves for the average and standard deviation of Vp estimated from Vs to get the Vp conversion functions from Vs, which can be applied to a wide Vs range. We also obtained the regression curves for two datasets with Vp lower and higher than 1,200 m/s. These regression curves can be applied when the groundwater level is known. In addition, we obtained the regression curves for density from Vs or Vp. An example of the application for those relationships in the velocity inversion is shown.


1976 ◽  
Vol 66 (2) ◽  
pp. 501-524
Author(s):  
Keiiti Aki ◽  
Anders Christoffersson ◽  
Eystein S. Husebye

abstract Using P-wave residuals for teleseismic events observed at the Montana Large Aperture Seismic Array (LASA), we have determined the three-dimensional seismic structure of the lithosphere under the array to a depth of 140 km. The root-mean-square velocity fluctuation was found to be at least 3.2 per cent which may be compared to estimate of ca. 2 per cent based on the Chernov random medium theory. The solutions are given by both the generalized inverse and stochastic inverse methods in order to demonstrate the relative merit of different inversion techniques. The most conspicuous feature of the lithosphere under LASA is a low-velocity anomaly in the central and northeast part of the array siting area with the N60°E trend and persisting from the upper crust to depths greater than 100 km. We interpret this low-velocity anomaly as a zone of weakness caused by faulting and shearing associated with the building of the Rocky Mountains.


1997 ◽  
Vol 40 (1) ◽  
Author(s):  
G. Drakatos ◽  
G. Karantonis ◽  
G. N. Stavrakakis

The three-dimensional velocity structure of the crust in the Aegean sea and the surrounding regions (34.0º-42.OºN, 19.0ºE-29.0ºE) is investigated by inversion of about 10000 residuals of arrival times of P-wave from local events. The resulting velocity structure shows strong horizontal variations due to the complicated crustal structure and the variations of crustal thickness. The northern part of the region generally shows high velocities. In the inner part of the volcanic arc (Southern Aegean area), relatively low velocities are observed, suggesting a large-scale absorption of seismic energy as confirmed by the low seismicity of the region. A low velocity zone was observed along the subduction zone of the region, up to a depth of 4 km. The existence of such a zone could be due to granitic or other intrusions in the crust during the uplift of the region during Alpidic orogenesis.


2021 ◽  
Vol 873 (1) ◽  
pp. 012056
Author(s):  
M F R Auly ◽  
A K Ilahi ◽  
I Madrinovella ◽  
S Widyanti ◽  
S K Suhardja ◽  
...  

Abstract The tectonic setting of Java island, located at southwestern edge of the Eurasia continent, is dominated by the subduction of Indo-Australia plate. One of the characteristics of active subduction is active seismicity, the generation of arc magmatism and volcanic activity. Mt. Merapi is one example of active volcano related with the subduction process. It is one of the most active volcanoes with location close to high population area. To better understand this area, we employed the Receiver Function technique, a method to image sub surface structure by removing the vertical component from horizontal component. First, we collected high magnitude events and processed RF with water level deconvolution method. Then, we constructed synthetic model with initial velocity input from previous tomography model. Note that we used reflectivity method in generating synthetic model with input parameters matched with parameters from real data processing. Next, we adjusted velocity inputs mainly on tops sediments (1-3 km) to include sediment layers and volcanic rocks, mid-depth low velocity zone that may be related with magma chamber and depth of crust-mantle boundary. Current forward velocity models show a relatively good agreement from 3 stations (ME25, ME32 and ME36). We estimate a thin layer of sediments followed a zone of velocity layer at a depth of 10-15 km and crust-mantle boundary ranging from 26-29 km. In this study, simulated that the signal of sediments layer and low velocity layers interfere main crust mantle boundary that supposed to be highest signal after the P wave in the typical receiver function study.


2021 ◽  
pp. M56-2020-19
Author(s):  
E. R. Ivins ◽  
W. van der Wal ◽  
D. A. Wiens ◽  
A. J. Lloyd ◽  
L. Caron

AbstractThe Antarctic mantle and lithosphere are known to have large lateral contrasts in seismic velocity and tectonic history. These contrasts suggest differences in the response time scale of mantle flow across the continent, similar to those documented between the northeastern and southwestern upper mantle of North America. Glacial isostatic adjustment and geodynamical modeling rely on independent estimates of lateral variability in effective viscosity. Recent improvements in imaging techniques and the distribution of seismic stations now allow resolution of both lateral and vertical variability of seismic velocity, making detailed inferences about lateral viscosity variations possible. Geodetic and paleo sea-level investigations of Antarctica provide quantitative ways of independently assessing the three-dimensional mantle viscosity structure. While observational and causal connections between inferred lateral viscosity variability and seismic velocity changes are qualitatively reconciled, significant improvements in the quantitative relations between effective viscosity anomalies and those imaged by P- and S-wave tomography have remained elusive. Here we describe several methods for estimating effective viscosity from S-wave velocity. We then present and compare maps of the viscosity variability beneath Antarctica based on the recent S-wave velocity model ANT-20 using three different approaches.


Geophysics ◽  
2021 ◽  
pp. 1-35
Author(s):  
M. Javad Khoshnavaz

Building an accurate velocity model plays a vital role in routine seismic imaging workflows. Normal-moveout-based seismic velocity analysis is a popular method to make the velocity models. However, traditional velocity analysis methodologies are not generally capable of handling amplitude variations across moveout curves, specifically polarity reversals caused by amplitude-versus-offset anomalies. I present a normal-moveout-based velocity analysis approach that circumvents this shortcoming by modifying the conventional semblance function to include polarity and amplitude correction terms computed using correlation coefficients of seismic traces in the velocity analysis scanning window with a reference trace. Thus, the proposed workflow is suitable for any class of amplitude-versus-offset effects. The approach is demonstrated to four synthetic data examples of different conditions and a field data consisting a common-midpoint gather. Lateral resolution enhancement using the proposed workflow is evaluated by comparison between the results from the workflow and the results obtained by the application of conventional semblance and three semblance-based velocity analysis algorithms developed to circumvent the challenges associated with amplitude variations across moveout curves, caused by seismic attenuation and class II amplitude-versus-offset anomalies. According to the obtained results, the proposed workflow is superior to all the presented workflows in handling such anomalies.


Sign in / Sign up

Export Citation Format

Share Document