Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

陕西关中盆地中部地下热水H、O同位素交换 及其影响因素

对关中盆地地下热水δ＆18O和8D数据的研究表明：盆地中部西安、咸阳深部的地压地热流体发生明显的δ^18O同位素交换,并出现^2H同位素交换,表明热储流体发生了强烈的水岩反应,盆地周边及中部的非地压地热流体^18O交换则不明显。根据研究区^18O同位素的交换程度（用^2H过量参数d表征）和水化学资料,可将关中盆地热储流体分为循环型和封闭型热储流体2类。地热水埋深越大、滞留时间越长、TDS和温度越高、地质环境越封闭,18^O交换程度就越大。西安和咸阳地下热水分属于不同的地热系统,具有不同的补给来源。     

双作用水泵泡沫增压装置应用及前景分析

双作用水泵泡沫增压装置是利用双缸双作用水泵的高压能力将低压空气和泡沫液二次增压至水泵的额定压力,并混合产生泡沫.采用BWZ-1200双作用水泵泡沫增压装置在北京地区深度为1 000 m以下的地热井钻探中,成功地克服了地热井钻进中井内漏失和洗井的难题;在极度缺水的宁夏南部地区水井钻探中,采用BWZ-1100双作用水泵泡沫增压装置解决了水井孔内严重漏失和施工用水的难题.该装置可代替价格昂贵的高压空压机实现气举反循环钻探,应用于深地热井、深水井和油气井钻探.     

鄂尔多斯盆地黄陵、东胜地区地温场对比

鄂尔多斯盆地黄陵、东胜铀矿区分别处于盆地南部渭北隆起的北侧边缘和盆地北部伊盟隆起的东部,赋矿层位都是中侏罗统直罗组。盆地南、北铀矿区在现今地温场及古地温场都存在明显差异,南部现今大地热值及热演化程度明显高于北部。对于下侏罗统延安组和石炭—二叠系煤层,黄陵地区镜质体反射率都高于东胜地区。通过镜质体反射率资料得出同一埋深的一套地层经历的最大古地温和对应的古地温梯度也有南部高于北部的现象。由于早白垩世后期盆地普遍整体抬升使得现今地温相对古地温降低,南部黄陵地区抬升剥蚀量大于北部东胜地区,导致古、今地温差异也大于后者。盆地南部庆阳—富县一带局部构造热运动,导致南部异常地温场的形成,使得南部热演化程度高于北部。     

鄂尔多斯高原地下水流系统的多层结构循环模式—来自深孔中PACKER系统分层水头测定的证据

地下水流系统和循环模式分析是研究地下水形成机理的基础,对正确认识和评价地下水资源具有重要意义。鄂尔多斯高原在含水系统和众多地下水排泄区的控制下,形成了多个不同规模、不同循环深度、相互独立的地下水流系统。PACKER系统分层试验测定的不同深度水头的数据证明,鄂尔多斯高原地下水流系统存在托斯多层水流模式,区域性水流系统一般包含浅循环、中间循环和深循环3个循环系统。浅循环系统的发育深度在200m以内,深循环系统的发育深度大于400m。     

提高管线探测水平分辨率的技术方法

在平顶山矿工路路面改造工程中,使用GX-2型及RD-4000型地下管线探测仪探测地下管网布设情况。根据地下管线种类、管线布设方位、埋置深度、探测难度大小进行了多次试验,在此基础上确定了电磁探测技术与方法：利用直连法提高被测管线中的交变电流,压制邻近平行管线和地下介质中的异常反映;在多管并存、且间距较小的情况下,应选择“梯度法”对磁场水平分量垂直梯度ΔHx进行观测,以得到最大清晰异常;另外,还应根据具体情况,合理选择诸如压制旁侧管线法、选择发射法、偏移感应法、动源发射法等发射方式,保证目标管线中有较强异常呈现。在管线密集区,应尽可能地降低工作频率,以减小旁侧管线中产生的二次电流及二次磁场;为防止信噪比下降,可适当减小收发距,以提高接收机的灵敏度。     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

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

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

 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?

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.