Changing grading of soil: effect on critical states   总被引：6，自引：3，他引：3  

David Muir Wood  Kenichi Maeda 《Acta Geotechnica》2008,3(1):3-14

Examples of situations are presented where the grading of a soil changes during its lifetime either by crushing of particles leading to an increase of fine material or by slow transport of fine particles with seepage leading to a decrease of fine material. Such grading changes influence the basic constitutive properties of the soil, in particular properties such as critical states which are dependent on the available range of densities of packing. Discrete element modelling is used to show the dependence of critical state conditions on grading and the way in which the particle assembly seeks out new critical state conditions as the grading changes.  相似文献   

Towards a nonperturbation transport theory in hetereogeneous aquifers     

Sergio E. Serrano 《Mathematical Geology》1996,28(6):701-721

A large body of existing theories of flow and contaminant transport in aquifers ignore the presence of recharge, eliminate the boundary conditions, neglect transient conditions in groundwater flow, conceive hydraulic gradients as linear, and require parameter variability to be stationary and Gaussian. The most outstanding and difficult to justify assumption is the subjective small size of the stochastic terms (i.e., small perturbation methods), which usually is forced by considering the logarithm of the hydraulic conductivity. Several problems in flow and contaminant subsurface hydrology, such as the enhanced dispersion parameters with plume size or time after injection, remain to be observed in the light of a stochastic theory that allows a more realistic consideration of physical and hydrologic properties. In this article, an attempt is made to reformulate a contaminant transport equation (the variable dispersion equation, VDE) with transport parameters in terms of regional hydrologic and aquifer hydraulic properties, such as recharge rate, spatially random transmissivity, hydraulic gradient, aquifer thickness, and soil porosity. Subsequently, a general analytic procedure, the method of decomposition, is used to derive a solution to the VDE. This procedure does not require small perturbation, logarithmic transformations, or specific probability law assumptions. Comparison tests with existing theoretical and field results are given. The tests illustrate the enhanced dispersion and shifting concentration effects produced by the variable dispersion equation. Finally a generalization of the method to nonstationary dispersion in three-dimensional domains is proposed.  相似文献   

黄河下游断面形态与水沙输移关系及数学模拟方法   总被引：4，自引：0，他引：4  

梁志勇  曾庆华 《地理研究》1993,12(2):37-46

以前人对黄河下游的实测资料分析为基础,本文讨论断面形态与来水来沙关系、断面形态对输水输沙的影响,并根据实测资料,提出黄河下游弯段与直段断面冲淤变化的两种模式,建立了主流摆动及坍岸影响断面形态变化以及纵向冲淤量在断面的分布两种数学模拟方法。  相似文献   

基于GIS城市洪水淹没模拟分析   总被引：10，自引：0，他引：10       下载免费PDF全文

向素玉  陈军 《地球科学》1995,20(5):575-580

城市洪水淹没范围动态模拟分析是城市防洪规划与防治决策的一个重要基础工作。本文主要研究洪水从翻堤口出发在地理空间蔓延、扩散、动态行进及确定淹没范围的数字模拟方法。为此,根据数学形态学及测地圆概念,研究设计了洪水扩散范围的“膨胀”模拟算法和淹没范围搜索算法,用于查询淹没通块中从翻堤处到任一点之间的淹没路径和ｔ时刻洪水扩散范围。  相似文献   

Using numerical modelling to evaluate tsunami hazard near the Kuril Island     

L. B. Chubarov  Yu. I. Shokin  K. V. Simonov 《Natural Hazards》1992,5(3):293-318

A sequence of computer experiments is used to study questions concerning the tsunami problem as a quantitative estimate of tsunami danger, detailed geographical tsunami classification, determination of the parameters of critical tsunami waves, and the conditions of their development. We call a wave critical, if its impact on the coast is most hazardous.Using the Middle Kuril Island as an example, we present the results of a computer experiment which includes determining the wavefields on the shelf and estimating the effects connected with the deep-water Bussol and Diana Straits.Numerical simulation of tsunami waves of different sources permits the assessment of the extent of tsunami danger in different areas of the coastal zone of Simushir Island, depending on the location of the focus zone and their geometry.The major singularities of the wavefield arise in the zones of the deep-water straits. The distribution of the amplification factors is determined by both the global parameters of the wavefields and the local properties of individual harbours. The results obtained for a particular harbour in the northern part of Simushir Island, formed the basis for the quantitative estimate of tsunami danger for this area.  相似文献   

Numerical study of droplet size dependent chemistry in oceanic,wintertime stratus cloud at southern mid-latitudes     

G. P. Ayers  T. V. Larson 《Journal of Atmospheric Chemistry》1990,11(1-2):143-167

Cloud droplet chemistry is modelled for the first 150 m of rise in a wintertime, mid-latitude, marine stratus cloud using observations made at and near the Cape Grim Baseline Station as a source of input parameters. The emphasis in this work was to study the variation in droplet chemistry as a function of both droplet size and nucleus composition, with a particular focus on the way in which oxidation of dissolved sulfur dioxide varied.At 150 m above the condensation level, solute concentration as a function of droplet size was found to increase by as much as 2 to 3 orders of magnitude for only a factor of 2 increase in droplet radius, primarily as a consequence of the 1/r dependence in the droplet growth equation. This type of size dependence exists at all levels in the model cloud, and has a significant influence on oxidation rate of sulfur dioxide in droplets growing on sulfate nuclei, oxidation by ozone being favoured in the smallest droplets, but oxidation by hydrogen peroxide being favoured in the larger droplets. Oxidation by ozone is favoured at all sizes in droplets formed on sea-salt nuclei as a result of the initially high alkalinity of these droplets, and in the cloud overall is calculated to be the more important oxidation pathway. Although based on a simplified chemical scheme, these results suggest that both size-dependent and nucleus-dependent chemistry of cloud droplets may need to be considered explicitly in cloud modelling work.Volume-weighted mean pH values in the range 5 to 6 were predicted from sensitivity studies in which input variables were varied over reasonable ranges, in agreement with two sets of bulk cloud-water pH data obtained by aircraft near Cape Grim.  相似文献   

