地震时间序列的周期图分析方法和极大熵谱分析方法

利用华北地震区的地震资料（1000～1977）,对地震时间序列进行了周期图分析和极大熵谱分析,提取了隐含在随机噪声中的隐含周期。结果表明,这两种方法所得结果很接近,地震应变能√Ei;释放的隐含周期利用周期图分析提取的主要周期是30年、245年和305年,极大熵谱分析方法提取的主要周期是30年、240年和300年。这三个主要周期存在“倍30”的关系,240年和300年周期是30年周期的整数倍。30年周期反映了各活跃幕中的幕式活动规律,300年周期同各活跃期高潮时段之间的时间差的平均值相符,华北自1000年以来的四个活跃期高潮时段之间的时间差平均值是300年。     

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

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

A wetting and drying algorithm with a combined pressure/free-surface formulation for non-hydrostatic models

A wetting and drying method for free-surface problems for the three-dimensional, non-hydrostatic Navier–Stokes equations is proposed. The key idea is to use a horizontally fixed mesh and to apply different boundary conditions on the free-surface in wet and dry zones. In wet areas a combined pressure/free-surface kinematic boundary condition is applied, while in dry areas a positive water level and a no-normal flow boundary condition are enforced. In addition, vertical mesh movement is performed to accurately represent the free-surface motion. Non-physical flow in the remaining thin layer in dry areas is naturally prevented if a Manning–Strickler bottom drag is used. The treatment of the wetting and drying processes applied through the boundary condition yields great flexibility to the discretisation used. Specifically, a fully unstructured mesh with any finite element choice and implicit time discretisation method can be applied. The resulting method is mass conservative, stable and accurate. It is implemented within Fluidity-ICOM [1] and verified against several idealized test cases and a laboratory experiment of the Okushiri tsunami.     

红板岩材料隐式本构模型的建立及其应用研究

针对红板岩材料在岩土工程中所表现的大量模糊的和不确定的因素等特点,基于人工神经网络的学习能力,借助于室内岩石力学试验,进行了对该材料的力学本构特性进行了神经网络模拟研究,提出了隐式本构模型的思想和方法,并通过该方法对该岩石的流变试验结果进行学习,获得了以网络权值结构保存的力学特性知识,由此得到了表征红板岩应力应变本构关系的隐式本构模型。应用结果表明,该方法对岩土类材料本构关系的模拟研究具有很好的应用前景。     

结合边界条件的截断高阶隐式有限差分方法(英文)

有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。     

隐式等权重粒子滤波在高维准地转模式中的特性研究

粒子滤波自从被引入资料同化领域以来,对于高维系统存在的粒子衰退问题一直困扰着资料同化领域的研究。隐式等权重粒子滤波（Implicit Equal-Weights Particle Filter,IEWPF）通过在高维的状态空间维数的前提下,隐式从每个粒子都具有特殊协方差的提议密度中进行采样,构建等权重的粒子集合,从而解决高维系统的粒子衰退问题。通过在高维准地转模式中应用IEWPF方法,验证了IEWPF的系统一致性和资料同化效果。通过对水平动能谱的检验,验证了IEWPF可以保持系统的原始平衡特性。通过IEWPF与等权重粒子滤波（Equivalent Weights Particle Filter,EWPF）的对比试验发现,两者的资料同化分析场非常接近,但在运行效率上,IEWPF远优于EWPF。同时,IEWPF也为解决一系列的资料同化问题,比如参数估计,提供了新的解决途径。      

Method and analysis for the upscaling of structural data

3D geological models are created to integrate a set of input measurements into a single geological model. There are many problems with this approach, as there is uncertainty in all stages of the modelling process, from initial data collection to the approach used in the modelling scheme itself to calculate the geological model. This study looks at the uncertainty inherent in geological models due to data density and introduces a novel method to upscale geological data that optimises the information in the initial dataset. This method also provides the ability for the dominant trend of a geological dataset to be determined at different scales. By using self-organizing maps (SOM's) to examine the different metrics used to quantify a geological model, we allow for a larger range of metrics to be used compared to traditional statistical methods, due to the SOM's ability to deal with incomplete datasets. The classification of the models into clusters based on the geological metrics using k-means clustering provides a useful insight into the models that are most similar and models that are statistical outliers. Our approach is guided and can be calculated on any input dataset of this type to determine the effect that data density will have on a resultant model. These models are all statistical derivations that represent simplifications and different scales of the initial dataset and can be used to interrogate the scale of observations.     

Stochastic simulations of fault networks in 3D structural modeling

3D structural modeling is a major instrument in geosciences, e.g. for the assessment of groundwater and energy resources or nuclear waste underground storage. Fault network modeling is a particularly crucial step during this task, for faults compartmentalize rock units and plays a key role in subsurface flow, whether faults are sealing barriers or drains. Whereas most structural uncertainty modeling techniques only allow for geometrical changes and keep the topology fixed, we propose a new method for creating realistic stochastic fault networks with different topologies. The idea is to combine an implicit representation of geological surfaces which provides new perspectives for handling topological changes with a stochastic binary tree to represent the spatial regions. Each node of the tree is a fault, separating the space in two fault blocks. Changes in this binary tree modify the fault relations and therefore the topology of the model.     

Formulation of discontinuous deformation analysis (DDA) — an implicit discrete element model for block systems

Continuum-based numerical methods have played a leading role in the numerical solution of problems in rock mechanics and engineering geology. However, for fractured rocks, a continuum assumption often leads to difficult parameters to define and over-simplified geometry to be realistic. In such case, discrete representations of fractures and individual blocks must be adopted. In this paper, a newly emerged member in the family of discrete element methods (DEM), the discontinuous deformation analysis (DDA), is presented, including its variational principle, governing equations, solution techniques and contact representation and detection algorithms. Its relative advantages and shortcomings are compared with the explicit distinct element method and the finite element method. An example of the analysis of tunnel stability is provided to demonstrate the capability of this new method.     

越洋海啸波传播数学模型的研究

应用内外解匹配的方法和驻相法推导了柱坐标系下地震引起的水面波动方程的解,即推导了地震海啸波生成与传播的理论方程,并对越洋地震海啸理论方程进行了求解。基于Boussinesq方程出发,建立了二维越洋海啸波传播数学模型,并对越洋海啸进行了数值计算,计算方法采用有限差分方法,差分格式采用交替方向隐格式(Alternating Direction Implicit即ADI方法)。利用越洋海啸计算模式对发生在大洋中的地震海啸进行了模拟,将数值模拟结果与地震海啸波理论方程的计算结果进行了比较,两种计算结果吻合较好。