Traction image method for irregular free surface boundaries in finite difference seismic wave simulation   总被引：16，自引：1，他引：16  

Wei Zhang  Xiaofei Chen 《Geophysical Journal International》2006,167(1):337-353

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.  相似文献   

Three types of pyroclastic currents and their deposits: examples from the Vulsini Volcanoes, Italy     

Danilo M. Palladino  Silvia Simei 《Journal of Volcanology and Geothermal Research》2002,116(1-2)

This paper deals with ground-hugging, gas–pyroclast currents from explosive volcanic eruptions and their deposits. Key field observations and laboratory determinations are proposed to relate specific deposit types with flow regimes and particle concentration in the transport and depositional systems. Three relevant flow scenarios and corresponding deposit types have been recognized from a survey of pyroclastic successions of the Vulsini Volcanic District (central Italy): (1) dilute, turbulent, pyroclastic currents producing normally or multiply graded beds by direct suspension sedimentation; (2) concentrated bedload regions beneath suspension currents, depositing inversely graded beds by traction carpet sedimentation; (3) self-sustained, high particle concentration, laminar, mass flows developing massive, poorly sorted bodies, with opposite grading of coarse lithic and pumice clasts, overlying fine-grained, inversely graded, basal layers. Main distinguishing criteria include the occurrence and pattern of clast grading, clast–thickness relationships, grain size, ash matrix componentry and pyroclast size–density relationships. Downcurrent and temporal transitions among identified flow scenarios are likely to occur for changing energy conditions and gas–pyroclast ratio both on regional and local scales. The nature and efficiency of magma fragmentation, volatile content, conduit geometry (which determine the characteristics of the erupted mixture and possible lateral blast component at the vent), and the angle of incidence of the column collapse, are suggested as the main factors controlling the generation of one type over the other at flow inception. Dilute, fine-grained, overpressured eruption clouds are thought to favor the formation of low particle concentration turbulent currents. Column collapse over slightly inclined volcano slopes, causing a high degree of compression of the collapsing mixture and of gas expulsion, would favor the generation of high particle concentration pyroclastic currents.  相似文献   

甘孜州色达县典型牵引式滑坡地质灾害形成机理分析--以甘孜州甲学乡甲热滑坡为例          下载免费PDF全文

罗雲丰  魏良帅 《探矿工程》2018,45(8):93-97

甲热滑坡位于甘孜州色达县甲学乡境内,属于典型的牵引式滑坡,滑坡体平均坡度约35°,主要由碎石土组成,滑坡长约195 m,宽约280 m,平均厚约15 m,滑坡体总方量约60万m3,属中型滑坡。滑坡坡体中上部变形现象较明显,若发生滑动,将直接威胁当地41户160人居民的生命财产安全。本文在进行野外现场基础地质调查后,对滑坡体的变形破坏机制进行了定性分析,并利用三位有限元数值模拟进行了定量分析,对以后研究同类滑坡的稳定性具有重要意义。  相似文献   

Depth–transport functions and erosion–deposition diagrams as indicators of slope inclination and time-averaged traction forces: applications in tropical reef environments     

Johann Hohenegger  & Elza Yordanova 《Sedimentology》2001,48(5):1025-1046

The trigonometric relationship between slope inclination, the horizontally acting time‐averaged traction force and the vertical depth of transport allows the estimation of one factor, when both others are known. Depth–transport functions can be deduced by comparing the depth distributions of living organisms and their skeletal remains, and this paper simplifies this comparison using foraminifera in which a single test represents an individual. Differences in distribution parameters between living individuals and empty tests allow depth–transport functions to be determined; these functions differ between species at a single transect according to the varying buoyancies of the tests. Within a single species, differences in depth–transport functions between locations are based on either slope inclination or traction intensities. After establishing a mean depth–transport function by averaging species‐characteristic functions, the time‐averaged traction force acting on the studied transect can be calculated. Transport intensities are also estimated using an erosion–deposition diagram that combines the relative frequency distributions of living individuals and empty tests. The proportion of ‘eroded’, ‘parautochthonous’ and ‘allochthonous’ tests mirrors the influence of both slope inclination and traction force for the deposition of empty tests. To test the model, six species of symbiont‐bearing benthic foraminifers were investigated at two transects in front of a NW Pacific coral reef. One transect is distinguished by a strong slope flattening below the steep reef slope (30 m), whereas further steepening characterizes the equivalent part in the other transect. These differences are mirrored in the depth–transport functions as well as in the erosion–deposition diagrams of all species. The time‐averaged traction forces differ in intensities between transects, because of the position of the reef front with respect to the predominant wind direction. However, the form of the functions is identical and distinguished by an increase from the surface to 35 m depth, followed by a decrease down to 105 m. This can be explained by successive onshore and offshore forces acting on the shallow slope, such as the tropical cyclones that cross the region every summer.  相似文献   

天津轨道交通对地磁观测干扰的影响研究   总被引：5，自引：0，他引：5       下载免费PDF全文

谢凡  滕云田  徐学恭  邢西淳  胡星星 《地震学报》2011,33(2):252-261

