Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Multiple overpressure systems and their relationship with hydrocarbon accumulation in the Qikou Depression of the Bohai Bay Basin,China

Pressure measurements using drill stem tests and estimates from log data calculation indicate that three vertically stacked regional pressure compartments exist in the Qikou Depression of Bohai Bay Basin, N. China. The compartments comprise hydrostatic, upper weak, and lower overpressure systems. Laterally, overpressure (pressure coefficient > 1.2) occurs in the deeper areas and weakens gradually from the centre to the margin of the depression. The accumulation of oil and gas exhibits the interesting characteristics of oil‐bearing layers above gas‐bearing layers in the Qikou Depression. The pattern can be accounted for by the evolution of overpressure system, the maturity process of the source rock and the main fault activity. In the late Dongying Formation (Ed, 30 Ma), the lower overpressure system began to form shape, and the hydrocarbon sources generated a large volume of oil. However, because there was no migration pathway, the oil only accumulated in the original strata. In the late Guantao Formation (Ng, 12 Ma), the gas was generated, the upper overpressure system formed gradually, and the activity of the main fault gradually increased. Then, the overpressure pushed the early gathered oil to flow from the lower overpressure system into the upper overpressure system. Afterwards, the activity of the main fault decreased again and remains weak until now. Thus, later generated natural gas cannot keep migrating along the main fault and can only accumulate in the lower overpressure system. Copyright © 2014 John Wiley & Sons, Ltd.     

基于FME的国家基本比例尺图幅坐标系转换研究

我国测绘成果长期存在坐标系统不统一的现状,使得不同坐标系测绘成果的使用存在一定的局限性,无法充分地利用测绘成果。本文以1∶10 000比例尺成果数据从西安80系到国家2000系转换为例,详细地介绍了基于FME软件零编程实现坐标系转换的方法。     

EMD-Wavelet-ICA耦合模型及其在GPS坐标序列降噪中的应用

在分析比较经验模态分解(EMD)、小波变换(Wavelet)和独立分量分析(ICA)优缺点的基础上,提出一种新的EMD-Wavelet-ICA耦合模型。该模型充分利用了EMD的自适应性,对原始信号进行分解获得不同频率的模态函数(IMF),采用标准化模量的累计均值对IMF进行尺度划分;进而分别采用Wavelet和ICA对高频和低频IMF进行降噪,将降噪后的IMF进行多尺度重构,获得降噪后的信号;采用信噪比、标准差、偏差和相关系数等指标对降噪效果进行评价。仿真数据和GPS坐标序列的处理结果表明:与EMD模型和EMD-ICA模型相比,新模型的标准差、偏差均有不同程度的减小;信噪比和相关系数有一定程度的增大,可以获得更好的降噪效果。     

局部区域模型重力异常快速算法

地球重力场位系数模型可以用于计算局部重力扰动场元。然而随着地球重力场模型阶次的提高、局域重力场计算范围的增大,其计算速度往往不能满足工程需求。针对这一问题,在对位系数模型泰勒级数展开的基础上提出了采用向量运算、混合编程的方法,同时对连带勒让德函数Belikov递推方法中与经纬度无关的量进行了预先计算,有效提高了计算速度。提出的方法对于利用超高阶次重力场模型快速解算大范围、高分辨率重力场元数据以及累加求和计算具有一定的参考与借鉴意义。     

Provenance of Holocene beach sand in the Western Iberian margin: the use of the Kolmogorov–Smirnov test for the deciphering of sediment recycling in a modern coastal system

