康滇地区重力场分离、密度反演与地壳构造

康滇地区位于我国著名的南北向活动构造带的南端，是中蒙经向中轴强地震带和巨大的地壳厚度变化带的一部分．由于该地区是地震活动的预测危险区，且资源非常丰富，引起了广大地学工作者的关注．本文采用对地壳内部密度分布最为敏感的重力场资料进行场分离、平差、反演等数学处理，求出地壳平均密度和深度10～100km各层的密度分布；综合其它地球物理资料进行分析．结果表明，康滇菱形块体与其东侧和西侧不论在构造特征还是地球物理场特征方面皆存在差异．      

基于最小反演拟合差的重磁场源深度计算方法

以等效源及位场物性反演为基础,本文提出一种新的求取重磁场源深度的方法.该方法将一层等效源以一定的间隔从浅部向深部移动,并将等效源作为初始模型进行反演,当反演拟合差最小时,停止反演,此时的等效源底深即为所求场源的中心深度.由于仅需要反演一层等效源,比传统的物性反演计算时间大大减少,并且不需要进行深度加权约束.理论模型数据处理结果表明该方法能够获得较准确的场源深度:以长宽比为7.5的薄板模型为例,深度计算误差约为1个点距(25 m);以长宽比为0.5~1.5的厚板模型为例,深度计算误差小于1个点距(25m).将该方法应用于实测航磁梯度数据,计算的磁源中心深度在200~250m之间,钻井资料显示该异常由埋藏深度在200~300m的闪长岩引起,计算结果与钻井资料较吻合.     

普光气田重力异常的视密度反演

普光气田是我国近年发现的最大气田．根据气田储层和盖层的岩石密度,估算普光气田上约有零点几毫伽量级的重力异常．但由于区域场和其他岩石的影响,气田的负异常受到严重的干扰．我们使用基于位场切割分离和位场大深度向下延拓的视密度反演新方法,在普光气田上得到明显的负密度异常．在川东北其他一些气田上,也得到负密度异常．根据视密度反演结果,对川东北的含气远景区做了预测．     

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

四川地区流动重力资料的位场 分离与异常特征提取

用于监测四川地区重力场随时间变化的流动重力测网于2008年5月12日汶川M_S8.0地震后进行了全面改造, 形成了新的闭合环线并对整网点位每年实施2次重力观测工作. 该文以2010—2012年获取的5期流动重力资料为基础, 分别采用小波多尺度分解和三维视密度反演等位场分离方法, 对测区内重力场的动态变化进行多层次分解及异常特征提取, 以探求四川测区范围内重力异常动态变化特征以及地震前后不同深度范围内密度变化情况与重力场响应之间的关系. 结果表明, 对同一期资料进行小波分解得到不同阶次的细节图像, 揭示了不同深度处重力异常动态变化特征; 分解结果同时也表明了四川测区范围内重力场的变化形态与川内特有的断裂构造格架和深部结构环境及过渡带特征密切相关. 通过对不同时期的差分重力异常资料进行三维视密度反演对比研究发现, 在地震前后四川测区范围内由浅到深密度变化呈现出不同程度的差异, 震区附近出现较为明显的横向密度变化特征且浅层物质密度局部性的异常变化差异较为突出, 测区范围内中上地壳深度层密度反演切片显示出的趋势性差异基本上反映了地震前后测区内地下物质的动态变化趋势, 因此, 深刻认识不同深度范围内重力场的动态演化特征, 对研究壳内物质密度变化以及与地震的孕育和发展有关的深部介质环境具有积极的意义.      

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

地震折射法中波场延拓技术的应用

从波动方程出发,推导出平面波场传播基本公式,阐述了利用地表观测波场反演地球内部波速结构的理论关系.观测波场可通过-p 变换分解为地表平面波场,其最大振幅的轨迹能稳定地反映出地球内部波速随深度变化的趋势,这一特征可用来对反演解空间进行约束.波场延拓技术可以充分地利用观测波场中包含的丰富信息,在迭代反演中不仅能简单快速地得到反演解,而且所得解具有良好的稳定性,较少受主观因素影响,这是一种很好的反演技术.对波场延拓进行了理论分析与数值模拟,采用了同态反褶积等改善资料信噪比的措施,使解的分辨能力得到了提高.文中对南海北部一个声纳折射剖面进行了分析和计算.结果表明:该区1.4km 深处,存在一个从1.76km/s 到2.21km/s 的速度间断面.间断面上下两层的速度梯度分别为0.54kms-1/km,0.63kms-1/km.最后,从构造演化的角度对浅海构造特征进行了探讨.      

视密度反演在东海及邻区重力异常解释中的应用

简要介绍了基于位场分离的视密度反演法的原理,并利用该技术对东海的重力异常进行处理并进行了区域地质解释应用.视密度反演的方法在对研究区内的区域性大断裂如郯庐断裂、江绍断裂以及东引—海礁隐伏断裂有很好的分辨能力.不同深度的视密度切片信息能够帮助我们了解这些大断裂的发育深度和规模.在沟弧盆等地形以及构造复杂地区利用视密度反演法研究深部地质构造,可以揭示不同深度上补偿质量分布情况,从而推测莫霍面的起伏情况.     

位场数据反演中深度分辨率及改进措施

