生物傳熱方程中灌注率函數(shù)的數(shù)值反演算法

曹慶發(fā)，胡 彬，萬 殊，王澤文

生物傳熱方程中灌注率函數(shù)的數(shù)值反演算法

曹慶發(fā)，胡 彬，萬 殊，*王澤文

（東華理工大學理學院，江西，南昌 330013）

本文研究了一類生物傳熱方程的灌注率函數(shù)反演問題?；诟郊拥姆蔷植織l件和有限差分的Crank-Nicolson方法，構(gòu)造了重建灌注率函數(shù)的迭代算法；經(jīng)進一步簡化后，得到了反演灌注率的一個顯格式。為克服計算的不穩(wěn)定性，引入移動平均濾波方法對誤差數(shù)據(jù)進行去噪，算例結(jié)果表明結(jié)合移動平均濾波去噪的數(shù)值反演算法是可行的，能有效反演出灌注率函數(shù)。

生物傳熱方程，灌注率，反問題，有限差分，移動平均

0 引言

生物醫(yī)學傳熱研究不僅在理論上很重要，而且具有重要的實際應(yīng)用價值[1]。例如，生物傳熱方程已被應(yīng)用于模擬高熱、血栓形成和血管硬化等研究中[2]。科學家提出了若干不同生物組織內(nèi)的傳熱模型，例如最常用的是由Pennes提出的生物傳熱模型[3]，以及文獻[4-6]提出連續(xù)型生物傳熱模型。

由Pennes提出的生物傳熱方程[3]為

它反映了血液灌注率，故稱其為血液灌注率函數(shù)。

和附加的非局部條件

不同于文獻[11-12]，本文受文獻[16]中研究的啟發(fā)，提出一種有限差分的反演算法，該方法無需事先將反問題轉(zhuǎn)化為源項反演。本文接下來安排如下：第二小節(jié)基于Crank-Nicolson格式給出兩種有限差分的數(shù)值反演算法；第三小節(jié)給出反問題的數(shù)值算例。

2 反問題的有限差分解法

2.1 反問題的數(shù)值解法

利用有限差分的Crank-Nicolson方法，將方程（2）離散為

由邊界條件（3）-（5），有

對于附加條件（6），利用數(shù)值積分的復(fù)化梯形公式得

將（13）-（14）改寫成矩陣形式為：

Step 3. 計算

上述差分格式的矩陣形式為

2.2 噪聲處理的移動平均濾波

3 數(shù)值算例

算例1 考慮生物傳熱反問題：

其中精確解為

圖1 算例1的反演結(jié)果對比

算例2 考慮生物傳熱反問題：

其中精確解為

圖2 算例2的反演結(jié)果對比

數(shù)值算例的結(jié)果圖1，圖2表明所給出的算法是可行的，且迭代算法的數(shù)值反演效果更佳，特別是在第二個數(shù)值算例中算法體現(xiàn)了較強的抗噪能力，這可能是移動平均濾波對算例2的數(shù)據(jù)去噪效果更佳的緣故。該方法也可以推廣到求解高維生物傳熱方程的相關(guān)反問題。

[1] Chato J C. Fundamentals of Bioheat Transfer[M]. Berlin: Springer-Verlag, 1990.

[2] Liu J, Xu L X. Boundary information based diagnostics on the thermal states of biological bodies[J]. International Journal of Heat and Mass Transfer, 2000,43(16): 2827-2839.

[3] Pennes H H. Analysis of tissue and arterial blood temperatures in the resting human forearm[J]. Journal of applied physiology, 1948, 1(2): 93-122.

[4] Wulff W. The energy conservation equation for living tissue[J]. IEEE transactions on biomedical engineering, 1974 (6): 494-495.

[5] Klinger H G. Heat transfer in perfused biological tissue—I: General theory[J]. Bulletin of Mathematical Biology, 1974, 36: 403-415.

[6] Chen M M, Holmes K R. Microvascular contributions in tissue heat transfer[J]. Annals of the New York Academy of Sciences, 1980, 335(1): 137-150.

