首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
Micropiping processes and biancana evolution in southeast Tuscany, Italy   总被引:1,自引:0,他引:1  
Biancane badlands consisting of small domes dissected by rills and micropipes, with rough disordered microrehef, can be found along the Apennines in Italy. The dominant processes forming biancane differ from those of badlands formed on smectite-rich mudrocks, as micropipes associated with pseudokarstic enlargement of pores and cracks predominate and form the main routes for evacuation of eroded material.Biancana evolution is controlled by water infiltration into intact bedrock, producing an erodible weathering ‘rind’ which is more porous than intact rock. This rind is easily removed by rill or micropipe flow, and erosion is therefore ‘weathering-controlled’, depending on rind production by infiltrating water. Infiltration is initially slow and stepped, due to slow water movement through very small capillary pores in intact rock alternating with rapid filling of macropores and cracks. This occurs due to rapid matrix pore enlargement by dispersion and/or dissolution. The infiltration pattern is accurately reproduced by a model built on progressive development of weathering layers by moisture penetration. Model results are consistent with weathering rind depths and erosion observed in the field, and show that a pipe network can be generated on newly exposed rock by the rainfall of one year.Propagation of the pipe network diverts a progressively larger proportion of runoff into micropipes, expanding weathering rind production within the biancana as well as on the surface. Internal weathering and flow progressively dominate with few unweathered corestones, and the biancana gradually collapses into a penultimate “soufflé-like” form.  相似文献   

2.
We report the results of an analytical investigation into the deformation behaviour of rate-dependent granular material as a refinement of previous studies on seepage phenomena during shear. The rheology has two components—a compliant part of the constitutive law associated with grain contacts as deformation takes place (dilatancy), and a rate-dependent viscous force transmitted by the melt phase. This formulation allows intermediate, time-dependent behaviour to be assessed for the dilatant porous medium. A key result is that during shear, the magnitude of the excess pore pressure first decreases then increases back to its initial value. Two characteristic timescales are identified that control the rate-dependent dilatancy of the mixture, τ1, the time constant that rules the increase of the magnitude of the excess pore pressure, and τ0 that controls its decline. We consider the dilatant effect to be an internal constraint in deforming magmas in the lithosphere and other porous (partially molten) regions in the solid earth. When such regions are exposed to external loading, secular pressure changes should drive fluid flow independent of local buoyancy forces, for the duration of the governing rate-dependent timescales. The accumulated heave of the process is also estimated.  相似文献   

3.
Numerical simulation of the propagation of P waves in fractured media   总被引:1,自引:0,他引:1  
We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q −1 of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q −1 of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) −1 in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.  相似文献   

4.
McKenzie's model of sedimentary basin evolution and its modification, widely used in geophysics, sometimes fails to explain discrepancies between predicted and observed values of extension, thinning and subsidence of the Earth's crust, as for the North Sea. We develop a numerical model of sedimentary basin evolution based on the mechanism suggested by Lobkovsky. In the course of rifting, accompanied by thinning of lower parts of the lithosphere, the roof of the underlying asthenosphere moves upward. the material of the mantle lifts and partially melts owing to the reduction of pressure. the density difference between the melt and the crystalline skeleton results in the filtration of the lighter melt and its accumulation in the form of a magmatic lens. Due to changed P-T conditions, the material of the lens undergoes the gabbro-eclogite phase transformation. the resultant anomalously heavy eclogite lens sinks in the surrounding material. This induces a viscous flow, changing the surface topography and forming a sedimentary basin. We construct a 2-D numerical model describing a viscous flow induced by subsidence of a heavy body and compute changes of surface topography. to compute the flow we employ the Galerkin-spline approach, with modifications allowing for density discontinuities and time dependence of the phase transformation. We apply the model to the cases of the Illinois, Michigan and Williston basins. the computed and tectonic subsidence curves agree well for these cases. the proposed model is compatible with the seismic structure of the crust and upper mantle below these basins. the model is also consistent with gravity data. the approach is applicable to other intracratonic basins.  相似文献   

