模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre (GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性.

     

Estimating groundwater recharge for a freshwater lens in an arid region: Formative and stability assessment

The formation of subsurface freshwater lenses on top of brackish groundwater is a fascinating hydrologic phenomenon that creates groundwater supplies of great potential value in arid regions. Information on the recharge quantity and mechanism of these lenses is both scarce and uncertain. This study examines the formation and macroscale stability of the Rawdatain freshwater lens in Kuwait, for which significant pre-development data are available. The Rawdatain is a large (150 million m³) subsurface freshwater lens overlying brackish groundwater compared to the other freshwater lenses in the Arabian Peninsula. In this study, a three-dimensional (3D) density-dependent groundwater flow model is tested against the following data targets to estimate long-term diffuse and focused groundwater recharge: (i) groundwater head, (ii) total dissolved solids (TDS) groundwater concentration, (iii) volume and vertical thickness of stored groundwater of three different water quality TDS ranges (0–700, 700–1000 and 1000–2000 mg/L) and (iv) geometrical shape features of the lens along cross-sections. To better represent the spatial variation in TDS, six different recharge zones were assigned to allocate diffuse and focused recharge conditions. Twelve recharge rate scenarios, encompassing a wide range of feasible long-term average annual recharge values (200,000–5,000,000 m³/year), were tested against the multiple targets and compared with the groundwater age of the Rawdatain lens. Based on comparison with data targets, the long-term average annual recharge is estimated to be 500,000 m³/year. Scenarios of reduced recharge, which may occur due to changes in land-use or climate, demonstrate the extremely slow response of the lens, which is in agreement with the slow development and formation of the lens (>2,000 years). Within a 100-year time frame, a 50% reduction in annual recharge reduces the lens volumes by 21, 17 and 9% for the three water quality categories, respectively. This study demonstrates the stability of freshwater lenses in arid regions and also provides methodology for similar focused rainfall recharge freshwater lenses.     

适用于声波方程数值模拟的时间-空间域隐式有限差分算子优化方法

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.     

无单元法用于地震波波动方程模拟与成像

波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二乘的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的.     

储层建模中几种原型模型的建立

储层地质知识库的建立是储层建模过程中极为重要的一步，它为建模提供基本的条件数据及各种统计参数。本文总结了建立地质知识库的主要步骤及基本内容。在资料较少的情况下，需要借助原型模型来提供一些无法直接获得但对建模十分重要的参数，如砂体宽厚比、水平方向变差函数等。文中用3个实例较详细的介绍了建立原型模型的3种主要方法，分别是利用露头、成熟油田和沉积模拟实验建立原型模型，并比较了各自的优缺点。     

空间波数混合域磁异常场积分解三维数值模拟

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证.

     

东亚地区区域气候模拟的研究进展

由于大气环流模式对东亚地区区域气候特征的模拟存在很大不足,采用区域气候模式模拟该地区特殊的季风气候成为目前发展的一个重要研究方向,并取得了显著的研究进展。通过回顾当前东亚地区区域气候模拟的现状,表明大部分区域气候模式都不同程度地模拟出了东亚季风区持续性洪涝现象的大尺度环流特征和演变过程,再现了东亚各主要气候区降水的年际、季节和季节内变化及主要雨带的季节进退和降水的时空演变特征。但是,大部分模式没能很好地模拟出大尺度特征的强度和量值,模拟的温度和降水存在系统性偏差。原因分析表明,进一步完善和改进区域气候模式的物理过程参数化方案及动力框架可以改善模拟效果。最后对区域气候模式未来的发展给出展望。     

THE APPLICATION OF IMAGE-BASED MODELING IN PALEOEARTHQUAKE TRENCH STUDY

The geological structure exposed by paleoearthquake trenches is the key material to the right cognition of fault activity and paleoearthquake. However, paleoearthquake trenching inevitably destroys active tectonic geomorphic evidence and trench exposures are usually difficult to reserve. The conventional process of recording the delicate geological information, manually constructing photomosaics by image-editing software, is time-consuming and produces undesirable artificial distortions. Herein, we explored the process of constructing trench orthophotomosaics and the 3D image model using the Image-based Modeling technology and applied it to the Liutiaohe trench across the Tianqiaogou-Huangyangchuan Fault, Gansu Province. Based on the 3D image modeling and orthophotomosaic, we firstly constructed the control points and scale bars on cleaned trench walls and collected photos of all sections of the trench with a digital camera in the field, and then reconstructed the 3D model of the trench through the Agisoft PhotoScan, an efficient image-based modeling software, and finally yielded the 3D image model of the trench and othophotomasaics of the trench exposures. The results show that the automated workflow can produce seamless, sub-millimeter-level high-resolution photomosaics more quickly, with precision in the centimeter range, and the 3D image model is of great help to identify strata and geological structures in trenches with much lower capital and labor costs and low expertise levels compared with LiDAR, meanwhile, the 3D archive benefits the share and communication and even allows future reinterpreting the site using new insights.     

Modeling uncertainties of reservoir flushing simulations

Every dam or barrage construction affects the watercourse and the retention of sediment that previously was carried by the river, which can lead to siltation of the reservoir and obstruction of water intakes over time, reducing their capacities. However, the information available regarding the effect of sediment and drawdown parameters, sediment management at reservoirs, as well as different equational approaches, is scarce. The current research aims to evaluate the effect of parameters associated with numerical modeling of sediment management in reservoirs considering scenarios with different drawdowns, transport equations, sediment size distributions, and thickness of the initial sediment layer. The case study of the Aimorés Hydropower Plant (HPP) is used, applying the Delft3D-FLOW model for two-dimensional modeling. All parameters influenced the volume of mobilized sediment, among which the initial layer thickness was the parameter that resulted in the greatest changes in simulated results. In general, the results show that the uncertainties in the input parameters outweigh the uncertainties between the techniques, which found large variations in results when evaluating the use of different transport equations. These results indicate the importance of proper estimation of model parameters for predicting effects with accuracy and the need for such studies before planning and management operations are evaluated to avoid environmental harm and energy waste.     

一次湖北暴雪天气的诊断与模拟

利用NCEP GFS资料分析了2007年1月15—16日鄂东南地区降雪过程,对造成暴雪过程的天气系统发生、发展背景场进行分析。并利用中尺度数值天气模式WRF模拟了这次暴雪过程,探讨了其发生发展的机制。天气系统的背景分析表明,这次暴雪过程主要是受700 hPa西南急流和地面冷空气的共同影响而产生的,降水过程与西南急流的变化密切联系。WRF模式较好地再现了此次暴雪的过程。模拟结果表明西南急流的减弱和移出,对应着降雪的开始和停止;在西南急流的左侧,由于低层涡度的增加,使低空辐合、高空辐散,在连续性原理和动力机制约束下导致上升运动的加强是该次暴雪的形成机制。模式结果说明,产生暴雪的上升运动要远小于产生暴雨的上升运动,且在暴雪过程中,中层为上升运动,近地层和高层伴随着下沉运动。