拟饱含水核磁共振自旋回波信号的分解与合成方法技术

拟饱含水核磁共振自旋回波信号的分解与合成方法。根据已知的饱和度数据和对应的弛豫信号，优化拟合确定不同油样饱和度对弛豫过程的影响系数。对于不同含油饱和度的弛豫信号，用同体积的盐水取代其中的油样，通过引入所述的影响系数，表征不同含油饱和度油样对弛豫信号的影响，据此写出不同含油饱和度下拟饱含水核磁共振弛豫信号。反演拟饱含水岩心的自旋回波弛豫信号，可以在一定程度上消除不同流体类型和不同饱和度的影响，提高核磁共振信号对岩石孔隙结构和流体的识别能力。本发明专利技术实现拟饱含水岩心核磁共振自旋回波信号的合成与分解，达到提高核磁共振数据评价岩石孔隙结构与流体识别的精度与可靠性。

【技术实现步骤摘要】
拟饱含水核磁共振自旋回波信号的分解与合成方法
本专利技术属于石油勘探
，应用于地层岩石孔隙流体的识别和孔隙结构分析。技术背景核磁共振自旋回波信号是核磁共振弛豫谱反演的基础数据，核磁共振弛豫谱可以用于地层岩石孔隙流体的识别和孔隙结构分析。在实际应用中，地层孔隙流体与孔隙结构对自旋回波信号的影响相互叠加，所以，目前直接反演自旋回波信号得到的弛豫谱分布，实际上包含了岩石孔隙流体和孔隙结构相互叠加的信息，导致弛豫谱数据在孔隙流体识别和孔隙结构分析上都出现困难。现有的自旋回波信号反演方法，都是不加区别的反演包含孔隙流体信息和孔隙结构信息的总回波信号。然后，根据反演弛豫谱的分布形式和数值特征分析孔隙流体类型和孔隙结构特征。尽管通过改进反演算法，可以提高弛豫谱分布的精度和稳定性，但是，由于没有实现孔隙流体信息和孔隙结构信息的分离反演，所以，无法从根本上提高识别孔隙流体和分析孔隙结构的精度。
技术实现思路
本专利技术目的在于公开一种核磁共振弛豫信号的分解与合成方法，这种方法可以提高核磁共振数据对岩石孔隙结构和流体类型的识别能力。概括的说，首先在实验流程中分别测量研究流体样品(油样和盐水)的横向自由弛豫时间，配置不同饱和度的油样，并测量在不同油样饱和度下岩样的横向弛豫时间；实验获得数据后，在数据处理中采用“挖补法”，将不同含油饱和度下岩样的弛豫信号与完全饱含水岩样的弛豫信号进行优化拟合，得到不同饱和度油样对岩石总弛豫信号的影响系数(即流体分离系数)，如此，通过分解合成得到不同饱和度下岩石拟完全含水弛豫信号。本专利技术获得的岩石拟完全含水弛豫信号，在一定程度上消除了不同流体类型和饱和度对总弛豫信号的影响，提高岩石弛豫信号对岩石孔隙结构和流体的识别能力。本专利技术通过以下技术方案实现：一种拟饱含水核磁共振自旋回波信号的分解与合成方法，其特征在于，该方法包括实验流程和数据分析两个部分，其中，(一)所述实验流程，即进行的核磁实验包括步骤有：1.1、准备好若干实验岩芯、一定矿化度的盐水和实验油样，并对岩芯进行必要的预处理，完成实验前的准备。(该部分属本领域普通技术人员具备的公知常识以及掌握的实验常规操作，属于现有技术)1.2、测量或者查询出实验配置盐水和实验油样(比如变压器油，或煤油)的自由横向弛豫时间。1.3、分不同饱和度阶段，向岩心中注入油样，驱替原先占据孔隙的部分盐水。在各个不同的饱和度阶段，记录下它们饱和度数据，同时测量对应的不同盐水饱和度岩芯的弛豫回波信号。继而，(二)所述数据分析包括步骤有：2.1、“挖补法”用以消除油相驰豫信号使用挖补法，去掉含油的体积，用相同的体积的水替代，将饱含油样和盐水的弛豫信号恢复成拟完全饱含盐水的弛豫信号，即方程(2)，式中，M(t)，挖补法恢复的饱含水总信号，也即拟完全含水弛豫信号，无单位；T2w，水的体积弛豫时间，ms；Mo用第一个回波信号代替；在水相和油相信号的前面加上修正系数Vc，方程2变成，2.2、确定不同流体弛豫信号的分离系数Vc方程3中的M(t)相当于挖补法确定的岩心饱含水的弛豫信号，设M(t)与实测饱含水岩心弛豫信号Mw(t)之间的均方差为Qa(Vc)，方程中分子部分是Vc的函数。在分子中对Vc求导，并令导数等于零，得到确定Vc的方程(方程5)，2.3、根据不同饱和度So的核磁共振实验数据，计算Vc值，得到含油饱和度(So)与Vc的关系：2.4、计算岩心拟完全含水的弛豫信号。2.5、反演岩样拟饱含水弛豫信号，计算岩样的拟饱含水T2谱分布，根据这个T2谱分布计算岩样的孔隙结构参数。上述技术方案实现了拟饱含水岩心核磁共振自旋回波信号的合成与分解，达到提高核磁共振数据评价岩石孔隙结构与流体识别的精度与可靠性。其中公开的实用化的算法，分离核磁共振自旋回波信号中的包含的孔隙流体信息与孔隙结构信息，分别得到孔隙流体的自旋回波信号与孔隙结构的自旋回波信号，并分别反演孔隙流体自旋回波信号和孔隙结构自旋回波信号，得到孔隙流体的弛豫谱分布和孔隙结构的弛豫谱。在此基础上，分别使用孔隙流体的弛豫谱和孔隙结构的弛豫谱，研究岩石的孔隙流体和孔隙结构，提高对岩石孔隙识别和孔隙结构分析的精度。附图说明图1本专利技术整个方法实施流程示意图图2为“挖补法”示意图图3为研究流体的自由横向弛豫时间测量图图4为Zhen53-240样品饱含盐水的弛豫回波及T2谱反演结果(A：弛豫回波，B：反弹T2谱分布)图5为Zhen53240盐岩不同含水饱和度的T2谱分布图6为不同饱和度弛豫回波信号分离后的反演结果对比(Zhen-53-241)图7为信号分离算法的软件化编程框图具体实施方式实施例1(理论依据)氢核在岩石孔隙中发生的弛豫过程，就是在外加磁场的作用下，氢核与孔隙壁之间相互作用而发生的能量传递和转化的过程。氢核在孔隙中的弛豫过程，与氢核弛豫的边界条件有关，而这个边界条件由孔隙内壁构成，所以，氢核的弛豫过程与岩石的孔隙结构有关，这是核磁共振弛豫信号评价岩石孔隙结构的基础。另一方面，岩石弛豫信号显著受到岩石中所有流体的影响，包括盐水和油气。而能够反映岩石孔隙结构的弛豫信号，主要来源于孔隙壁附近的流体。实际的计算结果表明，当岩石孔隙中存在多种流体时，各种流体的影响叠加在岩石样品的总弛豫信号上，会引起岩石弛豫信号的严重畸变，甚至完全掩盖孔隙结构对弛豫信号的影响。如果不剔除多种流体叠加在弛豫信号上的影响，将会严重降低弛豫信号及相应的反演结果对岩石孔隙结构的评价精度。本次专利技术对不同含油饱和度下的弛豫信号进行对比分析，得到不同含油饱和度下油样对总弛豫信号的影响程度，确定油样弛豫信号的分离系数，采用“挖补法”，恢复岩样拟完全含水弛豫信号，在一定程度上消除不同流体类型和不同含油饱和度对总弛豫信号的影响，提高弛豫信号对岩样孔隙结构和流体的识别精度。图1为本专利技术的实施流程图。步骤1，分别测量和查询研究流体的自由横向弛豫时间(图2)。对于不同矿化度的盐水，在必要时，需要采用杯状容器，测量盐水在各种矿化度下的自由横向弛豫时间。核磁共振实验仪器采用低场核磁共振分析仪器，例如，牛津公司的MARAN2S核磁分析仪器，采用CPMG(Carr-Purcel1-Meiboom-Gil1)脉冲序列，测量流体的自旋回波信号，并用SVD(SingularvalueDecomposition，奇异值分解)方法反演回波信号，获得研究自由流体样品的T2谱分布。除非明确说明，以下的核磁共振测量方法和数据反演方法与此类似。本次研究作了煤油本文档来自技高网...

