吳鑫垚 李發(fā)琪 李成海

基金項目：重慶市英才計劃（cstc2021ycjh-bgzxm0068）；重慶市自然科學基金創(chuàng)新發(fā)展聯(lián)合基金項目（CSTB2022NSCQ-LZX0028）；國家藥監(jiān)局超聲手術(shù)設備質(zhì)量評價重點實驗室開放課題（SMDTKL-2023-2-01）

第一作者簡介：吳鑫垚（1998-），男，碩士研究生。研究方向為生物醫(yī)學超聲。

*通信作者：李成海（1987-），男，博士，副研究員。研究方向為生物醫(yī)學超聲。

DOI：10.19981/j.CN23-1581/G3.2024.13.039

摘? 要：熱療在各種良惡性腫瘤以及非腫瘤疾病的治療中具有重要臨床應用價值，其原理為將物理能量轉(zhuǎn)化為靶區(qū)的熱能并體現(xiàn)為溫度升高進而實現(xiàn)疾病治療的目的。為保證熱療的安全性和有效性，治療中的無創(chuàng)溫度監(jiān)控一直是該領域的核心關鍵技術(shù)。該文對熱療中所使用到的核磁共振（Magnetic Resonance，MR）測溫、超聲（Ultrasound，US）測溫以及其他無創(chuàng)測溫技術(shù)及其相關研究進展進行綜述，旨在為無創(chuàng)測溫技術(shù)應用于臨床熱療監(jiān)控提供參考。

關鍵詞：無創(chuàng)測溫；熱療；核磁共振測溫；超聲測溫；紅外測溫；微波測溫

中圖分類號：TB559? ? ? 文獻標志碼：A? ? ? ? ? 文章編號：2095-2945（2024）13-0161-06

Abstract： Thermotherapy has important clinical applications in the treatment of various benign and malignant tumors as well as non-tumor diseases， and its principle is to convert physical energy into thermal energy in the target area and manifest it as a temperature increase to achieve the purpose of disease treatment. In order to ensure the safety and effectiveness of thermotherapy， noninvasive temperature monitoring has been the core technology in this field. In this paper， we summarize Magnetic Resonance （MR） thermometry， Ultrasound （US） thermometry， and other non-invasive thermometry techniques used in thermotherapy and their related research progress， aiming to provide reference for the application of non-invasive thermometry techniques in clinical thermotherapy monitoring.

Keywords： noninvasive thermometry; thermotherapy; MR thermometry; US thermometry; infrared temperature measurement; microwave thermometry

熱療在各種良惡性腫瘤以及非腫瘤疾病的治療中具有重要臨床應用價值，目前常用的治療方式主要有超聲、微波、射頻和紅外等，將物理能量轉(zhuǎn)化為靶區(qū)的熱能并體現(xiàn)為溫度的升高，產(chǎn)生熱療的治療效果。熱療過程中，為了保證治療的安全性和有效性，溫度的監(jiān)控是熱療必不可少的關鍵技術(shù)。

考慮到生物組織的復雜性，熱療的溫度監(jiān)控一直是一個令人關注的問題。熱療是熱能進入或者排出體內(nèi)的一種治療方式。熱療將熱源溫度控制在不損害正常組織范圍內(nèi)而獲得有效的治療結(jié)果，但是如果溫度超過閾值就會產(chǎn)生細胞不可逆的損傷。在臨床治療中需要對熱療進行監(jiān)控，防止過高的溫度累積造成正常細胞的凋亡?，F(xiàn)有的測溫技術(shù)可以分為有創(chuàng)和無創(chuàng)測溫，但是為了防止造成二次傷害，貼合熱療微無創(chuàng)的治療理念，熱療中的溫度監(jiān)控主要基于無創(chuàng)測溫技術(shù)。目前的無創(chuàng)測溫技術(shù)主要有核磁共振（Magnetic Resonance，MR）測溫[1]和超聲（Ultrasound，US）測溫[2]，紅外測溫和微波測溫也展示了一定的應用前景。本文對熱療中所使用到的無創(chuàng)測溫技術(shù)及其相關研究進展進行綜述，旨在為無創(chuàng)測溫技術(shù)應用于臨床熱療監(jiān)控提供參考。

1? MR測溫

