CN115824939B - 一种基于核磁共振及三维切片的溶蚀改造量化研究方法 - Google Patents

一种基于核磁共振及三维切片的溶蚀改造量化研究方法 Download PDF

Info

Publication number
CN115824939B
CN115824939B CN202211620990.XA CN202211620990A CN115824939B CN 115824939 B CN115824939 B CN 115824939B CN 202211620990 A CN202211620990 A CN 202211620990A CN 115824939 B CN115824939 B CN 115824939B
Authority
CN
China
Prior art keywords
sample
plunger sample
plunger
corrosion
magnetic resonance
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
CN202211620990.XA
Other languages
English (en)
Other versions
CN115824939A (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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202211620990.XA priority Critical patent/CN115824939B/zh
Publication of CN115824939A publication Critical patent/CN115824939A/zh
Application granted granted Critical
Publication of CN115824939B publication Critical patent/CN115824939B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

本发明公开了一种基于核磁共振及三维切片的溶蚀改造量化研究方法,包括如下步骤:制备白云岩的柱塞样;基于模拟条件对所述柱塞样进行溶蚀模拟实验,收集反应后的生产液;溶蚀模拟实验后再次测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;对实验前后的柱塞样分别进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样溶蚀模拟实验后的孔隙结构特征;本发明中通过构建模拟条件,并结合有损切片以及实验的特点,针对性的构建镜像对比样品,从而在有损实验前后均可对样品的孔隙特征进行三维比对,在进行三维比对时通过优化切片的方式从而形成真实的三维结构,便于真实观察和测定样品的物性变化。

Description

一种基于核磁共振及三维切片的溶蚀改造量化研究方法
技术领域
本发明涉及石油地质勘探技术领域,具体涉及一种基于核磁共振及三维切片的溶蚀改造量化研究方法。
背景技术
白云岩储层因其在世界范围内可观的油气储量而日益受到关注,且近年来,我国深层-超深层碳酸盐岩储层出现重大的油气勘探突破。
深层-超深层优质规模性储层的形成、分布规律成为近期研究的关注焦点。前人研究表明:
埋藏阶段的有机酸、TSR、热液等作用可以形成溶蚀孔洞,然而,由于白云岩储层复杂的成岩演化和强烈的非均质性,使得储层的开发和预测变得困难。深层-超深层溶蚀作用对三维孔隙结构的影响需要进一步研究。
核磁共振T2图谱测量作为一种样品孔隙结构分析方法,可以在不破坏样品的情况下分析样品的孔径分布、孔隙连通性等特征。分形维数(D)作为描述物体自相似性和粗糙度的重要参数,在储层孔隙结构的非均质性研究中得到广泛应用。
由于在溶蚀实验以及切片过程中均会对样品产生损伤,一旦进行前述任一操作后,同一个样品将无法进行后续的实验或者切片。因此,在现有技术中,样品无法进行实验前后的比对,而单纯的核磁共振数据并不足以清晰的反应实验前后物性的变化。基于此,提出了一种利用白云岩溶蚀模拟实验结合基于核磁共振的分形维数分析的方法,通过正演模拟埋藏环境下的温压、流体条件,结合孔隙结构的非均质性表征参数分析,从而为白云岩储层溶蚀孔隙演化规律研究提供依据。
发明内容
本发明的目的在于提供一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其通过溶蚀模拟实验,结合岩心核磁共振T2谱测量、分形维数分析等手段,定量研究白云岩微纳米孔溶蚀改造效果,以解决现有技术中由于样品有损无法进行前后比对的技术问题。
为解决上述技术问题,本发明具体提供下述技术方案:
一种基于核磁共振及三维切片的溶蚀改造量化研究方法,包括如下步骤:
制备白云岩的柱塞样,并对柱塞样进行预处理后称重,测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;
构建模拟条件,基于所述模拟条件对所述柱塞样进行溶蚀模拟实验,收集反应后的生产液;
溶蚀模拟实验后对所述柱塞样进行预处理后称重,并再次测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;
对所述柱塞样进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样溶蚀模拟实验后的孔隙结构特征。
进一步地,所述柱塞样至少包括两组平行样品,所述平行样品的制备具体包括如下步骤:
在岩样相邻位置钻取至少两个呈圆柱状的柱塞样,沿其中一个所述柱塞样的中轴线切分为柱塞样a和柱塞样b,其中所述柱塞样a和所述柱塞样b为一组平行样品,所述柱塞样a和所述柱塞样b的镜像相似度不低于90%。
进一步地,在进行溶蚀模拟实验前,对所述柱塞样a进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样a在溶蚀模拟实验前的孔隙结构特征;
在进行溶蚀模拟实验后,对所述柱塞样b进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样b在溶蚀模拟实验后的孔隙结构特征。
进一步地,在进行溶蚀模拟实验前后对所述柱塞样进行核磁共振测定分析所述柱塞样孔隙结构特征的具体方法为:
对所述柱塞样进行抽真空处理,并且真空条件下对所述柱塞样充注饱和蒸馏水直至压力大于常压,并维持该压力指定时间后测定所述柱塞样的T2谱;
对所述柱塞样进行多离心参数的离心处理后再次测定所述柱塞样的T2谱;
其中:
所述离心参数包括离心力和离心时间,在离心处理时通过耦合拟合的方式确定所述柱塞样的最佳离心参数;
基于分形维数分析所述柱塞样在核磁共振中三维孔隙结构的复杂程度,所述分形维数记为D:
logVp=(3-D)log(T2)+(D-3)log(T2max);
其中:
T2max是T2的最大值,Vp是由等于或小于T2所对应孔径所占孔隙体积占总孔隙体积的体积分数。
进一步地,计算所述柱塞样的T2截止值,以T2截止值为分界计算分段的分形维数。
进一步地,构建模拟条件包括模拟地层水类型、模拟实验温压条件和模拟实验流速。
进一步地,模拟地层水类型具体包括如下步骤:
采用苏林分类法计算原生水型的特性系数,获得研究区的地层水类型,对应配制不同酸性介质混合流体;
模拟实验温压条件具体包括如下步骤:
对研究区的地层埋藏史、温压场进行分析,获得研究区实际地质条件的温压梯度,确定溶蚀模拟实验所需的温压条件;
其中:
温度通过地温梯度计算法推算所取样品对应深度的温度;
压力采用静水压力公式计算法推算所取样品对应深度的压力。
进一步地,当研究区地层中有碳酸盐岩矿物包裹体数据时,以碳酸盐岩矿物包裹体为对象模拟包裹体成分数据,对应配制不同酸性介质混合流体。
进一步地,模拟实验流速具体包括如下步骤:
设定溶蚀模拟实验所需的地层水类型以及温压条件不变,按照溶蚀模拟实验的条件改变实验流速进行多次验证实验,在验证实验过程中以指定的时间间隔收集水样,并对水样进行离子成分分析,当水样的离子成分不再发生变化时,判断反应体系达到平衡状态,对应的流速即为溶蚀模拟实验的最佳反应流速。
本发明与现有技术相比具有如下有益效果:
本发明中通过构建模拟条件,并结合有损切片以及实验的特点,针对性的构建镜像对比样品,从而在有损实验前后均可对样品的孔隙特征进行三维比对,在进行三维比对时通过优化切片的方式从而形成真实的三维结构,便于真实观察和测定样品的物性变化。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种基于核磁共振及三维切片的溶蚀改造量化研究方法,包括如下步骤:
制备白云岩的柱塞样,并对柱塞样进行预处理后称重,测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;
构建模拟条件,基于所述模拟条件对所述柱塞样进行溶蚀模拟实验,收集反应后的生产液;
溶蚀模拟实验后对所述柱塞样进行预处理后称重,并再次测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;
对所述柱塞样进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样溶蚀模拟实验后的孔隙结构特征。
所述柱塞样至少包括两组平行样品,所述平行样品的制备具体包括如下步骤:
在岩样相邻位置钻取至少两个呈圆柱状的柱塞样(如每个柱塞样直径为2.5cm,高为6cm),沿其中一个所述柱塞样的中轴线切分为柱塞样a和柱塞样b,其中所述柱塞样a和所述柱塞样b为一组平行样品。
需要说明的是,在本发明中根据溶蚀实验的需求可以确定多组平行样品,其中每个平行样品中均包括柱塞样a和柱塞样b,其中柱塞样a和柱塞样b均为直径为2.5cm、高为3cm的圆柱体。
每组平行样品对应两组实验条件完全相同的重复验证实验,这是因为在切薄片时需要考虑到对样品的损坏,即在切薄片后无法再使用该样品进行实验。为了保证平行样品在溶蚀模拟实验前后可以进行比对,所述柱塞样a和所述柱塞样b的镜像相似度不低于90%。
因此,在本发明中:
由于切片(特制制备铸体薄片)和溶蚀模拟实验均是对样品产生损伤的,因此,不论是在进行切片还是溶蚀模拟实验时,一旦进行前述任一操作将无法再次使用该样品进行下一步操作。然而,在实际研究对比的过程中,真实的前后对比往往是最直观且最有说服力的证据,因此,在本实施方式中基于有损制样或实验的具体缺陷构建平行样品来实现制样、实验的同步进行和比对,从而克服前述缺陷。
在本实施方式中,采用相似度高的样品,设置反应条件完全相同的两组平行实验,防止产生偶然性误差。
另外,由于柱塞样a和柱塞样b是镜面对称(至少包括一个面是镜面对称的,因为共面)的,由此可以保证柱塞样b在溶蚀模拟实验前后可重复测定孔隙度和渗透率,而柱塞样a镜面剖开侧的岩片可以用来制薄片,柱塞样b镜面剖开侧的岩片可以用来进行溶蚀实验,并在实验后制作薄片以此来进行对比。由于两个岩片位置镜面对应,可以对比溶蚀反应前后的矿物。
在本发明中,为了更加全面的观察,因此在对样品进行切片制备铸体薄片时,采用三维切片的方式进行取样,从而获得不同平面上的孔隙结构,由此可以根据需要对其孔隙结构进行三维模拟构建来完成三维重建,从而在三维上来比对溶蚀过程中对于孔隙的作用。
需要强调的是,在三维重建的过程中需要配合核磁共振的数据,关于核磁共振的特征在后续中进一步说明。
基于前述,由于切片是有损样品的,因此在本实施方式中,具体的切片方式为:
在进行溶蚀模拟实验前,对所述柱塞样a进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样a在溶蚀模拟实验前的孔隙结构特征;
在进行溶蚀模拟实验后,对所述柱塞样b进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样b在溶蚀模拟实验后的孔隙结构特征。
特别的是,在本实施方式中,对于三维方向上的切片,每个维度上的切片数量不限于一个,切片的数量根据样品矿物或/和孔隙(吼道)等的大小决定。由于在每个方向进行切片后,样品的长度会有所缺损,因此在切片时需要对其位置进行标定,从而便于后续的三维重建。
对样品的预处理包括清洗和烘干,在烘干后进行称重。清洗使用超声清洗器配合去离子水清洗样品,将清洗后的样品置于烘箱内,100-110℃下持续烘干24h后称重,可以去除样品表面杂质离子对实验反应进程的影响,减小称重误差。
构建模拟条件,基于所述模拟条件进行开放式的溶蚀模拟实验,并收集反应后的生成液。
构建模拟条件包括模拟地层水类型、模拟实验温压条件和模拟实验流速。
模拟地层水类型具体包括如下步骤:
采用苏林分类法计算原生水型的特性系数,获得研究区的地层水类型,对应配制不同酸性介质混合流体。
其中,
当有研究区地层有碳酸盐岩矿物包裹体数据时,优先模拟包裹体成分数据,配制对应离子成分的不同酸性介质混合流体。相较于研究目标区域内地层水数据,包裹体数据能更好地反应成矿及成岩过程中的流体特点。
溶剂采用去离子水,盐类物质选择试剂纯的NaCl、CaCl2、MgCl2或Na2SO4,酸性介质选择CO2水溶液或油田水中常见的乙酸。
模拟实验温压条件具体包括如下步骤:
对研究区的地层埋藏史、温压场进行分析,获得研究区实际地质条件的温压梯度,确定溶蚀模拟实验所需的温压条件;
其中:
温度主要通过地温梯度计算法确定,即基于地热基本方法及理论,通过标准或近标准的地温梯度推算所取样品对应深度的温度;
压力主要采用静水压力公式计算法,确定所取样品对应深度的压力。
溶蚀模拟实验反应温度与压力的调整具体包括三种方式:
(1)溶蚀模拟实验温度维持不变,压力按压力梯度变化;
(2)溶蚀模拟实验压力维持不变,温度按地温梯度变化;
(3)溶蚀模拟实验温度压力按温压梯度同时变化。
在上述参数调整的过程中,通过控制变量法(1)(2)调整实验温压条件,实现温度/压力单因素条件控制下的溶蚀增孔效应分析,排除干扰;通过调整体系温压条件参数同时按地层深度变化(3),实现温度压力同时控制的溶蚀增孔效应分析,更接近真实埋藏环境。
埋藏环境下的流体-岩石作用具有孔隙空间狭小、岩石固体比面积大的特点,反应进程由化学热力学作用主导。通过设置流速验证实验,获得溶蚀模拟实验的最佳流速,可以排除流速变化对反应的干扰。
因此,一方面需要通过验证实验获得流速数据,另外在该过程中需要根据溶蚀的进度对流速数据进行适应性的调整,从而减少固定流速对反应的干扰。
在本实施方式中,模拟实验流速具体包括如下步骤:
设定溶蚀模拟实验所需的地层水类型以及温压条件不变,按照溶蚀模拟实验的条件改变实验流速进行多次验证实验,在验证实验过程中以指定的时间间隔收集水样,并对水样进行离子成分分析,当水样的离子成分不再发生变化时,判断反应体系达到平衡状态,对应的流速即为溶蚀模拟实验的最佳反应流速。
本发明中采用开放体系作为溶蚀模拟实验反应体系,在开放体系中,大量地层侵蚀性不饱和流体与岩石表面接触,导致水岩反应体系向着溶解的方向进行,形成一定规模溶孔。而封闭体系水岩比较低,流体迅速饱和,导致水岩反应体系达到平衡,因此埋藏环境主要以孔隙的保存为主,对储层物性的改变较小。含膏白云岩储层溶蚀具有事件性的特征,埋藏环境下开放体系和封闭体系交替出现,发生大规模溶蚀时的体系近似开放。为了定量表征含膏白云岩储层的溶蚀增孔效应,需要确保反应已经达到动态平衡,获得一定条件下的储层饱和溶蚀量,因此选择开放体系作为溶蚀模拟实验的体系。
在开放体系中,模拟条件的动态调整时最为重要的,由于在动态调整的过程中没有任何的依据可以参考,因此,在调整的过程中势必会造成一定的模拟误差,在本实施例中,依据具体的模拟条件,通过参数响应以及设定约束条件的方式来优化参数的设计值,从而使得整个模拟过程可以更符合要求。
在对本实施方式中的技术方案进行具体说明之前,首先需要强调的是开放式体系和闭合式体系之间的区别,以及开放式体系模拟条件设定的困难。
在闭合式体系中,由于模拟条件比较固定,因此起参数的设定是有据可循的,而在开放式体系中,由于模拟条件多而且是动态的,无法提前设计相关的数据,需要根据模拟过程进行动态的变化,
而且多个参数响应之间存在高度的相关性或具有不同的方差结构时可能难以获得合理的参数设计值,并且在模拟参数时可以会由于参数响应之间的相互干涉、冲突时难以通过相应的处理方式获得具体的参数设计。
具体溶蚀模拟实验步骤为:
实验开始前进行0.5h的常温常压预反应,增温增压至反应所需条件时开始正式实验,反应时长为12h;
实验过程中共收集两份水样各10mL,当生成液体积与柱塞样体积相等时收集第一份水样,反应进行至12h时收集第二份水样;
采用ICP-OES对实验过程中所取水样进行离子成分、浓度分析,与所配制的溶蚀流体进行比较,定量评价溶蚀能力。
溶蚀模拟实验后对样品进行预处理后称重,并对矿物表面形态进行观察,测定样品的孔隙度和渗透性,测试生成液的化学成分。
此处对样品进行预处理同样包括对样品的清洗和烘干,即将清洗后的样品b置于烘箱内,100-110℃下持续烘干24h后称重。
进一步地,利用全自动覆压渗透率-孔隙度联测仪测定柱塞样b的孔隙度和渗透率。
在进行溶蚀模拟实验前,对所述柱塞样a或所述柱塞样b进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样a或所述柱塞样b在溶蚀模拟实验前的孔隙结构特征;
在进行溶蚀模拟实验后,对所述柱塞样b或所述柱塞样a进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样b或所述柱塞样a在溶蚀模拟实验前的孔隙结构特征。
在进行溶蚀模拟实验前后对所述柱塞样进行核磁共振测定分析所述柱塞样孔隙结构特征的具体方法为:
对所述柱塞样进行抽真空处理,并且真空条件下对所述柱塞样充注饱和蒸馏水直至压力大于常压,如15MPa,并维持该压力指定时间后测定所述柱塞样的T2谱;
对所述柱塞样进行多离心参数的离心处理后再次测定所述柱塞样的T2谱;
其中:
所述离心参数包括离心力和离心时间,在离心处理时通过耦合拟合的方式确定所述柱塞样的最佳离心参数;
基于分形维数分析所述柱塞样在核磁共振中三维孔隙结构的复杂程度,所述分形维数记为D:
logVp=(3-D)log(T2)+(D-3)log(T2max);
其中:
T2max是T2的最大值,Vp是由等于或小于T2所对应孔径所占孔隙体积占总孔隙体积的体积分数。
优选地,计算所述柱塞样的T2截止值,以T2截止值为分界计算分段的分形维数。
对比反应前后参数变化,定量研究白云岩三维溶蚀改造效果。其具体方法为:
根据样品离子浓度数据,与所配制初始离子浓度进行比较,计算样品饱和溶蚀量,分析平衡体系离子组分;根据样品重量变化数据,计算样品重量变化率;对比溶蚀前后样品核磁共振T2图谱信号幅度变化,分析白云岩样品受溶蚀改造后的孔径分布、孔喉连通性变化;对比溶蚀前后样品的整体分形维数、分段分形维数变化,分析白云岩受溶蚀改造后的孔隙结构复杂程度演化方向。
结合在一定埋藏深度下白云岩样品受溶蚀改造后的饱和溶蚀量、质量变化率、核磁共振T2信号幅度、分型维数等参数变化,定量研究白云岩三维溶蚀改造效果。
综上所述,本发明具有以下优点:
本发明在研究区不同地层深度条件下,通过溶蚀模拟实验,结合薄片观察、核磁共振T2图谱测量等测试分析手段,实现了白云岩三维孔隙结构溶蚀改造效果定量表征和影响因素分析,从而为实际储层开发提供理论依据。
本发明利用分形维数理论,定量分析了白云岩三维孔隙结构复杂程度受溶蚀改造后的演化方向。
本发明能够准确解析溶蚀孔隙形成过程,认识溶蚀孔隙发育规律,有助于提升储层预测可靠性。
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。

Claims (7)

1.一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,包括如下步骤:
制备白云岩的柱塞样,并对柱塞样进行预处理后称重,测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;
构建模拟条件,基于所述模拟条件对所述柱塞样进行溶蚀模拟实验,收集反应后的生产液;
溶蚀模拟实验后对所述柱塞样进行预处理后称重,并再次测定所述柱塞样的孔渗物性,且通过核磁共振测定分析所述柱塞样的孔隙结构特征;
对所述柱塞样进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样溶蚀模拟实验后的孔隙结构特征;
其中:
所述柱塞样至少包括两组平行样品,所述平行样品的制备具体包括如下步骤:
在岩样相邻位置钻取至少两个呈圆柱状的柱塞样,沿其中一个所述柱塞样的中轴线切分为柱塞样a和柱塞样b,其中所述柱塞样a和所述柱塞样b为一组平行样品,所述柱塞样a和所述柱塞样b的镜像相似度不低于90%;
在进行溶蚀模拟实验前,对所述柱塞样a进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样a在溶蚀模拟实验前的孔隙结构特征;
在进行溶蚀模拟实验后,对所述柱塞样b进行三维切片形成多个铸体薄片,并基于所述铸体薄片观察所述柱塞样b在溶蚀模拟实验后的孔隙结构特征。
2.根据权利要求1所述的一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,在进行溶蚀模拟实验前后对所述柱塞样进行核磁共振测定分析所述柱塞样孔隙结构特征的具体方法为:
对所述柱塞样进行抽真空处理,并且真空条件下对所述柱塞样充注饱和蒸馏水直至压力大于常压,并维持该压力指定时间后测定所述柱塞样的T2谱;
对所述柱塞样进行多离心参数的离心处理后再次测定所述柱塞样的T2谱;
其中:
所述离心参数包括离心力和离心时间,在离心处理时通过耦合拟合的方式确定所述柱塞样的最佳离心参数;
基于分形维数分析所述柱塞样在核磁共振中三维孔隙结构的复杂程度,所述分形维数记为D:
logVp=(3-D)log(T2)+(D-3)log(T2max);
其中:
T2max是T2的最大值,Vp是由等于或小于T2所对应孔径所占孔隙体积占总孔隙体积的体积分数。
3.根据权利要求2所述的一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,计算所述柱塞样的T2截止值,以T2截止值为分界计算分段的分形维数。
4.根据权利要求1-3任一项所述的一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,构建模拟条件包括模拟地层水类型、模拟实验温压条件和模拟实验流速。
5.根据权利要求4所述的一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,模拟地层水类型具体包括如下步骤:
采用苏林分类法计算原生水型的特性系数,获得研究区的地层水类型,对应配制不同酸性介质混合流体;
模拟实验温压条件具体包括如下步骤:
对研究区的地层埋藏史、温压场进行分析,获得研究区实际地质条件的温压梯度,确定溶蚀模拟实验所需的温压条件;
其中:
温度通过地温梯度计算法推算所取样品对应深度的温度;
压力采用静水压力公式计算法推算所取样品对应深度的压力。
6.根据权利要求5所述的一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,当研究区地层中有碳酸盐岩矿物包裹体数据时,以碳酸盐岩矿物包裹体为对象模拟包裹体成分数据,对应配制不同酸性介质混合流体。
7.根据权利要求5或6所述的一种基于核磁共振及三维切片的溶蚀改造量化研究方法,其特征在于,模拟实验流速具体包括如下步骤:
设定溶蚀模拟实验所需的地层水类型以及温压条件不变,按照溶蚀模拟实验的条件改变实验流速进行多次验证实验,在验证实验过程中以指定的时间间隔收集水样,并对水样进行离子成分分析,当水样的离子成分不再发生变化时,判断反应体系达到平衡状态,对应的流速即为溶蚀模拟实验的最佳反应流速。
CN202211620990.XA 2022-12-16 2022-12-16 一种基于核磁共振及三维切片的溶蚀改造量化研究方法 Active CN115824939B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211620990.XA CN115824939B (zh) 2022-12-16 2022-12-16 一种基于核磁共振及三维切片的溶蚀改造量化研究方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211620990.XA CN115824939B (zh) 2022-12-16 2022-12-16 一种基于核磁共振及三维切片的溶蚀改造量化研究方法

Publications (2)

Publication Number Publication Date
CN115824939A CN115824939A (zh) 2023-03-21
CN115824939B true CN115824939B (zh) 2023-10-13

Family

ID=85516168

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211620990.XA Active CN115824939B (zh) 2022-12-16 2022-12-16 一种基于核磁共振及三维切片的溶蚀改造量化研究方法

Country Status (1)

Country Link
CN (1) CN115824939B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107727679A (zh) * 2017-11-03 2018-02-23 中国科学院地质与地球物理研究所 一种表征深层碳酸盐岩岩石物理学特征方法
CN109709301A (zh) * 2018-11-30 2019-05-03 中国石油天然气股份有限公司 一种裂缝孔隙型致密砂岩储层分布确定方法、装置及系统
CN110687612A (zh) * 2019-09-17 2020-01-14 中国石油天然气股份有限公司 吸附油和游离油含量连续表征的页岩油分析方法及装置
CN111595764A (zh) * 2020-04-26 2020-08-28 中国石油天然气股份有限公司 一种模拟硫酸根离子浓度对碳酸盐岩溶蚀影响的实验方法
CN112992283A (zh) * 2021-02-07 2021-06-18 中国石油天然气股份有限公司 晶体尺度的白云岩溶蚀孔隙形成演化模拟方法及装置
CN113624799A (zh) * 2021-08-20 2021-11-09 西南石油大学 一种基于核磁共振和分形维数的岩石渗透率预测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107727679A (zh) * 2017-11-03 2018-02-23 中国科学院地质与地球物理研究所 一种表征深层碳酸盐岩岩石物理学特征方法
CN109709301A (zh) * 2018-11-30 2019-05-03 中国石油天然气股份有限公司 一种裂缝孔隙型致密砂岩储层分布确定方法、装置及系统
CN110687612A (zh) * 2019-09-17 2020-01-14 中国石油天然气股份有限公司 吸附油和游离油含量连续表征的页岩油分析方法及装置
CN111595764A (zh) * 2020-04-26 2020-08-28 中国石油天然气股份有限公司 一种模拟硫酸根离子浓度对碳酸盐岩溶蚀影响的实验方法
CN112992283A (zh) * 2021-02-07 2021-06-18 中国石油天然气股份有限公司 晶体尺度的白云岩溶蚀孔隙形成演化模拟方法及装置
CN113624799A (zh) * 2021-08-20 2021-11-09 西南石油大学 一种基于核磁共振和分形维数的岩石渗透率预测方法

Also Published As

Publication number Publication date
CN115824939A (zh) 2023-03-21

Similar Documents

Publication Publication Date Title
US10161891B1 (en) Method for characterizing rock physical characteristics of deeply buried carbonate rocks
SA516370672B1 (ar) جهاز ثلاثي المحاور لاختبار الرنين النووي المغناطيسي (nmr)
US20180335494A1 (en) Systems and methods for tri-axial nmr testing
CN108627533A (zh) 一种测定多孔介质中流体动用特征的核磁共振实验方法及装置
CN105866002B (zh) 一种精确的含油页岩核磁共振孔隙度测试方法
CN106093299A (zh) 一种致密气储层钻井液伤害评价实验方法
EA029719B1 (ru) Система и компьютерно-реализуемый способ определения свойств смачиваемости флюидсодержащей пористой среды
CN103257151A (zh) 一种定量评价油气二次运移过程中孔喉动用规律的方法
CN113075102B (zh) 一种建立多孔介质自发渗吸量与时间关系数学模型的方法
CN106525681A (zh) 一种确定页岩储层孔隙直径大小的方法
US20240027379A1 (en) Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids
CN209821099U (zh) 基于核磁共振的多功能致密气储层动态参数联测装置
CN115824939B (zh) 一种基于核磁共振及三维切片的溶蚀改造量化研究方法
CN109901238A (zh) 一种基于应力差电阻率实验的高应力地层电阻率校正方法
CN115389387A (zh) 一种评价岩心伤害的实验方法
WO2021206822A9 (en) Systems and methods for slice selective nuclear magnetic resonance testing of fractured core plugs to determine in-situ pore volume
CN113484153A (zh) 室内真三轴水压致裂地应力测试模拟方法及装置
CN109709130B (zh) 一种测试全油基钻井液地层含油率的方法
JP2021503593A (ja) 全偏波レーダーに基づく氷床の内部組織構造と氷流れ場の分布特徴の識別方法
CN111693427B (zh) 油气藏流体可动性的分析方法
CN110410058B (zh) 一种校正岩心实验结果刻度二维核磁测井的方法
CN110441204B (zh) 一种基于数字岩心模拟的致密储层压裂液伤害数字化评价方法
Falcon-Suarez et al. Experimental rig to improve the geophysical and geomechanical understanding of CO2 reservoirs
CN116046657A (zh) 一种定量研究含膏白云岩储层溶蚀增孔能力的实验方法
CN112924360B (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