基于二维经验模式分解的重力资料多尺度分析

重力资料是从地表到深部所有密度不均匀体引起的重力效应的叠加。针对不同的工作任务,需要利用不同尺度的重力资料进行研究。二维经验模式分解能够自适应地对非线性、非平稳信号实现多尺度分解。文章研究将二维经验模式分解应用于重力资料的多尺度分析,对重力数据进行二维经验模式分解以得到各级固有模态函数和剩余分量,利用径向对数功率谱分析方法估计各级模态函数所反映的地下场源的近似深度,定性地解释不同深度下的场源产生的异常,较好地揭示场源赋存的地质信息。模型数据和实际数据的测试佐证了该方法技术的可行性。     

中国地区地磁正常场模型研究

基于2000年中国地区的实测地磁3分量数据,分别选用Taylor模型和第四代地磁场综合模型--CM4计算并绘制出正常场分布图.由于2种正常场分布较为相似,为了进一步比较研究,分别将实测值减去2种正常场模型值以得到2种磁异常值,并利用曲面Spline模型绘制各地磁分量磁异常图,结果显示:2种异常图分布基本一致,以CM4模型作为正常场与国界外的地磁分布结合得较好,适合进行大尺度的磁场研究,计算的磁异常分布在国内以负异常为主;以Taylor模型作为正常场适合进行局部区域的磁场研究,获得的磁异常图正负分布较为均匀、合理.     

经验模态分解和独立分量分解在航磁数据处理中的应用

以东天山1∶20万航磁数据为例,利用经验模态分解(EMD)方法对航磁数据进行分解,然后利用独立分量分析(ICA)方法对分解出的固有模态函数(IMF)数据进行重构。分解出来的从高频到低频IMF与地质体的分布具有一定的相关性,不同频率的IMF代表了一定的地质意义。经过ICA重构后的独立分量(IC)能够对异常IMF函数从盲源分离的角度进行有效的重构,其结果可以在一定程度上解释不同的地质事件对磁异常的贡献,并对构造识别和构造区域划分具有一定的指示意义。     

煤矿深部开采动力灾害区域应力效应监测与前兆特征研究

动力灾害的孕育及诱发与区域应力场的力学行为密切相关，煤矿开采扰动作用下区域局部应力场不断调整变化。为了研究区域应力场变化特征对动力灾害的前兆响应，采用变分模态分解和希尔伯特变换的方法对鲍店煤矿采集数据进行时频分析，识别和提取钻孔应变数据中反映动力异常变化的固有模态分量。研究结果表明，钻孔应变观测数据能够有效记录开采扰动相关区域岩体内部的微小变化，解析原始数据滤去干扰信号后获得表征动力异常的特定频段，分解变换后固有模态分量异常波动特征早于动压事件前2～7天出现，表现出形变稳定变化阶段向快速变化阶段“突跳”、形变快速变化阶段“震荡”以及失稳阶段“跌落”的显著特征。依据地震形变三阶段理论，判别动力灾害发生前形变异常开始时刻和“突跳”—“震荡”—“跌落”作为煤矿动力灾害发生的前兆判据。基于钻孔应变仪对区域应力场的精确观测，构建了应用于煤矿动力灾害超前预警的判据方法，可为煤矿相似开采环境下工作面的安全高效回采提供借鉴和参考。     

二维EMD方法在物化探数据处理中的应用研究——以云南个旧地区为例

文章采用二维EMD方法对云南个旧地区物化探数据进行非线性的多尺度分析与分解,有效分离了物化探异常,揭示深层次找矿信息.证明了二维EMD分解得到的所有分量都具有实际地质意义,它将是一种新的物化探数据处理方法.     

基于二维经验模态分解的重磁异常分离

由于地质过程的复杂性及成矿过程的多期次叠加性,原始重磁异常往往是多种地质因素的混合信息,既包含区域背景异常信息,也包含与矿床(体)、矿化蚀变带以及隐伏岩体等与找矿密切相关的地质要素所引起的局部重磁异常.如何从复杂的叠加重磁异常中分离出具有找矿意义的局部异常,是当前矿产勘查和资源潜力评价工作中面临的难题之一.采用经验模态分解(EMD)方法来分解重磁异常,为提高分解的稳健性提出了用双调和样条插值(BSI)进行包罗面插值的新方法,并以云南个旧地区重磁数据为例,对其进行非线性多尺度分解,实现对区域异常与局部重磁异常的分离,揭示了深层次找矿信息并拓宽了经验模态分解方法的应用领域.      

改进的完备经验模态分解和FWEO能量在南海油气识别中的应用