MR測溫技術(shù)是在MR成像基礎上，利用不同MR成像參數(shù)，包括擴散系數(shù)[3]、質(zhì)子共振頻率漂移（Proton Resonance Frequency shift，PRF）[4]、縱向弛豫時間（T1）[5]和橫向弛豫時間（T2）[6]等，在多個平面或體積上進行溫度成像。

1.1? 基于PRF的MR測溫

1995年，在前人提出相關理念的基礎上，Poorter等[7]首次基于PRF測溫開展了人體小腿肌肉的溫度檢測實驗。PRF測溫方法的原理主要基于溫度變化？駐T與相位差？駐？覫有關

？駐T=-■，

式中：Ψ為旋磁比，α為PRF的位移系數(shù)（-0.01 ppm·℃-1），B0為磁場，TE為回波時間[8]?；赑RF測溫的優(yōu)勢在于在各個組織中有均勻且相對較高的溫度敏感性，常用于臨床治療的溫度監(jiān)測。Jeanmonod等[9]已經(jīng)應用新的經(jīng)顱磁共振成像引導聚焦超聲技術(shù)進行非侵入性中央外側(cè)丘腦切除術(shù)作為慢性神經(jīng)病理性疼痛的治療。臨床驗證說明了MR測溫監(jiān)測熱療的安全性和有效性。雖然PRF已經(jīng)成熟應用于熱療的監(jiān)控，但是在脂肪中的效果不佳，原因是脂肪中缺乏氫鍵，導致PRF溫度測量造成誤差。為了解決這一問題，Zong等[10]提出了一種通過估計溫度不敏感的脂肪向量來消除脂肪的相位貢獻的方法，結(jié)果發(fā)現(xiàn)該方法用于監(jiān)控聚焦超聲消融離體脂肪組織時，與傳統(tǒng)的PRF方法之間的溫度差異高達28%。臨床上患者呼吸和其他身體運動或組織變化引起的磁場漂移和磁化率變化所帶來的誤差是影響PRF測溫準確性的另一個關鍵因素。Hensen等[11]在微波消融的條件下，建立了一種2D-PRF測溫法，創(chuàng)建了二進制掩碼S（x，y，t），可以用于矯正磁化率從而提高患者組織位置變化引起的測溫誤差。實現(xiàn)熱療中的實時溫度監(jiān)測，MR測溫需要更高的時間分辨率和空間分辨率。為了提高MR測溫的速度，Peng等[12]基于射頻擾相梯度回波序列和同步多層圖像成像序列在高強度聚焦超聲治療中進行快速的MR測溫實驗，在離體實驗中驗證了該序列可以提高時間的分辨率以及成像的覆蓋范圍。

1.2? 基于縱向弛豫時間T1的MR測溫

MR測溫的另一種重要測溫技術(shù)為基于縱向弛豫時間T1的測溫方式。1983年，Parker等[13]發(fā)現(xiàn)二維核磁共振中T1與溫度呈線性變化。基于T1對溫度的依賴性，McLean等[14]建立并評價了一種基于T1的校正方法，對乳腺脂肪組織加熱引起的非局域場效應引起的PRF測溫誤差進行校正，減少了由于脂肪組織加熱引起的非局部效應造成的MR測溫誤差。為了將T1在脂肪中良好的測溫效果與PRF進行結(jié)合，Zhang等[15]開發(fā)了可變翻轉(zhuǎn)角黃金角有序3D徑向堆疊MR技術(shù)，分別在水組織和脂肪組織中同時進行PRF和基于T1的測溫。在離體高強度聚焦超聲消融過程中，PRF和T1測量的溫度絕對平均差異在2 ℃以內(nèi)，可以獲得相較于單PRF測溫更準確的溫度圖像。

MR技術(shù)有著非侵入性和非電離的特點，并且對于臨床治療干擾較小，測溫精確度較高，已經(jīng)成為臨床測溫的金標準，在聚焦超聲、激光熱療等臨床治療中得到一定應用[16]。但是MR設備昂貴，在臨床使用中對于醫(yī)院的設備環(huán)境要求較為苛刻，是限制其普及的主要原因。

2? US測溫

