基于交错网格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.     

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

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

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

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

Receiver function imaging in strongly laterally heterogeneous crust: Synthetic modeling of BOLIVAR data

We resolve a large (~20 km) discrepancy in Moho depth determined from PdS receiver functions (RFs) and from active source seismic profiling in the complex Caribbean-South American plate boundary zone in eastern Venezuela. As part of the BOLIVAR experiment 20 broadband stations were deployed along an active source profile to record teleseisms. Using the extremely heterogeneous crustal model obtained from active source data, we generated 2D finite-difference elastic wave synthetics and from them calculated receiver functions and CCP stacks. We compare the observations with synthetic sections that have been spatially sampled at 0.25 km to 40 km. The densely sampled synthetics show that several events in the field data that were originally interpreted as the Moho are multiple reflections within sedimentary basins. Where the Moho has the steepest dip under the plate boundary the CCP stacks fail to image the Moho well, regardless of the density of spatial sampling. A suitable spatial sampling criterion for clearly imaging the lower crust and Moho is to overlap Fresnel zones by 50% at Moho depth, which for the 1 Hz receiver functions examined here, requires an instrument spacing of 15–20 km, with the actual field data density ranging from 20 km to 100 km.     

基于ObjectARX技术的管道参数化建模

基于AutoCAD 2008,结合ObjectARX 2008软件包,在.NET 2005环境下开发了管道参数化建模系统。阐述了此系统的数据文件的组织,所采用的关键技术以及实现的基本功能。     

复杂系统理论与Agent模型在土地变化科学中的研究进展

基于复杂系统理论与Agent 模型,分别从理论与方法的角度梳理并总结了当前土地变化科学的研究进展。复杂系统理论为研究土地变化所处的“耦合的人类-自然系统”及其复杂性问题提供理论支撑;而基于Agent 的模拟作为研究复杂系统的重要工具,为传统土地变化模拟提供新的方法支持。当前,ABMs与土地变化模型相结合(ABM/LUCC),无论是模型构建还是模型应用均取得了显著进展。然而,其也存在诸多问题,尤其是理论与应用脱节,使很多建模者容易陷入“为建模而建模”的误区。此外,现有的ABM模型设计较为随意,模型间差异太大,不利于跨区域比较研究的开展。因此,虽然ABM是认识复杂系统的有效方法,但是建模之前离不开对系统复杂性的深入研究;基于Agent 的土地变化模拟的意义在于充分表达“人类-自然”系统的综合复杂关系,从而合理的预测土地变化(结构与功能) 的动态过程,解释并评价土地变化的可能影响。     

Fault type analysis along the San Andreas Fault zone: A numerical approach

Finite Element (FE) modeling under plane stress condition is used to analyze the fault type variation with depth along and around the San Andreas Fault (SAF) zone. In this simulation elastic rheology was used and was thought justifiable as the variation in depth from 0.5 km to 20 km was considered. Series of calculations were performed with the variation in domain properties. Three types of models were created based on simple geological map of California, namely, 1) single domain model considering whole California as one homogeneous domain, 2) three domains model including the North American plate, Pacific plate, and SAF zone as separate domains, and 3) Four domains model including the three above plus the Garlock Fault zone. Mohr-Coulomb failure criterion and Byerlee's law were used for the calculation of failure state. All the models were driven by displacement boundary condition imposing the fixed North American plate and Pacific plate motion along N34°W vector up to the northern terminus of SAF and N50°E vector motion for the subducting the Gorda and Juan de Fuca plates. Our simulated results revealed that as the depth increased, the fault types were generally normal, and at shallow depth greater strike slip and some thrust faults were formed. It is concluded that SAF may be terminated as normal fault at depth although the surface expression is clearly strike slip.     

金融危机下国家大规模投资政策模拟

构建了一个中国多区域可计算一般均衡模型(MReCGE-C),利用该模型分析了我国应对国际金融危机的大规模投资政策。MReCGE-C根据城乡属性和年龄结构对居民进行了分组,放松了要素流动的限制,使劳动力和资本可以同时跨区域跨部门流动,并设计了一种基于居民终生效用的区域均衡机制,通过区域差距变量可以调控和测度区域差距。最后利用MReCGE-C模型针对国家大规模投资计划下的3种典型投资方案进行了模拟。结果发现在等量新增投资下,针对中部地区投资方案的GDP增长率高于针对中西部和全国投资方案的GDP增长率,同时针对中部地区投资方案对于城乡就业的拉动也要优于其余两种投资方案,而且可以实现缩小区域差距的目标。综合比较,目前国家大规模投资的优势应该指向中部地区。     

Assessing the natural component of the external load due to organic matter and biogenic elements on the water bodies in Northwestern Russia

We suggest the method to calculate the natural component of the external loads due to organic matter and biogenic elements on the unexplored lakes using the removal rate moduli for these elements, and indices of the lakes’ specific catchments. The schematic map has been generated for the distribution of natural moduli of removal on the catchment of the Gulf of Finland. We calculated the external natural loads for 94 lakes in landscapes with a high percentage of lakes.