南海重磁异常特征及火成岩分布(英文)

南海火成岩油气藏具有广阔的勘探前景,综合利用地球物理方法圈划与识别火成岩体、研究火成岩分布是火成岩油气藏研究的基础。针对南海重磁场特征,采用低纬度、变倾角化极技术进行了磁异常化极处理,利用优选延拓方法实现重磁异常分离并提取南海海域浅部火成岩重磁异常信息,利用磁异常三维相关成像给出南海火成岩的三维空间等效分布,在重磁梯度突出局部异常边界信息的基础上,通过梯度加权的重磁相关分析勾画不同类型火成岩的平面展布,火成岩的分布特征显示出受地壳深部结构及断裂构造的控制与影响。     

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

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

向上延拓对航磁梯度数据的影响及其规律研究

本文首次开展了实测梯度数据的延拓转换处理方法的研究,首先利用理论公式推导证明实测的梯度数据进行延拓转换处理的可行性,并设计了单一理论模型,绘制了梯度上延数据及ΔT上延后求导的理论平面剖面曲线图,定性分析两者的一致性程度;同时设计了复杂模型体群,开展了梯度数据进行延拓转换处理后的数据解释规律的研究,梯度数据的主要目的是突出浅部异常体,其上延处理的意义则是相对于较深部异常体更能突出浅部异常体,同时实测的梯度数据进行延拓后的结果相比于ΔT异常求导后再上延的结果更接近于真实的情况等结论.同时利用实测的两个不同飞行高度的梯度数据,针对利用模型得出的规律性结论进行实例验证,并分析该规律的应用思路.     

21世纪重力与磁法勘探的展望

对21世纪重力与磁法勘探的仪器，数据处理技术，解释理论与方法，应用领域等方面的发展方向进行了分析与展望，发展航空标量，矢量，梯度重力测量和航空全梯度磁力测量，三分量磁力测量，提高综合信息采集能力；开展卫星重磁测量，综合卫星，航空，地面重磁测量资料研究地球结构与构造；发展高精度数据处理技术；重磁异常弱信号的提取，不同深度重磁异常的划分，低纬底变倾角化磁极以及位场面延拓；发展复杂条件下三维重磁场多参数综合反演可视化技术以及快速自动反演技术；探索磁性多参数的应用新领域，充分发挥磁法在环境污染调查中的作用并开拓应用新领域。     

基于重磁资料研究贵州省深部构造特征

贵州省在大地构造位置上属于羌塘—扬子—华南板块下的扬子陆块,受不同构造旋回期的发展演化控制,早期地质构造被掩盖,而研究深部构造有利于构造单元划分、成矿作用等研究.本文利用区域重磁资料,结合其他地质资料,利用优选向上延拓分离不同尺度重磁异常,并对不同尺度重磁异常进行归一化总水平导数垂向导数处理,提取线性构造信息,分析区内主要断裂构造特征,进行定性分析;根据分离出的大尺度重磁异常反演研究区莫霍面、居里面起伏特征,莫霍面起伏介于34～46 km之间,总体上具有西高东低的特征;居里面深度介于22～28.5 km之间,西部以坳陷区为主,东部以隆起区为主,以区内两条深大断裂为界,深部构造具有南北分带、东西分块的特征,反映出贵州省上扬子地块与江南造山带两个三级构造分区的深部构造特征差异.     

位场方法在非常规油气勘探中的应用

在常规油气勘探领域,高精度的重磁资料能够非常好的完成油气资源的远景评价,为地震勘探和钻井提供指导.在非常规油气勘探领域,重磁勘探方法显得尤为重要.因为非常规油气资源的地质结构、圈闭条件、储层评价、成藏机理、渗流特征等方面更加复杂和未知,地震勘探和钻井难以直接开展,亟需非地震方法的优选规划.由于重磁位场方法并不能"特异性"的揭示地下页岩分布和页岩气的资源潜力,从某种程度上说,重磁位场方法在常规油气和非常规油气的勘探技术是没有多大区分的,核心都是"综合"二字.其中,对重磁资料的处理要经过常规的延拓、剩余重力异常提取和磁异常化极等,以及重力异常剥皮、park反演、梯度反演等关键处理,能够获取特定地层的厚度与埋深和断裂信息.本文以页岩气资源的重磁位场勘探技术为中心,通过调研国内外重磁方法在非常规油气勘探阶段的应用实例,总结规律,构建利用位场方法开展页岩气等资源勘探的技术体系与流程.     

重磁对应分析技术在缅甸X区块的应用研究

缅甸X区块属于陆上低成熟勘探区块,地震采集覆盖面积较少.为了预测研究区构造特征及火山岩分布情况,研究分析该区重磁异常的同源性,利用位场数据处理中的重磁对应分析技术、最大梯度模等方法,对该区重磁资料进行了处理和认识.综合分析了研究区重力场、磁力场、重磁对应分析中相关系数的基本特征,结合研究区已有地质资料及地震解释,进行构造特征认识,并综合分析了断裂划分及圈定勘探远景区.     