[7] Grabski J K, Lesnic D, Johansson B T. Identification of a time-dependent bio-heat blood perfusion coefficient[J]. International Communications in Heat and Mass Transfer, 2016, 75: 218-222.

[8] Lin Y. An inverse problem for a class of quasilinear parabolic equations[J]. SIAM journal on mathematical analysis, 1991, 22(1): 146-156.

[9] Cannon J R, Lin Y, Wang S. Determination of a control parameter in a parabolic partial differential equation[J]. Journal of the Australian Mathematical Society, Series B Applied mathematics, 1991, 33: 149-163.

[10] Kerimov N B, Ismailov M I. An inverse coefficient problem for the heat equation in the case of nonlocal boundary conditions[J]. Journal of Mathematical Analysis & Applications, 2012, 396(2): 546-554.

[11] Wang S, Lin Y. A finite-difference solution to an inverse problem for determining a control function in a parabolic partial differential equation[J]. Inverse Problems, 1989, 5(4): 631-640.

[12] Ismailov M I, Kanca F. An inverse coefficient problem for a parabolic equation in the case of nonlocal boundary and overdetermination conditions[J]. Mathematical Methods in the Applied Sciences, 2011, 34(6): 692-702.

[13] Ismailov M I, Bazán F S V, Bedin L. Time-dependent perfusion coefficient estimation in a bioheat transfer problem[J]. Computer Physics Communications, 2018, 230: 50-58.

[14] Hazanee A, Lesnic D. Determination of a time-dependent coefficient in the bioheat equation[J].International Journal of Mechanical Sciences, 2014, 88: 259-266.

[15] Trucu D, Ingham D, Lesnic D. An inverse coefficient identification problem for the bio-heat equation[J]. Inverse Problems in Science and Engineering, 2009, 17(1): 65-83.

[16] Wang Z, Ruan Z, Huang H, et al. Determination of an unknown time-dependent heat source from A nonlocal measurement by finite difference method[J]. Acta Mathematicae Applicatae Sinica, English Series, 2020, 36(1): 151-165.

[17] 黃何露,王澤文,阮周生,等. 一類擴散方程尋源反問題的有限差分法[J].贛南師范大學學報,2018,39(3):20-23.

[18] 邱淑芳,王澤文,曾祥龍,等. 一類時間分數(shù)階擴散方程中的源項反演解法[J]. 江西師范大學學報:自然科學版, 2018, 42(6): 610-615.

NUMERICAL METHOD FOR RECOVERING PERFUSION COEFFICIENT IN A BIOLOGICAL HEAT TRANSFER EQUATION

CAO Qing-fa, HU Bin, WAN Shu,*WANG Ze-wen

（School of Science, East China University of Technology, Nanchang, Jiangxi 330013, China）

The inversion of the perfusion coefficient function of a class of bioheat transfer equations is studied in this paper. Based on the additional non-local conditions and the Crank-Nicolson method of finite difference, an iterative algorithm for reconstructing the perfusion coefficient function is constructed; after further simplification, an explicit scheme for retrieving perfusion coefficient is obtained. In order to overcome the instability of calculation, the moving average filtering method is introduced to denoise the error data. The results of calculation examples show that numerical inversion algorithms combined with the moving average filtering denoising are feasible and effective for retrieving perfusion coefficient function.

bioheat transfer equation; perfusion coefficient; inverse problem; finite difference, moving average

1674-8085（2022）02-0022-06

O29

A

10.3969/j.issn.1674-8085.2022.02.004

2021-08-01；

2021-09-18

國家自然科學基金項目(11961002，11761007)；江西省教育廳科技計劃項目(GJJ170444)；東華理工大學大學生科技創(chuàng)新基金項目

曹慶發(fā)(1996-)，男，江西贛州人，碩士生，主要從事一般反問題的計算方法研究(E-mail:cqingfa58@163.com)；

*王澤文(1974-)，男，江西上饒人，教授，博士，主要從事一般反問題的計算方法研究(E-mail:zwwang6@163.com).