超聲測溫是一種低成本的測溫方法，其中最為常用的方法就是利用超聲波與溫度相關的聲學參數(shù)以及超聲成像信號的改變獲取目標區(qū)域溫度的信息。其中基于組織聲學特性變化測溫主要有聲速、聲衰減系數(shù)、非線性參數(shù)B/A等，基于超聲成像的測溫主要包括超聲圖像灰度變化、Nakagami成像、光聲成像、剪切波成像以及人工智能技術(shù)結(jié)合射頻信號預測組織溫度等。

2.1? 基于組織聲學特性變化的US測溫

聲速作為典型的聲學參數(shù)，在溫度改變時也會隨之變化。Engrand等[17]實驗得出了在10～30 ℃范圍內(nèi)縱波速度與心肌溫度之間的實驗線性關系。當組織中的溫度變化，使聲速局部產(chǎn)生變化，傳播介質(zhì)中發(fā)生熱膨脹從而觀察到回波時移。軸向溫升？駐T（z）與軸向的位移相關

？駐T（z）=■·■，

式中：t（z）為在組織深度z處的估計回波偏移，c0為環(huán)境下的聲速，α為熱膨脹系數(shù)，β為隨溫度改變的聲速[18]。聲波回波時移可以作為溫度判斷的標準，F(xiàn)an等[19]使用了改進的動態(tài)幀選擇算法和改進的自適應濾波方法與熱膨脹模型相結(jié)合，比較US的回波時移和熱電偶測得的溫度，確定了恒定溫度評價系數(shù)。當組織比基礎溫度高出6.4、9.8和19.3 ℃時，均方根評估誤差分別約為0.3、0.5和0.8 ℃。研究發(fā)現(xiàn)，溫度的改變造成的回波的偏移除了體現(xiàn)在時域上，也會體現(xiàn)在頻域上。Amini等[20]使用高分辨率譜估計方法來跟蹤與溫度變化相關的多個諧波頻率的頻移，發(fā)現(xiàn)基頻變化與溫度變化之間存在著線性關系。

熱應變成像（Thermal Strain Imaging，TSI）是利用熱膨脹和聲速的變化進行溫度成像的典型應用。Foiret等[21]在TSI基礎上，提出了一種運動補償方法，使用超聲測溫作為反饋，對存在運動的熱療進行了控制測試，實驗證明溫度可以保持在要求值的0.3 ℃以內(nèi)。為了減少TSI對于采集設備的性能要求，Yin等[22]驗證基于“無限小回聲應變?yōu)V波器”模型的內(nèi)插IQ圖像進行TSI的可行性，減少熱應變成像的誤差。他們在呼吸分離TSI的基礎上，提出了多線程TSI方法，在多線程中對每個線程進行單獨的計算處理，獲得的TSI結(jié)果進行平均以獲得合并的輸出。

衰減系數(shù)的改變與溫度的改變也有著直接的關聯(lián)。背向散射能量（Changes in Backscattered Energy，CBE）是在衰減系數(shù)基礎上提出的超聲測溫方式

CBE=■·■·■，

式中：TR為參考溫度，α（T）、η（T）和x分別為組織的衰減系數(shù)、反向散射系數(shù)和路徑長度[23]。US二維測溫可以通過CBE呈現(xiàn)，Shaswary等[24]采用基于二次諧波CBE的US測溫方法生成二維溫度圖，檢測低強度聚焦超聲產(chǎn)生的局部加熱區(qū)。結(jié)果表明，在37 ℃（基線）到47℃的溫度范圍內(nèi)，聲學諧波的CBE與病灶組織溫度之間存在直接相關性。說明基于組織衰減系數(shù)的變化也可反映組織中溫度的變化。

聲波的傳播過程是非線性的過程，溫度改變會使高次諧波成分發(fā)生改變。為了證明非線性參數(shù)B/A與溫度的關聯(lián)性和測溫的可行性，Van等[25]用水為實驗對象，對比了B/A以及16個模擬測量的組合。研究證實了B/A作為非侵入性測溫基礎的可行性，并發(fā)現(xiàn)在噪聲水平不超過40 dB的情況下，B/A對溫度變化比其他實際聲學參數(shù)更敏感。為了在組織實驗中證明非線性參數(shù)B/A與溫度的關聯(lián)，Dong等[26]進行了高強度聚焦超聲輻照離體豬肌肉和肝實驗，研究了非線性諧波幅值與溫度的關系，發(fā)現(xiàn)基波、二次諧波和三次諧波隨溫度的幅度變化與溫度變化密切相關。

