首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 406 毫秒
1.
Experimental measurements of fracture-induced seismic waves velocity variations at frequencies ~ 1 kHz, ~ 40 kHz and ~ 1 MHz were performed directly in the field at the rocky outcrop and in the laboratory on specific rock samples collected from the outcrops. The peridotite–lherzolite outcrop appeared macroscopically uniform and contained three systems of visible parallel sub-vertical fractures. This rock has substantial bulk density and higher than average value of seismic wave velocity. The presence of fracture systems gives rise to its velocity anisotropy. The seismic waves passing through the rock fractures are subject to velocity dispersion and frequency dependent attenuation. Our data, obtained from field and laboratory measurements, were compared with theoretical model predictions. In this model we successfully used displacement discontinuity approach. For the velocity dispersion evaluation we used multi-frequency measurements. The a priori observation of orientations and densities of fracture sets allowed evaluation of their stiffness. Our approach revealed that the first arrivals of seismic waves can be used for evaluation of P-wave group velocities, the specific case, in which we expect anomalous velocity dispersion. Our observations contribute to the issue of up-scaling of well-log derived velocities in fractured rock to the scale of standard seismic exploration frequencies.  相似文献   

2.
We use geodynamic models with imposed plate velocities to test the forward-modeled history of subduction based on a particular plate motion model against alternative seismic tomography models. We utilize three alternative published reference frames: a hybrid moving hotspot-palaeomagnetic, a hybrid moving hotspot-true polar wander corrected-palaeomagnetic, and a Subduction Reference Frame, a plate model including longitudinal shifts of subduction zones by matching subduction volumes imaged by P-wave tomography, to assess which model best predicts present day mantle structure compared with seismic tomography and volumetrically derived subduction history. Geodynamic modeling suggests paleo-longitudinal corrections applied to the Subduction Reference Frame result in lower mantle slab material beneath North America and East Asia accumulating up to 10–15° westward of that imaged by tomography, whereas the hybrid models develop material offset by 2–9°. However, the Subduction Reference Frame geodynamic model produces slab material beneath the Tethyan Domain coinciding with slab volumes imaged by tomography, whereas the hybrid reference frame models do not, suggesting regional paleo-longitudinal corrections are required to constrain slab locations. We use our models to test inferred slab sinking rates in the mantle focusing on well-constrained regions. We derive a globally averaged slab-sinking rate of 13 ± 3 mm/yr by combining the ages of onset and cessation of subduction from geological data and kinematic reconstructions with images of subducted slabs in the mantle. Our global average slab-sinking rate overlaps with the 15–20 mm/yr rate implied by mantle convection models using a lower mantle viscosity 100 times higher than the upper mantle.  相似文献   

3.
The Forsmark area belongs to the Paleoproterozoic Svecokarelian orogen (c. 1.9–1.8 Ga), the principal geological entity inside the Fennoscandian Shield, and is the site where Sweden has proposed to store its spent nuclear fuel. Three major sub-vertical (at the surface), composite ductile and brittle deformation zones that strike in a WNW or NW direction are present in the area. In between these zones the bedrock is less deformed and considered suitable for a repository. We present reprocessed reflection seismic data from seven profiles in which we have focused on improving the images in the depth range 1–5 km by passing lower frequencies through the processing flow at the cost of poorer resolution in the near-surface realm. The new images indicate that sub-horizontal to moderately dipping structures are possibly more extensive at depth than previously thought. Three main deeper reflective zones have been identified, one that is sub-horizontal and two that dip moderately to the southwest. The sub-horizontal reflective zone may represent a 1.27–1.26 Ga dolerite sill at about 3 km depth. One of the moderately dipping reflective zones may originate either from another dolerite sill or from a brittle fault system. The other moderately dipping structure may be present throughout most of the area and could cut all three sub-vertical deformation zones at depth. The new images and corresponding interpretation do not require a re-evaluation of the Forsmark site for storage of spent fuel, but they do influence how to interpret the deeper structures and, as a consequence, the tectonic evolution of the area.  相似文献   