传统希尔伯特变换(Hilbert-Huang transform, HHT)是一种识别精度较差的时频分析方法，存在端点效应和模态混叠等问题。改进的完备经验模态分解(Improve Complete Ensemble Empirical Mode Docomposition, ICEEMD)可以将复杂的地震信号分解为一系列单分量信号，较好地解决模态混叠问题，但结合希尔伯特变换提取的瞬时振幅和瞬时频率，对处理实际地震数据仍然有严重的端点效应。FWEO(Frequency-weighted Energy Operator)是一种非负频率权重算子，其结合TK能量差分算法和Hirbert变换复分析思想，具有比Hilbert变换更高的时间分辨率。但由于算法本身的原因，FWEO能量只能应用于单分量信号，不能直接应用于复杂的地震数据。因此，这里结合改进的完备经验模态分解方法和FWEO能量分离算法的优点，提取南海某工区实际地震记录的瞬时振幅和频率信息，并将预测结果与测井数据对比，预测吻合程度好、识别精度高、证明该方法可以准确地反映储层特征。     

基于EEMD的地震瞬时属性在砂砾岩中的应用

集总模态分解(EEMD)是经验模态分解(EMD)的一种改进算法,它克服了EMD中的模态混叠现象并可以有效抑制干扰信号。文中对EEMD方法进行理论试算,证实了其方法的改进效果及抗干扰能力,然后使用EEMD方法处理了中国东部某区块的砂砾岩地震资料,得到了不同频率的固有模态函数(IMF)分量;结合地质情况,找出了对砂砾岩敏感程度较高的IMF分量,并使用Hilbert变换计算出瞬时属性。将计算出的IMF瞬时频率与原始数据的瞬时频率对比得出前者对流体界面变化的识别能力高于后者。     

基于平面地质图的地质体三维建模

文章针对大尺度研究区域勘探范围有限、不能直接基于勘探数据构建区域地质体三维结构模型的问题,提出了一种以平面地质图为研究对象的地质体三维建模方法。此方法以平面地质图为基础,利用Section软件在DEM模型(数字高程模型)约束下半自动绘制一系列图切剖面图,将二维剖面映射到三维空间,利用三维剖面上轮廓线的对应性和连续性构建地质体三维模型。通过广西东平地区地质体三维模型构建实例,验证了该方法的可行性。     

中国区域地磁异常的提取及分析

区域地磁异常的提取对于研究地壳区域的、深部的情况是很有益的。为此,采用地磁台站的多分量磁场强度,经过处理获得了区域地磁异常,提供了丰富的地质信息。     

地磁场模型的计算和评述

评述了计算地磁场区域模型的各种数学方法,如:利用多项式方法可以计算地磁场的Taylor多项式模型和Legendre多项式模型,利用球冠谐和分析方法可以计算出磁场的冠谐模型。本文指出了建立区域模型各种方法的优缺点,指出球冠谐和分析方法是计算地磁场三维模型的好方法。     

公路交通荷载作用下分层地基的三维动响应分析

根据高速公路路基的实际情况,建立了三维分层地基模型,分别考虑了刚性路面和柔性路面的情况,利用传递矩阵法获得了置于刚性基岩上的层状地基和分层半无限地基的刚度矩阵,将车辆荷载简化为矩形均布移动荷载,利用Fourier变换方法得到了层状地基系统在交通荷载作用下的三维振动解,并用快速傅立叶逆变换方法得到了数值结果,与模型试验结果进行了比较,显示出了较好的一致性,同时分析了荷载大小、土层性状对响应结果的影响,得出了一些有意义的结论。研究结果可为路面施工设计及路基动力响应分析提供参考。     

Delineation of long-wavelength magnetic anomalies over central India by rectangular harmonic analysis

The method of rectangular harmonic analysis is applied to the geomagnetic field data from central India to isolate long wavelength magnetic anomalies associated with largescale crustal structures. The long-wavelength anomalies have accounted for approximately 20 % of the spatial variability of the residual magnetic field over the International Geomagnetic Reference Field. On the magnetic anomaly map, reflecting the surface expression of longwavelength anomalies, the Tapi-Narmada-Son zone is characterized by a feeble positive anomaly bounded by a strong negative anomaly. The anomaly pattern is believed to be caused by the large-scale undulation in Moho and related variations in the thickness of the lower (basaltic) crust. The other two prominent anomalies, the magnetic low striking northwest and the magnetic high trending east-northeast, appear to be related to the deep structural feature of the Godavari graben and the eastern Rajasthan lineament respectively.     

地磁总强度梯度的计算与分析

不仅要知道地磁场的分布规律,而且要知道地磁场梯度的分布规律。本上根据地磁场的高斯理论和地磁场的球谐模型,计算出中国及其邻近地区地磁总强度的水平梯度(南北方向和东西方向)和垂直梯度,绘制出相应的等值线图,并分析了总强度梯度的分布规律。地磁总强度的水平梯度和垂直梯度,主要随纬度变化,与经度关系不大。值得注意的是在我国的中心地区,地磁总强度的水平梯度值最大。     

均衡重力异常计算方法

