CN106137200B - 从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法 - Google Patents

Abstract

一种从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法，包括步骤：(1)计算核磁共振射频发射点发射的总能量，减去系统反射回来的能量，得到人体组织内存在的电磁场的总能量；(2)根据B₁Mapping技术，得到射频发射产生的磁场B₁ ⁺场的幅度分布，根据得到B₁ ⁺场的能量；根据麦克斯韦方程组电磁互生理论，同时得到电场的能量；(3)通过比较内外圈各点之间的电磁场总能量差并进行计算得到各圈各点的局部比吸收率和得到各处损耗角正切值。本发明从电磁场能量角度求解组织电特性和组织局部比吸收率，具有计算简单，结果精确的特点。

Description

从电磁场能量传播角度求解组织电特性分布及局部比吸收率 的方法

技术领域

本发明涉及磁共振成像技术领域，特别是涉及一种从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法。

背景技术

电特性(简称EPs)又称介电特性，主要是指组织的电导率σ和电容率ε，该参数能够反映组织的生理、病理状态。已有研究证明正常组织和肿瘤组织存在明显的电特性差异。因此，通过比较活体组织与正常组织EPs分布图的差异，可以了解组织的生理、病理状态，对于疾病的诊断及肿瘤的超早期发现有十分重要的意义。

在成像的过程中，当磁场作用于人体组织时，活体组织会吸收射频能量。随着主磁场场强的增加，比吸收率也会增加。如果组织对射频能量吸收过多，超出人体自我调节的承受范围，就会出现局部灼伤，甚至体温调节紊乱，导致电解质代谢紊乱、神经系统紊乱等一系列问题。因此，量化测定超高场MRI全身组织Local SAR有重要的临床价值。

近年来，磁共振人体组织电特性断层成像(MR EPT)技术引起了科学家们极大的兴趣。该研究结合B₁mapping技术，通过麦克斯韦方程导出。到目前为止，科学家们已在3T和7T磁共振下取得了一系列研究成果。主要如下所示:2009年，Philips欧洲研发中心的Katscher等系统地研究了3T下的MR EPT技术，Katscher等提出了忽略RF磁场强度矢量中的Z方向分量，只利用另外两个分量来解析求解组织各处电导率和电容率的方法，并在鸟笼线圈中实现了电导率和电容率的活体测量，大大推动了3T MR EPT技术的发展。2011年，德国Karlsruhe生物医学工程学院的Voigt等在Katscher等的研究成果基础上，对3T RF场的幅度和相位做出了一定的假设，Voigt等基于一定的假设，提出了基于相位的电导率求解方法和基于磁场幅度的电容率求解方法。Voigt等提出的基于以上两种假设的MR EPT算法，经与Katscher等提出的算法比较，二者测量结果相差10％以内(但这并不意味着与真实值之间的误差也在10％以内)，但计算方法简化了很多，将MR EPT技术方法向临床实际应用推进了一步。GE全球研发中心的Bulumulla等在Voigt等人的工作基础上，亦在3T MR中实现了基于RF场幅度的组织电容率和基于RF场相位的组织电导率断层成像，并进一步提出了MR EPT算法的快速优化计算方法。美国明尼苏达大学磁共振中心的Bin He等学者在2013年利用人脑组织左右大致对称的特点，提出了RF发射磁场和接收磁场幅度大致相当的假设，从依托于质子密度分布的RF多线圈通道接收到的信号中抽取整个脑部的质子密度分布，并计算各通道接收场的幅度和相位，再结合采用7T超高场多通道RF线圈B₁Mapping技术得到的各通道的发射磁场，依照MR EPT二阶微分算法，实现了7T MR人体脑部组织无创EPs断层成像。另有美国纽约大学Bernard&Irene Schwartz生物医学成像中心的Sodickson等学者针对7T MR多通道线圈发射和接收场的相位无法直接测量的问题，提出了利用局部麦克斯韦断层成像(Local Maxwell Tomography,LMT)方法来求解各线圈通道的相位。他们将描述每一个通道的Maxwell电磁方程两边的实部和虚部分离，得到两组不同的方程。理论上，有5组以上发射和接收线圈通道的测量数据，就可以对质子密度分布、各发射通道和各接收通道的相位等未知量进行求解，并进而依照MR EPT二阶微分算法得到成像区域各像素点的组织EPs值分布。