Controls on low-pressure anatexis   总被引：2，自引：0，他引：2  

C. C. GERBI  S. E. JOHNSON  P. O. KOONS 《Journal of Metamorphic Geology》2006,24(2):107-118

Low-pressure anatexis, whereby rocks melt in place after passing through the andalusite stability field, develops under more restricted conditions than does low-pressure metamorphism. Our thermal modelling and review of published work indicate that the following mechanisms, operating alone, may induce anatexis in typical pelitic rocks without inducing wholesale melting in the lower crust: (i) magmatic advection by pervasive flow; (ii) crustal-scale detachment faulting; and (iii) the presence of a high heat-producing layer. Of these, only magmatic advection by pervasive flow and crustal-scale detachment faulting have been shown quantitatively to provide sufficient heat to cause widespread melting. Combinations of the above mechanisms with pluton-scale magmatic advection, shear heating, removal of the lithospheric mantle, or with each other provide additional means of developing suitable high temperatures at shallow crustal levels to generate low-pressure anatexis.  相似文献   

Three-dimensional model and late stage warping of the Plattengneis Shear Zone in the Eastern Alps     

Martin Putz  Kurt Stüwe  Mark Jessell  Philippe Calcagno   《Tectonophysics》2006,412(1-2):87-103

The Plattengneis shear zone is a 250–600 m thick, flat lying, Cretaceous, eclogite facies, mylonitic shear zone, with north-over-south transport direction, that is exposed over almost 1000 km² in the Koralpe region along the eastern margin of the Alps. Although the shear zone is one of the largest in the Alps, its role in the Eoalpine metamorphic evolution and the subsequent exhumation of the region, remain enigmatic and its large-scale geometry is not well understood. The outcrop pattern suggests that the shear zone is made up of a single sheet that is folded into a series of open syn- and antiforms with wavelengths of about 10 km. Eclogite bodies occur above, within and below the shear zone and there is no metamorphic grade change across the shear zone. In the south, the fold axes strike east–west and plunge shallowly to the east. In the north, the fold axes are oriented in north–south direction and form a dome shaped structure of the shear zone. Total shortening during this late stage warping event was of the order of 5%. Indirect evidence constrains this folding event to have occurred between 80 and 50 Ma and the fold geometry implies that the final exhumation in the Koralpe occurred somewhat later than further north. Interestingly, the shear zone appears to strike out of the topography in the south and dip into the topography in the north, so that north of the shear zone only hanging-wall rocks are exposed and south of it only foot-wall rocks. Possibilities for the geometric relationship of the Plattengneis shear zone with the surrounding south dipping detachments are discussed.  相似文献   

An investigation of upper mantle heterogeneity beneath the Archaean and Proterozoic crust of western Canada from Lithoprobe controlled-source seismic experiments     

Andrew R. Gorman  Balzs Nmeth  Ron M. Clowes  Zoltan Hajnal 《Tectonophysics》2006,416(1-4):187

Observations of upper mantle reflectivity at numerous locations around the world have been linked to the presence of a heterogeneous distribution of rock types within a broad layer of the upper mantle. This phenomenon is observed in wide-angle reflection data from Lithoprobe's Alberta Basement Transect [the SAREX and Deep Probe experiments of 1995] and Trans-Hudson Orogen Transect [the THoRE experiment of 1993]. SAREX and Deep Probe image the Archaean lithosphere of the Hearne and Wyoming Provinces, whereas THoRE images the Archaean and Proterozoic lithosphere of the Trans-Hudson Orogen and neighbouring areas.Finite-difference synthetic seismograms are used to constrain the position and physical properties of the reflective layer. SAREX/Deep Probe modelling uses a 2-D visco-elastic finite-difference routine; THoRE modelling uses a pseudospectral algorithm. In both cases, the upper mantle is parameterized in terms of two media. One medium is the background matrix; the other is statistically distributed within the first as a series of elliptical bodies. Such a scheme is suitable for modelling: (1) variations in lithology (e.g., a peridotite matrix with eclogite lenses) or (2) variations in rheology (e.g., lenses of increased strain within a less strained background).The synthetic seismograms show that the properties of heterogeneities in the upper mantle do not change significantly between the two Lithoprobe transects. Beneath the Trans-Hudson Orogen in Saskatchewan, the layer is best modelled to lie at depths between 80 and 150 km. Based on observations from perpendicular profiles, anisotropy of the heterogeneities is inferred. Beneath the Precambrian domains of Alberta, 400 km to the west, upper mantle heterogeneities are modelled to occur between depths of 90 and 140 km. In both cases the heterogeneous bodies within the model have cross-sectional lengths of tens of kilometers, vertical thicknesses less than 1 km, and velocity contrasts from the background of − 0.3 to − 0.4 km/s. Based on consistency with complementary data and other results, the heterogeneous layer is inferred to be part of the continental lithosphere and may have formed through lateral flow or deformation within the upper mantle.  相似文献   

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model     

lvaro Gonzlez  Miguel Vzquez-Prada  Javier B. Gmez  Amalio F. Pacheco 《Tectonophysics》2006,424(3-4):319

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.  相似文献