CN112326710A - 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 - Google Patents

一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 Download PDF

Info

Publication number
CN112326710A
CN112326710A CN202011242806.3A CN202011242806A CN112326710A CN 112326710 A CN112326710 A CN 112326710A CN 202011242806 A CN202011242806 A CN 202011242806A CN 112326710 A CN112326710 A CN 112326710A
Authority
CN
China
Prior art keywords
magnetic resonance
nuclear magnetic
spectrum
centrifugal
state
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.)
Granted
Application number
CN202011242806.3A
Other languages
English (en)
Other versions
CN112326710B (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 Petroleum East China
Original Assignee
China University of Petroleum East China
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 Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202011242806.3A priority Critical patent/CN112326710B/zh
Publication of CN112326710A publication Critical patent/CN112326710A/zh
Application granted granted Critical
Publication of CN112326710B publication Critical patent/CN112326710B/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
    • 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
    • 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)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Developing Agents For Electrophotography (AREA)
  • Water Treatment By Electricity Or Magnetism (AREA)

Abstract

一种泥页岩饱和水后各离心状态核磁共振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谱校正方法，其特征在于：
步骤1：对泥页岩样品进行烘干、称重，确定泥页岩储层样品的干样质量，通过核磁共振实验分别测量该样品在饱和水状态和不同离心转速状态的核磁共振T2谱响应信号，将饱和水状态和不同离心转速状态的核磁共振T2谱响应信号除以泥页岩样品的干样质量，分别获得单位质量该页岩样品在饱和水状态和不同离心转速状态核磁共振T2谱响应信号，页岩储层样品的干样质量单位为克，核磁共振T2谱响应信号为无量纲参数，离心转速的单位为转/分钟；
步骤2：依据单位质量该页岩样品在饱和水状态和不同离心转速状态核磁共振T2谱响应信号，确定单位质量该泥页岩样品饱和水状态与各离心转速状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号和核磁共振T2谱响应信号幅度差均为无量纲参数，离心转速的单位为转/分钟；
步骤3：按照单位质量泥页岩样品饱和水状态核磁共振T2谱响应信号与各离心状态核磁共振T2谱左侧响应信号幅度差符合正态分布特征，先确定饱和水状态与各离心状态核磁共振T2谱响应信号幅度差左侧第一极大值，然后对该极大值左侧各弛豫时间所对应的响应信号幅度差依据下列公式进行校正，
Figure 632676DEST_PATH_IMAGE001
式中，ΔFj, i’是第i个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，弛豫时间的单位为毫秒（ms），ΔFj, i是第i个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号幅度差为无量纲参数，j=1，2，…n0，是离心转速由低至高的离心状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，i=1，2，…，n1-1，是核磁共振T2谱弛豫时间编号，弛豫时间编号为无量纲参数，其中n1是饱和水状态与离心状态核磁共振T2谱响应信号幅度差左侧第一个极大值所对应的弛豫时间编号，n2是饱和水状态与离心状态核磁共振核磁共振T2谱响应信号幅度差左侧第一个极小值所对应的弛豫时间编号，离心转速的单位为转/分钟；
步骤4：由低离心转速向高离心转速，分别确定饱和水状态与离心转速相邻两个核磁共振T2响应信号幅度差中相对应的极值之比，如果极大值的个数大于等于2，则从T2谱左侧的第一个极大值开始，利用下列公式向右侧计算响应信号幅度差极大值与响应信号幅度差极小值之间不同弛豫时间所对应的响应信号幅度差校正系数，并利用响应信号幅度差校正系数对饱和水状态与非最低离心转速状态核磁共振T2谱响应信号幅度差进行校正，
Figure 708079DEST_PATH_IMAGE002
Figure 655438DEST_PATH_IMAGE003
式中，λjj, k是饱和水状态与第jj个离心状态核磁共振核磁共振T2谱响应信号幅度差的校正系数，该校正系数是无量纲参数，其中k是响应信号幅度差极大值与响应信号幅度差极小值所对应的弛豫时间段之间的弛豫时间编号，弛豫时间编号为无量纲参数，弛豫时间的单位为毫秒（ms），ΔFmax,jj,x是饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差中的极大值，核磁共振T2谱响应信号幅度差中的极大值是无量纲参数，其中x是该极大值的核磁共振弛豫时间编号，弛豫时间编号为无量纲参数，ΔFmax,jj-1,x是饱和水状态与第jj-1个离心状态核磁共振T2谱响应信号幅度差中弛豫时间编号x附近的极大值，核磁共振T2谱响应信号幅度差中的极大值是无量纲参数，ΔFmin,jj,y是饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差中的极小值，核磁共振T2谱响应信号幅度差中的极小值是无量纲参数，其中y是该极小值的核磁共振弛豫时间编号，弛豫时间编号为无量纲参数，ΔFmin,jj-1,y是饱和水状态与第jj-1个离心状态核磁共振T2谱响应信号幅度差中弛豫时间编号y附近的极小值，核磁共振T2谱响应信号幅度差中的极小值是无量纲参数，ΔFjj, k’是第k个弛豫时间所对应的饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，ΔFjj, k是第k个弛豫时间所对应的饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号幅度差为无量纲参数，jj=2，…n0，是离心转速由低至高且非最低离心转速状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，离心转速的单位为转/分钟；
步骤5：利用饱和水状态与各离心状态核磁共振T2谱响应信号幅度差校正结果，结合饱和水状态核磁共振T2谱，利用下列公式计算校正后泥页岩样品饱和水之后的各离心状态核磁共振T2谱，
Figure 613029DEST_PATH_IMAGE004
式中，Fj, z’是第j个离心状态核磁共振T2谱中第z个弛豫时间所对应校正后响应值，核磁共振T2谱校正后响应值为无量纲参数，Fj, z是第j个离心状态核磁共振T2谱中第z个弛豫时间的响应值，核磁共振T2谱响应值为无量纲参数，ΔFj, z’是第z个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，F0, z是饱和水状态核磁共振T2谱中第z个弛豫时间的响应值，核磁共振T2谱响应值为无量纲参数，j=1，2，…n0，是离心转速由低至高的离心状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，z=1，2，…，10000，是核磁共振T2谱弛豫时间编号，弛豫时间编号为无量纲参数，弛豫时间的单位为毫秒（ms），n1是饱和水状态与离心状态核磁共振T2谱响应信号幅度差左侧第一个极大值所对应的弛豫时间编号，弛豫时间编号为无量纲参数，离心转速的单位为转/分钟。
本发明的有益效果：本发明计算页岩气储层赋存吸附气含量的方法，实现了对各温度、压力条件下页岩气储层赋存吸附气含量的定量计算，而且该评价方法易于操作、费用低廉，所评价的各温度、压力条件下页岩气储层赋存吸附气含量是页岩气勘探和开发中所必需的重要参数。
附图说明
图1是本发明的流程图。
图2是本发明实施实例1中单位质量泥页岩样品在饱和水状态、2000转/分钟离心转速状态和8000转/分钟离心转速状态的核磁共振T2谱。
图3是本发明实施实例1中单位质量泥页岩样品在饱和水状态与2000转/分钟离心转速状态、8000转/分钟离心转速状态的核磁共振T2谱响应信号幅度差。
图4是本发明实施实例1中单位质量泥页岩样品在饱和水状态与2000转/分钟离心转速状态、8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差的第一个极大值左侧幅度差校正结果。
图5是本发明实施实例1中单位质量泥页岩样品在饱和水状态与8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差的第一个极大值右侧幅度差校正结果。
图6是本发明实施实例1中单位质量泥页岩样品在饱和水之后2000转/分钟离心转速状态、8000转/分钟离心转速状态核磁共振T2谱响应信号校正结果。
具体实施方式：
实施例1：如图1所述，一种计算页岩气储层赋存吸附气含量的方法，含有以下步骤：
步骤1：对泥页岩样品进行烘干、称重，确定该泥页岩储层样品的干样质量为25g。通过核磁共振实验分别测量该样品在饱和水状态、2000转/分钟离心转速状态和8000转/分钟离心转速状态的核磁共振T2谱响应信号。将饱和水状态和不同离心转速状态的核磁共振T2谱响应信号除以泥页岩样品的干样质量，分别获得单位质量该页岩样品在饱和水状态、2000转/分钟离心转速状态和8000转/分钟离心转速状态核磁共振T2谱响应信号，核磁共振T2谱响应信号为无量纲参数。单位质量该样品在饱和水状态、2000转/分钟离心转速状态和8000转/分钟离心转速状态的核磁共振T2谱响应信号结果见图2。
步骤2：依据单位质量该页岩样品在饱和水状态、2000转/分钟离心转速状态和8000转/分钟离心转速状态核磁共振T2谱响应信号，确定单位质量该泥页岩样品饱和水状态与2000转/分钟核磁共振T2谱响应信号幅度差、饱和水状态与8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号和核磁共振T2谱响应信号幅度差均为无量纲参数。单位质量该泥页岩样品饱和水状态与2000转/分钟、8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差结果见图3。
步骤3：按照饱和水状态核磁共振T2谱响应信号与离心状态核磁共振T2谱左侧响应信号幅度差符合正态分布特征，先确定饱和水状态与2000转/分钟离心转速状态和8000转/分钟离心状态核磁共振T2谱响应信号幅度差左侧第一极大值，然后对该极大值左侧各弛豫时间所对应的响应信号幅度差依据下列公式进行校正。对单位质量泥页岩样品饱和水状态与2000转/分钟、8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差左侧第一个极大值左侧各弛豫时间所对应的响应信号幅度差校正结果见图4。
Figure 428670DEST_PATH_IMAGE005
式中，ΔFj, i’是第i个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，ΔFj, i是第i个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号幅度差为无量纲参数，j=1，2，是离心转速由低至高的离心状态编号，离心状态编号为无量纲参数，i=1，2，…，n1-1，是核磁共振T2谱弛豫时间编号，弛豫时间编号为无量纲参数，其中n1是饱和水状态与离心状态核磁共振T2谱响应信号幅度差左侧第一个极大值所对应的弛豫时间编号，当j为1时n1值均为87，当j为2时n1值均为85，n2是饱和水状态与离心状态核磁共振核磁共振T2谱响应信号幅度差左侧第一个极小值所对应的弛豫时间编号，当j为1和2时n2值均为60。
步骤4：由低离心转速向高离心转速，分别确定饱和水状态与离心转速相邻两个核磁共振T2响应信号幅度差中相对应的极值之比，如果极大值的个数大于等于2，则从T2谱左侧的第一个极大值开始，利用下列公式向右侧计算响应信号幅度差极大值与响应信号幅度差极小值之间不同弛豫时间所对应的响应信号幅度差校正系数，并利用响应信号幅度差校正系数对饱和水状态与8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差进行校正。对饱和水状态与8000转/分钟离心转速状态核磁共振T2谱响应信号幅度差第一个极大值右侧响应信号幅度差的校正结果见图5。
Figure 453520DEST_PATH_IMAGE002
Figure 837359DEST_PATH_IMAGE003
式中，λjj, k是饱和水状态与第jj个离心状态核磁共振核磁共振T2谱响应信号幅度差的校正系数，该校正系数是无量纲参数。其中k是响应信号幅度差极大值与响应信号幅度差极小值所对应的弛豫时间段之间的弛豫时间编号，弛豫时间编号为无量纲参数，k的范围为60至85。ΔFmax,jj,x是饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差中的极大值，该极大值为2.1409。其中x是该极大值的核磁共振弛豫时间编号，弛豫时间编号为无量纲参数，x值为85。ΔFmax,jj-1,x是饱和水状态与第jj-1个离心状态核磁共振T2谱响应信号幅度差中弛豫时间编号x附近的极大值，该极大值为0.9809。ΔFmin,jj,y是饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差中的极小值，该极小值为-0.0669。其中y是该极小值的核磁共振弛豫时间编号，弛豫时间编号为无量纲参数，y值为103。ΔFmin,jj-1,y是饱和水状态与第jj-1个离心状态核磁共振T2谱响应信号幅度差中弛豫时间编号y附近的极小值，该极小值为0.3306。ΔFjj, k’是第k个弛豫时间所对应的饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，ΔFjj, k是第k个弛豫时间所对应的饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号幅度差为无量纲参数，jj=2，是离心转速由低至高且非最低离心转速状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数。
步骤5：利用饱和水状态与各离心状态核磁共振T2谱响应信号幅度差校正结果，结合饱和水状态核磁共振T2谱，利用下列公式计算校正后泥页岩样品饱和水之后的各离心状态核磁共振T2谱。
Figure 485509DEST_PATH_IMAGE004
式中，Fj, z’是第j个离心状态核磁共振T2谱中第z个弛豫时间所对应校正后响应值，核磁共振T2谱校正后响应值为无量纲参数，Fj, z是第j个离心状态核磁共振T2谱中第z个弛豫时间的响应值，核磁共振T2谱响应值为无量纲参数，ΔFj, z’是第z个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，F0, z是饱和水状态核磁共振T2谱中第z个弛豫时间的响应值，核磁共振T2谱响应值为无量纲参数，j=1，2，…n0，是离心状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，z=1，2，…，10000，是核磁共振T2谱弛豫时间编号，弛豫时间编号为无量纲参数，n1是饱和水状态与离心状态核磁共振T2谱响应信号幅度差左侧第一个极大值所对应的弛豫时间编号，弛豫时间编号为无量纲参数。

