Abstract
The present application provides a method of designing a high shielding gradient coil for a planar superconducting magnetic resonance imaging (MRI) system and a gradient coil thereof, the method determines a shielding area according to an outer profile of a metal conductor around the position of the gradient coil in the planar superconducting MRI system, and performs partitioned shielding of a stray field. The constraint values of stray fields at different partitioned zones of the shielding area are adjusted according to the shielding requirements. The primary coils of both the transverse gradient coil and the longitudinal gradient coil optimized by the design method of the high shielding gradient coil contain a reverse coil, which generates a magnetic field that offsets leakage magnetic field of other coils, thus achieving the purpose of reducing the stray field of the gradient coil.
Claims
1. A method of designing a high shielding gradient coil for a planar superconducting magnetic resonance imaging (MRI) system, wherein, the method of designing the high shielding gradient coil comprises the following steps: Step 1, determining a shielding area to be an outer profile of a conductor structure on a magnet close to a gradient coil; Step 2, calculating a stray field distribution of a non-shielding gradient coil in the shielding area; Step 3, partitioning the shielding area according to the stray field distribution of the non-shielding gradient coil obtained in Step 2, and determining a stray field strength constraint at each partitioned zone, wherein, for a partitioned zone S.sub.i of the shielding area, an absolute value of the stray field strength does not exceed the determined stray field strength constraint value B.sub.si, wherein i is a serial number of the partitioned zone of the shielding area; Step 4, setting the stray field strength constraint determined according to each partitioned zone of the shielding area into a constraint condition of an optimization function, namely:
−B.sub.s.sub.i≤s.sub.Bz.sub.i≤B.sub.s.sub.i (1) wherein, s.sub.Bz.sub.i is the stray field strength in the partitioned zone i of the shielding area; Step 5, determining a target function for optimization design of the gradient coil, and setting a control deviation of a target magnetic field of an imaging area into the constraint condition of the optimization function; and establishing the optimization function according to the target function, the control deviation constraint of the target magnetic field of the imaging area and the stray field strength constraint in Step 4; Step 6, establishing the optimization function according to the target function determined in Step 5, and executing the optimization design to minimize the target function and simultaneously the magnetic field conditions meet the requirements for the stray field strength and the control deviation of the target magnetic field of the imaging area; obtaining a distribution of a stream function on the gradient coil according to a calculation result of a stream function X, and obtaining a structure and parameters of the gradient coil by a method of stream function equipotential line; and judging whether the structure and parameters of the gradient coil meet the design requirements, and if not, adjusting the stray field strength constraint and conducting re-optimization until the obtained structure and parameters of the gradient coil meet the design requirements.
2. The method of designing a high shielding gradient coil for a planar superconducting MRI system of claim 1, wherein, the target function for optimization design of the gradient coil in Step 5 is the following formula (2), and the deviation constraint of the target magnetic field of the gradient coil in the imaging area is the following formula (3), in the formula (2), a.sub.ij is an element in a magnetic field calculation matrix from the imaging area to a coil plane, b.sub.ij is an element in a magnetic field calculation matrix from the shielding area to the coil plane; x.sub.j is an element in a stream function vector at boundary element nodes when the gradient coil is designed using a boundary element method, t.sub.i is an element in a target magnetic field vector, X is the stream function vector at the boundary element nodes, and X.sup.T is a transposed vector of the stream function vector at the boundary element nodes, L is an inductance matrix, R is a resistance matrix, M.sub.1 is the number of sampling points in the imaging area, M.sub.2 is the number of sampling points in the shielding area, i and j are element position variables, w.sub.1 and w.sub.2 are weight factors, andfis the target function; and in the formula (3), A is a magnetic field calculation matrix from the imaging area to the coil plane, T.sub.Bz is a target magnetic field vector, and ε is a maximum deviation of magnetic field gradient linearity.
3. A high shielding gradient coil system for a planar superconducting MRI system, wherein, the gradient coil system comprises a transverse coil and a longitudinal coil, wherein the transverse coil includes an x gradient coil and a y gradient coil, and the longitudinal coil is a z gradient coil; a half of a single pole face of primary coils of the x gradient coil and they gradient coil contain both a forward coil and a reverse coil, wherein the reverse coil has a current running direction opposite to that of the forward coil; a single pole face of a primary coil of the longitudinal coil, which is the z gradient coil, contains both a forward coil and a reverse coil, wherein the reverse coil has a current running direction opposite to that of the forward coil; and the shielding coils of the transverse gradient coil and the longitudinal gradient coil only contain coils in the same direction.
4. The gradient coil system of claim 3, wherein at least the z gradient coil is a high shielding gradient coil, and both the x gradient coil and the y gradient coil are high shielding gradient coils, or ordinary shielding gradient coils or non-shielding gradient coils; the stray field shielding area of the high shielding gradient coil includes both a planar area in a magnet groove and a side area, so that the outer profile of the conductor structure on the magnet close to the gradient coil are fully shielding from the stray field; the ordinary shielding gradient coil refers to that the stray field shielding area of the gradient coil only comprises the planar area in the magnet groove; the non-shielding gradient coil refers to a gradient coil containing a primary coil and no shielding coil; wherein the primary coil is located in an inner layer and the shielding coil is located in an outer layer, wherein x primary coil, y primary coil and z primary coil in the inner layer are not required to be arranged in a certain order, and x shielding coil, y shielding coil or z shielding coil in the outer layer are also not required to be arranged in a certain order.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0019] FIG. 1 is a schematic structural diagram of a planar superconducting magnet and gradient device according to the present application, in which 1 denotes a planar superconducting magnet, 2 denotes a planar gradient device with a gradient coil arranged therein, 3 denotes an anti-eddy current plate, 4 denotes a steel band ring, and 5 denotes a shimming ring;
[0020] FIG. 2 shows the position arrangement of gradient coils in the planar gradient assembly;
[0021] FIG. 3 shows the design steps of the high shielding gradient coil in the present application;
[0022] FIG. 4a is a schematic diagram of partitioning of the shielding area according to an outer profile of a conductor structure around the gradient coil and the corresponding stray field strength constraint when designing the high shielding gradient coil;
[0023] FIG. 4b is a schematic diagram of the stray field strength constraint only considering the planar area in a magnet groove when designing an ordinary shielding gradient coil;
[0024] FIG. 4c is a schematic diagram showing the case comprising only a primary coil without any shielding coil when designing a non-shielding gradient coil;
[0025] FIG. 5 is a schematic diagram of mesh generation, target point sampling and shielding area sampling when designing high shielding planar gradient coils by a boundary element method;
[0026] FIG. 6a is a two-dimensional structure diagram of an x primary coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application; FIG. 6b is a two-dimensional structure diagram of an x shielding coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application; FIG. 6c is a three-dimensional structure diagram of an x primary coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application; and FIG. 6d is a three-dimensional structure diagram of an x shielding coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application.
[0027] FIG. 7a is a two-dimensional structure diagram of a z primary coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application; FIG. 7b is a two-dimensional structure diagram of a z shielding coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application; FIG. 7c is a three-dimensional structure diagram of a z primary coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application; and FIG. 7d is a three-dimensional structure diagram of a z shielding coil of a gradient coil for a planar superconducting MRI system obtained according to the design method of the present application.
DETAILED DESCRIPTION
[0028] The present application will be further explained with reference to the accompanying drawings and specific embodiments.
[0029] As shown in FIG. 1, the high shielding gradient coil and its surrounding components for the planar superconducting MRI system of the present application include a planar superconducting magnet 1, a planar gradient assembly 2, an anti-eddy current plate 3, a steel band ring 4 and a shimming ring 5. The anti-eddy current plate 3 is laid to be closely fitted to the planar superconducting magnet 1, the steel band ring 4 is located between the anti-eddy current plate 3 and the shimming ring 5, the planar gradient assembly 2 is located in a groove surrounded by the anti-eddy current plate 3, the steel band ring 4 and the shimming ring 5, and the gradient coil is placed in the planar gradient assembly 2. In FIG. 4a and FIG. 4b, the area with a stray field constraint strength of B.sub.s.sub.1 corresponds to an outer profile of the anti-eddy current plate. In FIG. 4a, the areas with stray field constraint strengths of B.sub.s.sub.2 and B.sub.s.sub.3 correspond to an outer profile of the steel band ring 4, and the area with a stray field constraint strength of B.sub.s.sub.4 corresponds to an outer profile of the shimming ring 5.
[0030] According to actual demand, the gradient coil of the present application can include the following three modes of combination:
[0031] mode (1), including a high shielding x gradient coil, a high shielding y gradient coil and a high shielding z gradient coil;
[0032] mode (2), including an ordinary shielding x gradient coil, an ordinary shielding y gradient coil and a high shielding z gradient coil;
[0033] mode (3), including a non-shielding x gradient coil, a non-shielding y gradient coil and a high shielding z gradient coil.
[0034] Among them, the ordinary shielding gradient coils and non-shielding gradient coils can be designed according to existing methods. For example, in the structure of the planar superconducting MRI system of FIG. 1, if the anti-eddy current plate can effectively resist the eddy current generated by the transverse gradient coil, the combination mode (3) of the gradient coil can be adopted; if there is no space for installing the anti-eddy current plate in the planar superconducting MRI system, which is limited by practical engineering, the transverse ordinary shielding gradient coil with low complexity can be adopted, and the combination mode (2) of the gradient coil can be adopted; if there is no space for installing the anti-eddy current plate in the planar superconducting MRI system, and the transverse and longitudinal gradient coils both have very serious magnetic leakage on the lateral sides thereof, the combination mode (1) of the gradient coil can be adopted.
[0035] The positional relationship of gradient coils in the present application is shown in FIG. 2: the primary coil is located in the inner layer and the shielding coil is located in the outer layer, wherein the x primary coil, the y primary coil and the z primary coil in the inner layer are not required to be arranged in a certain order, and the x shielding coil, the y shielding coil or the z shielding coil in the outer layer are also not required to be arranged in a certain order.
[0036] As shown in FIG. 3, according to an embodiment of the present application, the high shielding gradient coil is designed by the following steps:
[0037] Step 1, determining a shielding area to be an outer profile of a conductor structure around the gradient coil. As shown in FIG. 1, the anti-eddy current plate 3, the steel band ring 4 and the shimming ring 5 are conductor structures around the gradient coil, and the corresponding outer profile of the conductor structures is shown as the shielding area in FIG. 4a.
[0038] Step 2, calculating a stray field distribution of a non-shielding gradient coil in the above shielding area, that is, designing a non-shielding gradient coil that meets the requirements for a magnetic field of the imaging area by adopting the design method shown in FIG. 4c which only comprises a primary coil and imposes no constraint on the shielding area, and calculating the stray field of the non-shielding gradient coil in the shielding area shown in FIG. 4a;
[0039] Step 3, partitioning the shielding area according to the stray field distribution of the non-shielding gradient coil, and determining a stray field strength constraint at each partitioned zone, wherein, for a partitioned zone S.sub.i of the shielding area, an absolute value of the corresponding stray field strength does not exceed B.sub.si, wherein i is a serial number of the partitioned zone of the shielding area;
[0040] Step 4, setting the stray field strength constraint determined according to each partitioned zone of the shielding area into a constraint condition of an optimization function. In the numerical solution of the gradient coil, it is necessary to sample the shielding points in the shielding area. Taking FIG. 4a as an example, the magnetic field values of the shielding points need to meet the following requirements:
|s.sub.Bz.sub.1|≤B.sub.s.sub.1
|s.sub.Bz.sub.2|≤B.sub.s.sub.2
|s.sub.Bz.sub.3|≤B.sub.s.sub.3
|s.sub.Bz.sub.4|≤B.sub.s.sub.4 (4)
[0041] wherein s.sub.Bz.sub.1, s.sub.Bz.sub.2, s.sub.Bz.sub.3 and S.sub.Bz.sub.4 are stray field values at corresponding sampling points of the zones of the shielding area partitioned according to the outer profile of the conductor structure around the gradient coil, where the stray field values are obtained according to current density of boundary elements, and B.sub.s.sub.1, B.sub.s.sub.2, B.sub.s.sub.3 and B.sub.s.sub.4 are stray field strength constraint values corresponding to the partitioned zones of the shielding area.
[0042] Step 5, determining a target function for optimization design of the gradient coil, and setting a control deviation of a target magnetic field of an imaging area into the constraint condition of the optimization function. FIG. 5 is a schematic diagram of a numerical calculation method for designing high shielding planar gradient coils by a boundary element method. When solving the gradient coil numerically, it is necessary to sample the target points in the imaging area, and the magnetic field values at the target points need to meet the following requirements:
|t.sub.Bz−T.sub.Bz|≤ε|T.sub.Bz| (5)
[0043] wherein t.sub.Bz is the magnetic field value of the target point obtained according to current density of boundary elements, T.sub.Bz is a pre-set target magnetic field, and cε is the maximum deviation of the gradient linearity of the magnetic field.
[0044] Step 6, executing the optimization design to obtain the structure and parameter information of the gradient coil, and judging whether the information meets the requirements for configuration, if not, adjusting the stray field strength constraint, and performing optimization design until the obtained structure and parameter information of the gradient coil meets the design requirements.
[0045] A gradient coil system designed according to the design method of the present application comprises an x gradient coil, a y gradient coil and a z gradient coil, wherein at least the z gradient coil is a high shielding gradient coil, and both the x gradient coil and they gradient coil are high shielding gradient coils, or ordinary shielding gradient coils or non-shielding gradient coils.
[0046] According to Step 1 to Step 3 of the design method of the present application, a high shielding gradient coil that meets the requirements can be designed by adjusting the stray field strength constraint value. The stray field strength constraint value can be set in consideration of the stray field distribution of the non-shielding gradient coils, and the structures of high shielding gradient coils designed based on different stray field strength constraint values are varied.
[0047] The following Table 1 shows the dimensions of a gradient coil designed according to the design method of the present application, for a 0.7 T planar superconducting MRI system. The target gradient strength of the gradient coil is 15 mT/m, and the deviation of target magnetic field in the imaging area is less than 5%. The constraint values of stray field in shielding area are shown in Table 2, and the corresponding shielding area dimensions are shown in FIG. 4a, FIG. 4b and FIG. 4c. A high shielding x gradient coil and a high shielding z gradient coil are designed respectively according to the constraint parameters of the gradient coils. Since the coil structure of the high shielding y gradient coil is similar to that of the high shielding x gradient coil, except for a spatial phase difference of 90 degrees, the structure of the high shielding y gradient coil will not be described repeatedly herein.
TABLE-US-00001 TABLE 1 Design dimensions of a high shielding gradient coil Gradient coil Radius (r/m) Height (z/m) X gradient coil main coil 0.480 ±0.238 shielding coil 0.500 ±0.283 Y gradient coil main coil 0.480 ±0.242 shielding coil 0.500 ±0.287 Z gradient coil main coil 0.480 ±0.246 shielding coil 0.500 ±0.291
TABLE-US-00002 TABLE 2 Constraint values of the stray field in the shielding area x gradient coil z gradient coil B.sub.s1 (Gauss) 3.0 3.0 B.sub.s2 (Gauss) 4.0 3.5 B.sub.s3 (Gauss) 5.0 4.0 B.sub.s4 (Gauss) 5.0 5.0
[0048] FIGS. 6a, 6b, 6c and 6d show the structure of the high shielding x gradient coil, in which FIG. 6a shows a two-dimensional structure of the primary coil, FIG. 6b shows a two-dimensional structure of the shielding coil, FIG. 6c shows a three-dimensional structure of the primary coil and FIG. 6d shows a three-dimensional structure of the shielding coil. The structures of the primary coil in FIG. 6a and FIG. 6c comprise a reverse coil, which generates a magnetic field that offsets the leakage magnetic field of other coils, thereby eventually reducing the stray field of the high shielding x gradient coil.
[0049] FIGS. 7a, 7b, 7c and 7d show the structure of the high shielding z gradient coil, in which FIG. 7a shows a two-dimensional structure of the primary coil, FIG. 7b shows a two-dimensional structure of the shielding coil, FIG. 7c shows a three-dimensional structure of the primary coil and FIG. 7d shows a three-dimensional structure of the shielding coil. Similar to the high shielding x gradient coil, the structures of the primary coil of the high shielding z gradient coil in FIGS. 7a and 7c also comprise a reverse coil, which generates a magnetic field that offsets the leakage magnetic field of the other coils, thereby eventually reducing the stray field of the high shielding z gradient coil.
