磁共振成像原理(理論)4:自由進動和弛豫 (Free Precession and Relaxation)

當磁化自旋系統被射頻脈沖擾動而偏離其熱平衡態后,一旦移除外部激勵并給予足夠時間,系統將根據熱力學定律返回平衡態。這一過程包含三個特征現象:
(a) 自由進動——宏觀磁化矢量 (M?\vec{M}M) 繞( B0?\vec {B_0}B0?? )場的進動;
(b) 縱向弛豫——縱向磁化分量(Mz{M_z}Mz?) 的恢復;
? 橫向弛豫——橫向磁化分量 (Mxy{M_{xy}}Mxy?) 的衰減。

弛豫過程源于原子核周圍由隨機熱運動產生的時變微觀磁場,但具體機制因自旋系統而異,過于復雜。實際上不管是激發還是弛豫,只有量子力學才能完全精確的給出數學解釋,此處僅通過布洛赫方程進行唯象描述,它提供足夠的精度,對MR成像而言也足夠了

忽略射頻激勵項后,在拉莫爾旋轉坐標系中,弛豫過程由以下方程描述:
{dMz′dt=?Mz′?Mz0T1dMx′y′dt=?Mx′y′T2(3.121) \left\{ \begin{array}{l} \frac{d M_{z'}}{d t} = -\frac{M_{z'} - M_z^0}{T_1} \\ \frac{d M_{x'y'}}{d t} = -\frac{M_{x'y'}}{T_2} \end{array} \right. \tag {3.121} {dtdMz??=?T1?Mz??Mz0??dtdMxy??=?T2?Mxy???(3.121)

其中 Mz0M_z^0Mz0? 為熱平衡態縱向磁化強度,T1T_1T1?T2T_2T2? 分別為描述縱向和橫向弛豫過程的時間常數。解方程 (3.121) 得:
{Mx′y′(t)=Mx′y′(0+)e?t/T2Mz′(t)=Mz0(1?e?t/T1)+Mz′(0+)e?t/T1(3.122) \left\{ \begin{array}{l} M_{x'y'}(t) = M_{x'y'}(0_+) e^{-t/T_2} \\ M_{z'}(t) = M_z^0 (1 - e^{-t/T_1}) + M_{z'}(0_+) e^{-t/T_1} \end{array} \right. \tag {3.122} {Mxy?(t)=Mxy?(0+?)e?t/T2?Mz?(t)=Mz0?(1?e?t/T1?)+Mz?(0+?)e?t/T1??(3.122)
此處 Mx′y′(0+)M_{x'y'}(0_+)Mxy?(0+?)Mz′(0+)M_{z'}(0_+)Mz?(0+?) 分別表示射頻脈沖結束后,旋轉坐標系下觀測到的初始橫向和縱向磁化分量。
image.png

值得注意的是,T1T_1T1?T2T_2T2? 并不是縱向和橫向弛豫的完成時間。為了說得更加清楚,假設熱平衡態下產生一個 x′x'x 軸90度射頻脈沖,在射頻脈沖激發前(t=0?)(t=0_{-})(t=0??)和射頻脈沖激發后(t=0+)(t=0_{+})(t=0+?)的橫向和縱向磁化分量如下所示:
{Mxy′(0?)=0Mz′(0?)=Mz0→90x′{Mxy′(0+)=Mz0Mz′(0+)=0 \begin{cases} M_{xy^{\prime}}(0_{-}) = 0 \\ M_{z^{\prime}}(0_{-}) = M_{z}^{0} \end{cases} \xrightarrow{90_{x^{\prime}}} \begin{cases} M_{xy^{\prime}}(0_{+}) = M_{z}^{0} \\ M_{z^{\prime}}(0_{+}) = 0 \end{cases} {Mxy?(0??)=0Mz?(0??)=Mz0??90x??{Mxy?(0+?)=Mz0?Mz?(0+?)=0?
而在 t>0t>0t>0 的時的橫向和縱向磁化分量如下所示:
{Mx′y′(t)=Mz0e?t/T2Mz′(t)=Mz0(1?e?t/T1) \left\{ \begin{array}{l} M_{x'y'}(t) = M_{z}^{0} e^{-t/T_2} \\ M_{z'}(t) = M_z^0 (1 - e^{-t/T_1}) \end{array} \right. {Mxy?(t)=Mz0?e?t/T2?Mz?(t)=Mz0?(1?e?t/T1?)?
然后定義,當橫向磁化分量衰減為初始值Mx′y′(0+)M_{x'y'}(0_+)Mxy?(0+?)的37%的時刻就是 T2T_2T2? ,而當縱向磁化分量恢復為熱平衡態時的63%的時刻就是 T1T_1T1?: {Mx′y′(T2)≈37%Mx′y′(0+)Mz′(T1)≈63%Mz0\left\{\begin{array}{l}M_{x'y'}(T_2) \approx 37\% M_{x'y'}(0_+)\\M_{z'}(T_1) \approx 63\% M_z^0 \end{array}\right.{Mxy?(T2?)37%Mxy?(0+?)Mz?(T1?)63%Mz0??
生物組織典型值:

  • T1T_1T1?范圍:300–2000 ms
  • T2T_2T2?范圍:30–150 ms
  • 對一個給定的自旋系統, T1T_1T1? > T2T_2T2?

應用 (3.64)到(3.122) ,將旋轉坐標系下磁化矢量轉換到實驗室坐標系:
Mxy(t)=Mxy(0+)e?t/T2e?iω0tMz(t)=Mz0(1?e?t/T1)+Mz(0+)e?t/T1 \begin{align*} M_{xy}(t) &= M_{xy}(0_+) e^{-t/T_2} e^{-i\omega_0 t} \tag {3.124a} \\ M_z(t) &= M_z^0 (1 - e^{-t/T_1}) + M_z(0_+) e^{-t/T_1} \tag {3.124b} \end{align*} Mxy?(t)Mz?(t)?=Mxy?(0+?)e?t/T2?e?iω0?t=Mz0?(1?e?t/T1?)+Mz?(0+?)e?t/T1??(3.124a)(3.124b)?
其中
Mxy(0+)=Mx′y′(0+)e?iω0τp(3.125)M_{xy}(0_+) = M_{x'y'}(0_+) e^{-i\omega_0 \tau_p} \tag {3.125}Mxy?(0+?)=Mxy?(0+?)e?iω0?τp?(3.125)
式(3.125)這是在實驗室坐標系中由信號檢測系統觀測到的“初始”橫向磁化分量的強度。包含時間延遲項是因為我們在脈沖結束時(t=τpt=\tau_pt=τp??,其中τp\tau_pτp?為脈沖持續時間)將時間重置為零,以便描述弛豫效應。

式(3.124)給出了射頻脈沖后橫向與縱向磁化強度隨時間演化的"精確"唯象描述。具體而言,從方程(3.124a)可明確看出:在實驗室坐標系中,橫向磁化強度的演化特征表現為指數衰減e?t/T2e^{-t/T_{2}}e?t/T2?和繞B0B_{0}B0?場的進動e?iw0te^{-iw_{0}t}e?iw0?t。自由進動周期的時長取決于T2T_{2}T2?值。對于生物組織而言,T2T_{2}T2?值通常在數十毫秒量級,這使得在此期間能夠檢測到磁共振信號。

image.png

需要特別說明的是,當M?\vec{M}M因弛豫過程以螺旋運動方式”回歸“z軸時(如上圖所示),其模長并不守恒。該行為與激發期間M?\vec{M}M保持固定模長從z軸螺旋下行的運動特性存在本質區別。簡單的說,是因為:

  1. 射頻激發存在聚相的作用,這使得 z′z'z 軸宏觀磁化強分量 Mz′M_{z'}Mz? “全都”轉換為 x′y′x'y'xy 平面宏觀磁化強分量 Mx′y′M_{x'y'}Mxy?
  2. 弛豫過程中,聚相效果消失,因為失相,導致 x′y′x'y'xy 平面宏觀磁化強分量 Mx′y′M_{x'y'}Mxy? 快速衰減(相互抵消),但是這部分抵消并不能轉化為 z′z'z 軸的宏觀磁化強分量 Mz′M_{z'}Mz?
  3. Mz′M_{z'}Mz? 方向的恢復,實際上是高能級的原子躍遷回落到低能級

失相的唯象解釋(不是精確的量子力學角度的解釋):

  • 不同原子所處微觀磁場環境不同,使得不同的原子的進動頻率略有不同
  • 這使得原本相位同步的進動,開始慢慢分離,最后相互抵消,失相的過程請看下圖
    • 下圖是實驗室坐標下的 xyxyxy 平面,不同原子的磁矩 u?n\vec u_nun?xyxyxy 分量uxy,nu_{xy,n}uxy,n? (下圖黑線)和凈宏觀磁化矢量分量 MxyM_{xy}Mxy? (下圖紅線)
    • 下圖1,射頻激發剛結束,不同原子的磁矩 u?n\vec u_nun?xyxyxy 分量uxy,nu_{xy,n}uxy,n? 仍保持同相,其合成后的凈宏觀磁化矢量分量 MxyM_{xy}Mxy? 處于最大值,為Mxy(0+)M_{xy}(0_+)Mxy?(0+?)
    • 下圖2,不同原子的磁矩 u?n\vec u_nun?xyxyxy 分量uxy,nu_{xy,n}uxy,n? 因為進動頻率不同,在相同時間的旋轉角度不同,致使逐漸失去同相,他們的疊加后的MxyM_{xy}Mxy? 開始衰減
    • 下圖3,4,5,是一個失相逐漸擴大的橫向弛豫過程
    • 下圖6,弛豫完成,疊加后的凈宏觀磁化矢量分量MxyM_{xy}Mxy? 為0

image.png

本文來自互聯網用戶投稿,該文觀點僅代表作者本人,不代表本站立場。本站僅提供信息存儲空間服務,不擁有所有權,不承擔相關法律責任。
如若轉載,請注明出處:http://www.pswp.cn/pingmian/96015.shtml
繁體地址,請注明出處:http://hk.pswp.cn/pingmian/96015.shtml
英文地址,請注明出處:http://en.pswp.cn/pingmian/96015.shtml

如若內容造成侵權/違法違規/事實不符,請聯系多彩編程網進行投訴反饋email:809451989@qq.com,一經查實,立即刪除!

相關文章

ubuntu 20.04 安裝spark

安裝openjdk21 下載 wget https://download.java.net/openjdk/jdk21/ri/openjdk-2135_linux-x64_bin.tar.gz解壓 tar -xvf openjdk-2135_linux-x64_bin.tar.gzsudo mv jdk-21/ /opt/jdk-21/設置環境變量 echo export JAVA_HOME/opt/jdk-21 | sudo tee /etc/profile.d/java2…

第三方區塊鏈應用測評:【多簽錢包合約安全評估_閾值簽名機制與私鑰存儲安全性測試】

閾值簽名機制安全測試密碼學審計 采用門限簽名方案(TSS)的多簽錢包需驗證其閾值BLS簽名或ECDSA簽名算法的正確性。測試重點包括:分布式密鑰生成(DKG)過程的保密性(無密鑰信息泄露)、簽名碎片驗證…

大模型處理長文檔的挑戰和解決方案?

當前,AI 應用正處于極速發展階段,大語言模型(LLM)與檢索增強生成(RAG)系統已成為構建智能問答、知識管理等高階 AI 應用的核心引擎,被廣泛應用于金融分析、學術研究、企業合規等多個領域。然而&…

JavaWeb--day1--HTMLCSS

(以下內容全部來自上述課程及課件) web開發介紹 1. 什么是web? Web:全球廣域網,也稱為萬維網,能夠通過瀏覽器訪問的網站。 2. Web網站的工作流程 3. Web標準 Web標準也稱為網頁標準,由一系列的標準組成&#xf…

OpenEuler安裝gitlab,部署gitlab-runner

目錄 一、安裝gitlab 二、安裝部署docker設置源 三、下載部署runner ?編輯 四、構建CI/CD 一、安裝gitlab 1.查看OpenEuler版本 [rootlocalhost ~]# cat /etc/os-release NAME"openEuler" VERSION"24.03 (LTS-SP1)" ID"openEuler" VERSI…

實戰項目-----在圖片 hua.png 中,用紅色畫出花的外部輪廓,用綠色畫出其簡化輪廓(ε=周長×0.005),并在同一窗口顯示

實戰項目實現以下功能:對圖片 hua.png 進行輪廓提取,并在同一窗口中完成以下兩個繪制操作:用紅色畫出花的外部輪廓(即最外層輪廓)用綠色畫出該輪廓的近似多邊形,其中近似精度參數 ε 設置為輪廓周長的 0.00…

開源鴻蒙北向框架開發:系統服務理論詳解

系統服務的啟動 基本可以認為:OpenHarmony 的系統服務進程都是“由 init 直接或間接拉起”的。 直接方式: init 按 /system/etc/init/.cfg 啟動可執行(如 /system/bin/sa_main、/system/bin/samgr 等),這些進程的 PPid…

龍虎榜——20250909

上證指數今天縮量收陰線,跌破10日均線,目前日線總體處于高位寬幅震蕩中,小級別暫未明確方向,指數面臨方向選擇,需要注意高位股的風險。 深證指數今天縮量收陰線,跌破5日均線,接下來幾天方向的選…

基于dijkstra算法的WSN網絡MAC協議matlab仿真,分析網絡延遲與網絡開銷

目錄 1.程序功能描述 2.測試軟件版本以及運行結果展示 3.部分程序 4.算法理論概述 5.完整程序 1.程序功能描述 無線傳感器網絡(WSN, Wireless Sensor Network)是由大量低成本、低功耗的傳感器節點組成,通過無線通信實現數據采集、傳輸與…

ES數據庫啟動時可以設置1G一下內存嗎

可以,但強烈不建議在生產環境中這樣做。ES 7.0 版本確實可以設置小于 1GB 的堆內存,但這會帶來嚴重的性能問題和穩定性風險。 快速回答 # 最小化配置示例 - 僅用于測試或開發環境 export ES_JAVA_OPTS"-Xms512m -Xmx512m" ./bin/elasticsearch…

TI-92 Plus計算器:單位換算功能介紹

1 TI-92 Plus計算器:單位換算功能介紹 TI-92 Plus 內置了全面的單位換算功能,支持長度、質量、時間、溫度、面積、體積、速度等數十種物理量的單位轉換,操作直觀,無需手動輸入換算系數。以下是具體使用方法、示例及功能特點&#…

雪球科技Java開發工程師筆試題

單選 1.下列哪些語句關于內存回收的說明是正確的?( C ) A.內存回收程序允許程序員直接釋放內存 B.程序員必須創建一個線程來釋放內存 C.內存回收程序負責釋放無用內存 D.內存回收程序可以在指定的時間釋放內存對象 2.以下哪項不是Java基礎類型(A) A.String B.int C.b…

NV3041A-01芯片屏幕

1. 核心概覽這是一款集成了電源管理、顯示內存(RAM)、時序控制等多種功能的單片顯示驅動芯片(通常稱為Driver IC)。它采用COG(Chip-On-Glass) 工藝,直接將芯片綁定在玻璃基板上,使得…

aiagent知識點

一、MCP (Model Context Protocol) 1. 核心概念是什么:MCP是一個開放協議,用于在應用(如IDE、Agent) 和工具/數據源(如服務器、數據庫) 之間建立標準化的通信。目標:解決AI工具生態的碎片化問題…

第2節-過濾表中的行-WHERE

摘要:在本教程中,您將學習如何使用 PostgreSQL 的 WHERE 子句來篩選表中的行。 PostgreSQL WHERE 子句 SELECT FROM 語句從表中所有行的一個或多個列中查詢數據。實際上,你經常需要選擇滿足某個條件的行。 要根據條件從表中篩選行&#xf…

IACheck賦能AI環評報告審核,推動環保設備制造行業發展

在“雙碳目標”和綠色制造的背景下,環保設備制造行業正在迎來快速發展。然而,環評報告作為項目合規的“通行證”,卻一直是企業最頭疼的環節之一:編寫復雜、審核周期長、錯誤率高。傳統的審核模式不僅耗時耗力,還容易出…

常見的多態

一、核心概念多態(Polymorphism) 的字面意思是“多種形態”。在Java中,它指的是:同一個行為(方法)具有多個不同表現形式或形態的能力。更具體地說,它允許你:父類的引用指向子類的對象…

20250909_排查10.1.1.190檔案庫房綜合管理系統20250908備份缺失問題+優化scp腳本(把失敗原因記錄進日志)并測試腳本執行情況

目錄 一、現象(圖中現象) 二、可能原因 1、本地源文件沒生成(最可能原因) 2、清理邏輯誤刪 三、時間同步檢查 1、運行命令查看當前時間源:w32tm /query /source 2、修改為內網 NTP 服務器 10.1.10.251: 四、腳本優化 1、scp.bat vs roboscp.bat 對比表 2、robo…

OpenHarmony之有源NFC-connected_nfc_tag模塊詳解

1. 模塊概述 Connected NFC Tag模塊是OpenHarmony NFC子系統的一部分,專門用于處理有源NFC標簽的讀寫操作。與基礎NFC模塊不同,該模塊專注于與已連接的有源NFC標簽進行交互,提供更高層次的API供應用層使用。 設備必須具備有源NFC標簽芯片,才能使用有源NFC標簽的讀和寫服務…

MySql案例詳解之事務

下面我會從“事務是什么”→“為什么需要事務”→“事務的四大特性(ACID)”→“MySQL中怎么用事務”→“常見坑與調試技巧”→“完整實戰案例(含代碼、輸出、講解)”六個層次,給你一個“看完就能上手”的MySQL事務速查…