Claims (1)

1.一种泥页岩饱和水后各离心状态核磁共振T2谱校正方法，其特征在于：
步骤1：对泥页岩样品进行烘干、称重，确定泥页岩储层样品的干样质量，通过核磁共振实验分别测量该样品在饱和水状态和不同离心转速状态的核磁共振T2谱响应信号，将饱和水状态和不同离心转速状态的核磁共振T2谱响应信号除以泥页岩样品的干样质量，分别获得单位质量该页岩样品在饱和水状态和不同离心转速状态核磁共振T2谱响应信号，页岩储层样品的干样质量单位为克，核磁共振T2谱响应信号为无量纲参数，离心转速的单位为转/分钟；
步骤2：依据单位质量该页岩样品在饱和水状态和不同离心转速状态核磁共振T2谱响应信号，确定单位质量该泥页岩样品饱和水状态与各离心转速状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号和核磁共振T2谱响应信号幅度差均为无量纲参数，离心转速的单位为转/分钟；
步骤3：按照单位质量泥页岩样品饱和水状态核磁共振T2谱响应信号与各离心状态核磁共振T2谱左侧响应信号幅度差符合正态分布特征，先确定饱和水状态与各离心状态核磁共振T2谱响应信号幅度差左侧第一极大值，然后对该极大值左侧各弛豫时间所对应的响应信号幅度差依据下列公式进行校正，
Figure 112079DEST_PATH_IMAGE001
式中，ΔFj, i’是第i个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，弛豫时间的单位为毫秒（ms），ΔFj, i是第i个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号幅度差为无量纲参数，j=1，2，…n0，是离心转速由低至高的离心状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，i=1，2，…，n1-1，是核磁共振T2谱弛豫时间编号，弛豫时间编号为无量纲参数，其中n1是饱和水状态与离心状态核磁共振T2谱响应信号幅度差左侧第一个极大值所对应的弛豫时间编号，n2是饱和水状态与离心状态核磁共振核磁共振T2谱响应信号幅度差左侧第一个极小值所对应的弛豫时间编号，离心转速的单位为转/分钟；
步骤4：由低离心转速向高离心转速，分别确定饱和水状态与离心转速相邻两个核磁共振T2响应信号幅度差中相对应的极值之比，如果极大值的个数大于等于2，则从T2谱左侧的第一个极大值开始，利用下列公式向右侧计算响应信号幅度差极大值与响应信号幅度差极小值之间不同弛豫时间所对应的响应信号幅度差校正系数，并利用响应信号幅度差校正系数对饱和水状态与非最低离心转速状态核磁共振T2谱响应信号幅度差进行校正，
Figure 209611DEST_PATH_IMAGE002
Figure 280466DEST_PATH_IMAGE003
式中，λjj, k是饱和水状态与第jj个离心状态核磁共振核磁共振T2谱响应信号幅度差的校正系数，该校正系数是无量纲参数，其中k是响应信号幅度差极大值与响应信号幅度差极小值所对应的弛豫时间段之间的弛豫时间编号，弛豫时间编号为无量纲参数，弛豫时间的单位为毫秒（ms），ΔFmax,jj,x是饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差中的极大值，核磁共振T2谱响应信号幅度差中的极大值是无量纲参数，其中x是该极大值的核磁共振弛豫时间编号，弛豫时间编号为无量纲参数，ΔFmax,jj-1,x是饱和水状态与第jj-1个离心状态核磁共振T2谱响应信号幅度差中弛豫时间编号x附近的极大值，核磁共振T2谱响应信号幅度差中的极大值是无量纲参数，ΔFmin,jj,y是饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差中的极小值，核磁共振T2谱响应信号幅度差中的极小值是无量纲参数，其中y是该极小值的核磁共振弛豫时间编号，弛豫时间编号为无量纲参数，ΔFmin,jj-1,y是饱和水状态与第jj-1个离心状态核磁共振T2谱响应信号幅度差中弛豫时间编号y附近的极小值，核磁共振T2谱响应信号幅度差中的极小值是无量纲参数，ΔFjj, k’是第k个弛豫时间所对应的饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，ΔFjj, k是第k个弛豫时间所对应的饱和水状态与第jj个离心状态核磁共振T2谱响应信号幅度差，核磁共振T2谱响应信号幅度差为无量纲参数，jj=2，…n0，是离心转速由低至高且非最低离心转速状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，离心转速的单位为转/分钟；
步骤5：利用饱和水状态与各离心状态核磁共振T2谱响应信号幅度差校正结果，结合饱和水状态核磁共振T2谱，利用下列公式计算校正后泥页岩样品饱和水之后的各离心状态核磁共振T2谱，
Figure 772627DEST_PATH_IMAGE005
式中，Fj, z’是第j个离心状态核磁共振T2谱中第z个弛豫时间所对应校正后响应值，核磁共振T2谱校正后响应值为无量纲参数，Fj, z是第j个离心状态核磁共振T2谱中第z个弛豫时间的响应值，核磁共振T2谱响应值为无量纲参数，ΔFj, z’是第z个弛豫时间所对应的饱和水状态与第j个离心状态核磁共振T2谱响应信号幅度差校正结果，核磁共振T2谱响应信号幅度差为无量纲参数，F0, z是饱和水状态核磁共振T2谱中第z个弛豫时间的响应值，核磁共振T2谱响应值为无量纲参数，j=1，2，…n0，是离心转速由低至高的离心状态编号，离心状态编号为无量纲参数，其中n0是离心状态的数量，离心状态的数量为无量纲参数，z=1，2，…，10000，是核磁共振T2谱弛豫时间编号，弛豫时间编号为无量纲参数，弛豫时间的单位为毫秒（ms），n1是饱和水状态与离心状态核磁共振T2谱响应信号幅度差左侧第一个极大值所对应的弛豫时间编号，弛豫时间编号为无量纲参数，离心转速的单位为转/分钟。
CN202011242806.3A 2020-11-10 2020-11-10 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 Active CN112326710B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011242806.3A CN112326710B (zh) 2020-11-10 2020-11-10 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011242806.3A CN112326710B (zh) 2020-11-10 2020-11-10 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法

