A new model of edge-to-edge contact for three dimensional discontinuous deformation analysis

A robust contact theory can be regarded as key to three dimensional discontinuous deformation analyses (3D-DDA). Not only must this theory provide an efficient algorithm to judge the type and location of contacts but also it must be able to present comprehensive formulations for every kind of contact (open, sliding and locked contact). There are six types of contact in three dimensional discontinuous deformation analyses (vertex-to-vertex, vertex-to-edge, vertex-to-face, edge-to-edge, edge-to-face and face-to-face) that can be converted to vertex-to-face and edge-to-edge contacts. This paper presents a new model of edge-to-edge contact to three dimensional discontinuous deformation analyses (3D-DDA). This new model considers both kinds of edge-to-edge contact (cross-over and parallel edge-to-edge contact) and presents a criterion for inter-penetration. Sub matrices of normal and shear spring and friction force are derived by geometrical analysis and penalty method. This new model is implemented in a 3D-DDA computer programme, and the example results demonstrate the validity of the model.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

现代海底块状硫化物矿床储量估算方法探讨

海底块状硫化物矿床具有巨大的经济价值和良好的开发前景。但是由于勘探难度和技术手段的限制,在全球已经发现的600多个海底热液喷口和硫化物矿床中,不到5%的矿床进行了详细的勘探,进行资源量计算的矿床更是屈指可数。选择适当的方法估算储量对于深海资源开发意义重大。本文系统总结了前人估算海底硫化物储量的方法,并以勘探程度相对较高、数据资料相对丰富的Atlantis-II-Deep矿床和Solwara 1矿床为例,评析了前人采用不同储量计算方法得到的结果和各种方法计算的准确性;针对有一定钻孔数据的TAG硫化物堆丘,采用距离反比法重新计算了储量,并和块段法估算的结果进行对比。得到以下结论:(1)海底硫化物矿体形态是制约储量估算方法适用性的一个重要因素;(2)传统几何法可以在极低勘探条件下估算海底硫化物储量量级;(3)对于海底硫化物储量估算,地质统计学的精度高于传统几何法;(4)较高勘探条件下距离反比法可以与地质统计学相互对照、验证;(5)在一定钻孔数量的条件下,可以应用距离反比法进行海底硫化物储量估算。     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

聚合物驱组合段塞分注实验研究

萨尔图油田南二区东部140号断层以东地区葡11—4油层，以多段多韵律和正韵律沉积为主，平均孔隙度30．0％，平均空气渗透率1980×10^-3μm2。水驱结束时该油层平均含水92％，采出程度39．5％。计划对该区块实施聚合物驱开采。以往经验认为，组合段塞调整注聚技术能够有效改善流度控制作用。然而对组合注入情况下如何进一步提高驱油效果、注入井分注时机、层段注入强度对开发效果有何影响却很少探讨。通过室内驱油实验，在聚驱的组合段塞注入不同时机（水驱空白阶段、含水下降阶段、含水稳定阶段、含水回升阶段）下调整注入强度实施分注以及在含水回升分注基础上进一步调整注入强度，同时对比不分注实验，评价各方案驱油效果。实验表明：不实施分注时聚驱采收率在水驱基础上提高了19．07个百分点，而组合段塞注入时实施分注均能够在不分注基础上进一步提高驱油效果；空白水驱时实施分注采收率提高幅度最大，为26．55个百分点；随着分注时机的滞后采收率提高幅度增加值依次递减，分注时机越早越好。     

大气中一类非线性方程的非线性波解及周期的近似计算法

大气中许多模式可化为Ｆ＇＇＝ＡＦ／Ｆ＋Ｂ（ＡＢ≠０）一类非线性常微分方程。利用常微分方程定性理论说明了该方程存在非线性波解，讨论了其周期的计算问题，并求得了波解及其周期的近似表达式。     

利用原生岩浆定量反演原岩微量元素丰度的方法

利用同源不同程度平衡部分熔融作用形成的两种原生岩浆岩的微量元素浓度，提出了一种定量反演原岩微量元素丰度的方法。反演前不需要做任何人为假设，为定量研究源区的微量元素地球化学提供了一个有力的手段。对汉诺坝新生代橄榄拉斑玄武岩和碱性玄武岩的成因及地幔源区特征作了讨论。反演结果显示本区地幔相对球粒陨石有过明显的ＲＥＥ富集过程，并且不同的ＲＥＥ富集程度存在显著的差异。     

Forward Calculation of 2-D and 3-D Structures with a Cover by the Boundary Element Method while Using Electrical Methods

This paper develops the boundary element method,the authors employtwo-layered earth Green’s functions as the weighting functions of residual andderive boundary integral equations.The forward problems of point sources on2-D and 3-D structures with an influencing cover are solved by this method.The results show that this method markedly improves the original boundaryelement method.The features of the improved method are greater numericalaccuracy and much smaller systems of equations and thus considerable savingsfor the storage capacity of computers,allowing us to solve the above problemswith only ordinary microcomputers.The results in this paper extend the scopeof applying the boundary element method while using electrical methods forgeophysical prospecting.     

裂隙岩石介质中核素迁移的模型及其算子分裂、迎风、均衡格式

建立了放射性核素在裂隙岩石介质中迁移的双重介质模型,对模型的求解提出了一种新的数值方法—Galerkin有限元法与算子分裂、迎风、均衡格式相结合的新方法,给出了水质模型算子分裂、迎风、均衡格式的稳定性条件,且所得到的计算格式是非负的。最后通过对核素90Sr 100年、99Tc 1000年的预测计算,验证了本文所提方法的有效性和稳定性,并得出了一些有重要意义的结论。