电磁—地震联合研究深层火山岩储层——以辽东凹陷火山岩勘探为例

电磁—地震联合研究深层火山岩储层——以辽东凹陷火山岩勘探为例徐桂芬;何展翔;石艳玲;董卫斌【摘要】深层火山岩地层地震反射杂乱,使火山岩油藏的探测异常困难.针对火山岩储层制定了针对性的技术路线,采用电磁—地震联合模拟退火约束反演、利用钻井标定和地震构建浅层模型,提高深层的反演分辨率;通过岩性分析,揭示工区火山岩电性变化规律,提出以电性变化节律编码技术解释储层,有效识别出沙三段多期次火山岩及相带分布,精细刻画粗面岩的分布特征,为探测类似的复杂目标提供了可行的研究思路.【期刊名称】《石油地球物理勘探》【年(卷),期】2019(054)004【总页数】10页(P937-946)【关键词】电磁—地震联合约束反演;火山岩;节律编码;辽东凹陷【作者】徐桂芬;何展翔;石艳玲;董卫斌【作者单位】东方地球物理公司综合物化探处,河北涿州072751;河北工业大学电气工程学院,天津300401;南方科技大学前沿与交叉科学研究院,广东深圳518055;南方科技大学地球与空间科学系,广东深圳518055;东方地球物理公司综合物化探处,河北涿州072751;东方地球物理公司综合物化探处,河北涿州072751【正文语种】中文【中图分类】P6310 引言近年来，深层火山岩勘探屡获突破，已成为油气储量的重要增长点，火山岩勘探成为石油勘探领域的热点。

火山岩勘探实践表明，综合物探方法预测深层火山岩较为有效，可从重力、磁性、电性不同侧面获取火山岩的地球物理信息，研究火山岩地层的岩性并进行岩相划分、火山口识别、火山岩厚度及分布规律的研究，辅助地震勘探增强火山岩识别的可靠性[1-2]。

学者们在重磁电震综合研究火山岩方面进行了大量的探索。

刘云祥等[3-4]提出利用重力异常剥离、磁异常分离与弱信号处理、建场测深反演与异常提取等关键技术，从高精度重磁力资料和建场测深资料中分离、提取与火山岩有关的重磁电异常，综合解释火山岩岩性，进行岩相划分；杨辉等[5]研究松辽盆地徐家围子断陷火山岩岩性时，提出利用重力—地震剥层技术宏观预测火山岩岩性；王玉华等[6]在研究松辽盆地古龙断陷深层火山岩气藏时，采用重磁力弱信息提取、重力异常剥离、重磁震联合反演解释等技术，分离重磁异常、提取并增强与火山岩有关的重磁异常信息，实现了火山岩的综合地质解释。

徐礼贵等[7]提出了利用综合地球物理资料解释叠合盆地深层火山岩的思路：基于高精度重磁勘探和电法资料研究盆地深部结构和深大断裂，缩小并锁定深层火山岩地震勘探的有利靶区；利用地震相分析、地震属性解释、地震反演和综合评价等主要技术，实现深层火山岩从定性到定量的识别和预测，落实深层火山岩有利勘探目标。

索孝东等[8-9]提出利用重磁电震异常综合进行岩性识别的方法：根据岩石物性差异，建立不同岩性的重磁电异常和地震相特征模板，基于高精度重磁资料预测火山岩的平面分布，利用电法测深电阻率剖面和地震剖面对深层岩性垂直定深的优势，综合确定火山岩的空间位置、形态，利用重、磁、电、震多参数精细识别不同类型火山岩体。

时频电磁(TFEM)勘探方法是中国石油集团东方地球物理公司自主研发的一项电磁勘探技术[10]，它将可控源音频大地电磁(CSAMT)与瞬变电磁法(TEM)合并在一套系统中，一次采集即可获得时间域和频率域信号，改变了常规电磁法时、频域分离作业的模式，提高了采集效率和探测效果，形成了激发脉冲多样、多分量采集、多参数研究、勘探目标多样化的新型高精度大功率电磁测深方法。

TFEM法与其他电磁法相比，优势主要体现于采集方法和采集参数。

随着TFEM勘探技术的快速发展，勘探领域已从地质普查发展到深层目标勘探，例如深层火山岩及裂谷的分布等[11-12]。