Publications (2)

Publication Number Publication Date
CN112326710A true CN112326710A (zh) 2021-02-05
CN112326710B CN112326710B (zh) 2022-06-03

Family

ID=74316556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011242806.3A Active CN112326710B (zh) 2020-11-10 2020-11-10 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法

Country Status (1)

Country Link
CN (1) CN112326710B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130113479A1 (en) * 2011-11-04 2013-05-09 Baker Hughes Incorporated Method for correcting nmr and nuclear logs in formate mud filtrate invaded formatons
CN104634804A (zh) * 2013-11-08 2015-05-20 中国石油天然气股份有限公司 一种利用核磁共振t2谱确定储层相对渗透率的方法
CN105464654A (zh) * 2015-12-15 2016-04-06 中国石油天然气股份有限公司 确定致密砂岩饱和度指数的方法及装置
CN106249306A (zh) * 2016-10-12 2016-12-21 贵州大学 基于核磁共振的页岩孔隙结构检测方法
CN106351652A (zh) * 2016-11-25 2017-01-25 中国地质大学（北京） 一种含烃储集层核磁共振测井t2谱形态校正方法
US20170139027A1 (en) * 2014-06-19 2017-05-18 Total Sa Nmr method for determining non-oil volume of a rock sample
CN107843611A (zh) * 2016-09-20 2018-03-27 中国石油化工股份有限公司 低渗透砂岩储层可动流体饱和度核磁共振参数表征新方法
CN108458960A (zh) * 2018-03-27 2018-08-28 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN108956678A (zh) * 2018-06-11 2018-12-07 西南石油大学 一种基于核磁共振测井的t2谱敏感参数提取方法
CN109725016A (zh) * 2018-11-29 2019-05-07 中国石油天然气集团有限公司 一种用于含重油、沥青质岩心的核磁共振实验室测量方法
US20190257733A1 (en) * 2018-06-09 2019-08-22 China University Of Petroleum (East China) Isotope nuclear magnetic method for analyzing ineffective water absorption of rock pores
CN110306977A (zh) * 2018-03-27 2019-10-08 中国石油化工股份有限公司 核磁共振t2谱含烃校正方法及系统

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130113479A1 (en) * 2011-11-04 2013-05-09 Baker Hughes Incorporated Method for correcting nmr and nuclear logs in formate mud filtrate invaded formatons
CN104634804A (zh) * 2013-11-08 2015-05-20 中国石油天然气股份有限公司 一种利用核磁共振t2谱确定储层相对渗透率的方法
US20170139027A1 (en) * 2014-06-19 2017-05-18 Total Sa Nmr method for determining non-oil volume of a rock sample
CN105464654A (zh) * 2015-12-15 2016-04-06 中国石油天然气股份有限公司 确定致密砂岩饱和度指数的方法及装置
CN107843611A (zh) * 2016-09-20 2018-03-27 中国石油化工股份有限公司 低渗透砂岩储层可动流体饱和度核磁共振参数表征新方法
CN106249306A (zh) * 2016-10-12 2016-12-21 贵州大学 基于核磁共振的页岩孔隙结构检测方法
CN106351652A (zh) * 2016-11-25 2017-01-25 中国地质大学（北京） 一种含烃储集层核磁共振测井t2谱形态校正方法
CN108458960A (zh) * 2018-03-27 2018-08-28 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN110306977A (zh) * 2018-03-27 2019-10-08 中国石油化工股份有限公司 核磁共振t2谱含烃校正方法及系统
US20190257733A1 (en) * 2018-06-09 2019-08-22 China University Of Petroleum (East China) Isotope nuclear magnetic method for analyzing ineffective water absorption of rock pores
CN108956678A (zh) * 2018-06-11 2018-12-07 西南石油大学 一种基于核磁共振测井的t2谱敏感参数提取方法
CN109725016A (zh) * 2018-11-29 2019-05-07 中国石油天然气集团有限公司 一种用于含重油、沥青质岩心的核磁共振实验室测量方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
FANGWEN CHEN 等: "Evaluation of the density and thickness of adsorbed methane in differently sized pores contributed by various components in a shale gas reservoir: A case study of the Longmaxi Shale in Southeast Chongqing, China", 《CHEMICAL ENGINEERING JOURNAL》 *
FANGWEN CHEN 等: "Evaluation of the density and thickness of adsorbed methane in differently sized pores contributed by various components in a shale gas reservoir: A case study of the Longmaxi Shale in Southeast Chongqing, China", 《CHEMICAL ENGINEERING JOURNAL》, no. 367, 1 July 2019 (2019-07-01) *
JIJUN LI等: "A comparison of experimental methods for describing shale pore features — A case study in the Bohai Bay Basin of eastern China", 《INTERNATIONAL JOURNAL OF COAL GEOLOGY》 *
JIJUN LI等: "A comparison of experimental methods for describing shale pore features — A case study in the Bohai Bay Basin of eastern China", 《INTERNATIONAL JOURNAL OF COAL GEOLOGY》, no. 152, 1 December 2015 (2015-12-01) *
宁传祥等: "泥页岩核磁共振T2谱换算孔隙半径方法", 《科学技术与工程》 *
宁传祥等: "泥页岩核磁共振T2谱换算孔隙半径方法", 《科学技术与工程》, vol. 16, no. 27, 30 September 2016 (2016-09-30) *
林君等: "基于不等式约束的磁共振信号T2谱多指数分解法及算法改进", 《吉林大学学报（地球科学版）》 *
林君等: "基于不等式约束的磁共振信号T2谱多指数分解法及算法改进", 《吉林大学学报（地球科学版）》, vol. 43, no. 6, 30 November 2013 (2013-11-30) *
陈方文等: "渤海湾盆地冀中坳陷饶阳凹陷沙一下亚段页岩油可动量评价", 《石油与天然气地质》 *
陈方文等: "渤海湾盆地冀中坳陷饶阳凹陷沙一下亚段页岩油可动量评价", 《石油与天然气地质》, vol. 40, no. 3, 30 June 2019 (2019-06-30) *