5.
A theoretical approach to the propagation of interacting cracks   总被引:2,自引:0,他引:2  
We propose a scheme to compute interaction effects between two randomly oriented cracks under compressive stresses and we discuss the role crack interactions play in the crack coalescence process. Stress intensity factors are computed by using an iterative technique based on the method of successive approximations. Once crack propagation occurs, curved wing cracks grow from the initial crack tips. The stress intensity factors at the wing crack tips are calculated as the sum of two terms: a component for a single wing crack subjected to both the applied stresses and the interaction effect, and a component due to the sliding of the initial crack. We have applied our procedure to various crack geometries. Our results show that interaction effects act on the crack propagation path. For cracks under tension, our approach correctly predicts the curving, hook-shaped paths of interacting cracks that have been observed in various materials. For en echelon compressive cracks, interaction effects depend on the geometry of stepping. For right-stepping cracks, no mode I crack coalescence occurs. A mixedmode propagation criterion may be introduced to check whether coalescing secondary shear fractures initiate. For left-stepping cracks, depending on whether or not there is overlapping, crack coalescence is achieved by tension wing cracks at the inner crack tips. Without overlapping, the growing wing cracks delimit a region where a tensile secondary fracture may develop and lead to coalescence. These results are consistent with previous work and show that our procedure may be now extended to a population of cracks.  相似文献   

6.
Laboratory measurements of ultrasonic wave propagation in tuffaceous sandstone (Kimachi, Japan) and granite (Iidate, Japan) were performed during increasing fracturing of the samples. The fracturing was achieved by unconfined uniaxial compression up to and beyond the point of macrofracture of the specimen using a constant low strain rate. The observed variation of wave velocity (up to 40 per cent) due to the development of micro- and macrofractures in the rock is interpreted by rock models relating velocity changes to damage and crack density. The calculated density of the newly formed cracks reaches higher values for the sandstone than for the granite. Using the estimated crack densities, the attenuation behaviour is interpreted in terms of different attenuation mechanisms; that is, friction and scattering. Rayleigh scattering as described by the model of Hudson (1981 ) may explain the attenuation qualitatively if the largest plausible crack dimensions are assumed in modelling.  相似文献   

7.
陶波  李锋  马威  刘建雄  易守勇 《热带地理》2022,42(10):1761-1770
采用工程地质钻探、物探、地质测绘及室内试验等技术方法探讨飞鹅山Ⅲ号滑坡形成机理与防治技术。结果表明:1)滑坡体主要岩性为泥质粉砂岩,飞鹅山滑坡属于新形成的深层中型牵引式滑坡,在平面上呈圈椅状。2)滑坡属于双层滑面滑坡,主滑面以中型深层滑坡为主,主滑体上部发育中型中厚层滑坡。3)滑坡产生的原因为:①泥质粉砂岩倾向与坡向基本一致,且岩层倾角为中等倾角;②人工开挖使坡脚形成高陡临空面,抗滑力大为降低;③雨水沿层面及节理裂隙入渗至坡体深部,大大增加岩土体容重,同时泥质粉砂岩遇水软化,抗剪强度显著降低。4)结合该滑坡区地质环境条件,采用坡面削坡+锚杆(索)+格构梁+双排预应力锚拉抗滑桩+三维网植草绿化+截排水+毛石挡墙的综合治理方法进行防治,监测结果显示该滑坡变形及位移已得到有效控制,整治效果良好。  相似文献   