与此同时，活体组织射频热点分布测定问题随着B₀场的增高变得不可回避且日益重要，目前并未提出完善的解决方案，所以当前国际上相关研究非常活跃。现有的射频热点分布量化测定方法主要有两大类，第一类是：基于人体电磁模型的仿真计算方法：该方法应用一定的人体电磁模型，将其置于特定的MR射频线圈下，设置一定的射频扫描参数，计算B₁场在人体组织内的分布，参考离体组织复介电参数，得到组织的射频热点分布。第二类是：基于线性MR人体组织复介电参数断层测量(MR EPT)技术的MR扫描测量计算方法：Katscher等在2009年提出了确定组织复电导率和局部射频热点的方法；Voigt等在2011年提出了脑部组织复介电参数的线性测量方法,即线性MR EPT，并在此研究基础上于2012年提出了针对活体组织的基于MR EPT的射频热点分布计算方，但是仅在脑部的活体MR扫描中得到验证，还不能应用到全身，并且计算结果误差在20％左右，准确度不够理想，还有待提高。知名教授Bin He领导的研究团队在2013年报道了在7T下测量脑部组织复介电参数的研究，为进一步采用MR EPT技术计算射频热点分布提供支持。

关于MR EPT，现有的MR EPT技术中存在以下几个问题：1.运算过程和重建结果受噪声影响较大；2.无法解决不同组织交界处存在的成像误差较大等问题；3.成像分辨率较差，有待进一步提高。导致以上问题的原因有以下三个方面：1.现有技术均采用了近似处理的方法，基于组织复电特性梯度为零的假设；2.重建核心算法基于麦克斯韦电磁方程二阶微分运算，这就使得噪声在计算过程中被放大，导致计算结果对噪声非常敏感；3.需要得到B₁ ⁺场的幅度和相位信息，但B₁ ⁺场的相位信息是无法通过测量得到的，只能估计得到，因此基于B₁ ⁺场EPT重建误差较大。

关于Local SAR，第一类基于人体电磁模型的射频热点分布仿真计算方法的最大缺陷是：针对单个病人个体建立个性化全身电磁模型的建模工作量巨大，无法在实践中推广，因此只能采用通用人体电磁模型，这就导致个体差异无法测量。除此之外，计算电磁模型时参考离体组织的复介电参数，而离体组织的复介电参数与活体组织的复介电参数之间存在一定差异，这种差异也会导致计算结果出现偏差。第二类基于线性MR EPT的MR扫描计算射频热点分布的方法，均要求测量的目标区域内组织复介电参数的变化远远小于电场的变化。这种要求对以软组织为主的脑部组织而言，基本满足；但是对全身组织而言，就不能满足，因为全身范围内组织的复介电参数相互差异较大。

因此，针对现有技术不足，提供一种从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法以克服现有技术不足甚为必要。

发明内容

本发明的目的在于避免现有技术的不足之处而提供一种从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法，该从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法具有计算接单、结果准确特点。

本发明的上述目的通过如下技术手段实现。

提供一种从电磁场能量传播角度求解组织电特性分布的方法，该方法基于以下物理事实：电磁能量在射频线圈的作用下，从四周向体内传播，电场和磁场能量交替转换，且能量不可能突变，在电场和磁场能量相互转换过程中，一部分电场的能量会转换成热能消耗掉；

核磁共振射频发射线圈发射出的电磁场，在到达人体后，各点电磁场的总能量的表达式为：

其中，ω为角频率，μ是磁导率，ε是电容率，σ是电导率，B为磁感应强度，E为电场强度，“*”表示取共轭复数，表示磁场储能，表示电场储能，表示焦耳热损耗；

具体通过如下步骤进行：

(1)根据射频功放的功率大小和射频发射的时间，以及系统的工作效率，计算核磁共振射频线圈发射的总能量，减去系统反射回来的能量，得到人体组织内存在的电磁场的总能量；

