CN113253354A - 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 - Google Patents
一种核磁t2谱截止值的确定方法、装置、介质及电子设备 Download PDFInfo
- Publication number
- CN113253354A CN113253354A CN202110517574.6A CN202110517574A CN113253354A CN 113253354 A CN113253354 A CN 113253354A CN 202110517574 A CN202110517574 A CN 202110517574A CN 113253354 A CN113253354 A CN 113253354A
- Authority
- CN
- China
- Prior art keywords
- value
- nuclear magnetic
- peak
- spectrum
- spectral
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 150
- 239000011148 porous material Substances 0.000 claims abstract description 191
- 238000000034 method Methods 0.000 claims abstract description 53
- 230000003595 spectral effect Effects 0.000 claims abstract description 44
- 230000005311 nuclear magnetism Effects 0.000 claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 230000002902 bimodal effect Effects 0.000 claims description 42
- 238000005481 NMR spectroscopy Methods 0.000 claims description 21
- 238000002474 experimental method Methods 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 19
- 238000004364 calculation method Methods 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000005119 centrifugation Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000013307 optical fiber Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
-
- 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
- E21B47/00—Survey of boreholes or wells
-
- 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
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请实施例公开了一种核磁T2谱截止值的确定方法、装置、介质及电子设备。其中，方法包括：获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。本技术方案，可以实现提高核磁录井资料质量，解决核磁录井不能准确落实核磁T2谱截止值的难题，实现准确及时确定核磁T2谱截止值的目的。
Description
技术领域
本申请实施例涉及油气勘探中录井技术领域，尤其涉及一种核磁T2谱截止值的确定方法、装置、介质及电子设备。
背景技术
随着能源使用量的逐渐提高，油气勘探显得更加重要。在油气勘探过程中需要进行录井处理，核磁录井分析能够获得孔隙度、渗透率、含气饱和度、含水饱和度、可动流体饱和度、可动水饱和度以及束缚水饱和度等参数，以上数据是现场能够获得含气饱和度、可动水饱和度及束缚水饱和度等参数唯一的录井技术，是现场储层评价重要依据。而核磁T2谱截止值是核磁录井计算样品渗透率、可动流体饱和度、可动水饱和度及束缚水饱和度重要依据。所以及时快速确定核磁T2谱截止值意义重大。而现有技术中，核磁T2谱截止值分析方法最准确的是实验室离心法，该方法需要通过实验进行确定，在时效上和经济性上，无法达到油气勘探工程中的处理要求。
发明内容
本申请实施例提供一种核磁T2谱截止值的确定方法、装置、介质及电子设备，可以实现提高核磁录井资料质量，解决核磁录井不能准确落实核磁T2谱截止值的难题，实现准确及时确定核磁T2谱截止值的目的。
第一方面，本申请实施例提供了一种核磁T2谱截止值的确定方法，所述方法包括：
获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
进一步的，根据所述核磁录井数据确定核磁T2谱的形状类型，包括：
获取核磁录井数据中的T2谱数据点和孔隙信号；
根据所述孔隙信号的波动情况，确定核磁T2谱的波峰个数以及凹点个数；
根据所述波峰个数以及凹点个数，确定核磁T2谱的形状类型为单峰类型或者为双峰类型；其中，所述单峰类型包括单峰类型和包括至少三个波峰的多峰类型。
进一步的，所述核磁T2谱孔隙结构参数，包括：
最小孔喉孔隙半径值X起；频数最高的孔喉孔隙半径值X峰；最高峰对应的T2谱数据点参数值Y峰；频数最高的孔隙半径值的对数值Ln(X峰)；单峰类型起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数值积分面积S前峰；双峰类型第一个峰对应的T2谱数据点参数值X峰1；双峰类型第二个峰对应的T2谱数据点参数值X峰2；双峰类型凹点到第二个峰之间的T2谱数据点参数值积分面积S后峰；双峰类型第二个峰对应的T2谱数据点参数值Y峰2。
进一步的，根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数，包括：
利用Excel表格中的CORREL公式计算各核磁T2谱孔隙结构参数与T2谱截止值的相关系数；
根据所述相关系数的大小关系，确定目标孔隙结构参数。
进一步的，根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式，包括：
构建初始多参数T2谱截止值拟合公式；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，确定所述初始多参数T2谱截止值拟合公式中各目标孔隙结构参数的权重系数；
根据所述各目标孔隙结构参数的权重系数，确定当前形状类型对应T2谱截止值拟合公式。
进一步的，采用Excel表格中LINEST、INDEX函数组合返回初始多参数T2谱截止值拟合公式中各目标孔隙结构参数的权重系数。
进一步的，单峰类型的各目标孔隙结构参数包括：
单峰类型起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数值积分面积S前峰；最高峰对应的T2谱数据点参数值Y峰；最小孔喉孔隙半径值X起；频数最高的孔隙半径值的对数值Ln(X峰)；频数最高的孔喉孔隙半径值X峰；
双峰类型的各目标孔隙结构参数包括：
双峰类型第二个峰对应的T2谱数据点参数值Y峰2；双峰类型凹点到第二个峰之间的T2谱数据点参数值积分面积S后峰；双峰类型第二个峰对应的T2谱数据点参数值X峰2；最小孔喉孔隙半径值X起；双峰类型第一个峰对应的T2谱数据点参数值X峰1。
第二方面，本申请实施例提供了一种核磁T2谱截止值的确定装置，该装置包括：
形状类型确定模块，用于获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
孔隙结构参数确定模块，用于根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
目标孔隙结构参数筛选模块，用于根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
拟合公式确定模块，用于根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
第三方面，本申请实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本申请实施例所述的核磁T2谱截止值的确定方法。
第四方面，本申请实施例提供了一种电子设备，包括存储器，处理器及存储在存储器上并可在处理器运行的计算机程序，所述处理器执行所述计算机程序时实现如本申请实施例所述的核磁T2谱截止值的确定方法。
本申请实施例提供的技术方案，通过获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。通过执行本方案，可以实现提高核磁录井资料质量，解决核磁录井不能准确落实核磁T2谱截止值的难题，实现准确及时确定核磁T2谱截止值的目的。
附图说明
图1是本申请实施例一提供的核磁T2谱截止值的确定方法的流程图；
图2是本申请实施例一提供的正态单峰核磁T2谱与截止值的对应示意图；
图3是本申请实施例一提供的前单峰核磁T2谱与截止值的对应示意图；
图4是本申请实施例一提供的后单峰核磁T2谱与截止值的对应示意图；
图5是本申请实施例一提供的对称双峰核磁T2谱与截止值的对应示意图；
图6是本申请实施例一提供的前主双峰核磁T2谱与截止值的对应示意图；
图7是本申请实施例一提供的后主双峰核磁T2谱与截止值的对应示意图；
图8为本发明实施例二提供的核磁T2谱截止值的确定装置的结构框图；
图9是本申请实施例四提供的电子设备的结构示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本申请，而非对本申请的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本申请相关的部分而非全部结构。
在更加详细地讨论示例性实施例之前应当提到的是，一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各步骤描述成顺序的处理，但是其中的许多步骤可以被并行地、并发地或者同时实施。此外，各步骤的顺序可以被重新安排。当其操作完成时所述处理可以被终止，但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。
实施例一
图1是本申请实施例一提供的核磁T2谱截止值的确定方法的流程图，本实施例可适用于核磁录井场景，该方法可以由本申请实施例所提供的核磁T2谱截止值的确定装置执行，该装置可以由软件和/或硬件的方式来实现，并可集成于用于核磁录井的电子设备中。
本实施例为了解决核磁录井T2谱截止值动态确定问题，在收集以往离心法确定截止值参数的基础上，同时开展核磁T2谱截止值的离心实验，随机选取各地区及各层位砂岩样品202块，确定各种形状核磁T2谱截止值与孔隙结构参数的相关性。
通过对所分析的核磁T2谱与截止值对比发现，截止值与核磁T2谱的谱峰的起落位置及峰高都具有一定的关联，谱峰的起落位置代表储层的孔隙半径，峰高代表孔隙的大小，因此，截止值与核磁的孔隙结构具有较好的对应关系，并且核磁谱图与截止值具有一个相对固定的对应关系。
应用这些核磁T2谱表征储层孔隙结构参数，与实验室离心法确定的T2截止值进行相关性分析，优选参数，建立优选参数与实验室离心法确定的T2截止值的对应关系，用数学方法建立拟合公式，实现核磁录井T2谱截止值的动态确定。
如图1所示，所述核磁T2谱截止值的确定方法包括：
S110，获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型。
其中，核磁录井的核磁T2谱可以是由分析软件分析后，导出到Excel中。导出的文档格式最全时为三部分，最左部分为T2谱数据点。包括T2谱数据点、干样信号、孔隙信号及油信号。T2谱数据点的单位为ms。孔隙信号无单位。当一次分析或二次分析时没有含油信号或干样信号。中间为核磁T2驰豫谱图。谱图纵坐标为信号强度，横坐标为弛豫时间。右边为回波串数据点。
本方法涉及文档部分是最左部分的T2谱数据点，这部分涉及的只有T2谱数据点及孔隙信号两列数据。其它部分及其它数据均不涉及。
核磁T2谱图形状是根据孔隙信号进行划分的，分为以下二种：第一种就是双峰，第二种是单峰；其中，如果核磁T2谱为多峰时，实验确认可以按单峰处理。
在本实施例中，可选的，根据所述核磁录井数据确定核磁T2谱的形状类型，包括：
获取核磁录井数据中的T2谱数据点和孔隙信号；
根据所述孔隙信号的波动情况，确定核磁T2谱的波峰个数以及凹点个数；
根据所述波峰个数以及凹点个数，确定核磁T2谱的形状类型为单峰类型或者为双峰类型；其中，所述单峰类型包括单峰类型和包括至少三个波峰的多峰类型。
具体的，核磁T2谱图形状确认方法如下：对孔隙信号逐个相减，就是后一个减前一个，如果值是正的，说明孔隙信号在增加，如果某一孔隙信号后的差值是负的，那这一孔隙信号就是本井深核磁数据的第一个峰，如果在差值在负的情况下，到某一点之后的差值为正的，则这一点的孔隙信号是本井深核磁数据和第一点凹点，如果核磁谱图只有一个峰，那么这个谱图也不会出现凹点，如果核磁谱图有两个峰，则会有一个凹点，如果有n个峰，则有n-1个凹点。如果核磁T2谱为多峰时，按单峰处理，最高峰相当于单峰的最高峰。
S120，根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数。
其中，在确定形状类型之后，可以针对不同的形状类型，确定不同的核磁T2谱孔隙结构参数。
本方案中，可选的，所述核磁T2谱孔隙结构参数，包括：
最小孔喉孔隙半径值X起；频数最高的孔喉孔隙半径值X峰；最高峰对应的T2谱数据点参数值Y峰；频数最高的孔隙半径值的对数值Ln(X峰)；单峰类型起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数值积分面积S前峰；双峰类型第一个峰对应的T2谱数据点参数值X峰1；双峰类型第二个峰对应的T2谱数据点参数值X峰2；双峰类型凹点到第二个峰之间的T2谱数据点参数值积分面积S后峰；双峰类型第二个峰对应的T2谱数据点参数值Y峰2。
具体的，可以通过压汞实验确认，核磁谱图的参数与储层的孔隙结构具有较好的相关性，现介绍一下核磁谱图表征孔隙结构的主要参数。
X起：T2谱数据点参数，指最小孔隙半径值。是孔隙信号起始部分最小值对应的T2谱数据点值。例如孔隙信号起始部分最小值为0.09，它对应的T2谱数据点为0.1ms。
X峰：T2谱数据点参数，是指频数最高的孔隙半径值。是数据表中单峰或多峰的孔隙信号最高点处对应的T2谱数据点值。例如孔隙信号最高点值为85.47，它对应的T2谱数据点为5.57ms。
Y峰：孔隙信号参数，单峰或多峰核磁谱图中最高峰对应的T2谱数据点参数(孔隙半径)值。例如孔隙信号最高峰为85.47。
Ln(X峰)：频数最高的孔隙半径值的对数值。例如是ln(5.57)。
S前峰：孔隙信号参数，单峰或多峰谱图中起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数(孔隙半径)值积分面积。例如孔隙信号从0.09到85.47之间的孔隙信号积分值，为849.37。
X峰1：T2谱数据点参数，双峰核磁谱图中第一个峰对应的T2谱数据点参数(孔隙半径)值。
X峰2：T2谱数据点参数，双峰核磁谱图中第二个峰对应的T2谱数据点参数(孔隙半径)值。
S后峰：孔隙信号参数，双峰核磁谱图中双峰凹点到最后一个孔隙信号值之间的T2谱数据点参数(孔隙半径)值积分面积。
Y峰2：孔隙信号参数，双峰核磁谱图中第二个峰对应的T2谱数据点参数(孔隙半径)值。
其中，通过与实验，总结核磁谱图与离心法截止值对应分析，核磁T2谱与截止值具有以下特征。
一，单峰类型的核磁T2谱截止值：
单峰共分为三种，为正态单峰、前单峰及后单峰。正态单峰的截止值在峰顶左右，前单峰在右侧峰上，后单峰截止值在左侧峰上。
图2是本申请实施例一提供的正态单峰核磁T2谱与截止值的对应示意图，图3是本申请实施例一提供的前单峰核磁T2谱与截止值的对应示意图，图4是本申请实施例一提供的后单峰核磁T2谱与截止值的对应示意图，如图2、图3和图4所示，正态单峰的截止值在峰顶左右，前单峰在右侧峰上，后单峰截止值在左侧峰上。
二，双峰型的核磁T2谱截止值：
双峰共分为三种，为对称双峰、前主双峰及后主双峰。对称双峰的截止值在双峰凹点左右，前主双峰在前峰右侧峰上，后主双峰截止值在后峰的左侧峰上。
图5是本申请实施例一提供的对称双峰核磁T2谱与截止值的对应示意图，图6是本申请实施例一提供的前主双峰核磁T2谱与截止值的对应示意图，图7是本申请实施例一提供的后主双峰核磁T2谱与截止值的对应示意图，如图5、图6和图7所示，对称双峰的截止值在双峰凹点左右，前主双峰在前峰右侧峰上，后主双峰截止值在后峰的左侧峰上。
基于此，可以得到核磁T2谱与截止值分析数据表，如下表1所示：
表1
S130，根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数。
具体的，可以把核磁T2谱孔隙结构参数与实验获取的T2截止值进行对比，确定其相关性。目标孔隙结构参数可以是经过筛选得到的优选的孔隙结构参数。本方案中，可以根据相关性的大小进行筛选。
本方案中，可选的，根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数，包括：
利用Excel表格中的CORREL公式计算各核磁T2谱孔隙结构参数与T2谱截止值的相关系数；
根据所述相关系数的大小关系，确定目标孔隙结构参数。
此处选用Excel表格中的CORREL公式计算相关性。
相关系数计算公式为：ρxy＝CORREL($A$2:$A$65000,F$2:F$65000)
其中，A-E是X起等一系列数据，F一列是T2谱截止值，ρxy即为相关性的数值。
本方案使用的CORREL函数是一个统计函数，这个函数是用来统计单元格区域相关性的函数。这样设置的好处是计算过程更加准确和便捷。
单峰类型样品的孔隙结构参数与截止值汇总结果，可以参照下表2：
表2
| 截止值 | X<sub>起</sub> | X<sub>峰</sub> | Y<sub>峰</sub> | S<sub>前峰</sub> | lnX<sub>峰</sub> |
| 5.83 | 0.1000 | 1.55 | 796.2 | 4989.19 | 2.740840024 |
| 53.01 | 0.1000 | 215.44 | 501.38 | 4934.01 | 7.67 267541 |
| 31.79 | 0.1000 | 179 46 | 1998.8 | 18023.09 | 7.492537435 |
| 15.06 | 0.1000 | 71.97 | 2181.2 | 15852.15 | 6.57883445 |
| 5.29 | 0.1000 | 1.55 | 1167.33 | 7306.52 | 2.740840024 |
| 18.19 | 0.1000 | 86.4 | 1148.21 | 14792.92 | 6.761572769 |
| 7.51 | 0.1000 | 1.29 | 607.8 | 3590.98 | 2.557227311 |
| 9.53 | 0.1000 | 1.55 | 328.96 | 2019.82 | 2.740840024 |
| 28.47 | 0.1000 | 124.52 | 1408.82 | 17127.77 | 7.127051439 |
| 14.86 | 0.1000 | 59.95 | 1236.9 | 15232.95 | 6.396095974 |
| 6.18 | 0.1000 | 2.23 | 417.87 | 3337.92 | 3.104586678 |
| 16.26 | 0.1000 | 59.95 | 1492.9 | 18402.15 | 6.396095974 |
| … | … | … | … | … | … |
通过把上述的单峰类型孔隙结构参数与截止值带入相关系数计算公式中，就可以得到相应的相关性表，如下表3：
表3
| X<sub>峰</sub> | lnX<sub>峰</sub> | X<sub>起</sub> | Y<sub>峰</sub> | S<sub>前峰</sub> |
| 0.9599 | 0.73 | 0.6376 | 0.5535 | 0.5455 |
针对双峰类型样品的孔隙结构参数与截止值汇总结果，可参见表4：
表4
把双峰类型孔隙结构参数与截止值带入相关系数计算公式中，可以得到相应的相关性表，见表5：
表5
| X<sub>峰1</sub> | X<sub>起</sub> | X<sub>峰2</sub> | S<sub>后峰</sub> | Y<sub>峰2</sub> |
| 0.9253 | 0.7291 | 0.7079 | 0.6203 | 0.6146 |
在上述各技术方案中，可选的，单峰类型的各目标孔隙结构参数包括：
单峰类型起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数值积分面积S前峰；最高峰对应的T2谱数据点参数值Y峰；最小孔喉孔隙半径值X起；频数最高的孔隙半径值的对数值Ln(X峰)；频数最高的孔喉孔隙半径值X峰；
双峰类型的各目标孔隙结构参数包括：
双峰类型第二个峰对应的T2谱数据点参数值Y峰2；双峰类型凹点到第二个峰之间的T2谱数据点参数值积分面积S后峰；双峰类型第二个峰对应的T2谱数据点参数值X峰2；最小孔喉孔隙半径值X起；双峰类型第一个峰对应的T2谱数据点参数值X峰1。
其中，相应的目标孔隙结构参数的确定可以是基于相关性的大小来确定的，这样设置的好处是得到的你和公式能够更加准确的反映出对核磁T2谱截止值影响较大的因素，使得核磁T2谱截止值的拟合结果更加准确。
S140，根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
本方案中，可以预先确定一个拟合公式，并通过将上述的相关性数据代入的方式，确定拟合公式中的具体的参数数值。在得到参数数值之后，可以根据拟合公式，对后续的目标测试对象进行T2谱截止值的计算。从而，可以实现快速并准确的得到T2谱截止值的计算结果的目的，为后续的测试对象的岩性判断提供数据基础。
在本实施例中，可选的，根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式，包括：
构建初始多参数T2谱截止值拟合公式；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，确定所述初始多参数T2谱截止值拟合公式中各目标孔隙结构参数的权重系数；
根据所述各目标孔隙结构参数的权重系数，确定当前形状类型对应T2谱截止值拟合公式。
具体的，初始多参数T2谱截止值拟合公式，可以是采用多参数截止值拟合公式为：y＝m1×x1+m2×x2+m3×x3+…+mn×xn+b，其中的x1，x2，x3等可以是影响因子，即可以是前面确定的目标目标孔隙结构参数，m1，m2，m3等可以是影响因子的权重系数，即为各目标孔隙结构参数的权重系数，其中的b可以是一个常数。
本方案中，可选的，采用Excel表格中LINEST、INDEX函数组合返回初始多参数T2谱截止值拟合公式中各目标孔隙结构参数的权重系数。
具体的，采用Excel表格中LINEST、INDEX函数组合返回多重回归公式的系数。多重回归公式的系数计算公式为：
mn＝INDEX(LINEST($F$2:$F$103,$A$2:$E$103,1,FALSE),1,n)；
b＝INDEX(LINEST($F$2:$F$103,$A$2:$E$103,1,FALSE),1,n+1)。
通过这样的计算，可以准确并快速的得到当前形状类型对应T2谱截止值拟合公式。
具体的，单峰类型样品的孔隙结构参数与截止值拟合结果，可以见下表6：
表6
通过拟合，可以获得各目标孔隙结构参数的权重系数，如下表7：
表7
| 参数 | m1 | m2 | m3 | m4 | m5 | b |
| 拟合值 | -0.00068 | 0.005116 | 11.61061 | 1.556617 | 0.149799 | 0.962681 |
由此，可以得到单峰类型对应T2谱截止值拟合公式，如下所示：
单峰类型T2谱截止值＝
－0.000679785029174963×S前峰+0.00511586740746969×
Y峰+11.6106125967902×X起+1.55661701157281×
lnX峰+0.149798707164228×X峰+0.962680568749449；
与此对应的，也可以对双峰型样品的孔隙结构参数与截止值进行拟合，得到拟合结果，如下表8：
表8
通过拟合，可以获得各目标孔隙结构参数的权重系数，如下表9：
表9
| 参数 | m1 | m2 | m3 | m4 | m5 | b |
| 拟合值 | -0.00806 | 0.00075 | 0.02120 | -8.26244 | 2.61709 | 2.19610 |
由此，可以得到双峰类型对应T2谱截止值拟合公式，如下所示：
单峰类型T2谱截止值＝
－0.0080562646341124×Y峰2+0.000747591975362114×
S后峰+0.021202464778358×X峰2－8.26243834612535×
X起+2.61709028805596×X峰1+2.19610294055626；
本方案通过利用核磁T2谱截止值与核磁孔隙结构参数相对应，其值对相应谱图来说是固定的。通过大数据实现了优选核磁谱图的孔隙结构参数，并与实验室核磁T2截止值建立对应关系并拟合，实现及时准确计算截止值，打破行业目前使用区域固定经验值(14ms)的及采用谱图法大致确定截止值的做法，对核磁录井具有重要意义；而且，本方案能够获得准确的核磁T2截止值，确保了核磁录井能够准确判别储层物性及流体性质，利于油田高效勘探开发；另外，该方法适用于核磁录井样品分析，应用前景广阔。
本申请实施例所提供的技术方案，通过获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。通过执行本方案，可以实现提高核磁录井资料质量，解决核磁录井不能准确落实核磁T2谱截止值的难题，实现准确及时确定核磁T2谱截止值的目的。
实施例二
图8为本发明实施例二提供的核磁T2谱截止值的确定装置的结构框图，该装置可执行本发明任意实施例所提供的核磁T2谱截止值的确定方法，具备执行方法相应的功能模块和有益效果。如图8所示，该装置可以包括：
形状类型确定模块810，用于获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
孔隙结构参数确定模块820，用于根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
目标孔隙结构参数筛选模块830，用于根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
拟合公式确定模块840，用于根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
上述产品可执行本申请实施例所提供的核磁T2谱截止值的确定方法，具备执行方法相应的功能模块和有益效果。
实施例三
本发明实施例三提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本申请所有发明实施例提供的核磁T2谱截止值的确定方法：
获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括——但不限于——电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括——但不限于——无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
实施例四
本申请实施例四提供了一种电子设备，该电子设备中可集成本申请实施例提供的核磁T2谱截止值的确定方法。图9是本申请实施例四提供的电子设备的结构示意图。如图9所示，本实施例提供了一种电子设备900，其包括：一个或多个处理器920；存储装置910，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器920执行，使得所述一个或多个处理器920实现本申请实施例所提供的核磁T2谱截止值的确定方法，该方法包括：
获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
当然，本领域技术人员可以理解，处理器920还实现本申请任意实施例所提供的核磁T2谱截止值的确定方法的技术方案。
图9显示的电子设备900仅仅是一个示例，不应对本申请实施例的功能和使用范围带来任何限制。
如图9所示，该电子设备900包括处理器920、存储装置910、输入装置930和输出装置940；电子设备中处理器920的数量可以是一个或多个，图9中以一个处理器920为例；电子设备中的处理器920、存储装置910、输入装置930和输出装置940可以通过总线或其他方式连接，图9中以通过总线950连接为例。
存储装置910作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块单元，如本申请实施例中的核磁T2谱截止值的确定方法对应的程序指令。
存储装置910可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储装置910可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储装置910可进一步包括相对于处理器920远程设置的存储器，这些远程存储器可以通过网络连接。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置930可用于接收输入的数字、字符信息或语音信息，以及产生与电子设备的用户设置以及功能控制有关的键信号输入。输出装置940可包括显示屏、扬声器等电子设备。
本申请实施例提供的电子设备，可以实现提高核磁录井资料质量，解决核磁录井不能准确落实核磁T2谱截止值的难题，实现准确及时确定核磁T2谱截止值的目的。
上述实施例中提供的核磁T2谱截止值的确定装置、介质及电子设备可执行本申请任意实施例所提供的核磁T2谱截止值的确定方法，具备执行该方法相应的功能模块和有益效果。未在上述实施例中详尽描述的技术细节，可参见本申请任意实施例所提供的核磁T2谱截止值的确定方法。
注意，上述仅为本申请的较佳实施例及所运用技术原理。本领域技术人员会理解，本申请不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本申请的保护范围。因此，虽然通过以上实施例对本申请进行了较为详细的说明，但是本申请不仅仅限于以上实施例，在不脱离本申请构思的情况下，还可以包括更多其他等效实施例，而本申请的范围由所附的权利要求范围决定。
Claims (10)
1.一种核磁T2谱截止值的确定方法，其特征在于，所述方法包括：
获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
2.根据权利要求1所述的方法，其特征在于，根据所述核磁录井数据确定核磁T2谱的形状类型，包括：
获取核磁录井数据中的T2谱数据点和孔隙信号；
根据所述孔隙信号的波动情况，确定核磁T2谱的波峰个数以及凹点个数；
根据所述波峰个数以及凹点个数，确定核磁T2谱的形状类型为单峰类型或者为双峰类型；其中，所述单峰类型包括单峰类型和包括至少三个波峰的多峰类型。
3.根据权利要求2所述的方法，其特征在于，所述核磁T2谱孔隙结构参数，包括：
最小孔喉孔隙半径值X起；频数最高的孔喉孔隙半径值X峰；最高峰对应的T2谱数据点参数值Y峰；频数最高的孔隙半径值的对数值Ln(X峰)；单峰类型起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数值积分面积S前峰；双峰类型第一个峰对应的T2谱数据点参数值X峰1；双峰类型第二个峰对应的T2谱数据点参数值X峰2；双峰类型凹点到第二个峰之间的T2谱数据点参数值积分面积S后峰；双峰类型第二个峰对应的T2谱数据点参数值Y峰2。
4.根据权利要求1所述的方法，其特征在于，根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数，包括：
利用Excel表格中的CORREL公式计算各核磁T2谱孔隙结构参数与T2谱截止值的相关系数；
根据所述相关系数的大小关系，确定目标孔隙结构参数。
5.根据权利要求1所述的方法，其特征在于，根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式，包括：
构建初始多参数T2谱截止值拟合公式；
根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，确定所述初始多参数T2谱截止值拟合公式中各目标孔隙结构参数的权重系数；
根据所述各目标孔隙结构参数的权重系数，确定当前形状类型对应T2谱截止值拟合公式。
6.根据权利要求5所述的方法，其特征在于，采用Excel表格中LINEST、INDEX函数组合返回初始多参数T2谱截止值拟合公式中各目标孔隙结构参数的权重系数。
7.根据权利要求5所述的方法，其特征在于，单峰类型的各目标孔隙结构参数包括：
单峰类型起始孔隙信号到最高孔隙信号值之间的T2谱数据点参数值积分面积S前峰；最高峰对应的T2谱数据点参数值Y峰；最小孔喉孔隙半径值X起；频数最高的孔隙半径值的对数值Ln(X峰)；频数最高的孔喉孔隙半径值X峰；
双峰类型的各目标孔隙结构参数包括：
双峰类型第二个峰对应的T2谱数据点参数值Y峰2；双峰类型凹点到第二个峰之间的T2谱数据点参数值积分面积S后峰；双峰类型第二个峰对应的T2谱数据点参数值X峰2；最小孔喉孔隙半径值X起；双峰类型第一个峰对应的T2谱数据点参数值X峰1。
8.一种核磁T2谱截止值的确定装置，其特征在于，所述装置包括：
形状类型确定模块，用于获取核磁录井数据，根据所述核磁录井数据确定核磁T2谱的形状类型；
孔隙结构参数确定模块，用于根据所述形状类型，确定与当前形状类型对应的至少一个核磁T2谱孔隙结构参数；
目标孔隙结构参数筛选模块，用于根据各核磁T2谱孔隙结构参数与T2谱截止值的相关系数，确定目标孔隙结构参数；
拟合公式确定模块，用于根据各目标孔隙结构参数与实验获取的标准T2谱截止值进行拟合，得到当前形状类型对应T2谱截止值拟合公式；并根据所述T2谱截止值拟合公式，确定目标测试对象的T2谱截止值。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-7中任一项所述的核磁T2谱截止值的确定方法。
10.一种电子设备，包括存储器，处理器及存储在存储器上并可在处理器运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1-7中任一项所述的核磁T2谱截止值的确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110517574.6A CN113253354A (zh) | 2021-05-12 | 2021-05-12 | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110517574.6A CN113253354A (zh) | 2021-05-12 | 2021-05-12 | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113253354A true CN113253354A (zh) | 2021-08-13 |
Family
ID=77223160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110517574.6A Pending CN113253354A (zh) | 2021-05-12 | 2021-05-12 | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113253354A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116859473A (zh) * | 2023-05-23 | 2023-10-10 | 贵州省油气勘查开发工程研究院 | 面向非常规储层含水饱和度的形态学和分形学预测方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO984262L (no) * | 1997-09-16 | 1999-03-17 | Halliburton Energy Serv Inc | FremgangsmÕte og anordning for evaluering av geologiske formasjoner |
| CN1869733A (zh) * | 2005-05-27 | 2006-11-29 | 中国石油天然气股份有限公司 | 核磁共振测井t2谱t2截止值的确定方法 |
| US20130103319A1 (en) * | 2011-10-21 | 2013-04-25 | Saudi Arabian Oil Company | Methods for determining well characteristics and pore architecture utilizing conventional well logs |
| US20140035574A1 (en) * | 2012-08-03 | 2014-02-06 | Chevron U.S.A. Inc. | System and method for estimating a nuclear magnetic resonance relaxation time cutoff |
| CN104777181A (zh) * | 2014-12-31 | 2015-07-15 | 中国石油天然气股份有限公司 | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 |
| CN109580689A (zh) * | 2018-10-16 | 2019-04-05 | 中国石油天然气集团有限公司 | 一种核磁共振测井t2截止值的逐点计算方法 |
| CN110231272A (zh) * | 2019-07-09 | 2019-09-13 | 中国地质大学（北京） | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 |
| WO2019232299A1 (en) * | 2018-06-01 | 2019-12-05 | Board Of Regents, The University Of Texas System | Method for assessment of pore-throat size distribution and permeability in porous media |
| CN111289552A (zh) * | 2020-01-19 | 2020-06-16 | 中国地质大学（北京） | 煤的核磁共振截止值确定方法以及装置 |
| CN112067645A (zh) * | 2020-07-31 | 2020-12-11 | 中国石油大学(北京) | 核磁共振弛豫时间与孔喉半径的转换方法 |
-
2021
- 2021-05-12 CN CN202110517574.6A patent/CN113253354A/zh active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO984262L (no) * | 1997-09-16 | 1999-03-17 | Halliburton Energy Serv Inc | FremgangsmÕte og anordning for evaluering av geologiske formasjoner |
| CN1869733A (zh) * | 2005-05-27 | 2006-11-29 | 中国石油天然气股份有限公司 | 核磁共振测井t2谱t2截止值的确定方法 |
| US20130103319A1 (en) * | 2011-10-21 | 2013-04-25 | Saudi Arabian Oil Company | Methods for determining well characteristics and pore architecture utilizing conventional well logs |
| US20140035574A1 (en) * | 2012-08-03 | 2014-02-06 | Chevron U.S.A. Inc. | System and method for estimating a nuclear magnetic resonance relaxation time cutoff |
| CN104777181A (zh) * | 2014-12-31 | 2015-07-15 | 中国石油天然气股份有限公司 | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 |
| WO2019232299A1 (en) * | 2018-06-01 | 2019-12-05 | Board Of Regents, The University Of Texas System | Method for assessment of pore-throat size distribution and permeability in porous media |
| CN109580689A (zh) * | 2018-10-16 | 2019-04-05 | 中国石油天然气集团有限公司 | 一种核磁共振测井t2截止值的逐点计算方法 |
| CN110231272A (zh) * | 2019-07-09 | 2019-09-13 | 中国地质大学（北京） | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 |
| CN111289552A (zh) * | 2020-01-19 | 2020-06-16 | 中国地质大学（北京） | 煤的核磁共振截止值确定方法以及装置 |
| CN112067645A (zh) * | 2020-07-31 | 2020-12-11 | 中国石油大学(北京) | 核磁共振弛豫时间与孔喉半径的转换方法 |
Non-Patent Citations (7)
| Title |
|---|
| YANG ZHAO等: "Multifractal analysis of coal pore structure based on NMR experiment: A new method for predicting T2 cutoff value", 《FUEL》, vol. 283, pages 1 - 12 * |
| 宋明会: "核磁共振录井可动流体t_2截止值确定方法", 《录井工程》, vol. 18, no. 03, pages 5 - 8 * |
| 李建玉;赵辉;张凤生: "岩心核磁共振实验分析与参数评价", 《内蒙古石油化工》, vol. 36, no. 17, pages 1 - 4 * |
| 范宜仁;刘建宇;葛新民;邓少贵;刘洪亮;顾定娜: "基于核磁共振双截止值的致密砂岩渗透率评价新方法", 《地球物理学报》, vol. 61, no. 04, pages 1628 - 1638 * |
| 邵维志;丁娱娇;肖斐;慕红武;刘洪波;蔡浩;董旭锋: "利用T_2谱形态确定T_2截止值的方法探索", 《测井技术》, vol. 33, no. 05, pages 430 - 435 * |
| 郑可;徐怀民;陈建文;王永宏: "低渗储层可动流体核磁共振研究", 《现代地质》, vol. 27, no. 03, pages 710 - 718 * |
| 黄兴;李天太;王香增;高辉;倪军;赵金省;王琛: "致密砂岩储层可动流体分布特征及影响因素――以鄂尔多斯盆地姬塬油田延长组长8油层组为例", 《石油学报》, vol. 40, no. 05, pages 557 - 567 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116859473A (zh) * | 2023-05-23 | 2023-10-10 | 贵州省油气勘查开发工程研究院 | 面向非常规储层含水饱和度的形态学和分形学预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110717249B (zh) | 页岩气储层测井孔隙度快速预测方法及系统 | |
| CN109932297B (zh) | 一种致密砂岩储层渗透率的计算方法 | |
| CN109163997B (zh) | 一种基于声谱图深度学习的岩石表面强度测定方法 | |
| CN111425193A (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
| CN104932027A (zh) | 基于核磁共振测井的储层分类方法 | |
| CN110070220A (zh) | 一种基于神经网络算法的氨氮指标软测量方法 | |
| CN113253354A (zh) | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 | |
| CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
| CN110222981B (zh) | 一种基于参数二次选取的储层分类评价方法 | |
| CN106872669B (zh) | 一种花岗岩潜山裂缝储层电性分类方法 | |
| CN103117823A (zh) | 一种短波信道模型建模方法 | |
| CN105095552B (zh) | 一种通用概率分布参数估计流程与分析方法 | |
| CN115906624A (zh) | 危化品气体光谱生成方法、终端设备及存储介质 | |
| CN113818867B (zh) | 一种构建伪毛管压力曲线方法、系统、介质、设备及应用 | |
| CN105528657A (zh) | 基于北斗和向量机的建筑物震害预测方法 | |
| CN108898321B (zh) | 一种基于语义模板的制造技术问题标准冲突参数获取方法 | |
| CN114428043A (zh) | 多孔介质孔径分布表征方法及电子设备 | |
| CN114062190B (zh) | 表征页岩油微观赋存特征的方法、装置、终端及存储介质 | |
| CN110567905A (zh) | 一种水质检测方法及装置 | |
| CN103514288A (zh) | 客户端类别识别方法和系统 | |
| CN109342469B (zh) | 一种检测混凝土脱空深度的方法、装置和终端设备 | |
| CN112924360B (zh) | 一种定量评价岩心均质程度的方法 | |
| CN117113853B (zh) | 一种提高活塞机运行稳定性的方法及系统 | |
| CN114352257B (zh) | 岩心渗透率的表征方法及装置 | |
| CN113642656B (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 |