Albian extrusion evidences of the Triassic salt and clues of the beginning of the Eocene atlasic phase from the example of the Chitana-Ed Djebs structure (N.Tunisia): Implication in the North African Tethyan margin recorded events, comparisons

In the northern part of Tunisia, close to Testour/Slouguia, new observations and updated biostratigraphy make it possible to highlight the relation between the Triassic saliferous mass and the surrounding Mesozoic beds (T. M.). Near the (T. M.) boundary, the formations observed consist dominantly of Triassic evaporites reworked in the Early and Late Albian deep-water sedimentary deposits. Throughout the studied area, Jurassic rocks are absent. We propose to interpret the Chitana-Ed Djebs structure originally emplaced as gravitational stretch masses in a passive margin in the same way as the salt bodies of widespread salt province in the Gulf of Mexico. A reconstructed schematic position of the Chitana-Ed Djebs salt body displays a scenario of setting of the salt mass on a submarine palaeo-slope. Moreover, the starting clues of the paroxysmal event of the Late Mesozoic tectonic inversion clearly fossilized through the discordance of the Middle Eocene–Early Lutetian limestone on the Albian series.     

Identification of diffusion parameters in a nonlinear convection–diffusion equation using the augmented Lagrangian method

Numerical identification of diffusion parameters in a nonlinear convection–diffusion equation is studied. This partial differential equation arises as the saturation equation in the fractional flow formulation of the two-phase porous media flow equations. The forward problem is discretized with the finite difference method, and the identification problem is formulated as a constrained minimization problem. We utilize the augmented Lagrangian method and transform the minimization problem into a coupled system of nonlinear algebraic equations, which is solved efficiently with the nonlinear conjugate gradient method. Numerical experiments are presented and discussed. This work was partially supported by the Research Council of Norway (NFR), under grant 128224/431.     

基于PCA-GEP算法的边坡稳定性预测

提出一种基于主成分分析的基因表达式程序设计算法,并将其用于边坡稳定性预测。该算法先采用主成分分析法对样本数据进行预处理,有效地减少预测模型的输入量,消除输入数据间的相关性,再将得到的新样本数据输入基因表达式,构建边坡稳定性的预测模型。利用该预测模型对82个危险圆弧破坏边坡实例中的71个实例进行学习,对另外11个实例进行预测,取得了较好的效果。在保留传统的以误差值作为评判模型优劣标准的同时,引入AIC信息准则法,分别对v-SVR算法和GA-BP网络算法和PCA-GEP算法三种预测模型进行比较分析,结果表明,运用该算法可以获得更优的预测模型,其预测结果比v-SVR算法和GA-BP网络等其他算法得到的结果具有更高的预测精度。工程实例计算表明,该方法是合理、可行的。     

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.     

起伏地表条件下地震波场数值模拟有限积分变换有限差分方法

针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。     

岩土体渗流自由面问题的重心插值无网格方法

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。     

The earthquake foci and deep structure of the Earth’s crust and upper mantle along the southern Sakhalin-Sea of Okhotsk-Kamchatka profile

The Kuril-Kamchatka seismofocal zone was thought to be a single plate approximately 90 km wide and dipping to a depth of 700 km at an angle of 40°–45°. This concept reflects primarily the physical differences (elastic wave velocities, density, temperature, etc.) between the seismofocal zone and the mantle hosting it. Detailed investigations show that the seismofocal zone proper is also heterogeneous with earthquake hypocenters variably concentrated and clustered within this zone, where both seismogenic and aseismic strata, as well as subvertical zones, can be identified. The latter are reflected in the structure and faults of the Earth’s crust and upper mantle.     

西南山区矿调中的地球物理综合找矿技术

为了寻找隐伏矿、深部矿,以四川通安—小牛场地区为研究平台,开展了矿产远景调查,提出了将重力、航磁、地磁、激电等物探手段有机结合,综合应用于西南山区矿产调查的思路,并列举典型实例,逐步分析各种物探异常,最终实现从圈定找矿重点区到确定矿体平面分布形态的目的.以期这一重、磁、电相互结合的模式对今后西南山区的矿产勘查能起到一定的指导作用.     

转换波反射系数近似公式的比较及分析

多分量地震资料解释的首要问题就是多分量的标定和对比,其中最关键的是转换波反射系数的求取及子波的选取,然而Zoeppntz方程在数学表达上极其复杂,在物理意义上也不直观,因此很难得到直接的应用.为此,在详细介绍五种近似公式的基础上,利用六种模型分析了随入射角的变化各公式的适应性,最后得出结论,Aki和Richard的近似公式适应性普遍比较好.     

盆地超压层段非幕式突破期的地热场模型数值解法

盆地超压层段的地热场模拟对超压成因研究及油气生成，排放机理分析有重要意义，但求解超压层段的热传递方程的数值解法至今仍没有被很好地解决，问题的关键在于确定超压流体速度场，由于超压流体发生幕式突破之前的排出速率极低，其速度场v=vxi vzk近似满足条件Эux/Эx Эuz/Эz＝0，可将其近似地视为稳定的不可压缩的无源流体，利用这一条件及相应的边界压力条件，可使整个计算过程得到简化，在此基础上所建立的地热场模型有限元数值解法和模拟软件，能够实现对含油气盆地在幕式突破之前的超压层段地热场进行动态模拟。