CN104237966A - 测算地下有机质孔隙度的方法和装置 - Google Patents
测算地下有机质孔隙度的方法和装置 Download PDFInfo
- Publication number
- CN104237966A CN104237966A CN201410443833.5A CN201410443833A CN104237966A CN 104237966 A CN104237966 A CN 104237966A CN 201410443833 A CN201410443833 A CN 201410443833A CN 104237966 A CN104237966 A CN 104237966A
- Authority
- CN
- China
- Prior art keywords
- organic
- porosity
- unit
- well point
- represent
- 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.)
- Granted
Links
- 239000005416 organic matter Substances 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 40
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 164
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 157
- 239000011435 rock Substances 0.000 claims description 58
- 229930195733 hydrocarbon Natural products 0.000 claims description 34
- 150000002430 hydrocarbons Chemical class 0.000 claims description 34
- 239000004215 Carbon black (E152) Substances 0.000 claims description 33
- 238000013481 data capture Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract 2
- 238000011084 recovery Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 8
- 239000007789 gas Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000005498 polishing Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000003921 oil Substances 0.000 description 5
- 239000007787 solid Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000010779 crude oil Substances 0.000 description 4
- 239000011148 porous material Substances 0.000 description 3
- 238000000197 pyrolysis Methods 0.000 description 3
- 239000003079 shale oil Substances 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种测算地下有机质孔隙度的方法和装置,其中,该方法包括:在研究区中选取多个井点,获取所述研究区的有机碳恢复系数与有机质热演化程度之间的关系曲线;获取所述多个井点中各个井点的有机质热演化程度,对各个井点执行以下操作:计算该井点的原始有机碳重量百分含量;计算该井点的最大有机质孔隙度;计算得到该井点的观测有机质孔隙度;对该井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均,计算得到该井点的平均有机质孔隙度。本发明解决了现有技术中对有机质孔隙度测算结果不准确的技术问题,提高了测算结果的准确性。
Description
技术领域
本发明涉及页岩储层物性研究和页岩油气资源量技术领域,特别涉及一种测算地下有机质孔隙度的方法和装置。
背景技术
页岩油气是目前各种非常规油气类型中最受关注的一种,是各大油公司正在投入研究的潜在油气资源接替领域,测算地下状态下的有机质孔隙度,对页岩储层的物性研究和页岩油气资源量的计算非常重要。
有机质孔隙是指泥页岩中有机质在生烃过程中,固态干酪根转化成非固态烃类后残留的孔隙空间,有机质孔隙度是指有机质孔隙体积占岩石体积的百分比。
目前,国外有机质孔隙度的测算方法主要有以下三种:
1)利用测井曲线计算法
2012年,Modica和Lapierre提出一种计算PhiK模型,该模型认为:有机质孔隙度与原始有机碳、可转化碳比例和转化率成正比,与测井曲线解释出的岩石密度成正比,与测井曲线解释出的干酪根密度成反比。用该方法测算美国粉河盆地(PowderRiver Basin)Mowry页岩的有机质孔隙度,取得了良好的效果。
然而,PhiK模型计算结果为理论上最大有机质孔隙度,没有考虑后期有机质孔隙度被破坏的实际地质情况,因此其结果属于一种理想状态下的最大有机质孔隙度。
2)统计方法
2013年,郭秋麟等提出了一种利用Ar离子抛光薄片扫描电镜照片统计结果,拟合有机质面孔率与Ro的关系,然后,利用面孔率计算有机质孔隙度的方法。
然而,该方法实际上是一种预测地面条件下Ar离子抛光薄片观测孔隙度的,在有机质高演化阶段预测结果较好,在中低演化阶段(Ro<1.2)预测效果较差。这主要是因为地下高温高压与地面常温常压的环境差别较大,在地下被液态烃填充的孔隙,因轻质组分挥发在地面可能被沥青充填。
3)实验与Ar离子抛光薄片分析相结合的方法
2013年,陈方文等提出了一种泥页岩有机孔隙度检测方法,该方法包括:以代表性泥页岩样品和原油样品的热模拟实验为基础,利用化学动力学方法计算干酪根成油、干酪根成气和原油裂解成气的化学动力学参数,并结合目地层埋藏史和热史,确定研究层段泥页岩干酪根成油、干酪根成气和原油裂解成气的转化率;再利用目地层泥页岩残余氢指数和残余有机碳数据,结合干酪根成油、干酪根成气和原油裂解成气转化率,恢复目地层泥页岩原始氢指数和原始有机碳;利用目地层泥页岩样品的Ar离子抛光薄片分析泥页岩有机孔隙压缩系数。最终计算出目地层段泥页岩样品的有机孔隙度。
该方法实际是一种对PhiK模型的改进,该方法将理想状态下的最大有机质孔隙度折算成地面Ar离子抛光薄片观测孔隙度,也同样会出现在中低演化阶段(Ro<1.2),即生油窗阶段,预测效果较差的问题。
由上述分析可知,目前的几种有机质孔隙度的测算方法都存在预测结果偏大或偏小的问题,导致了有机质孔隙度的测算结果不准确。
发明内容
本发明实施例提供了一种测算地下有机质孔隙度的方法,以提高有机质孔隙度测算结果的准确性,该方法包括:
在研究区中选取多个井点,获取所述研究区的有机碳恢复系数与有机质热演化程度之间的关系曲线;
模拟所述研究区泥页岩的有机质热演化史,得到所述多个井点中各个井点的有机质热演化程度;
对各个井点执行以下操作:
获取该井点的有机碳重量百分含量、烃源岩密度和有机碳密度;
根据该井点的有机质热演化程度,从所述关系曲线中查找出该井点的碳恢复系数,根据获取的该井点的有机碳重量百分含量和查找到的该井点的碳恢复系数,计算该井点的原始有机碳重量百分含量;
根据获取的该井点的有机碳重量百分含量、原始有机碳重量百分含量、烃源岩密度和有机碳密度,计算该井点的最大有机质孔隙度;
根据该井点的有机质热演化程度确定面孔率,并根据确定的该井点的面孔率和获取的该井点的有机碳重量百分含量,计算得到该井点的观测有机质孔隙度;
对该井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均,计算得到该井点的平均有机质孔隙度。
在一个实施例中,按照以下公式计算原始有机碳重量百分含量:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,Cf表示碳恢复系数。
在一个实施例中,按照以下公式计算最大有机质孔隙度:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示最大有机质孔隙度,单位为%,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
在一个实施例中,按照以下公式确定面孔率:
psurf=a×Ro-b
其中,psurf表示面孔率,单位为%,a和b表示回归系数,Ro表示有机质热演化程度。
在一个实施例中,按照以下公式计算观测有机质孔隙度;
其中,Φmin表示观测有机质孔隙度,单位为%,psurf表示面孔率,单位为%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
在一个实施例中,按照以下公式计算平均有机质孔隙度:
Φom=(Φmin+Φmax)/2
其中,Φom表示各平均有机质孔隙度,单位为%,Φmin表示观测有机质孔隙度,单位为%,Φmax表示最大有机质孔隙度,单位为%。
本发明实施例还提供了一种测算地下有机质孔隙度的装置,以提高有机质孔隙度测算结果的准确性,该装置包括:
数据获取单元,用于在研究区中选取多个井点,获取所述研究区的有机碳恢复系数与有机质热演化程度之间的关系曲线;
模拟单元,用于模拟所述研究区泥页岩的有机质热演化史,得到所述多个井点中各个井点的有机质热演化程度;
执行操作单元,用于对各个井点执行以下操作:
获取该井点的有机碳重量百分含量、烃源岩密度和有机碳密度;
根据该井点的有机质热演化程度,从所述关系曲线中查找出该井点的碳恢复系数,根据获取的该井点的有机碳重量百分含量和查找到的该井点的碳恢复系数,计算该井点的原始有机碳重量百分含量;
根据获取的该井点的有机碳重量百分含量、原始有机碳重量百分含量、烃源岩密度和有机碳密度,计算该井点的最大有机质孔隙度;
根据该井点的有机质热演化程度确定面孔率,并根据确定的该井点的面孔率和获取的该井点的有机碳重量百分含量,计算得到该井点的观测有机质孔隙度;
对该井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均,计算得到该井点的平均有机质孔隙度。
在一个实施例中,所述执行操作单元具体用于按照以下公式计算原始有机碳重量百分含量:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,Cf表示碳恢复系数。
在一个实施例中,所述执行操作单元具体用于按照以下公式计算最大有机质孔隙度:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示最大有机质孔隙度,单位为%,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
在一个实施例中,所述执行操作单元具体用于按照以下公式确定面孔率:
psurf=a×Ro-b
其中,psurf表示面孔率,单位为%,a和b表示回归系数,Ro表示有机质热演化程度。
在一个实施例中,所述执行操作单元具体用于按照以下公式计算观测有机质孔隙度;
其中,Φmin表示观测有机质孔隙度,单位为%,psurf表示面孔率,单位为%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
在一个实施例中,所述执行操作单元具体用于按照以下公式计算平均有机质孔隙度:
Φom=(Φmin+Φmax)/2
其中,Φom表示各平均有机质孔隙度,单位为%,Φmin表示观测有机质孔隙度,单位为%,Φmax表示最大有机质孔隙度,单位为%。
在本发明实施例中,将各个井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均的结果,作为研究区中各个井点的平均有机质孔隙度,从而解决了现有技术中对有机质孔隙度测算结果不准确的技术问题,提高了测算结果的准确性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例的测算地下有机质孔隙度的方法流程图;
图2是本发明实施例的研究区泥页岩TOC分布示意图;
图3是本发明实施例的研究区有机碳恢复系数与Ro关系曲线图;
图4是本发明实施例的研究区Ro分布示意图;
图5是本发明实施例的研究区的泥页岩TOC与iTOC的对比图;
图6是本发明实施例的研究区的有机质孔隙度分布示意图;
图7是本发明实施例的测算地下有机质孔隙度的装置结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
考虑到现有的测算地下有机质孔隙度中,利用测井曲线计算法得到的结果是理想状态下的孔隙度,结果偏大,统计方法和实验与Ar离子抛光薄片分析相结合的方法得到的结果是观测孔隙度,在生油窗阶段(即,Ro<1.2时)结果偏小。
在本例中,提供了一种测算地下有机质孔隙度的方法,如图1所示,包括以下步骤:
步骤101:在研究区中选取多个井点,获取选取的多个井点中各个井点的有机碳重量百分含量、烃源岩密度和有机碳密度,并获取有机碳恢复系数与有机质热演化程度之间的关系曲线;
步骤102:模拟所述研究区泥页岩的有机质热演化史,得到所述多个井点中各个井点的有机质热演化程度;
步骤103:根据各个井点的有机质热演化程度,从所述关系曲线中查找出各个井点的碳恢复系数,根据获取的各个井点的有机碳重量百分含量和查找到的各个井点的碳恢复系数,计算所述各个井点的原始有机碳重量百分含量;
步骤104:根据获取的各个井点的有机碳重量百分含量、原始有机碳重量百分含量、烃源岩密度和有机碳密度,计算各个井点的最大有机质孔隙度;
步骤105:根据各个井点的有机质热演化程度确定面孔率,并根据确定的各个井点的面孔率和获取的各个井点的有机碳重量百分含量,计算得到各个井点的观测有机质孔隙度;
步骤106:对各个井点的最大有机质孔隙度,和各个井点的观测有机质孔隙度进行加权平均,计算得到所述研究区中各个井点的平均有机质孔隙度。
在本例中,将各个井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均的结果,作为研究区中各个井点的平均有机质孔隙度,从而解决了现有技术中对有机质孔隙度测算结果不准确的技术问题,提高了测算结果的准确性。
在上述步骤102中,可以利用通用盆地模拟软件系统模拟泥页岩有机质热演化史,这样相对于当天实测出的数据更为准确。
具体的,在步骤103中,可以按照以下公式计算各个井点的原始有机碳重量百分含量:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的各个井点的有机碳重量百分含量,单位为wt%,Cf表示各个井点的碳恢复系数。
具体的,在步骤104中,可以按照以下公式计算各个井点的最大有机质孔隙度:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示各个井点的最大有机质孔隙度,单位为%,iTOC表示各个井点的原始有机碳重量百分含量,单位为wt%,TOC表示获取的各个井点的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
具体实施是,可以先按照以下公式确定面孔率:
psurf=a×Ro-b
其中,psurf表示面孔率,单位为%,a和b表示回归系数,Ro表示各个井点的有机质热演化程度。
然后再按照以下公式计算各个井点的观测有机质孔隙度;
其中,Φmin表示观测有机质孔隙度,单位为%,psurf表示面孔率,单位为%,TOC表示获取的各个井点的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
具体的,可以按照以下公式计算各个井点的平均有机质孔隙度:
Φom=(Φmin+Φmax)/2
其中,Φom表示各个井点的平均有机质孔隙度,单位为%,Φmin表示各个井点的观测有机质孔隙度,单位为%,Φmax表示各个井点的最大有机质孔隙度,单位为%。
为了更好地说明本发明,下面结合一个具体实施例对本发明进行说明,然而值得注意的是,该具体实施例仅是为了更好地说明本发明,并不构成对本发明的不当限定。
在本例中提出了一种测算地下有机质孔隙度的方法,该方法具体包括:
步骤1:收集和准备相关数据,其中,该相关数据可以包括:烃源岩密度、有机碳重量百分含量TOC、有机质类型、有机碳密度、有机碳恢复系数与Ro之间的关系曲线等;
步骤2:利用通用盆地模拟软件系统模拟泥页岩有机质热演化史,得到有机质热演化程度Ro;
步骤3:计算原始有机碳重量百分含量。具体的,根据碳恢复系数与Ro之间的关系曲线,查找Ro-pd对应的碳恢复系数Cf,并计算原始有机碳重量百分含量,计算公式为:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示当前有机碳重量百分含量,单位为wt%,Cf表示碳恢复系数,用小数表示。
步骤4:根据物质守恒原理,计算地下泥页岩理论上的最大有机质孔隙度Φmax,原始有机碳重量百分含量与当前有机碳重量百分含量之差等于有机碳转化为烃的量,这部分固体碳转化成非固体的有机质后,原来所占的固体空间已转化成非固体的空间,即变为有机质孔隙。对应的孔隙度可称为理论上最大孔隙度,即未考虑孔隙度被破坏的理想孔隙度。计算公式为:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示各个井点的最大有机质孔隙度,单位为%,iTOC表示各个井点的原始有机碳重量百分含量,单位为wt%,TOC表示获取的各个井点的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
步骤5:利用面孔率计算地面有机质孔隙度Φmin,所谓面孔率是指图片中有机质孔隙面积占有机质面积的百分比,通过观测一组扫描电镜图片可以得到有机质面孔率及对应的Ro值,用回归方法建立两者之间的关系,然后根据面孔率和TOC计算有机质孔隙度。这种方法计算得到的孔隙度可称为观测孔隙度,观测孔隙度在生油窗阶段(Ro<1.2)比地下实际孔隙度小,这主要是因为地下高温高压与地面常温常压的环境差别较大,在地下被液态烃填充的孔隙,在地面可能被沥青充填。具体的,计算地面有机质孔隙度Φmin和面孔率的公式为:
其中,psurf表示面孔率,单位为%,TOC表示获取的有机碳重量百分含量,单位为%,Ro表示有机质热演化程度,单位为%,a和b表示回归系数,与有机质类型有关。
步骤6:计算地下平均有机质孔隙度Φom,具体的,可以将地下泥页岩理论上最大有机质孔隙度和地面观测有机质孔隙度进行加权平均,获得地下平均有机质孔隙度,计算公式为:
Φom=(Φmin+Φmax)/2
其中,Φom表示各个井点的平均有机质孔隙度,单位为%,Φmin表示各个井点的观测有机质孔隙度,单位为%,Φmax表示各个井点的最大有机质孔隙度,单位为%。
然而,值得注意是上,上述计算地下平均有机质孔隙度相当于将最大有机质孔隙度的权值和观测有机质孔隙度的权值都设定为0.5,这仅是一种实现方式,还可以采用其它的权值,例如,将地下泥页岩理论上最大有机质孔隙度的加权系数设定为0.4,将地面观测有机质孔隙度的加权系数设定为0.6。具体的,加权系数可以根据经验或者是实验数据的准确性选择。
下面以位于鄂尔多斯盆地西南部,面积约6.19×104km2,目的层为三叠系延长组长7-3泥页岩(简称T3C73)为例来说明通过本发明中提供的方法测算地下有机质孔隙度的方法,T3C73主体为烃源层,烃源岩厚度一般10~25m,最厚可达30m,优质烃源岩分布范围近5×104km2;有机碳含量主要分布于3%~9%之间;有机母质类型以Ⅱa型干酪根为主,有机质热演化程度Ro在介于生油窗之内,是该区中生界石油的主力油源。
在本例中,提供了一种测算长7-3泥页岩20个井点有机质孔隙度的方法,包括:
步骤1:收集研究区的相关数据,其中,主要包括:如图2所示的20个井点有机碳重量百分含量TOC,如图3所示的有机碳恢复系数与有机质热演化程度Ro之间的关系曲线等,其中,泥页岩的密度为2.5t/m3,有机碳密度为1.25t/m3。
步骤2:采用盆地模拟软件系统模拟泥页岩有机质热演化史,得到如图4所示的20个井点的有机质热演化程度Ro分布图。
步骤3:根据收集到的20个井点的有机碳重量百分含量计算这20个井点的原始有机碳重量百分含量,计算结果如图5所示。
步骤4:根据烃源岩密度、有机碳密度、收集到的有机碳重量百分含量和原始有机碳重量百分含量,计算得到如图6所示的地下泥页岩理论上最大有机质孔隙度。
步骤5:采用回归方法建立面孔率与Ro之间的关系;
步骤6:根据建立的面孔率和Ro之间的关系确定面孔率,根据面孔率和当前有机碳重量百分含量TOC,计算观测有机质孔隙度;
步骤7:将地下泥页岩理论上最大有机质孔隙度和观测有机质孔隙度加权平均,获得平均有机质孔隙度分布。如图6所示,是理论孔隙度(即最大孔隙度)、平均孔隙度和地面孔隙度(即观测孔隙度)的示意图。
相较于现有技术中的几种方式,本发明实施例可以达到以下效果:
1)在本例中,直接采用已有的有机碳恢复系数图版,计算原始有机碳和地下泥页岩理论上最大有机质孔隙度,不需要做实验,计算周期快,计算过程简便。
2)用一组观察数据进行拟合,形成了具有代表性的观测孔隙度计算模型,并用该计算模型来计算地面观测孔隙度,其结果比直接用观测结果更有代表性。
3)将理想状态下的最大有机质孔隙度和观测孔隙度的平均值作为地下孔隙度,结果更合理。
基于同一发明构思,本发明实施例中还提供了一种测算地下有机质孔隙度的装置,如下面的实施例所述。由于测算地下有机质孔隙度的装置解决问题的原理与测算地下有机质孔隙度的方法相似,因此测算地下有机质孔隙度的装置的实施可以参见测算地下有机质孔隙度的方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图7是本发明实施例的测算地下有机质孔隙度的装置的一种结构框图,如图7所示,包括:数据获取单元701、模拟单元702和执行操作单元703,下面对该结构进行说明。
数据获取单元701,用于在研究区中选取多个井点,获取所述研究区的有机碳恢复系数与有机质热演化程度之间的关系曲线;
模拟单元702,用于模拟所述研究区泥页岩的有机质热演化史,得到所述多个井点中各个井点的有机质热演化程度;
执行操作单元703,用于对各个井点执行以下操作:
获取该井点的有机碳重量百分含量、烃源岩密度和有机碳密度;
根据该井点的有机质热演化程度,从所述关系曲线中查找出该井点的碳恢复系数,根据获取的该井点的有机碳重量百分含量和查找到的该井点的碳恢复系数,计算该井点的原始有机碳重量百分含量;
根据获取的该井点的有机碳重量百分含量、原始有机碳重量百分含量、烃源岩密度和有机碳密度,计算该井点的最大有机质孔隙度;
根据该井点的有机质热演化程度确定面孔率,并根据确定的该井点的面孔率和获取的该井点的有机碳重量百分含量,计算得到该井点的观测有机质孔隙度;
对该井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均,计算得到该井点的平均有机质孔隙度。
在一个实施例中,执行操作单元703具体用于按照以下公式计算原始有机碳重量百分含量:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,Cf表示碳恢复系数。
在一个实施例中,执行操作单元703具体用于按照以下公式计算最大有机质孔隙度:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示最大有机质孔隙度,单位为%,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
在一个实施例中,执行操作单元703具体用于按照以下公式确定面孔率:
psurf=a×Ro-b
其中,psurf表示面孔率,单位为%,a和b表示回归系数,Ro表示有机质热演化程度。
在一个实施例中,所述执行操作单元具体用于按照以下公式计算观测有机质孔隙度;
其中,Φmin表示观测有机质孔隙度,单位为%,psurf表示面孔率,单位为%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
在一个实施例中,执行操作单元703具体用于按照以下公式计算平均有机质孔隙度:
Φom=(Φmin+Φmax)/2
其中,Φom表示各平均有机质孔隙度,单位为%,Φmin表示观测有机质孔隙度,单位为%,Φmax表示最大有机质孔隙度,单位为%。
在另外一个实施例中,还提供了一种软件,该软件用于执行上述实施例及优选实施方式中描述的技术方案。
在另外一个实施例中,还提供了一种存储介质,该存储介质中存储有上述软件,该存储介质包括但不限于:光盘、软盘、硬盘、可擦写存储器等。
从以上的描述中,可以看出,本发明实施例实现了如下技术效果:将各个井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均的结果,作为研究区中各个井点的平均有机质孔隙度,从而解决了现有技术中对有机质孔隙度测算结果不准确的技术问题,提高了测算结果的准确性。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种测算地下有机质孔隙度的方法,其特征在于,包括:
在研究区中选取多个井点,获取所述研究区的有机碳恢复系数与有机质热演化程度之间的关系曲线;
模拟所述研究区泥页岩的有机质热演化史,得到所述多个井点中各个井点的有机质热演化程度;
对各个井点执行以下操作:
获取该井点的有机碳重量百分含量、烃源岩密度和有机碳密度;
根据该井点的有机质热演化程度,从所述关系曲线中查找出该井点的碳恢复系数,根据获取的该井点的有机碳重量百分含量和查找到的该井点的碳恢复系数,计算该井点的原始有机碳重量百分含量;
根据获取的该井点的有机碳重量百分含量、原始有机碳重量百分含量、烃源岩密度和有机碳密度,计算该井点的最大有机质孔隙度;
根据该井点的有机质热演化程度确定面孔率,并根据确定的该井点的面孔率和获取的该井点的有机碳重量百分含量,计算得到该井点的观测有机质孔隙度;
对该井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均,计算得到该井点的平均有机质孔隙度。
2.如权利要求1所述的方法,其特征在于,按照以下公式计算原始有机碳重量百分含量:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,Cf表示碳恢复系数。
3.如权利要求1所述的方法,其特征在于,按照以下公式计算最大有机质孔隙度:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示最大有机质孔隙度,单位为%,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
4.如权利要求1所述的方法,其特征在于,按照以下公式确定面孔率:
psurf=a×Ro-b
其中,psurf表示面孔率,单位为%,a和b表示回归系数,Ro表示有机质热演化程度。
5.如权利要求4所述的方法,其特征在于,按照以下公式计算观测有机质孔隙度;
其中,Φmin表示观测有机质孔隙度,单位为%,psurf表示面孔率,单位为%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
6.如权利要求1至5中任一项所述的方法,其特征在于,按照以下公式计算平均有机质孔隙度:
Φom=(Φmin+Φmax)/2
其中,Φom表示各平均有机质孔隙度,单位为%,Φmin表示观测有机质孔隙度,单位为%,Φmax表示最大有机质孔隙度,单位为%。
7.一种测算地下有机质孔隙度的装置,其特征在于,包括:
数据获取单元,用于在研究区中选取多个井点,获取所述研究区的有机碳恢复系数与有机质热演化程度之间的关系曲线;
模拟单元,用于模拟所述研究区泥页岩的有机质热演化史,得到所述多个井点中各个井点的有机质热演化程度;
执行操作单元,用于对各个井点执行以下操作:
获取该井点的有机碳重量百分含量、烃源岩密度和有机碳密度;
根据该井点的有机质热演化程度,从所述关系曲线中查找出该井点的碳恢复系数,根据获取的该井点的有机碳重量百分含量和查找到的该井点的碳恢复系数,计算该井点的原始有机碳重量百分含量;
根据获取的该井点的有机碳重量百分含量、原始有机碳重量百分含量、烃源岩密度和有机碳密度,计算该井点的最大有机质孔隙度;
根据该井点的有机质热演化程度确定面孔率,并根据确定的该井点的面孔率和获取的该井点的有机碳重量百分含量,计算得到该井点的观测有机质孔隙度;
对该井点的最大有机质孔隙度和观测有机质孔隙度进行加权平均,计算得到该井点的平均有机质孔隙度。
8.如权利要求7所述的装置,其特征在于,所述执行操作单元具体用于按照以下公式计算原始有机碳重量百分含量:
iTOC=TOC×Cf
其中,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,Cf表示碳恢复系数。
9.如权利要求7所述的装置,其特征在于,所述执行操作单元具体用于按照以下公式计算最大有机质孔隙度:
Φmax=(iTOC-TOC)×ρrock/ρTOC
其中,Φmax表示最大有机质孔隙度,单位为%,iTOC表示原始有机碳重量百分含量,单位为wt%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
10.如权利要求7所述的装置,其特征在于,所述执行操作单元具体用于按照以下公式确定面孔率:
psurf=a×Ro-b
其中,psurf表示面孔率,单位为%,a和b表示回归系数,Ro表示有机质热演化程度。
11.如权利要求10所述的装置,其特征在于,所述执行操作单元具体用于按照以下公式计算观测有机质孔隙度;
其中,Φmin表示观测有机质孔隙度,单位为%,psurf表示面孔率,单位为%,TOC表示获取的有机碳重量百分含量,单位为wt%,ρrock表示烃源岩密度,单位为t/m3,ρTOC表示有机碳密度,单位为t/m3。
12.如权利要求7至11中任一项所述的装置,其特征在于,所述执行操作单元具体用于按照以下公式计算平均有机质孔隙度:
Φom=(Φmin+Φmax)/2
其中,Φom表示各平均有机质孔隙度,单位为%,Φmin表示观测有机质孔隙度,单位为%,Φmax表示最大有机质孔隙度,单位为%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410443833.5A CN104237966B (zh) | 2014-09-02 | 2014-09-02 | 测算地下有机质孔隙度的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410443833.5A CN104237966B (zh) | 2014-09-02 | 2014-09-02 | 测算地下有机质孔隙度的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104237966A true CN104237966A (zh) | 2014-12-24 |
CN104237966B CN104237966B (zh) | 2017-01-18 |
Family
ID=52226407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410443833.5A Active CN104237966B (zh) | 2014-09-02 | 2014-09-02 | 测算地下有机质孔隙度的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104237966B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104727813A (zh) * | 2015-02-12 | 2015-06-24 | 中国石油大学(华东) | 一种源储共生型地层的孔隙度测量方法 |
CN106223941A (zh) * | 2016-07-20 | 2016-12-14 | 中石化石油工程技术服务有限公司 | 基于测井资料的页岩气储层有机孔隙度确定方法 |
CN106326516A (zh) * | 2015-07-02 | 2017-01-11 | 中国石油天然气股份有限公司 | 烃源岩的排烃函数与排烃指数函数的构建方法 |
CN106370578A (zh) * | 2016-08-30 | 2017-02-01 | 中国石油天然气股份有限公司 | 一种确定页岩有机质孔隙度的方法及系统 |
WO2017024530A1 (zh) * | 2015-08-11 | 2017-02-16 | 深圳朝伟达科技有限公司 | 一种计算烃源岩中有机碳含量的方法 |
CN106680142A (zh) * | 2016-09-22 | 2017-05-17 | 中国矿业大学(北京) | 一种计算岩石中干酪根密度的方法 |
CN107643232A (zh) * | 2016-07-21 | 2018-01-30 | 中国石油化工股份有限公司 | 确定页岩中有机质密度的方法 |
CN108412488A (zh) * | 2018-03-13 | 2018-08-17 | 中石化石油工程技术服务有限公司 | 快速确定页岩气储层有机孔隙度的测井方法 |
CN113959919A (zh) * | 2020-07-21 | 2022-01-21 | 中国石油化工股份有限公司 | 一种页岩层段内有机质孔隙度的计算方法及其应用 |
CN114199911A (zh) * | 2020-08-27 | 2022-03-18 | 中国石油化工股份有限公司 | 泥质碳酸盐岩原始有机质丰度恢复的方法及其应用 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454198B (zh) * | 2013-04-24 | 2015-06-17 | 中国石油大学(华东) | 一种泥页岩有机孔隙度检测方法 |
CN103822866B (zh) * | 2014-03-21 | 2016-04-20 | 中国石油大学(华东) | 一种评价目的层段泥页岩孔隙度方法 |
-
2014
- 2014-09-02 CN CN201410443833.5A patent/CN104237966B/zh active Active
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104727813A (zh) * | 2015-02-12 | 2015-06-24 | 中国石油大学(华东) | 一种源储共生型地层的孔隙度测量方法 |
CN104727813B (zh) * | 2015-02-12 | 2017-08-29 | 中国石油大学(华东) | 一种源储共生型地层的孔隙度测量方法 |
CN106326516A (zh) * | 2015-07-02 | 2017-01-11 | 中国石油天然气股份有限公司 | 烃源岩的排烃函数与排烃指数函数的构建方法 |
WO2017024530A1 (zh) * | 2015-08-11 | 2017-02-16 | 深圳朝伟达科技有限公司 | 一种计算烃源岩中有机碳含量的方法 |
CN106223941B (zh) * | 2016-07-20 | 2019-06-14 | 中石化石油工程技术服务有限公司 | 基于测井资料的页岩气储层有机孔隙度确定方法 |
CN106223941A (zh) * | 2016-07-20 | 2016-12-14 | 中石化石油工程技术服务有限公司 | 基于测井资料的页岩气储层有机孔隙度确定方法 |
CN107643232A (zh) * | 2016-07-21 | 2018-01-30 | 中国石油化工股份有限公司 | 确定页岩中有机质密度的方法 |
CN107643232B (zh) * | 2016-07-21 | 2020-03-17 | 中国石油化工股份有限公司 | 确定页岩中有机质密度的方法 |
CN106370578A (zh) * | 2016-08-30 | 2017-02-01 | 中国石油天然气股份有限公司 | 一种确定页岩有机质孔隙度的方法及系统 |
CN106680142A (zh) * | 2016-09-22 | 2017-05-17 | 中国矿业大学(北京) | 一种计算岩石中干酪根密度的方法 |
CN108412488A (zh) * | 2018-03-13 | 2018-08-17 | 中石化石油工程技术服务有限公司 | 快速确定页岩气储层有机孔隙度的测井方法 |
CN113959919A (zh) * | 2020-07-21 | 2022-01-21 | 中国石油化工股份有限公司 | 一种页岩层段内有机质孔隙度的计算方法及其应用 |
CN113959919B (zh) * | 2020-07-21 | 2024-04-09 | 中国石油化工股份有限公司 | 一种页岩层段内有机质孔隙度的计算方法及其应用 |
CN114199911A (zh) * | 2020-08-27 | 2022-03-18 | 中国石油化工股份有限公司 | 泥质碳酸盐岩原始有机质丰度恢复的方法及其应用 |
CN114199911B (zh) * | 2020-08-27 | 2024-05-24 | 中国石油化工股份有限公司 | 泥质碳酸盐岩原始有机质丰度恢复的方法及其应用 |
Also Published As
Publication number | Publication date |
---|---|
CN104237966B (zh) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104237966A (zh) | 测算地下有机质孔隙度的方法和装置 | |
CN104237965A (zh) | 页岩油资源评价方法和装置 | |
Bai et al. | Distribution characteristics and oil mobility thresholds in lacustrine shale reservoir: Insights from N2 adsorption experiments on samples prior to and following hydrocarbon extraction | |
Herring et al. | Topological persistence for relating microstructure and capillary fluid trapping in sandstones | |
Swenson et al. | A GRACE‐based assessment of interannual groundwater dynamics in the C ommunity L and M odel | |
Partington et al. | Interpreting streamflow generation mechanisms from integrated surface‐subsurface flow models of a riparian wetland and catchment | |
Wise | Spatiotemporal variability of the precipitation dipole transition zone in the western United States | |
Baronas et al. | Mixing as a driver of temporal variations in river hydrochemistry: 2. Major and trace element concentration dynamics in the Andes‐Amazon transition | |
Allen et al. | A tree‐ring based reconstruction of Logan River streamflow, northern Utah | |
Santos et al. | Characterization of natural fracture systems: Analysis of uncertainty effects in linear scanline results | |
Xue et al. | Characterization of the reservoir in Lower Silurian and Lower Cambrian shale of south Sichuan Basin, China | |
Stone et al. | Unsaturated zone hydrostratigraphies: A novel archive of past climates in dryland continental regions | |
Mei et al. | The delivery of dissolved organic carbon from a forested hillslope to a headwater stream in southeastern Pennsylvania, USA | |
CN105138749A (zh) | 一种地质约束下高成熟烃源岩原始有机碳恢复方法 | |
Wang et al. | Anchoring multi‐scale models to micron‐scale imaging of multiphase flow in rocks | |
He et al. | Comparison of regression methods to compute atmospheric pressure and earth tidal coefficients in water level associated with Wenchuan Earthquake of 12 May 2008 | |
CN107092995A (zh) | 地热资源评价方法和装置 | |
Chen et al. | Capillary pressure curve determination based on a 2‐D cross‐section analysis via fractal geometry: a bridge between 2‐D and 3‐D pore structure of porous media | |
CN104950004A (zh) | 一种定量评价古老碳酸盐岩储层中沥青含量的方法 | |
Liao et al. | System reliability analysis of landslides subjected to fluctuation of reservoir water level: a case study in the Three Gorges Reservoir area, China | |
Bo et al. | Step loading compression of ultra-soft soil under radial drainage conditions | |
Becker | Development of efficient multiscale multiphysics models accounting for reversible flow at various subsurface energy storage sites | |
Smith et al. | Relative permeability of Williston Basin CO2 storage targets | |
Abdelhamed et al. | Challenges in Hydrologic‐Land Surface Modeling of Permafrost Signatures—A Canadian Perspective | |
CN105181933A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |