CN111965719B - 一种弛豫时间测量方法和装置 - Google Patents
一种弛豫时间测量方法和装置 Download PDFInfo
- Publication number
- CN111965719B CN111965719B CN202010707548.5A CN202010707548A CN111965719B CN 111965719 B CN111965719 B CN 111965719B CN 202010707548 A CN202010707548 A CN 202010707548A CN 111965719 B CN111965719 B CN 111965719B
- Authority
- CN
- China
- Prior art keywords
- echo
- frequency
- polarization
- relaxation time
- period
- 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
- 238000000691 measurement method Methods 0.000 title claims abstract description 20
- 230000010287 polarization Effects 0.000 claims abstract description 88
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000010586 diagram Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
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
-
- 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
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请实施例公开了一种弛豫时间测量方法和装置，该方法包括：在井下核磁测井仪器工作时，采用至少两个不同频率进行工作；两个不同频率包括第一频率和第二频率；在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集；所述第一极化时间段是第一频率对应的极化时间段；所述第二极化时间段是第二频率对应的极化时间段；根据第一回波和第二回波计算纵向弛豫时间T1和横向弛豫时间T2。通过该实施例方案，提高了T1‑T2测量效率，缩短了测量周期，提高了纵向分辨率。
Description
技术领域
本文涉及核磁测井技术，尤指一种弛豫时间测量方法和装置。
背景技术
目前，相关的井下T1-T2测量模式均是通过采用国外三大石油服务公司常用的测量模式。
相关测量模式中，一些测量模式虽然测量周期短，分辨率较高，但是反应T1-T2谱时因极化等待时间差异大，无法很好地形成T1-T2谱。还有一些测量模式虽然极化时长分布较大，但测量周期较长，纵向分辨率较低。
因此，相关测量模式测量精度有限，急需研发高分辨率T1-T2测量模式。
发明内容
本申请实施例提供了一种弛豫时间测量测量方法和装置，能够提高T1-T2测量效率，缩短测量周期，提高纵向分辨率。
本申请实施例提供了一种弛豫时间测量方法，所述方法可以包括：
在井下核磁测井仪器工作时，采用至少两个不同频率进行工作；所述两个不同频率包括第一频率和第二频率；
在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集；所述第一极化时间段是第一频率对应的极化时间段；所述第二极化时间段是第二频率对应的极化时间段；
根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2。
在本申请的示例性实施例中，所述第一极化时间段和所述第二极化时间段的长度范围均可以为0.02-10秒。
在本申请的示例性实施例中，所述在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集可以包括：循环执行以下步骤预设次数：
在所述第二极化时间段中对所述第一回波进行采集时停止对所述第二回波进行采集，所述第二频率处于极化等待状态；
在所述第一极化时间段中对所述第二回波进行采集时停止对所述第一回波进行采集，所述第一频率处于极化等待状态。
在本申请的示例性实施例中，所述在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集还可以包括：
在每个循环周期内，对所述第一回波进行一次采集，获取一个关于所述第一频率的回波串，对所述第二回波进行多次采集，获取多个关于所述第二频率的回波串；其中，一个循环周期包含一个第一极化时间段和一个第二极化时间段。
在本申请的示例性实施例中，所述方法还可以包括：
在井下核磁测井仪器工作时，采用多组频率进行工作；其中，每组频率中分别包括所述第一频率和所述第二频率。
在本申请的示例性实施例中，每组频率中的第一频率均相同，或至少两组中的第一频率不同；以及，
每组频率中的第二频率均相同，或至少两组中的第二频率不同。
在本申请的示例性实施例中，所述根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2可以包括：
对所述第一回波和所述第二回波进行反演计算，获取所述纵向弛豫时间T1和所述横向弛豫时间T2。
在本申请的示例性实施例中，所述方法还可以包括：构建T1-T2谱图；
所述构建T1-T2谱图可以包括：
根据获取的所述纵向弛豫时间T1和所述横向弛豫时间T2构建所述T1-T2谱图；或者，
采用关于所述第一回波和所述第二回波的第一回波信号图映射出所述T1-T2谱图。
在本申请的示例性实施例中，所述方法还可以包括：
对构建的T1-T2谱图进行回波正演，获取第二回波信号图；
根据所述第二回波信号图进行反演计算，获取模拟出的T1-T2谱图；
将所述模拟出的T1-T2谱图与所述构建的T1-T2谱图相比较，以根据比较结果检测所述模拟出的T1-T2谱图与所述构建的T1-T2谱图的相似率。
本申请实施例还提供了一种弛豫时间测量装置，可以包括处理器和计算机可读存储介质，所述计算机可读存储介质中存储有指令，当所述指令被所述处理器执行时，实现上述任意一项所述的弛豫时间测量方法。
本申请实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任意一项所述的弛豫时间测量方法。
本申请实施例可以包括在井下核磁测井仪器工作时，采用至少两个不同频率进行工作；所述两个不同频率包括第一频率和第二频率；在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集；所述第一极化时间段是第一频率对应的极化时间段；所述第二极化时间段是第二频率对应的极化时间段；根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2。通过该实施例方案，提高了T1-T2测量效率，缩短了测量周期，提高了纵向分辨率。
本申请的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本申请而了解。本申请的其他优点可通过在说明书以及附图中所描述的方案来实现和获得。
附图说明
附图用来提供对本申请技术方案的理解，并且构成说明书的一部分，与本申请的实施例一起用于解释本申请的技术方案，并不构成对本申请技术方案的限制。
图1为本申请实施例的核磁测井中的弛豫时间测量方法流程图；
图2为本申请实施例的高分辨率T1-T2模式序列顺序示意图；
图3为本申请实施例的高分辨率T1-T2模式序列采集示意图；
图4为本申请实施例的构建的T1-T2谱图示意图；
图5为本申请实施例的回波信号示意图；
图6为本申请实施例的高分辨率T1-T2模式仿真出的T1-T2谱图示意图；
图7为本申请实施例的核磁测井中的T1-T2测量装置组成框图。
具体实施方式
本申请描述了多个实施例，但是该描述是示例性的，而不是限制性的，并且对于本领域的普通技术人员来说显而易见的是，在本申请所描述的实施例包含的范围内可以有更多的实施例和实现方案。尽管在附图中示出了许多可能的特征组合，并在具体实施方式中进行了讨论，但是所公开的特征的许多其它组合方式也是可能的。除非特意加以限制的情况以外，任何实施例的任何特征或元件可以与任何其它实施例中的任何其他特征或元件结合使用，或可以替代任何其它实施例中的任何其他特征或元件。
本申请包括并设想了与本领域普通技术人员已知的特征和元件的组合。本申请已经公开的实施例、特征和元件也可以与任何常规特征或元件组合，以形成由权利要求限定的独特的发明方案。任何实施例的任何特征或元件也可以与来自其它发明方案的特征或元件组合，以形成另一个由权利要求限定的独特的发明方案。因此，应当理解，在本申请中示出和/或讨论的任何特征可以单独地或以任何适当的组合来实现。因此，除了根据所附权利要求及其等同替换所做的限制以外，实施例不受其它限制。此外，可以在所附权利要求的保护范围内进行各种修改和改变。
此外，在描述具有代表性的实施例时，说明书可能已经将方法和/或过程呈现为特定的步骤序列。然而，在该方法或过程不依赖于本文所述步骤的特定顺序的程度上，该方法或过程不应限于所述的特定顺序的步骤。如本领域普通技术人员将理解的，其它的步骤顺序也是可能的。因此，说明书中阐述的步骤的特定顺序不应被解释为对权利要求的限制。此外，针对该方法和/或过程的权利要求不应限于按照所写顺序执行它们的步骤，本领域技术人员可以容易地理解，这些顺序可以变化，并且仍然保持在本申请实施例的精神和范围内。
本申请实施例提供了一种弛豫时间测量方法，如图1所示，所述方法可以包括步骤S101-S103：
S101、在井下核磁测井仪器工作时，采用至少两个不同频率进行工作；所述两个频率包括第一频率和第二频率；
S102、在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集；所述第一极化时间段是第一频率对应的极化时间段；所述第二极化时间段是第二频率对应的极化时间段；
S103、根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2。
在本申请的示例性实施例中，所述在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集可以包括：循环执行以下步骤预设次数：
在所述第二极化时间段中对所述第一回波进行采集时停止对所述第二回波进行采集，所述第二频率处于极化等待状态(该极化等待状态，即处于等待底层中的氢原子被极化的状态，以使得氢原子有充分的极化时间)；
在所述第一极化时间段中对所述第二回波进行采集时停止对所述第一回波进行采集，所述第一频率处于极化等待状态。
在本申请的示例性实施例中，核磁共振测井是在井底条件下实现核磁共振测量。测量原理的核心之一是对地层施加外加磁场，使氢原子核磁化。氢核是一种磁性核，具有核磁矩。磁体放到井中，将在井周围地层产生磁场，使氢核的磁矩沿磁场方向取向，这个过程叫磁化或极化，极化时间常数用T1表示，称作纵向弛豫时间。T1与孔隙度的大小、孔隙直径的大小、孔隙中流体的性质以及地层的岩性等因素有关。核磁共振测井原理的核心之二是利用一个天线系统，向地层发射特定能量、特定频率和特定时间间隔的地磁波脉冲，产生所谓的自旋回波信号，并接收和采集这种回波信号，所采用的方法叫做自旋回波法。观测到的回波串为按指数规律衰减的信号，其衰减的时间常数用T2表示，叫做横向豫驰时间，他与地层孔隙度的大小。孔隙直径的大小、空隙中流体的性质、岩性以及采集参数等因素有关。
在本申请的示例性实施例中，依据核磁共振测井的原理，为了提高测井效率可以采用至少两个不同频率进行核磁测井，为了保证每个频率具有足够的极化时长，可以分别以第一频率和第二频率对应的极化时间段(第一极化时间段和第二极化时间段)为回波采集时段，在第一极化时间段内采集第二回波，第一频率处于极化等待状态，为第一频率留有足够的极化时长，在第二极化时间段内采集第一回波，第二频率处于极化等待状态，为第二频率留有足够的极化时长，通过两个频率的配合，在一个频率极化时对另一个频率的回波进行采集，因此缩短了采集周期，并保证了足够长的极化时长。
在本申请的示例性实施例中，所述第一极化时间段和所述第二极化时间段的长度范围均为0.02-10秒；例如，可以选择500毫秒。
在本申请的示例性实施例中，该极化时间段的时长分布范围较为广泛，能够很好地进行T1-T2谱的测量分析。
在本申请的示例性实施例中，所述在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集还可以包括：
在每个循环周期内，对所述第一回波进行一次采集，获取一个关于所述第一频率的回波串；对所述第二回波进行多次采集，获取多个关于所述第二频率的回波串；其中，一个循环周期包含一个第一极化时间段和一个第二极化时间段。
在本申请的示例性实施例中，对两个不同频率的回波信号(第一回波和第二回波)进行轮流采集时，每进行一次采集可以采集一组回波串，在每个采集周期内，可以仅采集一次第一回波(即采集一个回波串)，连续采集多次第二回波(即采集多个回波串)，如图2、图3所示，其中，A1为第一频率对应的第一回波的回波串，B1、B2、B3……等为第二频率对应的第二回波的回波串，其中的*2、*4等是指采集次数。
在本申请的示例性实施例中，采集的第一回波可以为主道回波(即主回波通道的回波)，该主道回波可以采集一次，获取一个回波串；采集的第二回波可以非主道回波，该非主道回波可以采集多组，获取多组回波串，每组回波串内可以包含多次采集的多个回波串。
在本申请的示例性实施例中，通过两个频率的配合，设置合理的回波采集方式，增强了测井时的纵向分辨率。
在本申请的示例性实施例中，在第一回波与第二回波之间进行采集切换时(即第一极化时间段和第二极化时间段切换时)，可以设置一定的延时，由于是对不同的回波进行采集，因此，这个延时可以越小越好，或者实现无缝切换(即没有延时)，以便进一步缩短采集时长。
在本申请的示例性实施例中，所述方法还可以包括：
在井下核磁测井仪器工作时，采用多组频率进行工作；其中，每组频率中包括所述第一频率和第二频率。
在本申请的示例性实施例中，可以不限于采用一组频率进行T1和T2测量采集，可以同时采用多组频率进行测量，以提高测量准确性。
在本申请的示例性实施例中，例如，可以采用一组为A频率和B频率，一组为C频率和D频率，一组为E频率和F频率同时测量，每一组测得的回波可以用以得出一个T1-T2图，各组之间的T1-T2图也可以比较相似性，从而确定测量准确率。
在本申请的示例性实施例中，每组频率中的第一频率均相同，或至少两组中的第一频率不同；以及，
每组频率中的第二频率均相同，或至少两组中的第二频率不同。
在本申请的示例性实施例中，所述根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2可以包括：
对所述第一回波和所述第二回波进行反演计算，获取所述纵向弛豫时间T1和所述横向弛豫时间T2。
在本申请的示例性实施例中，所述方法还可以包括：构建T1-T2谱图；
所述构建T1-T2谱图可以包括：
根据获取的所述纵向弛豫时间T1和所述横向弛豫时间T2构建所述T1-T2谱图；或者，
采用关于所述第一回波和所述第二回波的第一回波信号图映射出所述T1-T2谱图。
在本申请的示例性实施例中，所采用的反演算法可以为目前任意可以使用的反演算法，对于详细算法不做限制。
在本申请的示例性实施例中，根据图2和图3可以映射出T1-T2图谱。可以通过将图2和图3向左映射获得图4中所示的纵向弛豫时间T1，可以通过将图2和图3向下映射获得图4中所示的横向弛豫时间T2。在本申请的示例性实施例中，所述方法还可以包括：
对构建的T1-T2谱图进行回波正演，获取第二回波信号图；
根据所述第二回波信号图进行反演计算，获取模拟出的T1-T2谱图；
将所述模拟出的T1-T2谱图与所述构建的T1-T2谱图相比较，以根据比较结果检测所述模拟出的T1-T2谱图与所述构建的T1-T2谱图的相似率。
在本申请的示例性实施例中，所采用的正演算法和反演算法可以为目前任意可以使用的正演算法和反演算法，对于详细算法不做限制。
在本申请的示例性实施例中，对构建的T1-T2谱图进行回波正演可以获取如图5所示的回波信号图；反演结果在吻合的信噪比范围内。通过对回波信号图中的回波信号进行反演分析可以得到图6所示的高分辨率T1-T2模式仿真出的T1-T2谱图，通过将图4中仿真出的T1-T2谱图和图6中构建的T1-T2谱图相比较可知，两者的相似率极高。
在本申请的示例性实施例中，本申请实施例方案基于核磁常用的回波采集序列，设计了不同频率的不同极化及采集时间：
1、通过两个频率的配合，缩短了采集周期，增强了测井时的纵向分辨率；
2、预设的极化时长分布较为广泛，能够很好地进行T1-T2谱的测量分析。
本申请实施例还提供了一种弛豫时间测量装置1，如图7所示，可以包括处理器11和计算机可读存储介质12，所述计算机可读存储介质12中存储有指令，当所述指令被所述处理器11执行时，实现上述任意一项所述的弛豫时间测量方法。
本申请实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任意一项所述的弛豫时间测量方法。
在本申请的示例性实施例中，前述的方法实施例中的任意实施例均适用于该装置实施例和计算机可读存储介质实施例，在此不再一一赘述。
本领域普通技术人员可以理解，上文中所公开方法中的全部或某些步骤、系统、装置中的功能模块/单元可以被实施为软件、固件、硬件及其适当的组合。在硬件实施方式中，在以上描述中提及的功能模块/单元之间的划分不一定对应于物理组件的划分；例如，一个物理组件可以具有多个功能，或者一个功能或步骤可以由若干物理组件合作执行。某些组件或所有组件可以被实施为由处理器，如数字信号处理器或微处理器执行的软件，或者被实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读介质上，计算机可读介质可以包括计算机存储介质(或非暂时性介质)和通信介质(或暂时性介质)。如本领域普通技术人员公知的，术语计算机存储介质包括在用于存储信息(诸如计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机存储介质包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD-ROM、数字多功能盘(DVD)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。此外，本领域普通技术人员公知的是，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何信息递送介质。
Claims (8)
1.一种弛豫时间测量方法，其特征在于，所述方法包括：
在井下核磁测井仪器工作时，采用至少两个不同频率进行工作；所述两个不同频率包括第一频率和第二频率；
在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集；所述第一极化时间段是第一频率对应的极化时间段；所述第二极化时间段是第二频率对应的极化时间段；
根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2,
其中,所述在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集还包括：
在每个循环周期内，对所述第一回波进行一次采集，获取一个关于所述第一频率的回波串，对所述第二回波进行多次连续采集，获取多个关于所述第二频率的回波串；其中，一个循环周期包含一个第一极化时间段和一个第二极化时间段；
所述第一回波为主道回波，所述第二回波为非主道回波；
所述在预设的第二极化时间段中对所述第一频率对应的第一回波进行采集，在预设的第一极化时间段中对所述第二频率对应的第二回波进行采集包括：循环执行以下步骤预设次数：
在所述第二极化时间段中对所述第一回波进行采集时停止对所述第二回波进行采集，所述第二频率处于极化等待状态；
在所述第一极化时间段中对所述第二回波进行采集时停止对所述第一回波进行采集，所述第一频率处于极化等待状态。
2.根据权利要求1所述的弛豫时间测量方法，其特征在于，所述第一极化时间段和所述第二极化时间段的长度范围均为0.02-10秒。
3.根据权利要求1-2任意一项所述的弛豫时间测量方法，其特征在于，所述方法还包括：
在井下核磁测井仪器工作时，采用多组频率进行工作；其中，每组频率中分别包括所述第一频率和所述第二频率。
4.根据权利要求3所述的弛豫时间测量方法，其特征在于，
每组频率中的第一频率均相同，或至少两组中的第一频率不同；以及，
每组频率中的第二频率均相同，或至少两组中的第二频率不同。
5.根据权利要求1-2任意一项所述的弛豫时间测量方法，其特征在于，所述根据所述第一回波和所述第二回波计算纵向弛豫时间T1和横向弛豫时间T2包括：对所述第一回波和所述第二回波进行反演计算，获取所述纵向弛豫时间T1和所述横向弛豫时间T2。
6.根据权利要求1-2任意一项所述的弛豫时间测量方法，其特征在于，所述方法还包括：构建T1-T2谱图；
所述构建T1-T2谱图包括：
根据获取的所述纵向弛豫时间T1和所述横向弛豫时间T2构建所述T1-T2谱图；或者，
采用关于所述第一回波和所述第二回波的第一回波信号图映射出所述T1-T2谱图。
7.根据权利要求6所述的弛豫时间测量方法，其特征在于，所述方法还包括：
对构建的T1-T2谱图进行回波正演，获取第二回波信号图；
根据所述第二回波信号图进行反演计算，获取模拟出的T1-T2谱图；
将所述模拟出的T1-T2谱图与所述构建的T1-T2谱图相比较，以根据比较结果检测所述模拟出的T1-T2谱图与所述构建的T1-T2谱图的相似率。
8.一种弛豫时间测量装置，包括处理器和计算机可读存储介质，所述计算机可读存储介质中存储有指令，其特征在于，当所述指令被所述处理器执行时，实现如权利要求1-7任意一项所述的弛豫时间测量方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010707548.5A CN111965719B (zh) | 2020-07-21 | 2020-07-21 | 一种弛豫时间测量方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010707548.5A CN111965719B (zh) | 2020-07-21 | 2020-07-21 | 一种弛豫时间测量方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111965719A CN111965719A (zh) | 2020-11-20 |
| CN111965719B true CN111965719B (zh) | 2024-03-15 |
Family
ID=73364096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010707548.5A Active CN111965719B (zh) | 2020-07-21 | 2020-07-21 | 一种弛豫时间测量方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111965719B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111980663B (zh) * | 2020-07-21 | 2023-08-15 | 中海油田服务股份有限公司 | 一种多频多维核磁测井方法和装置 |
| CN114324441A (zh) * | 2021-12-31 | 2022-04-12 | 苏州纽迈分析仪器股份有限公司 | 二维反演图谱生成方法、设备及存储介质 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1345848A (en) * | 1972-05-09 | 1974-02-06 | Texaco Development Corp | Nuclear magnetism logging |
| GB9908393D0 (en) * | 1998-05-11 | 1999-06-09 | Schlumberger Holdings | NMR Well logging |
| US6049205A (en) * | 1997-10-01 | 2000-04-11 | Western Atlas International, Inc. | Multiple frequency method for nuclear magnetic resonance longitudinal relaxation measurement and pulsing sequence for power use optimization |
| US6331775B1 (en) * | 1999-09-15 | 2001-12-18 | Baker Hughes Incorporated | Gas zone evaluation by combining dual wait time NMR data with density data |
| CN1763563A (zh) * | 2004-10-21 | 2006-04-26 | 北京大学 | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 |
| GB201114159D0 (en) * | 2011-08-17 | 2011-10-05 | Schlumberger Holdings | NMR flowmwter with superconducting polarizer |
| DE102010022385A1 (de) * | 2010-06-01 | 2011-12-01 | Technische Universität Darmstadt | Field-Cycling NMR Relaxometrie |
| CN102608664A (zh) * | 2012-02-17 | 2012-07-25 | 中国石油大学(北京) | 深度维核磁共振反演获取横向弛豫时间谱的方法及装置 |
| CN105259198A (zh) * | 2015-09-28 | 2016-01-20 | 中国石油大学(北京) | 二维核磁共振弛豫时间的测量方法 |
| CN105486709A (zh) * | 2015-12-17 | 2016-04-13 | 中国石油大学(北京) | 井下稠油分子链长的测量方法和装置 |
| CN108426909A (zh) * | 2018-02-11 | 2018-08-21 | 中国石油大学(北京) | 核磁共振弛豫时间和流速同步快速测量的方法及装置 |
| CN110029990A (zh) * | 2019-04-15 | 2019-07-19 | 中国海洋石油集团有限公司 | 一种核磁共振测井方法和装置 |
| CN111980663A (zh) * | 2020-07-21 | 2020-11-24 | 中海油田服务股份有限公司 | 一种多频多维核磁测井方法和装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6956371B2 (en) * | 1995-10-12 | 2005-10-18 | Halliburton Energy Services, Inc. | Method and apparatus for detecting diffusion sensitive phases with estimation of residual error in NMR logs |
| US20020175682A1 (en) * | 2001-05-23 | 2002-11-28 | Songhua Chen | Rapid nmr multi-frequency t1 and t2 acquisition for earth formations evaluation with mwd or wireline tools |
| MXPA06003671A (es) * | 2003-10-03 | 2006-06-20 | Halliburton Energy Serv Inc | Sistema y metodos para diagrafias basadas en t1. |
| AU2007343899B2 (en) * | 2007-01-18 | 2011-09-15 | Halliburton Energy Services, Inc. | Simultaneous relaxation time inversion |
-
2020
- 2020-07-21 CN CN202010707548.5A patent/CN111965719B/zh active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1345848A (en) * | 1972-05-09 | 1974-02-06 | Texaco Development Corp | Nuclear magnetism logging |
| US6049205A (en) * | 1997-10-01 | 2000-04-11 | Western Atlas International, Inc. | Multiple frequency method for nuclear magnetic resonance longitudinal relaxation measurement and pulsing sequence for power use optimization |
| GB9908393D0 (en) * | 1998-05-11 | 1999-06-09 | Schlumberger Holdings | NMR Well logging |
| US6331775B1 (en) * | 1999-09-15 | 2001-12-18 | Baker Hughes Incorporated | Gas zone evaluation by combining dual wait time NMR data with density data |
| CN1763563A (zh) * | 2004-10-21 | 2006-04-26 | 北京大学 | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 |
| DE102010022385A1 (de) * | 2010-06-01 | 2011-12-01 | Technische Universität Darmstadt | Field-Cycling NMR Relaxometrie |
| GB201114159D0 (en) * | 2011-08-17 | 2011-10-05 | Schlumberger Holdings | NMR flowmwter with superconducting polarizer |
| CN102608664A (zh) * | 2012-02-17 | 2012-07-25 | 中国石油大学(北京) | 深度维核磁共振反演获取横向弛豫时间谱的方法及装置 |
| CN105259198A (zh) * | 2015-09-28 | 2016-01-20 | 中国石油大学(北京) | 二维核磁共振弛豫时间的测量方法 |
| CN105486709A (zh) * | 2015-12-17 | 2016-04-13 | 中国石油大学(北京) | 井下稠油分子链长的测量方法和装置 |
| CN108426909A (zh) * | 2018-02-11 | 2018-08-21 | 中国石油大学(北京) | 核磁共振弛豫时间和流速同步快速测量的方法及装置 |
| CN110029990A (zh) * | 2019-04-15 | 2019-07-19 | 中国海洋石油集团有限公司 | 一种核磁共振测井方法和装置 |
| CN111980663A (zh) * | 2020-07-21 | 2020-11-24 | 中海油田服务股份有限公司 | 一种多频多维核磁测井方法和装置 |
| WO2022017120A1 (zh) * | 2020-07-21 | 2022-01-27 | 中海油田服务股份有限公司 | 多频多维核磁测井方法、装置和存储介质 |
Non-Patent Citations (10)
| Title |
|---|
| A fast field-cycling device for high-resolution NMR：Design and application to spin relaxation and hyperpolarization experiments;Kiryutin A. S.等;Joumal of Magnetic Resonance;第第263卷卷;第79-91页 * |
| A Method of choosing the optimal number of singular values in the inverse laplace transform for the twodimensional NMR distribution function;Jiang Zhimin等;Chinese Physical letters;第30卷(第01期);第1-4页 * |
| D-T_2二维核磁共振脉冲序列及反演方法改进设计;范宜仁等;《物理学报》;20151231;第64卷(第09期);第524-539页 * |
| MRIL-P核磁共振成像测井技术及其应用;张彦;《石油仪器》;20090430;第23卷(第02期);第22-28页 * |
| T_2-D二维核磁共振测井评价在气水识别中的初步应用;罗利等;《测井技术》;20141031;第38卷(第05期);第564-568页 * |
| 二维核磁共振D―T_2测量数据处理算法及应用;姜志敏等;《测井技术》;20200229;第44卷(第01期);第27-31页 * |
| 二维核磁测井仪脉冲序列发生器的设计与实现;薛志波等;电子技术与软件工程(第12期);第153-155页 * |
| 基于FPGA的核磁共振测井仪控制逻辑设计;于慧俊等;《波谱学杂志》;20120630;第29卷(第02期);第201-208页 * |
| 尉中良等.《地球物理测井》.北京：地质出版社,2005,第202-203页. * |
| 快速孔渗核磁共振测井模式设计;姜志敏等;测井技术;第44卷(第02期);第129-134页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111965719A (zh) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101680936B (zh) | 用于高信噪比nmr测井的方法与装置 | |
| CN111965719B (zh) | 一种弛豫时间测量方法和装置 | |
| CN111980663B (zh) | 一种多频多维核磁测井方法和装置 | |
| CN104597420A (zh) | 基于多次激发的磁共振扩散成像方法 | |
| CN104345289A (zh) | 用于确定水的t1 时间和脂肪的t1 时间的方法和磁共振设备 | |
| RU2432587C2 (ru) | Способ анализа данных с общими и различающимися свойствами | |
| CN110029990B (zh) | 一种核磁共振测井方法和装置 | |
| CN105182431B (zh) | 石油储层岩石的组分识别及定量评价方法 | |
| CN110082381B (zh) | 一种利用dnp-mri测量油水分布方法 | |
| US6690167B2 (en) | Nuclear magnetic resonance pulse sequence with refocusing pulses having reduced tipping angle | |
| CN104330433A (zh) | 一种获得目的储层横向弛豫时间分布的方法及装置 | |
| CN103969610A (zh) | 采集磁共振数据和确定b1磁场的方法及磁共振设备 | |
| CN105301537A (zh) | 同时磁共振成像方法和用于同时多核磁共振成像的装置 | |
| AU2017204581A1 (en) | Nuclear magnetic resonance tool calibration | |
| NO20120654A1 (no) | Fremgangsmåte og anordning for inkorporering av intern gradient og begrenset diffusjon i NMR-inversjon ved brønnlogging | |
| CN109254255A (zh) | 核磁共振横向弛豫时间谱测量方法、装置及可读存储介质 | |
| CA2465809A1 (en) | Method and apparatus for subterranean formation flow imaging | |
| CN110596772B (zh) | 一种扩散t2弛豫二维谱的校正方法和装置 | |
| CN103635825A (zh) | 交错黑血和亮血动态对比增强（dce）mri | |
| US20130169275A1 (en) | Magnetic resonance imaging method and magnetic resonance imaging device | |
| CN107144885A (zh) | 一种基于地井模式的核磁共振探测系统及方法 | |
| WO2012140543A1 (en) | Mri of chemical species having different resonance frequencies using an ultra-short echo time sequence | |
| US10302733B2 (en) | NMR spin-echo amplitude estimation | |
| CN105785457B (zh) | 基于双极性脉冲的地面核磁共振横向弛豫时间的测量方法 | |
| CN115452874A (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 |