火山岩岩性识别方法
文献发布时间:2023-06-19 10:19:37
技术领域
本发明涉及油气勘探开发技术领域,具体而言,涉及一种火山岩岩性识别方法。
背景技术
火山岩岩性的正确识别是合理认识储层、综合评价地质特征的重要手段,对于科学开展火山岩相、火山岩储层预测研究,指导油气勘探、评价和开发部署具有重要的意义。历年来,基于常规测井刻度进行火山岩岩性识别难度大的问题是一个世界级难题。前人在火山岩岩性识别方面也做了大量的探索创新,如:
1)一种蚀变火山岩岩性识别方法及装置(专利号:CN102053268A)提出了一种利用钻、录、测井和化学分析数据关联开展火山岩岩性识别的方法和装置。其流程清晰,可操作性强,其难点在于钻、录、测井特征参数的提取及信息融合,对于经验要求较高,推广难。
2)一种油气勘探中蚀变火山岩有效储层确定方法及装置(专利号:CN103616731A)提出了一种通过岩石蚀变因子、储层品质因子,建立的蚀变火山岩有效储层的确定方法和原则,其难点是计算岩石的蚀变因子(即识别火山岩是否蚀变)。其过程较为简单,操作也相对容易,但局限于蚀变火山岩识别,且需要全直径岩心标定是否为有效储层,成本相对较高。
3)一种三维空间火山岩岩性识别方法(专利号为CN102071928A)提出了一种基于元素俘获测井、成像测井资料、岩石的结构信息,利用支持向量机的方法实现三维空间火山岩岩性识别的方法。其效果不错,但成本昂贵、推广性差。
4)一种石炭系火山岩油藏岩性识别方法(专利号为CN104929626A)提出了一种基于岩屑、测井交汇图、岩相划分、薄片定名的岩性识别方法。其方法可操作性强、结果可靠,但流程繁琐、耗时、工作量大。
5)一种基于电成像测井分形维数的火山岩岩性识别方法(专利号为CN108830140A)提出了一种基于取心、未取心井段成像测井数据并计算分维数进行岩性识别的方法。其精度相对可靠,但局限性大、成本高、对操作人员技术要求高、推广性不强。
成熟而有效的井下火山岩岩性识别方法主要依靠元素俘获能谱、化学测量(硅-碱法)、成像测井、钻井取心等手段,但是对于深层、大井距、相变快的区域,依靠上述手段的井下火山岩岩性识别方法存在成本高、推广难度大的问题,而目前为止基于常规测井资料还没有形成一套结构合理、经济高效、层次清晰、识别准确的交汇图版和识别方法。因此,在吸收前人认识、岩心露头样本标定、成像测井特征归纳和常规测井特征总结的基础上,需要提供一种低成本、高效且具有较高的推广适用性的识别火山岩岩性的方法。
发明内容
本发明的主要目的在于提供一种火山岩岩性识别方法,该方法能够利用常规测井资料来高效、准确地识别火山岩岩性,具有较高的推广适应性。
为了实现上述目的,本发明提供了一种火山岩岩性识别方法,火山岩岩性识别方法包括:步骤S10:收集目标地区的火山岩的常规测井资料,基于常规测井资料建立标准化数据库;步骤S20:从标准化数据库中,选择对火山岩岩性敏感的多条测井曲线;步骤S30:基于多条测井曲线,绘制能够识别不同火山岩岩性的测井曲线交汇图版,对火山岩岩性进行划分,确定火山碎屑岩和/或火山熔岩的分布范围。
进一步地,当火山岩岩性为火山碎屑岩时,在步骤S30之后,火山岩岩性识别方法还包括:步骤S40:从标准化数据库中,选择对火山碎屑岩岩性敏感的多条测井曲线;步骤S50:绘制能够识别不同火山碎屑岩岩性的测井曲线交汇图版,对火山碎屑岩岩性进行划分得到火山碎屑岩大类。
进一步地,在步骤S50之后,火山岩岩性识别方法还包括:步骤S60:从标准化数据库中,选择对火山碎屑岩大类的岩性敏感的多条测井曲线;步骤S70:绘制能够识别不同火山碎屑岩大类岩性的测井曲线交汇图版,对火山碎屑岩大类中的每一种岩性进行小类划分。
进一步地,能够识别不同火山碎屑岩岩性的测井曲线交汇图版为自然伽马与密度(GR-DEN)交汇图版。
进一步地,当火山岩岩性为火山熔岩时,在步骤S30之后,火山岩岩性识别方法还包括:步骤S45:从标准化数据库中,选择对火山熔岩岩性敏感的多条测井曲线;步骤S55:绘制能够识别不同火山熔岩岩性的测井曲线交汇图版,对火山熔岩岩性进行划分得到火山熔岩大类。
进一步地,在步骤S55之后,火山岩岩性识别方法还包括:步骤S65:从标准化数据库中,选择对火山熔岩大类的岩性敏感的多条测井曲线;步骤S75:绘制能够识别不同火山熔岩大类岩性的测井曲线交汇图版,对火山熔岩大类中的每一种岩性进行小类划分。
进一步地,能够识别不同火山熔岩岩性的测井曲线交汇图版为自然伽马与原状地层电阻率(GR-RT)的交汇图版。
进一步地,火山岩岩性识别方法还包括以下步骤:步骤S05:基于露头踏勘、岩心观察、实验分析和薄片鉴定方式中的一种或者多种得到的资料,为火山岩岩性进行定名,得到火山岩的岩性定名数据库;和/或,步骤S15:基于成像测井和/或元素俘获测井方式得到的火山岩岩性资料,建立火山岩的岩性资料数据库。
进一步地,火山岩岩性识别方法还包括:步骤S80:利用步骤S05中得到的火山岩的岩性定名数据库,和/或,步骤S15中建立的火山岩的岩性资料数据库,与划分出的火山岩岩性进行对比,判断火山岩岩性识别方法对火山岩岩性识别的准确性。
进一步地,能够识别不同火山岩岩性的测井曲线交汇图版为自然伽马与原状地层电阻率/声波时差(GR-RT/AC)的交汇图版。
应用本发明的技术方案,以常规测井资料为基础,通过绘制对火山岩岩性响应敏感的测井曲线的交汇图版,来识别火山岩岩性;可以解决在无取心、无特殊测井、有常规测井曲线的井段,通过建立火山岩岩性识别图版,快速对火山岩岩性进行识别,本发明的技术方案,操作步骤简单,成本低,排除了储层结构因素对火山岩岩性识别的干扰,提高了有效储层的岩性识别的准确率,提供了一种利用常规测井资料来快速、高效、低成本、准确地识别复杂火山岩岩性的方法,具有较高的推广适应性。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1示出了根据本发明的火山岩岩性识别方法的实施例的流程图;
图2示出了根据本发明的火山岩岩性识别方法的实施例的基础资料库的组成示意图;
图3示出了根据本发明的火山岩岩性识别方法的实施例的目标地区需要识别的火山岩岩性的层级分类图;
图4示出了根据本发明的火山岩岩性识别方法的实施例的火山碎屑岩-火山熔岩岩性识别的自然伽马与原状地层电阻率/声波时差(GR-RT/AC)交汇图版;
图5示出了根据本发明的火山岩岩性识别方法的实施例的火山碎屑岩岩性识别的自然伽马与密度(GR-DEN)交汇图版;
图6示出了根据本发明的火山岩岩性识别方法的实施例的火山熔岩大类识别的自然伽马与原状地层电阻率(GR-RT)的交汇图版;
图7示出了根据本发明的火山岩岩性识别方法的实施例的基性火山岩岩性识别的声波时差与补偿中子(AC-CNL)的交汇图版;
图8示出了根据本发明的火山岩岩性识别方法的实施例的中性火山岩岩性识别的自然伽马与密度(GR-DEN)交汇图版;
图9示出了根据本发明的火山岩岩性识别方法的实施例的酸性火山岩岩性识别的密度与声波时差(DEN-AC)交汇图版;以及
图10示出了根据本发明的火山岩岩性识别方法的实施例的次火山岩岩性识别的密度与补偿中子(DEN-CNL)交汇图版。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
需要指出的是,除非另有指明,本申请使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
在本发明中,在未作相反说明的情况下,使用的方位词如“上、下、顶、底”通常是针对附图所示的方向而言的,或者是针对部件本身在竖直、垂直或重力方向上而言的;同样地,为便于理解和描述,“内、外”是指相对于各部件本身的轮廓的内、外,但上述方位词并不用于限制本发明。
需要说明的是,本发明的火山岩岩性识别方法是一种基于常规测井资料进行火山岩岩性识别的方法。
需要说明的是,根据目标地区的不同,火山岩岩性存在差异;本发明的实施例选取新疆准噶尔盆地某地区作为目标地区,针对目标地区,应用本发明的火山岩岩性识别方法,其结果如图4至图10所示。
如图1至图4所示,本发明的实施例中,火山岩岩性识别方法包括:
步骤S10:收集目标地区的火山岩的常规测井资料,基于上述的常规测井资料建立标准化数据库;
步骤S20:从上述的标准化数据库中,选择对火山岩岩性敏感的多条测井曲线;
步骤S30:基于上述的多条测井曲线,绘制能够识别不同火山岩岩性的测井曲线交汇图版,对火山岩岩性进行划分,确定火山碎屑岩和/或火山熔岩的分布范围。
具体地,步骤S10为,收集、整理并标准化处理不同测井系列采集的常规测井数据,建立标准化测井资料数据库,即收集、整理目标地区的由不同测井系列采集的常规测井数据,并将上述的常规测井数据进行标准化处理,之后基于上述经过标准化处理的常规测井数据建立标准化测井资料数据库;通过对常规测井数据进行标准化处理,可以实现对由不同测井系列采集的常规测井数据进行修正,有助于消除系统误差,提高岩性识别的准确率;步骤S20为,在步骤S10的基础上,从上述的标准化测井资料数据库中,选择能够识别出火山岩岩性、并且识别效果较好的多条测井曲线;步骤S20利用步骤S10中建立的标准化测井资料数据库中的标准化的常规测井数据,确定了能够识别出火山岩岩性、且对火山岩岩性的识别效果相对较好的多条测井曲线,为保证后续步骤中的通过测井曲线准确识别火山岩岩性奠定了基础;步骤S30为,在步骤S20的基础上,绘制基于上述的多条测井曲线的交汇图版,用于识别火山岩岩性,分析绘制的测井曲线交汇图版,划分出不同火山岩岩性的分布区域,从而可以确定不同火山岩岩性的分布区域的上、下限值和/或建立不同火山岩岩性的线性关系式。
进一步地,在步骤S10中,常规测井资料为通过常规测井方法获得的测井数据,这里所说的常规测井方法主要是指在油气勘探领域所应用的9条常规测井曲线,包括自然伽马(GR)、声波时差(AC)、密度(DEN)、补偿中子(CNL)、自然电位(SP)、原状地层电阻率(RT)、侵入带电阻率(RI)、地层水电阻率(Ro)和井径(CAL)。
需要说明的是,对常规测井数据的标准化处理过程是本领域的常规技术手段,此处不再赘述。
需要说明的是,本发明的实施例中,对火山岩岩性响应敏感的测井曲线的选择是通过对火山岩岩性响应敏感的敏感参数来确定的,对上述的敏感参数的选择是在大量测井曲线及其组合的交汇图版试验基础上得来的,也就是,不断尝试做测井曲线的交汇图版,确定哪个曲线能够较准确区分火山岩岩性,能够准确区分和确定火山岩岩性的那个曲线即定为对火山岩岩性响应敏感的测井曲线。
具体地,如图4所示,本发明的实施例中,在步骤S20中,选择的对火山岩岩性响应敏感的多条测井曲线为自然伽马(GR)、原状地层电阻率(RT)和声波时差(AC)曲线,在步骤S20的基础上,步骤S30中绘制的能够识别不同火山岩岩性的测井曲线交汇图版为自然伽马与原状地层电阻率/声波时差(GR-RT/AC)的交汇图版。
如图4所示,建立自然伽马与原状地层电阻率/声波时差(GR-RT/AC)的交汇图版,其中,横坐标为自然伽马值,单位是API,纵坐标为原状地层电阻率与声波时差的比值,单位是(Ω·m)/(μs/ft);分析、划分火山碎屑岩和火山熔岩岩性分布区域;分别确定各自的分布区域的上、下限值,和/或根据需要建立线性关系式;通过上述步骤,对火山岩岩性进行大类划分,利用上述的自然伽马与原状地层电阻率/声波时差(GR-RT/AC)交汇图版将火山岩岩性区分为火山碎屑岩和火山熔岩,并划分火山碎屑岩和火山熔岩的分布区域;确定火山碎屑岩和火山熔岩的各自的分布区域的上、下限值,和/或根据需要建立火山碎屑岩和火山熔岩的各自的线性关系式。
上述技术方案,提供了一种基于常规测井资料识别火山岩岩性的方法,具体地,以常规测井资料为基础,通过绘制对火山岩岩性响应敏感的测井曲线的交汇图版,来识别火山岩岩性;可以解决在无取心、无特殊测井、有常规测井曲线的井段,通过建立火山岩岩性识别图版,快速对火山岩岩性进行识别,操作步骤简单,成本低,排除了储层结构因素对火山岩岩性识别的干扰,提高了有效储层的岩性识别的准确率,提供了一种利用常规测井资料来快速、高效、低成本、准确地识别复杂火山岩岩性的方法,具有较高的推广适应性。
需要说明的是,当确定火山岩岩性为火山碎屑岩和火山熔岩时,进一步地对火山碎屑岩岩性和火山熔岩岩性的识别不存在先后顺序,可以根据实际情况和实际需要,选择先对火山碎屑岩岩性或者火山熔岩岩性进行识别。
优选地,如图1和图2所示,本发明的实施例中,当火山岩岩性为火山碎屑岩时,在步骤S30之后,火山岩岩性识别方法还包括:
步骤S40:从标准化数据库中,选择对火山碎屑岩岩性敏感的多条测井曲线;
步骤S50:绘制能够识别不同火山碎屑岩岩性的测井曲线交汇图版,对火山碎屑岩岩性进行划分得到火山碎屑岩大类。
具体地,通过步骤S30,对火山岩岩性进行了划分,并确定了不同火山岩岩性的分布区域以及各自的分布区域的上、下限值;当火山岩岩性中含有火山碎屑岩时,通过步骤S30,确定了火山碎屑岩的分布区域的上、下限值,在步骤S30的基础上,建立测井曲线交汇图版,对火山碎屑岩岩性进行划分以得到火山碎屑岩大类;更具体地,在步骤S30的基础上,步骤S40为,从上述的标准化测井资料数据库中,选择能够识别出火山碎屑岩大类、并且识别效果较好的多条测井曲线;步骤S40利用标准化测井资料数据库中的标准化的常规测井数据,确定了能够识别出火山碎屑岩大类、且对火山碎屑岩大类的识别效果相对较好的多条测井曲线,为保证后续步骤中的通过测井曲线准确识别火山碎屑岩大类奠定了基础;步骤S50为,在步骤S40的基础上,结合通过步骤S30确定的火山碎屑岩的分布区域,绘制基于上述的多条测井曲线的交汇图版,用于识别火山碎屑岩岩性,通过分析绘制的测井曲线交汇图版,划分出不同火山碎屑岩岩性的分布区域,即划分出火山碎屑岩大类的分布区域,并进一步确定各自的分布区域的上、下限值和/或建立各自的线性关系式。
上述技术方案中,以常规测井资料为基础,结合在步骤S30中确定的火山岩各大类的分布区域,通过绘制对火山岩各大类的岩性响应敏感的测井曲线的交汇图版,来识别火山岩各大类的岩性,操作步骤简单,成本低;在步骤S10至步骤S30的基础上,本发明还提出了一种分级分类识别火山岩岩性的思路,通过绘制测井曲线交汇图版,先对火山岩岩性进行大类划分,根据划分出的火山岩各大类的分布区域,确定上述火山岩各大类的分布范围,在此基础上,针对确定的火山岩各大类,分别绘制火山岩各大类的测井曲线交汇图版,对上述各大类的岩性的分布区域进行划分,以确定上述各大类的岩性的分布范围,实现对火山岩各大类的岩性的快速识别;上述思路结构合理、层次清晰,通过以上步骤,可快速、高效、准确的对火山岩各大类的岩性进行识别,具有较高的推广适应性。
具体地,如图5所示,本发明的实施例中,在步骤S40中,选择的对火山碎屑岩岩性响应敏感的多条测井曲线为自然伽马(GR)和密度(DEN)曲线,在步骤S40的基础上,步骤S50中绘制的能够识别不同火山碎屑岩岩性的测井曲线交汇图版为自然伽马与密度(GR-DEN)的交汇图版。
如图5所示,根据步骤S30中划分出的火山碎屑岩的分布值域(即分布区域及分布范围),建立自然伽马与密度(GR-DEN)的交汇图版,其中,横坐标为自然伽马值,单位为API,纵坐标为密度值,单位为g/cm
上述技术方案,以常规测井资料为基础,结合在步骤S30中确定的火山碎屑岩的分布区域,通过绘制对火山碎屑岩岩性响应敏感的测井曲线的交汇图版,来识别火山碎屑岩岩性,上述操作步骤简单,成本低,且可对火山碎屑岩岩性进行快速、高效、准确的识别。
需要说明的是,当通过步骤S50划分出火山碎屑岩大类时,可根据划分出的火山碎屑岩大类的分布值域,判别是否需要建立新的交汇图版来进一步划分火山碎屑岩大类的更次一级的火山岩岩性。
优选地,如图1至图5所示,本发明的实施例中,在步骤S50之后,火山岩岩性识别方法还包括:
步骤S60:从标准化数据库中,选择对火山碎屑岩大类的岩性敏感的多条测井曲线;
步骤S70:绘制能够识别不同火山碎屑岩大类岩性的测井曲线交汇图版,对火山碎屑岩大类中的每一种岩性进行小类划分。
具体地,通过步骤S50,对火山碎屑岩岩性进行了划分,得到火山碎屑岩大类,并确定了火山碎屑岩各大类的分布区域以及各大类的分布区域的上、下限值;在步骤S50的基础上,步骤S60为,从上述的标准化测井资料数据库中,选择能够识别出火山碎屑岩各大类的岩性、并且识别效果较好的多条测井曲线;步骤S60利用标准化测井资料数据库中的常规测井数据,确定了能够识别出火山碎屑岩各大类的岩性、且对火山碎屑岩各大类的岩性的识别效果相对较好的多条测井曲线,为保证后续步骤中的通过测井曲线准确识别火山碎屑岩大类的岩性奠定了基础;步骤S70为,在步骤S60的基础上,结合通过步骤S50确定的火山碎屑岩各大类的分布区域,绘制基于上述的多条测井曲线的交汇图版,用于识别火山碎屑岩各大类的岩性,通过分析绘制的测井曲线交汇图版,划分出火山碎屑岩各大类的岩性的分布区域,并进一步确定分布区域的上、下限值和/或建立线性关系式。
上述技术方案中,以常规测井资料为基础,结合在步骤S50中确定的火山碎屑岩大类的分布区域,通过绘制对火山碎屑岩各大类的岩性响应敏感的测井曲线的交汇图版,来识别火山碎屑岩各大类的岩性,操作步骤简单,成本低;在步骤S40与步骤S50的基础上,本发明利用分级分类识别火山岩岩性的思路,通过绘制测井曲线交汇图版,快速对火山碎屑岩各大类的岩性的分布区域进行划分,以确定火山碎屑岩各大类的岩性的分布范围,实现对火山碎屑岩各大类的岩性的快速识别。通过上述步骤,可快速、高效、准确的对火山碎屑岩各大类的岩性进行识别,具有较高的推广适应性。
需要说明的是,对火山碎屑岩各大类的岩性的识别不存在先后顺序。
优选地,如图1至图4和图6所示,本发明的实施例中,当火山岩岩性为火山熔岩时,在步骤S30之后,火山岩岩性识别方法还包括:
步骤S45:从标准化数据库中,选择对火山熔岩岩性敏感的多条测井曲线;
步骤S55:绘制能够识别不同火山熔岩岩性的测井曲线交汇图版,对火山熔岩岩性进行划分得到火山熔岩大类。
具体地,通过步骤S30,对火山岩岩性进行了大类划分,并确定了不同火山岩岩性的分布区域以及分布区域的上、下限值;当火山岩岩性中含有火山熔岩时,通过步骤S30,确定了火山熔岩的分布区域的上、下限值,在步骤S30的基础上,建立测井曲线交汇图版,对火山熔岩的岩性进行大类划分以得到火山熔岩大类;更具体地,在步骤S30的基础上,步骤S45为,从上述的标准化测井资料数据库中,选择能够识别出火山熔岩大类、并且识别效果较好的多条测井曲线;步骤S45利用标准化测井资料数据库中的标准化的常规测井数据,确定了能够识别出火山熔岩大类、且对火山熔岩大类的识别效果相对较好的多条测井曲线,为保证后续步骤中的通过测井曲线准确识别火山熔岩大类奠定了基础;步骤S55为,在步骤S45的基础上,结合通过步骤S30确定的火山熔岩的分布区域,绘制基于上述的多条测井曲线的交汇图版,用于识别火山熔岩的岩性,分析绘制的测井曲线交汇图版,划分出不同火山熔岩岩性的分布区域,即划分出火山熔岩大类的分布区域以及分布区域的上、下限值和/或建立各自的线性关系式。
上述技术方案中,以常规测井资料为基础,结合在步骤S30中确定的火山岩的各大类的分布区域,通过绘制对火山岩的各大类的岩性响应敏感的测井曲线的交汇图版,来识别火山岩的各大类的岩性,操作步骤简单,成本低;在步骤S10至步骤S30的基础上,本发明还提出了一种分级分类识别火山岩岩性的思路,通过绘制测井曲线交汇图版,先对火山岩岩性进行大类划分,根据划分出的火山岩各大类的分布区域,确定上述火山岩各大类的分布范围,在此基础上,针对确定的火山岩各大类,分别绘制火山岩各大类的测井曲线交汇图版,对上述各大类的岩性的分布区域进行划分,以确定上述各大类的岩性的分布范围,实现对火山岩各大类的岩性的快速识别;上述思路结构合理、层次清晰,通过以上步骤,可快速、高效、准确的对火山岩各大类的岩性进行识别,具有较高的推广适应性。
具体地,如图6所示,本发明的实施例中,在步骤S45中,选择的对火山熔岩岩性响应敏感的多条测井曲线为自然伽马(GR)与原状地层电阻率(RT),在步骤S45的基础上,步骤S55中绘制的能够识别不同火山熔岩岩性的测井曲线交汇图版为自然伽马与原状地层电阻率(GR-RT)的交汇图版。
如图6所示,根据步骤S30中划分出的火山熔岩的分布值域,建立自然伽马与原状地层电阻率(GR-RT)的交汇图版,其中,横坐标为自然伽马值,单位是API,纵坐标为原状地层电阻率值,单位是Ω·m;分析、划分火山熔岩大类分布的区域(局部由于存在基质的过渡带而有重叠);分别确定火山熔岩大类的分布区域的上、下限值,和/或根据需要建立线性关系式;通过上述步骤,对火山熔岩岩性进行大类划分,利用上述的自然伽马与原状地层电阻率(GR-RT)的交汇图版,将火山熔岩岩性区分为基性火山岩、中性火山岩、酸性火山岩和次火山岩,并划分基性火山岩、中性火山岩、酸性火山岩和次火山岩的分布区域;确定基性火山岩、中性火山岩、酸性火山岩和次火山岩的各自的分布区域的上、下限值,和/或根据需要建立基性火山岩、中性火山岩、酸性火山岩和次火山岩的各自的线性关系式。
上述技术方案,以常规测井资料为基础,结合在步骤S30中确定的火山熔岩的分布区域,通过绘制对火山熔岩岩性响应敏感的测井曲线的交汇图版,来识别火山熔岩岩性,上述操作步骤简单,成本低,且可对火山熔岩岩性进行快速、高效、准确的识别。
需要说明的是,当通过步骤S55划分出火山熔岩大类(即基性火山岩、中性火山岩、酸性火山岩和次火山岩)时,可根据划分出的火山熔岩大类的分布值域,判别是否需要建立新的交汇图版来进一步划分基性火山岩、中性火山岩、酸性火山岩、次火山岩等的更次一级的火山岩岩性。
优选地,如图1至图4、图6至图10所示,本发明的实施例中,在步骤S55之后,火山岩岩性识别方法还包括:
步骤S65:从标准化数据库中,选择对火山熔岩大类的岩性敏感的多条测井曲线;
步骤S75:绘制能够识别不同火山熔岩大类岩性的测井曲线交汇图版,对火山熔岩大类中的每一种岩性进行小类划分。
具体地,通过步骤S55,对火山熔岩岩性进行了划分,得到火山熔岩大类(即基性火山岩、中性火山岩、酸性火山岩和次火山岩),并确定了火山熔岩各大类的分布区域以及各大类的分布区域的上、下限值;在步骤S55的基础上,步骤S65为,从标准化测井资料数据库中,选择能够识别出火山熔岩各大类的岩性、并且识别效果较好的多条测井曲线;步骤S65利用标准化测井资料数据库中的标准化的常规测井数据,确定了能够识别出火山熔岩大类的岩性、且对火山熔岩大类的岩性的识别效果相对较好的多条测井曲线,为保证后续步骤中的通过测井曲线准确识别火山熔岩大类的岩性奠定了基础;步骤S75为,在步骤S65的基础上,结合通过步骤S55确定的火山熔岩各大类的分布区域,绘制基于上述的多条测井曲线的交汇图版,用于识别火山熔岩各大类的岩性,通过分别绘制的测井曲线交汇图版,划分出火山熔岩各大类的岩性的分布区域,并进一步确定分布区域的上、下限值和/或建立线性关系式。
上述技术方案中,以常规测井资料为基础,结合在步骤S55中确定的火山熔岩大类的分布区域,通过绘制对火山熔岩各大类的岩性响应敏感的测井曲线的交汇图版,来识别火山熔岩各大类的岩性,操作步骤简单、成本低;在步骤S45与步骤S55的基础上,本发明利用分级分类识别火山岩岩性的思路,通过绘制测井曲线交汇图版,快速对火山熔岩各大类的岩性的分布区域进行划分,以确定火山熔岩各大类的岩性的分布范围,实现对火山熔岩各大类的岩性的快速识别。通过上述步骤,可快速、高效、准确地对火山熔岩各大类的岩性进行识别,具有较高的推广适应性。
需要说明的是,对火山熔岩各大类的岩性的识别不存在先后顺序。
具体地,当确定火山熔岩岩性(即确定火山熔岩大类为基性火山岩、中性火山岩、酸性火山岩和次火山岩)时,进一步对火山熔岩各大类的岩性的识别不存在先后顺序,可以根据实际情况和实际需要,选择先对基性火山岩或中性火山岩或酸性火山岩或次火山岩的岩性进行识别。
具体地,如图7所示,本发明的实施例中,在步骤S65中,选择的对基性火山岩岩性响应敏感的多条测井曲线为声波时差(AC)和补偿中子(CNL),在步骤S65的基础上,步骤S75中绘制的能够识别不同基性火山岩岩性的测井曲线交汇图版为声波时差与补偿中子(AC-CNL)交汇图版。
根据划分出的基性火山岩的分布值域,判断对基性火山岩进行进一步识别,建立声波时差与补偿中子(AC-CNL)交汇图版,其中,横坐标为声波时差值,单位为μs/ft,纵坐标为补偿中子值,补偿中子值为百分数;分析、划分基性火山岩岩性分布的区域;分别确定各种基性火山岩各自的分布的区域上、下限值,和/或根据需要建立线性关系式;通过上述步骤,对基性火山岩岩性进行小类划分,利用上述的声波时差与补偿中子(AC-CNL)交汇图版,将基性火山岩岩性区分为玄武岩、杏仁状玄武岩和玄武质火山角砾岩,并划分玄武岩、杏仁状玄武岩和玄武质火山角砾岩的分布区域,确定各自的分布区域的上、下限值,和/或根据需要建立各自的线性关系式。
上述技术方案,以常规测井资料为基础,结合在步骤S55中确定的基性火山岩的分布区域,通过绘制对基性火山岩岩性响应敏感的测井曲线的交汇图版,来识别基性火山岩岩性,上述操作步骤简单、成本低,且可对基性火山岩岩性进行快速、高效、准确的识别。
具体地,如图8所示,在本发明的实施例中,在步骤S65中,选择的对中性火山岩岩性响应敏感的多条测井曲线为自然伽马(GR)和密度(DEN),在步骤S65的基础上,步骤S75中绘制的能够识别不同中性火山岩岩性的测井曲线交汇图版为自然伽马与密度(GR-DEN)交汇图版。
根据划分出的中性火山岩的分布值域,判断对中性火山岩进行进一步识别,建立自然伽马与密度(GR-DEN)交汇图版,其中,横坐标为自然伽马值,单位为API,纵坐标为密度值,单位为g/cm
上述技术方案,以常规测井资料为基础,结合在步骤S55中确定的中性火山岩的分布区域,通过绘制对中性火山岩岩性响应敏感的测井曲线的交汇图版,来识别中性火山岩岩性,上述操作步骤简单、成本低,且可对中性火山岩岩性进行快速、高效、准确的识别。
具体地,如图9所示,本发明的实施例中,在步骤S65中,选择的对酸性火山岩岩性响应敏感的多条测井曲线为密度(DEN)和声波时差(AC),在步骤S65的基础上,步骤S75中绘制的能够识别不同酸性火山岩岩性的测井曲线交汇图版为密度与声波时差(DEN-AC)交汇图版。
根据划分出的酸性火山岩的分布值域,判断对酸性火山岩进行进一步识别,建立密度与声波时差(DEN-AC)交汇图版,其中,横坐标为密度值,单位为g/cm
上述技术方案,以常规测井资料为基础,结合在步骤S55中确定的酸性火山岩的分布区域,通过绘制对酸性火山岩岩性响应敏感的测井曲线的交汇图版,来识别酸性火山岩岩性,上述操作步骤简单、成本低,且可对酸性火山岩岩性进行快速、高效、准确的识别。
具体地,如图10所示,本发明的实施例中,在步骤S65中,选择的对次火山岩岩性响应敏感的多条测井曲线为补偿中子(CNL)和密度(DEN),在步骤S65的基础上,步骤S75中绘制的能够识别不同次火山岩岩性的测井曲线交汇图版为补偿中子与密度(CNL-DEN)交汇图版。
根据划分出的次火山岩的分布值域,判别对次火山岩进行进一步识别,建立补偿中子与密度(CNL-DEN)交汇图版,其中,横坐标为补偿中子值,补偿中子值为百分数,纵坐标为密度,单位为g/cm
上述技术方案,以常规测井资料为基础,结合在步骤S55中确定的次火山岩的分布区域,通过绘制对次火山岩岩性响应敏感的测井曲线的交汇图版,来识别次火山岩岩性,上述操作步骤简单、成本低,且可对次火山岩岩性进行快速、高效、准确的识别。
如图3、图4和图7至图9所示,本发明的实施例提供的基于常规测井资料识别火山岩岩性的火山岩岩性识别方法,对于基性火山岩、中性火山岩和酸性火山岩全序列岩性类别的识别效果较好,还能有效区分含基、中和酸性基质的火山角砾岩,排除了火山岩因含不同基质而产生的复杂定名的影响,提高了有效储层的岩性识别准确率。
本发明涉及一种基于常规测井刻度的火山岩岩性识别方法,利用常规测井资料建立一套快速、低成本、高效识别火山岩岩性的方法,对含基、中、酸性基质的火山角砾岩具有良好识别效果。
优选地,如图1和图2所示,本发明的实施例中,火山岩岩性识别方法还包括以下步骤:
步骤S05:基于露头踏勘、岩心观察、实验分析和薄片鉴定方式中的一种或多种得到的资料,为火山岩岩性进行定名,得到火山岩的岩性定名数据库;和
步骤S15:基于成像测井和/或元素俘获等特殊测井方法得到的火山岩岩性资料,建立火山岩的岩性资料数据库。
具体地,步骤S05为,收集、整理露头踏勘、岩心观察、实验分析、薄片鉴定后资料,建立基于以上基础研究信息(或者资料)的区域性(二级构造或者盆地范畴)标准火山岩岩性定名数据库。
即,步骤S05为,收集、整理目标地区的通过露头踏勘、岩心观察、实验分析和薄片鉴定方式中的一种或多种方式得到的资料(或者基础研究信息),之后基于上述的资料,为目标地区的火山岩岩性进行定名,建立标准化的火山岩的岩性定名数据库。
进一步地,在步骤S05中,对火山岩岩性的定名或者对火山岩岩性类别的划分要符合国家权威部门(如地矿部门、中石油行业协会等)发布的国家标准、行业标准。这样,根据国家权威部门发布的国家标准、行业标准,可以对目标地区的火山岩岩性进行统一定名,排除火山岩因含不同基质、不同内部结构特征而产生的复杂定名的影响,提高岩性识别准确率。
进一步地,本发明的实施例中,通过步骤S05,建立了基于上述的基础研究信息的区域性(二级构造或者盆地范畴)标准化的火山岩的岩性定名数据库。即,在国家权威部门发布的国家标准、行业标准中,选择部分区域(如二级构造或者盆地范畴)内的火山岩岩性的标准定名,作为目标地区的火山岩岩性定名的标准,这样,可以缩小对目标地区的火山岩的岩性进行定名的标准定名范围,减少工作量,提高工作效率,避免对目标地区的火山岩岩性定名的不规范。
具体地,步骤S15为,收集、整理特殊测井(主要是成像测井)拾取并识别出的火山岩岩性资料,建立基于特殊测井识别的火山岩信息库,即收集、整理目标地区的通过特殊测井方法(主要是成像测井)拾取并识别的火山岩岩性资料,建立基于特殊测井方法识别的火山岩的岩性资料数据库(即火山岩信息库)。
通过建立上述的火山岩的岩性定名数据库和火山岩的岩性资料数据库,可以实现,针对研究工区(即目标地区)的地质特征,确定出需要识别的岩性类别,也就是说,可以根据目标地区通过露头踏勘、岩心观察、实验分析和薄片鉴定方式及成像测井、元素俘获等特殊测井方法得到的火山岩岩性资料,确定目标地区的火山岩岩性类别,进一步确定出目标地区需要识别的火山岩岩性类别;并可根据国家权威部门发布的国家标准、行业标准,对目标地区需要识别的火山岩岩性类别进行统一定名,排除火山岩因含不同基质、不同内部结构特征而产生的复杂定名的影响,提高岩性识别准确率。
进一步地,如图2所示,本发明的实施例中,通过步骤S10、步骤S05和步骤S15,建立标准化的测井资料数据库、火山岩的岩性定名数据库和火山岩的岩性资料数据库,并建立由上述三个数据库形成的统一数据库(即图2中的基础资料库)。
进一步地,本发明的实施例中,通过上述三个数据库对比,优选对火山岩岩性响应敏感的测井项目(即测井曲线);具体地,针对目标地区同一井深的火山岩岩性,可以分别通过火山岩的岩性定名数据库和火山岩的岩性资料数据库进行识别,并同时利用标准化的测井资料数据库中测井曲线进行识别,将通过上述三个数据库分别识别的火山岩岩性进行对比,当利用标准化的测井资料数据库中测井曲线识别的火山岩岩性,与通过火山岩的岩性定名数据库和火山岩的岩性资料数据库识别的火山岩岩性相一致或相近时,上述标准化的测井资料数据库中测井曲线即为优选的对火山岩岩性响应敏感的测井曲线。
上述设置中,通过火山岩的岩性定名数据库和火山岩的岩性资料数据库能够确定出目标地区需要识别的火山岩岩性类别,并且对需要识别的火山岩岩性类别进行了符合国家权威部门发布的国家标准、行业标准的统一定名,可排除火山岩因含不同基质、不同内部结构特征而产生的复杂定名的影响;通过标准化的测井资料数据库,利用分级分类识别火山岩岩性的思路,建立火山岩岩性交汇图版,快速对火山岩岩性进行识别,排除了储层结构因素对火山岩岩性识别的干扰,提高了有效储层的岩性识别的准确率。
当然,在附图未示出的替代实施例中,还可以根据实际情况和实际需要,使本发明的火山岩岩性识别方法还包括以下步骤:
步骤S05:基于露头踏勘、岩心观察、实验分析和薄片鉴定方式中的一种或多种得到的资料,为火山岩岩性进行定名,得到火山岩的岩性定名数据库;或,
步骤S15:基于成像测井和/或元素俘获等特殊测井方法得到的火山岩岩性资料,建立火山岩的岩性资料数据库。
优选地,如图1和图2所示,本发明的实施例中,火山岩岩性识别方法还包括:
步骤S80:利用步骤S05中得到的火山岩的岩性定名数据库,和步骤S15中建立的火山岩的岩性资料数据库,与划分出的火山岩岩性进行对比,判断火山岩岩性识别方法对火山岩岩性识别的准确性。
具体地,如图2至图10所示,本发明的实施例中,图3为根据图2中基础资料库中的岩性定名数据库和岩性资料数据库中的资料确定的目标地区的火山岩岩性类别,即目标地区需要识别的火山岩岩性类别;图4至图10为根据图2中基础资料库中的测井资料数据库中的资料,通过本发明的火山岩岩性识别方法识别出的目标地区的火山岩岩性类别。
上述技术方案中,通过将图4至图10中识别出的目标地区的火山岩岩性类别与图3中已经确定的目标地区需要识别的火山岩岩性类别进行对比,可以获得本发明的火山岩岩性识别方法对火山岩岩性识别的准确性。通过本发明的实施例,可以得知,本发明的火山岩岩性识别方法对火山岩岩性的识别具有较高的准确性。该识别方法可以应用于其他需要进行岩性识别的地区。
进一步地,可以推断出,上述的标准化的测井资料数据库、火山岩的岩性定名数据库和火山岩的岩性资料数据库的建立以及基础资料库的建立,是提高本发明的火山岩岩性识别方法的识别精度适应性的必要条件之一,对于没有上述火山岩的岩性定名数据库和/或火山岩的岩性资料数据库的勘探新区本发明的火山岩岩性识别方法依然适用,本发明的火山岩岩性识别方法具有较高的推广适应性。
当然,在附图未示出的替代实施例中,根据实际情况和实际需要,本发明的火山岩岩性识别方法还包括:
步骤S80:利用步骤S05中得到的火山岩的岩性定名数据库,或,步骤S15中建立的火山岩的岩性资料数据库,与划分出的火山岩岩性进行对比,判断火山岩岩性识别方法对火山岩岩性识别的准确性。
下面,结合图1至图10详细说明本发明的操作方式。
本发明的火山岩岩性识别方法的设计思路如下:
1、针对研究工区的地质特征,确定出需要识别的岩性类别;
2、基于常规测井资料,分级分类建立各大类岩性交汇图版;
3、在大类岩性识别的基础上,建立各大类岩性的次一级岩性交汇图版。
进一步地,9条常规测井曲线的数据可能采集自多个测井系列并受不同环境的干扰,在应用本方法前需要对上述数据进行系统误差矫正,对上述数据进行系统误差矫正是保证本方法规范的必要条件之一;需要说明的是,系统误差矫正的技术要点及处理过程为本领域常规技术手段,此处不再赘述。
其中,大类岩性划分主要包括:火山碎屑岩类的凝灰质泥岩、凝灰质砂岩和火山熔岩类的基性、中性、酸性、次火山岩(侵入岩性)等。其中,岩性大类的划分如图3所示。
进一步地,次一级岩性交汇图版包括了火山碎屑岩、基性火山岩、中性火山岩、酸性火山岩、次火山岩(侵入岩性)类交汇图版,如图4~图10所示。需要说明的是,各次一级岩性交汇图版所确定的岩性的划分界限值因地域不同以及火山喷发的时间、期次、规模、持续时长和强度的差异而不同,各类火山岩性质有差异。
本发明的火山岩岩性识别方法具有以下有益效果:
1、本发明首先对各类岩性进行大类的划分,简化了岩性类别,解决了岩性类别复杂的问题;
2、本发明利用层级分类递进(即分级分类)的思路,将岩性先分大类进行交汇图版识别,再逐级做小类交汇图版识别,可识别基性、中性、酸性、次火山岩等火山熔岩及其小类岩性,包括对含基、中、酸性基质的火山角砾岩等复杂岩性的区分;
3、本发明所建立的岩性识别图版可快速、低成本、高效地识别出火山岩岩性,解决了岩性复杂、单一图版识别精度差等问题。
根据上述的设计思路,本发明的实施例的具体操作方式如下:
步骤S05:基于露头踏勘、岩心观察、实验分析和薄片鉴定方式中的一种或多种得到的资料,为火山岩岩性进行定名,得到火山岩的岩性定名数据库;
步骤S15:基于成像测井和/或元素俘获等特殊测井方法得到的火山岩岩性资料,建立火山岩的岩性资料数据库;
如图3所示,通过建立的火山岩的岩性定名数据库和火山岩的岩性资料数据库,可确定出目标地区需要识别的火山岩岩性类别,并且,火山岩岩性类别的划分要符合国家权威部门发布的国家标准、行业标准,根据上述的国家标准、行业标准,对目标地区需要识别的火山岩岩性类别进行统一定名。
步骤S10:收集目标地区的火山岩的常规测井资料,基于上述的常规测井资料建立标准化的测井资料数据库;
通过步骤S05、步骤S15和步骤S10,建立火山岩的岩性定名数据库、火山岩的岩性资料数据库和标准化的测井资料数据库,并建立由上述三个数据库形成的统一数据库(即图2中的基础资料库);
通过对上述的火山岩的岩性定名数据库、火山岩的岩性资料数据库和标准化的测井资料数据库进行对比,优选对火山岩岩性响应敏感的测井曲线。
步骤S20:从上述的标准化的测井资料数据库中,选择对火山岩岩性敏感的多条测井曲线;
步骤S30:基于上述的多条测井曲线,绘制能够识别不同火山岩岩性的测井曲线交汇图版,对火山岩岩性进行划分,确定火山碎屑岩和/或火山熔岩的分布范围;
如图4所示,建立自然伽马与原状地层电阻率/声波时差(GR-RT/AC)的交汇图版,对火山岩岩性进行划分,将火山岩岩性区分为火山碎屑岩和火山熔岩,并划分火山碎屑岩和火山熔岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
(需要说明的是,以下对火山碎屑岩和火山熔岩的岩性的识别不存在先后顺序)
步骤S40:从标准化数据库中,选择对火山碎屑岩岩性敏感的多条测井曲线;
步骤S50:绘制能够识别不同火山碎屑岩岩性的测井曲线交汇图版,对火山碎屑岩岩性进行划分得到火山碎屑岩大类;
如图5所示,建立自然伽马与密度(GR-DEN)的交汇图版,对火山碎屑岩岩性进行划分,将火山碎屑岩岩性区分为凝灰质泥岩和凝灰质砂岩,并划分凝灰质泥岩和凝灰质砂岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
步骤S45:从标准化数据库中,选择对火山熔岩岩性敏感的多条测井曲线;
步骤S55:绘制能够识别不同火山熔岩岩性的测井曲线交汇图版,对火山熔岩岩性进行划分得到火山熔岩大类。
如图6所示,建立自然伽马与原状地层电阻率(GR-RT)的交汇图版,对火山熔岩岩性进行划分,将火山熔岩岩性区分为基性火山岩、中性火山岩、酸性火山岩和次火山岩,并划分基性火山岩、中性火山岩、酸性火山岩和次火山岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
(需要说明的是,以下对火山熔岩各大类的岩性的识别不存在先后顺序)
步骤S65:从标准化数据库中,选择对火山熔岩大类的岩性敏感的多条测井曲线;
步骤S75:绘制能够识别不同火山熔岩大类岩性的测井曲线交汇图版,对火山熔岩大类中的每一种岩性进行小类划分;
如图7所示,建立声波时差与补偿中子(AC-CNL)交汇图版,对基性火山岩岩性进行小类划分,将基性火山岩岩性区分为玄武岩、杏仁状玄武岩和玄武质火山角砾岩,并划分玄武岩、杏仁状玄武岩和玄武质火山角砾岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
如图8所示,建立自然伽马与密度(GR-DEN)交汇图版,对中性火山岩岩性进行小类划分,将中性火山岩岩性区分为安山岩、熔结凝灰岩和安山质火山角砾岩,并划分安山岩、熔结凝灰岩和安山质火山角砾岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
如图9所示,建立密度与声波时差(DEN-AC)交汇图版,对酸性火山岩岩性进行小类划分,将酸性火山岩岩性区分为英安岩、流纹岩、碎裂流纹岩、酸性凝灰岩和流纹质火山角砾岩,并划分英安岩、流纹岩、碎裂流纹岩、酸性凝灰岩和流纹质火山角砾岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
如图10所示,建立密度与补偿中子(DEN-CNL)交汇图版,对次火山岩岩性进行小类划分,将次火山岩岩性区分为正长斑岩、二长玢岩、闪长玢岩和碎裂正长斑岩,并划分正长斑岩、二长玢岩、闪长玢岩和碎裂正长斑岩的分布区域,确定各自的分布区域的上、下限值,确定各自的分布值域。
步骤S80:利用步骤S05中得到的火山岩的岩性定名数据库,和步骤S15中建立的火山岩的岩性资料数据库,与划分出的火山岩岩性进行对比,判断火山岩岩性识别方法对火山岩岩性识别的准确性。
如图1至图10所示,将图4至图10中通过本发明的火山岩岩性识别方法识别出的目标地区的火山岩岩性类别与图3中根据图2中基础资料库确定的目标地区需要识别的火山岩岩性类别进行对比,可以得知,本发明的基于常规测井资料识别火山岩岩性的火山岩岩性识别方法,对火山岩岩性的识别具有较高的准确性,进一步地,可以推断出,对于没有上述的火山岩的岩性定名数据库和/或火山岩的岩性资料数据库的勘探新区,本发明的火山岩岩性识别方法依然适用,本发明的火山岩岩性识别方法具有较高的推广适应性。
从以上的描述中,可以看出,本发明上述的实施例实现了如下技术效果:本发明依据常规测井测量原理,通过分析对比已经确定的标准火山熔岩、火山碎屑岩测井响应特征及其差异,借助于已经完成定名的区域岩心和特殊测井的岩性刻度,优选对火山岩岩性响应敏感的测井曲线制作交汇图版,通过交汇图版首先区分火山碎屑岩类和火山熔岩类,其次细分火山碎屑岩类,第三步划分火山熔岩大类(基、中、酸、次火山岩),最后按照火山熔岩大类(基、中、酸、次火山岩)分步逐级细分出各大类中的诸多小类,特别是区分出含基、中、酸性基质的火山角砾岩。本发明以常规测井资料为基础,可以解决在无取心、无特殊测井、有常规测井曲线的井段,利用分级分类的思路建立火山岩岩性识别图版,快速对火山岩大类中的主要小类进行识别,排除火山岩因含不同基质、不同内部结构特征而产生的复杂定名的影响,以大类识别、区分酸碱性为主导,排除储层结构因素的干扰,提供了一种能够快速、高效、准确地识别火山岩岩性的火山岩岩性识别方法,同时能够有效区分含基、中、酸性基质的火山角砾岩,提高有效储层的岩性识别准确率,简化操作,降低成本,增加推广适应性,解决了基于常规测井进行火山岩岩性识别难度大的世界级难题。本发明的岩性识别准确率较高,操作简便,适于任何有测井曲线的录、测、钻井现场和科研工作中使用,能够用于辅助快速识别岩性,提高有效储层岩性识别精度,增加油气勘探开发效益。
显然,上述所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、工作、器件、组件和/或它们的组合。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本申请的实施方式能够以除了在这里图示或描述的那些以外的顺序实施。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
- 火山岩岩性识别方法
- 一种火山岩岩性识别方法