8.
The weathering of granodiorite porphyry is examined in the extremely cold and dry environment of the inland part of Antarctica. Weathering features include granular disintegration, rock varnish, and sheeting. Sheeting has gradually proceeded since the exposing of the nunatak, but the other types of weathering have not actively advanced in recent times.
Granular disintegration primarily comprises the release of individual crystals. The mechanism of this release is that cracks are created along crystal boundaries and cleavages in phenocrysts mainly due to differential thermal expansion. Plagioclase and hypersthene are released more readily than quartz. Disintegration of plagioclase has produced many prism-shaped holes, 0.5 to 2.0 mm in length, while fine-grained quartz crystals, 0.05 to 0.1 mm in diameter, are incidentally released without cracking.
The reddish brown (10R4/4) rock varnishes result from oxidation and consist of limonite, which fills cracks and penetrates into crystals. Ferrous iron in hypersthene, biotite, ilmenite and magnetite is transformed into limonite by oxidation. Manganese is not found in the varnishes. Sulfur, which is important for oxidation and which may have originated from adhered snow, is concentrated on the surface of the rock.
Sheeting has precipitated rock falls, and has gradually formed a gray mosaic on the varnished wall. The sheeting was caused by gravitational body force, which is internal stress of the rock body due to the mass of the overburden.  相似文献   

9.
The effect of cracks on the elastic properties of an isotropic elastic solid is studied when the cracks are saturated with a soft fluid. A polynomial equation in effective Poisson's ratio is obtained, whose coefficients are functions of Poisson's ratio of the uncracked solid, crack density and saturating fluid parameter. Elastic and dynamical constants used in Blot's theory of wave propagation in poroelastic solids are modified for the introduction of cracks. The effects of cracks on the velocities of three types of waves are observed numerically. The frequency equation is derived for the propagation of Rayleigh-type surface waves in a saturated poroelastic half-space lying under a uniform layer of liquid. Dispersion curves for a particular model of oceanic crust containing cracks are plotted. The effects of variations in crack density and saturation on the phase and group velocity are also analysed.  相似文献   

10.
This paper describes a method for determining Moho depth, lithosphere thinning factor (γ= 1 − 1/β) and the location of the ocean–continent transition at rifted continental margins using 3-D gravity inversion which includes a correction for the large negative lithosphere thermal gravity anomaly within continental margin lithosphere. The lateral density changes caused by the elevated geotherm in thinned continental margin and adjacent ocean basin lithosphere produce a significant lithosphere thermal gravity anomaly which may be in excess of −100 mGal, and for which a correction must be made in order to determine Moho depth accurately from gravity inversion. We describe a method of iteratively calculating the lithosphere thermal gravity anomaly using a lithosphere thermal model to give the present-day temperature field from which we calculate the lithosphere thermal density and gravity anomalies. For continental margin lithosphere, the lithosphere thermal perturbation is calculated from the lithosphere thinning factor (γ= 1 − 1/β) obtained from crustal thinning determined by gravity inversion and breakup age for thermal re-equilibration time. For oceanic lithosphere, the lithosphere thermal model used to predict the lithosphere thermal gravity anomaly may be conditioned using ocean isochrons from plate reconstruction models to provide the age and location of oceanic lithosphere. A correction is made for crustal melt addition due to decompression melting during continental breakup and seafloor spreading. We investigate the sensitivity of the lithosphere thermal gravity anomaly and the predicted Moho depth from gravity inversion at continental rifted margins to the methods used to calculate and condition the lithosphere thermal model using both synthetic models and examples from the North Atlantic.  相似文献   

11.
Summary . Expressions are available in the literature for the propagation constants of a 'mean' wave travelling in a material with a random distribution of cracks. These are approximations for small crack density and crack dimensions small compared with a wavelength. The formulae provided to second order in crack density by the method of smoothing are here extended to the case where the cracks consist of two or more sets, aligned in different directions.  相似文献   

12.
Summary. The polarizations of shear waves recorded by networks of digital three-component seismometers immediately above small earthquakes near the North Anatolian Fault in Turkey display shear-wave splitting on almost all shear-wave seismograms recorded within the shear-wave window. This splitting is incompatible with source radiation-patterns propagating through simple isotropic structures but is compatible with effective anisotropy of the internal structure of the rock along the ray paths. This paper interprets the phenomena in terms of widespread crack-induced anisotropy. Distributions of stress-induced cracks model many features of the observations, and synthetic polarization diagrams calculated for propagation through simulated cracked rock are similar to the observed patterns. This evidence for widespread crack-induced anisotropy lends strong support to the hypothesis of extensive-dilatancy anisotropy (EDA) suggested by laboratory experiments in subcritical crack-growth. The crucial evidence confirming some form of EDA would be observations of temporal changes in shear-wave splitting as the stress field alters the crack density and crack geometry. There is some weak evidence for such temporal changes at one site, but further analysis of suitable digital three-component seismometer networks in seismic areas is required to confirm EDA.  相似文献   