4.
Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average Vp/Vs ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.  相似文献   

5.
Cross-hole imaging method using Time Domain (TD) and Frequency Domain (FD) parts of cross-hole radar tomography data acquired using Step Frequency Ground Penetrating Radar (SFGPR) was implemented. This method was adopted for imaging foundation of a dam to check if the foundation was free of geological weak zones. The dam site is characterised by massive and jointed-phyllites associated with major and minor shears. The cross-hole radar tomography data was acquired in the frequency bandwidth of 250 MHz, from the deepest level gallery up to a depth of 40 m in the foundation. In TD, first arrival time and amplitudes of radio waves were inverted using Simultaneous Iterative Reconstruction Technique (SIRT) resulting in velocity and attenuation tomograms. The tomograms showed nearly uniform velocity or attenuation structure in the respective tomographic plane. Subsequently, cross-hole radar tomography data was analysed in FD for a variation of spectrum-amplitude at different frequencies. Amplitudes picked at each single frequency were then inverted using SIRT for obtaining frequency domain attenuation tomogram (FDAT). The FDAT clearly showed presence of anomalous high attenuation zones in the depth range of 23–33 m of the tomographic plane. The anomalous zones in the attenuation tomogram are weak zones in the foundation. To validate the above observations, cross-hole seismic tomography was also done in the same boreholes. Cross-hole seismic tomography results showed low velocity (p-wave) zones around the same location corresponding to the high attenuation zone in FDAT, bringing the dormant weak zone to light. This enabled fine-tuning of the reinforcement design and strengthening the weak zone. This paper discusses the cross-hole radar tomography imaging method, the results of its application in imaging weak zones in the foundation and the comparison of cross-hole radar tomography results (in TD and FD) with the cross-hole seismic tomography results.  相似文献   

6.
Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic VP/VS ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. VP varies between 8.9 km s? 1 and 5.6 km s? 1 (46% anisotropy), and 8.3 km s? 1 and 5.8 km s? 1 (37%), and VS between 5.1 km s? 1 and 2.5 km s? 1 (66%), and 4.7 km s? 1 and 2.9 km s? 1 (50%) for the single-crystal and aggregate, respectively. The VP/VS ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high VP/VS ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s?1 for VP and VS, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.  相似文献   

7.
We present a new method for the prediction of the discontinuities and lithological variations ahead of the tunnel face. The automatic procedure is applied to data collected by seismic reflection surveys, with the sources and sensors located along the tunnel. The method allows: i) to estimate an average value of the wave velocity; ii) to detect the discontinuities for each source point; and iii) to analyze and plot the number of superposing estimates for each node of the domain. The final result can be interpreted as the probability to detect a discontinuity at a certain distance from the tunnel face. The method automatically estimates the peaks in the seismograms that can be related to a reflection. On the base of this process, the method only requires the source–receiver geometry and the data acquisition parameters. The procedure has been tested on synthetic and real data coming from a seismic survey on a tunnel under construction. The results indicate that the method runs very fast and it is reliable in the identification of lithological changes and discontinuities, up to more than 100 m ahead of the tunnel face.  相似文献   

