内蕴大地边值问题

提出内蕴大地边值问题,使得有可能利用重力场边界观测研究地球重力场的内蕴结构。文中构造了椭球问题的迭代逼近求解程式,并给出了具体解式。     

华东地区具有短期预测意义的特殊显著性地震事件识别

依据华东地区中强地震前的地震活动特点,初步提出“特殊显著性地震事件”SPE的概念及其可能的识别方法.具有前兆指示意义的SPE可划分为三类,三类SPE分别表达了主震前中小地震与主震之间的时空位置关系.I类SPE距主震时、空距离均较近,在一定程度上与“直接前震”具有相似之处;II类SPE距主震时间较近、而空间距离可能较远,但与主震之间有一定的构造联系,具有时间上的“信号地震”的意义.Ⅲ类必须首先通过地震活动性方法对可能的SPE进行筛选,之后通过比较其与之前在同一位置发生的一系列中小地震的地震波参数之间的差异,在一定程度上进一步确认其短期预测意义.华东6次中强地震震例的研究结果表明,当可能的SPE的单台波速比小于此前该区域长时期的均值、归一化尾波持续时间大于此前该区域长时期的均值以及振幅衰减系数与此前该区域长时期均值之差大于长时期均值的30％时,可进一步认定其短期前兆意义.     

TTI介质qP波逆时偏移中伪横波噪声压制方法

在对地下复杂构造介质,特别是盐丘侧翼及岩下区域进行成像时,相对于传统的各向同性逆时偏移和VTI逆时偏移,具有倾斜对称轴的TTI逆时偏移成像效果最优.不仅反射同相轴更加的连续,而且能量得到了更好的聚焦.传统的各向异性介质全弹性波RTM的计算量大且计算效率低.由于目前仍以纵波勘探为主,因此TTI逆时偏移qP波波动方程的选取显得尤为重要.为了提高计算效率,采用将沿着对称轴方向的横波速度设为零的方法,简化得到qP波波动方程.然而,这样会引入一种严重影响成像效果的低速度、低振幅的qSV波人为干扰.本文建立了qP波方程的完全匹配层控制方程,而后借助于辅助波场采用一种高效的压制伪横波噪声传播的方法,通过模型测试验证了该方法的有效性.     

凸度量空间中渐近拟非扩张映象新的带误差的Ishikawa迭代逼近

在凸度量空间中，利用新的解析法，征明了由集合序列{On}产生的在多值渐近拟非扩张映象下的新的带误差的lshikawa迭代序列的收敛性。     

围岩参数Monte－Carlo有限元反分析

提出了应用隧道施工监测的随机值，采用蒙特一卡罗（Ｍｏｎｔｅ－Ｃａｒｌｏ）有限元随机反分析以反推地层参数和初始地应力分布的方法。反推参数不仅给出其均值，还给出参数的方差及分布类型，为地下结构稳定性和可靠性分析提供必要数据。为检验所提的反分析方法对工程的实用性，文中给出了一个实例。     

An enhanced ensemble Kalman filter scheme incorporating model error in sequential coupling between flow and geomechanics

In this work, we construct a new methodology for enhancing the predictive accuracy of sequential methods for coupling flow and geomechanics while preserving low computational cost. The new computational approach is developed within the framework of the fixed-stress split algorithm procedure in conjunction with data assimilation based on the ensemble Kalman filter (EnKF). In this context, we identify the high-fidelity model with the two-way formulation where additional source term appears in the flow equation containing the time derivative of total mean stress. The iterative scheme is then interlaced with data assimilation steps, which also incorporate the modeling error inherent to the EnKF framework. Such a procedure gives rise to an “enhanced one-way formulation,” exhibiting substantial improvement in accuracy compared with the classical one-way method. The governing equations are discretized by mixed finite elements, and numerical simulation of a 2D slab problem between injection and production wells illustrate the tremendous achievement of the method proposed herein.     

Algorithms for Stellar Perturbation Computations on Oort Cloud Comets

We investigate different approximate methods of computing the perturbations on the orbits of Oort cloud comets caused by passing stars, by checking them against an accurate numerical integration using Everhart’s RA15 code. The scenario under study is the one relevant for long-term simulations of the cloud’s response to a predefined set of stellar passages. Our sample of stellar encounters simulates those experienced by the Solar System currently, but extrapolated over a time of 10¹⁰ years. We measure the errors of perihelion distance perturbations for high-eccentricity orbits introduced by several estimators – including the classical impulse approximation and Dybczyński’s (1994, Celest. Mech. Dynam. Astron. 58, 1330–1338) method – and we study how they depend on the encounter parameters (approach distance and relative velocity). We introduce a sequential variant of Dybczyński’s approach, cutting the encounter into several steps whereby the heliocentric motion of the comet is taken into account. For the scenario at hand this is found to offer an efficient means to obtain accurate results for practically any domain of the parameter space.     

On the local multiscale determination of the Earth’s disturbing potential from discrete deflections of the vertical

As a first approximation, the Earth is a sphere; as a second approximation, it may be considered an ellipsoid of revolution. The deviations of the actual Earth’s gravity field from the ellipsoidal “normal” field are so small that they can be understood to be linear. The splitting of the Earth’s gravity field into a “normal” and a remaining small “disturbing” field considerably simplifies the problem of its determination. Under the assumption of an ellipsoidal Earth model, high observational accuracy is achievable only if the deviation (deflection of the vertical) of the physical plumb line, to which measurements refer, from the ellipsoidal normal is not ignored. Hence, the determination of the disturbing potential from known deflections of the vertical is a central problem of physical geodesy. In this paper, we propose a new, well-promising method for modelling the disturbing potential locally from the deflections of the vertical. Essential tools are integral formulae on the sphere based on Green’s function with respect to the Beltrami operator. The determination of the disturbing potential from deflections of the vertical is formulated as a multiscale procedure involving scale-dependent regularized versions of the surface gradient of the Green function. The modelling process is based on a multiscale framework by use of locally supported surface curl-free vector wavelets.      

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.