利用重磁资料解译川西地区深层断裂构造及预测火山岩顶界面深度.pdf

第 卷第期成都理工大学学报（自然科学版）年月 （）：文章编号 （）收稿日期 。基金项目中国石油西南油气田分公司科技重大专项（）。第一作者潘力（），男，助理工程师，研究方向：地球物理勘探，：。利用重磁资料解译川西地区深层断裂构造及预测火山岩顶界面深度潘力，何青林，陈康，梁生贤，朱家富，陈先洁，张富利，张晨，何昌兴悦（成都理工大学 地球物理学院，成都 ；中国石油西南油气田公司 勘探开发研究院，成都 ；中国地质调查局 成都地质调查中心，成都 ）摘要研究川西地区可能控制火山岩相带的深大断裂以及火山岩的分布范围、深度等基础问题。利用最新大比例尺重磁资料，对数据进行斜磁化、剩磁处理和场源分离，采用物性反演方法进行反演，利用互相关系数和深度加权函数同时进行自约束，引入迭代正则化的 法求解大型离散不适定问题以节省计算成本；采用正演模拟计算并确定了重磁资料主要解决的相关地质问题，即重力数据主要用于解译测区内具有密度差异的断裂构造，磁测数据主要用于预测火山岩；有针对性地采取响应的定性分析或三维反演成像，综合解译深层断裂构造、预测火山岩分布范围和埋深，并且根据已施工钻井对解译成果进行评价。火山岩埋深预测结果的最小相对均方误差仅 ，绝大多数情况下小于；重力资料能更好地反映深大基底断裂，而磁测资料对于火山岩的反应更为明显。关键词四川盆地；火山岩；重力；磁法；断裂构造 分类号 ；文献标志码犃 狆 狆 犾 犻 犮 犪 狋 犻 狅 狀狅 犳犵 狉 犪 狏 犻 狋 狔犪 狀 犱犿 犪 犵 狀 犲 狋 犻 犮犱 犪 狋 犪 狋 狅 狋 犺 犲 犻 狀 狋 犲 狉 狆 狉 犲 狋 犪 狋 犻 狅 狀狅 犳犱 犲 犲 狆犳 犪 狌 犾 狋 狊 狋 狉 狌 犮 狋 狌 狉 犲 狊犪 狀 犱狆 狉 犲 犱 犻 犮 狋 犻 狅 狀狅 犳狏 狅 犾 犮 犪 狀 犻 犮狉 狅 犮 犽狋 狅 狆 犻 狀 狋 犲 狉 犳 犪 犮 犲犱 犲 狆 狋 犺犻 狀狋 犺 犲狑 犲 狊 狋 犲 狉 狀犛 犻 犮 犺 狌 犪 狀犅 犪 狊 犻 狀，犆 犺 犻 狀 犪 ，犆 狅 犾 犾 犲 犵 犲 狅 犳犌 犲 狅 狆 犺 狔 狊 犻 犮 狊，犆 犺 犲 狀 犵 犱 狌犝 狀 犻 狏 犲 狉 狊 犻 狋 狔狅 犳犜 犲 犮 犺 狀 狅 犾 狅 犵 狔，犆 犺 犲 狀 犵 犱 狌 ，犆 犺 犻 狀 犪；犈 狓 狆 犾 狅 狉 犪 狋 犻 狅 狀犪 狀 犱犇 犲 狏 犲 犾 狅 狆 犿 犲 狀 狋犚 犲 狊 犲 犪 狉 犮 犺犐 狀 狊 狋 犻 狋 狌 狋 犲，犘 犲 狋 狉 狅 犆 犺 犻 狀 犪犛 狅 狌 狋 犺 狑 犲 狊 狋犗 犻 犾犪 狀 犱犌 犪 狊犉 犻 犲 犾 犱犆 狅 犿 狆 犪 狀 狔，犆 犺 犲 狀 犵 犱 狌 ，犆 犺 犻 狀 犪；犆 犺 犲 狀 犵 犱 狌犆 犲 狀 狋 犲 狉狅 犳犆 犺 犻 狀 犪犌 犲 狅 犾 狅 犵 犻 犮 犪 犾犛 狌 狉 狏 犲 狔，犆 犺 犲 狀 犵 犱 狌 ，犆 犺 犻 狀 犪犃 犫 狊 狋 狉 犪 犮 狋：，书书书 ，犓 犲 狔狑 狅 狉 犱 狊：；年川西地区 井取得重大突破，这是继 年 井以来再次在火山岩地层中获得高产工业气流，揭示了川西地区火山岩具有巨大的勘探潜力，并且各类相关研究成果呈井喷式发表。何青林等开展了火山岩地震反射特征研究，建立了火山岩不同岩相的地震反射模式；陆建林等基于 井对川西地区二叠系火山岩的形成发育特征、储层特征、天然气成因及成藏组合开展了研究；马新华等系统总结了火山岩气藏的勘探历程和 井区气藏的特征、天然气成矿的控制因素和模式；夏茂龙等、罗冰等综合钻井、重磁、地震等资料预测区域火山机构分布特征，认为简阳地区火山岩为峨眉山大火成岩省的特殊组成部分。目前的火山岩勘探手段也存在一定的难点，例如火山岩地震反射特征复杂，在剖面上难以直接识别火山岩和火山机构等。本文通过对川西地区最新的高精度重力和磁测资料进行精细处理和三维反演，探讨区内断裂构造和火山岩识别中重磁资料的作用及效果，从而为该地区火山岩天然气勘探提供地球物理证据。地质、地球物理概况 地质背景峨眉山玄武岩最早由赵亚曾先生于 年命名。峨眉山玄武岩喷发是峨眉地裂运动（东吴运动）短期强烈活动的表现，形成了巨厚的“峨眉山玄武岩”，呈一个长轴近南北向的菱形，面积达 ，是中国唯一一个被国际学术界认可的大火成岩省。其主要的喷发期是二叠纪茅口晚期龙潭期，底与下伏茅口组石灰岩不整合接触，顶与龙潭组陆相红色碎屑岩不整合接触。四川盆地位于扬子地块西北部，盆地呈北东向菱形展布，属于克拉通盆地，盆地四周被断裂围限，南西边缘以峨眉宜宾断裂带为界，北西以龙门山断裂带为界，南东以齐岳山断裂带为界，北缘以城口房县深断裂带为界。研究区位于四川盆地西部（图），面积约 ，区域地层除缺失泥盆系和石炭系以外，其余各个时代地层均有出露。区内发育的断层主要为 向的龙泉山三台巴中镇巴断裂带。峨眉山玄武岩受地形和距火山口远近影响整体有西南厚、向北东减薄的趋势 。主要地层密度、磁性特征重磁异常的特征取决于研究区岩石的密度和磁性参数，因此明确研究区的岩石密度、磁性参数特征是处理重磁数据的前提。前人研究表明，四川盆地存在个主要的密度界面，即须家河组底面和基底顶面，须家河组与下伏地层间、震旦系与基底之间的密度分别存在 、的差异。研究区内高磁化率层位为二叠系火山岩段以及基底的火山岩部分；其余的沉积岩磁化率在（）左右，而火山岩层位的磁化率普遍在（）之间，差异十分明显。数据处理 数据来源与噪声处理重力、磁法测量数据均为 年实测获得，面积 ，数据网格 ，共计 个重力和磁测点。测区内大小城镇与各种规模的工业企业较多，尤其在大型城镇与大规模工业生产场地附近，磁犜异常往往存在明显的局部极值。在对照第期潘 力，等：利用重磁资料解译川西地区深层断裂构造及预测火山岩顶界面深度图研究区位置及断裂分布图犉 犻 犵 （据马新华等）犜局部极值点、磁测原始记录和遥感图像后，将人文噪声影响所致的数据畸变点予以剔除。布格重力异常无明显的局部极值，其各类异常的等值线畸变形态均呈线状，分析认为，这种线性畸变应当为断裂构造的响应，也说明重力数据受人文干扰小，故无需做噪声处理。斜磁化、剩磁处理和场源分离由于受斜磁化和剩磁的影响，磁异常形态往往比较复杂，增加了解释的困难性。图 为假设磁化方向与地磁场方向一致时的化极结果，它与磁异常模量图（图 ）存在一定的差别，但总体上形态相近。由于后者受剩磁影响较小，因此在定性解释时，我们以磁异常模量为主；同时又考虑到后者的异常形态相对异常体，可能会存在一定的位置偏移。因此反演时的输入数据为犜，并假设磁化方向与现今地磁场方向一致，在获得扰动磁化强度后的定量（半定量）解释时，同时兼顾磁异常模量与反演结果。图测区磁异常模量与化极结果对比图犉 犻 犵 成都理工大学学报（自然科学版）第 卷 重磁数据反演重、磁数据反演就是根据重力数据或磁测数据获取地下地质体的密度或磁性。本次反演的数据为剩余重力异常和犜磁异常数据，反演前将地下空间剖分为若干个小单元，只求解各单元的扰动密度或扰动总磁化强度磁性，即物性反演方法。在反演中，我们利用互相关系数和深度加权函数同时进行自约束 ，以提升结果的分辨率、可靠性和稳定性；此外，引入迭代正则化方法 法求解大型离散不适定问题以节省计算成本 。解译方法 重磁数据解决的主要地质问题为明确测区内重磁数据特点能够反映出的主要地质问题，我们根据过井地震剖面和地层密度、磁性，设计了响应的模型（图）。正演模拟计算结果表明：在龙泉山断裂带附近，表现为明显的高重力加速度异常（重力勘探中，凡提到重力都是指重力加速度，所以下文中表述为重力异常），异常幅值比背景场大 （）以上，异常梯度陡，但对火山岩反应不够灵敏，在火山岩尖灭部位，甚至无明显的异常特征。断裂带不能引起明显的磁异常，但火山岩在平面上的变化能引起明显的磁异常特征，如在火山岩尖灭部位，磁异常幅值会陡然变低。根据上述模拟计算结果，认为重力数据主要用于解译测区内具有密度差异的断裂构造，磁测数据主要用于预测火山岩。为了进一步了解重磁数据对上述地质问题反应的灵敏度，我们分别就重力和磁测进行正演模拟计算。根据不同埋深、不同断距的基底断裂的重力正演结果（图）表明，当基底平均埋深为，基底断距大于 时，引起的重力异常幅值为 ，达 到 了 均 方 误 差（）的两倍以上，说明利用重力数据能够识别测区内断距大于 的基底断裂。图研究区重磁正演模型及模拟结果犉 犻 犵 （）重磁正演结果曲线；（）密度模型；（）磁力模型；（）对应地震剖面（据罗冰等）第期潘 力，等：利用重磁资料解译川西地区深层断裂构造及预测火山岩顶界面深度图不同埋深、不同断距的基底断裂的重力正演结果犉 犻 犵 通过不同埋深以及不同厚度磁测正演结果，并且根据结果用线性关系大致表示出火山岩厚度与对应磁异常幅值之间的关系（图），可以得到如下结论：在火山岩埋深相差 的情况下，其引起的磁异常幅值勉强满足大于倍均方误差（），说明利用磁异常判别火山岩埋深存在可能，但预测的误差可能较大。当火山厚度小于 时，引起的磁异常幅值小于倍均方误差，说明利用磁测数据能够识别的火山岩厚度必须大于。解译结果 断裂构造识别根据各类重力异常场特征共解译断裂 条，其中 向断裂条，其余均为 向；断裂构造全部分布于测区南西部。在测区北东部，无论各类重力异常还是各类磁异常图上，均无明显的断裂构造痕迹，说明在北东部断裂构造不发育。再根据剩余重力异常场不同上延高度的水平梯度模量，推测断裂规模、深度从大到小依次为：“、”“、”“、”“、”“、”“、”。其中区域断裂有：“、”“、”和“、”，并可能切穿上部沉积盖层，深达结晶基底；其余则为次级断裂（图）。火山岩识别采用自主研发的重磁反演算法 及软件，以获取测区内剩余磁化强度结构，然后定量化预测火山岩顶界面。剩余磁化强度成像结果显示，测区内 深度都无明显的磁性；在进入 深度后，高磁性异常突显；至 左右的深度，磁性最强。结合已施工钻孔情况，判断这种高磁异常就是火山岩的响应（图）。根据上述剩余磁化强度成像结果，尽管能判断出强磁性体就是火山岩的响应，但利用数字成像结果准确预测火山岩顶界面还存在一个阈值的问题，具体表现在：不同岩相的磁性强弱不同，理论上阈值不是某一个数字，而是随岩相变化的，但岩相又未知；不同的火山岩厚度、不同的岩相引起的体积效应不同，对应的阈值也会不同。可以说，火山岩顶界面预测埋深的影响因素有：岩相（岩性）、厚度以及反演结果的可靠性和准确度等，而各影响因素之间的关系又极其复杂，因此，目前还远未达到精确（准确）预测顶界面埋深的程度。本文根据实测总磁场强度、物性测试统计结果，并兼顾高磁性体的连续性，以某一阈值来大致成都理工大学学报（自然科学版）第 卷图磁测正演结果犉 犻 犵 预测火山岩顶界面埋深。由犜等值线形态与磁异常模量对比结果可知，测区内火山岩可能无强剩磁或古地磁方向与现今地磁方向差别不大，因此磁化强度与磁化率之间的转换关系式可写为 犜 犕 犜犕（）式中：为磁化率；犕为磁化强度；犜为磁感应强度。本文根据实测总磁场强度，取犜 ，根据物性测试统计结果，令磁化率 ，计算得出磁化强度犕约为。因此，本文以 为阈值预测火山岩顶界面（图）。结果评价前人研究表明四川盆地重力异常走向呈北东向的条带状分布，解译出川西基底断裂主要也为北东向展布，这与本文所识别的断裂构造特征具有高度的一致性。也就说明本文利用研究区内的重力资料识别出的深大断裂与实际基底断裂基本对应，处理结果可靠性高。同时，为判断本文中预测的火山岩顶界面深第期潘 力，等：利用重磁资料解译川西地区深层断裂构造及预测火山岩顶界面深度图重力水平模量图上的断裂分布犉 犻 犵 图三维高磁体切片图犉 犻 犵 图预测的火山岩埋深犉 犻 犵 度是否符合实际，根据搜集的研究区内钻遇火山岩测井资料统计出火山岩实际埋藏深度与预测深度的比对结果如表。其中钻遇火山岩的测井既有处于喷溢相，火山岩厚度较大的 井、井等；也有处于溢流相，火山岩厚度较小的 井等。从结果上看，火山岩埋藏深度预测最小相对均方误差仅为 ，且