位场勘探是一种较为成熟的地球物理勘探手段,然而其数据反演中存在的深度分辨率不足问题阻碍了其进一步发展.文章在深入研究位场数据反演理论的基础上,详细分析了位场反演中普遍存在的深度分辨率不足问题,分析了造成这种现象的主要客观原因.研究利用深度分辨率成图技术对给定位场反演问题的深度分辨率进行定量刻画.针对造成分辨率不足的原因提出了几点改进措施:利用离散Picard定理分析问题的适定水平,结合傅立叶系数与奇异值对比图进行合理的噪声剔除;反演计算之前结合DRP技术,给出较为合理的地下剖分方式;合理的数据加权.文章使用理论二维重力模型进行说明,最后将其应用于复杂三维重力模型,试验表明上述措施对于改进位场反演深度分辨率有一定的效果.     

3D inversion of gravity data by separation of sources and the method of local corrections: Kolarovo gravity high case study

We present a novel methodology for 3D gravity/magnetic data inversion. It combines two algorithms for preliminary separation of sources and an original approach to 3D inverse problem solution. The first algorithm is designed to separate sources in depth and to remove the shallow ones. It is based on subsequent upward and downward data continuation. For separation in the lateral sense, we approximate the given observed data by the field of several 3D line segments. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. It allows retrieving unknown 3D geometry of anomalous objects in terms of restricted bodies of arbitrary shape and contact surfaces. For restricted objects, we apply new integral equations of gravity and magnetic inverse problems. All steps of our methodology are demonstrated on the Kolarovo gravity anomaly in the Danube Basin of Slovakia.     

基于抛物线密度模型的频率域三维界面反演及其在川滇地区的应用

密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性.     

The inversion of density structure by graphic processing unit （GPU） and identification of igneous rocks in Xisha area

Organic reefs, the targets of deep-water petroleum exploration, developed widely in Xisha area. However, there are concealed igneous rocks undersea, to which organic rocks have nearly equal wave impedance. So the igneous rocks have become interference for future exploration by having similar seismic reflection characteristics. Yet, the density and magnetism of organic reefs are very different from igneous rocks. It has obvious advantages to identify organic reefs and igneous rocks by gravity and magnetic data. At first, frequency decomposition was applied to the free-air gravity anomaly in Xisha area to obtain the 2D subdivision of the gravity anomaly and magnetic anomaly in the vertical direction. Thus, the distribution of igneous rocks in the horizontal direction can be acquired according to high-frequency field, low-frequency field, and its physical properties. Then, 3D forward modeling of gravitational field was carried out to establish the density model of this area by reference to physical properties of rocks based on former researches. Furthermore, 3D inversion of gravity anomaly by genetic algorithm method of the graphic processing unit(GPU) parallel processing in Xisha target area was applied, and 3D density structure of this area was obtained. By this way, we can confine the igneous rocks to the certain depth according to the density of the igneous rocks. The frequency decomposition and 3D inversion of gravity anomaly by genetic algorithm method of the GPU parallel processing proved to be a useful method for recognizing igneous rocks to its 3D geological position. So organic reefs and igneous rocks can be identified, which provide a prescient information for further exploration.     

The inversion of density structure by graphic processing unit (GPU) and identification of igneous rocks in Xisha area

Organic reefs, the targets of deep-water petroleum exploration, developed widely in Xisha area. However, there are concealed igneous rocks undersea, to which organic rocks have nearly equal wave impedance. So the igneous rocks have become interference for future exploration by having similar seismic reflection characteristics. Yet, the density and magnetism of organic reefs are very different from igneous rocks. It has obvious advantages to identify organic reefs and igneous rocks by gravity and magnetic data. At first, frequency decomposition was applied to the free-air gravity anomaly in Xisha area to obtain the 2D subdivision of the gravity anomaly and magnetic anomaly in the vertical direction. Thus, the distribution of igneous rocks in the horizontal direction can be acquired according to high-frequency field, low-frequency field, and its physical properties. Then, 3D forward modeling of gravitational field was carried out to establish the density model of this area by reference to physical properties of rocks based on former researches. Furthermore, 3D inversion of gravity anomaly by genetic algorithm method of the graphic processing unit (GPU) parallel processing in Xisha target area was applied, and 3D density structure of this area was obtained. By this way, we can confine the igneous rocks to the certain depth according to the density of the igneous rocks. The frequency decomposition and 3D inversion of gravity anomaly by genetic algorithm method of the GPU parallel processing proved to be a useful method for recognizing igneous rocks to its 3D geological position. So organic reefs and igneous rocks can be identified, which provide a prescient information for further exploration.     

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.     

鄂尔多斯南缘地区重力变化场源特征

以鄂尔多斯南缘地区布格重力异常数据及2014—2017年相对重力观测数据为基础,采用欧拉反褶积方法,对引起重力变化的场源深度进行反演,并对空间分布规律予以解释。通过构建理论模型,反演得到最优模型参数,并对实际数据进行计算和分析。为了减弱和消除欧拉解的发散性,利用水平梯度滤波法优化反演结果。结果表明:①构造指数为1时,适合对相对重力数据进行反演,当滑动窗口选择4—8倍测点间距时,可以获得较为可靠的场源参数;②在主要断裂附近,不同时间尺度的场源呈较好的一致性特征;③引起重力变化的场源深度集中在10—30 km,与鄂尔多斯南缘地区的震源深度基本一致;④不同时间尺度的场源位置相对分散,表明引起重力变化的物质流动具有随机性和波动性。