Numerical 3D modeling of the effects of strain rate and confining pressure on the compressive behavior of Kuru granite

This paper deals with numerical modeling of the compressive behavior of granite rock under high strain rate dynamic loading and wide range of confining pressure. For this end, a constitutive model based on damage mechanics and viscoplasticity for rock is formulated and implemented in explicit dynamics FEM. The viscoplastic part is based on a simple power law type yield criterion that incorporates the rate-dependency with a linear viscosity term. Moreover, a Rankine type of tensile cut-off is employed. The damage part of the model is formulated with separate scalar damage variables in tension and compression. The model is calibrated for Kuru granite and validated with the experimental data from dynamic compression tests at the strain rate of 600 1/s up to 225 MPa of confining pressure. The numerical simulations demonstrate that, despite the underlying continuum modeling approach, the model captures the correct experimental failure modes, including the transition from single-to-multiple fragmentation, as well as the dynamic compressive strengths at different confining pressures.     

Cone penetration test in sand: A numerical-analytical approach

Separation of the effects of initial horizontal stress and relative density on cone tip resistance in sandy soils has been a complicated issue for many years. In order to overcome this problem, a numerical modeling of CPT which has been verified by calibration chamber tests, has been used in this paper to achieve a reliable analytical solution. The analytical solution has resulted in two relationships for sleeve friction and cone tip resistance in terms of the initial conditions of sandy soil. Based on the presented solution, the initial horizontal stress and relative density can be determined according to CPT measurements.     

塔里木盆地玉北地区埋藏史及热史分析

根据地层分层数据和古热流分析,对塔里木盆地玉北地区的5口代表性钻井进行了埋藏史、热史的模拟与分析,探讨了研究区烃源岩的热演化成熟过程,初步划分了烃源岩的成熟阶段和生烃期次。研究结果认为：玉北地区寒武系－泥盆系沉积时期,具有多次沉积埋藏和隆升剥蚀的交替;泥盆系沉积之后,石炭－二叠系呈现沉积埋藏的特征;中生代在玉北地区东部少量沉积三叠系,随后被剥蚀,而在中西部发生沉积间断;进入新生代,地层持续沉积,且沉积速率呈现增大趋势。研究区寒武系烃源岩自沉积以来成熟度演化相对连续,奥陶纪时进入生烃门限,成熟时期较早;志留－泥盆纪,受海西早期运动的影响,发生隆升剥蚀作用,烃源岩演化变缓;石炭－二叠纪受控于区域上持续稳定的地层沉积与广泛发育的火山活动,研究区寒武系烃源岩进入中高成熟阶段;由于喜马拉雅造山运动的远程效应,新生代研究区持续沉降,形成巨厚沉积,烃源岩成熟度达到高－过成熟,为天然气生成阶段     

基于谱元法的各向异性黏弹性介质地震波场模拟和分析

基于谱元法计算理论,通过坐标旋转算法实现任意各向异性黏弹性介质中的地震波场计算,研究复合介质中地震波场特征。根据计算结果,qSV波表现为剪切滞弹性形变,qP波表现为膨胀滞弹性形变;品质因子不仅影响地震波振幅,还影响其传播频率;黏弹性介质吸收高频快,在一定频率范围内,黏弹性介质没有频散现象。     

交通地理信息系统中的时空动态分段模型研究

动态分段技术是交通地理信息系统(GIS-T)中一项重要的线性要素动态显示与分析技术,该文提出了一种时空动态分段模型,使动态分段系统的组成成分融入时态因素,并将属性的时态信息作为属性信息表的字段存储,通过时态地理信息系统技术完成物理段的时态变化。利用统一建模语言(UML)的类图设计方式开发了物理模型,描述了时空动态分段模型中对象之间的关系,并定义了必要的属性和操作;同时基于线性参照系统(LRS)表达空间实体的方法,概括了模型中所涉及到的主要拓扑关系。实验表明,时空动态分段模型弥补其他模型将时空参考分为时间参考和空间参考所产生的不足,而且更明确地将多重属性和物理实体有机结合起来。     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

A mathematical model for steady‐state regolith production at constant erosion rate

It has been hypothesized that many soil profiles reach a steady‐state thickness. In this work, such profiles were simulated using a one‐dimensional model of reaction with advective and diffusive solute transport. A model ‘rock’ is considered, consisting of albite that weathers to kaolinite in the presence of chemically inert quartz. The model yields three different steady‐state regimes of weathering. At the lowest erosion rates, a local‐equilibrium regime is established where albite is completely depleted in the weathering zone. This regime is equivalent to the transport‐limited regime described in the literature. With an increase in erosion rate, transition and kinetic regimes are established. In the transition regime, both albite and kaolinite are present in the weathering zone, but albite does not persist to the soil–air interface. In the weathering‐limited regime, here called the kinetic regime, albite persists to the soil–air interface. The steady‐state thickness of regolith decreases with increasing erosion rate in the local equilibrium and transition regimes, but in the kinetic regime, this thickness is independent of erosion rate. Analytical expressions derived from the model are used to show that regolith production rates decrease exponentially with regolith thickness. The steady‐state regolith thickness increases with the Darcy velocity of the pore fluid, and in the local equilibrium regime may vary markedly with small variations in this velocity and erosion rate. In the weathering‐limited regime, the temperature dependences for chemical weathering rates are related to the activation energy for the rate constant for mineral reaction and to the ΔH of dissolution, while for local equilibrium regimes they are related to the ΔH only. The model illustrates how geochemical and geomorphological observations are related for a simple compositional system. The insights provided will be useful in interpreting natural regolith profiles. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.     

模拟地震波场的伪谱和高阶有限差分混合方法

伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算.本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的伪谱/有限差分混合方法.该方法在一个空间坐标方向上利用交错网格高阶有限差分算子,在另外的空间坐标方向上利用交错网格伪谱法算子,既保留了后者的高效、高精度优势,又便于在PC集群上实现并行计算.对二维模型的计算显示,混合方法能有效处理介质不连续面,在保证伪谱法计算精度的情况下,提供了一种并行计算的可能途径.     

基于COMSOL MULTIPHYSICS的直流电法正演模拟

以点电源和线电源为例,叙述了基于COMSOL MULTIPHYSICS的直流电法正演模拟过程。通过对COMSOL MULTIPHYSICS正演模拟解与理论值的对比分析,验证了应用COMSOL MULTIPHYSICS进行直流电法正演模拟的正确性与可行性。鉴于COMSOL MULTIPHYSICS强大的剖分求解功能及丰富的后处理操作,可以快速、准确实现正演模拟,其在地球物理正演模拟中的应用的研究意义重大。