Accelerated initial stiffness schemes for elastoplasticity

Iterative methods for the solution of non‐linear finite element equations are generally based on variants of the Newton–Raphson method. When they are stable, full Newton–Raphson schemes usually converge rapidly but may be expensive for some types of problems (for example, when the tangent stiffness matrix is unsymmetric). Initial stiffness schemes, on the other hand, are extremely robust but may require large numbers of iterations for cases where the plastic zone is extensive. In most geomechanics applications it is generally preferable to use a tangent stiffness scheme, but there are situations in which initial stiffness schemes are very useful. These situations include problems where a nonassociated flow rule is used or where the zone of plastic yielding is highly localized. This paper surveys the performance of several single‐parameter techniques for accelerating the convergence of the initial stiffness scheme. Some simple but effective modifications to these procedures are also proposed. In particular, a modified version of Thomas' acceleration scheme is developed which has a good rate of convergence. Previously published results on the performance of various acceleration algorithms for initial stiffness iteration are rare and have been restricted to relatively simple yield criteria and simple problems. In this study, detailed numerical results are presented for the expansion of a thick cylinder, the collapse of a rigid strip footing, and the failure of a vertical cut. These analyses use the Mohr–Coulomb and Tresca yield criteria which are popular in soil mechanics. Copyright © 2000 John Wiley & Sons, Ltd.     

山东莱芜煤田潘西煤矿深部煤层特征及赋煤规律浅析

石炭-二叠纪是山东省的主要聚煤时期之一,通过对潘西煤矿深部煤层特征和赋煤规律的分析、研究,在箕斗状反倾向盆地,"震钻结合,地震先行,钻探验证"、"就煤找煤"寻找石炭-二叠纪煤层的思路,取得了较好的深部找煤效果。     

基于平面波加速的VTI介质最小二乘逆时偏移

地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-PLSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.     

淮南顾北矿F104断层两侧岩溶水化学形成机制及导隔水性评价

利用断层两盘岩溶水化学特征判断其导隔水性对于我国华北型煤田水害防治具有十分重要的实践意义。以顾北矿F₁₀₄断层两侧太原组岩溶地下水为研究对象,在分析断层两盘水文地质条件基础上,采用Piper三线图、离子组合比和主成分法,分析了主要组分来源及水-岩作用差异性,并采用PHREEQC软件对岩溶地下水进行反向水化学模拟。结果表明:南北区岩溶水均存在方解石和白云石的溶解和沉淀现象,南区阳离子交换吸附和脱硫酸作用程度强于北区,而北区黄铁矿氧化和岩盐溶解作用较南区明显,南北两区水化学环境及水-岩作用存在显著差异,进而推断F₁₀₄断层具有较好的阻水性,且影响了其两侧的氧化-还原环境及温度差异,控制着地下水径流方向和水-岩作用程度。      

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

基于CT三维重建与逆向工程技术的煤体数字模型的建立

为了建立能够表征煤体真实孔隙结构的数字模型,进而进行数值模拟研究,通过μCT225kVFCB高精度CT系统扫描得到了大柳塔煤矿长焰煤CT数据,可观测到的最小孔隙直径为1.94 μm。使用基于Matlab语言的三维重建程序结合逆向工程技术,提出了一种将煤体CT三维数据转换为CAD数字模型的方法。以Ansys建立的煤体有限元模型为例,对煤体瓦斯渗流进行了模拟,分析了煤体孔隙内的速度及压力分布规律,并计算了沿X、Y、Z方向的渗透系数。计算结果表明：在微观尺度下（< 100 μm）煤体渗透系数呈现各向异性,其受煤体结构的影响较明显。提出的基于CT三维重建结合逆向工程技术构建的煤体CAD数字模型不仅可进行有限元分析,同时也可以被EDEM等离散元分析软件所使用,拓宽了煤体CT三维数据的应用范围,丰富了煤体在微观尺度上的研究方法。     

基于贝叶斯理论的逐次迭代非线性AVA反演方法

常规AVA三参数反演方法均基于横波速度与纵波速度之比γ为常数这一假设条件,且常被近似地取为0.5。然而在许多情况下γ并不为常数,而是在横向与纵向都渐变。若一概假定γ等于0.5,反演出的岩性参数势必要偏离真实值,因此有必要合理地选择γ。笔者基于贝叶斯理论,提出逐次迭代非线性AVA的反演方法。该方法把γ看成横向与纵向都渐变的反演初始背景,通过给定初始模型计算初始背景γ,并采用逐次迭代的策略求解该反演问题,解决了关于γ的选取问题以及由于引入变γ值而带来的非线性问题,提高了AVA三参数反演结果的精确度。     

Development of experimental apparatus about reverse circulation bit

A set of experimental apparatus on reverse circulation bit are developed, in order to lucubrate the mechanism of the new type reverse circulation bits, and the structure of the bits influencing the ability of taking core and carrying powder.Both the major structure of the equipment and the procession of experiment are described.     

LDPC码及其研究进展

对低密度奇偶校验码(LDPC)的结构、编译码原理进行了介绍，并与turbo码进行了简单的比较，最后对目前LDPC码的几个研究方向及其进展作了较详细的论述，指出LDPC码具有很大的发展前景。