Micromechanical approach to unsaturated water flow in structured geomaterials by two-scale computations     

Adam Szymkiewicz  Jolanta Lewandowska 《Acta Geotechnica》2008,3(1):37-47

This article presents a micromechanical approach to the problem of unsaturated water flow in heterogeneous porous media in transient conditions. The numerical formulation is based on the two-scale model obtained previously by periodic homogenization. It allows for a coupled solution of the non-linear flow equations at macroscopic and microscopic scales and takes into account the macroscopic anisotropy of the medium and the local non-equilibrium of the capillary pressure. The model was applied to simulate two-dimensional water infiltration at constant flux into an initially dry medium containing inclusions of square and rectangular shapes. The numerical results showed the influence of the inclusion–matrix conductivity ratio and the local geometry on the macroscopic behavior. The influence of the conductivity ratio manifested itself by the acceleration or retardation of the onset of the macroscopic water flux at the outlet, while the local geometry (anisotropy) significantly affected the macroscopic spatial distribution of the water flux. Such type of approach can be extended to simulate coupled phenomena (for example hydro-mechanical problems) with evolving local geometry.  相似文献   

The dynamics and thermodynamics of large ash flows   总被引：6，自引：6，他引：0  

Marcus I. Bursik  Andrew W. Woods 《Bulletin of Volcanology》1996,58(2-3):175-193

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996  相似文献   

Accelerating ground warming in the Canadian Prairie provinces: Is it a result of global warming?     

Jacek A. Majorowicz 《Pure and Applied Geophysics》1996,147(1):1-24

The results of precision temperature logs made to depths of several hundred meters in some 80 wells in Western Canada, most of which are located in the Prairie Provinces, show evidence of warming at the ground surface in the 0.5 K to 3.5 K range (average=2.2±0.7 K, for 80 unevenly distributed sites). Modeling shows that this warming mostly pertains to this century and it has been most substantal in the last four decades if the ramp function of the linear increase of surface temperature is assumed. Using the step function model's increase of surface temperature (land clearing, forest fires, etc.) the calculated onset of warming would pertain mostly to the last two decades. Contour maps of ground temperatures currently and previously and a contour map of the ground warming magnitude dilineate a large regional character of the ground temperature change at the southern marigin of permafrost for the large area of the Prairie Provinces. In many cases however, the magnitude of ground warming is much larger than the magnitude of air warming. This is especially evident for the northern areas of Alberta in the boreal forest ecoprovince. The magnitude of ground warming is equal to the magnitude of surface air warming in southern Alberta in the grassland and aspen parkland ecoprovinces. The analysis of the temperature depth response to the surface warming from well data shows the integrated effect of surface air warming together with the increases in ground temperature due to natural terrain effects and other anthropogenical changes to the surface of the earth.  相似文献   

地下水非稳定流的LTFAM算法     

王文科  李俊亭 《地球科学与环境学报》1993,(4)

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。  相似文献   

南海海盆区莫霍面分布规律及其对深部钻探的意义     

秦绪文  张宝金  赵斌  路允乾  陈玺  王利杰  许振强  张如伟  耿明会  杨振  李建平  吕文超  尉建功 《地质学报》2022,96(8):2635-2646

钻遇莫霍面是人类一直以来的梦想。深海海底是地球上离莫霍面最近的地方，目前有研究推测南海是世界上莫霍面深度最浅的海域之一，但缺乏足够的直接证据。深反射地震探测可以直接揭示岩石圈的构造形态，是莫霍面探测的重要手段。本文基于长达15000 km的深反射多道地震剖面的解释、处理、制图和分析，结合前人的研究，形成了南海海盆区莫霍面反射特征和空间分布的初步认识。① 南海东部次海盆南部早期经历了较快速扩张，岩浆供应充足，受扩张停止后岩浆活动影响较小，基底平坦，地质构造相对简单，同时洋壳地震速度结构不存在异常，且有较强的广角莫霍面反射波和可识别的地幔顶部折射波，具备莫霍面钻探的基本条件。② 南海海盆不同区域的莫霍面反射强度存在较大差异。其中东部次海盆莫霍面反射最为强烈且清晰，西北次海盆次之，西南次海盆仅有零星出现的清晰莫霍面反射且可信度不高。③ 识别南海海盆区莫霍面地震反射长度超过3500 km，首次形成了海盆区深度域莫霍面地震反射空间分布图。与重力反演的莫霍面深度相比，利用深反射多道地震计算的莫霍面深度细节更为丰富，并且可以在垂向上清晰刻画莫霍面的结构。整体上，南海海盆区莫霍面地震反射强烈和可信度高的区域中，深度较浅的区域之一是东部次海盆南部，最浅处仅约9. 5 km，其中水深4. 01 km，洋壳厚度仅5. 54 km。综合判断，东部次海盆南部是南海重要的莫霍面钻探备选区，这对南海莫霍面钻探选址具有重要意义。  相似文献   

地铁施工对济南白泉泉群流量的影响     

束龙仓  王小博  李虎  倪寒茜  李罡  余亚飞  王鑫  张曼琦 《吉林大学学报(地球科学版)》2021,51(1):192-200

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。  相似文献   

含煤层地质环境下地震波场的数值模拟   总被引：3，自引：2，他引：1  

潘冬明  杨顶辉 《地球物理学进展》1997,12(1):74-83

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。  相似文献   

用气象站地温资料计算多年平均土壤热流的初步结果   总被引：11，自引：3，他引：8  

董文杰  汤懋苍 《高原气象》1992,11(2):115-125

  相似文献   

Scattering of elastic waves by a fracture zone containing randomly distributed cracks     

Jun Kawahara  Teruo Yamashita 《Pure and Applied Geophysics》1992,139(1):121-144

We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ ^–1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ ^–1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ ^–1 is proportional tok ^–1 independently of model parameters and the wave modes. If the seismological observation thatQ ^–1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.  相似文献   

Transmission and reflection of Rayleigh waves at a thin low velocity vertical layer: A laboratory and theoretical study     

P. F. Biagi  T. S. Gegechkory  E. N. Its  P. Manjgaladze  V. Sgrigna  D. Zilpimiani 《Pure and Applied Geophysics》1990,133(2):317-327

A three-dimensional ultrasonic modelling was used to study the Rayleigh waves reflected from a thin low velocity vertical layer and those transmitted through the same zone. The reflection and transmission coefficients of the Rayleigh waves have been calculated by an approximate method based on the Green function technique. Nonrigid contact boundary conditions were used for the low velocity layer modelling. A rather good correspondence between the theoretical and experimental absolute values of the reflection and the transmission coefficients was obtained. Some discrepancies between experimental and theoretical results, in particular with regard to phase shifts, may be explained by a certain inadequacy of the theoretical model to the low velocity layer investigated by the ultrasonic modelling.  相似文献