松辽盆地北部徐家围子断陷火山岩分布及天然气富集规律

提出一个新的带通滤波器——延拓回返垂直二次导数,该滤波器可通过调节延拓回返高度及次数来达到调频目的,以突出勘探目标层段的信息.利用该技术对松辽盆地深层勘探程度最高的徐家围子断陷的航磁资料进行目标处理,与钻井资料对比表明,经处理后的航磁资料可以较好地反映火山岩的分布特征.以航磁资料为主结合钻井、地震资料预测了徐家围子三套火山岩的平面分布.指出航磁异常梯度带是火山岩储层的发育区带,近烃源岩或沟通烃源岩的航磁异常梯度带是火山岩天然气的富集区带.该认识对今后松辽盆地的勘探部署有着十分重要的指导意义,对我国其他地区的火山岩勘探也具有借鉴意义.     

位场向下延拓的改进迭代维纳滤波法

根据维纳滤波理论导出的位场向下延拓滤波器为最佳下延滤波器,但因其实现需要已知待求位场和噪声的功率谱而在实际应用中受到限制.针对该问题,本文首先提出一种基于位场径向平均功率谱的位场噪声水平估计方法,进而利用偏差准则求取正则化参数,实现位场正则化向下延拓;然后将位场正则化下延结果的功率谱作为待求位场功率谱的估计初值,采用带修正项的迭代维纳滤波方法来更新对待求位场功率谱的估计,最后提出本文的位场向下延拓改进迭代维纳滤波方法.基于理论重力模型数据及航磁实测数据进行了向下延拓对比试验,结果表明,改进迭代法具有较好的收敛性,且下延精度优于Tikhonov正则化法和递增型维纳滤波法.     

位场向下延拓的迭代最小二乘法

将位场向下延拓作为一种反问题,结合正则化方法和快速傅里叶算法的优点,本文提出了一种解决位场向下延拓问题的迭代最小二乘法.新方法将位场向下延拓问题转化为最小二秉问题,采用梯度下降法对最小二秉问题进行迭代求解,在迭代过程中,利用快速傅里叶变换方法解决了计算量大的问题.新方法对延拓数据中的噪声具有很好的抑制作用.利用不合噪声...     

基于分形径向谱的位场向下延拓截止波数自动确定方法

向下延拓是重磁位场数据处理与解释的一项重要技术,并因其固有的不稳定性而成为研究的热点.为获得稳定向下延拓结果,波数域向下延拓一般通过附加低通滤波器或改造向下延拓因子来完成.由此,滤波器或改造因子的截止波数则是精确向下延拓的关键.本文基于分形修正径向平均功率谱的物理特性,提出一种位场向下延拓截止波数的自动确定方法.基于理论重力模型和航磁实测数据的向下延拓对比实验结果表明：(1)本文所提出的自动确定方法物理意义明确,能快速有效地确定截止波数并进而获得向下延拓正则参数;(2)基于本文方法的正则化向下延拓结果优于改进导数迭代法的向下延拓结果.     

位场向下延拓的最小曲率方法

针对位场向下延拓的不适定性,我们将位场向下延拓视为向上延拓的反问题,提出以位场最小曲率作为约束条件来求解稳定的下延位场.我们将剖面位场向上延拓表达式用傅里叶矩阵的形式表示,以矩阵乘法形式给出延拓的表达式,同时向待反演的下延位场引入最小曲率约束,得到向下延拓的最小曲率解,并利用正交变换给出了更为简洁的频率域解.随后,利用Kronecker积将上述全部结果拓展至三维位场,给出了三维位场向下延拓的最小曲率解.此外,我们将位场数据的填充、扩充问题与向下延拓问题统筹考虑,提出一种新的向下延拓迭代格式,该算法面向实际资料处理需求、无须预扩充或填补数据.下延迭代时,对原始数据直接向下延拓,而空白部分利用上一次下延位场估计的上延值替代其空白值并对其向下延拓,直至获得最小曲率约束下稳定的向下延拓结果.同时,我们也讨论了利用改进L曲线和广义交叉验证(GCV)计算正则参数最优估计的问题.对理论模型和实际航空重力资料进行了向下延拓检验,处理结果表明位场向下延拓的最小曲率方法解能满足实际位场资料对向下延拓处理的需求,具有较高的下延精度.     

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.     

利用微分进化法确定海洋磁场向下延拓中的最优参数

