CN114355263B - 一种高阶匀场线圈设计方法 - Google Patents
一种高阶匀场线圈设计方法 Download PDFInfo
- Publication number
- CN114355263B CN114355263B CN202011047105.4A CN202011047105A CN114355263B CN 114355263 B CN114355263 B CN 114355263B CN 202011047105 A CN202011047105 A CN 202011047105A CN 114355263 B CN114355263 B CN 114355263B
- Authority
- CN
- China
- Prior art keywords
- coil
- order
- layer
- coils
- framework
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 14
- 238000004804 winding Methods 0.000 claims abstract description 11
- 238000009826 distribution Methods 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000012938 design process Methods 0.000 claims description 2
- 239000010410 layer Substances 0.000 abstract description 41
- 230000008878 coupling Effects 0.000 abstract description 18
- 238000010168 coupling process Methods 0.000 abstract description 18
- 238000005859 coupling reaction Methods 0.000 abstract description 18
- 239000002356 single layer Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004323 axial length Effects 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种高阶匀场线圈设计方法,适用于三阶或三阶以上匀场线圈。本发明采用多层线圈骨架绕线,每层骨架间隔一定的距离,每层上的导线串联,并且匀场线圈与其他线圈的互感足够小。该方法的主要优势是:从三阶线圈开始,匀场线圈之间的耦合问题变突出;除此之外,高阶线圈为了产生足够强度的场需要的线圈匝数比低阶线圈更多更稠密。这两个问题采用单层绕线结构都很难克服。为了解决这两个问题,本发明采用多层骨架绕线,通过合理选择绕线骨架位置,并在线圈优化时将线圈耦合做为约束条件进行优化,这种做法能使高阶线圈与低阶线圈之间的互感大幅度降低,并且避免了导线绕线密集造成的绕线困难。
Description
技术领域
本发明涉及一种高阶匀场线圈设计方法,更具体地涉及一种匀场线圈分层布线方案与设计方法。
背景技术
在磁共振系统中,要求背景场具有很高的均匀度。这种均匀度一般通过两种方法保证:被动匀场与主动匀场。所谓主动匀场,即设计若干个匀场线圈,每个匀场线圈对应一个谐波系数。通过调节匀场线圈内的电流可以将背景场对应的谐波分量进行抵消,从而提高背景场的均匀度。因此,匀场线圈在核磁共振系统中是一个很重要的部件。但是当线圈阶数达到3阶或者3阶以上时,匀场线圈的设计变得非常困难。这种困难表现在两方面。其一是高阶匀场线圈与低阶匀场线圈或梯度线圈之间可能会存在比较大的耦合,其次线圈阶数越高,则导线越密集,工程制作难度加大。除此之外,还有一个问题,采用流函数法设计高阶匀场线圈时,优化后的电流方向通常变化很快,这时候将电流密度转换为导线时会存在较大误差。这几个问题都为高阶匀场线圈的设计与应用造成了困难。
发明内容
发明目的:为了解决上述现有技术存在的问题,本发明提出一种新的设计方案,采用两层或两层以上骨架对高阶匀场线圈进行绕线,并给出了一种高阶匀场线圈设计方法。采用该方法能够保证所设计的线圈与已知线圈之间的耦合足够的小,并且能够保证每层线圈绕线的稀疏度。
技术方案:为实现上述发明目的,本发明所述的一种高阶匀场线圈设计方法,已知低阶匀场线圈结构以及该高阶匀场线圈设计参数,在允许的布线空间内选择两层或两层以上线圈骨架,每层骨架间隔一定的距离,所有骨架上的线圈串联,在线圈优化时将该高阶匀场线圈与低阶匀场线圈的互感作为约束条件,保证线圈的场满足要求的同时,线圈与可能存在耦合的低阶匀场线圈的互感满足工程要求。
进一步,采用流函数法设计线圈结构,设计过程为:已知线圈骨架的层数与每层线圈骨架的尺寸,在每层线圈骨架上定义电流密度基函数,然后构造目标函数与约束条件,通过对优化模型求解得到每层骨架上的电流密度分布,进一步通过电流密度函数得到高阶匀场线圈的结构。
进一步,设计方案采用两个线圈骨架,一个骨架位于梯度线圈主线圈附近,另一个位于梯度线圈屏蔽线圈附近。
进一步,线圈优化的目标如下:
上式中,An,m、Bn,m为所设计的线圈产生的磁场的谐波系数,A'n,m、B'n,m为理想的谐波系数,为待设计线圈骨架上任一点处的矢量坐标,/>为第t层待设计线圈骨架上坐标点/>处的电流密度,αn,m、βn,m、ω为权重系数,N为谐波阶数;St为第t层待设计线圈骨架所在的面,T为线圈骨架层数,ds为待设计线圈骨架内的面积微元。
进一步,线圈优化需要满足如下约束条件:
上式中,M为已知的低阶匀场线圈个数,Lj表示第j个已知低阶匀场线圈绕线路径,为第j个已知低阶匀场线圈内任一点处的矢量坐标,Indj表示高阶匀场线圈与第j个已知低阶匀场线圈的互感,εj表示允许的最大互感值,I0为待设计线圈内的电流幅值,μ0为真空磁导率,/>为已知低阶匀场线圈内的线微元。
进一步,所述线圈为柱面Z3匀场线圈,在每层骨架上采用的电流密度基函数表达式如下:
上式中,P为每层骨架上的基函数个数,Lt为第t个骨架长度,为待设计线圈骨架上任一点处的矢量坐标。
进一步,采用内点法对上述优化问题进行优化。
进一步,采用如下公式定义谐波系数:
当m=0时
Bn0=0
当m=1,2,3…时
上式中,r0为人为选取的球形参考球面半径,为参考球面上任意一点的矢量坐标,/>为参考面上的球坐标,/>表示球面上/>处的z方向磁场分量,sinθ、cosθ分别表示正弦函数与余弦函数,Pnm(cosθ)为n阶m次勒让德多项式。
有益效果:本发明设计方法的优点是:采用多层线圈骨架布线,并合理选择每层匀场线圈骨架位置,然后通过约束条件约束该线圈与其他已知线圈之间的耦合度,可以实现所设计的线圈与其他线圈之间具有非常小的耦合,从而避免线圈之间相互干扰对成像造成影响,而且能够保证线圈导线的稀疏度,减小了绕线困难。
附图说明
图1为已知的轴向梯度线圈结构图。
图2为单层Z3匀场线圈结构图。
图3为两层Z3匀场线圈结构图。
具体实施方式
本发明是关于一种分层高阶匀场线圈设计方法。发明内容是提出了一种基于多层绕线骨架的线圈设计方案。在该方案中,首先在允许的布线空间内合理选择每层线圈绕线骨架,然后分别在每层骨架上定义一组电流密度表达式,将该区域内的电流密度函数表示为基函数与待求系数的组合,然后定义待优化的目标函数以及约束条件,将线圈优化问题转化为一个带约束的优化问题,然后采用合适的优化算法进行求解。下面我们以超导磁共振系统中的三阶轴向匀场线圈设计来说明该方法的具体实施方案。
在超导磁体系统中,梯度线圈与匀场线圈都位于梯度线圈部件内,梯度线圈与匀场线圈均为柱面结构,每个线圈分布在不同的圆柱面骨架上。在工程中,梯度线圈同时用来做1阶匀场线圈。因此,本发明中所指的匀场线圈包括梯度线圈。原则上,要求所有匀场线圈之间不存在耦合。但是高阶匀场线圈与低阶匀场线圈之间可能会存在很大的耦合度。这里提供一种高阶匀场线圈设计方法,能够很好的去除高阶匀场线圈与低阶匀场线圈之间的耦合。假定已知已设计好的M个低阶匀场线圈的结构信息,接下来设计一种新的匀场线圈,要求与已设计的匀场线圈的耦合满足一定要求。则设计方案如下:
首先给定设计所需的输入参数,包括匀场线圈骨架层数T,每层线圈骨架所在柱面的半径r1,r2,…rT,轴线方向布线长度L1,2L2,…LT,磁场强度G与纯度E。为了能够去掉线圈之间的耦合,每层线圈半径要选取合适的值。这里以Z3线圈设计为例进行说明。Z3匀场线圈与Z梯度线圈之间可能存在耦合,与其他低阶匀场线圈不存在耦合。因此设计Z3时只需要考虑与Z线圈之间的耦合度即可。而且Z梯度线圈一般分为两层,一层主线圈,一层屏蔽线圈,两层线圈上的电流方向相反。为了减小Z3线圈与Z线圈之间的耦合,Z3线圈选取两层骨架,一层骨架选取在靠近Z线圈主线圈位置,另一层靠近Z屏蔽线圈位置。
选定骨架层数与位置后,接下来在每层线圈骨架上构造电流密度基函数:
上式中,P为每层骨架上的基函数个数,Lt为第t个骨架长度,为骨架上任一点的矢量坐标。
每层上的电流密度可以表示为基函数与其系数的组合:
接下来构造待优化的目标函数:
上式中,An,m、Bn,m为所设计的线圈产生的磁场的谐波系数,A'n,m、B'n,m为理想的谐波系数,为第t层待设计线圈骨架上坐标点/>处的电流密度,αn,m、βn,m、ω为权重系数,N为谐波阶数;St为第t层待设计线圈骨架所在的面,ds为待设计线圈骨架内的面积微元。
本发明中,采用如下公式定义谐波系数:
当m=0时
Bn0=0
当m=1,2,3…时
上式中,r0为人为选取的球形参考球面半径,为参考球面上任意一点的矢量坐标,/>为参考面上的球坐标,/>表示球面上/>处的z方向磁场分量,sinθ、cosθ分别表示正弦函数与余弦函数,Pnm(cosθ)为n阶m次勒让德多项式。
为了减小匀场线圈与其他已知线圈之间的耦合,在优化上述函数时添加以下约束条件:
上式中,M为已知的低阶匀场线圈个数,Lj表示第j个已知低阶匀场线圈绕线路径,为第j个已知低阶匀场线圈内任一点处的坐标矢量,Indj表示高阶匀场线圈与第j个已知低阶匀场线圈的互感,εj表示允许的最大互感值,I0为待设计线圈内的电流幅值,μ0为真空磁导率,/>为已知低阶匀场线圈内的线微元。
在添加上述约束条件后,则线圈设计问题变为一个带约束的优化问题,该数学模型的求解可以采用多种算法,例如内点法等,这些方法在很多论文里都有介绍,这里不再赘述。
求目标函数在上述约束条件下的最小值,可得到骨架上期望的电流密度分布,进一步可得到骨架上的线圈结构。
下面根据具体算例说明该算法的效果。已知某一带屏蔽的柱面轴向梯度线圈的结构如图1所示。图中的虚线表示主线圈,实线表示屏蔽线圈。主线圈的半径为0.35m,屏蔽线圈半径为0.45m,线圈的轴向长度为1.28m,在直径为45cm的球形范围内,梯度场强为50μT/m/A,线性度为10%。下面设计一款与该梯度线圈配套的Z3匀场线圈。Z3线圈的位置选取在梯度线圈主线圈外侧,屏蔽线圈内侧。设计的指标为:在直径为45cm的球内,主要的谐波分量为10μT/m/A,其他谐波纯度不大于4%。
首先,我们设计时只选取一层线圈骨架,骨架半径为r=0.36m,长度为1.2m。设计得到的匀场线圈如图2所示。该匀场线圈与已知梯度线圈之间的互感为:36.7μH。而且经过我们大量的实验,该互感很难通过在优化模型中添加互感约束条件减小。在医用磁共振应用中,一般要求梯度线圈与匀场线圈之间的耦合小于1μH,互感越大对成像造成的干扰越严重。可以看出,该匀场线圈与梯度线圈的耦合度不符合工程应用要求。另外从图中我们可以看出,该线圈中的导线非常密集。这限制了导线的尺寸,并且不易于工程制作。
下面我们选取两层骨架重新设计该匀场线圈。两层骨架半径分别为r1=0.36m,r2=0.44m,轴向长度分别为1.2m,1.24m。设计得到的结果如图3所示。该线圈与轴向梯度线圈的互感为0.1μH。同时,通过双层布线的方式,导线密度比单层下降了。可以看出,通过本发明中的算法可以很好的实现线圈解耦。
本实施方案以梯度线圈与匀场线圈为例来说明本申请中的高阶匀场线圈设计方案。但是该方法不局限于柱面匀场线圈设计。任何采用本方法设计的线圈都属于本申请保护范围。
Claims (5)
1.一种高阶匀场线圈设计方法,其特征在于,已知低阶匀场线圈结构以及该高阶匀场线圈设计参数,在允许的布线空间内选择两层线圈骨架,每层骨架间隔一定的距离,一个骨架位于梯度线圈主线圈附近,另一个位于梯度线圈屏蔽线圈附近,两层骨架上的线圈串联,在线圈优化时将该高阶匀场线圈与低阶匀场线圈的互感作为约束条件,保证线圈的场满足要求的同时,线圈与可能存在耦合的低阶匀场线圈的互感满足工程要求;
线圈优化的目标如下:
上式中,An,m、Bn,m为所设计的线圈产生的磁场的谐波系数,A'n,m、B'n,m为理想的谐波系数,为待设计线圈骨架上任一点处的矢量坐标,/>为第t层待设计线圈骨架上坐标点/>处的电流密度,αn,m、βn,m、ω为权重系数,N为谐波阶数;St为第t层待设计线圈骨架所在的面,T为线圈骨架层数,ds为待设计线圈骨架内的面积微元;
线圈优化需要满足如下约束条件:
上式中,M为已知的低阶匀场线圈个数,Lj表示第j个已知低阶匀场线圈绕线路径,为待设计线圈骨架上任一点处的矢量坐标,/>为第j个已知低阶匀场线圈内任一点处的矢量坐标,Indj表示高阶匀场线圈与第j个已知低阶匀场线圈的互感,εj表示允许的最大互感值,为坐标点/>处的电流密度,I0为待设计线圈内的电流幅值,μ0为真空磁导率,/>为已知低阶匀场线圈内的线微元,St为第t层待设计线圈骨架所在的面,T为线圈骨架层数,ds为待设计线圈骨架内的面积微元。
2.根据权利要求1所述的高阶匀场线圈设计方法,其特征在于,采用流函数法设计高阶匀场线圈结构,设计过程为:已知线圈骨架的层数与每层线圈骨架的尺寸,在每层线圈骨架上定义电流密度基函数,然后构造目标函数与约束条件,通过对优化模型求解得到每层骨架上的电流密度分布,进一步通过电流密度函数得到高阶匀场线圈的结构。
3.根据权利要求2所述的高阶匀场线圈设计方法,其特征在于,所述高阶匀场线圈为柱面Z3匀场线圈,在每层骨架上采用的电流密度基函数表达式如下:
上式中,T为线圈骨架层数,P为每层骨架上的基函数个数,Lt为第t个骨架长度,为待设计线圈骨架上任一点处的矢量坐标。
4.根据权利要求2所述的高阶匀场线圈设计方法,其特征在于,采用内点法对优化问题进行优化。
5.根据权利要求1所述的高阶匀场线圈设计方法,其特征在于,采用如下公式定义谐波系数:
当m=0时
Bn0=0
当m=1,2,3…时
上式中,r0为人为选取的球形参考球面半径,为参考球面上任意一点的矢量坐标,为参考面上的球坐标,/>表示球面上/>处的z方向磁场分量,sinθ、cosθ分别表示正弦函数与余弦函数,Pnm(cosθ)为n阶m次勒让德多项式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011047105.4A CN114355263B (zh) | 2020-09-29 | 2020-09-29 | 一种高阶匀场线圈设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011047105.4A CN114355263B (zh) | 2020-09-29 | 2020-09-29 | 一种高阶匀场线圈设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114355263A CN114355263A (zh) | 2022-04-15 |
CN114355263B true CN114355263B (zh) | 2023-10-24 |
Family
ID=81090319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011047105.4A Active CN114355263B (zh) | 2020-09-29 | 2020-09-29 | 一种高阶匀场线圈设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114355263B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020429A (zh) * | 2014-06-06 | 2014-09-03 | 南京工程学院 | 一种梯度线圈并联分层的布线结构和布线方法 |
CN106556813A (zh) * | 2016-11-25 | 2017-04-05 | 上海辰光医疗科技股份有限公司 | 一种磁共振系统中主动匀场线圈的线性混合优化方法 |
CN106990373A (zh) * | 2017-03-28 | 2017-07-28 | 中国科学院电工研究所 | 一种磁共振系统的解耦轴向匀场线圈设计方法 |
CN108872896A (zh) * | 2018-05-29 | 2018-11-23 | 河北惠仁医疗设备科技有限公司 | 一种双平面磁共振成像系统梯度线圈的设计方法 |
CN110568390A (zh) * | 2019-10-14 | 2019-12-13 | 惠仁望都医疗设备科技有限公司 | 一种双平面匀场线圈的设计方法 |
-
2020
- 2020-09-29 CN CN202011047105.4A patent/CN114355263B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020429A (zh) * | 2014-06-06 | 2014-09-03 | 南京工程学院 | 一种梯度线圈并联分层的布线结构和布线方法 |
CN106556813A (zh) * | 2016-11-25 | 2017-04-05 | 上海辰光医疗科技股份有限公司 | 一种磁共振系统中主动匀场线圈的线性混合优化方法 |
CN106990373A (zh) * | 2017-03-28 | 2017-07-28 | 中国科学院电工研究所 | 一种磁共振系统的解耦轴向匀场线圈设计方法 |
CN108872896A (zh) * | 2018-05-29 | 2018-11-23 | 河北惠仁医疗设备科技有限公司 | 一种双平面磁共振成像系统梯度线圈的设计方法 |
CN110568390A (zh) * | 2019-10-14 | 2019-12-13 | 惠仁望都医疗设备科技有限公司 | 一种双平面匀场线圈的设计方法 |
Non-Patent Citations (2)
Title |
---|
A Stream Function Method for Gradient Coil Design;ROSTISLAV A;Concepts in Magnetic Resonance Part B Magnetic Resonance Engineering;67-80 * |
Michael Poole.Improved equipment and techniques for dynamic shimming in high field MRI.2007,59-67. * |
Also Published As
Publication number | Publication date |
---|---|
CN114355263A (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5695241B2 (ja) | 傾斜磁場コイルおよび核磁気共鳴撮像装置 | |
CN102176368B (zh) | 一种用于磁共振成像超导磁体的优化设计方法 | |
CN106990373B (zh) | 一种磁共振系统的解耦轴向匀场线圈设计方法 | |
CN110780245B (zh) | 用于平面超导磁共振系统的屏蔽梯度线圈设计方法及其梯度线圈 | |
CN104198969B (zh) | 一种梯度线圈设计方法 | |
US20190353725A1 (en) | Resistive electromagnet systems and methods | |
CN104020429A (zh) | 一种梯度线圈并联分层的布线结构和布线方法 | |
JP2013532495A (ja) | 小型超伝導磁石装置 | |
JP2009502031A (ja) | 多層磁石 | |
CN110568390A (zh) | 一种双平面匀场线圈的设计方法 | |
CN109765510A (zh) | 一种带有圆角的径向超导匀场线圈及其设计方法 | |
CN114355263B (zh) | 一种高阶匀场线圈设计方法 | |
CN104062613A (zh) | 一种有源屏蔽梯度线圈及其设计方法 | |
US7427908B1 (en) | Magnetic shimming configuration with optimized turn geometry and electrical circuitry | |
KR20190005709A (ko) | 일정한 피치를 가지는 수신 코일을 포함하는 무선 전력 전송 시스템 | |
JPS6387712A (ja) | 磁界コイル | |
JP2014000388A (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用の磁石 | |
JPWO2016084164A1 (ja) | 超電導マグネットおよび超電導マグネットの製造方法 | |
US10217555B2 (en) | Compact inductor | |
Niu et al. | A novel design method of independent zonal superconducting shim coil | |
JP2005294537A (ja) | パネル状コイル、一様磁界発生装置、勾配磁界発生装置および磁界キャンセリング装置 | |
CN112100794A (zh) | 一种基于电流密度展开的磁共振线圈去耦合设计方法 | |
Zhang et al. | A spiral, bi-planar gradient coil design for open magnetic resonance imaging | |
CN205194456U (zh) | 一种用于关节成像的磁共振超导磁体 | |
CN114519289A (zh) | 用于磁共振成像的锥形阵列梯度线圈及其设计方法和应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |