CN103485758A - 一种利用多种资料拟合反演复杂岩性解释方法及装置 - Google Patents
一种利用多种资料拟合反演复杂岩性解释方法及装置 Download PDFInfo
- Publication number
- CN103485758A CN103485758A CN201310349346.8A CN201310349346A CN103485758A CN 103485758 A CN103485758 A CN 103485758A CN 201310349346 A CN201310349346 A CN 201310349346A CN 103485758 A CN103485758 A CN 103485758A
- Authority
- CN
- China
- Prior art keywords
- response
- characteristic value
- logging
- response characteristic
- log
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000004044 response Effects 0.000 claims abstract description 198
- 239000011435 rock Substances 0.000 claims abstract description 86
- 239000012530 fluid Substances 0.000 claims abstract description 61
- 238000011160 research Methods 0.000 claims abstract description 56
- 238000001228 spectrum Methods 0.000 claims abstract description 45
- 238000004458 analytical method Methods 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims description 67
- 239000003921 oil Substances 0.000 claims description 37
- 239000000470 constituent Substances 0.000 claims description 31
- 230000000704 physical effect Effects 0.000 claims description 22
- 230000005251 gamma ray Effects 0.000 claims description 18
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 12
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 claims description 8
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 8
- 239000004411 aluminium Substances 0.000 claims description 8
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 8
- 229910052782 aluminium Inorganic materials 0.000 claims description 8
- 229910052791 calcium Inorganic materials 0.000 claims description 8
- 239000011575 calcium Substances 0.000 claims description 8
- 229910052742 iron Inorganic materials 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 229910052710 silicon Inorganic materials 0.000 claims description 8
- 239000010703 silicon Substances 0.000 claims description 8
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 7
- 239000005864 Sulphur Substances 0.000 claims description 7
- 238000009415 formwork Methods 0.000 claims description 7
- 238000012549 training Methods 0.000 claims description 5
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000011161 development Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 238000012797 qualification Methods 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- -1 shale Substances 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用多种资料拟合反演复杂岩性解释方法及装置,所述方法包括:根据研究工区的地质情况,将储层划分为流体组分、泥质组分及岩石骨架组分;给定流体组分响应特征值及泥质组分响应特征值;根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;利用测井仪器,获取研究工区地层的常规测井响应值及元素含量;联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;利用测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出研究工区的地层中的流体组分、泥质组分及岩石骨架组分的最优相对体积,绘制岩性分析结果。
Description
技术领域
本发明是针对复杂岩性储层的分析,尤指一种利用多种资料拟合反演复杂岩性解释方法及装置。
背景技术
测井技术作为一项定性、定量评价储层的主要方法,贯穿于油气勘探开发各个阶段,在碎屑岩储层评价中发挥了重要作用,随着勘探对象的日益复杂,测井技术越来越受重视,但还存在很多难题和挑战。火山岩储层评价是火山岩油气藏描述的重点,而火山岩岩性解释又是火山岩储层评价的基础。岩性解释不准确,将直接导致漏失油气层。因此准确进行火山岩岩性解释是开展火山岩储层评价的基础和关键。
前期的研究者大多单一的利用常规测井资料、地震资料等,采用参数交会图法、成像测井识别法、神经网络法、对应分析识别法、偶极声波测井、岩石强度参数识别法、地震法等方法,识别结果为简单的岩性剖面,主要停留在定性程度,不能为精细的勘探开发提供准确的岩性剖面资料,特别是在火山岩的油气勘探中,由于岩性对测井响应的影响往往超过储层流体的影响,因此岩性成分的精细解释将是勘探成功与否的关键。元素俘获能谱测井(Elemental Capture Spectroscopy,简称ECS测井)能够测量组成地层岩石的各种元素信息,提供了一种从化学组分的角度解释火山岩岩性的可能性,但由于氧化物闭合模型不易确定,使得该项技术难以快速掌握、推广应用较难。另一方面,在以往研究中,元素俘获能谱测井资料的应用之前要做很多复杂的工作,例如氧化物闭合模型、原始值转为相对含量等,而这些计算不容易掌握,地区不同,模型参数差异也很大,未能与常规测井资料形成一种易操作、有效的火山岩岩性定量解释方法。
发明内容
本发明的目的是利用常规测井资料,元素俘获能谱测井等测井方法,基于地层组分,提供了一种简单、有效的火山岩岩性定量方法,以提高复杂储层岩性解释的符合率,以此作为储层综合评价的基础和关键。
为达到上述目的,本发明提供了一种利用多种资料拟合反演复杂岩性解释方法,所述方法包括:根据研究工区的地质情况,将储层划分为流体组分、泥质组分以及岩石骨架组分;给定所述流体组分的流体组分响应特征值及所述泥质组分的泥质组分响应特征值;根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定所述研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;利用测井仪器,获取所述研究工区地层的常规测井响应值及元素含量;根据所述流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;利用所述测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出所述研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果。
为达到上述目的,本发明还提供了利用多种资料拟合反演复杂岩性解释装置,所述装置包括:储层划分模块,用于根据研究工区的地质情况,将储层划分为流体组分、泥质组分以及岩石骨架组分;流体组分及泥质组分响应特征值获取模块,用于获取所述流体组分的流体组分响应特征值及所述泥质组分的泥质组分响应特征值,其中,所述流体组分响应特征值以及所述泥质组分响应特征值通过给定直接获得;岩石骨架组分响应特征值获取模块,根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定所述研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;响应值获取模块,用于利用测井仪器,获取所述研究工区地层的常规测井响应值及元素含量;测井响应方程组建立模块,用于根据所述流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;岩性定量解释模块,用于利用所述测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出所述研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果。
本发明实施例公开的利用多种资料拟合反演复杂岩性解释方法及装置具备以下显著特点:1、不仅立足于现有成熟的常规测井系列,而且综合应用了较难推广应用的元素俘获能谱测井资料,方法原理简单,技术上容易实现;2、提出了元素俘获能谱测井资料应用的新途径,应用前景广阔;3、精确计算出了火成岩岩性成分相对含量,为岩石命名提供量化的可能;4、实现了火山岩、碳酸盐岩等复杂岩性识别精度从定性到定量化的发展,极大提高了岩性解释的符合率,在一定程度上可用来修正、校正录井剖面记录上的人为误差。利用该发明处理了吐哈油田、大庆油田资料齐全的15口井,处理结果与312块薄片鉴定结果对比,符合率达到93%,较以前提高了10%左右。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1为本发明一实施例的岩性定量解释方法流程图。
图2为本发明一实施例的岩性定量解释装置的结构示意图。
图3A为本发明一具体实施例的玄武岩物性分析孔隙度与密度曲线交会图。
图3B为本发明一具体实施例的玄武岩物性分析孔隙度与中子曲线交会图。
图3C为本发明一具体实施例的玄武岩物性分析孔隙度与声波曲线交会图。
图4为本发明一具体实施例的多种资料拟合反演岩性定量解释成果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
图1为本发明一实施例的岩性定量解释方法流程图。如图1所示,本方法包括:
步骤S101,根据研究工区的地质情况,将储层划分为流体组分、泥质组分以及岩石骨架组分;
步骤S102,给定流体组分的流体组分响应特征值及泥质组分的泥质组分响应特征值;
步骤S103,根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;
步骤S104,利用测井仪器,获取研究工区地层的常规测井响应值及元素含量;
步骤S105,根据流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;其中,元素俘获能谱测井曲线是利用元素俘获能谱测井获得的曲线。
步骤S106,利用测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果。
在步骤S101中,一般认为储层由多种岩石骨架、泥质、石油、水(或者还有天然气,但天然气成分不是必须的)组成,根据研究工区的不同,主要体现在了岩石骨架的多样性方面,对于火山岩、碳酸盐等复杂储层的岩石骨架就包括了多种岩性成分。
在本实施例步骤S101、S102中,流体组分至少包括了石油以及水,石油及水对应的流体组分相应特征值通过直接给定获得。
在本实施例中,岩石骨架组分的岩性至少包括:玄武岩、安山岩、凝灰岩、辉绿岩及砂岩。
上述步骤S103的方法更详细的说,还包括:
根据测井资料,利用油田现场取心、地质描述、物性分析和实验室薄片鉴定方法,确定岩石骨架组分的岩性;
利用声波测井、密度测井及中子测井,结合物性分析与实际测井曲线,以数据拟合方式来确定研究工区的各岩石骨架组分的声波测井响应特征值、密度测井响应特征值及中子测井响应特征值;
利用自然伽马测井及元素俘获能谱测井,结合求取实际测井曲线平均值的方法确定岩石骨架组分中各个岩性的自然伽马测井响应特征值及元素俘获能谱测井响应特征值。
其中,测井资料包括了成像测井资料、常规测井资料;声波测井、密度测井及中子测井,自然伽马测井资料都属于常规测井资料;元素俘获能谱测井资料能够测量组成地层岩石的各种元素信息,提供了一种从化学组分的角度解释火山岩岩性的方法,这种测井资料属于非常规的测井资料。
在本实施例的步骤S104中,测井仪器包括:密度测井仪器、中子测井仪器、声波测井仪器及自然伽马测井仪器,元素俘获能谱测井仪器;其中
利用密度测井仪器、中子测井仪器、声波测井仪器及自然伽马测井仪器获取研究工区地层的密度响应值、中子响应值、声波响应值、自然伽马测井响应值;
还利用元素俘获能谱测井仪器,获取研究工区地层中的铝、钙、铁、硫、硅的元素含量。
在本实施例的步骤S105中,建立的测井响应方程组为:
其中,n为研究工区中的地层组分的个数;如一具体实施例中,地层组分包括了玄武岩、安山岩、凝灰岩、辉绿岩、砂岩、泥质、水、石油八个组分。(前5部分是岩石骨架组分、泥质是泥质组分、水和石油是流体组分)
xi为研究工区中的地层组分的相对体积;
Aij为研究工区中的流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值;
m为测试仪器测量的常规测井响应值及元素含量的个数;
Bj为测试仪器测量的常规测井响应值或元素含量。
在本实施例的步骤S105中,测井响应方程组中元素俘获能谱测井曲线的应用是直接用于建立测井响应方程组,而不需要转换为各元素的相对含量。
在本实施例的步骤S106中,先用Householder变换求解一个最小二乘解,然后将每个最小二乘解作为共轭梯度优化计算的初始值进行优化迭代计算求解;以得出研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果,然后再绘制岩性分析结果。
图2为本发明一实施例的岩性定量解释装置的结构示意图。如图2所示,装置包括:
储层划分模块10,用于根据研究工区的地质情况,将储层划分为流体组分、泥质组分以及岩石骨架组分;
流体组分及泥质组分响应特征值获取模块20,用于获取流体组分的流体组分响应特征值及泥质组分的泥质组分响应特征值,其中,流体组分响应特征值以及泥质组分响应特征值通过给定直接获得;
岩石骨架组分响应特征值获取模块30,用于根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;
响应值获取模块40,用于利用测井仪器,获取研究工区地层的常规测井响应值及元素含量;
测井响应方程组建立模块50,用于根据流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;
岩性定量解释模块60,用于利用测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果。
在本实施例中,流体组分至少包括石油及水,石油及水对应的流体组分响应特征值通过直接给定获得。岩石骨架组分至少包括:玄武岩、安山岩、凝灰岩、辉绿岩及砂岩。
在本实施例中,岩石骨架组分响应特征值获取模块30还用于根据测井资料,利用油田现场取心、地质描述、物性分析和实验室薄片鉴定方法,确定岩石骨架组分的岩性;
利用声波测井、密度测井及中子测井,结合物性分析与实际测井曲线,以数据拟合方式来确定研究工区的各岩石骨架组分的声波测井响应特征值、密度测井响应特征值及中子测井响应特征值;
利用自然伽马测井及元素俘获能谱测井,结合求取实际测井曲线平均值的方法确定岩石骨架组分中各个岩性的自然伽马测井响应特征值及元素俘获能谱测井响应特征值。
在本实施例中,响应值获取模块40还用于利用密度测井仪器、中子测井仪器、声波测井仪器及自然伽马测井仪器获取研究工区地层的密度响应值、中子响应值、声波响应值、自然伽马测井响应值;
还利用元素俘获能谱测井仪器,获取研究工区地层中的铝、钙、铁、硫、硅的元素含量。
在本实施例中,测井响应方程组建立模块50建立的测井响应方程组包括:
其中,n为研究工区中的地层组分的个数;(前5部分是岩石骨架组分、泥质是泥质组分、水和石油是流体组分)
xi为研究工区中的地层组分的相对体积;
Aij为研究工区中的流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值;
m为测试仪器测量的常规测井响应值及元素含量的个数;
Bj为测试仪器测量的常规测井响应值或元素含量。
在本实施例中,测井响应方程组建立模块50建立的测井响应方程组中,元素俘获能谱测井曲线的应用是直接用于建立测井响应方程组,而不需要转换为各元素的相对含量。
在本实施例中,岩性定量解释模块60还包括:先用Householder变换求解一个最小二乘解,然后将每个最小二乘解作为共轭梯度优化计算的初始值进行优化迭代计算求解;以得出研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果,然后再绘制岩性分析结果。
本发明提出的利用多种资料拟合反演复杂岩性解释方法及装置包括的技术主要有三种:地层组分的测井响应特征值提取、常规测井与元素俘获能谱测井综合利用、最优化求解测井响应方程组。已有的资料中没有提及受孔隙度影响较大的测井项目的地层组分响应特征值的提取技术,未见常规测井与元素俘获能谱测井联合建立方程组来定量计算地层组分。
结合图1及图2,下面以一具体实施例来本发明的方法及装置。
结合步骤S101,S102,S103,对吐哈油田某15口井(有元素俘获能谱测井资料)进行储层的初步划分,包括了岩石骨架组分,流体组分,泥质组分;利用油田现场取心、录井数据、薄片鉴定资料,对吐哈油田15口有元素俘获能谱测井资料的井进行分析,并对各种组分岩性出现的频率进行统计,结果发现出现频率最高的典型岩性大致有玄武岩、安山岩、凝灰岩、辉绿岩、砂岩、泥质、石油、水八个组分。
在具体实施例中,需要获取玄武岩、安山岩、凝灰岩、辉绿岩、砂岩、泥质、石油、水八个组分的声波测井、密度测井、中子测井、自然伽马测井及元素俘获能谱测井的铝、钙、铁、硫、硅元素含量的响应特征值。
其中,石油、水等流体的性质地区差异较小,应用理论值完全可以满足要求,因此流体组分的各测井响应特征值可以直接有技术人员根据理论值给定。关于泥质的响应特征值,也可以通过直接给定获得。
对于受孔隙度影响较大的声波测井、密度测井、中子测井的组分响应特征值主要通过物性分析资料与测井资料综合确定,如图3A至3C所示为玄武岩物性孔隙度与对应深度的声波测井曲线、密度测井曲线、中子测井曲线平均值的交会图,线性拟合二者的相关关系及趋势线,趋势线延长线与纵坐标轴交点即为相应的测井响应特征值。同理可确定出安山岩、辉绿岩、凝灰岩、砂岩的密度测井、声波测井、中子测井的响应特征值。
其中,结合图3A、图3B及图3C所示,由物性分析资料与实际的声波测井曲线、密度测井曲线、中子测井曲线获得各交汇图,线性拟合测井曲线值与物性分析孔隙度的趋势关系,然后将孔隙度为0时的曲线值作为该岩性的声波测井、密度测井、中子测井响应特征值。即:玄武岩的密度测井响应特征值为2.6554,玄武岩的中子测井响应特征值为0.0996,玄武岩的声波测井响应特征值为57.387。
对于受孔隙度影响较小或几乎无影响的自然伽马测井、元素俘获能谱测井的组分响应特征值主要通过纯单一岩性段的测井值平均法确定。
结合步骤S104,利用利用测井仪器,获取油田地层的常规测井响应值及元素含量;
首先,利用密度测井仪器、中子测井仪器、声波测井仪器及自然伽马测井仪器获取油田地层的密度响应值、中子响应值、声波响应值、自然伽马测井响应值;
还利用元素俘获能谱测井仪器,获取研究油田地层中的铝、钙、铁、硫、硅的元素含量。
结合步骤S105,利用上述步骤获取的数据,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组。建立方程组最好使响应方程组的方程个数(常规测井响应值及元素含量的个数)比储层组分的个数多一个,使得求解方程组时能得到一个最优化解。
在本实施例中,建立的测井方程组如下:
其中:Vb、Va、Vt、Vd、Vs、Vm、Vo、Vw分别为玄武岩、安山岩、凝灰岩、辉绿岩、砂岩、泥岩、石油、水的相对体积含量;
2.6554为玄武岩的密度测井响应特征值为,0.0996为玄武岩的中子测井响应特征值,2.512为安山岩的密度测井响应特征值,其他系数均为某一地层组分对应的响应特征值;
DEN、CNL、AC、GR分别为利用测井仪器获得的地层中的密度、中子、声波、自然伽马测井曲线响应值;
AL、CA、FE、S、SI分别为元素俘获能谱测井得到的铝、钙、铁、硫、硅元素含量。
结合步骤S106,先用Householder变换求解一个最小二乘解,然后将每个最小二乘解作为共轭梯度优化计算的初始值进行优化迭代计算求解;为了便于利于计算机进行自动岩性解释,在实际资料处理中采用最优化算法进行计算机实现。图4为本发明一具体实施例的多种资料拟合反演岩性定量解释成果图。如图所示,右边第二道为组分拟合反演岩性定量解释后绘制的精细岩性剖面,图中Ⅰ部分为凝灰岩(Tuff)相对体积含量,Ⅱ部分为玄武岩(Basalt)相对体积含量,Ⅲ部分为安山岩(Andesite)相对体积含量,Ⅳ部分为辉绿岩(Diabases)相对体积含量;右边第一道为薄片鉴定成果;左边第一道中,曲线①为自然电位测井曲线、②为自然伽马测井曲线、③为井径测井曲线;左边第二道中,曲线④为密度测井曲线、⑤为中子测井曲线、⑥为声波测井曲线,左边第三道中,曲线⑦、⑧、⑨、⑩分别任为元素俘获能谱测井的钙、铁、铝、硅元素曲线;左边第四道为录井记录的岩性剖面。经过录井岩性剖面与精细岩性解释剖面及薄片鉴定结果的对比,清晰可见,地层组分拟合反演岩性定量解释剖面与薄片鉴定结果完全一致,而录井岩性剖面则显的较为粗略,甚至岩性记录错误。
本发明实施例公开的利用多种资料拟合反演复杂岩性解释方法及装置具备以下显著特点:1、不仅立足于现有成熟的常规测井系列,而且综合应用了较难推广应用的元素俘获能谱测井资料,方法原理简单,技术上容易实现;2、提出了元素俘获能谱测井资料应用的新途径,应用前景广阔;3、精确计算出了火成岩岩性成分相对含量,为岩石命名提供量化的可能;4、实现了火山岩、碳酸盐岩等复杂岩性识别精度从定性到定量化的发展,极大提高了岩性解释的符合率,在一定程度上可用来修正、校正录井剖面记录上的人为误差。利用该发明处理了吐哈油田、大庆油田资料齐全的15口井,处理结果与312块薄片鉴定结果对比,符合率达到93%,较以前提高了10%左右。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (16)
1.一种利用多种资料拟合反演复杂岩性解释方法,其特征在于,所述方法包括:
根据研究工区的地质情况,将储层划分为流体组分、泥质组分以及岩石骨架组分;
给定所述流体组分的流体组分响应特征值及所述泥质组分的泥质组分响应特征值;
根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定所述研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;
利用测井仪器,获取所述研究工区地层的常规测井响应值及元素含量;
根据所述流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;
利用所述测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出所述研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果。
2.根据权利要求1所述的方法,其特征在于,所述流体组分至少包括石油及水,所述石油及水对应的流体组分响应特征值通过直接给定获得。
3.根据权利要求2所述的方法,其特征在于,所述岩石骨架组分至少包括:玄武岩、安山岩、凝灰岩、辉绿岩及砂岩。
4.根据权利要求3所述的方法,其特征在于,所述根根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定所述研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值,包括:
根据测井资料,利用所述油田现场取心、地质描述、物性分析和实验室薄片鉴定方法,确定所述岩石骨架组分的岩性;
利用所述声波测井、密度测井及中子测井,结合物性分析与实际测井曲线,以数据拟合方式来确定所述研究工区的各岩石骨架组分的声波测井响应特征值、密度测井响应特征值及中子测井响应特征值;
利用所述自然伽马测井及元素俘获能谱测井,结合求取实际测井曲线平均值的方法确定所述岩石骨架组分中各个岩性的自然伽马测井响应特征值及元素俘获能谱测井响应特征值。
5.根据权利要求4所述的方法,其特征在于,利用测井仪器,获取所述研究工区地层的常规测井响应值及元素含量,包括:
利用密度测井仪器、中子测井仪器、声波测井仪器及自然伽马测井仪器获取所述研究工区地层的密度响应值、中子响应值、声波响应值、自然伽马测井响应值;
还利用元素俘获能谱测井仪器,获取所述研究工区地层中的铝、钙、铁、硫、硅的元素含量。
6.根据权利要求5所述的方法,其特征在于,所述根据所述流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组,所述测井响应方程组为:
其中,n为所述研究工区中的地层组分的个数;
xi为所述研究工区中的地层组分的相对体积;
Aij为所述研究工区中的流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值;
m为所述测试仪器测量的常规测井响应值及元素含量的个数;
Bj为所述测试仪器测量的常规测井响应值或元素含量。
7.根据权利要求6所述的方法,其特征在于,所述测井响应方程组中元素俘获能谱测井曲线的应用是直接用于建立测井响应方程组,而不需要转换为各元素的相对含量。
8.根据权利要求7所述的方法,其特征在于,所述采用Householder变换和共轭梯度优化算法联合反演出所述研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果,包括:先用Householder变换求解一个最小二乘解,然后将每个最小二乘解作为共轭梯度优化计算的初始值进行优化迭代计算求解。
9.一种利用多种资料拟合反演复杂岩性解释装置,其特征在于,所述装置包括:
储层划分模块,用于根据研究工区的地质情况,将储层划分为流体组分、泥质组分以及岩石骨架组分;
流体组分及泥质组分响应特征值获取模块,用于获取所述流体组分的流体组分响应特征值及所述泥质组分的泥质组分响应特征值,其中,所述流体组分响应特征值以及所述泥质组分响应特征值通过给定直接获得;
岩石骨架组分响应特征值获取模块,根据测井资料,利用油田现场取心、地质描述、物性分析及实验室薄片鉴定方法,确定所述研究工区内的岩石骨架组分类型及个数,确定岩石骨架组分响应特征值;
响应值获取模块,用于利用测井仪器,获取所述研究工区地层的常规测井响应值及元素含量;
测井响应方程组建立模块,用于根据所述流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值、常规测井响应值及元素含量,联合常规测井曲线与元素俘获能谱测井曲线建立测井响应方程组;
岩性定量解释模块,用于利用所述测井响应方程组,采用Householder变换和共轭梯度优化算法联合反演出所述研究工区的地层中的流体组分、泥质组分以及岩石骨架组分的最优相对体积,绘制岩性分析结果。
10.根据权利要求9所述的装置,其特征在于,所述流体组分至少包括石油及水,所述石油及水对应的流体组分响应特征值通过直接给定获得。
11.根据权利要求10所述的装置,其特征在于,所述岩石骨架组分至少包括:玄武岩、安山岩、凝灰岩、辉绿岩及砂岩。
12.根据权利要求11所述的装置,其特征在于,所述岩石骨架组分响应特征值获取模块还用于根据测井资料,利用所述油田现场取心、地质描述、物性分析和实验室薄片鉴定方法,确定所述岩石骨架组分的岩性;
利用所述声波测井、密度测井及中子测井,结合物性分析与实际测井曲线,以数据拟合方式来确定所述研究工区的各岩石骨架组分的声波测井响应特征值、密度测井响应特征值及中子测井响应特征值;
利用所述自然伽马测井及元素俘获能谱测井,结合求取实际测井曲线平均值的方法确定所述岩石骨架组分中各个岩性的自然伽马测井响应特征值及元素俘获能谱测井响应特征值。
13.根据权利要求12所述的装置,其特征在于,所述响应值获取模块还用于利用密度测井仪器、中子测井仪器、声波测井仪器及自然伽马测井仪器获取所述研究工区地层的密度响应值、中子响应值、声波响应值、自然伽马测井响应值;
还利用元素俘获能谱测井仪器,获取所述研究工区地层中的铝、钙、铁、硫、硅的元素含量。
14.根据权利要求13所述的装置,其特征在于,所述测井响应方程组建立模块建立的测井响应方程组包括:
其中,n为所述研究工区中的地层组分的个数;
xi为所述研究工区中的地层组分的相对体积;
Aij为所述研究工区中的流体组分响应特征值、岩石骨架组分响应特征值、泥质组分响应特征值;
m为所述测试仪器测量的常规测井响应值及元素含量的个数;
Bj为所述测试仪器测量的常规测井响应值或元素含量。
15.根据权利要求14所述的装置,其特征在于,所述测井响应方程组建立模块建立的所述测井响应方程组中,元素俘获能谱测井曲线的应用是直接用于建立所述测井响应方程组,而不需要转换为各元素的相对含量。
16.根据权利要求15所述的装置,其特征在于,所述岩性定量解释模块还包括:先用Householder变换求解一个最小二乘解,然后将每个最小二乘解作为共轭梯度优化计算的初始值进行优化迭代计算求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310349346.8A CN103485758A (zh) | 2013-08-12 | 2013-08-12 | 一种利用多种资料拟合反演复杂岩性解释方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310349346.8A CN103485758A (zh) | 2013-08-12 | 2013-08-12 | 一种利用多种资料拟合反演复杂岩性解释方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103485758A true CN103485758A (zh) | 2014-01-01 |
Family
ID=49826304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310349346.8A Pending CN103485758A (zh) | 2013-08-12 | 2013-08-12 | 一种利用多种资料拟合反演复杂岩性解释方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103485758A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105003249A (zh) * | 2015-08-06 | 2015-10-28 | 北京航空航天大学 | 一种基于总流量与电导探针阵列信号的水平井流型识别方法 |
CN105089662A (zh) * | 2015-07-24 | 2015-11-25 | 中国石油天然气股份有限公司 | 碳酸盐岩地层录井综合柱状图的校正方法及系统 |
CN105201490A (zh) * | 2014-06-05 | 2015-12-30 | 中国石油化工股份有限公司 | 一种页岩层段岩石相分析方法 |
CN105221142A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物质量含量的方法 |
CN105257284A (zh) * | 2015-09-23 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种利用元素俘获能谱测井确定凝灰质含量的方法及装置 |
CN105257286A (zh) * | 2015-11-02 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种获取地层岩石组分含量的方法及装置 |
CN105298476A (zh) * | 2014-07-04 | 2016-02-03 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
CN106250619A (zh) * | 2016-07-29 | 2016-12-21 | 中国石油天然气股份有限公司 | 一种确定地层矿物含量的方法和装置 |
CN107144868A (zh) * | 2017-01-24 | 2017-09-08 | 核工业二六大队 | 一种现场多参数快速调查评价煤田放射性辐射环境的方法 |
CN110873904A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
CN112901156A (zh) * | 2021-01-14 | 2021-06-04 | 思凡(上海)石油设备有限公司 | 一种岩屑录井检测系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5128866A (en) * | 1989-09-20 | 1992-07-07 | Chevron Corporation | Pore pressure prediction method |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101149439A (zh) * | 2007-11-13 | 2008-03-26 | 符力耘 | 高分辨率非线性储层物性反演方法 |
CN101329405A (zh) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | 一种简单的多参数地震反演方法 |
CN101787884A (zh) * | 2010-01-28 | 2010-07-28 | 中国石油集团川庆钻探工程有限公司 | 声波孔隙度和中子孔隙度差值储层流体类型判别方法 |
CN101930082A (zh) * | 2009-06-24 | 2010-12-29 | 中国石油集团川庆钻探工程有限公司 | 采用电阻率资料进行储层流体类型判别的方法 |
CN102096107A (zh) * | 2009-12-09 | 2011-06-15 | 中国石油天然气股份有限公司 | 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法 |
CN102454398A (zh) * | 2010-10-28 | 2012-05-16 | 中国石油化工股份有限公司 | 一种适用于低孔低渗储层的气、水层识别方法 |
-
2013
- 2013-08-12 CN CN201310349346.8A patent/CN103485758A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5128866A (en) * | 1989-09-20 | 1992-07-07 | Chevron Corporation | Pore pressure prediction method |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101329405A (zh) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | 一种简单的多参数地震反演方法 |
CN101149439A (zh) * | 2007-11-13 | 2008-03-26 | 符力耘 | 高分辨率非线性储层物性反演方法 |
CN101930082A (zh) * | 2009-06-24 | 2010-12-29 | 中国石油集团川庆钻探工程有限公司 | 采用电阻率资料进行储层流体类型判别的方法 |
CN102096107A (zh) * | 2009-12-09 | 2011-06-15 | 中国石油天然气股份有限公司 | 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法 |
CN101787884A (zh) * | 2010-01-28 | 2010-07-28 | 中国石油集团川庆钻探工程有限公司 | 声波孔隙度和中子孔隙度差值储层流体类型判别方法 |
CN102454398A (zh) * | 2010-10-28 | 2012-05-16 | 中国石油化工股份有限公司 | 一种适用于低孔低渗储层的气、水层识别方法 |
Non-Patent Citations (3)
Title |
---|
崔博文: "基于Householder变换的复参数递推最小二乘估计方法", 《安徽大学学报(自然科学版)》 * |
张兆辉等: "基于地层组分分析的储层孔隙度计算方法研究", 《岩性油气藏》 * |
张兆辉等: "基于地层组分分析的火山岩岩性识别", 《测井技术》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105201490B (zh) * | 2014-06-05 | 2018-04-10 | 中国石油化工股份有限公司 | 一种页岩层段岩石相分析方法 |
CN105201490A (zh) * | 2014-06-05 | 2015-12-30 | 中国石油化工股份有限公司 | 一种页岩层段岩石相分析方法 |
CN105221142A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物质量含量的方法 |
CN105298476A (zh) * | 2014-07-04 | 2016-02-03 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
CN105298476B (zh) * | 2014-07-04 | 2018-06-22 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
CN105089662A (zh) * | 2015-07-24 | 2015-11-25 | 中国石油天然气股份有限公司 | 碳酸盐岩地层录井综合柱状图的校正方法及系统 |
CN105089662B (zh) * | 2015-07-24 | 2018-04-06 | 中国石油天然气股份有限公司 | 碳酸盐岩地层录井综合柱状图的校正方法及系统 |
CN105003249A (zh) * | 2015-08-06 | 2015-10-28 | 北京航空航天大学 | 一种基于总流量与电导探针阵列信号的水平井流型识别方法 |
CN105003249B (zh) * | 2015-08-06 | 2020-09-25 | 北京航空航天大学 | 一种基于总流量与电导探针阵列信号的水平井流型识别方法 |
CN105257284A (zh) * | 2015-09-23 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种利用元素俘获能谱测井确定凝灰质含量的方法及装置 |
CN105257284B (zh) * | 2015-09-23 | 2018-09-04 | 中国石油天然气股份有限公司 | 一种利用元素俘获能谱测井确定凝灰质含量的方法及装置 |
CN105257286A (zh) * | 2015-11-02 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种获取地层岩石组分含量的方法及装置 |
WO2017076151A1 (zh) * | 2015-11-02 | 2017-05-11 | 中国石油天然气股份有限公司 | 一种获取地层岩石组分含量的方法及装置 |
EP3333360A4 (en) * | 2015-11-02 | 2018-10-10 | Petrochina Company Limited | Method of acquiring rock component content of stratum |
CN105257286B (zh) * | 2015-11-02 | 2018-12-25 | 中国石油天然气股份有限公司 | 一种获取地层岩石组分含量的方法及装置 |
US10962679B2 (en) | 2015-11-02 | 2021-03-30 | Petrochina Company Limited | Method of acquiring rock component content of stratum |
CN106250619A (zh) * | 2016-07-29 | 2016-12-21 | 中国石油天然气股份有限公司 | 一种确定地层矿物含量的方法和装置 |
CN107144868A (zh) * | 2017-01-24 | 2017-09-08 | 核工业二六大队 | 一种现场多参数快速调查评价煤田放射性辐射环境的方法 |
CN110873904A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
CN110873904B (zh) * | 2018-09-03 | 2022-02-01 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
CN112901156A (zh) * | 2021-01-14 | 2021-06-04 | 思凡(上海)石油设备有限公司 | 一种岩屑录井检测系统 |
CN112901156B (zh) * | 2021-01-14 | 2022-11-04 | 思凡(上海)石油设备有限公司 | 一种岩屑录井检测系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103485758A (zh) | 一种利用多种资料拟合反演复杂岩性解释方法及装置 | |
CN105257286B (zh) | 一种获取地层岩石组分含量的方法及装置 | |
CN104564041B (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
CN101930082B (zh) | 采用电阻率资料进行储层流体类型判别的方法 | |
CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
CN105221142B (zh) | 一种识别页岩地层矿物质量含量的方法 | |
CN105651966A (zh) | 一种页岩油气优质储层评价方法及参数确定方法 | |
CN103744109B (zh) | 一种无井覆盖区碎屑岩风化壳结构的识别方法 | |
CN104698492B (zh) | 一种计算异常地层压力的方法 | |
US20090248309A1 (en) | Evaluating a reservoir formation | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN104727813B (zh) | 一种源储共生型地层的孔隙度测量方法 | |
CN105426612A (zh) | 一种地层组分最优化确定方法及装置 | |
CN104677734A (zh) | 一种多尺度信息融合的致密地层脆性指数测量方法 | |
CN107703561A (zh) | 一种页岩气水平井中计算有机碳含量的测井方法 | |
CN103867197A (zh) | 复杂岩性天然气层声波时差判别法 | |
CN103867198B (zh) | 碳酸盐岩天然气层地层密度判别法 | |
CN102012526A (zh) | 采用电阻率资料进行储层流体类型判别的方法 | |
CN104615843B (zh) | 一种识别隐伏火山岩型铀矿化信息的物化探综合方法 | |
CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
CN104573150A (zh) | 一种针对富含有机质页岩的各向异性岩石物理建模方法 | |
CN107436452A (zh) | 基于概率神经网络算法的烃源岩预测方法及装置 | |
CN105298476B (zh) | 一种识别页岩地层矿物体积含量的方法 | |
CN103678778B (zh) | 一种放射性物化探信息综合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140101 |
|
RJ01 | Rejection of invention patent application after publication |