位场的解析延拓是实现不同高度海洋地磁场相互转换的主要途径,是构建海洋三维磁空间背景场模型的关键技术.针对位场向下延拓迭代法中最优正则化参数及最佳迭代次数难以确定问题,尝试引入微分进化法,以正则化参数及迭代次数为种群变量,以延拓结果的熵值为目标函数,以目标函数最小化为搜索准则,实现两种参数的并行全局寻优.采用实测数据对微分进化法在几种常用的迭代法中最优正则化参数及最佳迭代次数的确定进行了分析,与传统L-曲线准则确定的最优正则化参数及多次试验确定的最佳迭代次数进行对比,结果表明：微分进化法确定的最优参数能使三种迭代法取得最佳迭代效果,延拓结果与真实地磁场最为接近,并且该法计算稳定、自适应强,建议在海洋磁场数据向下延拓中应用.     

基于向上延拓的航空重力向下解析延拓解

位场向下与向上延拓之间存在固有的内在联系,向上延拓解算具有稳定可靠的优良特性,本文据此提出了借助向上延拓信息实现航空重力向下延拓稳定解算的两种方法,分别建立了点对点向下解析延拓模型和最小二乘向下解析延拓模型.其核心思想是,依据泰勒级数展开模型,将位场向下延拓解算过程转换为向上延拓计算和垂向偏导数解算两个步骤,通过第一步的处理有效抑制数据观测噪声对解算结果的干扰,通过第二步的处理成功实现向下延拓反问题的稳定解算,较好地解决了向下延拓解算固有的不适定性问题.分析研究了两种解析延拓模型的计算精度及适用条件,利用超高阶位模型EGM2008建立的模拟标准场数据对两种模型解算结果的合理性和有效性进行了数值验证,证明本文新方法实用易行,具有较高的应用价值.     

航空重力测量数据向下延拓方法比较

首先介绍了空中重力异常向下延拓的直接代表法、正则化法、点质量法和基于球内Dirichlet问题的调和解法,分析了它们各自的特点.然后利用这些方法对中国大同航空重力测量数据进行向下延拓计算.通过对延拓结果与地面实测数据的比较分析,比较客观地评价了不同延拓方法的计算精度、可靠性及适应性.实际计算结果表明,直接代表法和正则化法的延拓结果具有较高的精度和可靠性,而球内Dirichlet问题调和解法的精度较差.     

等效源法三维随机点位场数据处理和转换

为了实现曲面随机点位场数据的曲面延拓和转换,以磁异常位场数据为例,采用一组磁偶极子作为等效源,置于观测面下方的一个曲面上,把观测磁异常作为这组磁偶极子所产生磁异常的边界条件,通过求解线性方程组的方法反演磁偶极子磁矩的大小,再根据反演结果正演所要计算的磁异常数据,实现了曲面随机点磁异常位场数据的向上延拓、向下延拓、求导以及化极处理.在数据量较大时,为了提高反演计算的速度,把磁异常数据和磁偶极子分成若干小块,再利用各块磁异常数据分别反演该块数据下方磁偶极子的磁矩,并通过迭代计算来逐步取得更准确的反演结果.模型试验表明,磁异常位场数据向上延拓的均方根误差小于±2nT,向下延拓和化极也可以取得较高的精度,所提出的分块处理方法提高了延拓和转换的速度,实际资料处理给出了曲面随机点航磁异常数据向下延拓和化极的一个例子.     

Gravity and magnetic investigation on the distribution of volcanic rocks in the Qinggelidi area,north‐eastern Junggar Basin (north‐west China)

A detailed investigation on the location of magmatic intrusions in the Carboniferous strata of the Qinggelidi area, north‐eastern Junggar Basin, is presented based on the interpretation of gravity and magnetic data constrained by petrophysical data, seismics and surface geology. The wavelet multi‐resolution analysis based on the discrete wavelet transform is adopted to the regional‐residual separation of gravity and magnetic anomalies. A power spectrum analysis is applied to estimate the source depths corresponding to different scales. A comparative analysis on the characteristics of local gravity and magnetic anomalies improved our understanding of volcanic rock distribution in the Carboniferous strata. Generally speaking, in total 75 anomalies are recognized, among which 23 are inferred to be the responses of basalts, diabases and andesites with high density and strong magnetization. Twelve anomalies are assumed to be caused by andesites, rhyolites and volcanic breccias with medium‐low density and high magnetization. There are still five anomalies that are believed to be generated by volcanic tuffs with low density and weak magnetization. Lastly, four cross‐sections in 3D gravity and magnetic modelling are displayed to provide a more thorough image of volcanic rocks in our study area.     

Gravity and magnetic field characteristics and hydrocarbon prospects of the Tobago Basin