(2)根据B₁Mapping技术，得到射频发射产生的磁场B₁ ⁺场在人体内各点的幅度分布，根据得到B₁ ⁺场在体内各点处的磁场能量；根据麦克斯韦方程组电磁互生理论，同时得到该处电场的能量；

(3)能量由外圈向内圈传递的过程中，由于热能的损耗，总能量不断减少，每圈对应各点之间减少的能量即是两圈之间各对应点之间B₁ ⁺场的能量的差，即：

其中,B_out表示外圈的磁感应强度,B_in表示内圈的磁感应强度,Ε_out表示外圈电场强度；B_in和B_out通过B₁Mapping技术得到；

外圈各点的电场储能转化为内圈对应各点的磁场能量，即：

用式(Ⅱ)除以式(Ⅲ)，得到损耗角正切，损耗角正切为外圈各点所在部位组织的电导率与相对介电常数的比值，相对介电常数为电容率与角频率之积，损耗角正切的表达式为

通过上述计算，得到成像区域内各处组织的电导率与相对介电常数的比值

本发明同时提供一种从电磁场能量传播角度求解组织局部比吸收率的方法，

该方法基于以下物理事实：电磁能量在射频线圈的作用下，从四周向体内传播，电场和磁场能量交替转换，且能量不可能突变，在电场和磁场能量相互转换过程中，一部分电场的能量会转换成热能消耗掉；

核磁共振射频发射线圈发射出的电磁场，在到达人体后，各点电磁场的总能量的表达式为：

其中，ω为角频率，μ是磁导率，ε是电容率，σ是电导率，B为磁感应强度，E为电场强度，“*”表示取共轭复数，表示磁场储能，表示电场储能，表示焦耳热损耗；

具体通过如下步骤进行：

(1)根据射频功放的功率大小和射频发射的时间，以及系统的工作效率，计算核磁共振射频线圈发射的总能量，减去系统反射回来的能量，得到人体组织内存在的电磁场的总能量；

(2)根据B₁Mapping技术，得到射频发射产生的磁场B₁ ⁺场在人体内各点的幅度分布，根据得到B₁ ⁺场在体内各点处的磁场能量；根据麦克斯韦方程组电磁互生理论，同时得到该处电场的能量；

(3)能量由外圈向内圈传递的过程中，由于热能的损耗，总能量不断减少，每圈对应各点之间减少的能量即是两圈之间各对应点之间B₁ ⁺场的能量的差，即：

其中，B_out表示外圈的磁感应强度，B_in表示内圈的磁感应强度，Ε_out表示外圈电场强度；B_in和B_out通过B₁Mapping技术得到；

外圈各点的电场储能转化为内圈对应各点的磁场能量，即：

对各圈而言，电场和磁场能量相互转化，即：

通过比较内外圈各点之间的电磁场总能量差，即得到各圈各点的热能损耗的值，即各圈各点的再利用人体组织各点处的组织密度ρ，得到各圈各点的局部比吸收率Local SAR，即：

本发明的从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法，从磁共振射频发射电磁场能量在体内传输和分布的角度求解EPs，计算过程中只用到了可以直接测量得到的射频场的幅度信息，不需要射频场的相位信息，规避了射频场的相位无法直接测量的问题；同时新方法避免了现有方法采用的麦克斯韦方程的二阶微分运算的方法，可以有效的提高运算结果的精度。

本发明的从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法，从电磁场能量——一个全新的角度求解组织电特性和组织局部比吸收率，具有计算简单，结果精确的特点。

附图说明

利用附图对本发明作进一步的说明，但附图中的内容不构成对本发明的任何限制。

图1是本发明实施例3的仿真实验中射频场幅度和相位的分布图。

图2是本发明实施例3的仿真实验中得到的SAR值分布图。

图3是本发明实施例3的仿真实验中计算得到的电导率和相对电容率的分布图。

图4是本发明实施例3的仿真实验中电场能量、磁场能量、热损耗和总能量的分布图。

具体实施方式

结合以下实施例对本发明作进一步描述。

实施例1。

一种从电磁场能量传播角度求解组织电特性分布的方法，该方法基于以下物理事实：电磁能量在射频线圈的作用下，从四周向体内传播，电场和磁场能量交替转换，且能量不可能突变，在电场和磁场能量相互转换过程中，一部分电场的能量会转换成热能消耗掉；

核磁共振射频发射点磁场的总能量的表达式为：

其中，ω为角频率，μ是磁导率，ε是电容率，σ是电导率，B为磁感应强度，E为电场强度，“*”表示取共轭复数，表示磁场储能，表示电场储能，表示焦耳热损耗。

具体通过如下步骤进行：

(1)根据射频功放的功率大小和射频发射的时间，以及系统的工作效率，计算核磁共振射频线圈发射的总能量，减去系统反射回来的能量，得到人体组织内存在的电磁场的总能量；其中，反射回来的能量可通过测量得到。

(2)根据B₁Mapping技术，得到射频发射产生的磁场B₁ ⁺场在人体内各点的幅度分布，根据得到B₁ ⁺场在体内各点处的磁场能量；根据麦克斯韦方程组电磁互生理论，同时得到该处电场的能量。

(3)能量由外圈向内圈传递的过程中，由于热能的损耗，总能量不断减少，每圈对应各点之间减少的能量即是两圈之间各对应点之间B₁ ⁺场的能量的差，即：

外圈的电场能量转化为内圈的磁场能量，即：

用式(Ⅱ)除以式(Ⅲ)，得到损耗角正切，损耗角正切为外圈各点所在部位组织的电导率与相对介电常数的比值，相对介电常数为电容率与角频率之积，损耗角正切的表达式为

通过上述计算，得到成像区域内各处组织的电导率与相对介电常数的比值

由于正常人体组织的介电特性值是已知的，故其比值即损耗角正切的值，也是已知的。将上述计算得到的分布成像的结果与正常组织的值进行对比，就可以发现值异常所在部位，为疾病的早期发现，如早期可疑癌变病灶等，提供有价值的空间定位信息，再辅以活检等技术手段，就可以实现肿瘤的早期发现，具有重要的临床应用价值。

本发明的从电磁场能量传播角度求解组织电特性分布的方法，从磁共振射频发射电磁场能量在体内传输和分布的角度求解EPs，计算过程中只用到了可以直接测量得到的射频场的幅度信息，不需要射频场的相位信息，规避了射频场的相位无法直接测量的问题；同时新方法避免了麦克斯韦方程的二阶微分运算，可以有效的提高运算结果的精度。

相比现有的MR EPT的算法中，在计算MR射频电磁场时，普遍采用的处理方法是提出局部组织电特性是均匀的(即假设局部组织的电特性的梯度为零)，并且假设忽略射频电场的横截面分量和射频磁场的纵轴分量，即假设E_x,y(r)≈0和B_z(r)≈0。现有的MR EPT算法这样处理的目的是为了简化得到组织EPs的近似计算公式，形成现有的MR EPT算法理论，比如之前提到的Voigt等给出的基于射频磁场幅度的电容率求解方法。这种简化处理在复杂电磁边界条件存在区域会导致测量结果误差增加。

本发明的方法从电磁场能量的角度求解EPs，计算过程中只应用了磁场的幅度信息，避开了射频场的相位问题。相位问题是现有MR EPT算法中未解决的难题，因为射频发射场的相位利用现有的MR技术，是无法在工程上直接测量得到的，在求解的过程中，同样是在做了许多假设的情况下求取相位分布近似解，且现有的求解方法都存在较大的误差。故，本发明的方法同时避免了麦克斯韦方程的二阶微分运算，可以有效的提高运算结果的精度。

实施例2。

一种从电磁场能量传播角度求解组织局部比吸收率的方法，该方法基于以下物理事实：电磁能量在射频线圈的作用下，从四周向体内传播，电场和磁场能量交替转换，且能量不可能突变，在电场和磁场能量相互转换过程中，一部分电场的能量会转换成热能消耗掉。

核磁共振射频发射点磁场的总能量的表达式为：

其中，ω为角频率，μ是磁导率，ε是电容率，σ是电导率，B为磁感应强度，E为电场强度，“*”表示取共轭复数，表示磁场储能，表示电场储能，表示焦耳热损耗。

具体通过如下步骤进行：

(1)根据射频功放的功率大小和射频发射的时间，以及系统的工作效率，计算核磁共振射频线圈发射的总能量，减去系统反射回来的能量，得到人体组织内存在的电磁场的总能量。

