航空物探勘查金属矿产时提取弱信息的解释方法

摘 要　 结构逻辑法是一种提取航空物探局部弱信息的新的解释方法。在介绍该方法原理的 基础上‚探讨了利用这种方法提取航空物探局部弱信息的效果及其在矿产资源快速勘查评价 中的作用。     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.     

基于GRACE的空间约束方法监测华北平原地下水储量变化

华北平原作为我国重要的工农业基地和政治经济中心,面临着严重的水资源危机.因此,开展对华北平原地下水储量变化的监测工作具有重要现实意义与科学价值.本文基于GRACE重力卫星的空间约束方法,研究了华北平原地下水储量变化的时空分布规律,并与地面水井实测与地下水模型结果进行了综合比较和分析.结果表明:2002-2014年,华北平原地下水存在明显的长期亏损,GRACE估计的亏损速率为-7.4±0.9 km~3·a~(-1),而地面水井资料估计的浅层地下水亏损速率为-1.2 km~3·a~1,对比两者之间的差异可以发现,华北平原的地下水亏损以深层地下水为主.2002-2008年,GRACE估计的华北平原地下水亏损速率为-5.3±2.2 km~3·a~(-1),这与华北平原两个地下水模型得到的平均亏损速率-5.4 km~3·a~(-1)十分吻合.通过华北平原区域地下水模型的独立验证,说明GRACE可以有效评估华北平原的地下水储量变化趋势.除了长期亏损的趋势项之外,华北平原地下水还存在明显的年际变化特征,并与该地区年降雨量变化特征一致.在降雨偏少的2002年、2005-2009年和2014年,华北平原地下水储量显著减少.在空间分布上,GRACE结果表明,华北平原的地下水储量减少主要发生在山前平原和中部平原区,这也与水井实测资料和区域地下水模型结果较为吻合.与GRACE和区域地下水模型相比,目前的全球水文模型仍无法准确估计华北平原地下水变化的空间分布和亏损速率.上述研究表明,GRACE提供了评估华北平原地下水储量变化的重要监测手段.     

Late Holocene slip rate and average recurrence interval of great earthquakes on the Shangzhi segment of the Yilan‐Yitong Fault Zone,northeastern China: Constraints from paleo‐earthquakes and historical written records

The Yilan‐Yitong Fault Zone (YYFZ) is considered to be the key branch of the Tancheng‐Lujiang Fault Zone (TLFZ) in northeastern China. Although the Mesozoic and early Cenozoic deformation of the YYFZ has been studied intensively over the past century, few estimates of slip rate and recurrence interval of large earthquakes in the late Quaternary, which are the two most important parameters for understanding the potential seismic hazard of this crucial structure, were obtained. Based on integrated interpretations of high resolution satellite images and detailed geologic and geomorphic mapping, linear landforms were identified, including fault scarps and troughs, along the Shangzhi segment of the YYFZ, which exceeds 25 km in length. Synthesized results of trench excavations and differential GPS measurements of terrace surfaces indicate that two events (E1, E2) occurred along the Shangzhi segment during the late Holocene, which resulted in 3.2 ±0.1 m of total vertical co‐seismic displacement with clear features of thrust motion. ¹⁴C dating of samples suggests that event E1 occurred between 440 ±30 years BP and 180 ±30 years BP and that event E2 occurred between 4 090 ±30 years BP and 3 880 ±30 years BP, which indicates that the minimum vertical slip rate of the Shangzhi segment of the YYFZ has been approximately 0.8 ±0.03 mm/year during the late Holocene. Constraints from paleo events and the slip rate suggest that the average recurrence interval of major earthquakes on the YYFZ is 3 800 ±200 years. Historical documents in Korea show that event E1 possibly corresponds to the earthquake that occurred in AD 1810 (the Qing Dynasty in Chinese history) in the Ningguta area, which had surface‐wave magnitude (Ms) of 6.8–7.5. Studies of kinematics show that the right‐lateral strike‐slip with a reverse component has been dominant along the YYFZ during the late Holocene.     

强造斜致绳索取心钻具折断分析

在定向钻探固体矿产领域,通常在造斜段使用螺杆钻具,大部分的工作量采用取心回转钻进。造斜段使用螺杆钻具,钻杆不回转,受力较简单;如遇强造斜转换绳索取心钻具稳斜钻进,下部钻具组合容易出现回转困难、下不到底的技术难题,甚至出现丝扣折断事故。为避免强造斜对绳索取心钻具的影响,本文通过建立强造斜段对绳索取心钻具的几何约束变形模型,量化分析强造斜段岩心外管长度与造斜强度、环空间隙之间的约束关系,提供了可快速预判下部钻具参数能否满足强造斜安全钻进的方法。该方法在山东某铁矿定向孔的强造斜致下部钻具折断的事故处理中,通过对绳索取心钻具安全计算模型校核和调整钻具参数,实现了绳索取心钻具通过造斜率0.62°/m的强造斜段安全钻进至靶点,证明对强造斜绳索取心工具安全钻进预测和事故处理计算具有指导意义。     