辽河盆地是中国较早发现火山岩油气藏并获得大规模油气储量的地区[13]，自1996年在辽东凹陷欧利坨子地区古近系粗面岩中获高产油气流以来，探明了共计数千万吨的油气储量，其中绝大部分富集在粗面岩中[14]，表明辽东凹陷火山岩具有良好的勘探潜力。

目前针对辽东凹陷火山岩进行了较为深入的研究，主要涉及火山喷发旋回和期次划分[15]、火山岩岩相分类[16]、火山体结构研究[17-18]、粗面岩成因及喷发模式[19-20]、断裂与火山岩分布的关系[21]、成藏条件分析与成藏模式[22-25]等。

钻井揭示辽东凹陷古近系火山岩有11种类型，其中最主要的是火山碎屑岩、粗面岩和玄武岩，兼具喷出和侵入、陆上与水下堆积特征的多层系、多类型火山岩组合[19]。

由于火山岩多期次发育，总厚度大，对地震信号的屏蔽作用强，使深层火山岩地震资料成像精度低，多呈空白或杂乱反射，导致火山岩解释结果可靠性差，阻碍了深层火山岩储层的油气勘探进程。

常规重磁电联合能够解决火山岩的平面分布，但对于沙三内幕这类多期次叠置的复杂火山岩目标难以有效辅助地震进行高精度识别。

针对研究区沙三内幕多层叠置的火山岩期次划分和岩性、岩相识别等勘探难题，本文提出电磁—地震联合反演及多重约束反演解释等技术，取得良好效果，对类似复杂火山岩油气目标的研究具有一定的借鉴和指导意义。

1 电磁—地震联合探测的可行性分析辽河凹陷位于渤海湾盆地的东北部，属典型的中新生代发育的大陆裂谷盆地。

在郯庐断裂多期活动以及印度和太平洋板块斜向俯冲的叠加作用下，自中生代以来发生了多期岩浆活动，新生代火山活动尤为频繁，以辽东凹陷火山活动最为强烈，期次多、持续时间长、火山岩层分布广泛。

由于火山岩厚度大，对地震信号屏蔽作用强，深层火山岩地震信息不足，对期次划分和岩性、岩相分布不清楚。

根据已有的地质资料及地震解释成果，部署了7条TFEM测线，其中包括6条主测线和1条联络线，工区位置和测线分布见图1。

此次研究的目的是探索联合电磁、地震勘探解决深层多期叠置的复杂火山岩储层分布问题。

1.1 电性特征分析岩石物性是联系地质与地球物理之间的纽带，因此地层(岩石)的地球物理特征是研究地质问题的基础。

通过研究区24口井的测井曲线(图2)及三界泡地区的电性测试资料，综合分析沙三内幕不同岩性的电阻率特征，为TFEM资料处理和解释提供依据。

沙三段火山喷发旋回和期次研究表明：沙三段构成一个火山喷发旋回，内部可划分为5个期次，自下而上依次为：qc1、qc2、qc3、qc4、qc5。

根据钻井资料，期次界面表现为沉积夹层或岩性组合及岩相组合的突变。

根据钻井统计，沙三段主要岩性有粗面岩、玄武岩、蚀变玄武岩、凝灰岩、角砾岩和沉积岩，含少量闪长岩和辉绿岩。

针对主要岩性分别进行电测井数据统计，直方图统计结果(图3)显示：沉积岩电阻率最低，约为1～8Ω·m，主值为5Ω·m；火山岩电阻率整体偏大，不同岩性之间电阻率差异较大。

粗面岩电阻率最大，为30～90Ω·m，主值为65Ω·m；玄武岩和蚀变玄武岩电阻率为2～20Ω·m，主值为10Ω·m；火山角砾岩和凝灰岩电阻率为10～30Ω·m，主值为20Ω·m。

主要岩性的电阻率特征可归纳为：粗面岩>火山碎屑岩(角砾岩、凝灰岩)>玄武岩(包含蚀变玄武岩)>沉积岩。

上述统计结果为电法反演识别高阻粗面岩储层提供了物性基础。

这些电性特征可概括为：沉积岩为低阻，粗面岩为高阻，其余火山岩为中等，因此在横向和纵向上可形成明显的电性异常。