(2)根据B₁Mapping技术，得到射频发射产生的磁场B₁ ⁺场在人体内各点的幅度分布，根据得到B₁ ⁺场在体内各点处的磁场能量；根据麦克斯韦方程组电磁互生理论，同时得到该处电场的能量。

(3)能量由外圈向内圈传递的过程中，由于热能的损耗，总能量不断减少，每圈对应各点之间减少的能量即是两圈之间各对应点之间B₁ ⁺场的能量的差，即：

外圈的电场能量转化为内圈的磁场能量，即：

对外圈而言：

其中,B_out表示外圈的磁感应强度,B_in表示内圈的磁感应强度,Ε_out表示外圈电场强度。B_in和B_out通过B₁Mapping技术得到。

通过比较内外圈各点之间的电磁场总能量差，即可以得到各圈各点的热能损耗的值，即各圈各点的再利用人体组织各点处的组织密度ρ(已知值)，即可得到各圈各点的局部比吸收率(Local SAR)值，即：

本发明的从电磁场能量传播角度求解组织局部比吸收率的方法，从磁共振射频发射电磁场能量在体内传输和分布的角度求解EPs，计算过程中只用到了可以直接测量得到的射频场的幅度信息，不需要射频场的相位信息，规避了射频场的相位无法直接测量的问题；同时新方法避免了麦克斯韦方程的二阶微分运算，可以有效的提高运算结果的精度。在组织电特性分布的基础上，结合组织密度、组织处电场幅值、组织处电导率，按照Local SAR计算公式即可以得到Local SAR分布。

本发明的从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法，从电磁场能量——一个全新的角度求解组织电特性和组织局部比吸收率，具有计算简单，结果精确的特点。

实施例3。

为了验证本发明的准确性，通过仿真、体模实验和人体实际测量等不同角度，来验证本发明方法的正确性。

1)仿真

利用SEMCAD软件，设置核磁共振的频率为128MHZ，建立16通道的高通鸟笼线圈，将DUKE头部模型置于鸟笼线圈中。仿真获取头部射频场B₁ ⁺场的分布，如图1所示。建立磁场由外向内的传播模型,将所得的B代入公式(Ⅱ)中,得到组织的热损耗值,又由于组织的密度是已知的,带入Local SAR的求解公式即可求得Local SAR的分布图，如图2所示。由公式(Ⅰ)和(Ⅳ)可计算得组织的电容率，再根据损耗角正切值，则可求出对应点的电导率，如图3所示。图4是本发明实施例3的仿真实验中电场能量、磁场能量、热损耗和总能量的分布图。

2)体模实验

采用各种能模拟不同组织电特性的介电材料，填充特定结构空腔，形成类似人体不同组织交界面的电磁边界条件，制成专用的体模。在3T的磁共振中采用不同的扫描序列，得到相应的图像，应用B₁Mapping技术，求出射频场的幅度|B₁ ⁺|。然后，结合核磁共振的能量场求解电特性及局部比吸收率即可。

3)人体实测

通过仿真及体模实验验证基于核磁共振能量场求解电特性的准备性后，进行临床志愿者的人体测量实验。志愿者包含正常人群及肿瘤患者两类。

通过仿真、体模实验和人体实测，结果证明，本发明的方法能够简单、精确地得到组织的电特性分布和局部比吸收率。相比现有的MR EPT方法，本发明的方法从电磁场能量的角度求解EPs，计算过程中只应用了磁场的幅度信息，避开了射频场的相位问题，因为相位在求解的过程中，比较复杂，且现有的求解方法都存在较大的误差；同时避免了麦克斯韦方程的二阶微分运算，可以有效的提高运算结果的精度。在组织电特性分布的基础上，结合组织密度、组织处电场幅值、组织处电导率，按照Local SAR计算公式即可以得到Local SAR分布。

本发明的从电磁场能量传播角度求解组织电特性分布及局部比吸收率的方法，从电磁场能量——一个全新的角度求解组织电特性和组织局部比吸收率，具有计算简单，结果精确的特点。

