发明名称 一种多孔介质内流体速度与温度同步测量的磁共振成像方法
摘要 一种多孔介质内流体速度与温度同步测量的磁共振成像方法,属于磁共振成像测量领域。本方法采用核磁共振翻转恢复标签方法,首先用自旋回波序列两点法测某一温度下被测流体纵向弛豫时间T1;在不同温度下重复步骤一得到多个温度下的T1值,拟合得到T1与温度T的关系;用翻转恢复标签序列对静态时该流体进行测量,再对任意温度流速下该流体进行测量;用matlab读入图像,得到图像信号强度SI分布,根据SI与T的关系计算温度T;通过比较该流速下标签坐标和静态时标签的坐标求出翻转恢复时间ti内的标签位移,进而求流体流速。本发明具有测量速度快,能有效补偿传统翻转恢复方法测量流场温度时由于流动导致的信号降低。
申请公布号 CN103776490A 申请公布日期 2014.05.07
申请号 CN201410014755.7 申请日期 2014.01.13
申请人 大连理工大学 发明人 刘瑜;周欣欢;宋永臣;蒋兰兰;赵越超;杨明军;赵佳飞;张毅;王大勇
分类号 G01D21/02(2006.01)I 主分类号 G01D21/02(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉
主权项 1.一种多孔介质和大体积内流体速度与温度同步测量的磁共振成像方法,其特征在于:(1)自旋回波序列两点法测量某一温度下流体T<sub>1</sub>值使用两个多层自旋回波序列SEMS测量某一温度下流体的信号强度值,两个自旋回波序列的重复时间TR分别为TR<sub>1</sub>和TR<sub>2</sub>,其他参数设置完全相同,对应视野区域内平均信号强度分别为SI<sub>1</sub>,SI<sub>2</sub>;自旋回波序列信号强度SI为:<maths num="0001"><![CDATA[<math><mrow><mi>SI</mi><mo>=</mo><msub><mi>M</mi><mn>0</mn></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>TE</mi><mo>/</mo><msub><mi>T</mi><mn>2</mn></msub></mrow></msup><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mrow><mn>2</mn><mi>e</mi></mrow><mrow><mo>-</mo><mrow><mo>(</mo><mi>TR</mi><mo>-</mo><mi>TE</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></msup><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>TR</mi><mo>/</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中M<sub>0</sub>为纵向平衡磁化强度,TE为回波时间,T<sub>2</sub>为横向弛豫时间;由于TE&lt;&lt;TR及TE&lt;&lt;T<sub>2</sub>,上式化简为<img file="FDA0000456234170000013.GIF" wi="444" he="64" />对两个不同TR值下的自旋回波信号强度做商:<maths num="0002"><![CDATA[<math><mrow><mfrac><msub><mi>SI</mi><mn>1</mn></msub><msub><mi>SI</mi><mn>2</mn></msub></mfrac><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><msub><mrow><mo>-</mo><mi>TR</mi></mrow><mn>1</mn></msub><mo>/</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><msub><mrow><mo>-</mo><mi>TR</mi></mrow><mn>2</mn></msub><mo>/</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>已知该温度下的两个信号强度SI<sub>1</sub>,SI<sub>2</sub>和两个序列的重复时间TR<sub>1</sub>,TR<sub>2</sub>,求解(2)式得到该温度下的T<sub>1</sub>值;(2)获得T<sub>1</sub>和温度T关系在不同温度下,重复步骤(1)测量计算多个温度点下的T<sub>1</sub>值;然后用不同温度下的T<sub>1</sub>值与温度T利用公式(3)进行数据拟合,从而得到常数A,B,C的值:T<sub>1</sub>=Ae<sup>-B/T</sup>+C    (3)(3)使用翻转恢复标签序列测量流体速度和温度用翻转恢复标签序列对感兴趣流动区域内静态待测流体进行测量,再进行某一温度和流速下该流体的测量,得到两幅图的信号强度即灰度级分布图;翻转恢复标签序列得到的灰度图上标签处温度与信号强度关系为:<maths num="0003"><![CDATA[<math><mrow><mi>T</mi><mo>=</mo><mfrac><mi>B</mi><mrow><mi>ln</mi><mo>[</mo><mo>-</mo><mi>A</mi><mi>ln</mi><mi>&alpha;</mi><mo>/</mo><mrow><mo>(</mo><mi>ti</mi><mo>+</mo><mi>C</mi><mi>ln</mi><mi>&alpha;</mi><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mfrac><msub><mi>SI</mi><mi>Z</mi></msub><msub><mi>SI</mi><mn>0</mn></msub></mfrac><mo>|</mo></mrow></math>]]></maths>上式中SI<sub>Z</sub>为标签上的信号强度,SI<sub>0</sub>为平衡信号强度,ti为翻转恢复标签序列中的参数恢复时间;使用图像处理工具得到步骤(3)中静态和某一流速下的两幅图的信号强度分布,得到某一流速下翻转恢复标签序列方法中的标签位置坐标,与静态得到的翻转恢复标签序列中标签位置坐标进行对比,根据ti内的标签变形,计算出流体流速。
地址 116024 辽宁省大连市凌工路2号