CN113153284B - 束缚水饱和参数的确定方法、装置、设备和存储介质 - Google Patents
束缚水饱和参数的确定方法、装置、设备和存储介质 Download PDFInfo
- Publication number
- CN113153284B CN113153284B CN202110485397.8A CN202110485397A CN113153284B CN 113153284 B CN113153284 B CN 113153284B CN 202110485397 A CN202110485397 A CN 202110485397A CN 113153284 B CN113153284 B CN 113153284B
- Authority
- CN
- China
- Prior art keywords
- time point
- parameter
- determining
- target
- parameters
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 115
- 238000000034 method Methods 0.000 title claims abstract description 55
- 239000011148 porous material Substances 0.000 claims abstract description 94
- 239000011435 rock Substances 0.000 claims abstract description 75
- 230000010354 integration Effects 0.000 claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims description 80
- 230000015654 memory Effects 0.000 claims description 12
- 230000001419 dependent effect Effects 0.000 claims description 10
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 28
- 238000004590 computer program Methods 0.000 description 17
- 238000012545 processing Methods 0.000 description 11
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供了一种束缚水饱和参数的确定方法、装置、设备和存储介质，属于储层评价技术领域。所述方法包括：确定待研究地区的目标储层；对于目标储层的每个深度，获取目标储层在深度的测井图谱和总孔隙参数，测井图谱为孔隙参数与时间之间的关系图谱；基于测井图谱，确定目标储层的测井参数和辅助系数；获取目标储层的多个岩样的离心图谱，离心图谱为孔隙参数与时间之间的关系图谱；基于测井图谱和多个岩样的离心图谱，确定目标储层的积分参数和岩样参数；基于测井参数、辅助系数、积分参数和岩样参数，确定束缚水体积；基于束缚水体积和总孔隙参数，确定目标储层的束缚水饱和参数。该方法提高了确定束缚水饱和参数的效率。
Description
技术领域
本申请涉及储层评价技术领域，特别涉及一种束缚水饱和参数的确定方法、装置、设备和存储介质。
背景技术
在储层开发过程中，储层的束缚水饱和度是评价储层的核心内容，关系到储层参数计算、储量计算、油藏描述和开发方案制定。因此，确定储层的束缚水饱和度，对于高效的开发储层十分重要。
相关技术中，一般采用T2(弛豫时间)谱截止值法确定储层饱和度；该过程为：事先测得大量的代表性岩样的核磁共振T2谱，基于核磁共振T2谱确定核磁共振T2谱截止值，进而建立T2谱截止值和T2谱几何形态的关系库。在实际测井时，将实际测得的T2谱几何形态与关系库对比，确定出实际测得的T2谱几何形态对应的T2谱截止值，进而根据T2谱截止值确定束缚水饱和度。由于该方法建立关系库时需要事先测得大量的代表性岩样的核磁共振T2谱，费事费力，进而降低了确定束缚水饱和度的效率。
发明内容
本申请实施例提供了一种束缚水饱和参数的确定方法、装置、设备和存储介质，能够提高确定束缚水饱和度的效率。所述技术方案如下：
一方面，提供了一种束缚水饱和参数的确定方法，所述方法包括：
确定待研究地区的目标储层；
对于所述目标储层的每个深度，获取所述目标储层在所述深度的测井图谱和总孔隙参数，所述测井图谱为孔隙参数与时间之间的关系图谱；
基于所述测井图谱，确定所述目标储层的测井参数和辅助系数；
获取所述目标储层的多个岩样的离心图谱，所述离心图谱为孔隙参数与时间之间的关系图谱；
基于所述测井图谱和所述多个岩样的离心图谱，确定所述目标储层的积分参数和岩样参数；
基于所述测井参数、所述辅助系数、所述积分参数和所述岩样参数，确定束缚水体积；
基于所述束缚水体积和所述总孔隙参数，确定所述目标储层的束缚水饱和参数。
在一种可能的实现方式中，所述测井图谱包括在第一回波间隔时间下的第一目标图谱和在第二回波间隔时间下的第二目标图谱，所述第一回波间隔时间小于所述第二回波间隔时间；
所述基于所述测井图谱，确定所述目标储层的测井参数和辅助系数，包括：
基于所述第一目标图谱和所述第二目标图谱，确定差图谱，所述差图谱为所述第一目标图谱和所述第二目标图谱中相同时间点的第一孔隙参数和第二孔隙参数之间的差值组成的曲线；
基于所述差图谱，确定第一时间点和第二时间点，所述第一时间点和所述第二时间点分别为所述差图谱的最大峰值和最小峰值对应的时间点，将所述第一时间点和所述第二时间点作为所述测井参数；
基于所述差图谱和所述测井参数，确定所述辅助系数。
在一种可能的实现方式中，所述基于所述差图谱和所述测井参数，确定所述辅助系数，包括：
基于所述第一时间点和所述第二时间点，生成第一关系数据，所述第一关系数据的两个参数分别为所述第一时间点和所述第二时间点，且所述第一关系数据为以时间为自变量、积分值为因变量、所述辅助系数为未知量的关系数据；
基于所述差图谱，确定第三时间点，所述第三时间点为所述差图谱的最小时间点；
对于所述差图谱的每个时间点，确定所述差图谱的第三孔隙参数在所述时间点与所述第三时间点之间的积分值；
基于所述差图谱的每个时间点和每个时间点对应的积分值，对所述第一关系数据进行拟合，得到所述辅助系数。
在一种可能的实现方式中，所述测井图谱包括在第一回波间隔下的第一目标图谱；
基于所述测井图谱和多个所述离心图谱，确定所述目标储层的积分参数和岩样参数，包括：
对于每个离心图谱，确定所述离心图谱的第四时间点，所述第四时间点为所述离心图谱的最大峰值对应的时间点；
基于多个所述离心图谱的第四时间点的平均值，确定第五时间点，将所述第五时间点作为所述岩样参数；
基于所述第一目标图谱和所述第五时间点，确定所述积分参数。
在一种可能的实现方式中，所述基于所述第一目标图谱和所述第五时间点，确定所述积分参数，包括：
基于所述第一目标图谱，确定第六时间点，所述第六时间点为所述第一目标图谱的最小时间点；
确定所述第一目标图谱的第一孔隙参数在所述第六时间点与所述第五时间点之间的积分值，将所述积分值作为所述积分参数。
在一种可能的实现方式中，所述基于所述束缚水体积和所述总孔隙参数，确定所述目标储层的束缚水饱和参数，包括：
确定所述束缚水体积与所述总孔隙参数之商，得到所述目标储层的束缚水饱和参数。
在一种可能的实现方式中，所述基于所述测井参数、所述辅助系数、所述积分参数和所述岩样参数，确定束缚水体积，包括：
从所述测井参数中确定目标测井值；
基于所述目标测井值、所述辅助系数和所述岩样参数，通过第二关系数据，确定拟合参数，所述第二关系数据的因变量为所述拟合参数，自变量为所述目标测井值、所述辅助系数和所述岩样参数；
基于所述积分参数与所述拟合参数的差值，确定所述束缚水体积。
另一方面，提供了一种束缚水饱和参数的确定装置，所述装置包括：
第一确定模块，用于确定待研究地区的目标储层；
第一获取模块，用于对于所述目标储层的每个深度，获取所述目标储层在所述深度的测井图谱和总孔隙参数，所述测井图谱为孔隙参数与时间之间的关系图谱；
第二确定模块，用于基于所述测井图谱，确定所述目标储层的测井参数和辅助系数；
第二获取模块，用于获取所述目标储层的多个岩样的离心图谱，所述离心图谱为孔隙参数与时间之间的关系图谱；
第三确定模块，用于基于所述测井图谱和所述多个岩样的离心图谱，确定所述目标储层的积分参数和岩样参数；
第四确定模块，用于基于所述测井参数、所述辅助系数、所述积分参数和所述岩样参数，确定束缚水体积；
第五确定模块，用于基于所述束缚水体积和所述总孔隙参数，确定所述目标储层的束缚水饱和参数。
在一种可能的实现方式中，所述测井图谱包括在第一回波间隔时间下的第一目标图谱和在第二回波间隔时间下的第二目标图谱，所述第一回波间隔时间小于所述第二回波间隔时间；
所述第二确定模块，包括：
第一确定单元，用于基于所述第一目标图谱和所述第二目标图谱，确定差图谱，所述差图谱为所述第一目标图谱和所述第二目标图谱中相同时间点的第一孔隙参数和第二孔隙参数之间的差值组成的曲线；
第二确定单元，用于基于所述差图谱，确定第一时间点和第二时间点，所述第一时间点和所述第二时间点分别为所述差图谱的最大峰值和最小峰值对应的时间点，将所述第一时间点和所述第二时间点作为所述测井参数；
第三确定单元，用于基于所述差图谱和所述测井参数，确定所述辅助系数。
在一种可能的实现方式中，所述第三确定单元，包括：
生成子单元，用于基于所述第一时间点和所述第二时间点，生成第一关系数据，所述第一关系数据的两个参数分别为所述第一时间点和所述第二时间点，且所述第一关系数据为以时间为自变量、积分值为因变量、所述辅助系数为未知量的关系数据；
第一确定子单元，用于基于所述差图谱，确定第三时间点，所述第三时间点为所述差图谱的最小时间点；
第二确定子单元，用于对于所述差图谱的每个时间点，确定所述差图谱的第三孔隙参数在所述时间点与所述第三时间点之间的积分值；
拟合子单元，用于基于所述差图谱的每个时间点和每个时间点对应的积分值，对所述第一关系数据进行拟合，得到所述辅助系数。
在一种可能的实现方式中，所述测井图谱包括在第一回波间隔下的第一目标图谱；
所述第三确定模块，包括：
第四确定单元，用于对于每个离心图谱，确定所述离心图谱的第四时间点，所述第四时间点为所述离心图谱的最大峰值对应的时间点；
第五确定单元，用于基于多个所述离心图谱的第四时间点的平均值，确定第五时间点，将所述第五时间点作为所述岩样参数；
第六确定单元，用于基于所述第一目标图谱和所述第五时间点，确定所述积分参数。
在一种可能的实现方式中，所述第六确定单元，包括：
第三确定子单元，用于基于所述第一目标图谱，确定第六时间点，所述第六时间点为所述第一目标图谱的最小时间点；
第四确定子单元，用于确定所述第一目标图谱的第一孔隙参数在所述第六时间点与所述第五时间点之间的积分值，将所述积分值作为所述积分参数。
在一种可能的实现方式中，所述第五确定模块，包括：
第七确定单元，用于确定所述束缚水体积与所述总孔隙参数之商，得到所述目标储层的束缚水饱和参数。
在一种可能的实现方式中，所述第四确定模块，包括：
第八确定单元，用于从所述测井参数中确定目标测井值；
第九确定单元，用于基于所述目标测井值、所述辅助系数和所述岩样参数，通过第二关系数据，确定拟合参数，所述第二关系数据的因变量为所述拟合参数，自变量为所述目标测井值、所述辅助系数和所述岩样参数；
第十确定单元，用于基于所述积分参数与所述拟合参数的差值，确定所述束缚水体积。
另一方面，提供了一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器中存储至少一条程序代码，所述至少一条程序代码由所述处理器加载并执行，以实现上述任一实现方式所述的束缚水饱和参数的确定方法的指令。
另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储至少一条程序代码，所述至少一条程序代码由处理器加载并执行，以实现上述任一实现方式所述的束缚水饱和参数的确定方法中的步骤。
另一方面，提供了一种计算机程序产品或计算机程序，所述计算机程序产品或所述计算机程序包括计算机程序代码，所述计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取所述计算机程序代码，处理器执行所述计算机程序代码，使得所述计算机设备执行上述的束缚水饱和参数的确定方法所执行的操作。
本申请实施例提供的技术方案的有益效果至少包括：
本申请实施例提供了一种束缚水饱和参数的确定方法，该方法通过目标储层的每个深度的测井图谱的测井参数和辅助系数，以及离心图谱的积分参数和岩样参数，确定目标储层的束缚水饱和参数；由此可见，该方法仅需获取目标储层的每个深度的测井图谱和岩样图谱，就能够确定该深度的束缚水饱和参数，从而避免了事先测量大量的核磁共振T2谱来确定束缚水饱和参数的情况，省事省力，从而提高了确定束缚水饱和参数的效率。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种束缚水饱和参数的确定方法的流程图；
图2是本申请实施例提供的一种测井图谱；
图3是本申请实施例提供的一种差图谱；
图4是本申请实施例提供的一种束缚水饱和参数的确定结果对比图；
图5是本申请实施例提供的一种束缚水饱和参数的拟合图；
图6是本申请实施例提供的一种束缚水饱和参数的确定装置的示意图；
图7是本申请实施例提供的一种计算机设备的框图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。
本申请的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们的任意变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
本申请实施例提供了一种束缚水饱和参数的确定方法，参见图1，方法包括：
步骤101：计算机设备确定待研究地区的目标储层。
其中，待研究地区为油气田区，目标储层为待研究地区中富含油气的地层，包括砂岩储层；例如，X油田的NP凹陷砂岩储层。
步骤102：计算机设备对于目标储层的每个深度，获取目标储层在该深度的测井图谱和总孔隙参数，测井图谱为孔隙参数与时间之间的关系图谱。
其中，计算机设备对于目标储层的每个深度，基于待研究地区的在该深度的P型核磁共振测井资料，得到目标储层在该深度的测井图谱。P型核磁共振测井资料通过P型核磁共振仪器进行测井工程得到，测井图谱为通过P型核磁共振仪器测井得到的T2谱，T2谱为根据核磁共振测井资料中的回波串信号反演得到的反映孔隙流体分布的阵列图谱。
其中，P型核磁共振测井按照《核磁共振成像测井仪刻度规范SY/T 6593-2016》和《核磁共振成像测井作业技术规范SY/T 6752-2016》标准规定的流程进行。P型核磁共振测井资料处理按照《核磁共振测井资料处理及解释规范SY/T6617-2016》标准规定的流程进行标准化处理。
其中，测井图谱为孔隙参数与时间之间的关系图谱，测井图谱的纵坐标为孔隙参数，横坐标为时间；孔隙参数可以为任一表示孔隙的参数；例如，孔隙参数为孔隙分量；横坐标的时间可以为弛豫时间，弛豫时间为核磁共振测井时横向磁化强度消失的时间。总孔隙参数为该深度地层的总孔隙度，通过对测井图谱上的孔隙分量积分得到。
其中，测井图谱包括在第一回波间隔时间下的第一目标图谱和在第二回波间隔时间下的第二目标图谱，第一回波间隔时间小于第二回波间隔时间。
其中，第一目标图谱为根据核磁共振测井资料中的回波串信号反演得到的在第一回波间隔时间下的T2谱，第二目标图谱为根据核磁共振测井资料中的回波串信号反演得到的在第二回波间隔时间下的T2谱。
参见图2，图2中的实线图谱为一种第一目标图谱，虚线图谱为一种第二目标图谱。
步骤103：计算机设备基于测井图谱，确定目标储层的测井参数和辅助系数。
其中，测井参数包括测井值。该步骤可以通过以下步骤(1)-(3)实现：
(1)计算机设备基于第一目标图谱和第二目标图谱，确定差图谱，差图谱为第一目标图谱和第二目标图谱中相同时间点的第一孔隙参数和第二孔隙参数之间的差值组成的曲线。
其中，第一孔隙参数包括第一目标图谱中的每个时间点的孔隙分量，第二孔隙参数包括第二目标图谱中的每个时间点的孔隙分量。计算机设备将第一目标图谱中每个时间点的孔隙分量减去第二目标图谱中相同时间点的孔隙分量，得到每个时间点的孔隙分量差值，基于该孔隙分量差值，绘制差图谱。
其中，孔隙分量差值通过以下公式一得到：
公式一：subT2(t)＝SshortTET2(t)-SlongTET2(t)
其中，subT2(t)为时间点t对应的孔隙分量差值，SshortTET2(t)为时间点t对应的第一目标图谱中的孔隙分量，SlongrET2(t)为时间点t对应的第二目标图谱中的孔隙分量。
参见图3，图3中的虚线图谱为一种差图谱，图3中的实线为一种第一目标图谱。
(2)计算机设备基于差图谱，确定第一时间点和第二时间点，第一时间点和第二时间点分别为差图谱的最大峰值和最小峰值对应的时间点，将第一时间点和第二时间点作为测井参数。
需要说明的一点是，计算机设备可以将第一时间点作为第一测井值，第二时间点作为第二测井值；相应的，测井参数包括第一测井值和第二测井值。
继续以图3中的差图谱为例，图3中的差图谱的最大峰值对应的时间点为T2maxsub，最小峰值对应的时间点为T2minsub，将T2maxsub作为第一测井值，将T2minsub作为第二测井值。
(3)计算机设备基于差图谱和测井参数，确定辅助系数。
其中，辅助系数包括第一辅助系数、第二辅助系数和第三辅助系数。该步骤可以通过以下步骤A1-A4实现：
A1：计算机设备基于第一时间点和第二时间点，生成第一关系数据，第一关系数据的两个参数分别为第一时间点和第二时间点，且第一关系数据为以时间为自变量、积分参数为因变量、辅助系数为未知量的关系数据。
其中，第一关系数据为：
其中，int_subT2(t)为时间点t对应的积分值，φwf为第一辅助系数，k1为第二辅助系数，k2为第三辅助系数，T2maxsub为第一时间点，T2minsub为第二时间点。
A2：计算机设备基于差图谱，确定第三时间点，第三时间点为差图谱的最小时间点。
其中，最小时间点为差图谱的横坐标起始点；继续以图3的差图谱为例，最小时间点为0.3ms，即第三时间点为0.3ms。
A3：计算机设备对于差图谱的每个时间点，确定差图谱的第三孔隙参数在时间点与第三时间点之间的积分值。
其中，第三孔隙参数为差图谱上的孔隙分量。每个时间点与第三时间点之间的积分值通过对该时间点和该时间点之前的每个时间点的孔隙分量累加得到，可通过以下公式二得到：
公式二：
其中，int_subT2(t)为时间点t对应的积分值，T2_min为第三时间点，subT2(t)为时间点t对应的孔隙分量。
A4：计算机设备基于差图谱的每个时间点和每个时间点对应的积分值，对第一关系数据进行拟合，得到辅助系数。
第一关系数据为：
其中，int_subT2(t)为时间点t对应的积分值，φwf为第一辅助系数，k1为第二辅助系数，k2为第三辅助系数，T2maxsub为第一时间点，T2minsub为第二时间点。
计算机设备基于差图谱的每个时间点和每个时间点对应的积分值，对第一关系数据进行非线性最小二乘法拟合，得到第一辅助系数，第二辅助系数和第三辅助系数。
步骤104：计算机设备获取目标储层的多个岩样的离心图谱，离心图谱为孔隙参数与时间之间的关系图谱。
其中，多个岩样为与测井图谱处于同一深度的地层的岩样，多个岩样的数量至少为两个，选取岩样按照《探井化验项目取样及成果要求SY/T6028-94》标准规定的流程进行。
其中，离心图谱为通过岩样的核磁共振实验数据得到的T2谱，为第一回波间隔时间下的图谱。岩样的核磁共振实验数据为通过岩样的核磁共振实验得到的。岩样的核磁共振实验数据包括第一回波间隔时间下的离心核磁共振实验数据，离心核磁共振实验数据为将被水饱和的岩样进行离心后，进行核磁共振实验得到的数据，离心后的岩样中包括束缚水。
第一回波间隔时间下的离心核磁共振实验数据包括第一回波间隔时间下的多个孔隙参数和与多个孔隙参数对应的时间点，计算机设备基于该多个孔隙参数和与多个孔隙参数对应的时间点绘制离心图谱。
其中，岩样的核磁共振实验按照《岩样核磁共振参数实验室测量规范SY/T6490-2014》标准规定的流程进行。
步骤105：计算机设备基于测井图谱和多个岩样的离心图谱，确定目标储层的积分参数和岩样参数。
该步骤可以通过以下步骤(1)-(3)实现：
(1)计算机设备对于每个离心图谱，确定离心图谱的第四时间点，第四时间点为离心图谱的最大峰值对应的时间点。
(2)计算机设备基于多个离心图谱的第四时间点的平均值，确定第五时间点，将第五时间点作为岩样参数。
其中，计算机设备对多个第四时间点求平均值，将该平均值作为第五时间点，且将该第五时间点作为岩样参数。
(3)计算机设备基于第一目标图谱和第五时间点，确定积分参数。
该步骤可以通过以下步骤A1-A4实现：
A1：计算机设备基于第一目标图谱，确定第六时间点，第六时间点为第一目标图谱的最小时间点。
其中，最小时间点为第一目标图谱的横坐标起始点；继续参见图2或图3，图中的第一目标图谱的最小时间点为0.3ms，即第六时间点为0.3ms。
A2：计算机设备确定第一目标图谱的第一孔隙参数在第六时间点与第五时间点之间的积分值，将积分值作为积分参数。
其中，第一孔隙参数为孔隙分量，第六时间点与第五时间点之间的积分值通过对第六时间点和第五时间点之间的每个时间点的孔隙分量的加权得到，可通过以下公式三得到：
公式三：
其中，int_T2wi(T2sirr)为积分参数，T2_min为第六时间点，T2sirr为第五时间点，SshortTET2(t)为时间点t对应的孔隙分量。
步骤106：计算机设备基于测井参数、辅助系数、积分参数和岩样参数，确定束缚水体积。
其中，测井参数包括测井值。辅助系数包括第一辅助系数和第二辅助系数。
该步骤可以通过以下步骤(1)-(3)实现：
(1)计算机设备从测井参数中确定目标测井值。
其中，测井参数包括测井值，计算机设备从测井值中确定目标测井值，目标测井值为第一测井值，第一测井值为差图谱的最大峰值对应的时间点。
(2)计算机设备基于目标测井值、辅助系数和岩样参数，通过第二关系数据，确定拟合参数。
其中，第二关系数据的因变量为拟合参数，自变量为目标测井值、辅助系数和岩样参数。
第二关系数据为：
其中，W为拟合参数，φwf为第一辅助系数，k1为第二辅助系数，T2sirr为岩样参数，T2maxsub为目标测井值。
其中，计算机设备将目标测井值、第一辅助系数、第二辅助系数和岩样参数代入公式四，得到拟合参数。
(3)计算机设备基于积分参数与拟合参数的差值，确定束缚水体积。
其中，束缚水体积通过以下公式四得到：
公式四：
其中，φwi为束缚水体积，int_T2wi(T2sirr)为积分参数，W为拟合参数。
其中，计算机设备确定积分参数与拟合参数的差值，将该差值的2倍作为束缚水体积。
步骤107：计算机设备基于束缚水体积和总孔隙参数，确定目标储层的束缚水饱和参数。
其中，计算机设备确定束缚水体积与总孔隙参数之商，得到目标储层的束缚水饱和参数。
其中，束缚水饱和参数为束缚水饱和度，束缚水饱和度通过公式五得到：
公式五：
其中，Swi为束缚水饱和度，φwi为束缚水体积，
为总孔隙参数；总孔隙参数通过对第一目标图谱上的孔隙分量积分得到。T2_max为第一目标图谱上的最大时间点，T2_min为第一目标图谱上的最小时间点，SshortTET2(t)为第一目标图谱上时间点t对应的孔隙分量。
在本申请实施例中通过步骤101到步骤107的方法逐一确定目标储层的每个深度的束缚水饱和参数，进而获得整个目标储层不同深度的束缚水饱和参数。
以X油田的NP凹陷砂岩储层为例，参加图4，图4为通过本申请实施例提供的方法确定的束缚水饱和参数的结果对比图。图4中第一道为岩性曲线道，其中GR为自然伽玛曲线，单位为API，SP为自然电位曲线，单位为mv。第二道为深度道，深度单位为m。第三道为电阻率道，其中RLLD为深侧向电阻率曲线，单位为Ω.m，RLLS为浅侧向电阻率曲线，单位为Ω.m。第四道为孔隙度道，其中AC为声波曲线，单位为μs/ft，CNL为补偿中子曲线，单位为％，DEN为补偿密度曲线，单位为g/cm3。第五道和第六道为束缚水饱和度道，TASPEC为第一回波间隔时间下的T2谱，TBSPEC为第二回波间隔时间下的T2谱，T2CUT为计算的T2谱截止值曲线，单位为ms。第七道为束缚水饱和度比较道，包括岩样实验确定的束缚水饱和度、T2CUT下以33ms为固定截止值确定的束缚水饱和度、通过本申请实施例的方法确定的束缚水饱和度，通过对比可以看出，通过本申请实施例的方法确定的束缚水饱和度与岩样实验确定的束缚水饱和度吻合度更高、说明通过本申请实施例的方法确定的束缚水饱和度的准确度高于T2谱截止值法确定的束缚水饱和度。
参见图5，图5为将通过本申请实施例的方法确定的束缚水饱和度与岩样实验确定的束缚水饱和度值进行对比的拟合图；横坐标为岩心实验确定的束缚水饱和度，纵坐标为同一深度本申请实施例的方法确定的束缚水饱和度，拟合的曲线方程为y＝0.9841x，方程的相关系数R＝0.81，说明本申请实施例的方法确定的束缚水饱和度的平均误差小，与岩心实验确定的束缚水饱和度高度吻合，进一步说明了本申请实施例的方法确定的束缚水饱和度的准确度高。
本申请实施例提供了一种束缚水饱和参数的确定方法，该方法通过目标储层的每个深度的测井图谱的测井参数和辅助系数，以及离心图谱的积分参数和岩样参数，确定目标储层的束缚水饱和参数；由此可见，该方法仅需获取目标储层的每个深度的测井图谱和岩样图谱，就能够确定该深度的束缚水饱和参数，从而避免了事先测量大量的核磁共振T2谱来确定束缚水饱和参数的情况，省事省力，从而提高了确定束缚水饱和参数的效率。
本申请实施例提供了一种束缚水饱和参数的确定装置，参见图6，装置包括：
第一确定模块601，用于确定待研究地区的目标储层；
第一获取模块602，用于对于目标储层的每个深度，获取目标储层在深度的测井图谱和总孔隙参数，测井图谱为孔隙参数与时间之间的关系图谱；
第二确定模块603，用于基于测井图谱，确定目标储层的测井参数和辅助系数；
第二获取模块604，用于获取目标储层的多个岩样的离心图谱，离心图谱为孔隙参数与时间之间的关系图谱；
第三确定模块605，用于基于测井图谱和多个岩样的离心图谱，确定目标储层的积分参数和岩样参数；
第四确定模块606，用于基于测井参数、辅助系数、积分参数和岩样参数，确定束缚水体积；
第五确定模块607，用于基于束缚水体积和总孔隙参数，确定目标储层的束缚水饱和参数。
在一种可能的实现方式中，测井图谱包括在第一回波间隔时间下的第一目标图谱和在第二回波间隔时间下的第二目标图谱，第一回波间隔时间小于第二回波间隔时间；
第二确定模块603，包括：
第一确定单元，用于基于第一目标图谱和第二目标图谱，确定差图谱，差图谱为第一目标图谱和第二目标图谱中相同时间点的第一孔隙参数和第二孔隙参数之间的差值组成的曲线；
第二确定单元，用于基于差图谱，确定第一时间点和第二时间点，第一时间点和第二时间点分别为差图谱的最大峰值和最小峰值对应的时间点，将第一时间点和第二时间点作为测井参数；
第三确定单元，用于基于差图谱和测井参数，确定辅助系数。
在一种可能的实现方式中，第三确定单元，包括：
生成子单元，用于基于第一时间点和第二时间点，生成第一关系数据，第一关系数据的两个参数分别为第一时间点和第二时间点，且第一关系数据为以时间为自变量、积分值为因变量、辅助系数为未知量的关系数据；
第一确定子单元，用于基于差图谱，确定第三时间点，第三时间点为差图谱的最小时间点；
第二确定子单元，用于对于差图谱的每个时间点，确定差图谱的第三孔隙参数在时间点与第三时间点之间的积分值；
拟合子单元，用于基于差图谱的每个时间点和每个时间点对应的积分值，对第一关系数据进行拟合，得到辅助系数。
在一种可能的实现方式中，测井图谱包括在第一回波间隔下的第一目标图谱；
第三确定模块605，包括：
第四确定单元，用于对于每个离心图谱，确定离心图谱的第四时间点，第四时间点为离心图谱的最大峰值对应的时间点；
第五确定单元，用于基于多个离心图谱的第四时间点的平均值，确定第五时间点，将第五时间点作为岩样参数；
第六确定单元，用于基于第一目标图谱和第五时间点，确定积分参数。
在一种可能的实现方式中，第六确定单元，包括：
第三确定子单元，用于基于第一目标图谱，确定第六时间点，第六时间点为第一目标图谱的最小时间点；
第四确定子单元，用于确定第一目标图谱的第一孔隙参数在第六时间点与第五时间点之间的积分值，将积分值作为积分参数。
在一种可能的实现方式中，第五确定模块607，包括：
第七确定单元，用于确定束缚水体积与总孔隙参数之商，得到目标储层的束缚水饱和参数。
在一种可能的实现方式中，第四确定模块606，包括：
第八确定单元，用于从测井参数中确定目标测井值；
第九确定单元，用于基于目标测井值、辅助系数和岩样参数，通过第二关系数据，确定拟合参数，第二关系数据的因变量为拟合参数，自变量为目标测井值、辅助系数和岩样参数；
第十确定单元，用于基于积分参数与拟合参数的差值，确定束缚水体积。
图7示出了本申请一个示例性实施例提供的计算机设备700的结构框图。该计算机设备700可以是便携式移动计算机设备，比如：智能手机、平板电脑、MP3播放器(MovingPicture Experts Group Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts Group Audio Layer IV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。计算机设备700还可能被称为用户设备、便携式计算机设备、膝上型计算机设备、台式计算机设备等其他名称。
通常，计算机设备700包括有：处理器701和存储器702。
处理器701可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器701可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器701也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器701可以集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器701还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。
存储器702可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器702还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器702中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器701所执行以实现本申请中方法实施例提供的束缚水饱和参数的确定方法。
在一些实施例中，计算机设备700还可选包括有：外围设备接口703和至少一个外围设备。处理器701、存储器702和外围设备接口703之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口703相连。具体地，外围设备包括：射频电路704、显示屏705、摄像头组件706、音频电路707、定位组件708和电源709中的至少一种。
外围设备接口703可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器701和存储器702。在一些实施例中，处理器701、存储器702和外围设备接口703被集成在同一芯片或电路板上；在一些其他实施例中，处理器701、存储器702和外围设备接口703中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
射频电路704用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路704通过电磁信号与通信网络以及其他通信设备进行通信。射频电路704将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路704包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路704可以通过至少一种无线通信协议来与其它计算机设备进行通信。该无线通信协议包括但不限于：万维网、城域网、内联网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路704还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。
显示屏705用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏705是触摸显示屏时，显示屏705还具有采集在显示屏705的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器701进行处理。此时，显示屏705还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏705可以为一个，设置在计算机设备700的前面板；在另一些实施例中，显示屏705可以为至少两个，分别设置在计算机设备700的不同表面或呈折叠设计；在另一些实施例中，显示屏705可以是柔性显示屏，设置在计算机设备700的弯曲表面上或折叠面上。甚至，显示屏705还可以设置成非矩形的不规则图形，也即异形屏。显示屏705可以采用LCD(Liquid Crystal Display，液晶显示屏)、OLED(Organic Light-EmittingDiode,有机发光二极管)等材质制备。
摄像头组件706用于采集图像或视频。可选地，摄像头组件706包括前置摄像头和后置摄像头。通常，前置摄像头设置在计算机设备的前面板，后置摄像头设置在计算机设备的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件706还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。
音频电路707可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器701进行处理，或者输入至射频电路704以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在计算机设备700的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器701或射频电路704的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路707还可以包括耳机插孔。
定位组件708用于定位计算机设备700的当前地理位置，以实现导航或LBS(Location Based Service，基于位置的服务)。定位组件708可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统或俄罗斯的伽利略系统的定位组件。
电源709用于为计算机设备700中的各个组件进行供电。电源709可以是交流电、直流电、一次性电池或可充电电池。当电源709包括可充电电池时，该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池，无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。
在一些实施例中，计算机设备700还包括有一个或多个传感器710。该一个或多个传感器710包括但不限于：加速度传感器711、陀螺仪传感器712、压力传感器713、指纹传感器714、光学传感器715以及接近传感器716。
加速度传感器711可以检测以计算机设备700建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器711可以用于检测重力加速度在三个坐标轴上的分量。处理器701可以根据加速度传感器711采集的重力加速度信号，控制显示屏705以横向视图或纵向视图进行用户界面的显示。加速度传感器711还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器712可以检测计算机设备700的机体方向及转动角度，陀螺仪传感器712可以与加速度传感器711协同采集用户对计算机设备700的3D动作。处理器701根据陀螺仪传感器712采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器713可以设置在计算机设备700的侧边框和/或显示屏705的下层。当压力传感器713设置在计算机设备700的侧边框时，可以检测用户对计算机设备700的握持信号，由处理器701根据压力传感器713采集的握持信号进行左右手识别或快捷操作。当压力传感器713设置在显示屏705的下层时，由处理器701根据用户对显示屏705的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器714用于采集用户的指纹，由处理器701根据指纹传感器714采集到的指纹识别用户的身份，或者，由指纹传感器714根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器701授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器714可以被设置在计算机设备700的正面、背面或侧面。当计算机设备700上设置有物理按键或厂商Logo时，指纹传感器714可以与物理按键或厂商Logo集成在一起。
光学传感器715用于采集环境光强度。在一个实施例中，处理器701可以根据光学传感器715采集的环境光强度，控制显示屏705的显示亮度。具体地，当环境光强度较高时，调高显示屏705的显示亮度；当环境光强度较低时，调低显示屏705的显示亮度。在另一个实施例中，处理器701还可以根据光学传感器715采集的环境光强度，动态调整摄像头组件706的拍摄参数。
接近传感器716，也称距离传感器，通常设置在计算机设备700的前面板。接近传感器716用于采集用户与计算机设备700的正面之间的距离。在一个实施例中，当接近传感器716检测到用户与计算机设备700的正面之间的距离逐渐变小时，由处理器701控制显示屏705从亮屏状态切换为息屏状态；当接近传感器716检测到用户与计算机设备700的正面之间的距离逐渐变大时，由处理器701控制显示屏705从息屏状态切换为亮屏状态。
本领域技术人员可以理解，图7中示出的结构并不构成对计算机设备700的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储至少一条程序代码，所述至少一条程序代码由处理器加载并执行，以实现上述任一实现方式所述的束缚水饱和参数的确定方法中的步骤。
另一方面，提供了一种计算机程序产品或计算机程序，所述计算机程序产品或所述计算机程序包括计算机程序代码，所述计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取所述计算机程序代码，处理器执行所述计算机程序代码，使得所述计算机设备执行上述的束缚水饱和参数的确定方法所执行的操作。
在一些实施例中，本申请实施例所涉及的计算机程序可被部署在一个计算机设备上执行，或者在位于一个地点的多个计算机设备上执行，又或者，在分布在多个地点且通过通信网络互连的多个计算机设备上执行，分布在多个地点且通过通信网络互连的多个计算机设备可以组成区块链系统。
本申请实施例提供了一种束缚水饱和参数的确定方法，该方法通过目标储层的每个深度的测井图谱的测井参数和辅助系数，以及离心图谱的积分参数和岩样参数，确定目标储层的束缚水饱和参数；由此可见，该方法仅需获取目标储层的每个深度的测井图谱和岩样图谱，就能够确定该深度的束缚水饱和参数，从而避免了事先测量大量的核磁共振T2谱来确定束缚水饱和参数的情况，省事省力，从而提高了确定束缚水饱和参数的效率。
以上所述仅为本申请的可选实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (8)
1.一种束缚水饱和参数的确定方法，其特征在于，所述方法包括：
确定待研究地区的目标储层；
对于所述目标储层的每个深度，获取所述目标储层在所述深度的测井图谱和总孔隙参数，所述测井图谱为孔隙参数与时间之间的关系图谱，所述测井图谱包括在第一回波间隔时间下的第一目标图谱和在第二回波间隔时间下的第二目标图谱，所述第一回波间隔时间小于所述第二回波间隔时间；
基于所述第一目标图谱和所述第二目标图谱，确定差图谱，所述差图谱为所述第一目标图谱和所述第二目标图谱中相同时间点的第一孔隙参数和第二孔隙参数之间的差值组成的曲线；
基于所述差图谱，确定第一时间点和第二时间点，所述第一时间点和所述第二时间点分别为所述差图谱的最大峰值和最小峰值对应的时间点，将所述第一时间点和所述第二时间点作为测井参数；
基于所述第一时间点和所述第二时间点，生成第一关系数据，所述第一关系数据的两个参数分别为所述第一时间点和所述第二时间点，且所述第一关系数据为以时间为自变量、积分值为因变量、辅助系数为未知量的关系数据；
基于所述差图谱，确定第三时间点，所述第三时间点为所述差图谱的最小时间点；
对于所述差图谱的每个时间点，确定所述差图谱的第三孔隙参数在所述差图谱的每个时间点与所述第三时间点之间的积分值；
基于所述差图谱的每个时间点和每个时间点对应的积分值，对所述第一关系数据进行拟合，得到所述辅助系数；
获取所述目标储层的多个岩样的离心图谱，所述离心图谱为孔隙参数与时间之间的关系图谱；
基于所述测井图谱和所述多个岩样的离心图谱，确定所述目标储层的积分参数和岩样参数；
基于所述测井参数、所述辅助系数、所述积分参数和所述岩样参数，确定束缚水体积；
基于所述束缚水体积和所述总孔隙参数，确定所述目标储层的束缚水饱和参数。
2.根据权利要求1所述的束缚水饱和参数的确定方法，其特征在于，所述测井图谱包括在第一回波间隔下的第一目标图谱；
所述基于所述测井图谱和所述多个岩样的离心图谱，确定所述目标储层的积分参数和岩样参数，包括：
对于每个离心图谱，确定所述离心图谱的第四时间点，所述第四时间点为所述离心图谱的最大峰值对应的时间点；
基于多个所述离心图谱的第四时间点的平均值，确定第五时间点，将所述第五时间点作为所述岩样参数；
基于所述第一目标图谱和所述第五时间点，确定所述积分参数。
3.根据权利要求2所述的束缚水饱和参数的确定方法，其特征在于，所述基于所述第一目标图谱和所述第五时间点，确定所述积分参数，包括：
基于所述第一目标图谱，确定第六时间点，所述第六时间点为所述第一目标图谱的最小时间点；
确定所述第一目标图谱的第一孔隙参数在所述第六时间点与所述第五时间点之间的积分值，将所述积分值作为所述积分参数。
4.根据权利要求1所述的束缚水饱和参数的确定方法，其特征在于，所述基于所述束缚水体积和所述总孔隙参数，确定所述目标储层的束缚水饱和参数，包括：
确定所述束缚水体积与所述总孔隙参数之商，得到所述目标储层的束缚水饱和参数。
5.根据权利要求1所述的束缚水饱和参数的确定方法，其特征在于，所述基于所述测井参数、所述辅助系数、所述积分参数和所述岩样参数，确定束缚水体积，包括：
从所述测井参数中确定目标测井值；
基于所述目标测井值、所述辅助系数和所述岩样参数，通过第二关系数据，确定拟合参数，所述第二关系数据的因变量为所述拟合参数，自变量为所述目标测井值、所述辅助系数和所述岩样参数；
基于所述积分参数与所述拟合参数的差值，确定所述束缚水体积。
6.一种束缚水饱和参数的确定装置，其特征在于，所述装置包括：
第一确定模块，用于确定待研究地区的目标储层；
第一获取模块，用于对于所述目标储层的每个深度，获取所述目标储层在所述深度的测井图谱和总孔隙参数，所述测井图谱为孔隙参数与时间之间的关系图谱，所述测井图谱包括在第一回波间隔时间下的第一目标图谱和在第二回波间隔时间下的第二目标图谱，所述第一回波间隔时间小于所述第二回波间隔时间；
第二确定模块，用于基于所述第一目标图谱和所述第二目标图谱，确定差图谱，所述差图谱为所述第一目标图谱和所述第二目标图谱中相同时间点的第一孔隙参数和第二孔隙参数之间的差值组成的曲线；基于所述差图谱，确定第一时间点和第二时间点，所述第一时间点和所述第二时间点分别为所述差图谱的最大峰值和最小峰值对应的时间点，将所述第一时间点和所述第二时间点作为测井参数；基于所述第一时间点和所述第二时间点，生成第一关系数据，所述第一关系数据的两个参数分别为所述第一时间点和所述第二时间点，且所述第一关系数据为以时间为自变量、积分值为因变量、辅助系数为未知量的关系数据；基于所述差图谱，确定第三时间点，所述第三时间点为所述差图谱的最小时间点；对于所述差图谱的每个时间点，确定所述差图谱的第三孔隙参数在所述差图谱的每个时间点与所述第三时间点之间的积分值；基于所述差图谱的每个时间点和每个时间点对应的积分值，对所述第一关系数据进行拟合，得到所述辅助系数；
第二获取模块，用于获取所述目标储层的多个岩样的离心图谱，所述离心图谱为孔隙参数与时间之间的关系图谱；
第三确定模块，用于基于所述测井图谱和所述多个岩样的离心图谱，确定所述目标储层的积分参数和岩样参数；
第四确定模块，用于基于所述测井参数、所述辅助系数、所述积分参数和所述岩样参数，确定束缚水体积；
第五确定模块，用于基于所述束缚水体积和所述总孔隙参数，确定所述目标储层的束缚水饱和参数。
7.一种计算机设备，其特征在于，所述计算机设备包括一个或多个处理器和一个或多个存储器，所述一个或多个存储器中存储有至少一条指令，所述至少一条指令由所述一个或多个处理器加载并执行以实现如权利要求1至权利要求5任一项所述的束缚水饱和参数的确定方法所执行的操作。
8.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令，所述至少一条指令由处理器加载并执行以实现如权利要求1至权利要求5任一项所述的束缚水饱和参数的确定方法所执行的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110485397.8A CN113153284B (zh) | 2021-04-30 | 2021-04-30 | 束缚水饱和参数的确定方法、装置、设备和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110485397.8A CN113153284B (zh) | 2021-04-30 | 2021-04-30 | 束缚水饱和参数的确定方法、装置、设备和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113153284A CN113153284A (zh) | 2021-07-23 |
| CN113153284B true CN113153284B (zh) | 2023-06-30 |
Family
ID=76873111
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110485397.8A Active CN113153284B (zh) | 2021-04-30 | 2021-04-30 | 束缚水饱和参数的确定方法、装置、设备和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113153284B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4752882A (en) * | 1986-05-05 | 1988-06-21 | Mobil Oil Corporation | Method for determining the effective water saturation in a low-resistivity hydrocarbon-bearing rock formation based upon rock matrix conductance |
| CN101413388A (zh) * | 2008-12-02 | 2009-04-22 | 大庆油田有限责任公司 | 获得油水同层原始含油饱和度的方法及估算未试油油水同层原始含油饱和度的方法 |
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN103678861A (zh) * | 2013-05-15 | 2014-03-26 | 中国石油大学(北京) | 一种确定致密砂岩气藏成藏终止深度及成藏范围的方法 |
| CN105114064A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
| CN105223116A (zh) * | 2015-08-28 | 2016-01-06 | 中国石油天然气集团公司 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
| CN106777515A (zh) * | 2016-11-23 | 2017-05-31 | 中国海洋石油总公司 | 一种基于岩心实验资料分析致密气井产能的方法 |
| CN106991223A (zh) * | 2017-03-27 | 2017-07-28 | 中国石油天然气股份有限公司 | 确定多层油藏水驱波及系数的方法及装置 |
| CN107194180A (zh) * | 2017-05-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 油藏水驱体积波及系数的确定方法、装置及系统 |
| CN112049623A (zh) * | 2019-06-05 | 2020-12-08 | 中国石油天然气股份有限公司 | 一种中基性火山岩束缚水出水的判定方法和装置 |
| CN112147172A (zh) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | 基于核磁共振t2谱评估含水饱和度的方法和装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2443480A1 (en) * | 2009-06-19 | 2012-04-25 | Conocophillips Company | Source rock volumetric analysis |
-
2021
- 2021-04-30 CN CN202110485397.8A patent/CN113153284B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4752882A (en) * | 1986-05-05 | 1988-06-21 | Mobil Oil Corporation | Method for determining the effective water saturation in a low-resistivity hydrocarbon-bearing rock formation based upon rock matrix conductance |
| CN101413388A (zh) * | 2008-12-02 | 2009-04-22 | 大庆油田有限责任公司 | 获得油水同层原始含油饱和度的方法及估算未试油油水同层原始含油饱和度的方法 |
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN103678861A (zh) * | 2013-05-15 | 2014-03-26 | 中国石油大学(北京) | 一种确定致密砂岩气藏成藏终止深度及成藏范围的方法 |
| CN105114064A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
| CN105223116A (zh) * | 2015-08-28 | 2016-01-06 | 中国石油天然气集团公司 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
| CN106777515A (zh) * | 2016-11-23 | 2017-05-31 | 中国海洋石油总公司 | 一种基于岩心实验资料分析致密气井产能的方法 |
| CN106991223A (zh) * | 2017-03-27 | 2017-07-28 | 中国石油天然气股份有限公司 | 确定多层油藏水驱波及系数的方法及装置 |
| CN107194180A (zh) * | 2017-05-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 油藏水驱体积波及系数的确定方法、装置及系统 |
| CN112049623A (zh) * | 2019-06-05 | 2020-12-08 | 中国石油天然气股份有限公司 | 一种中基性火山岩束缚水出水的判定方法和装置 |
| CN112147172A (zh) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | 基于核磁共振t2谱评估含水饱和度的方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113153284A (zh) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111077570B (zh) | 确定油藏中辉绿岩的侵入位置的方法、装置和存储介质 | |
| CN111255434A (zh) | 气井的试井方法、装置和计算机存储介质 | |
| CN110673214B (zh) | 预测水平井的入口靶点和末端点的深度的方法和装置 | |
| CN109281664B (zh) | 水平井产油量的预测方法和装置 | |
| CN116796503A (zh) | 一种页岩有机物碳含量计算模型建立方法和装置 | |
| CN108733902B (zh) | 油藏渗透率的确定方法、装置及存储介质 | |
| CN113153284B (zh) | 束缚水饱和参数的确定方法、装置、设备和存储介质 | |
| CN112964612B (zh) | 页岩储层分形维数的确定方法和装置 | |
| CN111999140B (zh) | 制备声发射样品、确定压裂泵的方法、装置及存储介质 | |
| CN113175321B (zh) | 流体饱和参数的确定方法、装置和计算机设备 | |
| CN113586043B (zh) | 束缚水饱和参数的确定方法、装置和计算机设备 | |
| CN111812744B (zh) | 地层可压性确定方法、装置及计算机存储介质 | |
| CN111425189B (zh) | 超深层缝洞碳酸盐岩储层品质评价方法、装置及存储介质 | |
| CN112096373B (zh) | 选定页岩气区域的方法、装置及计算机存储介质 | |
| CN113139877B (zh) | 页岩气井的筛选方法及装置 | |
| CN112412444B (zh) | 注采连通强度确定方法、装置、计算机设备和存储介质 | |
| CN115263244B (zh) | 气水层的排水控制方法、装置及计算机存储介质 | |
| CN112698391B (zh) | 砂体刻画方法、装置、设备及存储介质 | |
| CN113418848B (zh) | 岩心的渗透参数的确定方法、装置、设备及存储介质 | |
| CN110058306B (zh) | 确定三维速度体的方法、装置及计算机存储介质 | |
| CN115163056A (zh) | 动态储量的确定方法、装置、计算机设备及存储介质 | |
| CN116359997A (zh) | 声波速度的确定方法、装置、设备、存储介质及产品 | |
| CN115478836A (zh) | 孔隙参数的阈值信息确定方法、装置、设备及存储介质 | |
| CN115711123A (zh) | 页岩储层的孔隙度确定方法、装置及计算机存储介质 | |
| CN110927788A (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 |