13.
In the southwest U.S., cracks in alluvial fan surface clasts have a preferred orientation independent of rock fabric and shape. In this paper, we show that differential insolation of incipient cracks of random orientations predicts a distribution of crack orientations consistent with field observations. In this model, crack growth by hydration and/or thermal weathering is primarily a function of local water content at the crack tip. Crack tips that experience minimal solar insolation maintain a greater average moisture and, hence, weather more rapidly than cracks that experience greater solar insolation. To show this, we used a numerical radiative transfer code to quantify the solar insolation of rectangular cracks at 35° N. latitude with a range of depths and orientations. The amount of solar energy reaching the bottom of each crack was calculated at 5-min intervals over the day for several days of the year to determine hourly, daily, seasonal, and annual energy deposition as a function of crack depth and orientation. By assuming that only crack orientations that effectively shield their interiors and minimize their water loss are able to grow, the pattern of cracks produced by the model is consistent with field observations. The annual average insolation, which controls water retention, is associated with the two primary modes of crack orientation. The effect of daily recharge by summer rains of the North American monsoon system is consistent with the observed deviations from these primary modes. Model results suggest that both the annual average insolation and the daily pattern of rainfall is recorded in the preferred crack orientations of surface clasts in the southwest U.S.  相似文献   

14.
Wave speeds and attenuation of elastic waves in material containing cracks   总被引:38,自引:0,他引:38  
Summary. Expressions now exist from which may be calculated the propagation constants of elastic waves travelling through material containing a distribution of cracks. The cracks are randomly distributed in position and may be randomly orientated. The wavelengths involved are assumed to be large compared with the size of the cracks and with their separation distances so that the formulae, based on the mean taken over a statistical ensemble, may reasonably be used to predict the properties of a single sample. The results are valid only for small concentrations of cracks.
Explicit expressions, correct to lowest order in the ratio of the crack size to a wavelength, are derived here for the overall elastic parameters and the overall wave speeds and attenuation of elastic waves in cracked materials where the mean crack is circular, and the cracks are either aligned or randomly orientated. The cracks may be empty or filled with solid or fluid material. These results are achieved on the basis of simply the static solution for an ellipsoidal inclusion under stress.
The extension to different distributions of orientation or to mixtures of different types of crack is quite straightforward.  相似文献   

15.
Summary. It is known that flow in the mantle can produce preferred orientation in olivine crystals with seismic anisotropy as a consequence. Flow in the subcrustal lithosphere is unlikely because of the high viscosity. Lenses of high temperature and low-viscosity ( anomalous mantle ) are located under the crust in many tectonically active regions, and viscous flow can easily arise in such material resulting in seismic anisotropy. After cooling, such anomalous mantle acquires high viscosity and becomes incorporated into the lithospheric layer preserving the anisotropy produced by the flows which existed previously. The interaction of the stresses with cracks in the upper crust can be one of the causes of anisotropy in this layer.  相似文献   

