CN105653854A - 计算碳酸盐岩烃源岩有机碳含量的方法和装置 - Google Patents
计算碳酸盐岩烃源岩有机碳含量的方法和装置 Download PDFInfo
- Publication number
- CN105653854A CN105653854A CN201511009933.8A CN201511009933A CN105653854A CN 105653854 A CN105653854 A CN 105653854A CN 201511009933 A CN201511009933 A CN 201511009933A CN 105653854 A CN105653854 A CN 105653854A
- Authority
- CN
- China
- Prior art keywords
- organic carbon
- carbon content
- gamma ray
- correlation formula
- sample
- 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
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 111
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 111
- 239000011435 rock Substances 0.000 title claims abstract description 78
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 20
- 230000005251 gamma ray Effects 0.000 claims description 103
- 238000005259 measurement Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 13
- 229930195733 hydrocarbon Natural products 0.000 description 13
- 150000002430 hydrocarbons Chemical class 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 239000003086 colorant Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 239000012857 radioactive material Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明提供了一种计算碳酸盐岩烃源岩有机碳含量的方法和装置,其中,该方法包括:对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;检测得到多个样品中各个样品的有机碳含量;对测井的多个取样深度的自然伽马曲线进行取值,得到多个取样深度的自然伽马测井数据;对多个样品中各个样品的有机碳含量和多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;根据拟合关系式计算待测区域除测井之外的井,在目的层的任意井段的有机碳含量。本发明解决了现有技术中无法通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值的技术问题,达到了快速连续计算碳酸盐岩的TOC值、评价烃源岩质量的目的。
Description
技术领域
本发明涉及地质勘探技术领域,特别涉及一种计算碳酸盐岩烃源岩有机碳含量的方法和装置。
背景技术
有机碳含量(TotalOrganicCarbon,TOC)是烃源岩评价中的重要指标,利用岩心分析数据可获得精确的烃源岩有机碳含量数据,但是由于分析时间周期长、成本较高,并且取样数量有限,不能获得单井连续的有机碳分布曲线,对于目的层烃源岩的评价造成很大困难。基于测井数据评价烃源岩有机碳的分析技术,由于其固有的连续、快速、纵向分辨率高等特点在烃源岩有机碳含量的评价中被广泛应用。
有机质具有较强的吸附放射性物质的特性,富有机质地层(烃源岩)常常具有高放射性强度,因此可以利用二者间的正相关关系进行烃源岩有机碳含量计算。自然伽马测井曲线可在任何井眼环境汇总测量,且测井结果易于环境校正,因此通过自然伽马测井曲线进行有机碳含量技术具有较大优势。利用自然伽马测井进行有机碳含量计算多侧重于泥质烃源岩,然而,很少有通过自然伽马测井方法计算碳酸盐岩烃源岩的有机碳含量。
针对如有通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种计算碳酸盐岩烃源岩有机碳含量的方法,以解决现有技术中无法通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值的技术问题,该方法包括:
对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;检测得到所述多个样品中各个样品的有机碳含量;对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
在一个实施方式中,对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据,包括:对所述自然伽马曲线进行标准化和归一化处理,统一单位为API;对标准化和归一化处理的后自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;相应的,对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式,包括:将所述多个取样深度的自然伽马测井数据以预定个数个API为间隔,从小到大分为多个数据段;获取所述多个数据段中各个数据段的自然伽马测井数据和对应取样深度的有机碳含量;计算所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值;对所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值进行线性拟合,得到所述拟合关系式。
在一个实施方式中,对所述自然伽马曲线进行标准化和归一化处理,包括:
利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理。
在一个实施方式中,根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量,包括:确定所述拟合关系式是否满足预设要求;如果满足所述预设要求,则根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量;如果不满足,则重新拟合所述拟合关系式,直至拟合得到的拟合关系式满足所述预设要求,按照满足所述预设要求的拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
在一个实施方式中,确定所述拟合关系式是否满足预设要求,包括:确定根据所述拟合关系式计算得到的有机碳含量与实测的有机碳含量之间的吻合程度;如有吻合程度满足预设的吻合要求,则确定所述拟合关系式满足预设要求。
在一个实施方式中,在多个取样深度进行取样得到多个样品,包括:以2m为间隔取样得到多个样品。
本发明实施例还提供了一种计算碳酸盐岩烃源岩有机碳含量的装置,以解决现有技术中无法通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值的技术问题,该装置包括:
取样模块,用于对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;
检测模块,用于检测得到所述多个样品中各个样品的有机碳含量;
取值模块,用于对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
拟合模块,用于对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;
计算模块,用于根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
在一个实施方式中,所述取值模块包括:标准化和归一化单元,用于对所述自然伽马曲线进行标准化和归一化处理,统一单位为API;取值单元,用于对标准化和归一化处理的后自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;拟合模块包括:划分单元,用于将所述多个取样深度的自然伽马测井数据以预定个数个API为间隔,从小到大分为多个数据段;获取单元,用于获取所述多个数据段中各个数据段的自然伽马测井数据和对应取样深度的有机碳含量;第一计算单元,用于计算所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值;拟合单元,用于对所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值进行线性拟合,得到所述拟合关系式。
在一个实施方式中,所述归一化单元具体用于利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理。
在一个实施方式中,所述计算模块包括:确定单元,用于确定所述拟合关系式是否满足预设要求;第二计算单元,用于在确定满足所述预设要求的情况下,根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量;第三计算单元,用于在确定不满足所述预设要求的情况下,重新拟合所述拟合关系式,直至拟合得到的拟合关系式满足所述预设要求,按照满足所述预设要求的拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
在本发明实施例中,通过获取待测区域目的层段碳酸盐岩烃源岩的自然伽马测井数据和有机碳含量,并对自然伽马测井数据和有机碳含量进行拟合,得到拟合关系式,从而建立有机碳含量与自然伽马测井数据之间的关系,基于该关系便可以实现目的层的任意井段的有机碳含量的计算,通过上述方式解决了现有技术中无法通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值的技术问题,达到了快速连续计算碳酸盐岩的TOC值、评价烃源岩质量、分析区域有机碳丰度分布特征、进行油气资源评价的目的。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是根据本发明实施例的计算碳酸盐岩烃源岩有机碳含量的方法流程图;
图2是根据本发明实施例的利用自然伽马测井计算碳酸盐岩烃源岩有机碳含量的具体实例示意图;
图3是根据本发明实施例的中四川盆地东北地区下二叠统栖霞-茅口组GR-实测TOC回归关系图;
图4是根据本发明实施例的四川盆地东北地区张16井下二叠统实测TOC与计算TOC对比图;
图5是根据本发明实施例的计算碳酸盐岩烃源岩有机碳含量的装置结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
为了通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值,在本例中提供了一种计算碳酸盐岩烃源岩有机碳含量的方法,用以评价碳酸盐岩烃源岩的质量,寻找优质烃源岩层段,分析研究区域内有机碳平面分布特征,为油气资源评价提供详实数据支持,如图1所示,该方法包括以下步骤:
步骤101:对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;
在选取目的层碳酸盐岩烃源岩样品时,应该尽量选取岩心较为连续、完整的井,选择颜色较暗的岩石样品,并且取样间隔尽量均匀,例如,可以以2m为间隔取样,如果该井岩心段较短,可适当加大取样密度。在取样时,可以同时记录下取样井号、深度和岩性,还可以对典型的样品进行拍照并记录照片序号。
在取样过程中,尽量选择多口井的岩心样品,以便在面上覆盖、控制住研究区范围,并且如果分析样本数量比较多可以使得所取样品具有较好的代表性。
步骤102:检测得到所述多个样品中各个样品的有机碳含量;
在实现的过程中,可以将取得的碳酸盐岩样品,送到有资质的实验室,进行总有机碳(TOC)的测定。
步骤103:对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
步骤104:对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;
具体地,在上述步骤103和步骤104中,可以获取相同井对应取样深度的自然伽马(GR)数据;将GR数据以5个API单位为间隔,分成若干数据段,挑出每个数据段中的GR和对应的TOC数据点,分别计算得到各数据段GR的平均值(X)和TOC数据的平均值(Y),然后,得到Y和X二者之间的关系式:
Y=aX+b
并求取X和Y之间的相关系数R2。
即,在一个实施方式中,可以利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理,统一单位为API,其中,API为自然伽马曲线的常规单位,然后对归一化处理的后自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据,进一步的,将所述多个取样深度的自然伽马测井数据以预定个数个API为间隔,从小到大分为多个数据段;获取所述多个数据段中各个数据段的自然伽马测井数据和对应取样深度的有机碳含量;计算所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值;对所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值进行线性拟合,得到所述拟合关系式。
例如:按照深度顺序,将测试分析得到的TOC数据和对应深度的自然伽马数据值排列好在excel数据表中。将excel数据表中的自然伽马数据从小到大排列,以5个API单位为间隔,将自然伽马数据分成若干个数据段1,2,3,……,同时相对应的TOC数据也分文若干数据段1,2,3……。求取自然伽马和TOC各数据段的平均值X1,X2,X3,……和Y1,Y2,Y3,……。利用excel数据表自带的公式拟合工具,对Y1,Y2,Y3,……和X1,X2,X3,……,进行线性关系拟合,得到公式Y=aX+b,并计算方差。
然而,上述以5个API单位为间隔、利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理等,仅是为了更好的说明本发明,并不构成对本发明的不当限定,在具体实现的时候,还可以采用其它数值API的间隔,且可以采用其它的方式进行标准化和归一化处理。
步骤105:根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
考虑到有时候因为数据的偏差或者选取样品的特殊性,使得拟合得到的关系式并不能满足一般地层的要求,为此为了使得拟合得到关系式可以更好更精确地进行目的层任意井段的有机碳含量的计算,可以在得到拟合关系式,进行有机碳含量计算前,确定拟合关系式是否满足预设的要求,具体地,可以确定拟合关系式是否满足预设要求;如果满足所述预设要求,则根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量;如果不满足,则重新拟合所述拟合关系式,直至拟合得到的拟合关系式满足所述预设要求,按照满足所述预设要求的拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。即,在获得自然伽马值与TOC的关系式后,利用该关系式对同研究区域的其他井进行TOC计算,并且与该井统一深度的样品实测TOC值进行比较,两者值相近,则表明拟合关系式可靠,如果两者值相差较大,则表明拟合关系式不可靠。
在上例中,通过获取待测区域目的层段碳酸盐岩烃源岩的自然伽马测井数据和有机碳含量,并对自然伽马测井数据和有机碳含量进行拟合,得到拟合关系式,从而建立有机碳含量与自然伽马测井数据之间的关系,基于该关系便可以实现目的层的任意井段的有机碳含量的计算,通过上述方式解决了现有技术中无法通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值的技术问题,达到了快速连续计算碳酸盐岩的TOC值、评价烃源岩质量、分析区域有机碳丰度分布特征、进行油气资源评价的目的。
下面结合一个具体实施例对基于自然伽马测井数据计算碳酸盐岩烃源岩有机碳含量的方法进行说明,然而,值得注意的是,该具体实施例仅是为了更好地说明本发明,并不构成对本发明的不当限定。
为了利用自然伽马测井数据计算出碳酸盐岩烃源岩的有机碳含量,寻找优质烃源岩发育层段,在本例中提供了一种利用自然伽马测井数据计算碳酸盐岩烃源岩的有机碳含量的方法,如图2所示,包括以下步骤:
步骤201:选取有碳酸盐岩烃源岩岩心的井,对目的层样品进行取样,记录岩性、深度,并编好序号,在实验室进行有机碳(TOC)测试,且在选取烃源岩岩心样品的时候,应当尽量多选取样品,保证样品数量充足且具有代表性,数据越多,越具有代表性和普遍性;
即,首先选取待测区域有目的层取心的井,对这些井的目的层碳酸盐岩烃源岩岩心样品进行采集,并记录下井号、层位、深度并对采集的样品依次编上序号。在采集完碳酸盐岩样品后,将样品进行有机碳分析测试,具体的分析测试标准可以参考国标GB/T19145-2003,以实现沉积岩中总有机碳的测定。在具体操作的时候,可以选取颜色较暗的烃源岩样品,适当加大采样密度,例如可以以2m作为间隔采样。
步骤202:依据取样的深度,对自然伽马曲线进行取值,获取取样井对应取样深度的自然伽马(GR)测井数据;
步骤203:将获取的GR数据以5个API单位为间隔,从小到大,分成若干个数据段,挑出每个数据段中的GR值和对应深度的TOC值,分别计算各数据段的GR(X)和对应TOC平均值(Y);
步骤204:拟合Y与X二者的关系,得到拟合关系式,Y=aX+b,并求取两者之间的相关系数(R2);
即,可以利用测井软件CIFLog或Geolog等对自然伽马曲线进行标准化和归一化处理,统一单位为API,然后依据选取烃源岩岩心的深度,对自然伽马曲线进行取值。按照深度顺序,将测试分析得到的TOC数据和对应深度的自然伽马数据值排列好在excel数据表中。将excel数据表中的自然伽马数据从小到大排列,以5个API单位为间隔,将自然伽马数据分成若干个数据段1,2,3,……,同时相对应的TOC数据也分文若干数据段1,2,3……。求取自然伽马和TOC各数据段的平均值X1,X2,X3,……和Y1,Y2,Y3,……。利用excel数据表自带的公式拟合工具,对Y1,Y2,Y3,……和X1,X2,X3,……,进行线性关系拟合,得到公式Y=aX+b,并计算方差。
如图3所示,为实测的TOC(Y)值与自然伽马数据(X)的相关关系示意图,具体拟合得到的公式为Y=0.0028x-0.3343,R2=0.8789。
步骤205:选取本地区新井目的层碳酸盐岩烃源岩样品,采样并测试TOC值,比较计算得到的TOC值和实测的TOC值,如果两者接近,则证明上述拟合得到的关系式可靠。
即,在具体实施时,为了验证拟合得到的公式的可靠性,可以重新在研究区采集新井的目的层碳酸盐岩烃源岩样品,在实验室分析测试其TOC含量,再利用公式和其对应深度的自然伽马值,计算其TOC值,比较实测TOC和计算TOC,二者相近,则证明拟合的公式可靠,依据此公式即能对目的层段任意深度的碳酸盐岩烃源岩进行TOC计算。如果吻合度较差,则需要返回到步骤203,对公式进行重新拟合,或者再次从步骤201做起,通过选取更多的烃源岩样品,进行实测,到步骤202、步骤203、步骤204,再次拟合关系式,直至吻合度较高为止。
如图4所示,为川东地区张16井采用上述方式计算得到的TOC与实测的TOC的对比示意图,从图4中可以看出,两者的吻合度关系比较高。
步骤206:依据上述拟合关系式和检测到的GR值,计算待检测区域其他井目的层任意井段的TOC值。
在上例中,利用自然伽马曲线对碳酸盐岩烃源岩的有机碳进行计算,这对于快速、便捷统计烃源岩厚度、评价烃源岩质量,预测资源量具有重要的理论和实践意义。
基于同一发明构思,本发明实施例中还提供了一种计算碳酸盐岩烃源岩有机碳含量的装置,如下面的实施例所述。由于计算碳酸盐岩烃源岩有机碳含量的装置解决问题的原理与计算碳酸盐岩烃源岩有机碳含量的方法相似,因此计算碳酸盐岩烃源岩有机碳含量的装置的实施可以参见计算碳酸盐岩烃源岩有机碳含量的方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图5是本发明实施例的计算碳酸盐岩烃源岩有机碳含量的装置的一种结构框图,如图5所示,包括:取样模块501、检测模块502、取值模块503、拟合模块504和计算模块505,下面对该结构进行说明。
取样模块501,用于对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;
检测模块502,用于检测得到所述多个样品中各个样品的有机碳含量;
取值模块503,用于对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
拟合模块504,用于对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;
计算模块505,用于根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
在一个实施方式中,取值模块503可以包括:标准化和归一化单元,用于对所述自然伽马曲线进行标准化和归一化处理,统一单位为API;取值单元,用于对标准化和归一化处理的后自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;拟合模块504可以包括:划分单元,用于将所述多个取样深度的自然伽马测井数据以预定个数个API为间隔,从小到大分为多个数据段;获取单元,用于获取所述多个数据段中各个数据段的自然伽马测井数据和对应取样深度的有机碳含量;第一计算单元,用于计算所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值;拟合单元,用于对所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值进行线性拟合,得到所述拟合关系式。
在一个实施方式中,标准化和归一化单元具体用于可以利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理。
在一个实施方式中,计算模块505可以包括:确定单元,用于确定所述拟合关系式是否满足预设要求;第二计算单元,用于在确定满足所述预设要求的情况下,根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量;第三计算单元,用于在确定不满足所述预设要求的情况下,重新拟合所述拟合关系式,直至拟合得到的拟合关系式满足所述预设要求,按照满足所述预设要求的拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
从以上的描述中,可以看出,本发明实施例实现了如下技术效果:通过获取待测区域目的层段碳酸盐岩烃源岩的自然伽马测井数据和有机碳含量,并对自然伽马测井数据和有机碳含量进行拟合,得到拟合关系式,从而建立有机碳含量与自然伽马测井数据之间的关系,基于该关系便可以实现目的层的任意井段的有机碳含量的计算,通过上述方式解决了现有技术中无法通过直接的自然伽马曲线计算碳酸盐岩烃源岩的TOC值的技术问题,达到了快速连续计算碳酸盐岩的TOC值、评价烃源岩质量、分析区域有机碳丰度分布特征、进行油气资源评价的目的。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种计算碳酸盐岩烃源岩有机碳含量的方法,其特征在于,包括:
对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;
检测得到所述多个样品中各个样品的有机碳含量;
对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;
根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
2.如权利要求1所述的方法,其特征在于,
对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据,包括:
对所述自然伽马曲线进行标准化和归一化处理,统一单位为API;
对标准化和归一化处理的后自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
相应的,对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式,包括:
将所述多个取样深度的自然伽马测井数据以预定个数个API为间隔,从小到大分为多个数据段;
获取所述多个数据段中各个数据段的自然伽马测井数据和对应取样深度的有机碳含量;
计算所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值;
对所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值进行线性拟合,得到所述拟合关系式。
3.如权利要求2所述的方法,其特征在于,对所述自然伽马曲线进行标准化和归一化处理,包括:
利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理。
4.如权利要求1所述的方法,其特征在于,根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量,包括:
确定所述拟合关系式是否满足预设要求;
如果满足所述预设要求,则根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量;
如果不满足,则重新拟合所述拟合关系式,直至拟合得到的拟合关系式满足所述预设要求,按照满足所述预设要求的拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
5.如权利要求4所述的方法,其特征在于,确定所述拟合关系式是否满足预设要求,包括:
确定根据所述拟合关系式计算得到的有机碳含量与实测的有机碳含量之间的吻合程度;
如有吻合程度满足预设的吻合要求,则确定所述拟合关系式满足预设要求。
6.如权利要求1至5中任一项所述的方法,其特征在于,在多个取样深度进行取样得到多个样品,包括:
以2m为间隔取样得到多个样品。
7.一种计算碳酸盐岩烃源岩有机碳含量的装置,其特征在于,包括:
取样模块,用于对待测区域内测井的目的层的碳酸盐岩烃源岩,在多个取样深度进行取样得到多个样品;
检测模块,用于检测得到所述多个样品中各个样品的有机碳含量;
取值模块,用于对所述测井的所述多个取样深度的自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
拟合模块,用于对所述多个样品中各个样品的有机碳含量和所述多个取样深度的自然伽马测井数据进行线性拟合,得到拟合关系式;
计算模块,用于根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
8.如权利要求7所述的装置,其特征在于,
所述取值模块包括:
标准化和归一化单元,用于对所述自然伽马曲线进行标准化和归一化处理,统一单位为API;
取值单元,用于对标准化和归一化处理的后自然伽马曲线进行取值,得到所述多个取样深度的自然伽马测井数据;
拟合模块包括:
划分单元,用于将所述多个取样深度的自然伽马测井数据以预定个数个API为间隔,从小到大分为多个数据段;
获取单元,用于获取所述多个数据段中各个数据段的自然伽马测井数据和对应取样深度的有机碳含量;
第一计算单元,用于计算所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值;
拟合单元,用于对所述多个数据段中各个数据段的自然伽马测井数据的平均值和对应的有机碳含量的平均值进行线性拟合,得到所述拟合关系式。
9.如权利要求8所述的装置,其特征在于,所述归一化单元具体用于利用测井软件CIFLog或Geolog对所述自然伽马曲线进行标准化和归一化处理。
10.如权利要求7所述的装置,其特征在于,所述计算模块包括:
确定单元,用于确定所述拟合关系式是否满足预设要求;
第二计算单元,用于在确定满足所述预设要求的情况下,根据所述拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量;
第三计算单元,用于在确定不满足所述预设要求的情况下,重新拟合所述拟合关系式,直至拟合得到的拟合关系式满足所述预设要求,按照满足所述预设要求的拟合关系式计算所述待测区域除所述测井之外的井,在所述目的层的任意井段的有机碳含量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009933.8A CN105653854A (zh) | 2015-12-29 | 2015-12-29 | 计算碳酸盐岩烃源岩有机碳含量的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009933.8A CN105653854A (zh) | 2015-12-29 | 2015-12-29 | 计算碳酸盐岩烃源岩有机碳含量的方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105653854A true CN105653854A (zh) | 2016-06-08 |
Family
ID=56477742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511009933.8A Pending CN105653854A (zh) | 2015-12-29 | 2015-12-29 | 计算碳酸盐岩烃源岩有机碳含量的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105653854A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627381A (zh) * | 2017-03-23 | 2018-10-09 | 中国石油化工股份有限公司 | 一种钻井蒸发盐岩地层中碳酸盐岩样品的筛选方法 |
CN110716024A (zh) * | 2018-07-11 | 2020-01-21 | 中国石油天然气股份有限公司 | 物性曲线的确定方法及装置 |
CN112782778A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油天然气股份有限公司 | 碳酸盐岩台地咸化环境烃源岩发育确定方法及装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011133421A2 (en) * | 2010-04-21 | 2011-10-27 | Baker Hughes Incorporated | Method of predicting source rock thermal maturity from log responses |
CN103114840A (zh) * | 2013-01-09 | 2013-05-22 | 中国石油天然气股份有限公司 | 一种高-过高成熟页岩有机碳含量计算方法及装置 |
CN104278991A (zh) * | 2013-07-09 | 2015-01-14 | 中国石油化工股份有限公司 | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 |
US20150094960A1 (en) * | 2013-10-02 | 2015-04-02 | Schlumberger Technology Corporation | Quality metrics for tight oil reservoirs |
CN104533405A (zh) * | 2014-12-19 | 2015-04-22 | 中国石油天然气集团公司 | 泥页岩储层有机碳含量的检测方法和装置 |
CN104636588A (zh) * | 2014-08-25 | 2015-05-20 | 中国石油天然气股份有限公司 | 计算烃源岩中有机碳含量的方法及装置 |
CN104656162A (zh) * | 2015-02-13 | 2015-05-27 | 西安石油大学 | 一种富铀烃源岩铀含量及有机碳含量确定方法 |
CN103792338B (zh) * | 2014-01-28 | 2015-11-18 | 中国石油天然气股份有限公司 | 一种烃源岩有机碳含量的确定方法和装置 |
CN105158796A (zh) * | 2015-09-16 | 2015-12-16 | 中国石油天然气股份有限公司 | 确定toc含量的方法和装置 |
-
2015
- 2015-12-29 CN CN201511009933.8A patent/CN105653854A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011133421A2 (en) * | 2010-04-21 | 2011-10-27 | Baker Hughes Incorporated | Method of predicting source rock thermal maturity from log responses |
CN103114840A (zh) * | 2013-01-09 | 2013-05-22 | 中国石油天然气股份有限公司 | 一种高-过高成熟页岩有机碳含量计算方法及装置 |
CN104278991A (zh) * | 2013-07-09 | 2015-01-14 | 中国石油化工股份有限公司 | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 |
US20150094960A1 (en) * | 2013-10-02 | 2015-04-02 | Schlumberger Technology Corporation | Quality metrics for tight oil reservoirs |
CN103792338B (zh) * | 2014-01-28 | 2015-11-18 | 中国石油天然气股份有限公司 | 一种烃源岩有机碳含量的确定方法和装置 |
CN104636588A (zh) * | 2014-08-25 | 2015-05-20 | 中国石油天然气股份有限公司 | 计算烃源岩中有机碳含量的方法及装置 |
CN104533405A (zh) * | 2014-12-19 | 2015-04-22 | 中国石油天然气集团公司 | 泥页岩储层有机碳含量的检测方法和装置 |
CN104656162A (zh) * | 2015-02-13 | 2015-05-27 | 西安石油大学 | 一种富铀烃源岩铀含量及有机碳含量确定方法 |
CN105158796A (zh) * | 2015-09-16 | 2015-12-16 | 中国石油天然气股份有限公司 | 确定toc含量的方法和装置 |
Non-Patent Citations (7)
Title |
---|
BL MEYER ,ET AL.,: "Identification of source rocks on wireline logs by density-resistivity and sonic transit time-resistivity crossplots", 《 AAPG BULL》 * |
孟召平,等;: "页岩气储层有机碳含量与测井参数的关系及预测模型", 《煤炭学报》 * |
张立鹏,等;: "用测井资料识别烃源岩", 《测井技术》 * |
李松峰,等;: "普光地区须家河组烃源岩地球物理预测", 《断块油气田》 * |
杨惠连,等,编著;: "《误差理论与数据处理 1992年8月第1版》", 31 August 1992 * |
袁超,等;: "地层有机碳含量测井评价方法综述", 《地球物理学进展》 * |
郭彦省,: "基于非线性学习理论的非常规储层基本参数测井评价", 《中国博士学位论文全文数据库基础科学辑(月刊)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627381A (zh) * | 2017-03-23 | 2018-10-09 | 中国石油化工股份有限公司 | 一种钻井蒸发盐岩地层中碳酸盐岩样品的筛选方法 |
CN110716024A (zh) * | 2018-07-11 | 2020-01-21 | 中国石油天然气股份有限公司 | 物性曲线的确定方法及装置 |
CN112782778A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油天然气股份有限公司 | 碳酸盐岩台地咸化环境烃源岩发育确定方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
CN105221133A (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
AU2010263041A1 (en) | Source rock volumetric analysis | |
NO339596B1 (no) | Kalibrering av en permeabilitetsmodell ved bruk av målinger av reservoarfluidstrømning | |
CN105114067A (zh) | 一种岩性测井相的方法 | |
CN105653854A (zh) | 计算碳酸盐岩烃源岩有机碳含量的方法和装置 | |
Naylor et al. | Bias in detrital fission track grain-age populations: Implications for reconstructing changing erosion rates | |
Guo et al. | Identification of homogeneous region boundaries of fractured rock masses in candidate sites for Chinese HLW repository | |
Selzer et al. | RaDeCC Reader: Fast, accurate and automated data processing for Radium Delayed Coincidence Counting systems | |
Francesconi et al. | Reservoir rock types application-Kashagan | |
Martin et al. | Developing a quantitative mudrock calibration for a handheld energy dispersive X-ray fluorescence spectrometer | |
CN110685676B (zh) | 一种定量识别优质页岩段的方法 | |
CN114993759B (zh) | 用于放射性矿产勘查的岩粉取样方法 | |
Hussain et al. | Linking geochemical and mechanical properties of rock samples using new non-destructive techniques | |
Zohreh et al. | The evaluation of borehole imaging result comparing with cores in Sarvak fractured and non-fractured reservoir | |
Skupio et al. | Archival gamma ray logs standardization by nondestructive core measurements of the low-radioactivity rocks | |
Descalzi et al. | Synergetic log and core data treatment through Cluster Analysis: a methodology to improve reservoir description | |
ASFAHANI | Phosphate prospecting using natural gamma ray well logging in the Khneifiss Mine, Syria | |
Slatt et al. | Workflow for stratigraphic characterization of unconventional gas shales | |
Olesen et al. | Dual-burst thermal decay time logging overview and examples | |
Ighodalo et al. | Integrated NMR Fluid Characterization Guides Stimulation in Tight Sand Reservoirs | |
Eyvazzadeh et al. | Optimizing Petrophysics: The Ghawar Field, Saudi Arabia | |
Menon et al. | Novel Approach to Matrix Corrections in Complex Khuff Deep Gas Reservoirs, offshore Abu Dhabi | |
Alward et al. | Integrated Workflow for Petrophysical Evaluation in Heterogeneous Carbonate Reservoirs | |
Eltaher et al. | In-Situ Oil Carbon Density Characterization for Enhanced Reservoir Saturation Monitoring Using Carbon-Oxygen Logs |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160608 |