非线性结构时程分析的高阶单步法

本文提出了一种新的求解运动方程的直接积分法－高阶单步β法。     

Annotated bibliography of methodology for assessment of undiscovered oil and gas resources

An annotated bibliography of methodology of assessment of undiscovered oil and gas resources is presented as a useful reference for those engaged in resource assessment. The articles that are included deal only with quantitative assessment of undiscovered or inferred resources. the articles in this bibliography are classified largely according to the major assessment method that was applied in each situation. Major assessment methods include areal and volumetric yield methods, field size distributions, historical extrapolation, deposit modeling, organic geochemical mass balance methods, and direct expert assessment. Other categories include mathematical tools, reserve growth/confirmation, quantitative characterization of undiscovered resources, and general topics. For the purpose of future updates, we solicit contributions of articles that may have been missed in the preparation of this bibliography.     

二维控制场中控制点像点坐标的自动提取与识别

本文给出了一种二维控制场中控制点像点坐标自动提取与识别的原理以及实现方法，并编制程序试验，取得了较高的精度和良好的识别效果。主要实现方法包括直接线性变换、Hough变换、最小二乘拟合等，特别适用于控制点数量较多的大型二维控制场。     

用交错网格有限差分法计算三维频率域电磁响应

用交错网格有限差分法(SFD)，实现了三维频率域电磁场响应 的数值模拟. 该方法适用于任何方向的磁偶极子源. 经与解析方法、积分方程等 其他方法的计算结果对比表明，交错网格有限差分法结合散度校正和不完全乔累斯基分解预 处理的双共轭梯度迭代方法进行正演计算，速度快、精度高、结果稳定，能适应三维复杂介 质的数值模拟，为三维电磁反演奠定了基础.     

Microscopic observation and contact stress analysis of granite under compression

To study the damage process of microscale to macroscale in coarse-grained granite specimen under uniaxial compressive stress, we have observed micro-damage localization and propagation by using a newly developed experimental system that allows us to observe the damaging process continuously.

The results showed that pre-existing microcracks lead to macroscopic shear fracture through the damage development process. The mechanism of micro-damage initiation in a granite specimen under uniaxial compressive stress may be considered for two cases. One is that two grains such as quartz and feldspar contact each other in the same direction as the axial stress, and the other is that a biotite grain inclined to the axial stress direction is surrounded by feldspar grains. The homogenization theory was applied to verify numerically the micromechanics of stress-induced damage in the mineral contacts. Local stress distribution in the periodic-micro structure was also calculated by the homogenization theory. It is shown that this analysis, which takes into account the initial state of the specimen, is well adapted to the behavior of two grains for which microcracking is the fundamental mechanism of damage.     

基于Metropolis优化的叠前全局迭代地质统计学反演方法

叠前地质统计学反演将随机模拟与叠前反演相结合,不仅可以反演各种储层弹性参数,还提高了反演结果的分辨率.基于联合概率分布的直接序贯协模拟方法可以在原始数据域对数据进行模拟,不需要对数据进行高斯变换,拓展了地质统计学反演的应用范围;而联合概率分布的应用确保了反演参数之间相关性,提高了反演的精度.本文将基于联合概率分布的直接序贯协模拟方法与蒙特卡洛抽样算法相结合,参考全局随机反演策略,提出了基于蒙特卡洛优化算法的全局迭代地质统计学反演方法.为了提高反演的稳定性,我们修改了局部相关系数的计算公式,提出了一种新的基于目标函数的优化局部相关系数计算公式并应用到协模拟之中.模型测试及实际数据应用表明,该方法可以很好的应用于叠前反演之中.     

Forward plane-wave electromagnetic model in three dimensions using hybrid finite volume–integral equation scheme

We present a concept of the hybrid finite volume–integral equation technique for solving Maxwell's equation in a quasi-static form. The divergence correction was incorporated to improve the convergence and stability of the governing linear system equations which pose a challenge on the discretization of the curl–curl Helmholtz equation. A staggered finite volume approach is applied for discretizing the system of equations on a structured mesh and solved in a secondary field technique. The bi-conjugate gradient stabilizer was utilized with block incomplete lower-upper factorization preconditioner to solve the system of equation. To obtain the electric and magnetic fields at the receivers, we use the integral Green tensor scheme. We verify the strength of our hybrid technique with benchmark models relative to other numerical algorithms. Importantly, from the tested models, our scheme was in close agreement with the semi-analytical solution. It also revealed that the use of a quasi-analytical boundary condition helps to minimize the runtime for the linear system equation. Furthermore, the integral Green tensor approach to compute at the receivers demonstrates better accuracy compared with the conventional interpolation method. This adopted technique can be applied efficiently to the inversion procedure.     

利用微震数据对震源下方构造直接成像

现有的微震记录直接成像方法是将微震记录既当作入射记录,也当作散射记录,从而实现偏移成像.但此方法并不能突出透射波所携带的来自震源下方的深层散射波信息.本文在假设已知微震位置与子波的前提下,提出了对微震下方构造进行逆时偏移的成像方法.该方法类似于常规的逆时偏移,只是震源位置在地下.这使得在成像时,地下更深部的入射波场相比震源在地表时会更为精确,因此能够获得更加准确的成像结果.该方法会给成像结果带来一种尾波高频干扰：地下的震源发出的上行波与上方介质作用后,所产生的多级散射波会干扰反传波场.对此,在成像过程中,对入射场和散射场都进行左右行波分离,以压制该噪声.而在子波信息未知,无法重构入射场时,使用了激发时间成像条件,也能够实现同等效果的偏移成像,且不会出现尾波高频干扰.利用数值实验验证了本文方法的有效性.     

“一带一路”背景下巴基斯坦投资环境及区位选择

“一带一路”倡议以及中巴经济走廊建设为中国同巴基斯坦的经贸合作带来新的机遇,从发展趋势看,中国企业在巴基斯坦投资空间呈逐渐扩大之势;中国企业在巴基斯坦如何进行有效的区位选择,是文章要解决的问题。首先,从基础竞争力、核心竞争力、辅助竞争力3个维度,计算了巴基斯坦2000—2018年整体投资环境。结果显示：整体投资环境年均增长率36.05%,其中基础竞争力年均增长率28.32%,最为显著;表明巴基斯坦投资环境不断改善,特别是基础设施领域适合中国企业投资。其次,梳理了现有中国企业进入巴基斯坦投资的区位选择。结果表明,目前中国企业在信德省和旁遮普省的投资项目较多。最后,分别从趋于经济中心区、边界地区、社会联系密切区和安全形势良好区4个角度,分析了该国的主要省区和城市投资区位。结果表明,从趋于经济中心区的角度,信德省和旁遮普省存在经济上的优势;从趋于边界地区的角度,开伯尔-普什图省和瓜达尔港具有优势;从趋于社会联系密切地区角度,卡拉奇、拉合尔和伊斯兰堡具有优势;从趋于形势安全的角度,旁遮普省和信德省存在投资安全的优势。