8.
Seismic analysis and geochemical interpretations provide evidence that two separate hydrothermal cells circulate within the greater Lassen hydrothermal system. One cell originates south to SW of Lassen Peak and within the Brokeoff Volcano depression where it forms a reservoir of hot fluid (235–270 °C) that boils to feed steam to the high-temperature fumarolic areas, and has a plume of degassed reservoir liquid that flows southward to emerge at Growler and Morgan Hot Springs. The second cell originates SSE to SE of Lassen Peak and flows southeastward along inferred faults of the Walker Lane belt (WLB) where it forms a reservoir of hot fluid (220–240 °C) that boils beneath Devils Kitchen and Boiling Springs Lake, and has an outflow plume of degassed liquid that boils again beneath Terminal Geyser. Three distinct seismogenic zones (identified as the West, Middle, and East seismic clusters) occur at shallow depths (< 6 km) in Lassen Volcanic National Park, SW to SSE of Lassen Peak and adjacent to areas of high-temperature (≤ 161 °C) fumarolic activity (Sulphur Works, Pilot Pinnacle, Little Hot Springs Valley, and Bumpass Hell) and an area of cold, weak gas emissions (Cold Boiling Lake). The three zones are located within the inferred Rockland caldera in response to interactions between deeply circulating meteoric water and hot brittle rock that overlies residual magma associated with the Lassen Volcanic Center. Earthquake focal mechanisms and stress inversions indicate primarily N–S oriented normal faulting and E–W extension, with some oblique faulting and right lateral shear in the East cluster. The different focal mechanisms as well as spatial and temporal earthquake patterns for the East cluster indicate a greater influence by regional tectonics and inferred faults within the WLB. A fourth, deeper (5–10 km) seismogenic zone (the Devils Kitchen seismic cluster) occurs SE of the East cluster and trends NNW from Sifford Mountain toward the Devils Kitchen thermal area where fumarolic temperatures are ≤ 123 °C. Lassen fumaroles discharge geothermal gases that indicate mixing between a N2-rich, arc-type component and gases derived from air-saturated meteoric recharge water. Most gases have relatively weak isotopic indicators of upper mantle or volcanic components, except for gas from Sulphur Works where δ13C–CO2, δ34S–H2S, and δ15N–N2 values indicate a contribution from the mantle and a subducted sediment source in an arc volcanic setting.  相似文献   

9.
Reflection seismic data were acquired along a c. 23 km long profile over the Pärvie Fault system with a nominal receiver and source spacing of 20 m. An hydraulic breaking hammer was used as a source, generating signals with a penetration depth of about 5–6 km. Steeply dipping reflections from the end-glacial faults are observed, as well as sub-horizontal reflections. The location and orientation of the reflections from the faults agree well with surface geological observations of fault geometries. Reflections from a potential fourth end-glacial fault is observed further to the east along the profile. The more sub-horizontal reflections may originate from gabbroic bodies within the granitic basement or from deeper lying greenstones. Our results indicate that the end-glacial faults dip at moderate to steep dips down to at least 2–3 km depth, and possibly continue at this dip to depths of 6 km. This result has significant implications for determining the state of stress required to activate the faults in the past and in the future.  相似文献   

10.
We have studied the dependency between incoming plate structure, bending-related faulting, lithospheric hydration, and outer rise seismic activity offshore Maule, Chile. We derived a 2D Poisson's ratio distribution from P- and S-wave seismic wide angle data collected in the trench-outer rise. High values of Poisson's ratio in the uppermost mantle suggest that the oceanic lithosphere is highly hydrated due to the water infiltration through bending-related normal faults outcropping at the seafloor. This process is presumably facilitated by the presence of a seamount in the area. We conclude that water infiltrates deep into the lithosphere, when it approaches the Chile trench, producing a reduction of crustal and upper mantle velocities, supporting serpentinization of the upper mantle. Further, we observed a mantle Vp anisotropy of 8%, with the fast velocity axis running normal to the abyssal hill fabric and hence in spreading direction, indicating that outer rise processes have yet not affected anisotropy.The first weeks following the megatrust Mw = 8.8 Maule earthquake in 2010 were characterized by a sudden increase of the outer rise seismic activity, located between 34° S and 35°30′ S. We concluded that this phenomenon is a result of an intensification of the water infiltration process in the outer rise, presumably triggered by the main shock, whose epicenter was located some 100 km to the south east of the cluster.  相似文献   

