CN110108616B - 一种对核磁共振实验中煤样离心过程信号损失补偿方法 - Google Patents

一种对核磁共振实验中煤样离心过程信号损失补偿方法 Download PDF

Info

Publication number
CN110108616B
CN110108616B CN201910333621.4A CN201910333621A CN110108616B CN 110108616 B CN110108616 B CN 110108616B CN 201910333621 A CN201910333621 A CN 201910333621A CN 110108616 B CN110108616 B CN 110108616B
Authority
CN
China
Prior art keywords
nuclear magnetic
magnetic resonance
coal sample
experiment
pore
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
Application number
CN201910333621.4A
Other languages
English (en)
Other versions
CN110108616A (zh
Inventor
刘大锰
李夏伟
蔡益栋
王颖晋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN201910333621.4A priority Critical patent/CN110108616B/zh
Publication of CN110108616A publication Critical patent/CN110108616A/zh
Application granted granted Critical
Publication of CN110108616B publication Critical patent/CN110108616B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G01N15/0893—Investigating volume, surface area, size or distribution of pores; Porosimetry by measuring weight or volume of sorbed fluid, e.g. B.E.T. method
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Dispersion Chemistry (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Abstract

本发明实施例公开了一种对核磁共振实验中煤样离心过程信号损失补偿方法，包括如下步骤：步骤100、选取不同变质程度的煤样，对饱和水煤样进行核磁共振实验测试，对所述饱和水煤样进行离心实验后获得离心煤样，对离心煤样再次进行核磁共振实验测试，计算弛豫时间截止值并作出前后两次核磁共振的孔径分布曲线；步骤200、对原始的煤样进行低温氮吸附实验测试，计算得到比表面积S和孔体积V；步骤300、基于表面弛豫率的计算，对离心实验后的煤样在核磁共振实验中的信号损失量进行补偿；本发明对离心过程产生的核磁信号量损失进行补偿，提高了分析结果的精度。

Description

一种对核磁共振实验中煤样离心过程信号损失补偿方法
技术领域
本发明实施例涉及核磁共振实验技术领域，具体涉及一种对核磁共振实验中煤样离心过程信号损失补偿方法。
背景技术
煤层气是一种非常规清洁能源，它的开发可以缓解常规化石能源的不足，且能够降低煤炭直接开采带来的一些安全隐患。煤层气储量的评价体系中，煤岩孔隙物性的研究为一重要因素，包括孔隙度、孔径分布及渗透性等。目前煤岩孔隙研究方法众多，常规实验包括压汞法、低温氮吸附实验，核磁共振实验是一种新兴的实验手段，因具有无损、快速的特点在科研、生产中被得到广泛应用。
核磁共振实验尤其是在对煤岩孔隙表征方面因测点较多而更加精确，相比其他实验方法适用性更强，不会对孔隙结构产生损坏，但是在研究孔径分布中为确保准确，必须结合离心实验进行操作，但由于煤岩较脆的特性，在离心实验中离心力过大会对不同煤体结构发育煤样产生差异性破坏，造成核磁共振实验中信号量的损失，为确保表征结果的准确性，需要采取一定措施对损失部分的核磁共振信号量进行补偿。
发明内容
为此，本发明实施例提供一种对核磁共振实验中煤样离心过程信号损失补偿方法，通过对煤样离心产生的损失进行补偿以提高测量的精度，以解决现有技术中的问题。
为了实现上述目的，本发明实施例提供如下技术方案：一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，包括如下步骤：
步骤100、选取不同变质程度的煤样，对饱和水煤样进行核磁共振实验测试，对所述饱和水煤样进行离心实验后获得离心煤样，对离心煤样再次进行核磁共振实验测试，计算弛豫时间截止值并作出前后两次核磁共振的孔径分布曲线；
步骤200、对原始的煤样进行低温氮吸附实验测试，计算得到比表面积S和孔体积V；
步骤300、基于表面弛豫率的计算，对离心实验后的煤样在核磁共振实验中的信号损失量进行补偿。
进一步地，在步骤100中，饱和水煤样和离心煤样的核磁共振实验具体步骤为：
将饱和水煤样在2.38MHz的Rec Core 2500低场核磁共振测量仪中进行实验，反演T2弛豫时间谱，实验设置参数为：回波间隔时间0.240ms，等待时间6s，回波个数8000，扫描次数64，实验温度25℃；
将离心后的煤样在相同的测试环境下再次进行核磁共振实验。
进一步地，还包括所述饱和水煤样进行离心实验前确定最佳离心力，其具体步骤为：将饱和水煤样在恒定温度25℃，离心力分别为0.69、0.92、1.15、1.38、1.61MPa的条件下依次进行离心，每次离心时间均为1h，通过在不同离心转速作用下获取的煤样核磁共振信号量来确定最佳离心力。
进一步地，核磁共振实验测试弛豫时间截止值的具体计算步骤为：
在核磁共振实验测试中直接获得饱和水煤样和离心后煤样核磁共振的弛豫时间谱图；
对离心前后的弛豫时间谱图分别做核磁共振信号累计幅度曲线，并从离心后的累计幅度曲线最大值处作X轴的平行线，直至与饱和水累计幅度曲线相交，并由交点处再作垂线与X轴相交，即所述垂线的横坐标即为弛豫时间截止值。
进一步地，在步骤200中，煤样低温氮吸附实验测试的具体步骤为：
将煤样破碎并筛选60～80目的煤粉，将筛选后的煤样进行脱气处理后，在77K温度、相对压力介于0.01～0.995范围内进行煤样的低温氮吸/脱附实验，记录每个相对压力点下的吸附量，并依据BET模型或BJH模型计算对应比表面积S以及孔体积V。
进一步地，所述表面弛豫率的计算方法具体为：将横向弛豫时间简化为只有表面弛豫的等式：
Figure BDA0002038334110000031
其中T2表示实验测得的孔隙流体的弛豫时间，ρ2为煤样横向表面弛豫率，S为孔隙的表面积，V表示孔隙体积。
进一步地，根据
Figure BDA0002038334110000032
在等式两边均取孔径范围内以孔体积为权重的对数平均值，得到：
Figure BDA0002038334110000033
其中，T2LM表示弛豫时间T2分布的对数平均值，单位ms；T2i代表孔径范围内第i个T2值，单位ms；
Figure BDA0002038334110000034
代表T2值对应的NMR信号量的比例；
计算液氮吸附实验得到的区间比表面积S和孔体积V之比的对数平均值：
Figure BDA0002038334110000035
其中，S和V分别代表液氮吸附实验得到的比表面积和孔体积，单位分别为m2/g、cm3/g；
Figure BDA0002038334110000036
代表液氮吸附实验得到的区间比表面积S和孔体积V之比的对数平均值，单位×104cm-1；Sj和Vj分别代表液氮吸附实验得到的第j个孔径段孔隙的比表面积和孔体积，单位分别为m2/g、cm3/g；
根据上述两个公式对煤样依次进行计算，得到煤样横向表面弛豫率。
进一步地，所述目标孔径计算范围为10～100nm。
进一步地，运用煤样横向表面弛豫率可以对补偿后的弛豫时间截止值进行计算，其具体步骤为,将弛豫时间的表达式转化为包含孔径和孔隙几何因子的表达式：
Figure BDA0002038334110000041
其中，Fs代表孔隙几何因子，研究储层孔隙为柱状孔，取值为2；r代表孔隙半径。进一步地，还包括核磁信号补偿的验证，其具体步骤为，
根据弛豫时间T2的值与煤样孔径r之间的对应关系：
Figure BDA0002038334110000042
其中，弛豫时间截止值T2C为弛豫时间T2的截止值；
根据补偿后T2C值的计算结果对核磁共振表征孔径分布结构进行重新计算，并将该结算的结果与压汞实验测试、低温氮吸附测试对应的孔径分布结果进行对比。
本发明实施例具有如下优点：
本发明在实验中基于液氮吸附实验所得的数据来进行以孔体积为权重的对数平均值的计算，同时借助离心实验得到的核磁共振实验中T2值和r之间关系的不确定性，考虑到由于液氮吸附孔径分布分度值远大于核磁共振实验中的T2分布分度值，而且相对于全孔径范围内取对数平均值已经很大程度的缩小了误差，所以不会影响结果的准确性，通过这种补偿的方式来克服饱和水煤样进行离心操作时可能对不同煤体结构的煤样造成差异性破坏而产生核磁信号量损失的缺陷，从而提高孔隙表征结果的准确性。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。
本说明书所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本发明可实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本发明所能产生的功效及所能达成的目的下，均应仍落在本发明所揭示的技术内容得能涵盖的范围内。
图1为本发明实验过程流程示意图；
图2为本发明实验样品常规测试信息图；
图3为本发明核磁信号补偿前后T2C对比图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明提供了一种对核磁共振实验中煤样离心过程信号损失补偿方法，包括如下步骤：
步骤100、选取不同变质程度的煤样，对饱和水煤样进行核磁共振实验测试，对所述饱和水煤样进行离心实验后获得离心煤样，对离心煤样再次进行核磁共振实验测试，计算弛豫时间截止值并作出前后两次核磁共振的孔径分布曲线；
步骤200、对原始的煤样进行低温氮吸附实验测试，计算得到比表面积S和孔体积V；
步骤300、基于表面弛豫率的计算，对离心实验后的煤样在核磁共振实验中的信号损失量进行补偿。
煤样弛豫时间与煤样孔隙半径呈正相关关系，弛豫时间越长代表的煤岩孔隙半径越大，因此在现有技术中直接借用核磁共振实验得到的弛豫时间对应的核磁共振信号量对孔隙分布进行直接反应。这种方式虽然是一种无损测量方法，但是由于离心实验中离心力过大会对不同煤体结构发育煤样产生差异性破坏，煤样的损失从而导致核磁共振实验过程必然会损失部分核磁信号量，造成T2C值计算结果偏小，而以损失了部分信号量的T2C值作为直接的反应依据，其结果的可靠性和精度必然会有所缺陷，具体表现为不同实验在孔隙孔径分布上产生差异、峰值错位等，正是基于此，在本发明中构建了如何补偿煤样离心过程造成的核磁信号损失。在本实施方式中，为了更好的说明上述技术方案，将结合具体的实例来做进一步的说明和阐述。
本实施例选用的16组煤岩样品分别采自准噶尔、鄂尔多斯和沁水三个不同盆地的若干矿区，样品范围从低变质的褐煤到中变质的长焰煤、肥煤、焦煤和瘦煤以及高变质的无烟煤均有，并对煤样进行了镜质组反射率(Ro,max)、显微组分、工业分析、气测孔隙度和渗透率等系列实验，实验样品常规测试信息如图2所示。
在本实施方式中，对煤样进行以下实验，以确定相关的参数和关系：
(1)低温氮吸附实验
本实施例对16块煤岩样品进行了低温氮吸/脱附实验(N2GA)，在步骤200中，煤样低温氮吸附实验测试的具体步骤为：
步骤201、将煤样破碎并筛选60～80目的煤粉，将筛选后的煤样进行脱气处理后，在77K温度、相对压力介于0.01～0.995范围内进行煤样的低温氮吸/脱附实验，记录每个相对压力点下的吸附量，并依据BET模型或BJH模型计算对应比表面积S以及孔体积V。
N2GA实验对孔径测试范围在1.7～200nm之间，在液氮温度(77K)下，实验中的吸附量完全由相对压力P/P0决定，其中，P表示氮气实际压力，P0为77K温度下氮气的饱和蒸气压。当P/P0在0.05～0.35时，可通过BET模型对孔隙比表面积进行计算，在P/P0≥0.40时，通过BJH模型可以分析煤孔隙体积及孔径分布状况。
(2)核磁共振实验
核磁共振实验(NMR实验)的分析基础在于岩石中所含流体的自旋氢核在均匀分布的静磁场和射频场中的核磁共振弛豫行为。所述表面弛豫率的计算方法具体为，将横向弛豫时间简化为只有表面弛豫的等式：
Figure BDA0002038334110000071
式中：T2表示实验测得的孔隙流体的弛豫时间；ρ2为煤样横向表面弛豫率，常数；S为孔隙的表面积；V表示孔隙体积。所以根据公式(2-1)，煤样T2值与孔隙半径呈正比，弛豫时间越长代表的煤岩孔隙半径越大。为了确定T2值与孔径之间的确定对应关系，借助离心技术将孔隙中可动流体离心出去，根据Washburn方程：
Figure BDA0002038334110000072
式中，Pc—毛管压力，此处为离心力，MPa；σ—表面张力，N/m；θ—润湿接触角，(°)；rc—毛管半径，μm。当煤中孔隙流体为水时，表面张力为0.076N/m，接触角为60°，迭代后可得：
Figure BDA0002038334110000073
离心力与孔隙半径之间满足公式(2-3)。通过NMR实验测定离心后煤样核磁信号量，显示的信号量即为孔径小于rc时流体的信号量，近似于饱和水核磁共振数据中T2值小于T2C范围内孔隙流体信号量，所以通过rc值与T2C的对应关系可用来反映煤样整体孔径分布特征。
本实施例对16块煤样进行了低场核磁共振(NMR)实验。实验具体的步骤为：首先在新鲜煤样上钻取直径2.5cm、长度5.0cm的煤柱，将煤柱在105℃下烘干至恒重，然后置于真空泵中抽真空48h，在其内注入蒸馏水，并在8MPa的压力下对煤样饱和72h，使煤样达到完全饱和水状态，完成样品预处理。
在步骤100中，煤样在饱和水状态和离心实验后的核磁共振实验的具体步骤为：
将饱和水煤样在2.38MHz的Rec Core 2500低场核磁共振测量仪中进行实验，反演T2弛豫时间谱，实验设置参数为：回波间隔时间0.240ms，等待时间6s，回波个数8000，扫描次数64，实验温度25℃；
将离心后的煤样在相同的测试环境下再次进行NMR实验。
进一步地，所述饱和水煤样进行离心实验时需要先确定最佳离心力，其具体步骤为：将饱和水煤样在温度25℃，离心转速分别为3000r/min、4000r/min、5000r/min、6000r/min、7000r/min下进行离心，不同的离心转速对应离心力分别为0.69、0.92、1.15、1.38、1.61MPa，离心时间均为1h，通过不同离心转速作用后煤样的核磁共振信号量变化量来确定最佳离心力。最佳离心力主要是用于对核磁共振T2C值进行划分，用来准确表征孔径分布状态，同时也有利于样品测量数据的稳定，因为离心过程会导致构造发育的煤柱表面出现裂缝甚至破碎，造成核磁共振实验中的信号量存在误差，最佳离心力的确定可以确保因实验条件造成的系统误差的稳定。
基于上述实验过程，利用饱和水煤样以及不同离心力作用后煤样的核磁共振信号量，计算并比较煤样含水饱和度的变化来反映最佳离心力。
实验结果显示：随离心力增大，煤样含水饱和度下降程度逐渐降低，当离心力从1.38MPa增至1.61MPa时，煤样含水饱和度下降范围在0.51％～1.01％之间，平均0.71％，其含水饱和度下降范围基本可以忽略，故可以认定对于不同煤阶煤样，1.38MPa均可认定为使得饱和水煤样达到束缚水状态的最佳离心力。
使用核磁共振技术表征煤储层孔隙特征，需要通过T2截止值(T2C)来进行T2值和孔径的代换。在现有技术中，为了精确计算核磁共振实验的T2C值，除了离心法外，为了直接计算T2C建立了很多模型，如直接将T2分布曲线两个峰的谷底作为T2C值，但其运用具有很大的局限性，且缺少理论依据，又如以平均弛豫时间参数模型来计算T2C，通过孔隙分形方法定义T2C，二者结果准确性的衡量是通过与离心得到的T2C对比拟合来说明的，由于未从机理方面对T2值进行分析，所以拟合结果存在一定的偶然性。本发明中通过离心方式来计算T2C值，核磁共振实验测试的T2C具体计算步骤为：
在核磁共振实验测试中直接获得饱和水煤样和离心后煤样核磁共振T2谱图；
对离心前后的T2谱图分别做核磁共振信号累计幅度曲线，并从离心后的累计幅度曲线最大值处作X轴的平行线，直至与饱和水累计幅度曲线相交，并由交点处再作垂线与X轴相交，即所述垂线的横坐标即为T2C值。
但是由于离心过程煤样损失会导致离心煤样信号量偏小，而且本次实验中由于进行了多次离心，所以信号损失量会更大，而在现有技术条件下还不能够有效的对损失的信号量进行补偿。为了解决上述问题，在本发明中引入横向表面弛豫率的概念。
横向表面弛豫率与横向表面弛豫时间相对应，在核磁共振测井方面被广泛应用于弛豫时间与孔隙孔径分布之间的代换，以此来进行地层岩石物性的基础反映。依据式(2-1)弛豫时间可以转化为孔径的关系式：
Figure BDA0002038334110000091
其中，Fs代表孔隙几何形状因子(平行孔，Fs＝1；柱状孔，Fs＝2；球形孔，Fs＝3)，r表示T2值对应的孔径，在本实施方式中均假设孔隙为柱状孔，Fs＝2，利用上式即可通过ρ2值实现T2和r之间的转化，进而求取准确的T2C。
在现有技术中，虽然计算表面弛豫率的方法有很多，如将T2分布的峰值与液氮和压汞曲线的某一个峰值进行对应，直接确定T2值与r之间的关系，根据上式即可计算ρ2，但此方法在泥岩矿物含量及进汞饱和度的范围定义上没有定量分析，无法划分具体的范围，而且由于液氮及压汞实验测得的均为某一孔径范围内对应的孔体积，其孔径范围明显大于NMR实验T2值对应的孔径范围，所以直接进行峰值对应产生的误差会引起ρ2的计算误差；利用以
Figure BDA0002038334110000092
为权重计算T2的对数平均值的方法将核磁共振实验得到的T2值进行转化，进而与液氮吸附实验得到的S和V进行ρ2值计算，此方法虽然考虑到了每个T2值代表的孔径段占有的权重，但使用对数平均数即最具有代表性的T2值，将其与液氮吸附测量全孔径范围的S和V进行对应计算，仍存在明显错误。
基于煤储层孔径分布曲线，液氮吸附得到的S和V主要代表的孔径范围在10～100nm之间。因此：基于上述，本发明创造性的提出，在等式(2-1)两边均取孔径范围内以孔体积为权重的对数平均值，得到：
Figure BDA0002038334110000101
其中，T2LM表示T2分布的对数平均值，单位ms；T2i代表孔径范围内第i个T2值，单位ms；
Figure BDA0002038334110000102
代表T2值对应的NMR信号量的比例。在上述中，虽然在核磁共振实验中借助离心实验得到的T2值与r之间的关系不准确，但是由于液氮吸附孔径分布值远大于核磁共振实验中T2所代表的孔径分度值，而且相对于全孔径范围内取对数平均值已经很大程度的缩小了误差，所以基本不会影响结果的准确性，通过这种方式将克服现有技术方案中不确定性的因素，减少误差传递导致的计算精度问题，而这也是本实施方式中的主要创新点。
计算液氮吸附实验得到的区间比表面积S和孔体积V之比的对数平均值：
Figure BDA0002038334110000103
其中，S和V分别代表液氮吸附实验得到的比表面积和孔体积，单位分别为m2/g、cm3/g；
Figure BDA0002038334110000104
代表液氮吸附实验得到的区间比表面积S和孔体积V之比的对数平均值，单位×104cm-1；Sj和Vj分别代表液氮吸附实验得到的第j个孔径段孔隙的比表面积和孔体积，单位分别为m2/g、cm3/；
根据上述两个公式对煤样依次进行计算，得到煤样横向表面弛豫率，利用煤样横向表面弛豫率的计算结果，结合弛豫时间和孔径之间的关系式，取Fs＝2，从而获得核磁信号补偿后的T2截止值。
基于上述，为了进一步验证核磁信号补偿结果的准确性，如图3所示，依据补偿T2C值的计算结果，结合T2值与孔径之间的转化关系，对核磁共振实验表征的孔径分布结果进行重新计算，其与压汞实验测试、低温氮吸附测试的对比结果显示：由于补偿T2C值大于离心T2C值，而该值均对应于100nm孔径段，所以在孔径分布结果中，核磁共振实验孔径分布曲线会向左平移一段距离；孔径分布曲线中，MIP孔径分布曲线覆盖孔径范围较大，而结果显示，利用液氮吸附实验得到的S、V进行补偿后，其与核磁共振实验分布曲线对比结果更好，峰值基本处于同一孔径段内，所以可基本反映出运用本实施方式核磁信号量补偿方法的计算准确性。
该验证的具体步骤为：
第一步，根据T2值与煤样孔径r之间的对应关系：
Figure BDA0002038334110000111
其中，T2C为T2截止值；
第二步，依据补偿后T2C值对核磁共振表征孔径分布结构进行重新计算，并将该结果与压汞实验测试、低温氮吸附测试对应的孔径分布结果进行对比。
本发明以NMR实验中离心对煤样的损坏为起点，通过引入横向表面弛豫率的概念，并与N2GA实验中得出的S、V相结合，建立T2与S、V以
Figure BDA0002038334110000112
为权重的对数平均值的方程，进而对核磁信号量进行补偿，并对补偿后的T2C值进行计算。通过上述过程克服了现有计算方式中模型不适用以及计算方法存在误差等缺陷，从根本上创造性的解决了核磁信号补偿的缺陷。
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (8)

1.一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，包括如下步骤：
步骤100、选取不同变质程度的煤样，对饱和水煤样进行核磁共振实验测试，对所述饱和水煤样进行离心实验后获得离心煤样，对离心煤样再次进行核磁共振实验测试，计算弛豫时间截止值并作出前后两次核磁共振的孔径分布曲线；
步骤200、对原始的煤样进行低温氮吸附实验测试，计算得到比表面积S和孔体积V；
步骤300、基于表面弛豫率的计算，对离心实验后的煤样在核磁共振实验中的信号损失量进行补偿；
所述表面弛豫率的计算方法具体为：将横向弛豫时间简化为只有表面弛豫的等式：
Figure FDA0002584742580000011
其中T2表示实验测得的孔隙流体的弛豫时间，ρ2为煤样横向表面弛豫率，S为孔隙的表面积，V表示孔隙体积；
根据
Figure FDA0002584742580000012
在等式两边均取孔径范围内以孔体积为权重的对数平均值，得到：
Figure FDA0002584742580000013
其中，T2LM表示弛豫时间T2分布的对数平均值，单位ms；T2i代表孔径范围内第i个T2值，单位ms；
Figure FDA0002584742580000014
代表T2值对应的NMR信号量的比例；
计算液氮吸附实验得到的区间比表面积S和孔体积V之比的对数平均值：
Figure FDA0002584742580000015
其中，S和V分别代表液氮吸附实验得到的比表面积和孔体积，单位分别为m2/g、cm3/g；
Figure FDA0002584742580000021
代表液氮吸附实验得到的区间比表面积S和孔体积V之比的对数平均值，单位×104cm-1；Sj和Vj分别代表液氮吸附实验得到的第j个孔径段孔隙的比表面积和孔体积，单位分别为m2/g、cm3/g；
根据上述两个公式对煤样依次进行计算，得到煤样横向表面弛豫率。
2.根据权利要求1所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，在步骤100中，饱和水煤样和离心煤样的核磁共振实验具体步骤为：
将饱和水煤样在2.38MHz的Rec Core 2500低场核磁共振测量仪中进行实验，反演T2弛豫时间谱，实验设置参数为：回波间隔时间0.240ms，等待时间6s，回波个数8000，扫描次数64，实验温度25℃；
将离心后的煤样在相同的测试环境下再次进行核磁共振实验。
3.根据权利要求1所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，还包括所述饱和水煤样进行离心实验前确定最佳离心力，其具体步骤为：将饱和水煤样在恒定温度25℃，离心力分别为0.69、0.92、1.15、1.38、1.61MPa的条件下依次进行离心，每次离心时间均为1h，通过在不同离心转速作用下获取的煤样核磁共振信号量来确定最佳离心力。
4.根据权利要求3所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，核磁共振实验测试弛豫时间截止值的具体计算步骤为：
在核磁共振实验测试中直接获得饱和水煤样和离心后煤样核磁共振的弛豫时间谱图；
对离心前后的弛豫时间谱图分别做核磁共振信号累计幅度曲线，并从离心后的累计幅度曲线最大值处作X轴的平行线，直至与饱和水累计幅度曲线相交，并由交点处再作垂线与X轴相交，即所述垂线的横坐标即为弛豫时间截止值。
5.根据权利要求1所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，在步骤200中，煤样低温液氮吸附实验测试的具体步骤为：
将煤样破碎并筛选60～80目的煤粉，将筛选后的煤样进行脱气处理后，在77K温度、相对压力介于0.01～0.995范围内进行煤样的低温氮吸/脱附实验，记录每个相对压力点下的吸附量，并依据BET模型或BJH模型计算对应比表面积S以及孔体积V。
6.根据权利要求1所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，所述孔径计算范围为10～100nm。
7.根据权利要求1所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，运用煤样横向表面弛豫率可以对补偿后的弛豫时间截止值进行计算，其具体步骤为,将弛豫时间的表达式转化为包含孔径和孔隙几何因子的表达式：
Figure FDA0002584742580000031
其中，Fs代表孔隙几何因子，研究储层孔隙为柱状孔，取值为2；r代表孔隙半径。
8.根据权利要求7所述一种对核磁共振实验中煤样离心过程信号损失补偿方法，其特征在于，还包括核磁信号补偿的验证，其具体步骤为，
根据弛豫时间T2的值与煤样孔径r之间的对应关系：
Figure FDA0002584742580000032
其中，弛豫时间截止值T2C为弛豫时间T2的截止值；
根据补偿后T2C值的计算结果对核磁共振表征孔径分布结构进行重新计算，并将该计算的结果与压汞实验测试、低温氮吸附测试对应的孔径分布结果进行对比。
CN201910333621.4A 2019-04-24 2019-04-24 一种对核磁共振实验中煤样离心过程信号损失补偿方法 Active CN110108616B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910333621.4A CN110108616B (zh) 2019-04-24 2019-04-24 一种对核磁共振实验中煤样离心过程信号损失补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910333621.4A CN110108616B (zh) 2019-04-24 2019-04-24 一种对核磁共振实验中煤样离心过程信号损失补偿方法

Publications (2)

Publication Number Publication Date
CN110108616A CN110108616A (zh) 2019-08-09
CN110108616B true CN110108616B (zh) 2020-09-25

Family

ID=67486516

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910333621.4A Active CN110108616B (zh) 2019-04-24 2019-04-24 一种对核磁共振实验中煤样离心过程信号损失补偿方法

Country Status (1)

Country Link
CN (1) CN110108616B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111311575A (zh) * 2020-02-14 2020-06-19 中国地质大学（北京） 煤的核磁共振表面弛豫率确定方法以及装置
CN113567485A (zh) * 2020-04-29 2021-10-29 中国石油天然气股份有限公司 确定致密砂岩表面弛豫率的方法及装置
CN112067645B (zh) * 2020-07-31 2022-07-05 中国石油大学(北京) 核磁共振弛豫时间与孔喉半径的转换方法
CN113008927A (zh) * 2021-04-01 2021-06-22 中海石油(中国)有限公司 一种煤层核磁共振t2谱与孔隙分布的转换方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106442600A (zh) * 2016-11-23 2017-02-22 中国石油大学(华东) 页岩束缚水含量确定方法
CN107014728A (zh) * 2016-01-27 2017-08-04 中国石油化工股份有限公司 一种孔隙测量方法
CN107250828A (zh) * 2015-01-16 2017-10-13 因格瑞恩股份有限公司 用于改进的井底nmr特性描述的切屑分析
CN208155804U (zh) * 2018-05-22 2018-11-27 西南石油大学 基于核磁共振的致密储层岩电测量装置
CN109001243A (zh) * 2018-08-30 2018-12-14 中国地质大学（北京） 一种采用低场核磁共振评价煤的动态水锁效应的方法与装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107250828A (zh) * 2015-01-16 2017-10-13 因格瑞恩股份有限公司 用于改进的井底nmr特性描述的切屑分析
CN107014728A (zh) * 2016-01-27 2017-08-04 中国石油化工股份有限公司 一种孔隙测量方法
CN106442600A (zh) * 2016-11-23 2017-02-22 中国石油大学(华东) 页岩束缚水含量确定方法
CN208155804U (zh) * 2018-05-22 2018-11-27 西南石油大学 基于核磁共振的致密储层岩电测量装置
CN109001243A (zh) * 2018-08-30 2018-12-14 中国地质大学（北京） 一种采用低场核磁共振评价煤的动态水锁效应的方法与装置

Also Published As

Publication number Publication date
CN110108616A (zh) 2019-08-09

Similar Documents

Publication Publication Date Title
CN110108616B (zh) 一种对核磁共振实验中煤样离心过程信号损失补偿方法
CN104697915B (zh) 一种页岩微观孔隙大小及流体分布的分析方法
CN105866002A (zh) 一种精确的含油页岩核磁共振孔隙度测试方法
CN110296931B (zh) 一种致密砂岩油水相对渗透率信息的表征方法及系统
Xu et al. Effective porosity in lignite using kerosene with low-field nuclear magnetic resonance
CN113075102B (zh) 一种建立多孔介质自发渗吸量与时间关系数学模型的方法
US20240027379A1 (en) Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids
CN104142290A (zh) 一种岩样总体积变密度测定装置及方法
CN115078210B (zh) 页岩孔隙结构测试方法
Tian et al. Influence of pore water on the gas storage of organic-rich shale
Zheng et al. A multifractal-based method for determination NMR dual T2 cutoffs in coals
Wang et al. Fractal wave velocity evolution model of ultrasonic propagation and attenuation in water-bearing coal
Zhao et al. Investigation of pore structure characteristics of marine organic-rich shales using low-pressure N 2 adsorption experiments and fractal theory
CN115389387A (zh) 一种评价岩心伤害的实验方法
CN116539655B (zh) 一种基于核磁共振技术评价致密砂岩储层水敏性的方法
Averbakh et al. Acoustic spectroscopy of fluid saturation effects in carbonate rock
Zhao et al. Pore size distribution of high volatile bituminous coal of the southern Junggar Basin: a full-scale characterization applying multiple methods
CN108846218A (zh) 预测泥页岩储层不同类型孔隙的内比表面积的方法及系统
CN114965539A (zh) 页岩气水置换检测方法及系统
Dong et al. Quantitative characterization and characteristic analysis of pore structure of shale-gas reservoir in the Sichuan Basin, China
Wang et al. Experimental Study on Spontaneous Imbibition of Coal Samples of Different Ranks Based on the NMR Relaxation Spectrum
CN111562206A (zh) 一种非常规油气储层含油岩石孔径分布特征的测量方法
CN111426617A (zh) 一种地质条件下陆相湖盆致密储层孔隙度的定量预测方法
Zhang et al. Characterization of pore structure and its relationship with methane adsorption on medium-high volatile bituminous coal: an experimental study using nuclear magnetic resonance
CN111965329A (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