The Tobago Basin, which is located offshore northern Venezuela with a southern margin close to Trinidad and Tobago, has an area of approximately 59,600 km². The Tobago Basin has relatively favourable hydrocarbon prospects, and to date, exploration work has mainly concentrated on small areas of the southwestern portion of the basin. To conduct a comprehensive study of the structural framework of the basin and the characteristics of the basement in order to identify prospective zones for hydrocarbon exploration, shipborne‐measured and satellite‐measured gravity data, shipborne‐measured magnetic data, and aeromagnetic survey data were analysed. A regularisation filtering method was used to separate and obtain regional and residual gravity and magnetic anomalies. Directional gradients of gravity and magnetic anomalies and the total horizontal gradient and vertical second derivative of gravity anomalies were employed to extract information about fault structures. Regression analysis methods were used to determine the basement depth. The geological significance of the gravity and magnetic fields was examined, the structural framework of the basin was assessed, the basement depth was estimated, and favourable hydrocarbon exploration prospects within the basin were identified. The results show that the Tobago Basin contains complex structures consisting mainly of two groups of faults trending in northeasterly and northwesterly directions and that the major northeasterly trending faults control the main structural configuration and depositional system within the basin. The basement of the Tobago Basin has deep rises and falls. It can be divided into the following four secondary tectonic units: the western sub‐basin, the central uplift area, the southern sub‐basin, and the northeastern sub‐basin. The central uplift area and northeastern sub‐basin are most likely to have developed hydrocarbon accumulations and should be targeted for further exploration.     

Characteristics of gravity anomalies and prediction of volcanosedimentary boron deposit distribution in the Xiongba area,Tibet

Volcanosedimentary boron deposits are present within Tertiary lacustrine sediments and volcanic rocks in Xiongba, Tibet. Boron deposits are characterized by low density relative to country rocks; thus, it is possible to locate them by gravity measurements. We conducted a 1:50000 high-precision gravity survey in the Xiongba area, Tibet, and obtained the Bouguer and residual gravity anomalies. We analyzed fault systems and the distribution of sedimentary and volcanic rocks and their relation to the volcanosedimentary boron deposits. The processing of the gravity data revealed local gravity variations and fault structures. We applied preferential downward continuation and wavelet transform to the gravity data, and in conjunction with geological data, we predicted the distribution of volcanosedimentary boron deposits.     

Interpretation of gravity and magnetic anomalies near Naples,Italy, using computer techniques

Some aspects of the analytical procedure for an automatic fitting of gravity and magnetic anomalies are discussed. An example of application relative to the Campania volcanic district near Naples, Italy, is reported. Gravity anomalies in that area mainly reflect the depth to the carbonate basement rocks. A small gravity high over Vesuvius may indicate buried volcanic layers. Magnetic anomalies near L. Patria and Castel Volturno are interpreted to be caused by buried intrusive rocks of higher susceptibility. Computer models of subsurface rock units with appropriate densities and susceptibilities have been generated to match the observed anomalies.     

九江-瑞昌矿集区的3D结构及对区域找矿的启示

本文分析了九瑞地区重、磁场的分布特征.运用最新的方法技术对已有的重磁资料重新进行处理:采用重磁多尺度边缘检测方法,对九瑞矿集区区域重力和航磁数据进行了边缘检测,并根据检测结果重新厘定了断裂系统的展布位置.在整理、分析九瑞地区地质、地层物性资料的基础上,对实测的1∶5万重磁数据进行较细致的准三维反演.同时将重磁三维物性反演应用到岩浆岩空间结构研究中,获得了矿集区地层结构及岩浆岩三维空间形态特征.根据反演所得磁化率强弱,分析了岩体的基性程度,为寻找与火山岩、侵入岩体有关的金属矿产提供了指示信息.最终建立的模型给出了地下地层的分布特征、控矿构造的展布规律、与成矿相关岩体的三维形态以及已知矿点的空间分布特征,为在九瑞矿集区的深部寻找隐伏矿体提供了新的信息.     

Interpretation of gravimetric and magnetic anomalies in the Cameros Basin (North Spain): combination of deep and shallow sources

Outstanding potential field anomalies (gravimetric and magnetic) in the Cameros Basin (N Spain) follow a WNW-ESE trend, parallel to the geological structures resulting from Mesozoic extension and Tertiary basin inversion. The positive Bouguer gravity anomaly (15 mGal) is interpreted as the result of a strong contrast between the density of Tertiary rocks of the foreland basin and the Paleozoic and Mesozoic rocks, combined with crustal thickening in the Iberian Chain with respect to the Ebro Basin. The dipolar magnetic anomaly, slightly shifted to the south with respect to the relative maximum of the Bouguer anomaly, can be interpreted as related to volcanic rocks within the basement, which are linked to Triassic rifting as witnessed by outcrops of basalts along the basin margins. An exhaustive analysis of rock properties (density, magnetic susceptibility and remanence) and basin geometry from other sources (seismic reflection profiles) allow to constrain variations in crustal thickness and the location of large-scale basement faults.     