在沙三段内部5个期次中，火山岩喷发特征的差异形成以沉积岩为背景的5套高阻层，同时，由于岩性的变化形成了电阻率节律特征差异，可见电性分析为TFEM 资料识别高阻粗面岩储层提供了有效的电性模型。

1.2 可行性论证针对沙三内部火山岩埋藏深(普遍大于4km)、多层叠置分布，岩性、岩相复杂，区分难度大等难点，根据该区地震解释及电性特征，在地质分析基础上，建立包括3套沉积层的火山岩模型(图4a)，火山岩电阻率为70Ω·m，最小厚度为50m，埋深为2300～3300m。

先进行正演，获得模型仿真数据，再进行反演。

图4b的反演结果很好地恢复了火山岩薄层模型，分辨出三套高阻目标层，同时也清晰地反映了地层电性纵横向上的变化规律，表明电磁勘探在该区可以区分不同期次、具有一定厚度的火山岩。

图1 研究区位置示意图(左)和TFEM测线部署图(右)图2 研究区内主要井的电测井曲线图3 研究区岩性—电阻率直方图图4 基于地震—测井资料建立的火成岩模型及反演结果(a)火山岩模型；(b)约束反演结果2 电磁—地震联合勘探技术2.1 电磁—地震综合研究技术路线在TFEM资料预处理和常规处理的基础上，依据浅层可靠的地震资料建立浅层模型，反演浅层地电模型，在此基础上固定浅层，重点反演深层目标；通过钻井、电阻率剖面和地震剖面的综合对比分析，互相约束，进行火山岩期次的划分和横向对比追踪；总结不同岩性、岩相的电性异常特征和地震响应特征，对各期火山岩进行岩性、岩相综合解释。

具体研究思路和技术流程见图5。

图5 深层火山岩电磁—地震联合研究技术流程具体来说，包括以下4个步骤。

(1)如果地震资料在浅层碎屑岩中具有较好的分辨率，应首先通过地震资料解释构建浅层构造模型，通过钻井的层位、岩性和电性曲线信息，获得不同岩性地层的电阻率，一方面可以建立火山岩模型，通过正演和反演进行电法可行性分析，另一方面可以构建浅层地电结构，通过迭代反演，获得浅层地电模型。

(2)通过TFEM资料的预处理，获得综合地电参数，固定浅层几何参数和电阻率数据，重点反演深层地电结构，经过迭代反演建立深层地电结构模型。

(3)通过钻井火山岩期次界面和岩相界面的识别，地震剖面和电阻率剖面相结合，进行剖面上火山岩期次的划分和对比追踪；以钻井岩性、岩相标定做约束，建立不同岩性、岩相的地震相和电阻率剖面特征识别模板，二者互相约束、互为补充，完成各期次火山岩结构、火山岩相的识别。

(4)进行剖面综合解释，完成不同期次火山岩岩相平面分布和不同期次粗面岩储层在平面上分布特征的研究。

上述步骤中电磁—地震联合处理和电性节律编码解释技术是解决问题的关键，下面重点介绍。

2.2 电磁—地震联合反演技术研究区的浅层地震资料具有较高的分辨率，为了获得高分辨率的电阻率成像，采用基于井震建模的模拟退火反演算法进行约束反演[26]。

其基本思想是将待反演模型的每个参数看作是熔化物体的每一个分子，将目标函数看作是熔化物体的能量函数，通过缓慢减小一个模拟温度的控制参数进行迭代反演，使目标函数最终达到全局极值点。

L2范数情况下非线性的TFEM地球物理反演问题的目标函数为φ(m)=[d-A(m)]H*[d-A(m)](1)式中： d表示数据矩阵；A表示模型正演矩阵；m表示模型参数向量；上标H表示共轭转置，当观测数据为实数时，共轭转置退化为转置；“*”表示褶积运算。

将TFEM约束反演问题的每一个模型参数向量mi等效为物体的某种状态ri，将目标函数φ(m)等效为物体的能量函数Ei，引入一个随迭代次数变化的控制参数T模拟物体的温度，就可以得到TFEM的非线性反演的Metropolis接受准则P(mi→mj)(2)式中： kb是Boltzmann常数，在时频电磁反演过程中可设为1；T表示温度。