本文提出了采用埃利均衡补偿模式计算5ˊ×5ˊ平均均衡重力异常的一种方认。计算是使用5ˊ×5ˊ平均高程值,将A～O带分为近、中、远三个区:近区用1ˊ40〞×1ˊ40〞的柱体,中区用5ˊ×5ˊ的柱体,远区用15ˊ×15ˊ的柱体,分别采用矩形柱体、质线和质点的正演公式计算地形一均衡校正值。从华南某地区的均衡重力异常图可以看出,根据均衡重力异常可以消除由地壳厚度变化而产生的区域场,突出局部异常场。     

垂直界面附近椭圆形夹杂对SH波的散射与地震动研究

利用复变函数法和Green函数法给出了SH波对垂直界面附近椭圆形夹杂散射问题的解析解答。首先,将待求的半空间模型沿垂直边界分割为区域Ⅰ和区域Ⅱ两个直角域。通过保角映射的方法将区域Ⅰ内椭圆形夹杂的外域映射为单位圆外域,并利用镜像方法构造出两个区域内满足直角域边界条件的散射波场及适用的Green函数;其次,利用界面契合思想,通过在界面处添加附加力系的方法建立起满足界面处位移和应力连续条件的无穷代数方程组,并截断有限项求解;最后,给出了求解地表位移幅值的具体算例。结果表明,入射波数、入射角度、夹杂位置、垂直界面以及材料参数都对地表位移幅值的分布有影响。     

Aeromagnetic data analysis using the 2D directional continuous wavelet transform (DCWT)

In this paper, we have analyzed the geomagnetic anomaly field reduced to the pole by the 2D directional continuous wavelet transform (DCWT). The objective is to identify geological contacts. The proposed idea has been applied firstly at a synthetic model, the obtained results showed the robustness of the proposed technique, after that the synthetic data has been noised; analysis shows that contacts identified by DCWT are affected by noise. To resolve this problem, we have proposed an algorithm to reduce the noise effect in the contacts model. Application on the real geomagnetic data of In Ouzzal area located in the West of Hoggar (Algeria) shows clearly the robustness of the proposed method. Comparison with the analytic signal solutions exhibits a difference between the two models of contacts. Our method gives results that are accurate with the structural geological map.     

Geomagnetic field anomalies over the Lombok Island region: an attempt to understand the local tectonic changes

During the last years, several investigations on the earthquakes and related tectonic structures along the Java Trench have been conducted. In this study, we focus on the Lombok Island—West Nusa Tenggara (Indonesia), which lies between the centres of the highest intensity of contiguous negative–positive geomagnetic anomalies in this region. The geological and tectonic structures, however, are not known in detail for this island, whereas a better knowledge of these conditions could lead to better estimate the potential earthquake risks and thus a possible improvement of the existing early warning system. We have performed a ground-based geomagnetic survey at 56 stations in the Lombok Island during October–November 2006 for a detailed mapping and interpretation of geomagnetic anomalies related to the geological and tectonic characteristics. The 2D and 3D magnetic maps show a general geomagnetic anomaly pattern in the Lombok Island which consists of repeated contiguous negative–positive anomalies. Two forward models have been proposed for a profile connecting a strongest apparent dipolar structure. The first model assumes a uniformly magnetized sphere as the source of the anomaly, and could be interpreted as a specific local structure composed by a quite large magnetic body. The second model considers several rocks with different susceptibilities and magnetizations, and could be interpreted as a discontinuity in the geological structures. This model agrees with the local geological surface conditions and the known large scale regional tectonic structure. Therefore, it is used to interpret our results in terms of tectonic characteristics, which suggests the potential existence of a new tectonic element (e.g. a local normal fault) in this region.     

F-layer formation in the outer core with asymmetric inner core growth

Numerical calculations of thermochemical convection in a rotating, electrically conducting fluid sphere with heterogeneous boundary conditions are used to model effects of asymmetric inner core growth. With heterogeneous inner core growth but no melting, outer core flow consists of intense convection where inner core buoyancy release is high, weak convection where inner core buoyancy release is low, and large scale, mostly westward flow in the form of spiraling gyres. With localized inner core melting, outer core flow includes a gravity current of dense fluid that spreads over the inner core boundary, analogous to the seismic F-layer. An analytical model for gravity currents on a sphere connects the structure of the dense layer to the distribution of inner core melting and solidification. Predictions for F-layer formation by asymmetric inner core growth include large-scale asymmetric gyres below the core-mantle boundary and eccentricity of the geomagnetic field.     

基于无人直升机平台的航磁系统集成与应用

笔者介绍了无人直升机航磁系统的研发与集成,采用无人直升机作为飞行平台,搭载高精度航空磁测系统,具有低成本、高效率、不受机场跑道限制、可夜航、按设计测线全自主导航飞行等特点。该系统分别完成了磁补偿试验飞行与实际应用工作,补偿精度达到0.046 9 n T,测量成果与测区内以往航磁成果对比,其反映的地磁场特征形态基本一致,验证了该系统的有效性,航磁异常等值线在细节上表现更细致。本系统为大比例尺、高精度、小面积的航磁测量工作提供了一种高效灵活的工作手段。