直流牵引城市轨道交通系统在运行过程中产生的磁场对我国大城市周边地磁台站的观测数据产生较为严重的干扰.城市轨道交通运行过程中产生泄露电流是产生干扰的重要原因.本文以真空中的线电流磁场计算模型为基础,通过简化的远场磁场干扰模型对天津地铁轨道交通干扰进行理论建模.该模型计算简洁高效.通过对模型的数值模拟,获得了干扰随距离衰减...  相似文献   

Exact solution of earth-flattening transformation for P-SV waves: Taking surface wave as an example     

黄慧  陈晓非 《地震学报(英文版)》2008,21(2):109-117

Taking surface wave as an example this paper proposes an exact solution of earth-flattening transformation for P-SV waves and discusses the applicability of the approximate methods. The results show that the transform parameter m has little influence on the final results, and on the condition of short wave approximation, approximate earth-flattening transformation is suitable. Moreover, the efficiency of approximate transformation is twice of that of exact transformation. For low frequency problems exact earth-flattening transformation should be used.  相似文献   

水牛蹄型深海集矿机履齿设计及牵引性能研究     

丁聪  陈旭光  刘学麟  刘茜茜  李宗晏 《海洋工程》2024,(2):26-41

作为深海采矿系统的重要组成部分,深海集矿机的技术研究和开发一直是许多国家的难题。为了提高深海集矿机在稀软沉积物上的牵引性能,基于水牛蹄的曲面结构特征和特殊几何结构,设计了一款仿水牛蹄蹄型履齿。以该款履齿的结构参数为研究对象,通过对朗肯被动土压力理论的优化,得到了仿水牛蹄蹄型履齿的牵引力解析解,建立了对应的履齿仿生参数对履齿牵引力的影响公式。结合正交试验方法进行了不同形式履齿在沉积物上运动的单、多履齿剪切试验,将仿水牛蹄蹄型履齿与直型履齿、仿水牛蹄轮廓履齿（另一种仿生履齿）进行对比试验,并通过模型车试验验证了该款履齿的可行性。结果表明：3种履齿中,仿水牛蹄蹄型履齿提供的最大牵引力最高。研究可为进一步优化仿水牛蹄蹄型履齿结构参数和提高深海稀软底质机械的牵引性能提供参考依据。  相似文献   

浅议沉积学中的流体问题     

李向东 《世界地质》2020,39(1):45-55

笔者从流体的角度,对近年来有关流体沉积问题的研究成果进行整理、总结、归纳和诠释。从沉积学的角度总体上可将流体分为牵引流、过渡流和重力流3大类:其中牵引流是研究的基础,可进一步分为单向流、双向流和振荡流3种基本类型以及叠置流、复合流和叠加流3种复合类型,可按照流体力学的相关方法展开研究;过渡流和重力流沉积研究的关键在于各种重力流类型之间(浊流、液化沉积物流、颗粒流和碎屑流)及其和牵引流之间的相互转化。流体的沉积机制可分为垂直降落沉积(静水)、底床阻碍沉积(牵引流)和能量减弱沉积(重力流)3大类,各种流体的交互作用最终可形成丘状(洼状)交错层理的水动力条件。通过综述可以发现,从流体角度进行沉积学研究尽管困难,但仍然可行,特别是对于复杂水动力条件下的沉积学研究具有一定的意义。  相似文献   

大连地区南北向剪切滑褶变形   总被引：4，自引：0，他引：4  

王宗秀  杨天南  唐哲民 《岩石矿物学杂志》2001,20(2):208-216

大连地区经历了漫长而又特殊的地质构造演化过程，自中生代以来，全区处于活动大陆边缘环境，并遭受多期构造变形事件的影响，其中大规模南北向剪切滑褶变形为资产发现并确定，其为区内众多变形事件中规模最大、影响最广的一期构造事件，本文在详细的地质调查基础上，初步研究分析了研究区南北向剪切滑褶变形的构造-岩石组合特征及变形模式，并就变形动力学背景及变形时代进行了探讨。  相似文献   

基于控制区间牵引算法的地下施工变形预测     

马恩临  赖金星  王立新  汪珂  雷升祥  李储军  邱军领 《岩土力学》2023,(2):577-594

变形是地下结构施工安全状态的重要评价指标。针对数据驱动变形预测方法中地下施工动态影响因素难以量化这一问题,提出了基于控制区间牵引算法的地下施工变形预测模型。以双向长短时记忆网络（Bi-LSTM）对变形监测时间序列进行预测,以数值模拟在关键施工阶段处的结果为牵引点。根据实测值随阶段更新牵引点,并基于注意力机制双向长短时记忆网络（Bi-LSTM-AM）以牵引点修正数据驱动模型在控制区间的预测结果,实现更准确、更智能的地下施工变形预测。设置牵引相对权重,使模型可自适应判断当前合理牵引程度,准确融合了Bi-LSTM与数值模拟的结果。通过相关历史案例与数据,验证了牵引预测模型的有效性。依托西安市西咸新区丰镐三路地下通道上跨地铁1号线自动化监测实例,采用牵引算法预测了地下结构变形。结果表明：在关键施工阶段处,牵引作用改善了数据驱动预测方法存在的滞后问题,各阶段误差平均降低了24.34%;不断优化的牵引点逐渐趋近真值,降低了数值模拟出现偏差时造成的影响。该方法可为地下施工变形预测问题提供新途径。  相似文献