利用重力、磁力、电磁、地震资料综合识别火成岩方法
【技术领域】
[0001]本发明涉石油勘探中地球物理资料的处理及地质综合解释技术领域，具体涉及一种利用重力、磁力、电磁、地震资料综合识别火成岩方法。
【背景技术】
[0002]大型叠合含油气盆地历经多期构造运动与演化，早期裂谷或断陷阶段一般发生过多期强烈的岩浆活动，保存了众多的深层火成岩体。越来越多的油气探井揭示，受特殊的成岩机制和岩石物理特性影响，有些深层火成岩具有优越的储集性能，在合适的石油地质条件下可能富集石油或天然气体，形成大规模的油气藏。
[0003]就目前国内火成岩油气藏勘探技术来看，勘探早期，火成岩岩体的识别技术主要应用重、磁，目前常规应用二维和三维地震解释和反演技术。
[0004]在火成岩岩性识别上，主要应用测井技术，近年来提出了不少新方法。火成岩岩相划分和内幕解剖技术上，主要应用三维地震解释和反演技术。火成岩储层预测技术主要应用测井分析和地震解释与反演。火成岩油气藏的预测技术上，利用AVO技术已获得成功，但是由于对火成岩储层预测地质模型认识的不足，火成岩油气藏的地震预测技术精度还很不够。
[0005]但是，在某些火成岩发育区，厚层火成岩对地震信号有较强的屏蔽和吸收作用，使得地震资料火成岩内幕多呈空白或杂乱反射，制约了石炭系火成岩的描述与目标优选。
[0006]重磁电勘探虽然精度和分辨率低于地震，但在火成岩识别方面具有特定的优势条件。火成岩具有高磁化率、高电阻、基性岩密度较大、酸性岩密度较小的物性特征，可产生一定规模的重磁电异常，是非地震技术最有利的勘探目标，因此，加强非地震、特别是磁力勘探在识别火成岩方面的优势，结合地震、地质、钻井等其它资料，对火成岩平面相带分布、岩性描述及内幕进行反演，是解决石炭系火成岩描述的有效途径，也是目前火成岩油气藏勘探的技术发展趋势。
【发明内容】
[0007]本发明的目的在于提供一种利用重力、磁力、电磁、地震多种资料综合进行火成岩识别的方法；该方法针对钻井资料少，地震资料空白或杂乱反射，无法进行火成岩有效识别的地区，发挥多种资料在火成岩识别方面的优势，并且通过技术手段将各种资料有机结合，能够有效识别火成岩发育区及岩性、岩相展布特征，提高了识别精度及可靠性，降低勘探风险。
[0008]为了达成上述目的，本发明采用了如下技术方案，利用重力、磁力、电磁、地震资料综合识别火成岩方法其步骤包括:
1)利用支持向量机(SVM)物性分类技术，建立不同类型火成岩与地层物性模型；
2)火成岩的重、磁、电、震正演模拟，明确不同岩性、岩相火成岩异常响应特征及规律；
3)火成岩平面重磁异常提取及弱异常增强处理，解释划分火成岩平面分布区带及范围；
4)电法资料反演及背景分离，明确剖面火成岩电性分布特征；
5)地震资料火成岩目标处理及属性分析，明确不同火成岩相地震属性特征；
6)重磁电震资料剖面联合反演，建立研究区火成岩纵向分布框架；
7)利用剖面联合反演成果及钻井资料作为约束条件，开展重磁异常火成岩物性(视密度、视磁化率)三维反演；
8)依据物性(密度、磁化率)三维分布结合火成岩平面、剖面重磁电震特征，综合识别火成岩空间分布及岩性、岩相展布规律。
[0009]进一步，建立不同类型火成岩与地层物性模型的具体步骤是:
1.1)野外采样和钻井岩心通过测试分析获取不同岩石的密度、磁化率、电阻率、速度数据；
1.2)利用钻测井资料得到地下不同岩石的电阻率、速度数据；
1.3)利用支持向量机(SVM)物性分类技术。对上述数据进行特征分析，统计出不同岩石及地层的物性模型；
进一步，利用重磁异常解释划分火成岩平面分布区带及范围的具体步骤是:
3.1)火成岩平面磁力异常利用化极处理消除玮度对磁力异常影响；
3.2)重磁数据利用水平方向导数极值位置及水平总梯度极大值，综合划分断裂位置及走向；
3.3)火成岩平面重力异常利用正演剥层法消除火成岩以外的影响；
3.4)利用倾斜角处理增强火成岩重磁弱异常，根据增强处理等值线O值位置确定火成岩平面分布范围；
3.5)火成岩重磁异常对应分析，结合断裂分布及不同火成岩性、岩相重磁异常响应特征，定性划分火成岩岩性及岩相平面分布。
[0010]进一步，电法资料反演及解释的步骤是:
4.1)参考磁力异常对电法资料进行飞点及静态位移校正等精细预处理；
4.2)进行TE、TM高精度二维定量反演；
4.3)反演结果进行背景异常分离，突出火成岩电性异常；
4.4)分离后异常进行垂向导数处理，提高电性异常火成岩分辨率；
4.5)依据电性分布特征及火成岩电性特征，确定剖面上火成岩分布；
进一步，地震资料火成岩目标处理及属性分析的步骤是:
5.1)地震资料火成岩目标处理包含剩余振幅补偿、高精度速度分析、多项式拟合弱信号加强处理；
5.2)对处理后的地震资料进行瞬时频率、瞬时相位、振幅、相干、倾角属性分析，确定不同类型火成岩地震相及地震属性特征；
5.3)地震、电法综合解释，确定剖面火成岩分布；
进一步，重磁电震资料剖面联合反演的步骤是:
6.1)剖面联合反演初始模型由地震、电法资料综合构建；
6.2)模型填充密度、磁化率由权利要求1步骤I)确定；
6.3)在模型相对简单情况下，进行2.5维反演； 6.4)反演过程中构造修改需参考地震、电法剖面资料；
6.5)若误差小于设定值，则接受当前火成岩及地层构造模型并输出，否则以人机交互方式修改物性或构造参数，再次进行重磁正演拟合，如此重复，直到找到合理的火成岩及地层构造模为止。
[0011]6.6)构建火成岩分布框架需要交叉剖面，剖面交叉处需吻合；
进一步，重磁异常火成岩物性三维反演的是视密度、视磁化率分布；
进一步，根据地层物性特征将反演的视密度、视磁化率转换为实际岩石与地层的密度、磁化率；
进一步中，综合识别的火成岩空间分布及岩性、岩相展布规律可根据钻井资料进行修改。
[0012]本发明优点在于实现了重力、磁力、电磁、地震多种资料综合火成岩识别，较单一资料火成岩识别可靠性更高，可广泛应用于火成岩发育区火成岩解释。
【附图说明】
[0013]图1为本发明利用重力、磁力、电磁、地震资料综合识别火成岩方法的流程框图； 图2为支持向量机(SVM)物性分类技术流程框图；
图3为电磁资料火成岩识别流程图；
图4为为本发明的一具体实施例中某地区电磁资料结合重磁资料火成岩岩相识别图； 图5为本发明的一具体实施例中某地区地震资料结合重磁资料火成岩联合反演解释图；
图6为本发明的一具体实施例中某地区石炭系火成岩岩相、岩性分布预测图。
【具体实施方式】
[0014]有关本发明的详细说明及技术内容，配合【附图说明】如下，然而附图仅提供参考与说明之用，并非用来对本发明加以限制。
[0015]如图1所示，利用重力、磁力、电磁、地震资料综合识别火成岩方法包括下列步骤: I)建立不同火成岩与地层物性模型。其方法是:根据野外采样、钻井岩心测试分析资料以及测井资料等，利用支持向量机物性(密度、磁化率、电阻率、速度)分类技术，建立不同类型火成岩、地层物性模型及多物性火成岩综合识别模板。
[0016]2)利用正演模拟技术，根据步骤I)得到的火成岩及地层物性模型，构建地质模型进行重力、磁力、电法、地震正演模拟，分析所得异常其响应规律，作为下步重磁电震资料处理解释的基础。
[0017]3)火成岩平面重磁异常提取及弱异常增强处理，解释划分火成岩平面分布范围及区带具体步骤是:
步骤1:通过实测重磁异常的剩余分析、水平方向导数处理、水平总梯度处理，结合重磁异常断裂识别特征，预测断裂系统展布；
步骤2:通过延拓、滤波处理及正演剥层处理，提取由火成岩引起的重磁异常。
[0018]步骤3:通过倾斜角处理、垂向导数处理对提取后的重磁异常进行弱异常增强，结合增强处理后的O等值线位置，定性识别火成岩平面分布边界； 步骤4:利用重磁异常对应分析结果，结合已识别断裂展布、火成岩正演异常特征对火成岩平面展布范围及岩性特征进行定性描述。
[0019]4)电法资料反演及背景分离，明确剖面火成岩电性分布特征的具体步骤是:
步骤1:参考磁力异常对电法资料进行飞点及静态位移校正预处理；
步骤2:对预处理后的电法资料进行TE、TM高精度二维定量反演；
步骤3:反演结果进行背景异常分离，突出火成岩电性异常；
步骤4:分离后异常进行垂向导数处理，提高电性异常火成岩分辨率；
步骤5:依据电性分布特征及火成岩电性特征，确定剖面上火成岩分布；
5)地震资料火成岩目标处理及属性分析