大兴安岭西坡德尔布干断裂地球物理特征与构造属性

对1:250万重力异常向上延拓不同高度并分别计算135°方向水平导数,德尔布干断裂表现为北东向延伸的重力场分界线,采用基于DCT法欧拉反褶积对1:10万高精度重磁数据向上延拓0～20 km追踪断裂轨迹,并进行1:100万大地电磁测深反演和海拉尔盆地地震剖面解释,识别出德尔布干断裂具有向南东倾斜,切割深度至下地壳的特征....     

准噶尔盆地克-百断裂带火山岩分布特征的重磁资料解释

文中从各类岩石磁化率、密度统计数据出发,对比分析不同地层和岩性的磁性、密度,指出在研究区内开展火山岩重磁研究具有有利的物性基础。在利用带通滤波方法对重磁异常进行位场分离,消除沉积盖层和区域场影响后,获得由火山岩引起的剩余重磁异常,位场分离结果显示,反映火山岩的局部重磁异常多呈串珠状NE向分布。同时采用水平总梯度和垂向导数等边界信息增强手段,圈定研究区内的断裂构造分布情况,划分出了4条向SE弧形突出的NE向主干断裂及其他次级断裂共计15条。通过大地电磁测深了解主要断裂的垂向延伸情况,进一步厘定出深、浅2个断裂系统,并结合构造背景对2个层次断裂的性质、特征和对区内火山活动的作用进行了讨论。在综合分析剩余重磁异常分布特征和断裂构造相互关系的基础上,圈定出3条火山活动带。最后,利用剩余重力异常进行视密度反演,获得古生界顶部的视密度异常,以此为基础,结合钻孔岩性、岩石实测密度值和剩余重力异常分布规律,对古生界顶部的岩性分布进行细致的划分,确定了中奥陶统科克萨依组、下石炭统太勒古拉组和下二叠统佳木河组的分布范围及岩性组成。     

扩展的倾斜角和倾斜角总水平导数方法

为了突出显示重磁场中地质体分布、断裂构造特征等信息,使图像显示错落有序、层次分明、清楚直观,有利于解释人员从不同角度认识重磁场异常特征.本文通过分析倾斜角概念,根据反余切函数特性,提出了扩展的倾斜角(iTdr)和倾斜角总水平导数(iTdr_Thdr)方法,来增强重磁异常信息,突出异常形状特点,以便于识别重磁异常中地质体边界和提取地质体分布、断裂构造平面展布特征等信息.理论模型和实际资料对比分析表明该方法简单、实用、具有较好的识别效果和较高的分辨能力,能够获取丰富的地质信息,对识别地质体边界、划分大地构造单元、确定断裂带和地质构造走向等具有实际意义.     

华北地区重力场变化特征与孕震模型研究

本文采用基于贝叶斯原理的重力平差方法对华北地区2009—2017年的16期陆地流动重力观测资料进行了统一处理,结合研究区内2009年以来4次M≥4.5 地震活动,分析区域重力场动态演化特征。基于华北地区均衡重力异常、布格重力异常、沉积层重力异常和莫霍面重力异常,研究华北地区三维密度结构特征,分析区域重力场变化与深部密度结构之间的关系。结果表明:研究区内近年来M≥4.5地震均位于区域差分重力场和累积重力场“零值”等值线附近、异常变化剧烈的梯度带及梯度带转弯部位,以及深部结构中理论均衡厚度与实际地壳厚度存在差异的非均衡区和三维密度结构的高低密度过渡带内。      

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.     

大别造山带东段重力异常多尺度分解及其构造意义

文中通过多源数据融合、模型构建、数据试验、二维离散小波变换和功率谱分析等方法获取了大别造山带东段深、浅部场源布格异常及其场源似深度,并结合地壳结构、地质构造、岩石圈有效弹性厚度和地震活动等资料,讨论了地壳深、浅部的结构特征及地震活动构造背景。结果表明,低频布格异常显示大别造山带东段与华北地块间深部构造缝合带在东部应位于青山-晓天断裂前缘,在落儿岭-土地岭断裂和商城-麻城断裂之间向N偏移至梅山-龙河口断裂之下,造山带南侧与扬子地块间深部构造缝合带位于襄樊-广济断裂以北约20km,造山带东侧与扬子地块间的深部构造转换带位于郯庐断裂带之下,造山带东段腹地显著的低频布格异常低值表明对应部位的莫霍面存在明显下凹,造山带内部的布格异常高梯度带表明其深部结构不完整;高频布格异常揭示肥中断裂、六安-合肥断裂、肥西-韩摆渡断裂和郯庐断裂带等主要断裂对地壳中上部密度结构的影响明显,落儿岭-土地岭断裂对地壳中上部密度结构的影响范围向N延伸至肥西-韩摆渡断裂前缘。结合地震活动资料进一步分析认为,大别造山带东段与华北地块在青山-晓天断裂前缘附近接触和相互作用,且大别造山带东段地壳深、浅部结构均不完整,不利于应力积累,趋向于在断裂交错的脆弱部位频繁释放应力,是霍山地区小地震活动频繁的主要原因。     