【技术保护点】
一种拟饱含水核磁共振自旋回波信号的分解与合成方法，其特征在于，该方法包括实验流程和数据分析两个部分，其中，(一)所述实验流程，即进行的核磁实验包括步骤有：1.1、准备好若干实验岩芯、一定矿化度的盐水和实验油样，并对岩芯进行预处理；1.2、测量或者查询出实验配置盐水和实验油样的自由横向弛豫时间；1.3、分不同饱和度阶段，向岩心中注入油样，驱替原先占据孔隙的部分盐水，在各个不同的饱和度阶段，记录下它们饱和度数据，同时测量对应的不同盐水饱和度岩芯的弛豫回波信号；继而，(二)所述数据分析包括步骤有：2.1、“挖补法”用以消除油相驰豫信号使用挖补法，去掉含油的体积，用相同的体积的水替代，将饱含油样和盐水的弛豫信号恢复成拟完全饱含盐水的弛豫信号，即方程(2)，M(t)=Mt(t)-SoMoexp(-tT2o)+SoMoexp(-tT2w)---(2)式中，M(t)，挖补法恢复的饱含水总信号，也即拟完全含水弛豫信号，无单位；T2w，水的体积弛豫时间，ms；Mo用第一个回波信号代替；在水相和油相信号的前面加上修正系数Vc，方程2变成，M(t)=Mt(t)-VcSoMoexp(-tT2o)+VcSoMoexp(-tT2w)---(3)2.2、确定不同流体弛豫信号的分离系数Vc方程3中的M(t)相当于挖补法确定的岩心饱含水的弛豫信号，设M(t)与实测饱含水岩心弛豫信号Mw(t)之间的均方差为Qa(Vc)，Qa(Vc)=Σi=12048(M(ti)-Mw(ti))22048-1---(4)方程中分子部分是Vc的函数，在分子中对Vc求导，并令导数等于零，得到确定Vc的方程(5)Vc=Σi=12048[Mt(ti)-Mw(ti)]SoMoΣi=12048[exp(-tiT2o)-exp(-tiT2w)---(5)2.3、根据不同饱和度So的核磁共振实验数据，计算Vc值，得到含油饱和度(So)与Vc的关系：Vc=3.943So2-4.689So+3.179R2=0.8478---(6)2.4、计算岩心拟完全含水的弛豫信号；2.5、反演岩样拟饱含水弛豫信号，计算岩样的拟饱含水T2谱分布，根据这个T2谱分布计算岩样的孔隙结构参数。...

【技术特征摘要】
1.一种拟完全饱含盐水核磁共振自旋回波信号的分解与合成方法，其特征在于，该方法包括实验流程和数据分析两个部分，其中，(一)所述实验流程，即进行的核磁共振实验包括步骤有：1.1、准备好3块实验岩心，分别进行预处理，包括洗油，洗盐，烘干，抽真空；1.2、配置矿化度为80000ppm的实验盐水；1.3、针对配置的实验盐水和实验油样，分别测量实验盐水和实验油样的自由横向弛豫时间，测得实验盐水和实验油样的自由横向弛豫时间分别为T2w＝773ms和T20＝875ms；1.4、用饱和实验盐水饱和经过预处理的岩心，经过抽真空处理后的岩心，排除了气泡的影响，岩心完全被实验盐水饱和；1.5、对饱和实验盐水的岩心测量自由横向弛豫回波信号，测量参数：自由横向弛豫回波时间间隔分别为Te＝0.2ms，0.4ms，核磁共振等待时间分别为Tw＝1，3，6，8ms，计算结果表明，Te＝0.2ms，Tw＝3ms的测量参数组合，可以测量到比较真实有效的弛豫回波信号，此时，岩心的饱含实验盐水的饱和度为100％；1.6、分不同实验盐水饱和度阶段，向岩心中注入实验油样，驱替原先占据孔隙的部分实验盐水，在各个不同的实验盐水饱和度阶段，记录下它们饱和度数据，同时测量对应的不同实验盐水饱和度岩心的自由横向弛豫回波信号；(二)所述数据分析包括步骤有：2.1、“挖补法”用以消除实验油样相自由横向驰豫回波信号使用挖补法，去掉含实验油样的体积，用相同的体积的实验盐水替代，将饱含实验油样和实验盐水的自由横向弛豫回波信号恢复成拟完全饱含实验盐水的自由横向弛豫回波信号，即方程(1)，式中，M(t)，挖补法恢复的饱含实验盐水总信号，即拟完全饱含实验盐水的自由横向弛豫回波信号，无单位；Mt(...

【专利技术属性】
技术研发人员：刘堂晏，邬龙，薛苗苗，王红涛，潘豪杰，王飞飞，
申请(专利权)人：同济大学，
类型：发明
国别省市：