11.
Many authors have proposed that the study of seismicity rates is an appropriate technique for evaluating how close a seismic gap may be to rupture. We designed an algorithm for identification of patterns of significant seismic quiescence by using the definition of seismic quiescence proposed by Schreider (1990). This algorithm shows the area of quiescence where an earthquake of great magnitude may probably occur. We have applied our algorithm to the earthquake catalog on the Mexican Pacific coast located between 14 and 21 degrees of North latitude and 94 and 106 degrees West longitude; with depths less than or equal to 60 km and magnitude greater than or equal to 4.3, which occurred from January, 1965 until December, 2014. We have found significant patterns of seismic quietude before the earthquakes of Oaxaca (November 1978, Mw = 7.8), Petatlán (March 1979, Mw = 7.6), Michoacán (September 1985, Mw = 8.0, and Mw = 7.6) and Colima (October 1995, Mw = 8.0). Fortunately, in this century earthquakes of great magnitude have not occurred in Mexico. However, we have identified well-defined seismic quiescences in the Guerrero seismic-gap, which are apparently correlated with the occurrence of silent earthquakes in 2002, 2006 and 2010 recently discovered by GPS technology.  相似文献   

12.
The strong earthquake (M = 7) that occurred in the Fucino basin (central Italy) on January 13, 1915 was followed by six earthquakes of M > 5.5 and several other shocks of M > 5 in the major seismic zones of the northern Apennines from 1916 to 1920. This seismicity pattern is consistent with the implications of the present tectonic setting in the study area, which suggests that strong decoupling earthquakes in the central Apennines cause a significant increase of tectonic load, and possibly of seismicity, in the northern Apennines. A numerical simulation, carried out by an elastic-viscous model, of the stress diffusion induced by the Fucino and successive largest earthquakes, shows that each of the above shocks occurred when the respective zone was reached by the highest values of the strain and strain rate perturbation triggered by the previous events. Furthermore, the computed strain regime at each earthquake site is consistent with the known faulting pattern. The results provide important insights into the physical mechanism that controls the interaction of seismic sources in the central and northern Apennines.  相似文献   

13.
Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m3) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m3. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.  相似文献   

14.
Although offset and age data from displaced landforms are essential for identifying earthquake clusters and thus testing whether faults slip at uniform or secularly varying rates, it is not clear how the uncertainties in such measurements should be propagated so as to yield a robust fault-slip history (i.e., record of fault displacement over time). Here we develop a Monte Carlo approach for estimating the distribution of geologically reasonable fault-slip histories that fit the offset and age data from a population of dated and displaced landforms. The model assumes that the landforms share common faulting histories, the offset and age constraints are correct, and the fault has not reversed shear sense. Analysis of the model results yields both a precise average slip rate, in the case where a linear fit is applied to the data, and a best-fit fault-slip history, in the case where the linear constraint is removed. The method can be used to test for secular variation in slip because the uncertainty on this best-fit history is quantified. By applying the method to previously published morphochronologic data from faulted late Quaternary terrace risers along the Kunlun fault in China and the Awatere fault in New Zealand, we have assessed the extent to which our modeled average slip rates match previously reported values and the data support previous interpretations of uniform slip rate. The Kunlun data set yields average slip rates of 8.7 + 3.6/?2.1 mm/yr and 5.1 + 1.6/?1.2 mm/yr (68.27% confidence), for the central and eastern reaches of the fault, respectively, both of which match previously published slip rates. Our analysis further indicates that these fault reaches have both slipped uniformly over the latest Quaternary. In contrast, analysis of data from the Saxton River site along the Awatere fault reveals a mid-Holocene deceleration in slip rate from 6.2 + 1.6/?1.4 mm/yr to 2.8 + 1.0/?0.6 mm/yr. This result contradicts previous interpretations of uniform slip along the Awatere fault. The Monte Carlo method we present here for quantifying fault-slip histories using the offset and age data from a population of faulted landforms provides an important tool for distinguishing temporally uniform from secularly varying fault slip.  相似文献   

