CN111307652A - 吸附气量的确定方法及装置 - Google Patents
吸附气量的确定方法及装置 Download PDFInfo
- Publication number
- CN111307652A CN111307652A CN201811509638.2A CN201811509638A CN111307652A CN 111307652 A CN111307652 A CN 111307652A CN 201811509638 A CN201811509638 A CN 201811509638A CN 111307652 A CN111307652 A CN 111307652A
- Authority
- CN
- China
- Prior art keywords
- geological
- parameter values
- function
- adsorbed gas
- gas amount
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000001179 sorption measurement Methods 0.000 title claims description 84
- 239000011435 rock Substances 0.000 claims abstract description 122
- 230000015572 biosynthetic process Effects 0.000 claims description 80
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 13
- 229910052799 carbon Inorganic materials 0.000 claims description 13
- 230000006870 function Effects 0.000 abstract description 206
- 239000007789 gas Substances 0.000 description 221
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000015654 memory Effects 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000003795 desorption Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000594592 Lanugo Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000013547 langmuir monolayer Substances 0.000 description 1
- 230000007246 mechanism Effects 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
- 239000011146 organic particle Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N7/00—Analysing materials by measuring the pressure or volume of a gas or vapour
- G01N7/02—Analysing materials by measuring the pressure or volume of a gas or vapour by absorption, adsorption, or combustion of components and measurement of the change in pressure or volume of the remainder
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种吸附气量的确定方法及装置,属于吸附气技术领域。所述方法包括:基于岩层中第一样本岩心的多个地层特征的参数值、多个地质特征的参数值以及吸附气量,获取该岩层的地质特征函数和吸附气量函数,并基于第二样本岩心的各参数值对上述两个函数进行验证,基于验证通过的两个函数,以及该岩层的待测岩心的各参数值来确定吸附气量。由于上述两个函数是通过第一样本岩心的多个地层特征的参数值、多个地质特征的参数值得到的,所以这两个函数中包含了地层特征和地质特征的参数值对吸附气量的影响,能较全面、准确的反映不同因素对该吸附气量的影响程度,所以基于上述两个函数所得的待测岩心的吸附气量的准确度较高。
Description
技术领域
本发明涉及吸附气技术领域,特别涉及一种吸附气量的确定方法及装置。
背景技术
吸附气是指以吸附状态保存在有机质颗粒表面的气体,在温度一定的条件下,该吸附气的含量与地层压力等因素有关。在页岩气资源的勘探开发过程中,页岩作为一种致密的非常规的天然气储层,其微纳尺度的孔隙发育程度较高,发育的微观孔隙使得吸附解吸效应成为影响页岩气储量计算、采出机理研究、开发方案设计和采收率预测的关键因素,而研究页岩的吸附气量就是研究吸附解吸效应的必要步骤。
目前常用的吸附气量的确定方法是等温吸附实验法,该方法包括:在一定温度条件下,将一定粒度的样本岩心置于密封容器中;当该容器中的气体分别在不同压力条件下达到吸附平衡时,测量该样本岩心所吸附的甲烷等试验气体的体积;根据朗格缪尔单分子层吸附理论,计算出表征该样本岩心对甲烷等试验气体吸附特性的兰氏体积、兰氏压力以及等温吸附曲线;基于待测岩心的压力值来确定待测岩心的吸附气量。
上述方法获取到的等温吸附曲线不能准确的反映该待测岩心的吸附特性,所以根据该等温吸附曲线所得的待测岩心的吸附气量也存在偏差,准确度较低。
发明内容
本发明实施例提供了一种吸附气量的确定方法及装置,能够解决目前常用的方法获取到的等温吸附曲线不能准确的反映该待测岩心的吸附特性,所以根据该等温吸附曲线所得的待测岩心的吸附气量也存在偏差,准确度较低的问题。所述技术方案如下:
一方面,提供了一种吸附气量的确定方法,该方法包括:
基于岩层中各第一样本岩心的多个地层特征的参数值、该各第一样本岩心的吸附气量以及吸附气量函数,获取各第一样本岩心的第一地质特征的参数值,该多个地层特征包括地层压力、兰氏体积和兰氏压力,该第一地质特征的参数值用于表示地质特征对该吸附气量的影响程度,该吸附气量函数为该多个地层特征的参数值、该第一地质特征的参数值以及该吸附气量之间的函数关系;
基于该第一地质特征的参数值和该各第一样本岩心的多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数,该多个第二地质特征包括有机碳含量和孔隙度;
基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度;
当该准确度高于预设阈值时,基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量。
在一种可能实现方式中,该吸附气量函数的获取过程包括:
建立初始吸附气量函数,该初始吸附气量函数为该多个地层特征的参数值与吸附气量之间的函数关系;
基于该初始吸附气量函数,以及该第一地质特征的参数值与该吸附气量之间的函数关系,获取吸附气量函数。
在一种可能实现方式中,该基于该第一地质特征的参数值和该各第一样本岩心的多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数包括:
建立初始地质特征函数;
基于该初始地质特征函数、该第一地质特征的参数值和该多个第二地质特征的参数值,获取该地质特征函数,该地质特征函数为该多个第二地质特征的参数值、多个系数值与该第一特征的参数值之间的函数关系。
在一种可能实现方式中,该基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度包括:
基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个第二样本岩心的理论吸附气量;
获取该各第二样本岩心的实际吸附气量与该各第二样本岩心的理论吸附气量之间的差值;
获取多个差值与该各第二样本岩心的实际吸附气量之间的多个比值;
获取该多个比值的绝对值的平均数,将100%与该平均数之差获取为该准确度。
在一种可能实现方式中,该基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量包括:
基于该每个待测岩心的多个第二地质特征的参数值以及该地质特征函数,获取该每个待测岩心的第一地质特征的参数值;
基于该每个待测岩心的多个地层特征的参数值、该每个待测岩心的第一地质特征的参数值以及该吸附气量函数,获取该每个待测岩心的吸附气量。
一方面,提供了一种吸附气量的确定装置,该装置包括:
第一地质特征的参数值获取模块,用于基于岩层中各第一样本岩心的多个地层特征的参数值、该各第一样本岩心的吸附气量以及吸附气量函数,获取各第一样本岩心的第一地质特征的参数值,该多个地层特征包括地层压力、兰氏体积和兰氏压力,该第一地质特征的参数值用于表示地质特征对该吸附气量的影响程度,该吸附气量函数为该多个地层特征的参数值、该第一地质特征的参数值以及该吸附气量之间的函数关系;
地质特征函数获取模块,用于基于该第一地质特征的参数值和该各第一样本岩心的多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数,该多个第二地质特征包括有机碳含量和孔隙度;
准确度获取模块,用于基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度;
吸附气量获取模块,用于当该准确度高于预设阈值时,基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量。
在一种可能实现方式中,该装置还包括吸附气量函数获取模块,该吸附气量函数获取模块用于:
建立初始吸附气量函数,该初始吸附气量函数为该多个地层特征的参数值与吸附气量之间的函数关系;
基于该初始吸附气量函数,以及该第一地质特征的参数值与该吸附气量之间的函数关系,获取吸附气量函数。
在一种可能实现方式中,该地质特征函数获取模块用于:
建立初始地质特征函数;
基于该初始地质特征函数、该第一地质特征的参数值和该多个第二地质特征的参数值,获取该地质特征函数,该地质特征函数为该多个第二地质特征的参数值、多个系数值与该第一特征的参数值之间的函数关系。
在一种可能实现方式中,该准确度获取模块用于:
基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个第二样本岩心的理论吸附气量;
获取该各第二样本岩心的实际吸附气量与该各第二样本岩心的理论吸附气量之间的差值;
获取多个差值与该各第二样本岩心的实际吸附气量之间的多个比值;
获取该多个比值的绝对值的平均数,将100%与该平均数之差获取为该准确度。
在一种可能实现方式中,该吸附气量获取模块用于:
基于该每个待测岩心的多个第二地质特征的参数值以及该地质特征函数,获取该每个待测岩心的第一地质特征的参数值;
基于该每个待测岩心的多个地层特征的参数值、该每个待测岩心的第一地质特征的参数值以及该吸附气量函数,获取该每个待测岩心的吸附气量。
本发明实施例提供的技术方案带来的有益效果是:
基于岩层中第一样本岩心的多个地层特征的参数值、多个地质特征的参数值以及吸附气量,获取该岩层的地质特征函数和吸附气量函数,并基于第二样本岩心的各参数值对上述两个函数进行验证,基于验证通过的两个函数,以及该岩层的待测岩心的各参数值来确定吸附气量。由于上述两个函数是通过第一样本岩心的多个地层特征的参数值、多个地质特征的参数值得到的,所以这两个函数中包含了地层特征和地质特征的参数值对吸附气量的影响,能较全面、准确的反映不同因素对该吸附气量的影响程度,所以基于上述两个函数所得的待测岩心的吸附气量的准确度较高。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种吸附气量的确定方法的流程图;
图2是本发明实施例提供的一种吸附气量的确定方法的流程图;
图3是本发明实施例提供的一种吸附气量的确定装置的结构示意图;
图4是本发明实施例提供的一种计算机设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种吸附气量的确定方法的流程图。该方法可以应用于计算机设备,参见图1,该实施例包括:
101、基于岩层中各第一样本岩心的多个地层特征的参数值、该各第一样本岩心的吸附气量以及吸附气量函数,获取各第一样本岩心的第一地质特征的参数值。
其中,该多个地层特征包括地层压力、兰氏体积和兰氏压力,该第一地质特征的参数值用于表示地质特征对该吸附气量的影响程度,该吸附气量函数为该多个地层特征的参数值、该第一地质特征的参数值以及该吸附气量之间的函数关系。
102、基于该第一地质特征的参数值和该各第一样本岩心的多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数。
其中,该多个第二地质特征包括有机碳含量和孔隙度。
103、基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度。
104、当该准确度高于预设阈值时,基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量。
本发明实施例提供的方法,基于岩层中第一样本岩心的多个地层特征的参数值、多个地质特征的参数值以及吸附气量,获取该岩层的地质特征函数和吸附气量函数,并基于第二样本岩心的各参数值对上述两个函数进行验证,基于验证通过的两个函数,以及该岩层的待测岩心的各参数值来确定吸附气量。由于上述两个函数是通过第一样本岩心的多个地层特征的参数值、多个地质特征的参数值得到的,所以这两个函数中包含了地层特征和地质特征的参数值对吸附气量的影响,能较全面、准确的反映不同因素对该吸附气量的影响程度,所以基于上述两个函数所得的待测岩心的吸附气量的准确度较高。
在一种可能实现方式中,该吸附气量函数的获取过程包括:建立初始吸附气量函数,该初始吸附气量函数为该多个地层特征的参数值与吸附气量之间的函数关系;基于该初始吸附气量函数,以及该第一地质特征的参数值与该吸附气量之间的函数关系,获取吸附气量函数。
在一种可能实现方式中,该基于该第一地质特征的参数值和该各第一样本岩心的多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数包括:建立初始地质特征函数;基于该初始地质特征函数、该第一地质特征的参数值和该多个第二地质特征的参数值,获取该地质特征函数,该地质特征函数为该多个第二地质特征的参数值、多个系数值与该第一特征的参数值之间的函数关系。
在一种可能实现方式中,该基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度包括:基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个第二样本岩心的理论吸附气量;获取该各第二样本岩心的实际吸附气量与该各第二样本岩心的理论吸附气量之间的差值;获取多个差值与该各第二样本岩心的实际吸附气量之间的多个比值;获取该多个比值的绝对值的平均数,将100%与该平均数之差获取为该准确度。
在一种可能实现方式中,该基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量包括:基于该每个待测岩心的多个第二地质特征的参数值以及该地质特征函数,获取该每个待测岩心的第一地质特征的参数值;基于该每个待测岩心的多个地层特征的参数值、该每个待测岩心的第一地质特征的参数值以及该吸附气量函数,获取该每个待测岩心的吸附气量。
图2是本发明实施例提供的一种吸附气量的确定方法的流程图。该方法可以应用于计算机设备,参见图2,该实施例包括:
201、获取岩层中各第一样本岩心的多个地层特征的参数值和吸附气量。
其中,多个地层特征包括地层压力、兰氏体积和兰氏压力。地层压力是指作用在单位面积的岩石孔隙内的流体上的压力,单位为MPa;兰氏体积又称为朗缪尔体积,是指岩层中各样本岩心在平衡水条件下测得的最大吸附气量,单位为m3/t;兰氏压力是指当岩层中各样本岩心的吸附量达到兰氏体积一半时,该样本岩心的压力,单位为MPa,该多个地层特征的参数值均能够影响吸附气量。以该岩层中7个不同深度位置的第一样本岩心为例,各第一样本岩心的地层压力、兰氏体积、兰氏压力和吸附气量可以是如表1所示的数据。
表1
202、建立初始吸附气量函数,该初始吸附气量函数为该多个地层特征的参数值与吸附气量之间的函数关系。
岩层中通常含有多种不同性质的岩石,不同性质的岩石的兰氏体积和兰氏压力也不相同,该初始吸附气量函数包含了该岩层中多种不同性质的岩石对应的吸附气量。例如,根据单一岩石吸附气量函数(1)得到n种不同性质的岩石的总吸附气量函数(2),该函数(2)即为该初始吸附气量函数。
其中,Vads为吸附气量,单位为m3/t;VL为兰氏体积,单位为m3/t;p为地层压力,单位为MPa;pL为兰氏压力,单位为MPa;n表示该岩层中岩石的种类总数;i为该岩层中任一种岩石的标识。
203、基于该初始吸附气量函数,以及该第一地质特征的参数值与该吸附气量之间的函数关系,获取吸附气量函数。
其中,该第一地质特征的参数值用于表示地质特征对该吸附气量的影响程度,该吸附气量函数为该多个地层特征的参数值、该第一地质特征的参数值以及该吸附气量之间的函数关系。
由于岩层的吸附气量的影响因素较多,包括多个地层特征的参数值和多个地质特征的参数值。为了使该吸附气量函数包括不同的影响因素对吸附气量的影响,通过引入多个初始参数的方式对初始吸附气量函数进行变换,例如,可以引入与地层特征相关的初始参数a和c,以及与地质特征相关的初始参数b,根据该参数a、b和c与吸附气量以及与多个地层特征的参数值、多个地质特征的参数值之间的函数关系,将函数(2)变换为函数(9),在一种可能实现方式中,该具体过程可以为下述步骤:
第一步,在函数(2)中引入与地层特征相关的初始参数a,得到函数(3);
第二步,在函数(3)中引入与地质特征相关的初始参数b以及与地层特征相关的初始参数c,也即是,将函数(4)和函数(5)输入函数(3),得到函数(6);
第三步,将n种不同性质的岩石的总吸附气量函数变换为用该岩层的吸附气量函数来表示的形式,也即是,将函数(7)输入函数(6),得到函数(8);
第四步,根据函数(8),当p趋近于无穷大时,可得Vads=a,根据兰氏体积的定义:兰氏体积是岩层中各第一样本岩石在平衡水条件下测得的最大吸附气量,得到a=VL;基于函数(8)与函数(1)的对应关系,得到c=pL;进而得到引入与地质特征相关的初始参数b后的吸附气量函数(9),其中,该与地质特征相关的初始参数b的值即为第一地质特征的参数值。
在上述函数(3)至函数(9)中,Vads为吸附气量,单位为m3/t;VL为兰氏体积,单位为m3/t;p为地层压力,单位为MPa;pL为兰氏压力,单位为MPa;n表示岩层中岩石的种类总数;i为岩层中任一种岩石的标识。
204、基于岩层中各第一样本岩心的多个地层特征的参数值、该各第一样本岩心的吸附气量以及吸附气量函数,获取该各第一样本岩心的第一地质特征的参数值。
具体地,将该岩层中各第一样本岩心的多个地层特征的参数值、该各第一样本岩心的吸附气量输入该吸附气量函数,分别得到每个样本岩心的第一地质特征的参数值。例如,将步骤201的表1中的数据输入吸附气量函数(9)中,可以得到与表1中的每一个样本岩心对应的第一地质特征的参数值,以便基于该第一地质特征的参数值,获取该第一地质特征的参数值与多个第二地质特征的参数值之间的函数关系。
205、获取该各第一样本岩心的多个第二地质特征的参数值。
其中,该多个第二地质特征包括有机碳含量和孔隙度。有机碳含量是指水体或岩石中溶解性和悬浮性有机物含碳的总量,无因次;孔隙度是指样本岩心中所有孔隙空间体积之和与该样本岩心体积的比值,无因次。同样以该岩层中的该7个不同深度位置的第一样本岩心为例,各第一样本岩心的有机碳含量和孔隙度可以是如表2所示的数据。
表2
深度(m) | 有机碳含量(%) | 孔隙度(%) |
3860 | 0.66 | 4.55 |
3866 | 0.88 | 1.17 |
3870 | 1.71 | 3.68 |
3875 | 2.78 | 4.66 |
3879 | 3.37 | 3.32 |
3882 | 1.59 | 2.12 |
3886 | 2.50 | 2.90 |
3893 | 2.34 | 0.74 |
3316 | 1.65 | 5.39 |
3336 | 3.71 | 5.56 |
3343 | 2.66 | 6.44 |
3346 | 2.73 | 4.46 |
3348 | 2.57 | 5.48 |
206、基于该各第一样本岩心的第一地质特征的参数值和多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数。
在一种可能实现方式中,该步骤206可以通过下述步骤实现:
第一步,建立初始地质特征函数。
该初始地质特征函数能够反映该第一地质特征的参数值与多个第二地质特征的参数值之间的函数关系,例如,用不同的字母d、e、f和g表示多个系数值,该初始地质特征函数(10)可以是:
第二步,基于该初始地质特征函数、该各第一样本岩心的第一地质特征的参数值和多个第二地质特征的参数值,获取该地质特征函数。
其中,该地质特征函数为该多个第二地质特征的参数值、多个系数值与该第一特征的参数值之间的函数关系。
具体地,根据上述步骤204中获取到的第一地质特征的参数值和步骤205中获取的多个第二地质特征的参数值,拟合得到该地质特征函数,以上述步骤204中的第一地质特征的参数值和步骤205中表2中的数据为例,所得的地质特征函数(11)可以是:
该地质特征函数(11)以及上述步骤203中获取的吸附气量函数(9),能够反映该岩层的吸附特性,由于上述两个函数(9)和函数(11)是基于岩层中样本岩心的实际数据拟合得到的,包含了各地层特征和各地质特征的参数值对吸附气量的影响,能较全面、准确的反映各种不同因素对该吸附气量的影响程度,所以基于上述两个函数所得的待测岩心的吸附气量的准确度较高。
207、获取岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量。
在岩层中,获取与该第一样本岩心不同的第二样本岩心,进而获取该第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,用于对基于第一样本岩心获取的地质特征函数和该吸附气量函数的准确度进行验证。
208、基于岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度。
在一种可能实现方式中,该步骤208可以通过下述步骤实现:
第一步,基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个第二样本岩心的理论吸附气量。
以上述步骤中获取的地质特征函数(11)和吸附气量函数(9)为例,将该岩层中各第二样本岩心的多个第二地质特征的参数值输入函数(11)中,得到各第二样本岩心的第一地质特征的参数值,将该各第二样本岩心的第一地质特征的参数值和各第二样本岩心的多个地层特征的参数值输入函数(9)中,得到该各第二样本岩心的理论吸附气量。
第二步,获取该各第二样本岩心的实际吸附气量与理论吸附气量之间的差值;获取多个差值与该多个实际吸附气量之间的多个比值;获取该多个比值的绝对值的平均数,将100%与该平均数之差获取为该准确度。
例如,基于上述步骤207中的岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及吸附气量数据对上述步骤203中获取的吸附气量函数(9)以及上述步骤206中获取的地质特征函数(11)进行验证,该准确度为92%。
当然,该准确度还可以通过其他方式获取,本实施例对此不作限定。
第三步,当该准确度高于预设阈值时,确定该地质特征函数和该吸附气量函数符合预设要求。
如果该准确度高于该预设阈值,则该地质特征函数和该吸附气量函数的准确度符合预设要求,能够在后续步骤中用于确定待测岩心的吸附气量;如果该准确度低于或等于该预设阈值,则该地质特征函数和该吸附气量函数的准确度不符合预设要求,不能在后续步骤中用于确定待测岩心的吸附气量。
例如,该预设阈值可以是90%,若以第三步中所得的准确度为92%为例,该准确度高于预设阈值,也即是该地质特征函数(11)和该吸附气量函数(9)符合预设要求。
209、获取该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值。
基于上述参数值、该地质特征函数以及该吸附气量函数,能够获取该岩层中每个待测岩心的吸附气量。
210、当该准确度高于预设阈值时,基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量。
在一种可能实现方式中,该步骤210可以通过下述步骤实现:
基于该每个待测岩心的多个第二地质特征的参数值以及该地质特征函数,获取该每个待测岩心的第一地质特征的参数值;基于该每个待测岩心的多个地层特征的参数值、该每个待测岩心的第一地质特征的参数值以及该吸附气量函数,获取该每个待测岩心的吸附气量。
该步骤的具体过程与上述步骤208中计算第二样本岩心的理论吸附气量的步骤相似,在此不再赘述。
上述所有可选技术方案,可以采用任意结合形成本公开的可选实施例,在此不再一一赘述。
本发明实施例提供的方法,基于岩层中第一样本岩心的多个地层特征的参数值、多个地质特征的参数值以及吸附气量,获取该岩层的地质特征函数和吸附气量函数,并基于第二样本岩心的各参数值对上述两个函数进行验证,基于验证通过的两个函数,以及该岩层的待测岩心的各参数值来确定吸附气量。由于上述两个函数是通过第一样本岩心的多个地层特征的参数值、多个地质特征的参数值得到的,所以这两个函数中包含了地层特征和地质特征的参数值对吸附气量的影响,能较全面、准确的反映不同因素对该吸附气量的影响程度,所以基于上述两个函数所得的待测岩心的吸附气量的准确度较高。进一步地,通过第一地质特征的参数值这一中间量,来获取多个第二地质特征的参数值与吸附气量之间的两个函数,简化了计算步骤。
图3是本发明实施例提供的一种吸附气量的确定装置的结构示意图,该装置包括:
第一地质特征的参数值获取模块301,用于基于岩层中各第一样本岩心的多个地层特征的参数值、该各第一样本岩心的吸附气量以及吸附气量函数,获取各第一样本岩心的第一地质特征的参数值,该多个地层特征包括地层压力、兰氏体积和兰氏压力,该第一地质特征的参数值用于表示地质特征对该吸附气量的影响程度,该吸附气量函数为该多个地层特征的参数值、该第一地质特征的参数值以及该吸附气量之间的函数关系。
地质特征函数获取模块302,用于基于该第一地质特征的参数值和该各第一样本岩心的多个第二地质特征的参数值,获取该第一地质特征的参数值与该多个第二地质特征的参数值之间的地质特征函数,该多个第二地质特征包括有机碳含量和孔隙度。
准确度获取模块303,用于基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取该地质特征函数和该吸附气量函数的准确度。
吸附气量获取模块304,用于当该准确度大于预设阈值时,基于该岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个待测岩心的吸附气量。
在一种可能实现方式中,该装置还包括吸附气量函数获取模块,该吸附气量函数获取模块用于:建立初始吸附气量函数,该初始吸附气量函数为该多个地层特征的参数值与吸附气量之间的函数关系;基于该初始吸附气量函数,以及该第一地质特征的参数值与该吸附气量之间的函数关系,获取吸附气量函数。
在一种可能实现方式中,该地质特征函数获取模块302用于:建立初始地质特征函数;基于该初始地质特征函数、该第一地质特征的参数值和该多个第二地质特征的参数值,获取该地质特征函数,该地质特征函数为该多个第二地质特征的参数值、多个系数值与该第一特征的参数值之间的函数关系。
在一种可能实现方式中,该准确度获取模块303用于:基于该岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、该地质特征函数以及该吸附气量函数,获取该岩层中每个第二样本岩心的理论吸附气量;获取该各第二样本岩心的实际吸附气量与该各第二样本岩心的理论吸附气量之间的差值;获取多个差值与该各第二样本岩心的实际吸附气量之间的多个比值;获取该多个比值的绝对值的平均数,将100%与该平均数之差获取为该准确度。
在一种可能实现方式中,该吸附气量获取模块304用于:基于该每个待测岩心的多个第二地质特征的参数值以及该地质特征函数,获取该每个待测岩心的第一地质特征的参数值;基于该每个待测岩心的多个地层特征的参数值、该每个待测岩心的第一地质特征的参数值以及该吸附气量函数,获取该每个待测岩心的吸附气量。
需要说明的是:上述实施例提供的吸附气量的确定装置在吸附气量的确定时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的吸附气量的确定装置与吸附气量的确定方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
本发明实施例提供的装置,基于岩层中第一样本岩心的多个地层特征的参数值、多个地质特征的参数值以及吸附气量,获取该岩层的地质特征函数和吸附气量函数,并基于第二样本岩心的各参数值对上述两个函数进行验证,基于验证通过的两个函数,以及该岩层的待测岩心的各参数值来确定吸附气量。由于上述两个函数是通过第一样本岩心的多个地层特征的参数值、多个地质特征的参数值得到的,所以这两个函数中包含了地层特征和地质特征的参数值对吸附气量的影响,能较全面、准确的反映不同因素对该吸附气量的影响程度,所以基于上述两个函数所得的待测岩心的吸附气量的准确度较高。
图4是本发明实施例提供的一种计算机设备的结构示意图,该计算机设备400可因配置或性能不同而产生比较大的差异,可以包括一个或一个以上处理器(centralprocessing units,CPU)401和一个或一个以上的存储器402,其中,所述存储器402中存储有至少一条指令,所述至少一条指令由所述处理器401加载并执行以实现上述各个方法实施例提供的方法。当然,该计算机设备还可以具有有线或无线网络接口、键盘以及输入输出接口等部件,以便进行输入输出,该计算机设备还可以包括其他用于实现设备功能的部件,在此不做赘述。
在示例性实施例中,还提供了一种计算机可读存储介质,例如包括指令的存储器,上述指令可由计算机设备中的处理器执行以完成上述实施例中吸附气量的确定方法。例如,所述计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种吸附气量的确定方法,其特征在于,所述方法包括:
基于岩层中各第一样本岩心的多个地层特征的参数值、所述各第一样本岩心的吸附气量以及吸附气量函数,获取所述各第一样本岩心的第一地质特征的参数值,所述多个地层特征包括地层压力、兰氏体积和兰氏压力,所述第一地质特征的参数值用于表示地质特征对所述吸附气量的影响程度,所述吸附气量函数为所述多个地层特征的参数值、所述第一地质特征的参数值以及所述吸附气量之间的函数关系;
基于所述第一地质特征的参数值和所述各第一样本岩心的多个第二地质特征的参数值,获取所述第一地质特征的参数值与所述多个第二地质特征的参数值之间的地质特征函数,所述多个第二地质特征包括有机碳含量和孔隙度;
基于所述岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取所述地质特征函数和所述吸附气量函数的准确度;
当所述准确度高于预设阈值时,基于所述岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、所述地质特征函数以及所述吸附气量函数,获取所述岩层中每个待测岩心的吸附气量。
2.根据权利要求1所述的方法,其特征在于,所述吸附气量函数的获取过程包括:
建立初始吸附气量函数,所述初始吸附气量函数为所述多个地层特征的参数值与吸附气量之间的函数关系;
基于所述初始吸附气量函数,以及所述第一地质特征的参数值与所述吸附气量之间的函数关系,获取所述吸附气量函数。
3.根据权利要求1所述的方法,其特征在于,所述基于所述第一地质特征的参数值和所述各第一样本岩心的多个第二地质特征的参数值,获取所述第一地质特征的参数值与所述多个第二地质特征的参数值之间的地质特征函数包括:
建立初始地质特征函数;
基于所述初始地质特征函数、所述第一地质特征的参数值和所述多个第二地质特征的参数值,获取所述地质特征函数,所述地质特征函数为所述多个第二地质特征的参数值、多个系数值与所述第一特征的参数值之间的函数关系。
4.根据权利要求1所述的方法,其特征在于,所述基于所述岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取所述地质特征函数和所述吸附气量函数的准确度包括:
基于所述岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、所述地质特征函数以及所述吸附气量函数,获取所述岩层中每个第二样本岩心的理论吸附气量;
获取所述各第二样本岩心的实际吸附气量与所述各第二样本岩心的理论吸附气量之间的差值;
获取多个所述差值与所述各第二样本岩心的实际吸附气量之间的多个比值;
获取所述多个比值的绝对值的平均数,将100%与所述平均数之差获取为所述准确度。
5.根据权利要求1所述的方法,其特征在于,所述基于所述岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、所述地质特征函数以及所述吸附气量函数,获取所述岩层中每个待测岩心的吸附气量包括:
基于所述每个待测岩心的多个第二地质特征的参数值以及所述地质特征函数,获取所述每个待测岩心的第一地质特征的参数值;
基于所述每个待测岩心的多个地层特征的参数值、所述每个待测岩心的第一地质特征的参数值以及所述吸附气量函数,获取所述每个待测岩心的吸附气量。
6.一种吸附气量的确定装置,其特征在于,所述装置包括:
第一地质特征的参数值获取模块,用于基于岩层中各第一样本岩心的多个地层特征的参数值、所述各第一样本岩心的吸附气量以及吸附气量函数,获取所述各第一样本岩心的第一地质特征的参数值,所述多个地层特征包括地层压力、兰氏体积和兰氏压力,所述第一地质特征的参数值用于表示地质特征对所述吸附气量的影响程度,所述吸附气量函数为所述多个地层特征的参数值、所述第一地质特征的参数值以及所述吸附气量之间的函数关系;
地质特征函数获取模块,用于基于所述第一地质特征的参数值和所述各第一样本岩心的多个第二地质特征的参数值,获取所述第一地质特征的参数值与所述多个第二地质特征的参数值之间的地质特征函数,所述多个第二地质特征包括有机碳含量和孔隙度;
准确度获取模块,用于基于所述岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值以及实际吸附气量,获取所述地质特征函数和所述吸附气量函数的准确度;
吸附气量获取模块,用于当所述准确度高于预设阈值时,基于所述岩层中每个待测岩心的多个地层特征的参数值、多个第二地质特征的参数值、所述地质特征函数以及所述吸附气量函数,获取所述岩层中每个待测岩心的吸附气量。
7.根据权利要求6所述的装置,其特征在于,所述装置还包括吸附气量函数获取模块,所述吸附气量函数获取模块用于:
建立初始吸附气量函数,所述初始吸附气量函数为所述多个地层特征的参数值与吸附气量之间的函数关系;
基于所述初始吸附气量函数,以及所述第一地质特征的参数值与所述吸附气量之间的函数关系,获取所述吸附气量函数。
8.根据权利要求6所述的装置,其特征在于,所述地质特征函数获取模块用于:
建立初始地质特征函数;
基于所述初始地质特征函数、所述第一地质特征的参数值和所述多个第二地质特征的参数值,获取所述地质特征函数,所述地质特征函数为所述多个第二地质特征的参数值、多个系数值与所述第一特征的参数值之间的函数关系。
9.根据权利要求6所述的装置,其特征在于,所述准确度获取模块用于:
基于所述岩层中各第二样本岩心的多个地层特征的参数值、多个第二地质特征的参数值、所述地质特征函数以及所述吸附气量函数,获取所述岩层中每个第二样本岩心的理论吸附气量;
获取所述各第二样本岩心的实际吸附气量与所述各第二样本岩心的理论吸附气量之间的差值;
获取多个所述差值与所述各第二样本岩心的实际吸附气量之间的多个比值;
获取所述多个比值的绝对值的平均数,将100%与所述平均数之差获取为所述准确度。
10.根据权利要求6所述的装置,其特征在于,所述吸附气量获取模块用于:
基于所述每个待测岩心的多个第二地质特征的参数值以及所述地质特征函数,获取所述每个待测岩心的第一地质特征的参数值;
基于所述每个待测岩心的多个地层特征的参数值、所述每个待测岩心的第一地质特征的参数值以及所述吸附气量函数,获取所述每个待测岩心的吸附气量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811509638.2A CN111307652B (zh) | 2018-12-11 | 2018-12-11 | 吸附气量的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811509638.2A CN111307652B (zh) | 2018-12-11 | 2018-12-11 | 吸附气量的确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111307652A true CN111307652A (zh) | 2020-06-19 |
CN111307652B CN111307652B (zh) | 2022-07-05 |
Family
ID=71144682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811509638.2A Active CN111307652B (zh) | 2018-12-11 | 2018-12-11 | 吸附气量的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111307652B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103592687A (zh) * | 2013-11-22 | 2014-02-19 | 中国石油化工集团公司 | 煤岩吸附气含量的定量计算方法 |
CN104458489A (zh) * | 2014-12-03 | 2015-03-25 | 中国石油大学(北京) | 泥页岩吸附气含量预测方法和装置 |
CN105184106A (zh) * | 2015-10-19 | 2015-12-23 | 中国石油大学(北京) | 一种确定泥页岩储层中吸附气含量的方法及装置 |
CN105842416A (zh) * | 2015-01-16 | 2016-08-10 | 中石化西南石油工程有限公司 | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 |
CN106153495A (zh) * | 2016-08-31 | 2016-11-23 | 中国石油天然气股份有限公司 | 岩石气体吸附相体积测量方法及装置 |
CN106251232A (zh) * | 2016-08-12 | 2016-12-21 | 中国石油天然气股份有限公司 | 确定页岩含气量的方法和装置 |
CN107818379A (zh) * | 2017-09-28 | 2018-03-20 | 中国石油天然气股份有限公司 | 吸附气量确定方法和装置 |
CN108535135A (zh) * | 2018-05-24 | 2018-09-14 | 中国地质大学(北京) | 用于测量气体吸附-扩散-置换的实验系统及方法 |
US20180306009A1 (en) * | 2017-04-20 | 2018-10-25 | Korea Gas Corporation | Modeling method for gas production of cbm reservoir rocks |
-
2018
- 2018-12-11 CN CN201811509638.2A patent/CN111307652B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103592687A (zh) * | 2013-11-22 | 2014-02-19 | 中国石油化工集团公司 | 煤岩吸附气含量的定量计算方法 |
CN104458489A (zh) * | 2014-12-03 | 2015-03-25 | 中国石油大学(北京) | 泥页岩吸附气含量预测方法和装置 |
CN105842416A (zh) * | 2015-01-16 | 2016-08-10 | 中石化西南石油工程有限公司 | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 |
CN105184106A (zh) * | 2015-10-19 | 2015-12-23 | 中国石油大学(北京) | 一种确定泥页岩储层中吸附气含量的方法及装置 |
CN106251232A (zh) * | 2016-08-12 | 2016-12-21 | 中国石油天然气股份有限公司 | 确定页岩含气量的方法和装置 |
CN106153495A (zh) * | 2016-08-31 | 2016-11-23 | 中国石油天然气股份有限公司 | 岩石气体吸附相体积测量方法及装置 |
US20180306009A1 (en) * | 2017-04-20 | 2018-10-25 | Korea Gas Corporation | Modeling method for gas production of cbm reservoir rocks |
CN107818379A (zh) * | 2017-09-28 | 2018-03-20 | 中国石油天然气股份有限公司 | 吸附气量确定方法和装置 |
CN108535135A (zh) * | 2018-05-24 | 2018-09-14 | 中国地质大学(北京) | 用于测量气体吸附-扩散-置换的实验系统及方法 |
Non-Patent Citations (1)
Title |
---|
钟光海等: "四川盆地页岩气储层含气量的测井评价方法", 《天然气工业》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111307652B (zh) | 2022-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107132170B (zh) | 一种储层应力敏感性的确定方法和装置 | |
US10909281B2 (en) | History matching of hydrocarbon production from heterogenous reservoirs | |
CN111911115B (zh) | 一种页岩气井动态配产方法 | |
Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
CN105487122B (zh) | 用于确定储层流体识别因子的敏感性的方法及系统 | |
CN103485770B (zh) | 一种基于人工神经网络求取含油饱和度的方法及系统 | |
CN103527172B (zh) | 可变岩电耦合指数含水饱和度计算方法 | |
CN112796738A (zh) | 一种阵列声波测井和常规测井结合的地层渗透率计算方法 | |
CN104360412A (zh) | 致密深盆气成藏预测方法和装置 | |
EP2678718A2 (en) | Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region | |
CN109799164B (zh) | 吸附气量的确定方法及装置 | |
CN111307652B (zh) | 吸附气量的确定方法及装置 | |
CN110863818A (zh) | 一种剩余油/气分布的描述方法及装置 | |
CN110344822A (zh) | 含钙质致密储层含水饱和度的确定方法 | |
CN110857625B (zh) | 碳酸盐岩有封闭水体缝洞单元储量获取方法、装置及设备 | |
CN115271367A (zh) | 高原隧道围岩分级方法、装置、设备以及存储介质 | |
CN112419493B (zh) | 页岩储层三维属性模型建立方法及装置 | |
CN113586042B (zh) | 孔隙型碳酸盐岩动态渗透率确定方法及装置 | |
CN117248895A (zh) | 一种强非均质气藏开发早期气井合理配产图版的建立方法 | |
CN114441401A (zh) | 覆压孔隙度的确定方法、装置、计算机设备及存储介质 | |
CN114076987B (zh) | 基于致密砂岩油气藏有效储层的测井识别方法及装置 | |
CN114565116A (zh) | 一种低渗储层物性参数的预测方法和装置 | |
CN112377184A (zh) | 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置 | |
CN109828098B (zh) | 非均质性储层聚集成藏系数的确定方法、装置及存储介质 | |
CN110879424A (zh) | 流体性质识别方法及装置 |
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 |