Numerical simulation of a stratigraphic model     

Véronique Gervais  Roland Masson 《Computational Geosciences》2008,12(2):163-179

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.  相似文献   

Strain localization due to structural in-homogeneities in the Central European Basin System     

Mauro Cacace  Ulf Bayer  Anna Maria Marotta 《International Journal of Earth Sciences》2008,97(5):899-913

The large-scale crustal deformations observed in the Central European Basin System (CEBS) are the result of the interplay between several controlling factors, among which lateral rheological heterogeneities play a key role. We present a finite-element integral thin sheet model of stress and strain distribution within the CEBS. Unlike many previous models, this study is based on thermo-mechanical data to quantify the impact of lateral contrasts on the tectonic deformation. Elasto-plastic material behaviour is used for both the mantle and the crust, and the effects of the sedimentary fill are also investigated. The consistency of model results is ensured through comparisons with observed data. The results resemble the present-day dynamics and kinematics when: (1) a weak granite-like lower crust below the Elbe Fault System is modelled in contrast to a stronger lower crust in the area extending north of the Elbe Line throughout the Baltic region; and (2) a transition domain in the upper mantle is considered between the shallow mantle of the Variscan domain and the deep mantle beneath the East European Craton (EEC), extending from the Elbe Line in the south till the Tornquist Zone. The strain localizations observed along these structural contrasts strongly enhance the dominant role played by large structural domains in stiffening the propagation of tectonic deformation and in controlling the basin formation and the evolution in the CEBS.  相似文献   

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula   总被引：7，自引：3，他引：4  

Vicente Medina  Marcel Hürlimann  Allen Bateman 《Landslides》2008,5(1):127-142

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.  相似文献   

分布式光纤传感技术在滑坡监测中的应用   总被引：5，自引：1，他引：4  

史彦新  张青  孟宪玮 《吉林大学学报(地球科学版)》2008,38(5)

将分布式光纤传感技术引入滑坡监测,既可得到整个滑坡体的概要特征,又能提高监测效率。FBG与BOTDR是两种最具代表性的分布式光纤传感技术,FBG通过测量其反射光波长的变化获得应变或温度值。BOTDR通过测定后向布里渊散射光的频移实现分布式温度、应变测量。FBG传感器灵敏度高,但只能实现离散点的准分布式测量,BOTDR可实现分布式、长距离、不间断测量,但其空间分辨率不高。笔者提出FBG与BOTDR联合监测滑坡的方案。在巫山残联滑坡,在整个滑坡体上铺设监测光纤,利用BOTDR获得整个滑坡体的概要信息;在滑坡体变形的关键部位———变形缝安装FBG传感器,获得某些关键部位的应变信息,从而实现由点到线再到面的监测,获得滑坡体较完整的应变信息。  相似文献   

努力提高三维地震勘探精度，为西部煤炭工业做出新贡献     

孙升林 《中国煤炭地质》2008,(6):1-5

以国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”项目的由来、意义和总体研究目标为引,概括的介绍了项目依托工程中各个专项技术研究完成情况,并对非均匀介质成像技术、高精度三维地震静校正技术、高密度采集技术、特观技术、岩性反演技术、属性体解释技术等六项重大关键技术取得的突破性进展进行了重点说明。指出随着我国煤炭生产重点的逐步西移,应加强诸如叠前、叠后深度偏移技术的研究,以解决复杂山区三维地震面元内地震反射波散射问题,提高其三维地震勘探精度,为西部煤炭工业做出新贡献!  相似文献   

应用全站仪进行三角高程测量方法的探讨     

吴昊  蒋建新 《现代测绘》2008,31(4)

传统的三角高程水准测量由于受地形起伏的限制,外业工作量大,施测精度和工作效率不高。本文采用全站仪对传统的三角高程水准测量方法进行了改进,并论述了全站仪三角高程测量方法的原理和实际作业方法。  相似文献   

三维地震小面元采集技术在晋城矿区的应用     

张胤彬 《中国煤炭地质》2008,(6):70-72

依托“西部煤炭资源高精度三维地震勘探技术”工程,对晋城矿区进行了旨在提高小断层,小陷落柱探测能力的高密度三维地震勘探。根据面元选择因素及该区地质任务,采用5m×5m网格进行野外数据采集;考虑炮检距、方位角、覆盖次数、排列片横纵比及煤层埋深（350～500m）等因素,采用中点放炮、60道接收,24次覆盖（横向4次,纵向6次）的8线16炮束状观测系统,基岩中激发。原始资料经同一处理流程后,获得5m×5m×1ms、5m×10m×1ms、10m×10m×1ms及2.5m×2.5m×1ms不同单元的三维数据体多个,通过对比可以发现小断层,小陷落柱在其小面元叠加时间剖面、顺层切片及相干切片都有清晰的反映。实例说明,小面元采集技术可以提高对小构造的纵、横向分辨能力,满足山区对三维地震精确勘探的要求。  相似文献   