Detrital zircons from Holocene beach sand and igneous zircons from the Cretaceous syenite forming Cape Sines (Western Iberian margin) were dated using laser ablation – inductively coupled plasma – mass spectrometry. The U–Pb ages obtained were used for comparison with previous radiometric data from Carboniferous greywacke, Pliocene–Pleistocene sand and Cretaceous syenite forming the sea cliff at Cape Sines and the contiguous coast. New U–Pb dating of igneous morphologically simple and complex zircons from the syenite of the Sines pluton suggests that the history of zircon crystallization was more extensive (ca 87 to 74 Ma), in contrast to the findings of previous geochronology studies (ca 76 to 74 Ma). The U–Pb ages obtained in Holocene sand revealed a wide interval, ranging from the Cretaceous to the Archean, with predominance of Cretaceous (37%), Palaeozoic (35%) and Neoproterozoic (19%) detrital‐zircon ages. The paucity of round to sub‐rounded grains seems to indicate a short transportation history for most of the Cretaceous zircons (ca 95 to 73 Ma) which are more abundant in the beach sand that was sampled south of Cape Sines. Comparative analysis using the Kolmogorov–Smirnov statistical method, analysing sub‐populations separately, suggests that the zircon populations of the Carboniferous and Cretaceous rocks forming the sea cliff were reproduced faithfully in Quaternary sand, indicating sediment recycling. The similarity of the pre‐Cretaceous ages (>ca 280 Ma) of detrital zircons found in Holocene sand, as compared with Carboniferous greywacke and Pliocene–Pleistocene sand, provides support for the hypothesis that detritus was reworked into the beach from older sedimentary rocks exposed along the sea cliff. The largest percentage of Cretaceous zircons (<ca 95 Ma) found in Holocene sand, as compared with Pliocene–Pleistocene sand (secondary recycled source), suggests that the Sines pluton was the one of the primary sources that became progressively more exposed to erosion during Quaternary uplift. This work highlights the application of the Kolmogorov–Smirnov method in comparison of zircon age populations used to identify provenance and sediment recycling in modern and ancient detrital sedimentary sequences.     

依兰-伊通断裂新构造活动规律分析

作为郯庐断裂带北段主干的依兰-伊通断裂, 其新构造活动性与活动规律仍然存在不同的认识.本次工作通过详细的野外调查, 发现该断裂内活断层广泛存在, 由东、西两支北东走向的主干活断层构成, 沿着古近纪地堑边界断层发育.这些活断层主要呈破碎型结构, 多为逆右行平移活动.通过对这些活断层一系列实测擦痕反演应力场, 显示它们多是在东西向挤压中活动的, 而现今应力场转变为北东东-南西西向区域性挤压.依据本次野外观察与¹⁴ C定年, 并结合前人定年结果与近代地震分布, 表明依兰-伊通西支活断层的最新活动时代为全新世与晚更新世相间, 而东支活断层的最新活动时代主要为早-中更新世.依兰-伊通断裂内活断层显示了明显的差异性活动, 表现为西支的活动强度明显大于东支, 西支的最新活动时代皆晚于东支, 沿走向上活动性强、弱相间与最新活动时代不断变化, 以及近代地震活动不均一分布.它们沿走向上的分段性、差异性活动主要是因为被一系列北西向断层切断所致.     

内蒙古敖汉旗七家金矿闪长玢岩锆石U-Pb年代学、地球化学及地质意义

内蒙古敖汉旗七家金矿构造位置上处于兴蒙造山带与华北克拉通北缘的结合部位, 靠近兴蒙造山带一侧, 是近年来新发现的石英脉型金矿床.区内岩脉以闪长玢岩为主, 且与石英脉型矿体空间关系密切.根据岩脉与矿体的穿插关系, 将其分为成矿前闪长玢岩和成矿后闪长玢岩.LA-ICP-MS锆石U-Pb年代学表明, 成矿前闪长玢岩加权平均年龄为166.3±3.4 Ma, 成矿后闪长玢岩加权平均年龄为128.6±4.5 Ma, 与区域上两期构造-岩浆活动时间(分别为160 Ma左右和128 Ma左右)近似吻合.主量元素特征显示: 二者均属于准铝质, 高钾钙碱性系列; 微量元素特征显示: 二者均富集轻稀土元素、大离子亲石元素(Rb、K)和活泼的不相容元素(U、Pb), 相对亏损重稀土元素和高场强元素(Nb、Ta、Ti), 反映出岩浆来源具壳源特征.主微量元素特征均显示出二者具有岛弧岩浆岩的地球化学属性.另外, 成矿后闪长玢岩有较高的Nb/Ta(18.66~20.27)和Zr/Hf(37.16~39.23)比值, 暗示其岩浆来源具明显的幔源特征.岩石地球化学特征表明, 成矿前闪长玢岩可能起源于太平洋板块俯冲导致的下地壳部分熔融岩浆; 成矿后闪长玢岩可能起源于俯冲流体交代的地幔楔熔融岩浆, 并在其演化侵位过程中伴有地壳物质的混染.结合闪长玢岩成岩时代、岩石成因及区域构造演化, 反映出七家金矿的形成与太平洋板块的西向俯冲密切相关, 矿床在地壳由挤压增厚向伸展减薄的转换过程中形成.