WO2013166810A1 - Système d'aimant supraconducteur d'imagerie à résonance magnétique et procédé et dispositif pour acquérir des paramètres structuraux de celui-ci - Google Patents
Système d'aimant supraconducteur d'imagerie à résonance magnétique et procédé et dispositif pour acquérir des paramètres structuraux de celui-ci Download PDFInfo
- Publication number
- WO2013166810A1 WO2013166810A1 PCT/CN2012/083287 CN2012083287W WO2013166810A1 WO 2013166810 A1 WO2013166810 A1 WO 2013166810A1 CN 2012083287 W CN2012083287 W CN 2012083287W WO 2013166810 A1 WO2013166810 A1 WO 2013166810A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- magnetic field
- current
- coil
- coil system
- value
- Prior art date
Links
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000009826 distribution Methods 0.000 claims abstract description 78
- 238000005457 optimization Methods 0.000 claims abstract description 43
- 230000005291 magnetic effect Effects 0.000 claims description 306
- 238000013461 design Methods 0.000 claims description 41
- 238000010586 diagram Methods 0.000 claims description 40
- 238000003384 imaging method Methods 0.000 claims description 23
- 230000000670 limiting effect Effects 0.000 claims description 11
- 230000002441 reversible effect Effects 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 12
- 230000005672 electromagnetic field Effects 0.000 description 10
- 239000010410 layer Substances 0.000 description 10
- 238000004519 manufacturing process Methods 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 5
- 239000000306 component Substances 0.000 description 4
- 238000004146 energy storage Methods 0.000 description 4
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 201000004569 Blindness Diseases 0.000 description 2
- 206010009244 Claustrophobia Diseases 0.000 description 2
- 239000003302 ferromagnetic material Substances 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 208000019899 phobic disease Diseases 0.000 description 2
- 241000237858 Gastropoda Species 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/381—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using electromagnets
- G01R33/3815—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using electromagnets with superconducting coils, e.g. power supply therefor
Definitions
- the magnetic resonance imaging technology of the present invention especially relates to a magnetic resonance imaging superconducting magnet system and a method and a device for acquiring the structural parameters thereof.
- Magnetic Resonance Imaging is a high-tech image based on the magnetic properties of biological magnetic nuclei, such as hydrogen atoms, in magnetic fields. It has no electromagnetic radiation, high image contrast, and can be imaged in any direction. The advantages make it the first in medical imaging equipment.
- the core components of the RI system mainly include magnet systems (also known as: superconducting magnet systems), spectrometer systems, computer systems, and image display systems.
- the magnet system also known as the MRI superconducting magnet system
- the function of the magnet system is to produce a high field strength, high uniformity, and high stability spatial magnetic field distribution in a generally spherical imaging region, so that the hydrogen atoms inside the human body are magnetized in the magnetic field, and the RF is generated by the RF line.
- the signal causes the hydrogen atom to resonate and attract the RF excitation, so that the energy absorbed by the hydrogen atom is released to be captured by the signal receiving device. Capture, finally through the image processing after the image processing, the resulting image resolution is proportional to the magnetic field strength.
- Hpp peak-to-peak homogeneity
- the ratio, ie: Hpp (B max - B min ) / mean (B) xl0 6 , in parts per million, ppm, where B max represents the highest value of the magnetic field strength within the DSV, 8 criz ⁇ indicates the lowest value of the magnetic field strength in the DSV, mean(B) represents the average value of the magnetic field strength in the DSV, and the smaller the value of the peak-to-peak uniformity of the magnetic field, the better the peak-to-peak uniformity of the magnetic field.
- the DSV requires a high uniformity magnetic field distribution with a peak-to-peak uniformity of the magnetic field better than lOppni in an imaging spherical region with a diameter of 40 to 50 cm.
- the magnet system In the early MRI system, the magnet system was more than two meters long, and the patient performed magnetic resonance detection. Out of tension and worry, this phenomenon is called claustrophoscopic clinically. In order to reduce the phenomenon of claustrophobia, MRI superconducting magnet system design has been moving toward how to reduce the length of the magnet system without reducing the size and magnetic of the imaging area. In the direction of uniformity, for example, in a 1.5T MRI system in 1989, the magnet system was up to 2.4m in length and weighed 13t. In 2009, the same 1.5T MRI system, the length of the magnet system was reduced to 1.37m, and the weight was 3.2t. However, if the magnet system is too short, it will bring other problems, such as difficulty in achieving uniformity, large electromagnetic stress, and difficult construction due to compact structure. Therefore, the MRI superconducting magnet system is designed. The process of balancing multiple parameters.
- the design parameters of the MRI superconducting magnet system mainly include: 1 the spatial size of the line pattern to be arranged; 2 the size, shape, central magnetic field strength and peak-to-peak uniformity of the magnetic field; 3 screens Masking technique and 5 Gaussian stray field range; the highest value of magnetic field strength and current safety margin in 4-line graph. These parameters become the main indicators determining the performance of MRI superconducting magnets. among them:
- Space size to be laid out Generally, it is a hollow core solenoid space with a rectangular cross section.
- the size of the hollow core solenoid space is mainly determined by the inner diameter miograph, outer diameter and length of the hollow core solenoid.
- the inner diameter of the space determines the minimum inner diameter of the magnet system, determines the room temperature aperture of the final magnet system and the comfort level during the patient's diagnosis;
- the outer diameter of the space determines the outer diameter of the magnet system, outside the magnet system.
- the size of the path determines the overall system cost and floor space; and the length of the space determines the final length of the magnet system, which is the decisive factor in reducing the phenomenon of claustrophobia.
- Imaging area size, shape, central magnetic field strength and peak-to-peak value of the magnetic field is generally a spherical body with a diameter of 40cm ⁇ 50cm, which can meet the needs of whole body imaging; the central magnetic field strength determines the resolution of imaging, 1.5T and 3T is the mainstream product in the current market.
- the peak-to-peak uniformity of the magnetic field is generally better than lOppm to meet the imaging needs.
- Shielding technology is generally divided into passive shielding and active shielding technology.
- Magnetic shielding is used to shield the influence of the magnetic field generated by MRI on the surrounding environment.
- the magnetic field coupling of ferromagnetic substances around the magnetization affects the MRI magnetic field uniformity.
- Passive shielding technology is to form a magnetic field loop around the MR1 through a ferromagnetic material, and to control the magnetic field outside the ferromagnetic material within a certain range. This technology makes the MRI large-area, inconvenient to install and susceptible to environmental factors such as temperature.
- the structure of the magnet is relatively simple and the cost of the magnet is greatly reduced.
- Many early products shield the magnetic field in this way.
- one or more pairs of reverse current coils can be added to the outer layer of the magnet system by active shielding technology to shield the stray field.
- This design makes the magnet system footprint greatly reduced.
- the turns structure is more complicated than the passive shield magnet system.
- Stray fields generally require a magnetic field strength of less than 5 Gauss in an ellipsoid.
- the highest value of the magnetic field strength and the current safety margin in the 4-wire When the coil is energized, a magnetic field distribution is generated in the space, and the energized coil generates an electromagnetic force in the magnetic field. The larger the magnetic field is the electromagnetic force received by the coil. The greater the force, the greater the electromagnetic force will weaken the superconducting wire Performance, even destroying superconducting magnets. Therefore, in order to avoid excessive electromagnetic force on the coil, the maximum value of the magnetic field strength in the coil is usually limited to less than 8T.
- the working point of the coil including the running current (lop) of the coil and the highest value of the magnetic field strength (B max ) in the coil, combined with the critical characteristics of the selected superconducting wire, calculates the critical point corresponding to the working point (Ic, Bc), calculate the current safety margin of the magnet system by the ratio of the operating current (lop) at the operating point to the operating current (Ic) at the critical point.
- the current safety margin should be set to less than 80%.
- the difficulty in designing an MRI superconducting magnet system is to design a magnetic field distribution that produces a high magnetic field and high uniformity in a small space.
- the design parameters of the MRI superconducting magnet system are contradictory to each other. The shorter the length of the magnet, the more difficult it is to achieve the image area and the more complex the magnet structure. Therefore, for the magnet system designer to balance the advantages and disadvantages, select the appropriate design parameters, and design the MRI superconducting magnet system to meet the requirements.
- the design of the magnet system is essentially a solution to the electromagnetic field problem.
- the solution of the electromagnetic field problem can be divided into two categories, one is the positive problem of the electromagnetic field, that is: the electromagnetic field distribution of the space is calculated according to the current source distribution in the space; the other is the inverse problem of the electromagnetic field, ie: according to the space
- the required electromagnetic field distribution is used to solve the current source position of the space.
- magnet system designers used the method of solving the positive problem of electromagnetic field, and through the distribution characteristics of the magnetic field generated by the solenoid coil, patching multiple turns in space to improve the uniformity of the magnetic field generated by the magnet system, the workload is large. It is difficult to meet the requirements.
- designers continue to solve the electromagnetic field inverse problem through a variety of numerical optimization algorithms, such as Monte Carlo, simulated annealing algorithm, genetic algorithm, etc.
- the numerical optimization algorithm can be divided into global Optimization algorithm and local optimization algorithm.
- the global optimization algorithm can solve the global optimal solution in the whole space, but due to the slow calculation speed, it is difficult for the magnet designer to modify the program effectively and effectively, so that the magnet design efficiency is low.
- the local optimization algorithm can optimize the parameters of the coil according to the initial values of the parameters. However, the selection of the initial values determines whether the global optimal solution can be solved. For many years, magnet system designers have been trying to solve a reasonable initial value and combine a local optimization algorithm to design a more reasonable magnet system, so that the entire design process of the magnet system has higher computational efficiency and can achieve a global optimal solution.
- Dr. Huawei Zhao proposed a multilayer current density design method for designing a compact MI magnet system. He proposed dividing the space in the line to be divided into multiple layers. Structure, each source layer is divided into several source points, each source point represents an ideal current ring, and the current distribution curve at all source point positions satisfying the requirements of the magnet system is solved by the regularization numerical method, and the current distribution is obtained. The peak position of the curve is used as the initial position of the spiral pipeline, and combined with the nonlinear optimization algorithm, the magnet system that meets the system design requirements is finally solved.
- the design method can obtain the initial value of the position of the coil by a layered method, and then design a magnetic system that meets the requirements by a local optimization algorithm.
- the design of the multilayer current density has at least the following problems:
- the stratified position selection has blindness, making it difficult to ensure that the line graph structure is a global optimal solution; a magnet system design scheme needs to run for 4 to 10 hours, and the design efficiency is relatively high. low.
- Dr. Quang M. Tien proposed a global optimal design method based on the minimum energy storage of magnet systems.
- each grid represents a spiral line ⁇
- the current distribution map in all the grids is calculated under the condition that the magnet system requirements are met and the magnet system has the minimum energy storage; also in the current distribution At the peak position, the initial position of the solenoid coil is reasonably arranged, and the sequence quadratic programming algorithm is combined to solve the magnet system that meets the requirements.
- the design method is that the global optimum initial value is obtained in the entire layout space, and the coil structure is more reasonable.
- the global optimal design method based on the minimum energy storage of the magnet system has at least the following problems:
- the grid current is a continuously changing curved surface, and it is difficult to define each spiral pipeline. In the initial position, it is difficult to achieve an optimal design. Summary of the invention
- One of the technical problems to be solved by the embodiments of the present invention is to provide a magnetic resonance imaging superconducting magnet system and a method and device for acquiring the structural parameters thereof, so as to improve the acquisition efficiency of the structural parameters of the MRI superconducting magnet system, and
- the coil structure in the MRI superconducting magnet system manufactured by the structural parameters is small, so that the MRI superconducting magnet system manufactured based on the structural parameters has higher performance.
- a method for acquiring structural parameters of a magnetic resonance imaging superconducting magnet system includes:
- each grid obtained by the division is regarded as a current ring, and the current value of each current ring includes a positive value, a negative value or zero; And uniformly dividing the surface of the spherical imaging region of the superconducting magnet system and the surface of the ellipsoid of the 5 Gaussian stray field into a plurality of target points according to a preset standard;
- the magnetic field + value uniformity at each target point on the surface of the spherical imaging region is not greater than the first predetermined magnetic field peak-to-peak uniformity, and the limiting condition of the magnetic field strength value at each target point on the surface of the ellipsoid of the 5 Gaussian stray field is a constrained component, calculating a current distribution at each of the meshes that minimizes the volume of the coil system under a given operating current density condition by a linear programming algorithm, and obtaining a current distribution map;
- the first predetermined magnetic field value uniformity is greater than or Equal to the second preset magnetic field value uniformity in the design requirements of the superconducting magnet system;
- the number and spatial position of the non-zero current values in the current distribution diagram are respectively taken as the number and initial position of the spiral pipeline diagram in the coil system;
- the limiting condition and the current safety margin limit of the operation of the coil system are the second constraint condition.
- the nonlinear optimization algorithm is used to calculate the size and position parameters of each solenoid coil when the volume of the coil system is the smallest, and the final structural parameters of the coil system are obtained.
- the size and positional parameters of each solenoid coil include an inner radius, an outer radius, and an axial position of each end of each spiral line.
- the spatial extent of the coil system to be arranged is specifically a solenoid-shaped region having a rectangular cross-section, the extent of the region being determined by parameters of a rectangular cross-section, the parameters of the rectangular cross-section including the inner diameter of the rectangular cross-section , outer diameter and length;
- Continuous meshing of the spatial extent of the coil system to be arranged includes: successively meshing the radial cross sections in the radial and axial directions, respectively.
- the continuous meshing of the rectangular cross section in the radial direction and the axial direction respectively comprises: performing a two-dimensional continuous mesh division on the rectangular cross section in the radial direction and the axial direction, respectively
- the radial and axial directions are divided into several equal parts to form a two-dimensional continuous space grid.
- the current profile gives the magnitude and direction of the current at each grid, and non-zero current values at the grid are brought together to form a clear boundary non-zero current cluster.
- a forward current spiral pipeline map is placed at each positive position of the non-zero current cluster, and a reverse current spiral pipeline map is respectively placed at each negative position of the non-zero current cluster.
- the limiting conditions of the magnetic field strength values at the respective target points on the surface of the ellipsoid of the 5 Gaussian stray field include: 5
- the magnetic field strength value at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss.
- the second constraint comprises:
- the peak-to-peak uniformity of the magnetic field of the axial magnetic field at each target point on the surface of the spherical imaging region is not greater than the peak-to-peak uniformity of the second predetermined magnetic field;
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not greater than the preset magnetic field strength
- the current safety margin of the coil system operation is not greater than a preset current safety margin calculated from a highest magnetic field strength value in the coil system and a critical characteristic of the superconducting wire selected by the coil system;
- the spacing between adjacent solenoid coils in the turns system is greater than the preset spacing
- each of the solenoid lines ⁇ causes the line drawing system to be within the spatial extent of the line system to be placed.
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the 5 Gaussian stray field is not greater than the preset magnetic field strength including:
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss.
- the method further includes:
- the current safety margin of the coil system operation is obtained by the maximum magnetic field strength value in the preset coil system and the critical characteristics of the superconducting wire selected by the coil system.
- a receiving unit configured to receive spatial range information of the coil system to be arranged, a first predetermined magnetic field peak-to-peak uniformity, a limit condition of a magnetic field strength value at each target point on an ellipsoid surface of a 5 Gaussian stray field, and a running current density condition , the second preset magnetic field peak-to-peak uniformity limit condition in the design requirements of the superconducting magnet system, and the target points on the surface of the ellipsoid of the 5 Gaussian stray field
- each divided grid is regarded as a current ring, and the current value of each current ring includes a positive value, a negative value or zero;
- a first calculating unit configured to: the magnetic field peak-to-peak uniformity at each target point on the surface of the spherical imaging region is not greater than the first predetermined magnetic field peak-to-peak uniformity, and the magnetic field at each target point on the surface of the ellipsoid of the 5 Gaussian stray field
- the limiting condition of the intensity value is a first constraint component, and a current distribution map is obtained by a linear programming algorithm for calculating a current distribution at each grid at a given operating current density condition to minimize the volume of the coil system;
- a second calculating unit configured to respectively use the number and spatial position of the non-zero current values in the current distribution map as the number and initial position of the solenoid line ⁇ in the coil system; and the second predetermined magnetic field peak-to-peak uniformity Restricted condition, magnetic field strength value limit at each target point on the surface of the ellipsoid of the 5 Gaussian stray field, maximum magnetic field strength value limit condition in the coil system, current safety margin limit of the operation of the coil system is the second constraint
- the size and position parameters of each solenoid coil are calculated by a nonlinear optimization algorithm to obtain the final structural parameters of the coil system; wherein the size and position parameters of each spiral line include each spiral The inner radius, outer radius, and axial position of the two ends of the line.
- the spatial extent of the coil system to be arranged is specifically a solenoid-shaped region having a rectangular cross section, the range of the region being determined by parameters of a rectangular cross section, and the parameters of the rectangular cross section include a rectangular cross section.
- the rectangular cross section is divided into two dimensions in the radial direction and the axial direction, respectively, and the rectangular cross section is along the radial direction.
- the axial directions are divided into several equal parts to form a two-dimensional continuous space grid.
- the current distribution diagram gives the magnitude and direction of the current at each grid, and the non-zero current values at the grid are gathered together to form a clear-cut non-zero current cluster;
- Each positive position of the non-zero current cluster is used to arrange a forward current solenoid line ⁇ , and each negative value position is used to arrange a reverse current solenoid line ⁇ , respectively.
- the second constraint comprises:
- the peak-to-peak uniformity of the magnetic field of the axial magnetic field at each target point on the surface of the spherical imaging region is not greater than the uniformity of the second predetermined magnetic field value
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not greater than the preset magnetic field strength
- the current safety margin of the coil system operation is not greater than a preset current safety margin calculated from the highest magnetic field strength value in the line surrounding system and the critical characteristic of the superconducting wire selected by the coil system. ;
- the spacing between adjacent spiral conduits in the coil system is greater than a predetermined spacing
- each solenoid line ⁇ causes the coil system to be within the spatial extent of the coil system to be placed.
- the limitation of the magnetic field strength value at each target point on the surface of the ellipsoid of the 5 Gaussian stray field comprises: 5 The magnetic field strength value at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss;
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not greater than the preset magnetic field strength including: 5
- the axial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field The superimposed magnetic field strength with the radial magnetic field is not more than 5 Gauss.
- the second calculating unit is further configured to obtain a current safety margin of the coil system operation according to a highest magnetic field strength value in the preset coil system and a critical characteristic of the superconducting wire selected by the coil system.
- a magnetic resonance imaging superconducting magnet system provided by an embodiment of the present invention is obtained by acquiring a structural parameter of the magnetic resonance system.
- the method and device for acquiring structural parameters of the MRI superconducting magnet system wherein the peak-to-peak uniformity of the magnetic field at each target point on the surface of the DSV is not greater than the first predetermined magnetic field peak-to-peak uniformity, and the ellipsoid of the 5 Gaussian stray field
- the constraint condition of the magnetic field strength value at each target point on the surface is the first constraint condition, and the condition of the given operating current density is calculated by the linear programming algorithm.
- the current distribution at each grid that minimizes the volume of the coil system is obtained, and the current distribution map is obtained, and then the number and spatial position of the non-zero current values in the current distribution map are respectively taken as the number of spiral pipeline turns in the coil system.
- the initial position is the second constraint condition in the design requirements of the superconducting magnet system.
- the nonlinear structure optimization parameter is used to calculate the size and position parameters of each spiral pipeline diagram when the volume of the coil system is the smallest, and the final structural parameters of the coil system are obtained.
- the hybrid design method combining the linear programming algorithm and the nonlinear optimization algorithm is adopted in the embodiment of the present invention, the number and initial position of the spiral pipeline ⁇ in the coil system are obtained by the linear programming algorithm with higher efficiency, and the nonlinear optimization algorithm is combined.
- the MRI superconducting magnet system with the structural parameters obtained in the above embodiment has a smaller wire loop structure, requires the least amount of superconducting wire, has the advantages of simple and compact wire loop structure, low cost, and is easy to construct and install, thereby making it based on
- the MRI superconducting magnet system manufactured by the structural parameters has high performance, and overcomes the problem that the existing MRI superconducting magnet system is inefficient in electromagnetic design and difficult to realize optimal design.
- the current safety margin limit provides a robust guarantee for the stable operation of the entire MRI superconducting magnet system.
- FIG. 1 is a schematic flow chart of a method for obtaining structural parameters of an MRI superconducting magnet system according to an embodiment of the present invention
- FIG. 2 is a schematic flow chart of a method for acquiring structural parameters of an MRI superconducting magnet system according to another embodiment of the present invention
- FIG. 3 is a structural parameter of an MRI superconducting magnet system according to an application embodiment of the present invention. Schematic flow chart of the acquisition method
- FIG. 4 is a schematic diagram of a spatial extent, a DSV region, and a 5 Gaussian stray field region in which the turns are arranged in the application embodiment shown in FIG. 3;
- FIG. 5 is a current distribution diagram obtained by a linear programming algorithm in the application embodiment shown in FIG. 3;
- FIG. 6 is a magnetic field «value uniformity distribution diagram of an axial magnetic field at a target point on the surface of the DSV obtained in the application example shown in FIG. 3;
- FIG. 7 is a magnetic field intensity distribution map at each target point on the surface of an ellipsoid of a 5 Gaussian stray field obtained in the application example shown in FIG. 3;
- FIG. 8 is a schematic diagram of initial position and size parameters of a non-zero current cluster discrete into a spiral line ⁇ in the application embodiment shown in FIG. 3;
- FIG. 9 is a schematic diagram of final structural parameters of a line ⁇ system optimized by a nonlinear optimization algorithm in the application embodiment shown in FIG. 3;
- FIG. 10 is a magnetic field uniformity distribution diagram of a DSV surface optimized by a nonlinear optimization algorithm in the application embodiment shown in FIG. 3;
- FIG. 11 is a magnetic field intensity distribution diagram at each target point on the surface of an ellipsoid of a 5 Gaussian stray field obtained by a nonlinear optimization algorithm in the application embodiment shown in FIG. 3;
- Figure 12 is a current safety margin diagram of the MRI superconducting magnet system in the application embodiment shown in Figure 3;
- Figure 13 is a diagram showing the distribution of magnetic field strength generated by the coil system in space in the application embodiment shown in Figure 3;
- Figure 14 is a magnetic field intensity distribution diagram of a solenoid enthalpy in one end of the inner layer in the coil system of the application embodiment shown in Figure 3;
- Figure 15 is a schematic block diagram of an apparatus for acquiring structural parameters of an MRI superconducting magnet system according to an embodiment of the present invention. detailed description
- FIG. 1 is a schematic flow chart showing a method of obtaining structural parameters of an MRI superconducting magnet system according to an embodiment of the present invention. As shown in FIG. 1, the method for obtaining the structural parameters of the MRI superconducting magnet system of this embodiment includes:
- each divided grid is regarded as a current ring, and the current value of each current ring may be positive, negative or zero;
- the preset standard evenly divides the DSV surface of the superconducting magnet system and the ellipsoidal surface of the 5 Gaussian stray field into several target points.
- the space range in which the coil system is to be arranged is preset according to the production requirements of the MRI superconducting magnet system.
- the first constraint condition is obtained by calculating a current distribution at each grid at a given operating current density to minimize the volume of the coil system by a linear programming algorithm.
- the first predetermined magnetic field peak-to-peak degree is greater than or equal to the design of the superconducting magnet system Finding the second preset magnetic field value uniformity.
- the limitation of the magnetic field strength value at each target point on the surface of the ellipsoid of the 5 Gaussian stray field may include: 5
- the magnetic field intensity value at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss.
- the nonlinear structure optimization parameter is used to calculate the size and position parameters of each spiral line ⁇ when the volume of the coil system is the smallest, and the final structural parameters of the coil system are obtained.
- each spiral line includes the inner radius, the outer radius of each solenoid coil, and the axial position of the two ends.
- a hybrid design method combining a linear programming algorithm and a nonlinear optimization algorithm is adopted, and the number and initial position of the spiral pipelines in the coil system are obtained by a linear programming algorithm with high efficiency, and then combined with a nonlinear optimization algorithm.
- the acquisition of the entire structural parameter takes only about 20 minutes; and the coil structure of the MRI superconducting magnet system using the structural parameters obtained by the above embodiment of the present invention is small, and the required amount of superconducting wire is minimum.
- the utility model has the advantages of simple compactness and low cost, and is easy to construct and install, so that the MRI superconducting magnet system manufactured based on the structural parameter has high performance and overcomes the existing MRI superconducting magnet system in electromagnetic field. Designed to be less efficient and difficult to achieve optimal design issues.
- the current safety margin limits provide a robust guarantee for the stable operation of the entire MRI superconducting magnet system.
- the second constraint may specifically include:
- the peak-to-peak uniformity of the magnetic field of the axial magnetic field at each target point on the surface of the spherical imaging region is not greater than the peak-to-peak uniformity of the second predetermined magnetic field; 5
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not greater than the preset magnetic field strength, for example, the axial direction at each target point on the surface of the ellipsoid of the 5 Gaussian stray field
- the superimposed magnetic field strength of the magnetic field and the radial magnetic field is not more than 5 Gauss;
- the current safety margin of the coil system operation is not greater than the preset current safety margin, and the preset current safety margin is the highest magnetic field strength value and coil in the coil system
- the critical characteristics of the superconducting wire selected by the system are calculated;
- the spacing between adjacent solenoid coils in the turns system is greater than the preset spacing
- each solenoid line ⁇ place the coil system within the space within which the coil system is to be placed.
- the spatial extent of the wire loop system to be disposed in the embodiment of the present invention may specifically be a solenoid-shaped region having a rectangular cross section, the range of the region being determined by parameters of the rectangular cross section, wherein the parameters of the rectangular cross section include The inner diameter, outer diameter, and length of the rectangular section.
- the continuous meshing of the spatial extent of the system to be arranged may be : Continuous meshing of the rectangular section in the radial direction and the axial direction respectively.
- continuous meshing of the rectangular section in the radial direction and the axial direction respectively includes: two-dimensional continuous radial and axial directions respectively for the rectangular section Meshing divides the rectangular section into several equal parts along the radial and axial directions to form a two-dimensional continuous space grid.
- the current distribution diagram gives the magnitude and direction of the current at each grid, and the current values at the grid are not zero. Together form a clear, non-zero current cluster.
- a forward current spiral line ⁇ is respectively arranged at each positive value position of the non-zero current cluster, and a reverse current spiral line ⁇ is respectively arranged on the Each negative position of a non-zero current cluster.
- the preset current safety margin can be calculated according to the highest magnetic field strength value in the preset coil system and the critical characteristic of the superconducting wire selected by the coil system.
- the spatial range of the coil system to be arranged is preset according to the manufacturing requirements of the MRI superconducting magnet system, and the spatial range is specifically a rectangular cross section.
- 2 is a schematic flow chart of a method for acquiring structural parameters of an MRI superconducting magnet system according to another embodiment of the present invention. As shown in FIG. 2, the method for acquiring structural parameters of the MRI superconducting magnet system of the embodiment includes:
- 201 Perform a two-dimensional continuous mesh division on the rectangular section of the coil system to be arranged in the radial direction and the axial direction respectively, and divide the rectangular section into a plurality of equal parts along the radial direction and the axial direction to form a two-dimensional continuous space grid.
- Each grid is regarded as a current ring, and the current value of each current ring can be positive, negative or zero; and the DSV surface of the superconducting magnet system and the ellipsoid of the 5 Gaussian stray field according to preset criteria The surface is evenly divided into several target points.
- the current value of the current ring at each mesh obtained by dividing the rectangular section is used as a variable to be optimized, so that the peak-to-peak uniformity of the magnetic field at each target point on the DSV surface is not greater than the first predetermined magnetic field peak-to-peak uniformity, 5
- the height of the magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss is the first constraint.
- the linear programming algorithm is used to calculate the meshes that minimize the volume of the coil system under the given operating current density.
- a current distribution diagram is obtained.
- the current distribution diagram shows the magnitude and direction of the current at each grid.
- the current value at most of the grid is zero, and only a few current values are not zero, not Zero current values are grouped together to form a clear, non-zero current cluster.
- the first predetermined magnetic field peak-to-peak uniformity is greater than or equal to the second predetermined magnetic field value uniformity in the superconducting magnet system design requirements.
- a forward current spiral line ⁇ is respectively arranged at each positive value position of the non-zero current cluster, and a reverse current solenoid coil is respectively arranged in the non-zero current coil.
- Each negative position of the zero current cluster the arrangement of the spiral line ⁇ conforms to the requirement that the actual magnet system consists of a plurality of separate spiral line diagrams of the same current magnitude; and, by the position and total current of the non-zero current cluster
- the size can be used to construct the initial position and size parameters of each spiral line.
- Initial position and size parameters of the solenoid line As the initial value of the spiral pipeline ⁇ structural parameters in the subsequent nonlinear optimization algorithm, it compensates for the blindness of the local optimization algorithm in the initial value selection of the solenoid coil structure parameters.
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss;
- the current safety margin of the coil system operation is not greater than the preset current safety margin to ensure the safety of the MRI superconducting magnet system
- the spacing between adjacent solenoid coils in the coil system is greater than the preset spacing, and the axial dimension of each solenoid coil may be constrained so that the spacing of adjacent spiral pipelines is greater than a preset spacing, thereby avoiding Each of the solenoid coils overlap each other to facilitate the construction of the solenoid line ;;
- each solenoid line ⁇ place the coil system within the space within which the coil system is to be placed.
- N is the number of solenoids included in the coil system, ie: the number of solenoid coils required for the MRI superconducting magnet system
- r inner is the ith spiral line of the N solenoid coils
- z lert is the axial position of one end of the i-th solenoid coil
- z right is the axial position of the other end of the i-th solenoid line.
- the embodiment shown in FIG. 2 of the present invention greatly improves the acquisition efficiency of the structural parameters of the MRI superconducting magnet system.
- the acquisition of the entire structural parameter takes only about 20 minutes, wherein operation 201 - 202 takes about 5 minutes, and operation 203 - 203 takes about 15 minutes.
- the MRI superconducting magnet system using the structural parameters obtained by the embodiment shown in Fig. 2 above requires the least amount of superconducting wire, which reduces the i-price of the entire MRI superconducting magnet system; meanwhile, the wire loop system has a simple structure and is easy to construct. And installation, the limitation of current safety margin provides an effective guarantee for the stable operation of the entire MRI superconducting magnet system.
- FIG. 3 it is a schematic flow chart of a method for acquiring structural parameters of an MRI superconducting magnet system according to an application embodiment of the present invention.
- the MRI superconducting magnet system is required to generate a magnetic field distribution with a central magnetic field strength of 1.5T in a DSV having a diameter of 50 cm, and the peak-to-peak uniformity of the magnetic field needs to be better than 10 ppm; 5 Gaussian stray field is constrained to an ellipsoid In the region, the long semi-axis length of the ellipsoid is 5m, and the short semi-axis is 4m; the inner radius of the coil system is not less than 0,40m, the outer radius is not more than 0.80m, and the length is shorter than 1.15m, namely: the coil system The minimum inner radius is 0.40m, the maximum outer radius of the coil system is 0.80m, and the maximum length of the coil system is 1.15m; the highest magnetic field strength in the coil system is less than
- the spatial extent of the coil system to be arranged is set according to the production requirements of the MRI superconducting magnet system.
- the two-dimensional continuous mesh division of the rectangular cross-section of the coil system to be arranged is respectively performed in the axial direction and the radial direction, and the rectangular cross-section is respectively divided into 80 and 40 in the axial direction and the radial direction.
- Each grid is regarded as an ideal current ring, and the current value of each current ring can be positive, negative or zero. .
- FIG. 4 it is a schematic diagram of a spatial range, a DSV surface area, and a 5 Gaussian stray field area in which the turns are to be arranged in the application embodiment of the present invention.
- Fig. 4 it is a schematic diagram of a spatial range, a DSV surface area, and a 5 Gaussian stray field area in which the turns are to be arranged in the application embodiment of the present invention.
- the upper part is a schematic diagram of a two-dimensional continuous space grid; the lower part is a schematic diagram of a coil system of a spiral line arranged in a subsequent operation, and the outermost rectangle represents a section of the coil system, each of which is small
- the rectangle represents the cross section of a spiral line.
- the outer ellipse represents the cross section of the ellipsoidal surface of the 5 Gaussian stray field, and the inner circle represents the section of the DSV surface.
- Gap represents the spacing between adjacent solenoid coils, specifically the inter-turn spacing.
- the current value of the current ring at each grid obtained by dividing the rectangular section is used as a variable to be optimized, and the current distribution of each grid at the smallest volume of the line graph system is calculated by a linear programming algorithm to obtain a current distribution map. .
- the matrix of the contribution matrix A zs and the radial magnetic field contribution matrix A rs , the magnetic field contribution matrix A zd , the magnetic field contribution matrix A zs and the magnetic field contribution matrix ⁇ are both 51x3200; the variable to be optimized is the current value I at all the grids, .
- the ideal current loop can be represented by a spiral line ⁇ with a certain cross section, establishing a mathematical model of the linear programming algorithm, since the mesh with the smallest volume of the coil system is obtained.
- Current distribution at all grids The total volume of ⁇ is used as the objective function, and the minimum value of the objective function takes the smallest volume of the line ⁇ system.
- the minimum value of the peak-to-peak uniformity of the axial magnetic field at each target point on the DSV surface is set to each target point on the surface of the ellipsoid of 20 ppm, 5 Gaussian stray field.
- the maximum values of the axial magnetic field and the radial magnetic field strength are set to 5 Gauss; assuming that the cross-sectional area of the superconducting wire used in the spiral line is 4.5312 m 2 , the critical current of the superconducting wire at 9 T background magnetic field strength is 950 A, Set the operating current density of the solenoid line ⁇ to 148 ⁇ / ⁇
- the radial position of the i-th grid is ri
- the cross-sectional area of the spiral line ⁇ represented by the ⁇ grid is Ai
- the volume ⁇ at the i-th grid is 2 ⁇ > ⁇ , assuming the same operating current density, the total volume of the coil system is:
- the first constraint includes: A... xl ⁇ 5 Gauss
- max(A zd xI) represents the maximum value of the magnetic field strength of the target point on the DSV surface
- min(A zd xI) represents the minimum value of the magnetic field strength of the target point on the DSV surface.
- B. The value of the center magnetic field strength is 1.5T; J is the operating current density of the spiral line ;; A mesh is the area occupied by each grid.
- a positive value indicates that a spiral line ⁇ with a forward current is required at the corresponding grid
- a negative value indicates that a solenoid coil with a reverse current is required at the corresponding grid, and zero corresponds to There is no need to arrange a spiral line around the grid. From the position of the non-zero current cluster of the current profile, it can be clearly seen that the spatial position of the six pairs of solenoids and the arrangement of the solenoid coils needs to be arranged.
- FIG. 6 the magnetic field peak-to-peak uniformity distribution map of the axial magnetic field at each target point on the surface of the DSV obtained corresponding to the current profile shown in Fig. 5 is shown. It can be seen from Fig. 6 that the magnetic field peak-to-peak uniformity of the axial magnetic field at each target point on the surface of the DSV satisfies 20 ppm.
- the intensity values are not greater than 5 Gauss.
- the initial position and size parameters of each spiral line ⁇ can be obtained and used as the initial values of the spiral line group structure parameters in the subsequent nonlinear optimization algorithm. As shown in FIG. 8, the initial position and size parameters of the non-zero current clusters in the current distribution diagram obtained by the linear programming algorithm in this application embodiment are separated into spiral pipelines.
- the solenoid line obtained in operation 303 is performed by a nonlinear optimization algorithm
- the initial values of the structural parameters of the circle are optimized, and the size and position parameters of each spiral line ⁇ are calculated when the volume of the coil system is the smallest, and the final structural parameters of the coil system are obtained.
- the variables to be optimized in the nonlinear optimization algorithm are the structural parameters of each spiral line, including the following dimensions and spatial position parameters of each solenoid coil: r inner (i), r outer (i), z left ( i), z right (i), where 1,2 ..,6, r inner (i) is the inner radius of the ith spiral line ,, r.
- Uter (i) is the outer radius of the ith spiral line
- z, eft (i) is the end axial position of the ith spiral line
- z right (i) is the ith solenoid
- the other end of the coil is axially positioned.
- the minimum value of the objective function takes the minimum volume of the system.
- the magnetic field strength of the axial magnetic field and the radial magnetic field of the target point on the surface of the ellipsoid of the 5 Gaussian stray field is set to be no more than 5 Gauss;
- the highest magnetic field strength in the coil system is less than 8T, and the current safety is less than 80% according to the highest magnetic field strength in the coil system and the performance of the selected superconducting wire; the same layer of spiral pipelines are convenient for construction and It is avoided that the solenoid coils overlap each other, exemplarily preset that the inter-turn spacing of the spiral line turns is greater than 1 cm, and the radial different interlayer spacing is greater than 1 cm.
- the objective function is: ⁇ 2 x [_ r ou, er ( 2 - r imer ⁇ if" x [z right (i) - ⁇ 1 (0)
- the second constraint includes:
- 8 3 ⁇ 4 ( ⁇ is the magnetic field strength of the axial magnetic field generated by the solenoid coil at each target point on the DSV surface, 8 ⁇ is the ellipse of the spiral line ⁇ in the 5 Gaussian stray field
- B rstray is the magnetic field strength of the radial magnetic field generated at each target point on the surface of the 5 Gaussian stray field ellipsoid
- max(B zdsv ) Indicates the maximum value of the magnetic field strength at all target points
- max(B zdsv ) represents the minimum of the magnetic field strength at all target points
- mean(B zdsv ) represents the average of the magnetic field strength at all target points
- lop is the superconducting magnet
- the operating current of the system at the operating point, Ic(B max ) is the corresponding critical current value Ic in the superconducting wire calculated from the operating current lop
- FIG. 9 it is a schematic diagram of the final structural parameters of the coil system obtained by the nonlinear optimization algorithm in the application embodiment shown in FIG.
- the black area indicates the solenoid coil in which the energization direction is forward
- the white area indicates the negative solenoid line ⁇ in the energization direction.
- the line drawing system has two layers of six pairs of spiral lines, wherein the inner layer has five pairs of spiral lines and the outer layer has a pair of shielded spiral lines.
- Fig. 10 is a view showing the distribution of the peak-to-peak uniformity of the magnetic field of the axial magnetic field at each target point on the DSV surface optimized by the nonlinear optimization algorithm in the application embodiment.
- 11 is a magnetic field intensity distribution diagram at each target point on the surface of an ellipsoid of a 5 Gaussian stray field obtained by a nonlinear optimization algorithm in the application embodiment shown in FIG. 3.
- FIG. 12 is a current safety margin diagram of the MRI superconducting magnet system in the application embodiment shown in FIG. 3, wherein the thick curve (Ic, Be) represents the critical current characteristic Ic-Bc curve of the selected superconducting wire, and the thin line (B) Max lop) represents the relationship between the operating current lop of the coil system at the operating point and the highest magnetic field strength value B max in the coil system. The intersection of the two curves is the critical point corresponding to the operating point of the coil system.
- Current safety margin for this application example when Ic 953.3A Degree is 70.35%,
- FIG. 13 is a magnetic field intensity distribution diagram of a coil system generated in a space in the application embodiment shown in FIG. 3, wherein the magnetic field distribution in the coil system has symmetry
- FIG. 14 is an inner layer in the coil system of the application embodiment shown in FIG.
- Table 1 shows the structural parameters of each solenoid coil in the MRI superconducting magnet system corresponding to the application embodiment shown in Fig. 3.
- Table 2 shows the structural parameters of the MRI superconducting magnet system finally obtained in the application example:
- all or part of the operations of implementing the above method embodiments may be performed by hardware associated with the program instructions, the foregoing program including machine readable instructions for implementing the method according to the present invention, which may be stored in a
- Machine readable instructions Thus, the invention also covers a storage medium storing a program for performing the method according to the invention. Previous storage media include: ROM, RAM, disk or optical disk, and other media that can store program code.
- the method for acquiring the structural parameters of the MRI superconducting magnet system of the present invention is implemented by software, hardware, firmware or any combination of software, hardware, and firmware, the method can be specifically implemented by:
- the production requirements parameters of the MRI superconducting magnet system for example, including but not Limited to central magnetic field strength, DSV size, magnetic field value uniformity requirement, 5 Gaussian stray field requirements, spatial extent of the coil system, operating current density of a given coil system, current safety margin;
- the loop iterative calculation is performed by the linear programming algorithm and the nonlinear optimization algorithm, and the structural parameters of the coil system satisfying the production requirements of the MRI superconducting magnet system are obtained.
- the process is divided into linear programming. Process and nonlinear optimization process two steps.
- the operation of 101 - 102 or 201 ⁇ 202 is performed, and the variable to be optimized in the linear programming algorithm is the current value at each mesh in the spatial range of the system to be arranged, and the objective function is the snail at all the grids.
- the total volume of the line pipeline is cyclically iteratively calculated by a linear programming algorithm.
- the linear programming algorithm automatically adjusts all variables to be optimized (ie, the current values at each grid) and Perform the next iteration calculation until the magnetic field at all target points satisfies the first constraint and the total volume of the spiral line ⁇ represented by all the grids is the smallest, stop the iterative calculation, and obtain the optimal solution of the variable to be optimized.
- the optimal solution of the optimization variables forms a current distribution map at each grid.
- the iterative calculation of the preset number of times can be performed based on the convergence characteristic of the total volume of the solenoid coil. For example, iterative calculation is performed 150 times, and the first constraint is satisfied and all the networks are selected.
- the current value of the smallest volume of the solenoid coil represented by the grid forms a current profile.
- the non-zero current cluster is converted into a spiral pipeline ⁇ , and the total current and spatial position of each non-zero current cluster are respectively corresponding.
- the initial total current and initial position of the solenoid coil, performing 103 or 203 ⁇ 204 operation, the variables to be optimized in the nonlinear optimization algorithm are the inner and outer radii of each solenoid coil and the axial directions of the two ends Position, the objective function is the total volume of the spiral line ⁇ at all grids.
- the loop iterative calculation is performed by the nonlinear optimization algorithm. If the second constraint is not satisfied, the nonlinear optimization algorithm automatically adjusts the size and position parameters of each spiral pipeline within the spatial range of the system to be arranged for the next iterative calculation. , until the conditions in the second constraint are satisfied and the total volume of all the solenoid coils is the smallest, outputting the inside and outside of each solenoid coil The radius and the axial position of the two ends provide the final structural parameters of the coil system.
- the iterative calculation of the preset number of times can be performed based on the convergence characteristic of the total volume of the spiral pipeline circle, for example, iterative calculation
- an MRI superconducting magnet system can be constructed to meet the production requirements.
- Figure 15 is a schematic block diagram of an apparatus for acquiring structural parameters of an MRI superconducting magnet system according to an embodiment of the present invention.
- the obtaining device of this embodiment can be used to implement the flow of the embodiment of the method for acquiring structural parameters of each of the above MRI superconducting magnet systems of the present invention. As shown in FIG. 15, it includes a receiving unit 301, a dividing unit 302, a first calculating unit 303, and a second calculating unit 304.
- the receiving unit 301 is configured to receive spatial range information of the coil system to be arranged, a first predetermined magnetic field peak-to-peak uniformity, a limitation condition of a magnetic field strength value at each target point on an ellipsoid surface of a 5 Gaussian stray field, and an operating current density
- the first predetermined magnetic field peak-to-peak value uniformity is greater than or equal to the second preset magnetic field value uniformity.
- the dividing unit 302 is configured to perform continuous meshing on the spatial extent of the coil system to be arranged, and each divided grid is regarded as a current ring, and the current value of each current ring includes a positive value, a negative value or zero; And according to the preset standard, the DSV surface of the superconducting magnet system and the ellipsoid surface of the 5 Gaussian stray field are evenly divided into thousands of target points.
- the first calculating unit 303 is configured to: the magnetic field peak-to-peak uniformity at each target point on the surface of the DSV is not greater than the first predetermined magnetic field peak-to-peak uniformity, and the magnetic field strength at each target point on the surface of the ellipsoid of the 5 Gaussian stray field
- the value constraint is the first constraint.
- the current distribution map is obtained by linear programming algorithm to calculate the current distribution at each grid which minimizes the volume of the coil system under a given operating current density.
- a second calculating unit 304 configured to respectively use the number of current values that are not zero in the current distribution map Head and space position, as the number and initial position of the solenoid coil in the coil system; limit of the magnetic field strength at each target point on the surface of the ellipsoid with 5 Gaussian stray field with the second preset magnetic field peak-to-peak uniformity limit condition
- the condition of the highest magnetic field strength in the line graph system limits the current safety margin limit of the coil system operation as the second constraint condition.
- the size of each solenoid coil is calculated by the nonlinear optimization algorithm when the volume of the coil system is the smallest.
- the positional parameter, the final structural parameters of the coil system are obtained.
- the size and position parameters of each solenoid coil include an inner radius, an outer radius, and an axial position of the two end portions of each solenoid coil.
- the spatial extent of the coil system to be arranged is specifically a solenoid-shaped region having a rectangular cross section, the range of which is rectangular
- the parameters of the rectangular section determine the inner diameter, outer diameter and length of the rectangular section.
- the two-dimensional continuous meshing may be performed on the rectangular cross-section in the radial direction and the axial direction, respectively, and the rectangular cross-section is along the radial direction.
- the axial directions are divided into several equal parts to form a two-dimensional continuous space grid.
- the current distribution diagram gives the magnitude and direction of the current at each grid, and the non-zero current values at the grid are gathered together.
- each positive position of the non-zero current cluster is used to arrange a forward current solenoid line ⁇
- each negative value position is used to arrange a reverse current solenoid line ⁇ , respectively.
- the second constraint includes:
- the peak-to-peak uniformity of the magnetic field of the axial magnetic field at each target point on the surface of the DSV is not greater than the uniformity of the second predetermined magnetic field value
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not greater than the preset magnetic field strength
- the current safety margin of the coil system operation is not greater than the preset current safety margin, and the preset current safety margin is calculated from the highest magnetic field strength value in the coil system and the critical characteristic of the superconducting wire selected by the coil system;
- the spacing between adjacent spiral lines in the turns system is greater than the preset spacing;
- the size and positional parameters of each solenoid line ⁇ place the coil system within the spatial extent of the system in which the coil is to be placed.
- the limiting value of the magnetic field strength value at each target point on the surface of the ellipsoid of the 5 Gaussian stray field may be: 5
- the magnetic field strength value at each target point on the surface of the ellipsoid of the Gaussian stray field is not more than 5 Gauss;
- the superimposed magnetic field strength of the axial magnetic field and the radial magnetic field at each target point on the surface of the ellipsoid of the Gaussian stray field is not greater than the preset magnetic field strength, which may be: 5
- the axial direction at each target point on the surface of the ellipsoid of the Gaussian stray field The superposed magnetic field strength of the magnetic field and the radial magnetic field is not more than 5 Gauss.
- the second calculating unit 304 can also be used to calculate the current safety margin of the coil system operation according to the highest magnetic field strength value in the preset coil system and the critical characteristic of the superconducting wire selected by the coil system. degree.
- the embodiment of the present invention further provides an MRI superconducting magnet system, and the structural parameters of the MRI superconducting magnet system can be obtained by the method for acquiring structural parameters of any of the MRI superconducting magnet systems of the present invention.
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
L'invention concerne un système d'aimant supraconducteur d'imagerie à résonance magnétique ainsi qu'un procédé et un dispositif pour acquérir des paramètres structuraux de celui-ci. Le procédé comprend : la réalisation d'une division en grille continue sur une plage spatiale où il faut disposer un système de bobine, la considération de chaque grille obtenue par la division en tant que cercle de courant, et la division uniforme d'une surface DSV d'un système d'aimant supraconducteur et d'une surface ellipsoïdale d'un champ de dispersion de 5 Gauss en une pluralité de points cibles ; au moyen d'un algorithme de programmation linéaire, le calcul de la distribution du courant à l'endroit de chaque grille, ce qui rend le volume du système de bobine minimum, de manière à obtenir un schéma de distribution du courant ; et le relevé du nombre et des positions spatiales des valeurs du courant différentes de zéro dans le schéma de distribution du courant respectivement en tant que nombre et que positions initiales des bobines de solénoïde dans le système de bobine et, au moyen d'un algorithme d'optimisation non linéaire, le calcul des paramètres de taille et de position de chaque bobine de solénoïde lorsque le volume du système de bobine est minimum de manière à acquérir les paramètres structuraux finaux du système de bobine. Le procédé et le dispositif peuvent améliorer l'efficacité d'acquisition des paramètres structuraux d'un système d'aimant supraconducteur d'IRM et rendre relativement petite une structure de bobine dans le système d'aimant supraconducteur d'IRM.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210147319.8 | 2012-05-11 | ||
CN201210147319.8A CN102707250B (zh) | 2012-05-11 | 2012-05-11 | 一种磁共振成像超导磁体系统的设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013166810A1 true WO2013166810A1 (fr) | 2013-11-14 |
Family
ID=46900180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2012/083287 WO2013166810A1 (fr) | 2012-05-11 | 2012-10-22 | Système d'aimant supraconducteur d'imagerie à résonance magnétique et procédé et dispositif pour acquérir des paramètres structuraux de celui-ci |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN102707250B (fr) |
WO (1) | WO2013166810A1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113971349A (zh) * | 2021-12-22 | 2022-01-25 | 华中科技大学 | 一种螺线管形超导磁体线圈数目和初始位置的获取方法 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707250B (zh) * | 2012-05-11 | 2014-04-02 | 中国科学院电工研究所 | 一种磁共振成像超导磁体系统的设计方法 |
CN104685584B (zh) * | 2012-11-23 | 2017-07-18 | 中国科学院高能物理研究所 | 用于构造磁共振成像超导磁体的方法 |
CN103499797B (zh) * | 2013-09-09 | 2016-03-02 | 中国科学院电工研究所 | 磁共振成像超导磁体螺线管线圈数目和初始位置获取方法 |
CN103592610B (zh) * | 2013-11-26 | 2015-12-30 | 中国科学院武汉物理与数学研究所 | 一种基于吸收线形二次矩的自动搜索匀场方法 |
CN104899441B (zh) * | 2015-06-05 | 2018-02-27 | 河海大学 | 一种适用于计算磁场vrms均匀度的数值方法 |
CN107871024B (zh) * | 2016-09-23 | 2023-07-18 | 中国电力科学研究院 | 一种高温超导环形储能磁体的电磁优化方法和装置 |
CN106556813B (zh) * | 2016-11-25 | 2021-09-24 | 上海辰光医疗科技股份有限公司 | 一种磁共振系统中主动匀场线圈的线性混合优化方法 |
CN108802645B (zh) * | 2017-10-24 | 2021-06-11 | 中国计量大学 | 基于0-1整数规划的永磁型纵向梯度线圈设计方法 |
CN108107390A (zh) * | 2017-12-29 | 2018-06-01 | 鑫高益医疗设备股份有限公司 | 一种超导磁体外磁屏蔽线圈的优化设计方法 |
CN108491623B (zh) * | 2018-03-21 | 2023-05-12 | 南京磁升超导技术有限公司 | 一种降低高场超导磁体线圈应力的方法 |
CN110632541B (zh) * | 2018-06-22 | 2021-12-28 | 上海尚磁机电科技有限公司 | 一种磁共振系统的匀场方法和装置 |
CN109765510B (zh) * | 2019-02-20 | 2021-04-27 | 中国科学院电工研究所 | 一种带有圆角的径向超导匀场线圈设计方法 |
CN112395735B (zh) * | 2020-09-23 | 2023-01-24 | 核工业西南物理研究院 | 一种rebco导体载流效率的仿真方法 |
CN113343551B (zh) * | 2021-06-30 | 2022-10-28 | 北京航空航天大学 | 一种基于磁屏蔽筒耦合的圆柱形和球形匀强磁场线圈设计方法 |
CN114896939B (zh) * | 2022-07-12 | 2022-09-23 | 山东奥新医疗科技有限公司 | 超导磁体线圈设计方法、系统、装置及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6255929B1 (en) * | 1998-03-23 | 2001-07-03 | The Board Of Trustees Of The Leland Stanford Junior University | Method of making optimized, air-core electromagnets |
WO2002045582A1 (fr) * | 2000-12-07 | 2002-06-13 | New York University | Procedes d'optimisation de systemes d'imagerie par resonance magnetique |
WO2006065027A1 (fr) * | 2004-12-14 | 2006-06-22 | Korea Basic Science Institute | Procede de mise au point d'un aimant supraconducteur a champ magnetique fort |
CN101533077A (zh) * | 2009-04-17 | 2009-09-16 | 中国科学院电工研究所 | 一种磁共振成像装置超导磁体的优化设计方法 |
CN102707250A (zh) * | 2012-05-11 | 2012-10-03 | 中国科学院电工研究所 | 一种磁共振成像超导磁体系统的设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176368B (zh) * | 2011-01-24 | 2012-10-03 | 中国科学院高能物理研究所 | 一种用于磁共振成像超导磁体的优化设计方法 |
-
2012
- 2012-05-11 CN CN201210147319.8A patent/CN102707250B/zh not_active Expired - Fee Related
- 2012-10-22 WO PCT/CN2012/083287 patent/WO2013166810A1/fr active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6255929B1 (en) * | 1998-03-23 | 2001-07-03 | The Board Of Trustees Of The Leland Stanford Junior University | Method of making optimized, air-core electromagnets |
WO2002045582A1 (fr) * | 2000-12-07 | 2002-06-13 | New York University | Procedes d'optimisation de systemes d'imagerie par resonance magnetique |
WO2006065027A1 (fr) * | 2004-12-14 | 2006-06-22 | Korea Basic Science Institute | Procede de mise au point d'un aimant supraconducteur a champ magnetique fort |
CN101533077A (zh) * | 2009-04-17 | 2009-09-16 | 中国科学院电工研究所 | 一种磁共振成像装置超导磁体的优化设计方法 |
CN102707250A (zh) * | 2012-05-11 | 2012-10-03 | 中国科学院电工研究所 | 一种磁共振成像超导磁体系统的设计方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113971349A (zh) * | 2021-12-22 | 2022-01-25 | 华中科技大学 | 一种螺线管形超导磁体线圈数目和初始位置的获取方法 |
CN113971349B (zh) * | 2021-12-22 | 2022-04-08 | 华中科技大学 | 一种螺线管形超导磁体线圈数目和初始位置的获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102707250A (zh) | 2012-10-03 |
CN102707250B (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2013166810A1 (fr) | Système d'aimant supraconducteur d'imagerie à résonance magnétique et procédé et dispositif pour acquérir des paramètres structuraux de celui-ci | |
CN104685584B (zh) | 用于构造磁共振成像超导磁体的方法 | |
CN102376410B (zh) | 一种开放式核磁共振超导磁体的设计方法 | |
US5382904A (en) | Structured coil electromagnets for magnetic resonance imaging and method for fabricating the same | |
CN102176368B (zh) | 一种用于磁共振成像超导磁体的优化设计方法 | |
CN102150222B (zh) | 用于限定用于成像体积的磁场的磁组件和方法 | |
US9389291B2 (en) | Gradient coil, magnetic resonance imaging device, and method for designing coil pattern | |
CN101533077A (zh) | 一种磁共振成像装置超导磁体的优化设计方法 | |
CN112424625B (zh) | 电阻电磁体系统 | |
CN106556813A (zh) | 一种磁共振系统中主动匀场线圈的线性混合优化方法 | |
CN105445683B (zh) | 一种柱面横向自屏蔽梯度线圈设计方法 | |
CN110568390A (zh) | 一种双平面匀场线圈的设计方法 | |
Ni et al. | Globally optimal algorithm for design of 0.7 T actively shielded whole-body open MRI superconducting magnet system | |
CN104833930B (zh) | 开放式磁共振系统梯度线圈磁场强度的计算方法 | |
CN112307657A (zh) | 一种用于新生儿成像的超导磁共振磁体及设计方法 | |
US11740301B2 (en) | Eigenmode transmit array coil for magnetic resonance imaging | |
Tashiro et al. | Simple-Box-9 coil system: A novel approach to design of a square coil system for producing uniform magnetic fields | |
Siebold | Design optimization of main, gradient and RF field coils for MR imaging | |
Yue et al. | Design and magnetic field simulation of 3.0 T MRI superconducting magnet | |
CN103151136B (zh) | 一种非对称自屏蔽开放式磁共振成像超导磁体 | |
Sattarov et al. | High-field open MRI for breast cancer screening | |
CN118446077B (zh) | 一种低功率的三轴亥姆霍兹线圈模型的优化方法 | |
Dou et al. | Design of highly uniformity planar grid coils for obtaining near-zero magnetic field in magnetoencephalogram | |
Wang | Gradient coil design and acoustic noise control in magnetic resonance imaging systems | |
CN119471500A (zh) | 一种单端开口式高性能磁屏蔽系统及其设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12876555 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 12876555 Country of ref document: EP Kind code of ref document: A1 |