15.
Airborne electromagnetic (AEM) surveys, when regionally extensive, may sample a wide-range of geological formations. The majority of AEM surveys can provide estimates of apparent (half-space) conductivity and such derived data provide a mapping capability. Depth discrimination of the geophysical mapping information is controlled by the bandwidth of each particular system. The objective of this study is to assess the geological information contained in accumulated frequency-domain AEM survey data from the UK where existing geological mapping can be considered well-established. The methodology adopted involves a simple GIS-based, spatial join of AEM and geological databases. A lithology-based classification of bedrock is used to provide an inherent association with the petrophysical rock parameters controlling bulk conductivity. At a scale of 1:625k, the UK digital bedrock geological lexicon comprises just 86 lithological classifications compared with 244 standard lithostratigraphic assignments. The lowest common AEM survey frequency of 3 kHz is found to provide an 87% coverage (by area) of the UK formations. The conductivities of the unsampled classes have been assigned on the basis of inherent lithological associations between formations. The statistical analysis conducted uses over 8 M conductivity estimates and provides a new UK national scale digital map of near-surface bedrock conductivity. The new baseline map, formed from central moments of the statistical distributions, allows assessments/interpretations of data exhibiting departures from the norm. The digital conductivity map developed here is believed to be the first such UK geophysical map compilation for over 75 years. The methodology described can also be applied to many existing AEM data sets.  相似文献   

16.
Geodetically-detected episodes of slow slip appear in several subduction zones to be accompanied by bursts of low-frequency coherent noise known as seismic tremor, but whether a single physical process governs this association or even whether slow slip is invariably accompanied by tremor remains unresolved. Detailed analysis of broadband seismic data spanning a slow slip episode in the Hikurangi subduction zone, New Zealand, reveals that slow slip was accompanied by distinct reverse-faulting microearthquakes, rather than tremor. The timing, location, and faulting style of these earthquakes are consistent with stress triggering down-dip of the slow slip patch, either on the subduction interface or just below it. These results indicate that tremor is not ubiquitous during subduction zone slow slip, and that slow slip in subduction zone environments is capable of triggering high-frequency earthquakes near the base of the locked subduction thrust. In this and other locations (Hawaii, Boso Peninsula) where slow slip is accompanied by triggered microseismicity, the estimated upper extent of the slow slip is shallower (less than ~ 20 km) than in those locations from which tremor has been reported. This suggests that ambient temperature- or pressure-dependent factors govern the character of the seismic response to slow slip on subduction thrusts and other large faults, with rheological or lithological conditions at shallow depths triggering high-frequency microearthquakes and those at greater depths triggering seismic tremor.  相似文献   

17.
We have studied the scaling behavior of compressional-wave velocity and density logs from an exploration borehole that extends down to about 700 m depth in the Brunswick No. 6 mining area, Bathurst Mining Camp, Canada. Using statistical methods, vertical and horizontal scale lengths of heterogeneity were estimated. Vertical scale length estimates from the velocity, density and calculated acoustic impedance are 14 m, 33 m, and about 20 m, respectively. Although the estimated scale length for the acoustic impedance implies a weak scattering environment, elastic finite difference modeling of seismic wave propagation in 2D heterogeneous media demonstrates that even this weak scattering medium can mask seismic signals from small, but yet economically feasible, massive sulfide deposits. Further analysis of the synthetic seismic data suggests that in the presence of heterogeneity, lenticular-shaped targets may only exhibit incomplete diffraction signals whereby the down-dip tails of these diffractions are mainly visible on the stacked sections. Therefore, identification of orebody generated diffractions is much easier on the unmigrated stacked sections than on migrated stacked sections. The numerical seismic modeling in 2D heterogeneous media indicates that in the presence of large horizontal, but small vertical scale lengths (structural anisotropy), identification of massive sulfide deposits is possible, but their delineation at depth requires detailed velocity modeling and processing algorithms which can handle the anisotropy.  相似文献   