Strain analysis in Jurassic argillites of the Monte Sirino area (Lagonegro Zone,southern Apennines,Italy) and implications for deformation paths in pelitic rocks     

S. Mazzoli 《International Journal of Earth Sciences》1995,84(4):781-793

Analysis of strain in Jurassic argillites forming part of the folded and thrusted sedimentary succession of the Lagonegro basin (southern Italian Apennines) has been carried out using ellipsoid-shaped reduction spots as strain markers. Most of the determined finite strain ellipsoids are of oblate type and show a peculiar distribution of the maximum extension direction (X), with maxima either subparallel or subperpendicular to the local fold axes. Using the strain matrix method, two different deformation histories have been considered to assist the interpretation of the observed finite strain pattern. A first deformation history involved vertical compaction followed by horizontal shortening (occurring by a combination of true tectonic strain and volume loss), whereby all strain is coaxial and there is no change in the intermediate axis of the strain ellipsoid. By this type of deformation sequence, which produces a deformation path where total strain moves from the oblate to the prolate strain field and back to the oblate field, prolate strain ellipsoids can be generated and may be recorded where tectonic deformation has not been large enough to reverse pretectonic compaction. This type of deformation history may be of local importance within the study area (i.e. it may characterize some fold hinge regions) and, more generally, is probably of limited occurrence in deformed pelitic rocks. A second deformation sequence considered the superposition of pre-tectonic compaction and tectonic strain consisting of initial layer-parallel shortening followed by layer-parallel shear (related to flexural folding). Also in this instance, volume change during tectonic deformation and tectonic plane strain have been assumed. For geologically reasonable amounts of volume loss due to compaction and of initial layer-parallel shortening, this type of deformation history is capable of producing a deformation path entirely lying within the oblate strain field, but still characterized by a changeover, during deformation, of the maximum extension axis (X) from a position parallel to the fold axis to one perpendicular to it. This type of deformation sequence may explain the main strain features observed in the study area, where most of the measured finite strain ellipsoids, determined from the limb regions of flexural folds, display an oblate shape, irrespective of the orientation of their maximum extension direction (X) with respect to the local structural trends. More generally, this type of deformation history provides a mechanism to account for the predominance of oblate strains in deformed pelitic rocks.  相似文献   

Uniaxial compression and tension tests of anthracite and loading rate dependence of peak strength   总被引：2，自引：0，他引：2  

S. Okubo  K. Fukui  Qi Qingxin 《International Journal of Coal Geology》2006,68(3-4):196-204

The mechanical characteristics of coal under uniaxial compressive stress and uniaxial tensile stress are key factors in investigations of the stability of galleries and coal faces, and are vital considerations for the efficient design of coalmines, disaster prevention, and environmental preservation. These mechanical characteristics have long been a topic of research, but remain insufficiently understood. In this study we performed uniaxial compression tests and uniaxial tension tests on coal, with particular attention to two concerns. The first was to measure the loading rate dependence of the peak strength of coal. Coal shows high sample-to-sample scatter in strength; test methods that are successful in comparatively homogeneous rock types cannot be used for coal. We therefore used an alternative test method that was recently developed by the authors. Anthracite samples were subjected to alternating slow and fast strain rates. Measured variations in stress during this process were used to estimate the loading rate dependence of peak strength. The second objective was to obtain complete stress–strain curves for coal under uniaxial tensile stress. It is difficult to hold samples secure during a tensile test; consequently, such curves have yet to be obtained. This study presents a successful application of the authors' method to the analysis of coal samples, yielding a complete stress–strain curve under tensile stress. The two methods presented here hold promise for application not only to anthracite but also to a wide variety of coals. It is possible to derive the values of the constants used in constitutive equations from the obtained experimental results. Once these equations are determined, they can be incorporated into finite-element software to investigate various time-dependent behaviors of coal and aid in the efficient design of coalmines and the prevention of subsidence and mine disasters.  相似文献   

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D     

Peter Frolkovič  Jozef Kačur 《Computational Geosciences》2006,10(3):279-290

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).  相似文献