2.2? 基于超聲成像的測溫

考慮到組織聲學特性對溫度具有依賴性，憑借不同條件下組織聲學特性的測試可以間接實現(xiàn)組織溫度的測量，但除了溫度之外影響組織聲學特性的因素仍然較多，獲得單一的組織聲參數(shù)——溫度關系仍然較為困難。如何利用現(xiàn)有較為成熟的超聲成像手段實現(xiàn)組織溫度測量廣受研究者的關注，其中包括B超灰度變化、Nakagami成像、光聲成像和剪切波等。

B超圖像的灰度變化可以反映溫度的變化，Teixeira等[27]評估了超聲圖像的平均灰度級在非侵入性地估計溫度在時間和空間上變化的潛力，在2個均勻加熱/冷卻循環(huán)中，兩組實驗的平均溫差分別為0.19 ℃和0.06 ℃。說明通過研究灰度值的變化可以在一定程度上預測溫度的變化量。

考慮到基于超聲的包絡信號處理得到的Nakagami成像的參數(shù)在溫度改變時會相應發(fā)生改變，Zhang等[28]計算了熱病變中超聲波Nakagami參數(shù)的動態(tài)變化，同時重建了超聲波 B 型和Nakagami圖像。其研究發(fā)現(xiàn)在熱消融過程中，Nakagami參數(shù)的平均值呈上升趨勢，體外模型從0.72增加到1.01，體內(nèi)模型從0.54增加到0.72。進一步的研究發(fā)現(xiàn)，由于溫度升高而引起的Nakgami形狀參數(shù)m的變化隨著初始m值的增加而增加，經(jīng)過m值的修正成像，可以清楚地顯示軟組織內(nèi)的溫度分布[29]。

光聲成像（Photoacoustic imaging，PA）可以通過跟蹤溫度引起的光聲信號幅度的變化來測量溫度。2005年，Larina等[30]發(fā)現(xiàn)光聲測量的組織溫度與實際組織溫度之間有很好的一致性，溫度監(jiān)測精度優(yōu)于1 ℃，空間分辨率約為1 mm。PA與其他監(jiān)測方式結(jié)合進行多模態(tài)的溫度監(jiān)測。Zhang等[31]采用順序傳輸、PA成像、PA測溫和US成像等技術(shù)，顯示雙模式圖像并記錄局部溫度變化。實驗發(fā)現(xiàn)，對于皮膚等軟組織，如果局部溫度不低于50 ℃，則達到PA幅度與測量溫度之間的線性比例。如果局部組織溫度超過臨界點，生物組織的吸收光譜可能會因蛋白質(zhì)變性而發(fā)生變化。

剪切波的變化與組織的硬度有關，當溫度改變時，組織的剪切波也會產(chǎn)生相應變化?；诩羟胁▽τ跍囟鹊拿舾行裕珹rnal等[32]發(fā)現(xiàn)剪切波成像對局部組織彈性的估計可以反映HIFU治療期間溫度變化，在溫度小于45 ℃時，基于超聲的溫度估算與組織硬度變化高度相關（r2=0.91～0.97）。剪切波在不同溫度區(qū)間的變化是不同的，研究發(fā)現(xiàn)在55～60 ℃范圍內(nèi)，剪切波速對溫度的敏感性很低，在所有情況下，速度均小于5m·s-1，溫度大于55～60 ℃時，速度呈陡增趨勢[33]。