18.
The spontaneous imbibition of water and other liquids into gas-filled fractures in variably-saturated porous media is important in a variety of engineering and geological contexts. However, surprisingly few studies have investigated this phenomenon. We present a theoretical framework for predicting the 1-dimensional movement of water into air-filled fractures within a porous medium based on early-time capillary dynamics and spreading over the rough surfaces of fracture faces. The theory permits estimation of sorptivity values for the matrix and fracture zone, as well as a dispersion parameter which quantifies the extent of spreading of the wetting front. Quantitative data on spontaneous imbibition of water in unsaturated Berea sandstone cores were acquired to evaluate the proposed model. The cores with different permeability classes ranging from 50 to 500 mD and were fractured using the Brazilian method. Spontaneous imbibition in the fractured cores was measured by dynamic neutron radiography at the Neutron Imaging Prototype Facility (beam line CG-1D, HFIR), Oak Ridge National Laboratory. Water uptake into both the matrix and the fracture zone exhibited square-root-of-time behavior. The matrix sorptivities ranged from 2.9 to 4.6 mm s−0.5, and increased linearly as the permeability class increased. The sorptivities of the fracture zones ranged from 17.9 to 27.1 mm s−0.5, and increased linearly with increasing fracture aperture width. The dispersion coefficients ranged from 23.7 to 66.7 mm2 s−1 and increased linearly with increasing fracture aperture width and damage zone width. Both theory and observations indicate that fractures can significantly increase spontaneous imbibition in unsaturated sedimentary rock by capillary action and surface spreading on rough fracture faces. Fractures also increase the dispersion of the wetting front. Further research is needed to investigate this phenomenon in other natural and engineered porous media.  相似文献   