给定观测精度下的点(线)源模型重力与磁力垂向识别能力研究

重、磁勘探具有效率高、成本低、工作范围广等优点,已在地球物理勘探中得到了广泛应用.前人大多在不考虑重、磁勘探观测精度的条件下进行了垂向识别能力的研究,但在考虑重、磁观测精度条件下,重力（重力异常、重力张量）与磁力（磁力异常、磁力三分量、磁力张量）对孤立异常的垂向识别能力如何则需要进行深入的理论研究.本文从重、磁场正演理论出发,以球体（点源模型）和无限延伸水平圆柱体（线源模型）为例,考虑给定观测精度条件下,以重力和磁力幅值大小与观测精度的关系来研究垂向识别能力,从而消除了背景场的影响,提高了研究结果的可靠度.通过研究表明,对于孤立异常,重力张量在浅部一定深度内比重力异常的垂向识别能力强,该深度与重力异常和重力张量观测精度的比值成正比;垂直磁化磁力张量在浅部一定深度内比化极磁力异常的垂向识别能力强,该深度与磁力异常与磁力张量观测精度的比值成正比;磁力在浅部一定深度内比重力的垂向识别能力强,该深度与地质体的磁化强度和剩余密度比值、重力观测精度和磁力观测精度比值成正比.通过重力和磁力垂向识别能力的研究将为重、磁勘探的实际应用起到指导作用.     

二连盆地“三明治”电性结构特征及其对古生界的研究意义

由于长期以来东北地区晚古生代的构造属性一直被认为是海西期褶皱带或碰撞造山带,而一直被作为中、新生代盆地的变质结晶基底成为油气勘查的禁区,限制了古生界的油气勘探和开发."佳木斯—兴蒙地块"新认识及二连盆地赛51井取得的石炭纪碳酸盐岩油藏的发现,开辟了二连盆地古生界油气勘探的新领域.本文对二连盆地3条典型的CEMP电测量资料进行了二维电阻率反演,结合研究区地震反射资料、钻孔电测井资料、钻井岩芯资料、岩石露头资料,揭示了二连盆地具有"三明治"的电性结构特征,同时指出二连盆地古生界广泛分布,在古生界与中生界之间广泛分布一套以玄武岩、安山岩、凝灰岩等为主的火山岩.这种"三明治"电性结构特征不但给古生界的研究提供了线索和依据,而且火山岩(高阻层)与其下伏古生界地层(低阻层)之间可形成良好的储盖组合(或生储盖组合),其对于古生界油气勘探的研究具有重要的意义.同时,认为二连盆地中、新生界和古生界的构造沉积演化主要经历了4个重要阶段,即:泥盆纪的南北差异性升降阶段、石炭-二叠纪的地层稳定沉积阶段、三叠纪-晚侏罗世的古生代盆地被改造和中新生代盆地构造格架形成的重要演化阶段以及中新生代的地层沉积阶段.     

贺根山镁铁-超镁铁岩与缝合带的关系——来自地球物理的证据

长期以来,对内蒙古贺根山缝合带中的镁铁-超镁铁岩,有着"蛇绿岩"、"岩浆岩"、"幔源熔-流体"等不同认识.近年来在铬铁矿中发现了金刚石等深部地幔矿物,如何认识携带这些物质的幔源熔-流体的上涌?缝合带在镁铁-超镁铁岩的形成过程中起到什么作用?解决问题的关键之一是弄清镁铁-超镁铁岩的深部产状,这需要来自地球物理观测数据的依据.本文基于航磁和重力数据的研究表明,贺根山地区的蛇绿岩块均呈现高磁异常特征,其中贺根山岩块埋深明显大于其他岩块,表现为高磁异常与低重力异常,与地表出露的蛇纹石化的镁铁-超镁铁岩带对应.对航磁化极异常与布格重力异常数据做了向上延拓处理,进行了磁源形态及底界深度的估算,并利用基于相异度算子的边缘增强方法辅助识别断裂.结果表明,贺根山岩块贯穿地壳,且附近存在超壳断裂.地表至中地壳主要由蛇纹石化的镁铁-超镁铁岩组成,下地壳主要为超镁铁质岩组成,它们充填在一组宽约30km的NEE向断裂带中,大地电磁测深剖面揭示的壳-幔电性结构进一步证实控制缝合带的是岩石圈断裂,贺根山缝合带具有明显的根部.由此推测,在地幔底辟上涌的背景下,幔源熔-流体沿着岩石圈断裂持续上升到达浅表,暗示该区板块的拼合可能是通过深部幔源物质的侵入而成的.