超聲射頻信號中包含的信息復雜繁多，經(jīng)過這些年人工智能技術(shù)的發(fā)展，可以更有效地發(fā)掘信號中深層特征，進行溫度的預測研究。2006年，Teixeira等[34]就提出了從采集的射頻線中提取的光譜特征和時間特征，利用徑向基函數(shù)神經(jīng)網(wǎng)絡進行多點溫度的預測。最佳擬合的預測模型在預測范圍內(nèi)的最大絕對誤差小于0.4 ℃。2008年，Teixeira等[35]又做出新的突破，對凝膠體模采集的超聲回波信號的時移進行跟蹤，并將過去的溫度值作為徑向基函數(shù)神經(jīng)網(wǎng)絡的輸入信息。使用多目標遺傳算法選擇最好的神經(jīng)模型，預測結(jié)果的最大絕對誤差低于0.5 ℃，標準偏差為0.11 ℃。近些年來，隨著神經(jīng)網(wǎng)絡模型的性能提升，可以考慮更多信號的深層特征。Chen等[36]接收探頭中的超聲元件256個HIFU元件發(fā)出的超聲脈沖，使用卷積長短期記憶神經(jīng)網(wǎng)絡從超聲通道數(shù)據(jù)中生成溫度圖像。軸面平均差為0.57±0.33 ℃，最大差值為1.99±1.07 ℃。在冠狀平面上，平均差為0.33±0.19 ℃，最大差為1.54±1.04 ℃。Byra等[37]使用連體網(wǎng)絡對射頻數(shù)據(jù)的小型二維斑塊進行局部比較，得出相似性分數(shù)，從而評估組織溫度。組織樣本的結(jié)果顯示，CBE和深度學習技術(shù)的r2分別為0.95和0.90。

由于US的低成本特點，US測溫是熱療溫度監(jiān)控的理想方式。但是現(xiàn)有的研究都只能在較低溫度（小于50 ℃）下有較為精確的準確度，一旦超過閾值，就會出現(xiàn)較大的溫差，這也限制了US測溫在臨床的使用。

3? 其他測溫方式

紅外測溫是一種非破壞性、非接觸式溫度測量技術(shù)，可實現(xiàn)實時、全場成像。Geng等[38]在射頻消融實驗中，超聲回波時移法測得的溫度分布與紅外圖像測得的溫度分布一致，與紅外圖像相比，當溫度超過約45 ℃時，超聲波檢測到的溫度分布區(qū)域往往被低估。這一實驗展現(xiàn)了紅外測溫的潛能。

微波測溫可以通過測量物體或物體發(fā)出的自然熱電磁輻射的變化來表征內(nèi)部溫度，是一種非電離、非侵入性、廉價、快速和被動監(jiān)測的測溫方式。Villa等[39]利用偽相關輻射計反演物體或生物組織溫度的新方案，對多頻率輻射計的性能進行了評估，結(jié)果顯示在檢索未知溫度方面效果顯著。獲得相對于由噪聲源提供的溫度的誤差小于0.1%。

4? 結(jié)束語

本文總結(jié)了熱療中無創(chuàng)測溫的相關方法及研究進展。目前MR和US測溫是臨床上普遍看好的無創(chuàng)測溫方式，相比于成熟的MRI測溫，US測溫具有便攜、低成本的特點。盡管US測溫在50 ℃以下已經(jīng)取得令人矚目的研究成果，但缺乏相關研究證明50 ℃以上存在可靠的US測溫技術(shù)。US射頻回波信號蘊含很多與溫度相關特征，隨著人工智能技術(shù)的發(fā)展，射頻數(shù)據(jù)結(jié)合人工智能技術(shù)將讓US測溫進一步突破成為了可能。而不同測溫方式的結(jié)合也是未來US測溫的研究方向，通過針對不同溫升選擇合適的測溫方式可以提升測溫的準確性、實時性以及抗干擾能力。總而言之，熱療測溫未來將向著智能化、高實時性、高準確性和低成本的方向發(fā)展。

參考文獻：

[1] DESCLIDES M， OZENNE V， BOUR P， et al. Real-time automatic temperature regulation during in vivo MRI-guided laser-induced thermotherapy（MR-LITT）[J]. Scientific Reports， 2023，13（1）：3279.

[2] PENG K， ZHOU P， LIU W. Long-term efficacy of ultrasound-guided percutaneous laser ablation for low-risk papillary thyroid microcarcinoma： A 5-year follow-up study[J]. BioMed Research International， 2021.

[3] SPARACIA G， CANNELLA R， LO R V， et al. Brain-core temperature of patients before and after orthotopic liver transplantation assessed by DWI thermometry[J]. Japanese Journal of Radiology， 2018（36）：324-330.