Also Published As

Publication number Publication date
CN112326710B (zh) 2022-06-03

Similar Documents

Publication Publication Date Title
CN104697915B (zh) 一种页岩微观孔隙大小及流体分布的分析方法
CN105866002B (zh) 一种精确的含油页岩核磁共振孔隙度测试方法
Li et al. Effect of pore structure on shale oil accumulation in the lower third member of the Shahejie formation, Zhanhua Sag, eastern China: Evidence from gas adsorption and nuclear magnetic resonance
CN109001243B (zh) 一种采用低场核磁共振评价煤的动态水锁效应的方法与装置
CN107991335B (zh) 致密砂岩水锁损害评价试验方法
CN109725016B (zh) 一种用于含重油、沥青质岩心的核磁共振实验室测量方法
CN109884109B (zh) 一种页岩孔隙水吸附量、游离量及其分布的评价方法
CN110296931B (zh) 一种致密砂岩油水相对渗透率信息的表征方法及系统
CA2524993C (en) Methods and apparatus for measuring capillary pressure in a sample
CN110108616B (zh) 一种对核磁共振实验中煤样离心过程信号损失补偿方法
CN107843611B (zh) 低渗透砂岩储层可动流体饱和度核磁共振参数表征新方法
CN111537543B (zh) 低场核磁共振测定页岩黏土与脆性矿物相对含量的方法
CN115078210B (zh) 页岩孔隙结构测试方法
CN104142290A (zh) 一种岩样总体积变密度测定装置及方法
CN113075102B (zh) 一种建立多孔介质自发渗吸量与时间关系数学模型的方法
CN112378943A (zh) 一种页岩含油饱和度评价模型、评价方法、应用
CN109781765A (zh) 一种计算致密储层束缚水液膜厚度的新方法
CN106483057A (zh) 一种定量评价超深储层可动流体的方法及其应用
CN112326710B (zh) 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法
CN111441765A (zh) 一种含裂缝致密油藏空气重力驱潜力评价实验方法及装置
CN111693361A (zh) 一种基于全流触探的土体强度测量装置
US20180143148A1 (en) A Method For Analysis Of Hydrocarbons In Water
EP0586001A1 (en) Process for generating, by way of experimental tests, the capillary pressure curve in a porous medium
CN115032222A (zh) 致密岩石核磁共振t2截止值的拟合计算方法
CN109991675A (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