发明名称 拟饱含水核磁共振自旋回波信号的分解与合成方法
摘要 “拟饱含水核磁共振自旋回波信号的分解与合成方法”。根据已知的饱和度数据和对应的弛豫信号,优化拟合确定不同油样饱和度对弛豫过程的影响系数。对于不同含油饱和度的弛豫信号,用同体积的盐水取代其中的油样,通过引入所述的影响系数,表征不同含油饱和度油样对弛豫信号的影响,据此写出不同含油饱和度下拟饱含水核磁共振弛豫信号。反演拟饱含水岩心的自旋回波弛豫信号,可以在一定程度上消除不同流体类型和不同饱和度的影响,提高核磁共振信号对岩石孔隙结构和流体的识别能力。本发明实现拟饱含水岩心核磁共振自旋回波信号的合成与分解,达到提高核磁共振数据评价岩石孔隙结构与流体识别的精度与可靠性。
申请公布号 CN103487837B 申请公布日期 2016.11.23
申请号 CN201310418890.3 申请日期 2013.09.13
申请人 同济大学 发明人 刘堂晏;邬龙;薛苗苗;王红涛;潘豪杰;王飞飞
分类号 G01V3/14(2006.01)I;G01V3/38(2006.01)I 主分类号 G01V3/14(2006.01)I
代理机构 上海天协和诚知识产权代理事务所 31216 代理人 叶凤
主权项 一种拟完全饱含盐水核磁共振自旋回波信号的分解与合成方法,其特征在于,该方法包括实验流程和数据分析两个部分,其中,(一)所述实验流程,即进行的核磁共振实验包括步骤有:1.1、准备好3块实验岩心,分别进行预处理,包括洗油,洗盐,烘干,抽真空;1.2、配置矿化度为80000ppm的实验盐水;1.3、针对配置的实验盐水和实验油样,分别测量实验盐水和实验油样的自由横向弛豫时间,测得实验盐水和实验油样的自由横向弛豫时间分别为T<sub>2w</sub>=773ms和T<sub>20</sub>=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),<maths num="0001"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>M</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mi>o</mi></msub><msub><mi>M</mi><mi>o</mi></msub><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>t</mi><msub><mi>T</mi><mrow><mn>2</mn><mi>o</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>S</mi><mi>o</mi></msub><msub><mi>M</mi><mi>o</mi></msub><mi>exp</mi><mo>(</mo><mrow><mo>-</mo><mfrac><mi>t</mi><msub><mi>T</mi><mrow><mn>2</mn><mi>w</mi></mrow></msub></mfrac></mrow><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001056835990000011.GIF" wi="1550" he="146" /></maths>式中,M(t),挖补法恢复的饱含实验盐水总信号,即拟完全饱含实验盐水的自由横向弛豫回波信号,无单位;M<sub>t</sub>(t)饱含实验油样和实验盐水的自由横向弛豫回波总测量信号,无单位;S<sub>0</sub>为含实验油样饱和度;t为测量时间,ms;T<sub>2w</sub>,实验盐水的自由横向弛豫时间,ms;T<sub>20</sub>,实验油样的自由横向弛豫时间,ms;M<sub>o</sub>为岩心的初始自由横向弛豫回波强度,无单位,用第一个回波信号代替;在实验盐水相和实验油样相信号的前面加上分离系数V<sub>c</sub>,方程(1)变成,<maths num="0002"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>M</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>c</mi></msub><msub><mi>S</mi><mi>o</mi></msub><msub><mi>M</mi><mi>o</mi></msub><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>t</mi><msub><mi>T</mi><mrow><mn>2</mn><mi>o</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>V</mi><mi>c</mi></msub><msub><mi>S</mi><mi>o</mi></msub><msub><mi>M</mi><mi>o</mi></msub><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>t</mi><msub><mi>T</mi><mrow><mn>2</mn><mi>w</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001056835990000021.GIF" wi="1526" he="155" /></maths>2.2、确定不同流体自由横向弛豫回波信号的分离系数V<sub>c</sub>方程(2)中的M(t)相当于挖补法确定的岩心饱含实验盐水的自由横向弛豫回波信号,设M(t)与实测饱含实验盐水岩心的自由横向弛豫回波信号M<sub>w</sub>(t)之间的均方差为Q<sub>a</sub>(V<sub>c</sub>),<maths num="0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>V</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2048</mn></munderover><msup><mrow><mo>(</mo><mrow><mi>M</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>M</mi><mi>w</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mn>2048</mn><mo>-</mo><mn>1</mn></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001056835990000022.GIF" wi="1206" he="244" /></maths>t<sub>i</sub>表示采样的时间;方程中分子部分是V<sub>c</sub>的函数,在分子中对V<sub>c</sub>求导,并令导数等于零,得到确定V<sub>c</sub>的方程(4)<maths num="0004"><math><![CDATA[<mrow><msub><mi>V</mi><mi>c</mi></msub><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2048</mn></munderover><mo>&lsqb;</mo><msub><mi>M</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>M</mi><mi>w</mi></msub><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo><mo>&rsqb;</mo></mrow><mrow><msub><mi>S</mi><mi>o</mi></msub><msub><mi>M</mi><mi>o</mi></msub><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2048</mn></munderover><mo>&lsqb;</mo><mi>exp</mi><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>t</mi><mi>i</mi></msub></mrow><msub><mi>T</mi><mrow><mn>2</mn><mi>o</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>-</mo><mi>exp</mi><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>t</mi><mi>i</mi></msub></mrow><msub><mi>T</mi><mrow><mn>2</mn><mi>w</mi></mrow></msub></mfrac><mo>)</mo><mo>&rsqb;</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001056835990000023.GIF" wi="1400" he="360" /></maths>2.3、根据不同实验油样饱和度S<sub>o</sub>的核磁共振实验数据,计算V<sub>c</sub>值,得到实验油样饱和度S<sub>o</sub>与分离系数V<sub>c</sub>的关系:<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>V</mi><mi>c</mi></msub><mo>=</mo><mn>3.943</mn><msubsup><mi>S</mi><mi>o</mi><mn>2</mn></msubsup><mo>-</mo><mn>4.689</mn><msub><mi>S</mi><mi>o</mi></msub><mo>+</mo><mn>3.179</mn></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>=</mo><mn>0.8478</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001056835990000024.GIF" wi="1350" he="197" /></maths>其中,R为数据拟合的相关系数;2.4、计算岩心拟完全饱含实验盐水的自由横向弛豫回波信号;2.5、反演岩心拟饱含实验盐水的自由横向弛豫回波信号,计算岩心的拟饱含实验盐水T2谱分布,根据这个T2谱分布计算岩心的孔隙结构参数。
地址 200092 上海市杨浦区四平路1239号