[4] GUO R， WANG X， LIU J， et al. Experimental study on the accuracy of non-invasive temperature measurement by magnetic resonance in microwave ablation[J]. Journal of Cancer Research and Therapeutics， 2022，18（5）：1256-1260.

[5] QIAN E， POOJAR P， FUNG M， et al. Magnetic Resonance Fingerprinting based Thermometry （MRFT）： application to ex vivo imaging near DBS leads[J]. Physics in Medicine & Biology， 2023，68（17）：17NT01.

[6] BARON P， RIES M， DECKERS R， et al. In vivo T2-based MR thermometry in adipose tissue layers for high-intensity focused ultrasound near-field monitoring[J]. Magnetic resonance in medicine， 2014， 72（4）： 1057-1064.

[7] POORTER J D， WAGTER C D， DEENE Y D， et al. Noninvasive MRI thermometry with the proton resonance frequency （PRF） method： in vivo results in human muscle[J]. Magnetic resonance in medicine， 1995，33（1）：74-81.

[8] BING C， STARUCH R M， TILLANDER M， et al. Drift correction for accurate PRF-shift MR thermometry during mild hyperthermia treatments with MR-HIFU[J]. International Journal of Hyperthermia， 2016，32（6）：673-687.

[9] JEANMONOD D， WERNER B， MOREL A， et al. Transcranial magnetic resonance imaging-guided focused ultrasound： noninvasive central lateral thalamotomy for chronic neuropathic pain[J]. Neurosurgical focus， 2012，32（1）：E1.

[10] ZONG S， SHEN G， MEI C S. Proton resonance frequency-based thermometry for aqueous and adipose tissues[J]. Medical Physics， 2021， 48（10）： 5651-5660.

[11] HENSEN B， HELLMS S， WERLEIN C， et al. Correction of heat-induced susceptibility changes in respiratory-triggered 2D-PRF-based thermometry for monitoring of magnetic resonance-guided hepatic microwave ablation in a human-like in vivo porcine model[J]. International Journal of Hyperthermia， 2022，39（1）：1387-1396.

[12] PENG Y， ZOU C， QIAO Y， et al. Fast MR thermometry using an echo-shifted sequence with simultaneous multi-slice imaging[J]. Magnetic Resonance Materials in Physics， Biology and Medicine， 2018（31）：771-779.

[13] PARKER D L， SMITH V， SHELDON P， et al. Temperature distribution measurements in two-dimensional NMR imaging[J]. Medical physics， 1983， 10（3）： 321-325.

[14] MCLEAN M K， PARKER D L， ODON P， et al. A T1-based correction method for proton resonance frequency shift thermometry in breast tissue[J]. Medical physics， 2021，48（9）：4719-4729.

[15] ZHANG L， ARMSTRONG T， Li X， et al. A variable flip angle goldenethod for proton resonance frequency shiftnique for simultaneous proton resonant frequency shift and T1-based thermometry[J]. Magnetic resonance in medicine， 2019，82（6）：2062-2076.

[16] BLACKWELL J， KRA■NY M J， O'BRIEN A， et al. Proton resonance frequency shift thermometry： a review of modern clinical practices[J]. Journal of Magnetic Resonance Imaging， 2022，55（2）：389-403.

[17] ENGRAND C， LAUX D， FERRANDIS J Y， et al. Velocimetric ultrasound thermometry applied to myocardium protection monitoring[J]. Ultrasonics， 2018（87）：1-6.

[18] MAASSND C， L R， DAMIANOU C A. Noninvasive temperature estimation in tissue via ultrasound echo-shifts. Part I. Analytical model[J]. The Journal of the Acoustical Society of America， 1996，100（4）：2514-2521.

[19] FAN P， YIN C， XUE H， et al. In vivo evaluation of two-dimensional temperature variation in perirenal fat of pigs with B-mode ultrasound[J]. Journal of Applied Physics， 2019，126（8）.

[20] AMINI A N， EBBINI E S， GEORGIOU T T. Noninvasive estimation of tissue temperature via high-resolution spectral analysis techniques[J]. IEEE Transactions on Biomedical Engineering， 2005，52（2）：221-228.

[21] FOIRET J， FERRARA K W. Spatial and temporal control of hyperthermia using real time ultrasonic thermal strain imaging with motion compensation， phantom study[J]. PLoS One， 2015，10（8）：e0134938.

