本发明专利技术是利用重磁电异常的组合特征识别火成岩岩性的地球物理勘探技术方法。步骤是：采集记录岩芯物理特性；用ｋ－均值聚类对岩石的密度、磁化率、电阻率进行聚类分析，获得岩芯聚类编码；对重磁勘探结果进行重磁异常剥离获得异常，采用同样的聚类进行三维空间异常编码；将岩石物性聚类编码与重磁电异常编码相互结合反映探区每一处火成岩的岩性。本发明专利技术将重磁电综合识别火成岩岩性的方法从定性推测发展到定量统计，减少了人为的经验的主观的因素，也不会出现不同地质解释人员会出现不同的结果，提高了火成岩的识别效果。
【技术实现步骤摘要】
本专利技术涉及地球物理勘探技术，具体是一种利用重磁电异常的组合特征 识别火成岩岩性的方法。
技术介绍
目前，在油气勘探中，火山岩相关的油气藏相继被大量发现，与火成岩 相关的油气藏勘探和研究己引起了地球物理界的重视。但是在复杂的火成岩 发育区，单一的利用地震勘探识别火成岩岩性、岩相存在一定难度，主要是 由于火成岩特殊的物性和内部结构，如果深度比较浅，地震勘探能得到火成 岩顶部形态，但难以确定火成岩岩体展布及其岩性；如果火成岩比较深，上 覆地层存在高速层，则对火成岩的有效反射信息被大大消弱，甚至难以发现 火成岩目标，使地震方法难以有效识别火成岩。火成岩具有强磁性，利用高 精度磁力勘探能够发现火成岩，但难以确定具体深度，浅部火成岩与深部火 成岩需要利用其它信息来加以分离，但火成岩也有磁性强磁性弱的，还是难 以确定和有效识别；重力勘探也可以用于识别火成岩，火成岩与其它岩石具 有不同密度，而且不同火成岩之间密度也有差异，但是也需要利用其它资料 进行重力异常的分离，获得不同目标体的布格异常，从而研究目标体的密度， 但也不能明确分辨火成岩及岩性；电磁法具有一定测深能力，可以直接探测 到火成岩体，因为，火成岩与围岩的电阻率有时存在很大差别，也有差别很 小的，要识别火成岩也有一定难度。因此，人们采用综合的方法，把一个探 区的重磁电震资料进行联合处理、综合解释，利用地震可以作为重磁浅层模4型，剥离浅层异常，从而获得反映深部的重磁异常，再把反映目标的重磁电 震异常进行综合对比，定性的推断火成岩范围及火成岩岩性岩相。但是，目 前只能依靠人为的经验推断，存在人的主观因素，不同地质解释人员会出现 不同的结果，同一解释者两次解释也会产生相左的结果，使火成岩的识别效 果难以令人信服，进而影响到钻探验证成功率。
技术实现思路
本专利技术目的是提供一种综合利用重、磁、电、震等地球物理勘探技术和 岩石密度、磁化率、电阻率等性质准确无误圈定和识别火成岩的利用重磁电 异常的组合特征识别火成岩岩性的方法。本专利技术提供以下实现步骤1) 在探区钻井采集岩芯，测量并记录每块岩芯标本样本的物理特性； 步骤1所述的物理特性是密度、磁化率、电阻率。2) 采用k-均值聚类对岩石的密度d、磁化率m、电阻率r，数据进行聚 类分析；步骤2所述的聚类分析是先给出分类数，选出聚类种子数据核心，根据与聚类种子数据核心距离把所有点分类，再求均值作为新的种子，重新按照距离分类，如此反复，直到各类最后没有变化了为止。步骤2所述的聚类分析将距离最短的类别合并为一类。步骤2所述的聚类分析分别对密度、磁化率和电阻率岩芯测量数据进行聚类分析，分别获得对每块岩石标本三种物性的聚类，每种物性选择的聚类数要求一致。3) 对聚类结果按升序或降序进行编码，获得探区所有岩芯的聚类编码；步骤3所述的聚类编码以岩石标本密度、磁化率、电阻率的顺序，每一 块标本有一个编码，并用0和1的二进数表示。4) 采用常规的方法对探区进行重磁勘探，对勘探结果进行重磁异常剥离 处理获得反映探区目标的重力、磁力的异常分布；对电磁法采用二维反演获 得电阻率二维断面；步骤4所述的剥离处理是利用浅层地震和钻井建立密度模型，计算出浅 层重力响应，并从实测异常中剥离Gs=GfGq，式中G。为实测异常，Gq浅 层模型异常，对深层区域异常采用重力向上延拓方法获得Gy，将Gy从Gs中 减去，就获得较探区深层探区目标的重力异常Gm=Gs-Gy=G。-Gq-Gy。步骤4所述的剥离处理对磁力异常Mo直接向上延拓获得区域异常My，将 My从Mo中减去Ms=Mo_My，再对Ms求取垂直二次导数，即可得到反映探区 深层目标的磁力异常Mm二Ms"。步骤4所述的电磁法二维反演处理，采用广义逆二维反演获得分层电阻率 异常，直接获得反映深探区层目标的电阻率。5) 对由重磁电异常处理获得的反映探区空间的重力、磁力和电阻率异常， 采用步骤2)和3)同样的聚类处理得到三维空间异常编码；6) 模式识别解释将步骤3)得到的岩石物性聚类编码与步骤5)得到 的探区重磁电异常编码相互对应，就可一一对应反映探区每一处火成岩的岩 性。本专利技术将重磁电综合识别火成岩岩性的方法从定性推测发展到定量统 计，减少了人为的经验的主观的因素，也不会出现不同地质解释人员会出现 不同的结果，提高了火成岩的识别效果。本专利技术虽然是针对火山岩的，但完全适用所有岩石岩性的识别。 附图说明图1为本专利技术实施例古生界地层岩性编码及岩性解释图。 具体实施例方式采用本专利技术实施例是对某地区火成岩主要分布于古生界，中新生界基本 无火成岩，识别古生界岩石的物性特征。1、 采集该探区钻井岩芯，测量并记录每块岩芯标本样本的物理特性我们收集了该区55 口钻井的古生界岩芯1650块，对岩芯进行编号il i2 i3……in, n=1650，然后对1650块岩芯进行物性测量，即测量岩芯的密度、 磁化率和电阻率的测量工作，获得该区密度、磁化率和电阻率的岩芯物性数据岩芯密度dl d2 d3……dn;岩芯磁化率ml m2 m3 ......mn;电阻率rl r2r3 ......rn。2、 对岩石的密度d、磁化率m、电阻率r，数据进行k-均值聚类分析 对获得的岩石的三种物性数据，密度、磁化率、电阻率分别进行k-均值聚类分析。首先，我们对密度数据采用k-均值聚类，给出分类数为3类， 并挑选出3个密度数据点，d3,d601，dl23作为"聚类种子"，即初始聚类核 心；根据和这3个点的距离远近，把所有密度数据分成3类。再把这3类的 中心(均值)作为新的种子，重新按照距离分类。如此反复下去，直到各类 最后没有变化了为止，即完成了密度数据的聚类，把岩芯密度分为高密、中 密和低密三个等级。同样，我们也把岩芯磁化率分为强磁性、中磁性和弱磁 性三个等级；把岩芯电阻率分为高电阻率、中电阻率和低电阻率三个等级。3、对聚类结果按升序进行编码，获得一个探区所有岩芯的聚类编码编 码方法以密度、磁化率、电阻率的顺序，每一块标本有一个三种物性的组合 编码，对每一种物性均用0和1的二进数表示，这样每种岩石的三种物性的 编码位数为6位，前两位数代表密度，中间两位数是磁化率，后两位是电阻 率。下表列出了该探区岩石聚类编码<table>table see original document page 8</column></row><table>4、重磁电异常的分离与聚类处理对探区重力异常进行剥离浅层用地 震钻井建模，并计算浅层模型异常Gq，区域异常用向上延拓的方法，可直接 获得Gy，然后从实测异常中减去浅层和深部区域异常，即获得研究目标的重 力异常Gm-Gm二 G。-Gq-Gy。对探区磁力异常，采用延拓方法获得区域异常My，并从实测异常中减去 My获得剩余磁力异常Ms，然后求取剩余磁力异常的二次导数，即获得研究目 标的磁力异常Mn^Ms";对电磁法采用广义逆二维反演，获得从浅到深的电阻率断面，根据电阻 率异常特征我们可以获得研究目标古生界地层的电阻率异常分布。5、同样，对上面获得的古生界地层的重力、磁力异常和电阻率异常的平面数据进行k-均值聚类分析处理，同样选择聚类数为3，获得研究目标层的 平面重磁电异常编码，见附本文档来自技高网...
【技术保护点】
一种利用重磁电异常的组合特征识别火成岩岩性的方法，其特征在于：采用以下实现步骤： １）采用常规的方法在探区钻井采集岩芯，测量并记录每块岩芯标本样本的物理特性； ２）采用ｋ－均值聚类对岩石的密度ｄ、磁化率ｍ、电阻率ｒ，数据进行聚类分析； ３）对聚类结果按升序或降序进行编码，获得探区所有岩芯的聚类编码； ４）采用常规的方法对探区进行重磁勘探，对勘探结果进行重磁异常剥离处理获得反映探区目标的重力、磁力的异常分布，对电磁法采用二维反演获得电阻率二维断面；５）对由重磁电异常处理获得的反映探区空间的重力、磁力和电阻率异常，采用步骤２）和３）同样的聚类处理得到三维空间异常编码； ６）模式识别解释：将步骤３）得到的岩石物性聚类编码与探区重磁电异常编码相互对应，就可一一对应反映探区每一处火成岩的岩性。
【技术特征摘要】
【专利技术属性】
技术研发人员：何展翔，江汶波，王晓帆，刘云祥，杨战军，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：13[中国|河北]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。