发明名称 岩石微裂缝损伤变量的核磁共振定量分析方法
摘要 岩石微裂缝损伤变量的核磁共振定量分析方法,先加工试样,将试样浸泡在蒸馏水中,然后进行核磁共振检测,通过对试样进行CPMG脉冲序列测试,得到自旋回波串衰减信号,对自旋回波串衰减信号进行反演,得到T<sub>2</sub>谱的分布图、谱面积值;分析了弛豫时间T<sub>2</sub>和孔隙尺寸r的对应关系,根据谱峰面积正比于微裂缝孔隙数量,本发明建立由核磁共振横向弛豫时间T<sub>2</sub>谱表征的岩石损伤变量公式,为识别岩石材料的损伤提供了定量分析的方法。
申请公布号 CN105241913A 申请公布日期 2016.01.13
申请号 CN201510655268.3 申请日期 2015.10.10
申请人 西安石油大学 发明人 王萍;屈展;王亮;黄海;高辉
分类号 G01N24/08(2006.01)I 主分类号 G01N24/08(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贺建斌
主权项 岩石微裂缝损伤变量的核磁共振定量分析方法,其特征在于,包括以下步骤:第一步,按行业标准将试样加工成直径为25毫米、长为50毫米的圆柱体;将试样浸泡在蒸馏水中,分别在浸泡10min、30min、1h、2h、4h、6h、8h、1d、3d、5d不同时间后取出试样;对经过浸泡的试样进行核磁共振检测,通过对试样进行CPMG脉冲序列测试,得到自旋回波串衰减信号;利用核磁共振弛豫时间反演拟合软件对自旋回波串衰减信号进行反演,得到T<sub>2</sub>谱的分布图、谱面积值和各峰所占比例值,利用MiniMR核磁共振成像软件得到同一试样不同浸泡时间的核磁共振成像;第二步,根据核磁共振原理,T<sub>2</sub>弛豫由表面弛豫决定,表面驰豫与介质表面积有关,多孔介质孔隙表面积S与孔隙体积V之比为介质比表面,介质比表面越大,则驰豫越强,反之亦然,T<sub>2</sub>表面表示为:<img file="FDA0000818947390000011.GIF" wi="814" he="157" />式中ρ<sub>2</sub>为T<sub>2</sub>表面驰豫强度;(S/V)<sub>孔隙</sub>为孔隙表面积S与孔隙体积V之比,继而得到T<sub>2</sub>与孔径r的关系为:<img file="FDA0000818947390000012.GIF" wi="838" he="142" />Fs称为几何形状因子,由上式可见孔隙内流体的弛豫时间和孔隙空间大小有关,孔隙越小,比面积越大,表面相互作用的影响越强烈,T<sub>2</sub>时间也越短,弛豫时间T<sub>2</sub>和孔径r是一一对应的,利用T<sub>2</sub>分布来评价孔隙裂纹大小、分布;第三步,建立微裂纹密度分布函数n(r,t),对理想微裂纹系统进行统计描述,n表示在t时刻时,单位体积内尺度在r~r+dr范围内的裂纹数n(r,t)dr,在损伤初期,略去裂纹间相关效应,成核过程和扩展过程决定了n的演化规律,考察孔径r所张成的一维相空间中相体元(r,r+Δr)内微裂纹数密度变化,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mo>&Integral;</mo><mi>r</mi><mrow><mi>r</mi><mo>+</mo><mi>&Delta;</mi><mi>r</mi></mrow></msubsup><mi>n</mi><mrow><mo>(</mo><msup><mi>r</mi><mo>&prime;</mo></msup><mo>,</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>dr</mi><mo>&prime;</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818947390000021.GIF" wi="1130" he="134" /></maths>第四步,设连续损伤D用微损伤数密度来表示:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mi>h</mi><mi>&infin;</mi></msubsup><mi>r</mi><mi>n</mi><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mi>d</mi><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818947390000022.GIF" wi="1214" he="119" /></maths>将第二步中的式(2)中的r代入上式(4)则得<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>T</mi><mn>0</mn></msub><msub><mi>T</mi><mn>2</mn></msub></msubsup><msub><mi>F</mi><mi>S</mi></msub><msub><mi>T</mi><mn>2</mn></msub><mi>n</mi><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>T</mi><mo>)</mo></mrow><msub><mi>dT</mi><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818947390000023.GIF" wi="1211" he="125" /></maths>第五步,对第一步反演得到的T<sub>2</sub>谱面积进行分析,核磁共振横向弛豫时间T<sub>2</sub>谱的积分面积正比于试样中所含流体的多少,T<sub>2</sub>谱中X坐标表示弛豫时间T<sub>2</sub>值大小,谱峰的位置与孔径大小有关,由此得到微损伤扩展前沿在Δt内新扩展的面积S≈∫F<sub>S</sub>n(t,T<sub>2</sub>)dT<sub>2</sub>      (6)第六步,根据(5)和(6)两式,建立以T<sub>2</sub>谱表征的损伤变量为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>D</mi><mi>T</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>S</mi><msub><mi>T</mi><mn>2</mn></msub></msub><mo>-</mo><msub><mi>S</mi><msub><mi>T</mi><mn>0</mn></msub></msub></mrow><msub><mi>S</mi><msub><mi>T</mi><mn>2</mn></msub></msub></mfrac><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>S</mi><msub><mi>T</mi><mn>0</mn></msub></msub><msub><mi>S</mi><msub><mi>T</mi><mn>2</mn></msub></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818947390000024.GIF" wi="1069" he="143" /></maths>S<sub>T0</sub>由核磁共振测量初始损伤时的T<sub>2</sub>谱积分面积,S<sub>T2</sub>为经过水化浸泡后由核磁共振测量水化损伤的T<sub>2</sub>谱积分面积。
地址 710065 陕西省西安市电子二路东段18号