[22] YIN C， WANG G， YANG K， et al. Thermal strain imaging in vivo using interpolated IQ-images[J]. Ultrasonics， 2021（110）：106292.

[23] ZHOU Y. Noninvasive thermometry in high-intensity focused ultrasound ablation[J]. Ultrasound quarterly， 2017，33（4）：253-260.

[24] SHASWARY E， ASSI H， YANG C， et al. Noninvasive calibrated tissue temperature estimation using backscattered energy of acoustic harmonics[J]. Ultrasonics， 2021（114）：106406.

[25] VAN D K W A， VERWEIJ M D. A feasibility study for non-invasive thermometry using non-linear ultrasound[J]. International Journal of Hyperthermia， 2011，27（6）：612-624.

[26] DONG H， QIAN S Y， TAN Q， et al. Temperature measurement using passive harmonics during high intensity focused ultrasound exposures in porcine tissue[J]. Acta Phys. Pol. A， 2018（134）：524-528.

[27] TEIXEIRA C A， ALVARENGA A V， CORTELA G， et al. Feasibility of non-invasive temperature estimation by the assessment of the average gray-level content of B-mode images[J]. Ultrasonics， 2014，54（6）：1692-1702.

[28] ZHANG S， SHANG S， HAN Y， et al. Ex vivo and in vivo monitoring and characterization of thermal lesions by high-intensity focused ultrasound and microwave ablation using ultrasonic nakagami imaging[J]. IEEE transactions on medical imaging， 2018，37（7）：1701-1710.

[29] TAKEUCHI M， MATSUI Y， SAKAI T， et al. Investigation of initial value dependence in the statistical analysis of ultrasonic scattered echoes for the non-invasive estimation of temperature distribution in biological tissue[J]. Japanese Journal of Applied Physics， 2019，58（SG）：SGGE09.

[30] LARINA I V， LARIN K V， ESENALIEV R O. Real-time optoacoustic monitoring of temperature in tissues[J]. Journal of Physics D： Applied Physics， 2005，38（15）：2633.

[31] ZHANG Y， WANG L. Array-based high-intensity focused ultrasound therapy system integrated with real-time ultrasound and photoacoustic imaging[J]. Biomedical Optics Express， 2023，14（3）：1137-1145.

[32] ARNAL B， PERNOT M， TANTER M. Monitoring of thermal therapy based on shear modulus changes： I. shear wave thermometry[J]. IEEE transactions on ultrasonics， ferroelectrics， and frequency control， 2011，58（2）：369-378.

[33] GIURAZZA F， MASSARONI C， SILVESTRI S， et al. Preliminary analysis of ultrasound elastography imaging-based thermometry on non-perfused ex vivo swine liver[J]. Journal of Ultrasound， 2020（23）：69-75.

[34] TEIXEIRA C A， RUANO A E， RUANO M G， et al. Non-invasive temperature prediction of in vitro therapeutic ultrasound signals using neural networks[J]. Medical and Biological Engineering and Computing， 2006（44）：111-116.

[35] TEIXEIRA C A， RUANO M G， RUANO A E， et al. A soft-computing methodology for noninvasive time-spatial temperature estimation[J]. IEEE Transactions on Biomedical Engineering， 2008，55（2）：572-580.

[36] CHEN C， KIM Y， AUDIGIER C， et al. Ultrasound thermometry using an ultrasound element and deep learning for HIFU[C]//2019 IEEE International Ultrasonics Symposium （IUS）. IEEE， 2019：548-551.

[37] BYRA M， KLIMONDA Z， KRUGLENKO E， et al. Unsupervised deep learning based approach to temperature monitoring in focused ultrasound treatment[J]. Ultrasonics， 2022（122）：106689.

[38] GENG X， ZHOU Z， LI Q， et al. Comparison of ultrasound temperature imaging with infrared thermometry during radio frequency ablation[J]. Japanese Journal of Applied Physics， 2014，53（4）：047001.

[39] VILLA E， AJA B， DE LA FUENTE L， et al. Multifrequency Microwave Radiometry for Characterizing the Internal Temperature of Biological Tissues[J]. Biosensors， 2022，13（1）：25.