CN112858129B - 岩石储集空间测定方法及装置 - Google Patents
岩石储集空间测定方法及装置 Download PDFInfo
- Publication number
- CN112858129B CN112858129B CN201911178756.4A CN201911178756A CN112858129B CN 112858129 B CN112858129 B CN 112858129B CN 201911178756 A CN201911178756 A CN 201911178756A CN 112858129 B CN112858129 B CN 112858129B
- Authority
- CN
- China
- Prior art keywords
- fine
- characteristic information
- scanning
- rock
- pore
- 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.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 252
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000011148 porous material Substances 0.000 claims abstract description 269
- 230000035515 penetration Effects 0.000 claims abstract description 213
- 239000000945 filler Substances 0.000 claims abstract description 95
- 238000009826 distribution Methods 0.000 claims abstract description 64
- 238000005259 measurement Methods 0.000 claims abstract description 22
- 238000012937 correction Methods 0.000 claims description 68
- 238000004458 analytical method Methods 0.000 claims description 31
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 29
- 229910052753 mercury Inorganic materials 0.000 claims description 29
- 238000000691 measurement method Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 13
- 239000002245 particle Substances 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000003825 pressing Methods 0.000 claims description 7
- 238000012512 characterization method Methods 0.000 description 20
- 238000002591 computed tomography Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 239000007789 gas Substances 0.000 description 6
- 239000001307 helium Substances 0.000 description 6
- 229910052734 helium Inorganic materials 0.000 description 6
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 6
- 239000003350 kerosene Substances 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000005553 drilling Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000012216 screening Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000008094 contradictory effect Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012856 packing Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 241000208467 Macadamia Species 0.000 description 1
- 241000238633 Odonata Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000005510 radiation hardening Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
Landscapes
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Dispersion Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本申请实施例提供一种岩石储集空间测定方法及装置，方法包括：对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征；根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比；对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息；根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息；本申请能够得到精准的目标岩石的孔隙结构特征，进而提升目标岩石储集空间测定的准确度。
Description
技术领域
本申请涉及地质勘探领域，具体涉及一种岩石储集空间测定方法及装置。
背景技术
作为油气储层的重要类型，以砾岩为代表的粗粒碎屑岩成为地质学家研究的重要对象，特别是近期，伴随着准噶尔盆地环玛湖地区砾岩大油区及渤海湾盆地东营凹陷北坡砾岩规模储量区的发现，以砾岩为核心的粗粒沉积体系成为重要的研究领域。
孔隙结构分析是储层评价的基础，重点对孔隙的大小、分布及连通性进行评价，随着近期CT技术在石油工业的快速发展，三维孔隙结构分析成为研究者追求的对象，目前已广泛应用于泥页岩、砂岩等储层评价。在CT进行储集空间三维成像研究时，考虑到粗粒碎屑岩储层样品粒度差异大，而X射线CT扫描过程中无法兼顾样品代表性与分辨率，即样品尺寸与分辨率具有矛盾关系，样品尺寸除以CCD的像素数即得到像素点分辨率，像素点分辨率一般为样品尺寸的一千分之一或者是两千分之一，例如选择对全直径等较大粗粒碎屑岩进行CT扫描时，样品的分辨率只能达到几十微米，难以对其孔隙结构进行全面的表征，而如果选择较小的粗粒碎屑岩储层样品进行扫描，考虑到样品的非均质性，则样品的代表性受到质疑。
发明内容
针对现有技术中的问题，本申请提供一种岩石储集空间测定方法及装置，能够得到精准的目标岩石的孔隙结构特征，进而提升目标岩石储集空间测定的准确度。
为了解决上述问题中的至少一个，本申请提供以下技术方案：
第一方面，本申请提供一种岩石储集空间测定方法，包括：
对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征；
根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比；
对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息；
根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息。
进一步地，在所述确定所述目标岩石的整体孔隙特征信息之后，还包括：
对所述目标岩石进行高压压汞分析，得到所述目标岩石的精细孔隙特征信息；
根据所述整体孔隙特征信息中的孔隙度与预先测量得到的实际孔隙度的数值比较结果，确定所述整体孔隙特征信息的第一修正权重和所述精细孔隙特征信息的第二修正权重；
根据所述第一修正权重和所述第二修正权重分别对所述整体孔隙特征信息和所述精细孔隙特征信息进行加权修正，并将经过所述加权修正后的所述整体孔隙特征信息和所述精细孔隙特征信息进行合并处理，得到所述目标岩石的最终孔隙特征信息。
进一步地，所述根据所述整体孔隙特征信息中的孔隙度与预先测量得到的实际孔隙度的数值比较结果，确定所述整体孔隙特征信息的第一修正权重和所述精细孔隙特征信息的第二修正权重，包括：
根据所述整体孔隙特征信息中的孔隙度在预先测量得到的实际孔隙度中的数值占比，确定所述整体孔隙特征信息的第一修正权重；
根据所述孔隙度与所述实际孔隙度的孔隙度差值，确定所述孔隙度差值在所述实际孔隙度中的数值占比和所述精细孔隙特征信息的第二修正权重。
进一步地，所述确定在所述目标岩石中所述精细穿透扫描区域的区域体积占比，包括：
根据所述初步穿透扫描的结果中所述填充物在所述目标岩石中的第一体积占比和所述精细穿透扫描区域在所述填充物中的第二体积占比，确定所述精细穿透扫描区域在所述目标岩石中的区域体积占比。
进一步地，所述根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域，包括：
将所述填充物中符合预设粒度范围和分布代表性条件的区域确定为所述精细穿透扫描区域。
进一步地，在所述对目标岩石进行初步穿透扫描之前，还包括：
根据所述目标岩石的实际尺寸和预设填充物识别清晰度条件，确定所述初步穿透扫描的扫描工作配置参数。
第二方面，本申请提供一种岩石储集空间测定装置，包括：
初步扫描模块，用于对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征；
精细扫描区域确定模块，用于根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比；
精细扫描模块，用于对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息；
特征信息确定模块，用于根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息。
进一步地，还包括：
压汞分析模块，用于对所述目标岩石进行高压压汞分析，得到所述目标岩石的精细孔隙特征信息；
权重确定模块，用于根据所述整体孔隙特征信息中的孔隙度与预先测量得到的实际孔隙度的数值比较结果，确定所述整体孔隙特征信息的第一修正权重和所述精细孔隙特征信息的第二修正权重；
加权修正模块，用于根据所述第一修正权重和所述第二修正权重分别对所述整体孔隙特征信息和所述精细孔隙特征信息进行加权修正，并将经过所述加权修正后的所述整体孔隙特征信息和所述精细孔隙特征信息进行合并处理，得到所述目标岩石的最终孔隙特征信息。
进一步地，所述权重确定模块包括：
第一修正权重确定单元，用于根据所述整体孔隙特征信息中的孔隙度在预先测量得到的实际孔隙度中的数值占比，确定所述整体孔隙特征信息的第一修正权重；
第二修正权重确定单元，用于根据所述孔隙度与所述实际孔隙度的孔隙度差值，确定所述孔隙度差值在所述实际孔隙度中的数值占比和所述精细孔隙特征信息的第二修正权重。
进一步地，所述精细扫描区域确定模块包括：
区域体积占比确定单元，用于根据所述初步穿透扫描的结果中所述填充物在所述目标岩石中的第一体积占比和所述精细穿透扫描区域在所述填充物中的第二体积占比，确定所述精细穿透扫描区域在所述目标岩石中的区域体积占比。
进一步地，所述精细扫描区域确定模块包括：
区域确定条件判定单元，用于将所述填充物中符合预设粒度范围和分布代表性条件的区域确定为所述精细穿透扫描区域。
进一步地，还包括：
扫描配置单元，用于根据所述目标岩石的实际尺寸和预设填充物识别清晰度条件，确定所述初步穿透扫描的扫描工作配置参数。
第三方面，本申请提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的岩石储集空间测定方法的步骤。
第四方面，本申请提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的岩石储集空间测定方法的步骤。
由上述技术方案可知，本申请提供一种岩石储集空间测定方法及装置，通过对目标岩石进行常规分辨率的初步穿透扫描，以确定在整个目标岩石中填充物的基本构成，例如属性特征和分布特征，再基于填充物的属性特征和分布特征，在填充物中确定至少一个能够表征填充物整体特性的精细穿透扫描区域，同时确定该精细穿透扫描区域在目标岩石中的区域体积占比，然后，对该精细穿透扫描区域进行高分辨率的精细穿透扫描，以得到详细准确的区域孔隙特征信息，结合对应的精细穿透扫描区域在整个目标岩石中的区域体积占比对该区域孔隙特征信息进行比例换算，得到能够表征目标岩石整体孔隙特征的整体孔隙特征信息，由于相较于现有技术中依据常规分辨率进行穿透扫描得到的孔隙结构特征，本申请通过对能够存储油气的填充物中具有代表性的区域进行精细穿透扫描，并根据精细穿透扫描得到的区域孔隙特征信息依比例来表征目标岩石的整体孔隙特征信息，由此确定的目标岩石的孔隙结构特征更加精准，对目标岩石储集空间测定的准确度也得到进一步提升。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例中的岩石储集空间测定方法的流程示意图之一；
图2为本申请实施例中的岩石储集空间测定方法的流程示意图之二；
图3为本申请实施例中的岩石储集空间测定方法的流程示意图之三；
图4为本申请实施例中的岩石储集空间测定装置的结构之一；
图5为本申请实施例中的岩石储集空间测定装置的结构之二；
图6为本申请实施例中的岩石储集空间测定装置的结构之三；
图7为本申请实施例中的岩石储集空间测定装置的结构之四；
图8为本申请实施例中精细穿透扫描区域划分示意图；
图9为本申请实施例中精细穿透扫描区域1的精细穿透扫描结果示意图；
图10为本申请实施例中精细穿透扫描区域2的精细穿透扫描结果示意图；
图11为本申请实施例中精细穿透扫描区域3的精细穿透扫描结果示意图；
图12为本申请实施例中精细穿透扫描区域1的三维孔隙模型示意图；
图13为本申请实施例中精细穿透扫描区域2的三维孔隙模型示意图；
图14为本申请实施例中精细穿透扫描区域3的三维孔隙模型示意图；
图15为本申请实施例中精细穿透扫描区域1的孔径分布直方图；
图16为本申请实施例中精细穿透扫描区域2的孔径分布直方图；
图17为本申请实施例中精细穿透扫描区域3的孔径分布直方图；
图18为本申请实施例中由精细穿透扫描区域1、2和3的区域孔隙特征信息确定的整体孔隙特征信息的示意图；
图19为本申请实施例中压汞分析得到的精细孔隙特征信息示意图；
图20为本申请实施例中整体孔隙特征信息和精细孔隙特征信息合并后得到的最终孔隙特征信息的示意图；
图21为本申请实施例中的电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
考虑到现有技术中通过CT扫描目标岩石时无法兼顾样品代表性与分辨率，进而导致目标岩石尺寸较大时扫描得到的孔隙特征信息不准确的问题，本申请提供一种岩石储集空间测定方法及装置，通过对目标岩石进行常规分辨率的初步穿透扫描，以确定在整个目标岩石中填充物的基本构成，例如属性特征和分布特征，再基于填充物的属性特征和分布特征，在填充物中确定至少一个能够表征填充物整体特性的精细穿透扫描区域，同时确定该精细穿透扫描区域在目标岩石中的区域体积占比，然后，对该精细穿透扫描区域进行高分辨率的精细穿透扫描，以得到详细准确的区域孔隙特征信息，结合对应的精细穿透扫描区域在整个目标岩石中的区域体积占比对该区域孔隙特征信息进行比例换算，得到能够表征目标岩石整体孔隙特征的整体孔隙特征信息，由于相较于现有技术中依据常规分辨率进行穿透扫描得到的孔隙结构特征，本申请通过对能够存储油气的填充物中具有代表性的区域进行精细穿透扫描，并根据精细穿透扫描得到的区域孔隙特征信息依比例来表征目标岩石的整体孔隙特征信息，由此确定的目标岩石的孔隙结构特征更加精准，对目标岩石储集空间测定的准确度也得到进一步提升。
在现有技术中，通常采用CT扫描方式获取粗粒碎屑岩储层储集空间三维孔隙模型，即通过对特定尺寸砾岩样品进行X射线无损扫描，获得三维孔隙模型，进而得到孔隙度、孔径分布直方图等孔隙特征信息。由于砾岩粒度远大于砂岩或者泥岩，因此普遍采用大功率CT扫描，对全直径岩心或者较大的岩心进行CT扫描，主要步骤如下：
(1)制备砾岩样品，选择扫描区域；
(2)待CT设备稳定后，根据扫描区域确定CT设备扫描工作时的配置参数，获取三维孔隙模型；
(3)利用现有三维分析软件对三维孔隙模型进行分析，得到孔隙特征信息。
由实际生产经验可知，上述现有技术无法兼顾样品代表性与分辨率精度，样品尺寸与分辨率是固有矛盾，具体地，样品尺寸除以CCD的像素数即为像素点分辨率，像素点分辨率一般为样品尺寸的一千分之一或者是两千分之一。如选择对全直径等较大粗粒碎屑岩进行CT扫描时，样品的分辨率只能达到几十微米，难以对其孔隙结构进行全面的表征；而如果选择较小的粗粒碎屑岩储层样品进行扫描，考虑到样品的非均质性，则样品的代表性受到质疑。
为了能够得到精准的目标岩石的孔隙结构特征，进而提升目标岩石储集空间测定的准确度，本申请提供一种岩石储集空间测定方法的实施例，参见图1，所述岩石储集空间测定方法具体包含有如下内容：
步骤S101：对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征。
可以理解的是，所述目标岩石可以为粗粒碎屑岩中的砾岩，也可以为泥页岩或砂岩等主要有碎屑和填充物组合而成的岩体，根据实际生产经验可知，在油气勘探领域，油气往往存储于所述填充物的孔隙中，即储集空间。
可选地，可以通过现有穿透式扫描设备和方法对目标岩石进行初步穿透扫描，所述初步穿透扫描的结果可以为所述目标岩石的三维岩石模型，然后通过现有三维数据分析软件(例如Avizo Fire或Dragonfly Pro)对该三维岩石模型进行数据分析，从而确定所述目标岩石中填充物的属性特征和分布特征，例如先识别出目标岩石中的填充物及其具体位置和体积占比(即所述分布特征)，同时再进一步分析得到所述填充物的粒度(即所述属性特征)。
可选地，考虑到粗粒碎屑岩粒度较粗影响CT扫描的效果，因此所述目标岩石可以为全直径样品或直径为1.5英寸的圆柱形样品，砾径较小、胶结较好的小粗粒碎屑岩圆柱形样品直径可为1英寸，但其他粗粒碎屑岩样品要求最小为1.5英寸的圆柱。
可以理解的是，当目标岩石的尺寸较大时，对其进行初步穿透扫描只能得到目标岩石的基本构成特征，难以对其能够存储油气的孔隙特征进行测量采集，因此需要对至少一个能够基本表征目标岩石整体特性的特定区域单独进行精细穿透扫描。
步骤S102：根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比。
可选地，在上述步骤S101中确定目标岩石中的填充物的具体位置和体积占比(即所述分布特征)及其粒度(即所述属性特征)后，可以通过在系统中预设至少一个筛选条件(例如在填充物中占比较高的一粒度范围或一设定的粒度范围)，确定符合该筛选条件的至少一个对应的填充物的区域，并将其确定为精细穿透扫描区域，由于通过现有三维数据分析软件能够获知三维岩石模型的实际体积，因此也能确定该精细穿透扫描区域在目标岩石中的区域体积占比，即该精细穿透扫描区域的体积与目标岩石体积的比值。
在一个具体实例中，参见图8，可以根据预设的三个粒度范围筛选出三个对应的精细穿透扫描区域1、精细穿透扫描区域2和精细穿透扫描区域3，并根据现有三维数据分析软件对填充物整体粒度的分析，确定各精细穿透扫描区域在目标岩石中的区域体积占比，例如精细穿透扫描区域1的区域体积占比为25％，精细穿透扫描区域2的区域体积占比为60％，精细穿透扫描区域3的区域体积占比为15％。
在另一个具体实例中，参见图8，也可以人为自定义三个精细穿透扫描区域，并自定义各精细穿透扫描区域对应的在目标岩石中的区域体积占比，或通过现有三维数据分析软件对各精细穿透扫描区域的属性特征与整个填充物进行对比，进而确定对应的区域体积占比。
步骤S103：对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息。
可选地，在步骤S102确定好至少一个精细穿透扫描区域后，可以利用煤油作为介质，使用机械钻头在各所述精细穿透扫描区域钻取直径为5mm、高度为5mm的圆柱形样品，其中，钻取直径为5mm、高度为5mm的圆柱形样品主要考虑到两点：一方面粗粒碎屑岩砾径主体较大，泥质或凝灰质填充物胶结比较疏松，很难钻取直径更小的样品；另一方面小直径样品可以提高X射线的穿透性，有利于提高分辨率，进而得到该精细穿透扫描区域样品的区域孔隙特征信息。
可选地，可以将精细穿透扫描区域样品置于CT扫描设备中进行360°高分辨率扫描，且像素点分辨率设置在1μm～2μm，扫描张数不低于3000张，根据仪器性能设置扫描电压，一般电压应介于60～80KV，要求单张图像光子数不低于2500，样品透过率应介于25～50％之间，单张曝光时间不低于3秒，主要目的是提高CT表征精度与小孔隙刻画比例。
可以理解的是，参见图9、图10和图11，为经过所述精细穿透扫描后得到的扫描结果，由图中可知此时已经能够扫描得到对应的孔隙结构，参见图12、图13和图14，通过现有三维成像软件得到该精细穿透扫描区域样品的三维孔隙模型，通过现有三维数据分析软件对其进行进一步孔隙特征分析，能够得到对应的区域孔隙特征信息，所述区域孔隙特征信息包括但不限于：孔隙度和孔径分布数据(例如孔径分布直方图，参见图15、图16和图17)。
步骤S104：根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息。
可以理解的是，通过步骤S103取得了目标岩体上具有代表性的至少一个精细穿透扫描区域的精确的区域孔隙特征信息，且已知该精细穿透扫描区域在整个目标岩石中的体积占比，将各个所述区域孔隙特征信息分别依据各自的体积占比(各体积占比之和为1)进行比例换算，并将经过比例换算后的各个区域孔隙特征信息相加，得到该目标岩石的整体孔隙特征信息。
在一个具体实例中，参见图15，在精细穿透扫描区域1的孔径分布直方图中孔隙等效直径大于200μm的孔隙体积比例为69.71％，参见图16，在精细穿透扫描区域2的孔径分布直方图中孔隙等效直径大于200μm的孔隙体积比例为81.89％，参见图17，在精细穿透扫描区域3的孔径分布直方图中孔隙等效直径大于200μm的孔隙体积比例为0，此时在结合上述步骤S102中各精细穿透扫描区域的区域体积占比，将区域孔隙特征信息中的孔径分布直方图中孔隙等效直径大于200μm的孔隙体积比例进行比例换算，并将经过比例换算后的各个区域孔隙特征信息相加，得到该目标岩石的整体孔隙特征信息，例如，参见图18，在目标岩石的整体孔隙特征信息中的孔径分布直方图中孔隙等效直径大于200μm的孔隙体积比例的具体计算公式为：69.71×25％+81.89×60％+0×15％＝66.56。
从上述描述可知，本申请实施例提供的岩石储集空间测定方法，能够通过对目标岩石进行常规分辨率的初步穿透扫描，以确定在整个目标岩石中填充物的基本构成，例如属性特征和分布特征，再基于填充物的属性特征和分布特征，在填充物中确定至少一个能够表征填充物整体特性的精细穿透扫描区域，同时确定该精细穿透扫描区域在目标岩石中的区域体积占比，然后，对该精细穿透扫描区域进行高分辨率的精细穿透扫描，以得到详细准确的区域孔隙特征信息，结合对应的精细穿透扫描区域在整个目标岩石中的区域体积占比对该区域孔隙特征信息进行比例换算，得到能够表征目标岩石整体孔隙特征的整体孔隙特征信息，由于相较于现有技术中依据常规分辨率进行穿透扫描得到的孔隙结构特征，本申请通过对能够存储油气的填充物中具有代表性的区域进行精细穿透扫描，并根据精细穿透扫描得到的区域孔隙特征信息依比例来表征目标岩石的整体孔隙特征信息，由此确定的目标岩石的孔隙结构特征更加精准，对目标岩石储集空间测定的准确度也得到进一步提升。
为了能够对CT扫描采集不到的更小粒度的填充物进行进一步的孔隙特征采集，从而使最终得到的目标岩石的孔隙特征信息更加完善和精准，在本申请的岩石储集空间测定方法的一实施例中，在步骤S104之后，还可以具体包含有通过高压压汞进一步确定精细孔隙特征信息的步骤，参见图2，该步骤具体包含有如下内容：
步骤S201：对所述目标岩石进行高压压汞分析，得到所述目标岩石的精细孔隙特征信息。
可以理解的是，在实际生产应用中，由于CT扫描设备自身硬件参数的限制，其对目标岩石的孔隙结构扫描具有较为明显的局限性，例如，常规CT扫描设备难以获取孔隙等效直径2μm以下的孔隙结构(例如孔吼)，因此，本申请还进一步提供了一种确定孔隙等效直径2μm以下的精细孔隙特征信息，并将该精细孔隙特征信息与步骤S101至步骤S104得到的整体孔隙特征信息进行有效结合，且得到更加完整、更加精确的目标岩石的最终孔隙特征信息的技术方案。
可选地，可以通过现有的压汞仪设备对目标岩石进行高压压汞分析，控制最高进汞压力不低于200MPa，根据高压压汞的分析结果作为现有三维数据分析软件的输入数据，并输出该目标岩石的精细孔隙特征信息，所述精细孔隙特征信息包括但不限于：孔隙度和孔径分布数据(例如孔径分布直方图，参见图19)。
可选地，针对同一目标岩石，进行高压压汞分析和精细穿透扫描的先后顺序本申请并不进行具体限定，但若现进行精细穿透扫描，再进行高压压汞分析，则在进行高压压汞分析前可以对目标岩石进行样品制备处理，具体地，可以以煤油为介质，将目标岩石两端切平，并切除精细穿透扫描区域的钻样部分，剩余部分制备成符合高压压汞分析要求的样品，一般直径与高度均为为1英寸的圆柱样品或者直径与高度均为1.5英寸的圆柱。
步骤S202：根据所述整体孔隙特征信息中的孔隙度与预先测量得到的实际孔隙度的数值比较结果，确定所述整体孔隙特征信息的第一修正权重和所述精细孔隙特征信息的第二修正权重。
可选地，经过精细穿透扫描得到的整体孔隙特征信息虽能反映详细准确的孔隙结构且具有较强整体代表性，但其与目标岩石实际情况仍有难以弥补的微小差距，因此，可以通过类似孔隙度这种定量数据准确表征所述差距，并以此确定第一修正权重对该差距进行弥补，使整体孔隙特征信息更加精准；同理，根据该差距可以确定第二修正权重，以对高压压汞分析得到的精细孔隙特征信息进行修正，使其更加精准。
可选地，可以在精细穿透扫描和高压压汞分析工作开始之前，通过现有孔隙度测定仪测量目标岩石的实际孔隙度，然后与精细穿透扫描得到的整体孔隙特征信息中的孔隙度进行数值比较，根据数值比较结果，确定用于对整体孔隙特征信息进行修正的第一修正权重和用于对精细孔隙特征信息进行修正的第二修正权重。
例如，经孔隙度测定仪测量得到的目标岩石的实际孔隙度为ΦZ＝6.5％，经精细穿透扫描得到的整体孔隙特征信息中的孔隙度为ΦCT＝3％，二者相减得到一差值ΦS＝ΦZ-ΦCT＝3.5％，该差值与实际孔隙度的比值能够准确表征CT扫描设备难以扫描确定的孔隙部分在目标岩石中的占比，而该孔隙部分恰好能够被高压压汞分析得到，由此，所述整体孔隙特征信息的第一修正权重可以为孔隙度ΦCT与实际孔隙度ΦZ的比值所述精细孔隙特征信息的第二修正权重可以为该差值ΦS与实际孔隙度ΦZ的比值/>
步骤S203：根据所述第一修正权重和所述第二修正权重分别对所述整体孔隙特征信息和所述精细孔隙特征信息进行加权修正，并将经过所述加权修正后的所述整体孔隙特征信息和所述精细孔隙特征信息进行合并处理，得到所述目标岩石的最终孔隙特征信息。
可以理解的是，由于精细穿透扫描往往用于确定孔隙等效直径大于2μm的孔隙结构，而高压压汞分析往往用于确定孔隙等效直径小于2μm的孔隙结构，因此可以将二者进行有效结合，以使目标岩石的最终孔隙特征信息更加完善，在进行整体孔隙特征信息与精细孔隙特征信息的结合时，所述整体孔隙特征信息依据第一修正权重进行数值比例换算，所述精细孔隙特征信息依据第二修正权重进行数值比例换算，对于同一信息类型，二者换算后的结果相加即为该信息类型的最终数值。
在一个具体实例中，参见图18，所述整体孔隙特征信息中孔隙等效直径位于1μm～2μm的孔隙体积比例为0，参见图19，所述精细孔隙特征信息中孔隙等效直径位于1μm～2μm的孔隙体积比例为3.99％，且由上述步骤S202可知，所述整体孔隙特征信息的第一修正权重所述精细孔隙特征信息的第二修正权重/>则参见图20，将所述整体孔隙特征信息与所述精细孔隙特征信息有效结合后，在目标岩石的最终孔隙特征信息中孔隙等效直径位于1μm～2μm的孔隙体积比例的计算公式为：/>
为了能够对高压压汞确定的精细孔隙特征信息和精细穿透扫描确定的整体孔隙特征信息进行有效合并，从而使得到目标岩石完整的最终孔隙特征信息，在本申请的岩石储集空间测定方法的一实施例中，在步骤S202中，还可以具体包含有分别确定整体孔隙特征信息和精细孔隙特征信息在合并时各自修正权重的步骤，参见图3，该步骤具体包含有如下内容：
步骤S301：根据所述整体孔隙特征信息中的孔隙度在预先测量得到的实际孔隙度中的数值占比，确定所述整体孔隙特征信息的第一修正权重。
步骤S302：根据所述孔隙度与所述实际孔隙度的孔隙度差值，确定所述孔隙度差值在所述实际孔隙度中的数值占比和所述精细孔隙特征信息的第二修正权重。
例如，经孔隙度测定仪测量得到的目标岩石的实际孔隙度为ΦZ＝6.5％，经精细穿透扫描得到的整体孔隙特征信息中的孔隙度为ΦCT＝3％，二者相减得到一差值ΦS＝ΦZ-ΦCT＝3.5％，该差值与实际孔隙度的比值能够准确表征CT扫描设备难以扫描确定的孔隙部分在目标岩石中的占比，而该孔隙部分恰好能够被高压压汞分析得到，由此，所述整体孔隙特征信息的第一修正权重可以为孔隙度ΦCT与实际孔隙度ΦZ的比值所述精细孔隙特征信息的第二修正权重可以为该差值ΦS与实际孔隙度ΦZ的比值/>
为了能够准确确定各精细穿透扫描区域在目标岩石中的区域体积占比，以使其能够准确表征整个目标岩石的孔隙特征，在本申请的岩石储集空间测定方法的一实施例中，在步骤S102中，还可以具体包含有根据精细穿透扫描区域在填充物中的占比，间接确定其在目标岩石中的占比的步骤，该步骤具体包含有如下内容：
根据所述初步穿透扫描的结果中所述填充物在所述目标岩石中的第一体积占比和所述精细穿透扫描区域在所述填充物中的第二体积占比，确定所述精细穿透扫描区域在所述目标岩石中的区域体积占比。
可选地，为了使各精细穿透扫描区域的区域体积占比更加准确，除可以通过现有三维数据分析软件直接确定各精细穿透扫描区域在目标岩石中的区域体积占比外，还可以通过间接计算方式确定区域体积占比，例如，先通过现有三维数据分析软件确定填充物整体在所述目标岩石中的填充物体积占比PF，在确定了精细穿透扫描区域后，再通过现有三维数据分析软件确定各精细穿透扫描区域在填充物整体中的各自占比K1，K2，......，KN；以此确定各精细穿透扫描区域在所述目标岩石中的区域体积占比Ci＝Ki×PF，i＝1，2，......，N。
为了能够在填充物中准确确定至少一个精细穿透扫描区域，以使其能够准确表征填充物整体的孔隙特征，在本申请的岩石储集空间测定方法的一实施例中，在步骤S102中，还可以具体包含有根据填充物中碎屑的不同粒度和不同分布特征，间接确定其在目标岩石中的占比的步骤，该步骤具体包含有如下内容：
将所述填充物中符合预设粒度范围和分布代表性条件的区域确定为所述精细穿透扫描区域。
可选地，在上述步骤S101中确定目标岩石中的填充物的具体位置和体积占比(即所述分布特征)及其粒度(即所述属性特征)后，可以通过在系统中预设至少一个筛选条件(例如在填充物中占比较高的一粒度范围(即符合所述分布代表性条件)或一设定的粒度范围(即符合所述预设粒度范围))，确定符合该筛选条件的至少一个对应的填充物的区域，并将其确定为精细穿透扫描区域。
在一个具体实例中，参见图8，可以根据预设的三个粒度范围筛选出三个对应的精细穿透扫描区域1、精细穿透扫描区域2和精细穿透扫描区域3，并根据现有三维数据分析软件对填充物整体粒度的分析，确定各精细穿透扫描区域在目标岩石中的区域体积占比，例如精细穿透扫描区域1的区域体积占比为25％，精细穿透扫描区域2的区域体积占比为60％，精细穿透扫描区域3的区域体积占比为15％。
在另一个具体实例中，参见图8，也可以人为自定义三个精细穿透扫描区域，并自定义各精细穿透扫描区域对应的在目标岩石中的区域体积占比，或通过现有三维数据分析软件对各精细穿透扫描区域的属性特征与整个填充物进行对比，进而确定对应的区域体积占比。
为了能够根据扫描结果清晰划分出填充物，以使后续对填充物的体积计算更加准确，在本申请的岩石储集空间测定方法的一实施例中，在步骤S101之前，还可以具体包含有根据对扫描设备的工作参数进行配置的步骤，该步骤具体包含有如下内容：
根据所述目标岩石的实际尺寸和预设填充物识别清晰度条件，确定所述初步穿透扫描的扫描工作配置参数。
可以理解的是，若所述目标岩石的实际尺寸较大，则常规CT扫描设备为了能完整反映目标岩石的整体结构，其所呈现的扫描结果难以得到目标岩石内部精准的结构特征，为了能够清晰扫描出目标岩石内部填充物的结构特征，需对扫描设备的扫描工作配置参数进行动态调节，所述扫描工作配置参数包括但不限于：工作电压、曝光时间和像素点分辨率。
为了能够得到精准的目标岩石的孔隙结构特征，进而提升目标岩石储集空间测定的准确度，本申请提供一种用于实现所述岩石储集空间测定方法的全部或部分内容的岩石储集空间测定装置的实施例，参见图4，所述岩石储集空间测定装置具体包含有如下内容：
初步扫描模块10，用于对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征。
精细扫描区域确定模块20，用于根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比。
精细扫描模块30，用于对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息。
特征信息确定模块40，用于根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息。
从上述描述可知，本申请实施例提供的岩石储集空间测定装置，能够通过对目标岩石进行常规分辨率的初步穿透扫描，以确定在整个目标岩石中填充物的基本构成，例如属性特征和分布特征，再基于填充物的属性特征和分布特征，在填充物中确定至少一个能够表征填充物整体特性的精细穿透扫描区域，同时确定该精细穿透扫描区域在目标岩石中的区域体积占比，然后，对该精细穿透扫描区域进行高分辨率的精细穿透扫描，以得到详细准确的区域孔隙特征信息，结合对应的精细穿透扫描区域在整个目标岩石中的区域体积占比对该区域孔隙特征信息进行比例换算，得到能够表征目标岩石整体孔隙特征的整体孔隙特征信息，由于相较于现有技术中依据常规分辨率进行穿透扫描得到的孔隙结构特征，本申请通过对能够存储油气的填充物中具有代表性的区域进行精细穿透扫描，并根据精细穿透扫描得到的区域孔隙特征信息依比例来表征目标岩石的整体孔隙特征信息，由此确定的目标岩石的孔隙结构特征更加精准，对目标岩石储集空间测定的准确度也得到进一步提升。
为了能够对CT扫描采集不到的更小粒度的填充物进行进一步的孔隙特征采集，从而使最终得到的目标岩石的孔隙特征信息更加完善和精准，在本申请的岩石储集空间测定装置的一实施例中，参见图5，还具体包含有如下内容：
压汞分析模块50，用于对所述目标岩石进行高压压汞分析，得到所述目标岩石的精细孔隙特征信息。
权重确定模块60，用于根据所述整体孔隙特征信息中的孔隙度与预先测量得到的实际孔隙度的数值比较结果，确定所述整体孔隙特征信息的第一修正权重和所述精细孔隙特征信息的第二修正权重。
加权修正模块70，用于根据所述第一修正权重和所述第二修正权重分别对所述整体孔隙特征信息和所述精细孔隙特征信息进行加权修正，并将经过所述加权修正后的所述整体孔隙特征信息和所述精细孔隙特征信息进行合并处理，得到所述目标岩石的最终孔隙特征信息。
为了能够对高压压汞确定的精细孔隙特征信息和精细穿透扫描确定的整体孔隙特征信息进行有效合并，从而使得到目标岩石完整的最终孔隙特征信息，在本申请的岩石储集空间测定装置的一实施例中，参见图6，所述权重确定模块60包括：
第一修正权重确定单元61，用于根据所述整体孔隙特征信息中的孔隙度在预先测量得到的实际孔隙度中的数值占比，确定所述整体孔隙特征信息的第一修正权重。
第二修正权重确定单元62，用于根据所述孔隙度与所述实际孔隙度的孔隙度差值，确定所述孔隙度差值在所述实际孔隙度中的数值占比和所述精细孔隙特征信息的第二修正权重。
为了能够准确确定各精细穿透扫描区域在目标岩石中的区域体积占比，以使其能够准确表征整个目标岩石的孔隙特征，在本申请的岩石储集空间测定装置的一实施例中，参见图7，所述精细扫描区域确定模块20包括：
区域体积占比确定单元21，用于根据所述初步穿透扫描的结果中所述填充物在所述目标岩石中的第一体积占比和所述精细穿透扫描区域在所述填充物中的第二体积占比，确定所述精细穿透扫描区域在所述目标岩石中的区域体积占比。
为了能够在填充物中准确确定至少一个精细穿透扫描区域，以使其能够准确表征填充物整体的孔隙特征，在本申请的岩石储集空间测定装置的一实施例中，参见图7，所述精细扫描区域确定模块20包括：
区域确定条件判定单元22，用于将所述填充物中符合预设粒度范围和分布代表性条件的区域确定为所述精细穿透扫描区域。
为了更进一步说明本方案，本申请还提供一种应用上述岩石储集空间测定装置实现岩石储集空间测定方法的具体应用实例，具体包含有如下内容：
步骤1，选择待测定的粗粒碎屑岩样品，利用煤油和机械钻头，将其制备成圆柱形样品，要求圆柱高度不小于直径，圆柱形样品在旋转过程中，X射线穿透样品的厚度没有明显变化，因此可避免射线硬化，提升CT成像效果；考虑到粗粒碎屑岩粒度较粗，要求样品尽量为全直径样品或直径为1.5英寸的圆柱形样品，粒径较小、胶结较好的小粗粒碎屑岩圆柱形样品直径可为1英寸，但其他粗粒碎屑岩样品要求最小为1.5英寸的圆柱；
步骤2，样品制备完成后，将其置入现有的孔隙度测定仪中，测定样品的氦气孔隙度(ΦZ)；主要目的是了解和认识样品整体的孔隙空间，为后续孔隙结构精细刻画提供总体约束；
步骤3，选择合适的CT设备对样品进行整体扫描，根据样品的尺寸确定扫描参数，保证图像可清晰分辨砾石与填充物；全直径样品，选择医用CT，扫描电压介于150～220KV(不同设备可以略有差异)；1.5英寸或者1英寸直径样品，选择可调节扫描电压的微米CT，扫描电压120～140KV，扫描角度为360°，图像中可对砾石与填充物进行清晰的识别，主要目的是后续二者比例关系的定量计算和精细扫描单元的选取奠定良好的图像基础；
步骤4，扫描结束后，利用现有三维分析软件，计算砾石体积比例(PC)与填充物体积比例(PF)；根据二维图像中填充物的类型与分布特征，确定特征性精细扫描单元(F1，F2，......，FN)，计算每一个精细扫描单元占填充物的体积比例(K1，K2，......，KN)；特征性精细扫描单元主要位于填充物之中，主要的依据是已有研究成果揭示粗粒碎屑岩储集空间主体位于填充物内而非砾石内；此时，每一个精细扫描单元占岩石的体积比例(C1，C2，......，CN)可确定为：
Ci＝Ki×PF，i＝1，2，......，N；
步骤5，确定好精细扫描单元位置后，利用煤油作为介质，使用机械钻头钻取直径为5mm、高度为5mm的样品；钻取5mm的样品主要考虑到两点：一方面粗粒碎屑岩砾径主体较大，泥质或凝灰质填充物胶结比较疏松，很难钻取直径更小的样品；另一方面小直径样品可以提高X射线的穿透性，有利于提高分辨率，这与精细扫描样品需研究孔隙结构的需求相吻合；粗粒碎屑岩精细扫描单元的数量根据样品整体扫描的结果确定，一般一个端面上不超过3个样品点；
步骤6，将精细扫描单元样品置于CT中完成360°高分辨率扫描，要求像素点分辨率应在1μm～2μm，扫描张数不低于3000张，根据仪器性能设置扫描电压，一般电压应介于60～80KV，要求单张图像光子数不低于2500，样品透过率应介于25～50％之间，单张曝光时间不低于3秒；主要目的是提高CT表征精度与小孔隙刻画比例；
步骤7，利用现有三维分析软件完成每一个精细扫描单元的孔隙模型重构，计算得到孔径分布直方图与三维CT孔隙度(ΦCT-1，ΦCT-2，......，ΦCT-N)；
步骤8，根据精细扫描单元占岩石整体的体积比例，将CT三维模型计算的孔隙度与孔径分布进行加权平均，计算得到综合孔径分布图与综合孔隙度(ΦCT)，例如，所述综合孔隙度的计算公式为；
ΦCT＝ΦCT-1×C1+ΦCT-2×C2+……ΦCT-N×CN)；
步骤9，以煤油为介质，将整体样品两端切平，切除精细扫描单元钻样部分，剩余部分制备成符合高压压汞分析要求的样品，一般直径与高度均为为1英寸的圆柱样品或者直径与高度均为1.5英寸的圆柱；
步骤10，利用孔隙度仪测定样品氦气孔隙度(ΦM)，完成高压压汞分析，要求最高进汞压力不低于200MPa，根据高压压汞结果编制孔径直径分布直方图；测定氦气孔隙度ΦM的目的是用于校正整体样品的氦气孔隙度；
步骤11，将CT重构的三维孔隙模型与高压压汞计算的孔隙模型进行融合，得到样品整体的孔隙结构；主要的考虑是两种分析方法具有差异，CT受限于分辨率，只能表征部分孔隙空间，而高压压汞是间接的定量方法，无法提供直观的孔隙空间结构，因此将两种方法进行融合，CT分辨率以上的部分根据CT三维重构结构进行确定，CT分辨率之下的部分则由高压压汞对应部分确定比例关系，这两部分总体的权重根据样品氦气孔隙度与CT计算综合孔隙度的结果确定，具体如下：
假定CT精细扫描可将样品中分辨率以上的孔隙全部识别，那么将样品整体氦气孔隙度(ΦZ)与CT计算的综合孔隙度相减(ΦCT)，得到CT分辨率之下的粗粒碎屑岩孔隙度(ΦS)，进而得到CT计算的综合孔隙度(ΦCT)占总体孔隙度(ΦZ)的比例(KCT)与CT扫描未表征孔隙度(ΦS)占总体孔隙度(ΦZ)的相对比例(KCT-N)，KCT是CT计算的综合孔径分布直方图的权重、KCT-N是压汞分析中位于CT分辨率之下的孔径分布直方图权重，加权后将两图融合得到整体样品的孔径分布。
由上述描述可知，本申请还可以实现如下技术效果：
利用不同分辨率CT的多次扫描与高压压汞联测，可实现对粗粒碎屑岩储集空间三维的直观定量表征。有效填补了粗粒碎屑岩储集空间三维测定方法的空白，提高了结果的准确性，为粗粒碎屑岩储层评价、渗流机理研究与资源量计算提供重要的技术支持，并进一步推动相关实验室基础研究工作的快速发展。
本申请的实施例还提供能够实现上述实施例中的岩石储集空间测定方法中全部步骤的一种电子设备的具体实施方式，参见图21，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；所述通信接口603用于实现岩石储集空间测定装置、在线业务系统、客户端设备以及其他参与机构之间的信息传输；
所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的岩石储集空间测定方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤S101：对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征。
步骤S102：根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比。
步骤S103：对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息。
步骤S104：根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息。
从上述描述可知，本申请实施例提供的电子设备，能够通过对目标岩石进行常规分辨率的初步穿透扫描，以确定在整个目标岩石中填充物的基本构成，例如属性特征和分布特征，再基于填充物的属性特征和分布特征，在填充物中确定至少一个能够表征填充物整体特性的精细穿透扫描区域，同时确定该精细穿透扫描区域在目标岩石中的区域体积占比，然后，对该精细穿透扫描区域进行高分辨率的精细穿透扫描，以得到详细准确的区域孔隙特征信息，结合对应的精细穿透扫描区域在整个目标岩石中的区域体积占比对该区域孔隙特征信息进行比例换算，得到能够表征目标岩石整体孔隙特征的整体孔隙特征信息，由于相较于现有技术中依据常规分辨率进行穿透扫描得到的孔隙结构特征，本申请通过对能够存储油气的填充物中具有代表性的区域进行精细穿透扫描，并根据精细穿透扫描得到的区域孔隙特征信息依比例来表征目标岩石的整体孔隙特征信息，由此确定的目标岩石的孔隙结构特征更加精准，对目标岩石储集空间测定的准确度也得到进一步提升。
本申请的实施例还提供能够实现上述实施例中的岩石储集空间测定方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的岩石储集空间测定方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤S101：对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征。
步骤S102：根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比。
步骤S103：对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息。
步骤S104：根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息。
从上述描述可知，本申请实施例提供的计算机可读存储介质，能够通过对目标岩石进行常规分辨率的初步穿透扫描，以确定在整个目标岩石中填充物的基本构成，例如属性特征和分布特征，再基于填充物的属性特征和分布特征，在填充物中确定至少一个能够表征填充物整体特性的精细穿透扫描区域，同时确定该精细穿透扫描区域在目标岩石中的区域体积占比，然后，对该精细穿透扫描区域进行高分辨率的精细穿透扫描，以得到详细准确的区域孔隙特征信息，结合对应的精细穿透扫描区域在整个目标岩石中的区域体积占比对该区域孔隙特征信息进行比例换算，得到能够表征目标岩石整体孔隙特征的整体孔隙特征信息，由于相较于现有技术中依据常规分辨率进行穿透扫描得到的孔隙结构特征，本申请通过对能够存储油气的填充物中具有代表性的区域进行精细穿透扫描，并根据精细穿透扫描得到的区域孔隙特征信息依比例来表征目标岩石的整体孔隙特征信息，由此确定的目标岩石的孔隙结构特征更加精准，对目标岩石储集空间测定的准确度也得到进一步提升。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (10)
1.一种岩石储集空间测定方法，其特征在于，所述方法包括：
对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征；
根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比；
对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息；
根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息；
对所述目标岩石进行高压压汞分析，得到所述目标岩石的精细孔隙特征信息；
根据所述整体孔隙特征信息中的孔隙度在预先测量得到的实际孔隙度中的数值占比，确定所述整体孔隙特征信息的第一修正权重；
根据所述孔隙度与所述实际孔隙度的孔隙度差值，确定所述孔隙度差值在所述实际孔隙度中的数值占比和所述精细孔隙特征信息的第二修正权重；
根据所述第一修正权重和所述第二修正权重分别对所述整体孔隙特征信息和所述精细孔隙特征信息进行加权修正，并将经过所述加权修正后的所述整体孔隙特征信息和所述精细孔隙特征信息进行合并处理，得到所述目标岩石的最终孔隙特征信息。
2.根据权利要求1所述的岩石储集空间测定方法，其特征在于，所述确定在所述目标岩石中所述精细穿透扫描区域的区域体积占比，包括：
根据所述初步穿透扫描的结果中所述填充物在所述目标岩石中的第一体积占比和所述精细穿透扫描区域在所述填充物中的第二体积占比，确定所述精细穿透扫描区域在所述目标岩石中的区域体积占比。
3.根据权利要求1所述的岩石储集空间测定方法，其特征在于，所述根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域，包括：
将所述填充物中符合预设粒度范围和分布代表性条件的区域确定为所述精细穿透扫描区域。
4.根据权利要求1所述的岩石储集空间测定方法，其特征在于，在所述对目标岩石进行初步穿透扫描之前，还包括：
根据所述目标岩石的实际尺寸和预设填充物识别清晰度条件，确定所述初步穿透扫描的扫描工作配置参数。
5.一种岩石储集空间测定装置，其特征在于，包括：
初步扫描模块，用于对目标岩石进行初步穿透扫描，并根据所述初步穿透扫描的结果确定所述目标岩石中填充物的属性特征和分布特征；
精细扫描区域确定模块，用于根据所述填充物的属性特征和分布特征，确定至少一个精细穿透扫描区域以及在所述目标岩石中所述精细穿透扫描区域的区域体积占比；
精细扫描模块，用于对所述精细穿透扫描区域进行精细穿透扫描，并根据所述精细穿透扫描的结果得到与所述精细穿透扫描区域对应的区域孔隙特征信息；
特征信息确定模块，用于根据所述区域体积占比和所述区域孔隙特征信息，确定所述目标岩石的整体孔隙特征信息；
压汞分析模块，用于对所述目标岩石进行高压压汞分析，得到所述目标岩石的精细孔隙特征信息；
权重确定模块，包括：第一修正权重确定单元，用于根据所述整体孔隙特征信息中的孔隙度在预先测量得到的实际孔隙度中的数值占比，确定所述整体孔隙特征信息的第一修正权重；第二修正权重确定单元，用于根据所述孔隙度与所述实际孔隙度的孔隙度差值，确定所述孔隙度差值在所述实际孔隙度中的数值占比和所述精细孔隙特征信息的第二修正权重；
加权修正模块，用于根据所述第一修正权重和所述第二修正权重分别对所述整体孔隙特征信息和所述精细孔隙特征信息进行加权修正，并将经过所述加权修正后的所述整体孔隙特征信息和所述精细孔隙特征信息进行合并处理，得到所述目标岩石的最终孔隙特征信息。
6.根据权利要求5所述的岩石储集空间测定装置，其特征在于，所述精细扫描区域确定模块包括：
区域体积占比确定单元，用于根据所述初步穿透扫描的结果中所述填充物在所述目标岩石中的第一体积占比和所述精细穿透扫描区域在所述填充物中的第二体积占比，确定所述精细穿透扫描区域在所述目标岩石中的区域体积占比。
7.根据权利要求5所述的岩石储集空间测定装置，其特征在于，所述精细扫描区域确定模块包括：
区域确定条件判定单元，用于将所述填充物中符合预设粒度范围和分布代表性条件的区域确定为所述精细穿透扫描区域。
8.根据权利要求5所述的岩石储集空间测定装置，其特征在于，还包括：
扫描配置单元，用于根据所述目标岩石的实际尺寸和预设填充物识别清晰度条件，确定所述初步穿透扫描的扫描工作配置参数。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述的岩石储集空间测定方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至4任一项所述的岩石储集空间测定方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911178756.4A CN112858129B (zh) | 2019-11-27 | 2019-11-27 | 岩石储集空间测定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911178756.4A CN112858129B (zh) | 2019-11-27 | 2019-11-27 | 岩石储集空间测定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112858129A CN112858129A (zh) | 2021-05-28 |
| CN112858129B true CN112858129B (zh) | 2023-09-26 |
Family
ID=75985851
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911178756.4A Active CN112858129B (zh) | 2019-11-27 | 2019-11-27 | 岩石储集空间测定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112858129B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103822865A (zh) * | 2014-03-20 | 2014-05-28 | 中国石油大学(华东) | 一种高分辨率三维数字岩心建模方法 |
| CN104352246A (zh) * | 2014-12-02 | 2015-02-18 | 东南大学 | 基于可视化的锥束ct感兴趣区域的扫描方法 |
| CN106780744A (zh) * | 2016-12-27 | 2017-05-31 | 中国石油天然气集团公司 | 利用不同分辨率ct图像构建多尺度三维数字岩心的方法 |
| CN106918607A (zh) * | 2017-03-06 | 2017-07-04 | 西南石油大学 | 一种孔隙结构获取方法及装置 |
| CN107240154A (zh) * | 2017-05-02 | 2017-10-10 | 中国石油天然气股份有限公司 | 一种双孔隙三维数字岩心建模方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| MX2011001035A (es) * | 2011-01-27 | 2012-07-27 | Mexicano Inst Petrol | Procedimiento para determinar la porosidad efectiva y total de rocas sedimentarias carbonatadas, y caracterizacion morfologica de sus micro nanoporos. |
-
2019
- 2019-11-27 CN CN201911178756.4A patent/CN112858129B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103822865A (zh) * | 2014-03-20 | 2014-05-28 | 中国石油大学(华东) | 一种高分辨率三维数字岩心建模方法 |
| CN104352246A (zh) * | 2014-12-02 | 2015-02-18 | 东南大学 | 基于可视化的锥束ct感兴趣区域的扫描方法 |
| CN106780744A (zh) * | 2016-12-27 | 2017-05-31 | 中国石油天然气集团公司 | 利用不同分辨率ct图像构建多尺度三维数字岩心的方法 |
| CN106918607A (zh) * | 2017-03-06 | 2017-07-04 | 西南石油大学 | 一种孔隙结构获取方法及装置 |
| CN107240154A (zh) * | 2017-05-02 | 2017-10-10 | 中国石油天然气股份有限公司 | 一种双孔隙三维数字岩心建模方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112858129A (zh) | 2021-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ketcham et al. | Nondestructive high-resolution visualization and measurement of anisotropic effective porosity in complex lithologies using high-resolution X-ray computed tomography | |
| Artioli et al. | X-ray diffraction microtomography (XRD-CT), a novel tool for non-invasive mapping of phase development in cement materials | |
| Alshibli et al. | Quantifying morphology of sands using 3D imaging | |
| Oddershede et al. | Measuring the stress field around an evolving crack in tensile deformed Mg AZ31 using three-dimensional X-ray diffraction | |
| Zhang et al. | Computational investigation on mass diffusivity in Portland cement paste based on X-ray computed microtomography (μCT) image | |
| Ketcham et al. | Three-dimensional measurement of fractures in heterogeneous materials using high-resolution X-ray computed tomography | |
| US8938045B2 (en) | Method of determining reservoir properties and quality with multiple energy X-ray imaging | |
| Semnani et al. | Quantifying the heterogeneity of shale through statistical combination of imaging across scales | |
| Remeysen et al. | Application of microfocus computed tomography in carbonate reservoir characterization: Possibilities and limitations | |
| CN105628726B (zh) | 一种致密砂岩矿物组成的分析方法及系统 | |
| CN105510204A (zh) | 一种基于ct图像的渗透率预测方法 | |
| Chen et al. | The dynamic evolution of compaction bands in highly porous carbonates: the role of local heterogeneity for nucleation and propagation | |
| CN111366521B (zh) | 多尺度确定孔隙度的方法及相关设备 | |
| Wang et al. | Construction of dynamic three-dimensional microstructure for the hydration of cement using 3D image registration | |
| Borges et al. | Computed tomography to estimate the representative elementary area for soil porosity measurements | |
| Chen et al. | A synchrotron-based local computed tomography combined with data-constrained modelling approach for quantitative analysis of anthracite coal microstructure | |
| Pal et al. | Porosity estimation by digital image analysis | |
| CN112858129B (zh) | 岩石储集空间测定方法及装置 | |
| Miller et al. | A scale-consistent method for imaging porosity and micrite in dual-porosity carbonate rocks | |
| Halisch et al. | Pore-scale tomography and imaging: applications, techniques and recommended practice | |
| Yang et al. | X-ray dark-field contrast imaging of water transport during hydration and drying of early-age cement-based materials | |
| Qiu et al. | 3D microstructure simulation of reactive aggregate in concrete from 2D images as the basis for ASR simulation | |
| Sardini et al. | Mapping and quantifying the clay aggregate microporosity in medium-to coarse-grained sandstones using the 14C-PMMA method | |
| WO2019151889A1 (en) | A method for determining a three-dimensional spatial distribution of porosity in a sample of a heterogeneous porous medium | |
| Abutaha et al. | Evaluation of 3D small-scale lithological heterogeneities and pore distribution of the Boda Claystone Formation using X-ray computed tomography images (CT) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |