共查询到18条相似文献，搜索用时 171 毫秒
1.
五大连池火山群老黑山火烧山熔岩流流动特征及其形成动力过程 总被引：2，自引：0，他引：2
从火山灾害的预测预防角度出发，叙述了五大连池火山群老黑山、火烧山近代火山喷发所形成的熔岩流产物——熔岩的分布，表面特征，流动过程。讨论了漫流、溪谷熔岩流和潜流三种流动方式。对结壳熔岩、渣状熔岩、块状熔岩、喷气锥、熔岩遂道等的形成动力过程以及对古地形古环境条件的制约做了深入的讨论，提出新的形成模式和一些新的认识。 相似文献
2.
天池火山东北侧造盾玄武岩可划分出8个流动单元,熔岩流的流动距离主要集中在30~50km,熔岩流宽度以5km左右为主。通过由野外调查获得的天池火山东北侧不同熔岩流单元的地表坡度、熔岩流厚度等,结合温度、密度与黏度等物理参数,按照熔岩流速度公式恢复的头道组和早白山组0.5m厚晶体含量5%的玄武岩熔岩流流速集中在0~1m/s之间。晶体含量为30%、厚度为0.5m的晚白山组和老房子小山组玄武岩熔岩流的流动速度集中在0~0.12m/s之间。厚度增大至2m左右,晶体含量不变的头道组和早白山组的玄武岩熔岩流流动速度可加快至11m/s。天池火山2m厚的碱性熔岩流在12h内达到或接近了它的最远距离,而各组内2m厚拉斑玄武岩熔岩流在20h内接近了最远距离。0.5m厚的熔岩流在10d内接近最大距离。50km是预计的熔岩流长度,在未来制定减灾措施时,可将此长度作为重要依据之一。天池火山熔岩流灾害主要表现为熔岩流动时对房屋建筑、农田、道路、林地、电站的毁坏,火灾及大量的人口伤亡 相似文献
3.
湖泊模型为数值天气预报模型提供热量通量、水汽通量和动量通量等下边界条件,但是不同时间尺度上湖泊水热通量变化的控制因子不同,因此有必要对湖泊模型进行多时间尺度上的离线评估.本文利用2012-2016年太湖中尺度通量网避风港站的气象资料和辐射数据驱动CLM4-LISSS模型(Community Land Model version 4-Lake,Ice,Snow and Sediment Simulator),并与涡度相关观测(Eddy Covariance,EC)结果进行对比,以年平均潜热通量模拟结果最佳为目标调整了模式中的消光系数、粗糙度长度方案,研究了该模型从半小时到年尺度上对湖表温度和水热通量的模拟性能.结果表明:模型对湖表温度的模拟在各时间尺度上均比较理想,但是模拟的日较差较小;从半小时到年尺度上潜热通量的变化趋势都能被很好地模拟出来,但在季节尺度上,潜热通量的模拟出现了秋冬季偏高、春夏季偏低的情况,季节变化模拟不准确.湖表温度和潜热通量模拟偏差的原因可能是消光系数的参数化方案.相比之下,感热通量尽管年际变化趋势的模拟值与观测值一致,但是从半小时到年尺度均被高估.特别地,冷锋过境期间,模型能较好地模拟出潜热通量和感热通量的变化趋势,但对于高风速条件下的感热通量模拟效果不佳.本文的研究结果能为湖泊模式的应用与发展提供有用信息. 相似文献
4.
在山区进行MT勘探时,用规则网格有限元方法模拟起伏地形会受到限制.本文采用非结构三角网格可以有效地模拟任意二维地质结构,如起伏地形、倾斜岩层和多尺度构造等.正演引入自适应有限元方法,其在网格剖分过程中能根据单元误差自动细化网格,保证了正演结果的精度.将自适应有限元与Occam算法结合,且引用并行处理技术提高正反演计算速度.通过对比两个理论模型,讨论了地形对MT正演响应的影响;其次进行了不同地电模型带地形反演展示了本文算法的正确性和适用性;最后将该方法应用于实测MT数据处理,证明了自适应非结构有限元方法是复杂地形下处理MT数据的有力工具. 相似文献
5.
将大陆岩石层视为由幂指数律控制的一层薄层,它上伏在粘滞性较低的软流层之 上蠕变流动,其运动限制在与东亚大陆构造形态较相似的边界模型的梯形框架之中．设印度 板块以一恒定的速度向北推进,并被视为青藏高原挤压隆升的主要动力．用数值模拟的方法 研究了青藏高原的挤压隆升演化过程,并对数值模拟的隆升过程作了剥蚀修正．结果表明, 由挤压模型所产生的地形和现代青藏高原及其邻区的地形格局比较吻合．同时也表明,挤压 隆升过程受多种因素（如岩石层的力学特性、边界条件以及剥蚀作用）的制约,无论从空间还 是从时间上看,模拟所反映的高原隆升都是不均匀的演化过程． 相似文献
6.
本文利用优化的25点频率-空间域有限差分算法对基于BISQ模型双相各向同性介质中的地震波进行了数值模拟.通过与经典的Biot模型理论模拟结果进行对比,分析了Biot流动（宏观流体流动）和Squirt流动（微观流体流动）耦合作用对地震波在孔隙介质中传播特性的影响.数值模拟在地震频段进行,结果显示：在理想相界和黏滞相界情况下,Squirt流动机制都比Biot流动机制产生了更大的速度频散和能量衰减.其中,在Biot流动和Squirt流动耦合作用下的快P波的速度和振幅小于仅考虑Biot流动影响下快P波速度和振幅,而且慢P波的衰减也更加强烈.本文还研究了地震波在双层双相各向同性介质分界面处的反射和透射特征,双相介质中波的反射与透射现象类似于单相介质的情况.模拟结果表明,利用优化25点频率-空间域有限差分法模拟双相孔隙介质中的地震波场是可行的,这为开展双相孔隙介质全波形反演问题的研究提供了可能. 相似文献
7.
2008年5月12日四川汶川地区发生MW7.9地震，震中位置103.4°E，31.06°N.这次地震造成了以汶川、映秀为中心及其周边地域建筑物的严重破坏和人员的重大伤亡，且因为高山等地形复杂区域抢险救灾的艰巨性，为及时救援造成很大干扰.为更好理解地形因素对于强地面数值模拟结果的影响，建立了包含地形起伏影响及去除地形影响的两类模型.同时，依据震源破裂过程运动学反演结果，建立了包含障碍体破裂过程的震源滑动模型，实现断层分段、空间倾角以及滑移角的动态设定.基于动力学的地震动模拟方法，通过对地震波传播过程的数值计算和后处理分析，模拟由地震激发的区域强地面运动过程.结果显示：（1）强震动台站的断层距对地形效应具有放大或抑制作用，距离断层破裂带越近，地形效应越明显，反之，距离越远，则地形效应越微弱；（2）因为地形高差与障碍体的影响，地震造成的峰值可能出现在震中区域之外；（3）考虑地形影响模型的地表峰值速度（PGV）区域位于汶川与北川附近；而未考虑地形影响模型的PGV区域位于灌县—江油断层的后半段，处清平、安县附近；对汶川地震近实时强地面运动波场的模拟、峰值图谱的圈定及未来大地震强地面运动特征的预测都有重要指示意义. 相似文献
8.
9.
本文基于自适应非结构有限元算法实现海洋电磁起伏海底地形三维正演模拟.通过采用隐式后推欧拉时间离散技术,保证在较大的时间步长条件下获得正确结果.为获得多时间道海洋电磁正演模拟的有效网格,我们采用基于法向电流连续的后验误差估计的自适应方法和网格融合技术;同时为了控制网格数量和保证正演模拟稳定性,我们还在网格融合过程中应用了随机网格挑选技术.对于方程组求解我们使用MUMPUS直接求解器.当时间步长不变时,只需对系数矩阵进行一次分解,大大提高计算速度.将本文计算结果与半空间模型解析解进行对比,验证了本文算法精度.针对海洋电磁半拖曳式和双船拖曳式工作方式,我们通过典型模型的模拟计算,研究海底地形影响及海底高阻层识别特征. 相似文献
10.
地下热水资源定量评价是地热资源开发利用中的重要环节.在地下热水系统中,地下水流场和温度场相互影响,因而在研究地下热水运移时须同时考虑地下热水的流动和热量的运移.本文建立三维不稳定地下水流-热运移耦合模型,模拟陕西渭北东部岩溶地下热水运移特征.针对研究区多组断层纵横交错而且对岩溶水运移起着控制作用的特点,模型中将断层剖分为条带以更好地刻画其水力学特征.利用水化学测试结果估算地下水-地表水转化量,为模型边界流量的确定提供依据.结合地下水测年数据进行模型校正与参数识别,提高了模型仿真性. 相似文献
11.
We present a kinematic, self-adaptive, numerical model to describe the down-flow thermal and rheological evolution of channel-contained lava. As our control volume of lava advances down a channel it cools and crystallizes, an increasingly thick and extensive surface crust grows, and its heat budget and rheology evolve. By estimating down-flow heat and velocity loss, our model calculates the point at which the control volume becomes stationary, giving the maximum distance lava flowing in the channel can extend. Modeled effusion rates, velocities, widths, surface crust parameters, heat budget, cooling, temperature, crystallinity, viscosity, and yield strength all compare well with field data collected during eruptions at Mauna Loa, KOlauea, and Etna. Modeled lengths of 25-27, 2.5-5.7, and 0.59-0.83 km compare with measured lengths of 25-27, 4, and 0.75 km for the three flows, respectively. Over proximal flow portions we calculate cooling, crystallization, viscosity, and yield strength of 1-10°C km-1, 0.001-0.01 volume fraction km-1, 103-104 Pa s, and 10-3-102 Pa, respectively. At the flow front, cooling, crystallization, viscosity, and yield strength increase to >100°C km-1, 0.1 volume fraction km-1, 106-107 Pa s, and 103-104 Pa, respectively, all of which combine to cause the lava to stop flowing. Our model presents a means of (a) analyzing lava flow thermo-rheological relationships; (b) identifying important factors in determining how far a channel-fed flow can extend; (c) assessing lava flow hazard; and (d) reconstructing flow regimes at prehistoric, unobserved, or remote flows. 相似文献
12.
A dynamical model of lava flows cooling by radiation 总被引：1，自引：0，他引：1
Michele Dragoni 《Bulletin of Volcanology》1989,51(2):88-95
The behaviour of a lava flow is reproduced by a two-dimensional model of a Bingham liquid flowing down a uniform slope. Such a liquid is described by two rheological parameters, yield stress and viscosity, both of which are strongly temperature-dependent. Assuming a flow rate and an initial temperature of the liquid at the eruption vent, the temperature decrease due to heat radiation and the consequent change in the rheological parameters are computed along the flow. Both full thermal mixing and thermal unmixing are considered. The equations of motion are solved analytically in the approximation of a slow downslope change of the flow parameters. Flow height and velocity are obtained as functions of the distance from the eruption vent; the time required for a liquid element to reach a certain distance from the vent is also computed. The gross features of observed lava flows are reproduced by the model which allows us to estimate the sensitivity of flow dynamics to changes in the initial conditions, ground slope and rheological parameters. A pronounced increase in the rate of height increase and velocity decrease is found when the flow enters the Bingham regime. The results confirm the observation according to which lava flows show an initial rapid advance, followed by a marked deceleration, while the final length of a flow is such that the Graetz number is in the order of a few hundreds. 相似文献
13.
We introduce a 3D model for near-vent channelized lava flows. We assume the lava to be an isothermal Newtonian liquid flowing
in a rectangular channel down a constant slope. The flow velocity is calculated with an analytical steady-state solution of
the Navier-Stokes equation. The surface velocity and the flow rate are calculated as functions of the flow thickness for different
flow widths, and the results are compared with those of a 2D model. For typical Etna lava flow parameters, the influence of
levees on the flow dynamics is significant when the flow width is less than 25 m. The model predicts the volume flow rate
corresponding to the surface velocity, taking into account that both depend on flow thickness. The effusion rate is a critical
parameter to evaluate lava flow hazard. We propose a model to calculate the effusion rate given the lava flow width, the topograhic
slope, the lava density, the surface flow velocity, and either the lava viscosity or the flow thickness.
Received: 20 January 1998 / Accepted: 8 January 1999 相似文献
14.
The evolution of lava flows emplaced on Mount Etna (Italy) in September 2004 is examined in detail through the analysis of morphometric measurements of flow units. The growth of the main channelized flow is consistent with a layering of lava blankets, which maintains the initial geometry of the channel (although levees are widened and raised), and is here explicitly related to the repeated overflow of lava pulses. A simple analytical model is introduced describing the evolution of the flow level in a channelized flow unit fed by a fluctuating supply. The model, named FLOWPULSE, shows that a fluctuation in the velocity of lava extrusion at the vent triggers the formation of pulses, which become increasingly high the farther they are from the vent, and are invariably destined to overflow within a given distance. The FLOWPULSE simulations are in accordance with the observed morphology, characterized by a very flat initial profile followed by a massive increase in flow unit cross-section area between 600 and 700 m downflow. The modeled emplacement dynamics provides also an explanation for the observed substantial “loss” of the original flowing mass with increasing distance from the vent. 相似文献
15.
The dynamics of a lava flow is studied by a two-dimensional model describing a viscous fluid with Bingham rheology, flowing down a slope. The temperature in the flow is calculated assuming that heat is transferred through the plug by conduction and is lost by radiation to the atmosphere at the top of the flow. Taken into account is that the increasing crystallization takes place in the flow as a consequence of cooling. The lava viscosity and yield stress are expressed as a function of crystallization degree as well as of temperature: in particular it is assumed that yield stress reaches a maximum value above the solidus temperature, according to experimental data. Dynamical variables, such as velocity and thickness of the flow, are calculated for different values of the maximum crystallization degree and the flow rate. The model shows how the lava flow dynamics is affected by cooling and crystallization. The cooling of the flow is controlled by the increase of yield stress, which produces a thicker plug and makes the heat loss slower. The increasing crystallization has two opposing effects on viscosity: it produces an increase of viscosity, but at the same time produces an increase of yield stress and hence reduces the heat loss and keeps the internal temperature high. As a consequence, lava flows are significantly affected by the dependence of yield stress on temperature and scarcely by the maximum crystallization degree. 相似文献
16.
Volcanic geomorphologists have investigated various relationships between eruption rate and morphologic parameters of lava flows, particularly with regard to preferred statistical correlations and the conditions under which they are valid. Here we employ two simple models for lava flow heat loss by Stefan-Boltzmann radiation to derive eruption rate versus planimetric area relationships. Both of these models predict a linear relationship between eruption rate and planimetric area, modulated by distinct prefactors potentially sensitive to compositional and temperature differences among different flows. Regardless of any theoretical considerations, we show that eruption rate is highly correlated with planimetric area for the Hawaiian basaltic flows analyzed in this work. Moreover, this observed correlation is superior to those from other obvious combinations of eruption rate and flow dimensions. On the basis of the theoretical models for lava flow heat loss, the correlations obtained here suggest that the surfaces of Hawaiian flows radiate at an effective temperature much less than the inner parts of the flowing lava in agreement with numerous field observations. This work also indicates that eruption rate versus planimetric area correlations can be markedly degraded when data from different vents, volcanoes and epochs are combined. These previously unrecognized sensitivities identified by the thermal loss modeling may have contributed to past unresolved debates on relationships between eruption rates and morphologic dimensions. 相似文献
17.
由给定的随时间和空间变化的二维应变速率分布计算速度场,根据速度场计算二维热传导方程,得到岩石圈的热演化结构,然后计算由热扰动引起的密度变化,进而确定加载到岩石圈上的载荷,根据载荷计算挠曲方程,最终得到盆地的沉降.在此基础上讨论了热源对地表热流和沉降的影响以及岩石圈有效弹性厚度对盆地几何形态的影响.结果表明,热源对地表热流影响显著,对沉降影响不大,岩石圈有效弹性厚度则控制盆地的几何形态.二维应变速率正演为下一步的反演奠定了基础. 相似文献
18.
The initial cooling of pahoehoe flow lobes 总被引：1，自引：0，他引：1
In this paper we describe a new thermal model for the initial cooling of pahoehoe lava flows. The accurate modeling of this
initial cooling is important for understanding the formation of the distinctive surface textures on pahoehoe lava flows as
well as being the first step in modeling such key pahoehoe emplacement processes as lava flow inflation and lava tube formation.
This model is constructed from the physical phenomena observed to control the initial cooling of pahoehoe flows and is not
an empirical fit to field data. We find that the only significant processes are (a) heat loss by thermal radiation, (b) heat
loss by atmospheric convection, (c) heat transport within the flow by conduction with temperature and porosity-dependent thermal
properties, and (d) the release of latent heat during crystallization. The numerical model is better able to reproduce field
measurements made in Hawai'i between 1989 and 1993 than other published thermal models. By adjusting one parameter at a time,
the effect of each of the input parameters on the cooling rate was determined. We show that: (a) the surfaces of porous flows
cool more quickly than the surfaces of dense flows, (b) the surface cooling is very sensitive to the efficiency of atmospheric
convective cooling, and (c) changes in the glass forming tendency of the lava may have observable petrographic and thermal
signatures. These model results provide a quantitative explanation for the recently observed relationship between the surface
cooling rate of pahoehoe lobes and the porosity of those lobes (Jones 1992, 1993). The predicted sensitivity of cooling to
atmospheric convection suggests a simple field experiment for verification, and the model provides a tool to begin studies
of the dynamic crystallization of real lavas. Future versions of the model can also be made applicable to extraterrestrial,
submarine, silicic, and pyroclastic flows.
Received: 26 November 1994 / Accepted: 1 December 1995 相似文献