16.
Effects of fractures on seismic-wave velocity and attenuation   总被引:1,自引:0,他引:1  
The effects of fractures on the seismic velocity and attenuation of a rock are investigated using theoretical results and experimental data. Fractures in a rock mass influence the traveltimes and amplitudes of seismic waves that have propagated through them. The displacement discontinuity model, recently employed in fracture investigations, is modified to describe the effect of fractures on seismic-wave velocity and attenuation. This new model, the modified displacement discontinuity model (MDD), is formulated in a way analogous to transmission-line analysis. The fractures are treated as transmission lines for the passage of seismic waves. The MDD takes into consideration realistic fracture parameters which include the fracture length, the fractional area of a fracture surface in contact, and the nature of the infilling material. A single fracture of varying geometric and material properties is shown to affect dramatically the transmission properties of a propagating waveform, and hence the seismic velocity and attenuation. These effects have been shown to result in a frequency-dependent velocity and attenuation. The sensitivity of the fracture parameters to seismic-wave velocity and attenuation was investigated and interesting results were obtained. Fracture parameters used in designing experimental models consisting of synthetically manufactured cracks were fed into the MDD and a well-known crack model, Hudson's model, for comparison. Velocities as a function of the incident-wave angle were obtained from both numerical models and were compared with the results from the experimental modelling. For P waves, the MDD model results show better agreement with those of the experimental model for all crack densities investigated than those from Hudson's model.  相似文献   

17.
Summary . A two-space singular perturbation technique is employed to derive approximate governing equations for flow of a viscous heat-conducting fluid through a rigid porous solid. It is assumed that buoyancy forces are significant, and it is shown that standard approximations used in the study of flow through a porous medium are valid provided that Gr ≫ 1, where Gr is a Grashof number calculated using a typical pore radius as the length scale. Results previously derived in the literature for flow through an isotropic random array of spherical particles are used to show how the permeability and conductivity tensors can be calculated for a problem of interest in planetary science.  相似文献   

18.
青藏高原山地系统动力学数值模拟研究   总被引:1,自引:2,他引:1  
毕思文 《山地学报》2001,19(4):289-298
青藏高原山地系统地壳巨厚,约为正常地壳厚度的两倍,由于第四纪以来高原的急剧隆升而成为全球构造研究的热点,根据山地系统岩石圈运动的主要因素-构造带和断裂带,本文将青藏高原山地系统的地壳演变化简化成一个平面问题和一个典型剖面问题的力学模型,应用变形体模拟构造带,采用摩擦模型模拟断裂带,并根据周边位移情况给出边界条件,应用ABAQUS有限元分析软件和位移加载,变成了青藏高原山地系统岩石圈应力-位移场吻合较好,由此分析和推断了青藏高原山地系统演变的历史和现状。  相似文献   

19.
The McKenzie model proposed in 1978, which is widely used in calculating the thermal history of rift basins and other extensional basins, incorrectly assumes that all heat passing through the lithosphere originates below the lithosphere. In reality, heat from radiogenic sources within the lithosphere, especially in the upper crust, may represent more than half the heat flow at the top of basement. Thinning of the lithosphere during extension does indeed result in an increase of heat flowing from the asthenosphere, but this thinning also reduces the radiogenic heat from within the lithosphere. Because these two effects cancel to a large degree, the direct effects of lithospheric extension on heat flow at the top of basement are smaller than those predicted by the McKenzie model. Because of permanent loss of radiogenic material by lithospheric thinning, the heat flow at the top of basement long after rifting will be lower than the pre-rift heat flow.The McKenzie model predicts an instantaneous increase in heat flow during rifting. The Morgan model proposed in 1983, however, predicts a substantial time delay in the arrival of the higher heat flow from the asthenosphere at the top of basement or within sediments. Using the Morgan model, heat flow during the early stages of rifting will actually be lower than prior to rifting, because the time delay in the loss of radiogenic heat is less than the time delay in arrival of new heat from the asthenosphere.  相似文献   

20.
A model is suggested to simulate the physical aspect of diagenesis in porous rocks. A bidisperse ballistic deposition model with relaxation of deposited grains is used to generate the porous structure. Sedimentation and erosion are allowed to restructure the pore space as a fluid flows through the rock. The effect of this restructuring of the pore space on permeability is studied. The Navier–Stokes equation is solved numerically by the finite difference method to determine the pressure and velocity distributions in the pore space. We find that though deposition is the dominant process in our model of diagenesis, reducing the porosity, the permeability may increase dramatically in some cases. These are when the erosion takes place at a single narrow constriction in the pore channel.  相似文献   

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

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