最后应当说明的是，以上实施例仅用以说明本发明的技术方案而非对本发明保护范围的限制，尽管参照实施例对本发明作了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和范围。

Claims (2)

1.一种从电磁场能量传播角度求解组织电特性分布的方法，其特征在于，

该方法基于以下物理事实：电磁能量在射频线圈的作用下，从四周向体内传播，电场和磁场能量交替转换，且人体组织内存在的电磁场的总能量不可能突变，在电场和磁场能量相互转换过程中，一部分电场的能量会转换成热能消耗掉；

核磁共振射频线圈发射出的电磁场，在到达人体后，人体组织内各点电磁场的总能量的表达式为：

其中，ω为角频率，μ是磁导率，ε是电容率，σ是电导率，B为磁感应强度，E为电场强度，“*”表示取共轭复数，表示磁场储能，表示电场储能，表示焦耳热损耗；

具体通过如下步骤进行：

(1)根据射频功放的功率大小和射频发射的时间，以及MRI设备系统的工作效率，计算核磁共振射频线圈发射的能量，减去MRI设备系统的射频电缆反射回来的能量，得到人体组织内存在的电磁场的总能量；

(2)根据B₁ Mapping技术，得到射频发射产生的磁场B₁ ⁺场在人体内各点的幅度分布，根据得到B₁ ⁺场在体内各点处的磁场能量；根据麦克斯韦方程组电磁互生理论，同时得到体内各点对应的电场的能量；

(3)总能量由外圈向内圈传递的过程中，由于热能的损耗，总能量不断减少，两圈之间减少的能量即是两圈之间对应点之间B₁ ⁺场的能量的差，即：

其中,B_out表示外圈的磁感应强度,B_in表示内圈的磁感应强度,Ε_out表示外圈电场强度；B_in和B_out通过B₁ Mapping技术得到；

外圈各点的电场储能转化为内圈对应各点的磁场能量，即：

用式(Ⅱ)除以式(Ⅲ)，得到损耗角正切，损耗角正切为外圈各点所在部位组织的电导率与相对介电常数的比值，相对介电常数为电容率与角频率之积，损耗角正切的表达式为

通过上述计算，得到成像区域内各处组织的电导率与相对介电常数的比值

2.一种从电磁场能量传播角度求解组织局部比吸收率的方法，其特征在于，

该方法基于以下物理事实：电磁能量在射频线圈的作用下，从四周向体内传播，电场和磁场能量交替转换，且人体组织内存在的电磁场的总能量不可能突变，在电场和磁场能量相互转换过程中，一部分电场的能量会转换成热能消耗掉；

核磁共振射频线圈发射出的电磁场，在到达人体后，人体组织内各点电磁场的总能量的表达式为：

其中，ω为角频率，μ是磁导率，ε是电容率，σ是电导率，B为磁感应强度，E为电场强度，“*”表示取共轭复数，表示磁场储能，表示电场储能，表示焦耳热损耗；

具体通过如下步骤进行：

(1)根据射频功放的功率大小和射频发射的时间，以及MRI设备系统的工作效率，计算核磁共振射频线圈发射的能量，减去MRI设备系统的射频电缆反射回来的能量，得到人体组织内存在的电磁场的总能量；

(2)根据B₁ Mapping技术，得到射频发射产生的磁场B₁ ⁺场在人体内各点的幅度分布，根据得到B₁ ⁺场在体内各点处的磁场能量；根据麦克斯韦方程组电磁互生理论，同时得到体内各点对应的电场的能量；

(3)总能量由外圈向内圈传递的过程中，由于热能的损耗，总能量不断减少，两圈之间减少的能量即是两圈之间对应点之间B₁ ⁺场的能量的差，即：

其中，B_out表示外圈的磁感应强度，B_in表示内圈的磁感应强度，Ε_out表示外圈电场强度；B_in和B_out通过B₁ Mapping技术得到；

外圈各点的电场储能转化为内圈对应各点的磁场能量，即：

对各圈而言，电场和磁场能量相互转化，即：

通过比较内外圈各点之间的电磁场总能量差，即得到各圈各点的热能损耗的值，即各圈各点的再利用人体组织各点处的组织密度ρ，得到各圈各点的局部比吸收率Local SAR，即：