CN115421205A - 火成岩核磁共振测井t2谱束缚水截止值的确定方法及装置 - Google Patents
火成岩核磁共振测井t2谱束缚水截止值的确定方法及装置 Download PDFInfo
- Publication number
- CN115421205A CN115421205A CN202211047615.0A CN202211047615A CN115421205A CN 115421205 A CN115421205 A CN 115421205A CN 202211047615 A CN202211047615 A CN 202211047615A CN 115421205 A CN115421205 A CN 115421205A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- magnetic resonance
- nuclear magnetic
- igneous rock
- echo data
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 184
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 147
- 238000001228 spectrum Methods 0.000 title claims abstract description 135
- 238000000034 method Methods 0.000 title claims abstract description 71
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 67
- 239000011148 porous material Substances 0.000 claims abstract description 77
- 238000004088 simulation Methods 0.000 claims abstract description 39
- 238000012937 correction Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 11
- 239000012530 fluid Substances 0.000 claims description 23
- 238000004891 communication Methods 0.000 claims description 13
- 238000011156 evaluation Methods 0.000 abstract description 7
- 230000035699 permeability Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 4
- 239000010410 layer Substances 0.000 description 10
- 239000007789 gas Substances 0.000 description 9
- 230000007547 defect Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000005415 magnetization Effects 0.000 description 4
- 230000005298 paramagnetic effect Effects 0.000 description 4
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000003921 oil Substances 0.000 description 3
- 238000005295 random walk Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000005653 Brownian motion process Effects 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000005537 brownian motion Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000005291 magnetic effect Effects 0.000 description 2
- 239000002344 surface layer Substances 0.000 description 2
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000013589 supplement Substances 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/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种火成岩核磁共振测井T2谱束缚水截止值的确定方法及装置。该方法基于构建的火成岩孔隙网络模型，进行不同温压条件下核磁共振模拟，得到不同温压条件下核磁共振模拟回波数据，建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系，根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，对校正前后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。本发明可以连续获取井下核磁共振T2谱的束缚水截止值，提高了核磁共振在孔隙度、渗透率、饱和度评价等方面的计算精度。
Description
技术领域
本发明涉及油气勘探技术领域，具体涉及一种火成岩核磁共振测井T2谱束缚水截止值的确定方法及装置。
背景技术
T2截止值在核磁共振测井中是用于区分可动流体与束缚流体的界限，是束缚水饱和度和渗透率评价的重要参数。目前测井确定T2截止值的方法有两种：一种是固定T2截止值，砂岩是33ms，碳酸盐岩92ms；另一种是通过统计岩心核磁共振T2谱的参数、形态建立的变T2截止值法。固定T2截止法不考虑岩性、流体的影响，变T2截止值法与岩心分布区域有关，且人为因素影响大，在使用上均存在较大的限制。火成岩中含有大量的顺磁性物质Fe，T2分布受内部梯度磁场、温度的影响大，流体性质识别难度大，流体组分较砂岩、碳酸盐岩更加难以区分，目前火成岩确定T2截止值的方法仅限于核磁共振物理实验，利用离心前后测量岩心的核磁共振T2谱得到，费时费力，难以应用于全井段分析，而火成岩储层岩性非均质性强，应用效果亦不佳。
发明内容
鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的火成岩核磁共振测井T2谱束缚水截止值的确定方法及装置。
根据本发明的一个方面，提供了一种火成岩核磁共振测井T2谱束缚水截止值的确定方法，包括：
获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型；
基于火成岩孔隙网络模型进行不同温压条件下核磁共振数值模拟，得到不同温压条件下核磁共振模拟回波数据；
建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系，根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，得到校正后的火成岩核磁共振测井回波数据；
对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱；
将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
根据本发明的另一方面，提供了一种火成岩核磁共振测井T2谱束缚水截止值的确定装置，包括：
建模模块，适于获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型；
模拟模块，适于基于火成岩孔隙网络模型进行不同温压条件下核磁共振数值模拟，得到不同温压条件下核磁共振模拟回波数据；
刻度模块，适于建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系；
校正模块，适于根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，得到校正后的火成岩核磁共振测井回波数据；
反演模块，适于对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱；
截止值确定模块，适于将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
根据本发明的又一方面，提供了一种计算设备，包括：处理器、存储器、通信接口和通信总线，处理器、存储器和通信接口通过通信总线完成相互间的通信；
存储器用于存放至少一可执行指令，可执行指令使处理器执行上述火成岩核磁共振测井T2谱束缚水截止值的确定方法对应的操作。
根据本发明的再一方面，提供了一种计算机存储介质，存储介质中存储有至少一可执行指令，可执行指令使处理器执行如上述火成岩核磁共振测井T2谱束缚水截止值的确定方法对应的操作。
根据本发明提供的方案，依据的原理是顺磁物质主导的表面弛豫率不受温度、压力的影响，因此火成岩孔隙中束缚水的弛豫时间不受温压影响，而孔隙中可动流体的弛豫时间受温压影响较大。本发明采用数值模拟与实验测量相结合的方式，模拟不同火成岩在不同温压条件下的磁化衰减曲线，弥补了物理实验样本的不足，有效解决了因实验样品不足致使校正方法普适性低的问题。此外利用井下、地面不同温压条件下火成岩孔隙中束缚水、可动流体弛豫时间随温压变化的差异，来获取火成岩核磁共振测井T2谱截止值，有效解决了当前核磁共振测井确定T2截止值方法的局限性，本发明可以连续获取井下核磁共振T2谱的束缚水截止值，提高了核磁共振在孔隙度、渗透率、饱和度评价等方面的计算精度。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1A示出了根据本发明一个实施例的火成岩核磁共振测井T2谱束缚水截止值的确定方法的流程示意图；
图1B为火成岩核磁共振测井T2谱的束缚水截止值的示意图；
图2示出了根据本发明一个实施例的火成岩核磁共振测井T2谱束缚水截止值的确定装置的结构示意图；
图3示出了根据本发明一个实施例的计算设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
图1A示出了根据本发明一个实施例的火成岩核磁共振测井T2谱束缚水截止值的确定方法的流程示意图。如图1A所示，该方法包括以下步骤：
步骤S101，获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型。
通过以往采集的岩心和岩屑分析可知，火成岩岩性种类多样，垂向非均质性强。而采集的岩心往往集中于目的层，不能代表整个井段的岩性变化，因此还需要进行数值模拟研究。
数值模拟研究分为两步，第一步需要先构建符合实际地层的火成岩孔隙网络模型，具体地，获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型，其中，目标区域是火成岩研究目标区。
例如，对火成岩岩心进行CT扫描，CT技术可以无损检测岩石的内部组分及结构，提供岩石的高分辨率图像，将岩石灰度图像数值化，构建三维火成岩孔隙网络模型。从表现形式上看，数字岩心呈现的是孔隙空间和基质，孔隙网络模型反映的是孔隙空间的连通关系。以数字岩心及其孔隙中轴线为基础，可以快速进行孔隙和喉道的识别、确定孔隙大小和长度、确定喉道长度和喉道半径，基于此建立与数字岩心孔隙空间等价拓扑结构的孔隙网络模型。
在本发明一种可选实施方式中，可以通过如下方法来获取岩石孔喉参数：对目标区域的目标层位的火成岩岩心的组分及结构进行检测，得到岩石孔喉参数，其中，所述岩石孔喉参数用于表示目标区域岩石孔喉特征。
步骤S102，基于火成岩孔隙网络模型进行不同温压条件下核磁共振数值模拟，得到不同温压条件下核磁共振模拟回波数据。
数值模拟第二步是通过随机游走或有限元等方法实现核磁共振响应的模拟，具体地，在获得符合实际地层的火成岩孔隙网络模型后，通过调整模型的孔、喉等参数即可获取不同孔喉的孔隙网络模型；利用随机游走、有限元等方法，给孔隙网络模型赋予温度、压力、氢质子等不同的参数值，即可进行不同温压条件下核磁共振数值模拟实验，得到不同温压条件下核磁共振模拟回波数据。
例如，随机游走算法是一种基于对分子布朗运动现象的模拟方法。核磁共振测井检测的是地层中氢核的弛豫信号，该方法首先会将给定数量的氢质子随机分配在已经建立好的三维孔隙空间中(即孔隙网络模型)，给定氢质子扩散系数、弛豫率等参数，模拟氢质子在孔隙空间运动过程中的磁化衰减，从而得到一条类似核磁共振CPMG自旋回波串的弛豫衰减曲线，它是孔隙网络模型中氢质子的信号总和，将其反演得到核磁共振T2谱后，与真实火成岩岩心核磁共振实验得到的T2谱进行对照，即可验证该数值模拟方法能否呈现真实岩石孔隙空间分布及氢核信号弛豫衰减规律。
通过进行不同温压条件下核磁共振数值模拟，得到了大批量的核磁共振模拟回波数据，从而实现了数据补充。
步骤S103，建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系，根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，得到校正后的火成岩核磁共振测井回波数据。
由步骤S101、S102得到不同温压条件下核磁共振模拟回波数据，根据回波数据与温度、压力的对应关系，可建立不同温压条件下核磁共振回波数据与温度、压力的刻度关系。
在本发明一种可选实施方式中，火成岩孔隙流体性质包括：油、水、气，由于油、水的弛豫信号主要受温度影响，压力影响较小，而气的弛豫信号同时受温度、压力的影响，因此，针对不同的火成岩孔隙流体性质，可以建立对应的刻度关系，即，针对任一火成岩孔隙流体性质，建立不同温压条件下目标核磁共振模拟回波数据与温度压力的刻度关系，其中，所述刻度关系的校正系数由数据拟合得到，所述火成岩孔隙流体性质包括：油、水、气。
所述刻度关系具体包括：油、水层火成岩核磁共振温压校正公式及气层火成岩核磁共振温压校正公式。其中，油层、水层火成岩核磁共振温压校正公式为：
气层火成岩核磁共振温压校正公式为：
其中，M0为井下火成岩核磁共振测井回波数据、M为校正后的火成岩核磁共振测井回波数据、Tref为实验室温度、T为井下温度、Pref为实验室压力、P为井下压力，a、b、c为系数。
所述刻度关系是一种可将井下连续段的核磁共振测井回波数据校正至任一温压条件下的核磁共振回波数据的校正公式。
一般地，通常将井下温压条件校正至地面，即实验室温度压力为20℃、一个标准大气压。具体做法是：利用井下设备测量当前井深温压(或利用公式估算当前深度温压)，将当前深度温压的火成岩核磁共振回波数据输入至形成的温压校正公式中，即可计算得到地面温压条件下校正后的火成岩核磁共振测井回波数据。此时会存在井下火成岩核磁共振测井回波数据及校正后的火成岩核磁共振测井回波数据。
步骤S104，对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱。
火成岩井下核磁共振测井采集得到井下火成岩核磁共振测井回波数据，经过所建立的刻度关系(一种校正方法)校正后得到校正后的火成岩核磁共振测井回波数据，后续为了方便确定火成岩核磁共振测井T2谱的束缚水截止值(即，T2截止值)，这里可以通过反演算法对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱(其中，第一T2谱是对井下火成岩核磁共振测井回波数据进行反演处理得到的T2谱，又可称为井下核磁共振T2谱)、第二T2谱(其中，第二T2谱是对校正后的火成岩核磁共振测井回波数据反演处理得到的T2谱，又可称为校正核磁共振T2谱)，具体地反演过程这里不再赘述。
步骤S105，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
孔隙中流体分为束缚水和可动流体两部分，其中可动流体弛豫时间由自由弛豫主导，受温度压力的影响；束缚水弛豫时间由表面弛豫主导，而表面弛豫又分为两部分：一部分是表面层质子与电子异核偶极作用导致的弛豫，不受温度影响，另一部分是表面层质子与其他表面二氧化硅作用导致的弛豫，受温度影响。火成岩束缚水的弛豫时间主要由第一部分主导，表面弛豫基本不会随温度改变而改变(见图1B)。
具体地，利用上述理论获取火成岩核磁共振截止值的具体做法是：将第一T2谱与第二T2谱投影在同一坐标系中，将第一T2谱与第二T2谱进行比较，以确定第一T2谱与第二T2谱之间重合的短弛豫部分，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。如图1B所示。
由于流体分子的布朗运动以及仪器死时间的不同，校正前后所生成的第一T2谱与第二T2谱的束缚水部分的峰值可能不完全一致，因此，可以通过如下方法来判定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分是否重合，进而确定束缚水截止值：
判断所述第一T2谱的短弛豫部分的峰值与所述第二T2谱的短弛豫部分的峰值之间的差值是否小于或等于预设阈值；
若是，则判定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分之间重合，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
具体地，计算第一T2谱的短弛豫部分的峰值与第二T2谱的短弛豫部分的峰值之间的差值，将差值与预设阈值进行比较，其中，预设阈值可以根据实际经验设定，例如，设置为5％，因此，当差值小于或等于5％时，可以确定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分之间基本重合；若差值大于5％时，可以确定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分之间不重合，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
在本发明一种可选实施方式中，该方法还可以计算束缚水孔隙度，具体地，可以将重合的短弛豫部分的T2谱的面积进行累加，得到束缚水孔隙度，例如，可以利用积分方法来计算得到重合的短弛豫部分的T2谱的覆盖面积或者将重合的短弛豫部分的T2谱划分成单位面积为1的小矩形，然后，通过统计小矩形的数量来得到重合的短弛豫部分的T2谱的覆盖面积，该重合的短弛豫部分的T2谱的覆盖面积即为束缚水孔隙度。当存在确定全井段核磁共振测井的T2截止值的需求时，可以利用上述方法逐点获取全井段核磁共振测井T2截止值。
该方法得到的T2截止值与岩心核磁实验得到T2截止值较为接近，且用于全井段的T2截止值的计算，能够很好的解决火成岩储层非均质性强、T2截止值难以确定的难题，该方法在潜山几十口井中得到应用证实，在潜山储层参数评价、储层有效性评价和产能预测中均发挥了重要作用。
本发明提供的方案，依据的原理是顺磁物质主导的表面弛豫率不受温度、压力的影响，因此火成岩孔隙中束缚水的弛豫时间不受温压影响，而孔隙中可动流体的弛豫时间受温压影响较大。本发明采用数值模拟与实验测量相结合的方式，模拟不同火成岩在不同温压条件下的磁化衰减曲线，弥补了物理实验样本的不足，有效解决了因实验样品不足致使校正方法普适性低的问题。此外利用井下、地面不同温压条件下火成岩孔隙中束缚水、可动流体弛豫时间随温压变化的差异，来获取火成岩核磁共振测井T2谱截止值，有效解决了当前核磁共振测井确定T2截止值方法的局限性，本发明可以连续获取井下核磁共振T2谱的束缚水截止值，提高了核磁共振在孔隙度、渗透率、饱和度评价等方面的计算精度。
图2示出了根据本发明一个实施例的火成岩核磁共振测井T2谱束缚水截止值的确定装置的结构示意图。如图2所示，该装置包括：建模模块201、模拟模块202、刻度模块203、校正模块204、反演模块205、截止值确定模块206。
建模模块201，适于获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型；
模拟模块202，适于基于火成岩孔隙网络模型进行不同温压条件下核磁共振数值模拟，得到不同温压条件下核磁共振模拟回波数据；
刻度模块203，适于建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系；
校正模块204，适于根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，得到校正后的火成岩核磁共振测井回波数据；
反演模块205，适于对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱；
截止值确定模块206，适于将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
可选地，截止值模块进一步适于：判断所述第一T2谱的短弛豫部分的峰值与所述第二T2谱的短弛豫部分的峰值之间的差值是否小于或等于预设阈值；
若是，则判定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分之间重合，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
可选地，刻度模块进一步适于：针对任一火成岩孔隙流体性质，建立不同温压条件下目标核磁共振模拟回波数据与温度压力的刻度关系，其中，所述刻度关系的校正系数由数据拟合得到，所述火成岩孔隙流体性质包括：油、水、气。
可选地，所述刻度关系是一种将井下连续段的核磁共振测井回波数据校正至任一温压条件下的核磁共振回波数据的校正公式。
可选地，所述刻度关系具体包括：油、水层火成岩核磁共振温压校正公式及气层火成岩核磁共振温压校正公式；
其中，油层、水层火成岩核磁共振温压校正公式为：
气层火成岩核磁共振温压校正公式为：
其中，M0为井下火成岩核磁共振测井回波数据、M为校正后的火成岩核磁共振测井回波数据、Tref为实验室温度、T为井下温度、Pref为实验室压力、P为井下压力，a、b、c为系数。
可选地，建模模块进一步适于：对目标区域的目标层位的火成岩岩心的组分及结构进行检测，得到岩石孔喉参数，其中，所述岩石孔喉参数用于表示目标区域岩石孔喉特征。
本发明提供的方案，依据的原理是顺磁物质主导的表面弛豫率不受温度、压力的影响，因此火成岩孔隙中束缚水的弛豫时间不受温压影响，而孔隙中可动流体的弛豫时间受温压影响较大。本发明采用数值模拟与实验测量相结合的方式，模拟不同火成岩在不同温压条件下的磁化衰减曲线，弥补了物理实验样本的不足，有效解决了因实验样品不足致使校正方法普适性低的问题。此外利用井下、地面不同温压条件下火成岩孔隙中束缚水、可动流体弛豫时间随温压变化的差异，来获取火成岩核磁共振测井T2谱截止值，有效解决了当前核磁共振测井确定T2截止值方法的局限性，本发明可以连续获取井下核磁共振T2谱的束缚水截止值，提高了核磁共振在孔隙度、渗透率、饱和度评价等方面的计算精度。
本申请实施例还提供了一种非易失性计算机存储介质，所述计算机存储介质存储有至少一可执行指令，该计算机可执行指令可执行上述任意方法实施例中的火成岩核磁共振测井T2谱束缚水截止值的确定方法。
图3示出了根据本发明一个实施例的计算设备的结构示意图，本发明具体实施例并不对计算设备的具体实现做限定。
如图3所示，该计算设备可以包括：处理器(processor)302、通信接口(Communications Interface)304、存储器(memory)306、以及通信总线308。
其中：
处理器302、通信接口304、以及存储器306通过通信总线308完成相互间的通信。
通信接口304，用于与其它设备比如客户端或其它服务器等的网元通信。
处理器302，用于执行程序310，具体可以执行上述火成岩核磁共振测井T2谱束缚水截止值的确定方法实施例中的相关步骤。
具体地，程序310可以包括程序代码，该程序代码包括计算机操作指令。
处理器302可能是中央处理器CPU，或者是特定集成电路ASIC(ApplicationSpecific Integrated Circuit)，或者是被配置成实施本发明实施例的一个或多个集成电路。计算设备包括的一个或多个处理器，可以是同一类型的处理器，如一个或多个CPU；也可以是不同类型的处理器，如一个或多个CPU以及一个或多个ASIC。
存储器306，用于存放程序310。存储器306可能包含高速RAM存储器，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
程序310具体可以用于使得处理器302执行上述任意方法实施例中的火成岩核磁共振测井T2谱束缚水截止值的确定方法。程序310中各步骤的具体实现可以参见上述火成岩核磁共振测井T2谱束缚水截止值的确定实施例中的相应步骤和单元中对应的描述，在此不赘述。所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的设备和模块的具体工作过程，可以参考前述方法实施例中的对应过程描述，在此不再赘述。
在此提供的算法或显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述，构造这类系统所要求的结构是显而易见的。此外，本发明实施例也不针对任何特定编程语言。应当明白，可以利用各种编程语言实现在此描述的本发明的内容，并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
类似地，应当理解，为了精简本发明并帮助理解各个发明方面中的一个或多个，在上面对本发明的示例性实施例的描述中，本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该公开的方法解释成反映如下意图：即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如下面的权利要求书所反映的那样，发明方面在于少于前面公开的单个实施例的所有特征。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解，可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件，以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外，本领域的技术人员能够理解，尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在下面的权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现，或者以在一个或者多个处理器上运行的软件模块实现，或者以它们的组合实现。本领域的技术人员应当理解，可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如，计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上，或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到，或者在载体信号上提供，或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤，除有特殊说明外，不应理解为对执行顺序的限定。
Claims (10)
1.一种火成岩核磁共振测井T2谱束缚水截止值的确定方法，包括：
获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型；
基于火成岩孔隙网络模型进行不同温压条件下核磁共振数值模拟，得到不同温压条件下核磁共振模拟回波数据；
建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系，根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，得到校正后的火成岩核磁共振测井回波数据；
对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱；
将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
2.根据权利要求1所述的方法，其中，所述将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值进一步包括：
判断所述第一T2谱的短弛豫部分的峰值与所述第二T2谱的短弛豫部分的峰值之间的差值是否小于或等于预设阈值；
若是，则判定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分之间重合，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
3.根据权利要求1所述的方法，其中，所述建立不同温压条件下目标核磁共振模拟回波数据与温度压力的刻度关系进一步包括：
针对任一火成岩孔隙流体性质，建立不同温压条件下目标核磁共振模拟回波数据与温度压力的刻度关系，其中，所述刻度关系的校正系数由数据拟合得到，所述火成岩孔隙流体性质包括：油、水、气。
4.根据权利要求1-3中任一项所述的方法，其中，所述刻度关系是一种将井下连续段的核磁共振测井回波数据校正至任一温压条件下的核磁共振回波数据的校正公式。
5.根据权利要求1或3所述的方法，其中，所述刻度关系具体包括：油、水层火成岩核磁共振温压校正公式及气层火成岩核磁共振温压校正公式；
其中，油层、水层火成岩核磁共振温压校正公式为：
气层火成岩核磁共振温压校正公式为：
其中，M0为井下火成岩核磁共振测井回波数据、M为校正后的火成岩核磁共振测井回波数据、Tref为实验室温度、T为井下温度、Pref为实验室压力、P为井下压力，a、b、c为系数。
6.根据权利要求1-3中任一项所述的方法，其中，所述获取目标区域火成岩的岩石孔喉参数进一步包括：
对目标区域的目标层位的火成岩岩心的组分及结构进行检测，得到岩石孔喉参数，其中，所述岩石孔喉参数用于表示目标区域岩石孔喉特征。
7.一种火成岩核磁共振测井T2谱束缚水截止值的确定装置，包括：
建模模块，适于获取目标区域火成岩的岩石孔喉参数，根据岩石孔喉参数构建符合实际地层的火成岩孔隙网络模型；
模拟模块，适于基于火成岩孔隙网络模型进行不同温压条件下核磁共振数值模拟，得到不同温压条件下核磁共振模拟回波数据；
刻度模块，适于建立不同温压条件下核磁共振模拟回波数据与温度压力的刻度关系；
校正模块，适于根据所建立的刻度关系对井下火成岩核磁共振测井回波数据进行校正，得到校正后的火成岩核磁共振测井回波数据；
反演模块，适于对井下火成岩核磁共振测井回波数据、校正后的火成岩核磁共振测井回波数据分别进行反演处理，得到第一T2谱、第二T2谱；
截止值确定模块，适于将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
8.根据权利要求7所述的装置，其中，截止值模块进一步适于：判断所述第一T2谱的短弛豫部分的峰值与所述第二T2谱的短弛豫部分的峰值之间的差值是否小于或等于预设阈值；
若是，则判定第一T2谱的短弛豫部分与第二T2谱的短弛豫部分之间重合，将第一T2谱与第二T2谱中重合的短弛豫部分对应的最大弛豫时间确定为火成岩核磁共振测井T2谱的束缚水截止值。
9.一种计算设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如权利要求1-6中任一项所述的火成岩核磁共振测井T2谱束缚水截止值的确定方法对应的操作。
10.一种计算机存储介质，所述存储介质中存储有至少一可执行指令，所述可执行指令使处理器执行如权利要求1-6中任一项所述的火成岩核磁共振测井T2谱束缚水截止值的确定方法对应的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211047615.0A CN115421205A (zh) | 2022-08-29 | 2022-08-29 | 火成岩核磁共振测井t2谱束缚水截止值的确定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211047615.0A CN115421205A (zh) | 2022-08-29 | 2022-08-29 | 火成岩核磁共振测井t2谱束缚水截止值的确定方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115421205A true CN115421205A (zh) | 2022-12-02 |
Family
ID=84200086
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211047615.0A Pending CN115421205A (zh) | 2022-08-29 | 2022-08-29 | 火成岩核磁共振测井t2谱束缚水截止值的确定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115421205A (zh) |
-
2022
- 2022-08-29 CN CN202211047615.0A patent/CN115421205A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9070049B2 (en) | Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties | |
| US10909281B2 (en) | History matching of hydrocarbon production from heterogenous reservoirs | |
| WO2017028161A1 (en) | Nmr anaylysis system and method for porous media | |
| CN103513285A (zh) | 一种确定横向表面弛豫速率的方法及装置 | |
| CN113919219A (zh) | 基于测井大数据的地层评价方法及系统 | |
| CN107688037B (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
| CN106323836A (zh) | 一种井壁渗透率计算方法及装置 | |
| WO2022226332A1 (en) | Methods and systems for determining reservoir and fracture properties | |
| CN114283254B (zh) | 基于核磁共振数据的岩心数字化孔隙网络模型构建方法 | |
| CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
| NO20120654A1 (no) | Fremgangsmåte og anordning for inkorporering av intern gradient og begrenset diffusjon i NMR-inversjon ved brønnlogging | |
| WO2014142976A1 (en) | Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties | |
| CN116776513A (zh) | 地热资源量计算方法、装置、存储介质及电子设备 | |
| CN110232208A (zh) | 用于识别优势渗流段的储层参数模型构建方法及设备 | |
| CN110046414A (zh) | 用于识别优势渗流段的储层参数模型构建方法及设备 | |
| CN111485874B (zh) | 确定水平井储层含水饱和度的方法及装置 | |
| CN115421205A (zh) | 火成岩核磁共振测井t2谱束缚水截止值的确定方法及装置 | |
| CN114185108B (zh) | 裂缝孔隙度确定方法、存储介质和计算机设备 | |
| CN108979629B (zh) | 一种基于密度和中子和核磁共振测井的气层密度计算方法 | |
| CN115452874A (zh) | 基于核磁共振对火山岩进行测量的方法、设备及存储介质 | |
| CN113781443A (zh) | 一种基于数字图像颜色梯度的岩层裂隙分布规律获取方法、系统、终端及可读存储介质 | |
| CN107818185B (zh) | 一种获取致密砂岩地层因素的方法 | |
| CN116819644B (zh) | 页岩油储层最小水平主应力的确定方法及装置 | |
| CN112711076B (zh) | 提取石油钻井中泥浆侵入地层深度的方法和装置 | |
| CN117192607A (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 |