19.
The location in the Barberton Greenstone Belt (Kaapvaal Craton) of ∼3.26–3.24 Ga asteroid impact ejecta units at, and immediately above, a sharp break between a > 12 km-thick mafic–ultramafic volcanic crust (Onverwacht Group ∼3.55–3.26 Ga, including the ∼3.298 > 3.258 Ga Mendon Formation) and a turbidite–felsic volcanic rift-facies association (Fig Tree Group ∼3.258–3.225 Ga), potentially represents the first documented example of cause–effect relations between extraterrestrial bombardment and major tectonic and igneous events [D.R. Lowe, G.R. Byerly, F. Asaro, F.T. Kyte, Geological and geochemical record of 3400 Ma old terrestrial meteorite impacts, Science 245 (1989) 959–962; D.R. Lowe, G.R. Byerly, F.T. Kyte, A. Shukolyukov, F. Asaro, A. Krull, Spherule beds 3.47–3.34 Ga-old in the Barberton greenstone belt, South Africa: a record of large meteorite impacts and their influence on early crustal and biological evolution, Astrobiology 3 (2003) 7–48; A.Y. Glikson, The astronomical connection of terrestrial evolution: crustal effects of post-3.8 Ga mega-impact clusters and evidence for major 3.2 ± 0.1 Ga bombardment of the Earth–Moon system, J. Geodyn. 32 (2001) 205–229]. Here we correlate this boundary with a contemporaneous break and peak magmatic and faulting events in the Pilbara Craton, represented by the truncation of a 3.255–3.235 Ga-old volcanic sequence (Sulphur Springs Group—SSG) by a turbidite-banded iron formation–felsic volcanic association (Pincunah Hill Formation, basal Gorge Creek Group). These events are accompanied by ∼3.252–3.235 Ga granitoids (Cleland plutonic suite). The top of the komatiite–tholeiite–rhyolite sequence of the SSG is associated with a marker chert defined at 3.238 ± 3–3.235 ± 3 Ga, abruptly overlain by an olistostrome consisting of mega-clasts of felsic volcanics, chert and siltstone up to 250 × 150 m-large, intercalated with siliciclastic sedimentary rocks and felsic volcanics (Pincunah Hill Formation-basal Gorge Creek Group-GCG [R. M. Hill, Stratigraphy, structure and alteration of hanging wall sedimentary rocks at the Sulphur Springs volcanogenic massive sulphide (VMS) prospect, east Pilbara Craton, Western Australia. B.Sc Hon. Thesis, University of Western Australia (1997) 67 pp.; M.J. Van Kranendonk, A.H. Hickman, R.H. Smithies, D.R. Nelson, Geology and tectonic evolution of the Archaean north Pilbara terrain, Pilbara Craton, Western Australia, Econ. Geol. 97 (2002) 695–732; M.J. Van Kranendonk, Geology of the North Shaw 1 : 100 000 Sheet. Geological Survey Western Australia 1 : 100 000 Geological Series (2000) 86 pp., R. Buick, C.A.W. Brauhart, P. Morant, J.R. Thornett, J.G. Maniew, J.G. Archibald, M.G. Doepel, I.R. Fletcher, A.L. Pickard, J.B. Smith, M.B. Barley, N.J. McNaughton, D.I. Groves, Geochronology and stratigraphic relations of the Sulphur Springs Group and Strelley Granite: a temporally distinct igneous province in the Archaean Pilbara Craton, Australia, Precambrian Res. 114 (2002) 87–120]). The structure and scale of the olistostrome, not seen elsewhere in the Pilbara Craton, is interpreted in terms of intense faulting and rifting, supported by topographic relief represented by deep incision of overlying arenites (Corboy Formation) into underlying units [M.J. Van Kranendonk, Geology of the North Shaw 1 : 100 000 Sheet. Geological Survey Western Australia 1 : 100 000 Geological Series (2000) 86 pp.]. The age overlaps between (1) 3.255 ± 4–3.235 ± 3 Ga peak igneous activity represented by the SSG and the Cleland plutonic suite (Pilbara Craton) and the 3.258 ± 3 Ga S2 Barberton impact unit, and (2) 3.235 ± 3 Ga top SSG break and associated faulting and the 3.243 ± 4 S3–S4 Barberton impact units may not be accidental. Should correlations between the Barberton S2–S4 impact units and magmatic and tectonic events in the Pilbara Craton be confirmed, they would imply impact-triggered reactivation of mantle convection, crustal anatexis, faulting and strong vertical movements in Archaean granite–greenstone terrains associated with large asteroid impacts, culminating in transformation from sima-dominated crust to continental rift environments.  相似文献   

20.
We determined crustal structure along the latitude 30°N through the eastern Tibetan Plateau using a teleseismic receiver function analysis. The data came mostly from seismic stations deployed in eastern Tibet and western Sichuan region from 2004 to 2006. Crustal thickness and Vp/Vs ratio at each station were estimated by the H–k stacking method. On the profile, the mean crustal thickness and Vp/Vs ratio were found to be 62.3 km and 1.74 in the Lhasa block, 71.2 km and 1.79 near the Bangong–Nujiang suture, 66.3 km and 1.80 in the Qiangtang block, 59.8 km and 1.81 in the Songpan–Garze block, and 42.9 km and 1.76 in the Yangtze block, respectively. The estimated crustal thicknesses are consistent with predictions based on the topography and the Airy isostasy, except near the Bangong–Nujiang suture and in the Qiangtang block where the crust is 5–10 km thicker than predicted, indicating that the crust may be denser, possibly due to mafic underplating. We also inverted receiver functions for crustal velocity structure along the profile, which reveals a low S-wave velocity zone in the lower crust beneath the eastern Tibetan Plateau, although the extent of the low-velocity zone varies considerably. The low-velocity zone, together with previous results, suggests limited partial melting and localized crustal flow in the lower crust of the eastern Tibetan Plateau.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号