基于BP神经网络的薄互层储层预测

根据地震属性来进行储层预测的研究,对于寻找油气具有十分重要的意义。而由于地下地质情况千变万化,不确定的因素太多,对这种多自变量与因变量的复杂关系模拟问题,神经网络技术是目前较成熟、实际应用效果也较好的方法之一。用BP神经网络来预测薄互层储层厚度,它可以建立属性参数与预测目标之间的高度非线性映射,并应用于具体的实例,为薄互层储层预测提供了新的思路。     

Incorporating known petrophysical model in the seismic full-waveform inversion using the Gramian constraint

In this paper, we develop a general approach to integrating petrophysical models in three-dimensional seismic full-waveform inversion based on the Gramian constraints. In the framework of this approach, we present an example of the frequency-domain P-wave velocity inversion guided by an electrical conductivity model. In order to introduce a coupling between the two models, we minimize the corresponding Gramian functional, which is included in the Tikhonov parametric functional. We demonstrate that in the case of a single-physics inversion guided by a model of different physical type, the general expressions of the Gramian functional and its gradients become simple and easy to program. We also prove that the Gramian functional has a non-negative quadratic form, so it can be easily incorporated in a standard gradient-based minimization scheme. The developed new approach of seismic inversion guided by the known petrophysical model has been validated by three-dimensional inversion of synthetic seismic data generated for a realistic three-dimensional model of the subsurface.     

结点虚流量法理论基础阐释及改进算法

含自由面的无压渗流问题本质上是一类非线性自由边值问题,固定网格的结点虚流量法在全域范围内不断扣除虚域流量贡献,从而使该问题得到求解。它具有网格依赖性小、出逸点收敛快等优点,但其内在理论基础尚未被完全揭示。通过引入互补型约束条件建立了结点虚流量法和Signorini型变分不等式提法的等价性桥梁,在此基础上引入过渡区放大系数 对自由面判别准则进行优化,并以砂槽模型试验为例进行验证。对比结果表明,优化后算法数值稳定性更好,计算结果与试验数据吻合度更高。研究成果为超大规模网格的渗控结构优化设计提供了有效分析手段。     

厦门花岗岩残积土物理力学指标关联性定量表征初探

厦门地区花岗岩残积土的颗分试验表明其粒度呈“两头多,中间少”的特征,级配较独特。孔隙比等物性指标与压缩模量等设计中常用指标定量关联性不强,离散性大且没有明显规律。结合厦门地铁工程地质勘察资料,通过粒间状态参量、级配控制模式确定所研究土体的关键物理状态参量,分析该参量与岩土工程设计中常用指标之间的相关关系。结果表明,引入Thevanayagam提出的粒间孔隙比作为参变量,压缩模量随粒间孔隙比的增大而减小,建立了估计残积砾质黏性土、残积砂质黏性土压缩模量的经验公式。级配控制模式下的特征粒径比可以作为估计花岗岩残积土抗剪强度指标的关键参变量,黏聚力随特征粒径比的增大而增大,内摩擦角随特征粒径比的增大而减小,且线性相关性较好。文中建立的预测花岗岩残积土压缩模量、抗剪强度指标经验公式,可供厦门地区工程设计时参考。     

Midlittoral polychaete communities in the Eastern Mediterranean Sea: new information from the implementation of the Natural Geography in Shore Areas (NaGISA) protocol and comparisons at local and regional scales

The aim of this study was to attempt to fill current knowledge gaps on midlittoral Mediterranean biodiversity at local and regional scales, by using benthic polychaetes as a model taxon. Two different data sets were analysed: (i) a quantitative data set from the two Natural Geography in Shore Areas (NaGISA) study sites in Crete and (ii) a qualitative data set from multiple sites across the Mediterranean. At the local scale, the results provide evidence that (i) discrete species communities are formed in midlittoral Mediterranean habitats, which vary by geographical location and year, depending on the scale of observation; (ii) macrophyte coverage and Chl‐a are the only environmental variables associated, albeit weakly, with the above pattern; (iii) although naturally disturbed, the Cretan NaGISA sites do not seem to experience any anthropogenic stress; (iv) environmental heterogeneity and history seem to be much less important in shaping the polychaete communities than inter‐specific interactions; however, it is not possible to specify at this stage whether local or regional processes or even their interactions may shape the polychaete communities. At the regional scale, the results indicate that (i) the only factor that seems to be involved in the regional pattern is the identity of the study providing the data sets, which implies variability and bias in how research projects are carried out, from the sampling design through to data collection and analysis; (ii) the Cretan NaGISA sites may be considered as representative of the habitat in the Mediterranean, under certain conditions.