中国东北地区重力场研究——利用改进的三方向小子域滤波划分主构造线及大地构造单元

利用重力异常进行大地构造单元的划分是区域重力资料处理与解释的重要任务.为了精细地划分中国东北地区大地构造单元及其分界线位置,对三方向小子域滤波进行了改进,探讨了改进的三方向小子域滤波的合理性.利用小子域及改进的三方向小子域滤波对研究区布格重力异常进行处理,并与欧拉反褶积反演结果进行对比分析.结果表明,改进的三方向小子域滤波在识别平面断裂的应用中具有更好的效果.综合分析东北地区布格重力异常特征、小子域和改进的三方向小子域滤波结果,划分出延续性较好、规模较大的断裂18条.对研究区重力异常进行了分区,描述了一、二级异常的分区特点,给出了东北地区大地构造单元的基本格架.黑河-齐齐哈尔-白城等重要构造线平面位置的确定突破了传统的重力异常分区模式,指出东北地区NE向构造具有走滑特征,且走滑距离具有由西向东逐渐增大的特点.     

基于实测重力异常和地形数据确定重力梯度的研究

本文利用澳大利亚北领地West Arnhem Land 地区实测重力异常数据并联合DEM(9")和SRTM3(3")地形高程数据,使用移去-恢复技术和Stokes积分方法计算了该地区两条剖面的重力梯度及其功率谱密度,使用FFT方法解算了整个地区的重力梯度值,结果证明了联合重力异常数据和高分辨率地形高程数据能有效地提高重力梯度的解算精度;功率谱密度的计算结果与国外成熟的重力梯度功率谱密度模型相吻合,表明高于0.3 Hz频率范围的功率谱密度可看做噪声,为重力梯度数据处理中噪声的辨别和剔除提供了借鉴,另外对重力梯度辅助导航基准图的构建以及重力梯度测量系统的标定提供了有益的探索.     

各向异性标准化方差计算重磁源边界

在重磁源边界定位方法中,传统的梯度方法易受干扰的影响使计算的边界混乱,而且在弱异常处由于叠加异常的影响很难识别场源边界.本文首先利用坐标旋转构造了各向异性高斯函数,提出了各向异性标准化方差计算重磁源边界的方法.理论分析与模型实验详细阐明了该方法的数学含义,并通过干扰分析验证了方法的稳定性与有效性,结合中扬子地区航磁异常...     

黄河上游地区地热场初步分析

运用热流、地温梯度及温泉等资料,划分了黄河上游地区的地热异常区,讨论了地热与构造及地震的关系。认为断裂构造及深部热源是地热异常的主要成因;在空间上,强震分布与热异常区的分布呈相关关系。     

Topography of the Moho and earth crust structure beneath the East Vietnam Sea from 3D inversion of gravity field data

The Moho depth, crustal thickness and fault systems of the East Vietnam Sea (EVS) are determined by 3D interpretation of satellite gravity. The Moho depth is calculated by 3D Parker inversion from residual gravity anomaly that is obtained by removing the gravity effects of seafloor and Pre-Cenozoic sediment basement topographies from the free air anomaly. The 3D inversion solution is constrained by power density spectrum of gravity anomaly and seismic data. The calculated Moho depths in the EVS vary from 30–31 km near the coast to 9 km in the Central Basin. A map of the lithosphere extension factor in the Cenozoic is constructed from Moho and Pre-Cenozoic sediment basement depths. The fault systems constructed by the maximum horizontal gradient approach include NE-SW, NW-SE, and N-S oriented faults. Based on the interpretation results, the EVS is sub-divided into five structural zones which demonstrated the different characteristics of the crustal structure.     

新疆精河6.6级地震周边地区密度构造、均衡异常以及岩石圈挠曲机理

依据EIGEN-6C4重力模型和ETOPO1高程模型数据,围绕新疆精河6.6级地震展开岩石圈均衡与挠曲机理研究,得到如下结论：（1）震中附近的布格与自由空气重力异常分别为-221和-92mGal(10~(-5 )m·s~(-2)）,震中位于重力异常高梯度带上;（2）震中周边地区地壳厚度约为50km,密度结构总体变化平缓,东西方向地壳厚度变化较小,但自南向北地壳厚度逐渐变薄,精河6.6级地震初始破裂发生在上中地壳分界面附近;（3）震中附近岩石圈承载的垂向构造应力为20MPa左右,震中位于岩石圈垂向构造应力极大值附近的高梯度带上;（4）地震周边地区岩石圈有效弹性厚度最优解为26km,加载比最优解为F_1=1,F_2=F_3=0,表明该区域岩石圈相对坚硬,且导致岩石圈变形的初始加载全部来自地表．