CN109276811B - Setting device and setting method for zero magnetic field point position - Google Patents
Setting device and setting method for zero magnetic field point position Download PDFInfo
- Publication number
- CN109276811B CN109276811B CN201811170387.XA CN201811170387A CN109276811B CN 109276811 B CN109276811 B CN 109276811B CN 201811170387 A CN201811170387 A CN 201811170387A CN 109276811 B CN109276811 B CN 109276811B
- Authority
- CN
- China
- Prior art keywords
- zero
- tumor
- magnetic field
- coil
- setting device
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N2/00—Magnetotherapy
- A61N2/02—Magnetotherapy using magnetic fields produced by coils, including single turn loops or electromagnets
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N2/00—Magnetotherapy
- A61N2/004—Magnetotherapy specially adapted for a specific therapy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Abstract
The invention provides a setting device and a setting method for a zero magnetic field point position, wherein the setting device for the zero magnetic field point position comprises the following steps: the device comprises 6 groups of coil assemblies, wherein 6 groups of coil assemblies enclose a center circle, each group of coil assemblies comprises 3 coils, the circle center of each coil is uniformly distributed on the center circle according to every 120 degrees, a connecting line between the circle center of the center circle and the circle center of each coil divides the center circle into 18 areas in a clock clockwise direction, corresponding current values are loaded in each coil, and the current values in each coil are determined according to position information of a tumor in the setting device so as to control the magnetic field intensity at the position of the tumor to be zero or close to zero. Through the technical scheme of the invention, the whole device has reasonable structure, can accurately control the magnetic field at the position of the tumor in the device to be zero or close to zero, ensures other organs and cells in a non-focus area, and has higher practical and research values.
Description
Technical Field
The invention relates to the technical field of medical equipment, in particular to a device and a method for setting a zero magnetic field point position.
Background
Tumor thermotherapy uses the physiological phenomenon that cancer cells are thermolabile than normal cells to kill tumor cells after heating tissues, is a thermotherapy method called high-temperature thermotherapy, and is a green therapy following surgery, radiotherapy and chemotherapy. At present, thermal therapy can be divided into methods of microwave radiation, radio frequency radiation, ultrasonic focusing, resistance heating, alternating magnetic field heating magnetic particles, and the like according to heating methods.
Heating magnetic particles by an alternating magnetic field is a new thermotherapy mode in recent years, and a magnetic particle material is injected into pathological tissues, and the alternating magnetic field is applied to the outside, so that the magnetic material generates heat due to hysteresis loss. During the thermal treatment, the magnetic particles inevitably diffuse to the normal tissue due to their fluidity, and thus the normal tissue is also heated and even damaged. In order to protect normal tissues, a static magnetic field is established around the tumor to protect normal organs and cells around the lesion area. The static magnetic field of the tumor region is required to be as zero as possible, so that the injected magnetic nanoparticles are ensured to directly reach the tumor region by utilizing the medicinal property of the carrier, and other organs and cells in a non-focus region are protected by the static magnetic field, so that how to make the magnetic field close to zero at the tumor position in the constructed thermotherapy device becomes a technical problem to be solved urgently.
Disclosure of Invention
Based on at least one of the technical problems, the invention provides a new setting scheme of the zero magnetic field point position, the whole device has reasonable structure, the magnetic field at the position of the tumor in the device can be accurately controlled to be zero or close to zero, other organs and cells in a non-focus area are effectively ensured, and the device has higher practical and research values.
In view of the above, the present invention provides a setting device for a zero magnetic field point position, comprising: the device comprises 6 groups of coil assemblies, wherein 6 groups of coil assemblies enclose a center circle, each group of coil assemblies comprises 3 coils, the circle center of each coil is uniformly distributed on the center circle according to every 120 degrees, a connecting line between the circle center of the center circle and the circle center of each coil divides the center circle into 18 areas in a clock clockwise direction, corresponding current values are loaded in each coil, and the current values in each coil are determined according to position information of a tumor in the setting device so as to control the magnetic field intensity at the position of the tumor to be zero or close to zero.
In the technical scheme, the circle centers of all coils in the setting device are positioned on the same plane, the magnetic field strength in the vertical direction can be mutually offset on the plane, the calculated amount of the subsequent process is greatly reduced when the position of the tumor is positioned on the plane, in addition, the current value determined according to the position information of the tumor in the setting device is loaded in each coil, the magnetic field strength of the position of the whole setting device where the tumor is positioned can be controlled to be zero or close to zero, other organs and cells in a non-focus area are effectively ensured, the whole device is reasonable in structure, and the device has higher practical and research values.
In the above technical solution, preferably, the diameter of each coil is 80mm, and the diameter of the central circle is 800 mm.
In the technical scheme, the universality of the whole setting device can be improved by setting the diameter of the coil to be 80mm and the diameter of the central circle to be 800mm, and the actual use requirement is met.
According to a second aspect of the present invention, a method for setting a zero magnetic field point position is provided, which is based on the method for setting a zero magnetic field point position in the above technical solution, and includes: acquiring position information of a tumor in the setting device; bringing the acquired position information into a preset solving model, and carrying out solving analysis based on an improved bat algorithm to obtain 18 current values; and respectively loading the 18 current values obtained by the solution into corresponding coils in the setting device so as to control the magnetic field intensity at the position of the tumor to be zero or close to zero.
In the technical scheme, 18 current values are obtained by modeling the position information of the tumor and performing a bat algorithm, so that the error in the calculation process can be effectively reduced, the calculation efficiency can be improved, the accuracy of the calculation result can be ensured, and the magnetic field intensity of the position of the tumor can be ensured to be zero or close to zero after the obtained accurate current values are respectively loaded into the corresponding current loops.
In the above technical solution, preferably, the preset solution model is a mathematical model constructed based on a geometric shape of the setting device.
In any one of the above technical solutions, preferably, the improved bat algorithm is:
wherein the content of the first and second substances,
is the spatial position of the bat individual i in the t +1 th search; x is the number of * Is the current global optimal position; levy (λ) is a random perturbation vector whose step size μ obeys Levy distribution, λ -3/2 is a scale parameter;
representing a vector operation; t is max Is the maximum number of iterations and t is the current number of iterations.
In the technical scheme, the Levy disturbance amount is gradually reduced while the iteration times are increased, and multiple small disturbances and occasional large disturbances can increase the diversity of the population when the position of the bat is updated by using a Levy disturbance strategy, thereby expanding the search space of the bat individual, improving the optimizing capability of the bat algorithm in a high-dimensional space and avoiding the bat individual from falling into local optimization.
Through the technical scheme, the whole device is reasonable in structure, the magnetic field at the position where the tumor is located in the device can be accurately controlled to be zero or close to zero, other organs and cells in a non-focus area are effectively ensured, and the device has high practical and research values.
Drawings
FIG. 1 is a schematic three-dimensional structure diagram of a setting device for a zero-field point position according to an embodiment of the invention;
FIG. 2 shows a top view of a setup device for zero field point locations according to an embodiment of the invention;
FIG. 3 shows a flow diagram of a method of setting a zero field point position according to an embodiment of the invention;
FIG. 4 shows a flow diagram of an improved bat algorithm according to an embodiment of the present invention;
FIG. 5A shows a graph of simulated effect of a tumor location at (0,0,0) according to an embodiment of the invention;
FIG. 5B shows a graph of simulated effect of tumor location at (0,100,0) according to an embodiment of the invention;
fig. 5C shows a graph of simulated effect of tumor location at (200, 100,0) according to an embodiment of the invention.
Detailed Description
In order that the above objects, features and advantages of the present invention may be more clearly understood, there will now be described in connection with the following examples, without conflict, the present disclosure and features of the examples may be combined with each other.
The technical solution of the present invention is further explained with reference to fig. 1 to 5C as follows:
as shown in fig. 1, the setting device for the zero magnetic field point position includes: 6 groups of coil assemblies and 6 groups of coil assemblies enclose a central circle, wherein each group of coil assemblies comprises 3 coils, as shown in fig. 2, the circle center of each coil is uniformly distributed on the central circle every 120 degrees, a connecting line between the circle center of the central circle and the circle center of each coil divides the central circle into 18 areas in a clock clockwise direction, under the condition of no grouping, the included angle between every two adjacent coils is 20 degrees, in order to increase the universality of the whole setting device and meet the actual use requirement, the diameter of each coil can be set to be 80mm, the diameter of the central circle is set to be 800mm, in addition, a corresponding current value is loaded in each coil, and the current value in each coil is determined according to the position information of the tumor in the setting device so as to control the magnetic field intensity at the position of the tumor to be zero or close to zero.
In particular, a method of controlling the magnetic field strength at a location of a tumor to be zero or close to zero, as shown in particular in fig. 3, comprises the steps of:
Specifically, the thermal therapy object is placed in a circle surrounded by 18 coils (i.e. a central circle) in the device, and position information of the tumor is obtained under a coordinate system established as shown in fig. 2.
And step 304, bringing the acquired position information into a preset solving model, and carrying out solving analysis based on an improved bat algorithm to obtain 18 current values.
In the established predetermined solution model (which is a mathematical model constructed based on the geometry of the setup device), the coil is approximated to a circular current. According to the Biot-Saval law, the magnetic field component expressions excited by the circular current at the space point under a rectangular coordinate system established by taking the circle center of the circular current as the origin and the central axis as the z axis are respectively as follows:
wherein, mu 0 For vacuum permeability, R is the coil radius and I is the coil current. By means of mathematical software, the distribution of the magnetic induction intensity of the circular current magnetic field in the whole three-dimensional space can be obtained by utilizing the Simpson integration method for calculation. Magnetic field of Z-axis when tumor is located in XY planeThe induction strengths cancel each other to zero. From the above formula, the magnetic field strength generated by each coil at the tumor location can be calculated, and by vector addition, the modulus of the magnetic field strength at the tumor location can be obtained. Wherein the vector addition formula is as follows:
wherein, BX i For the magnetic induction produced BY coil i in the direction of the X-axis at the tumor location, BY i For the magnetic induction produced by coil i in the direction of the Y-axis at the tumor site, BZ i The magnetic induction intensity generated by the coil i in the Z-axis direction at the position of the tumor is obtained, the magnetic induction intensity at the position of the tumor is obtained, and the mode of the coil B is required to be zero or close to zero in order to meet the requirement that the magnetic field intensity at the position of the tumor is zero and close to zero.
Fig. 4 specifically shows a process of performing solution analysis based on the improved bat algorithm (wLevyOLBA), which includes:
before the algorithm starts iteration, each bat is reversely learned by the following formula to obtain a more excellent initial population.
x’ i ＝x max +x min -x i
In the formula, x max And x min Position information, x, when the fitness function takes the maximum and minimum values, respectively i As position information of the current particle, x i ' is its reverse position information.
The position updating in the iteration process is carried out by the following formula, so that the bat realizes Levy disturbance when the position is updated.
Wherein the content of the first and second substances,
is the spatial position of the bat individual i in the t +1 th search; x is the number of * Is the current global optimal position; levy (λ) is a random perturbation vector whose step size μ obeys Levy distribution, λ -3/2 is a scale parameter;
representing a vector operation; t is max The maximum iteration number is t, the current iteration number is t, the Levy disturbance amount is gradually reduced while the iteration number is increased, multiple small disturbances and occasional large disturbances are generated when the position of the bat is updated by using a Levy disturbance strategy, the diversity of the population is increased, the search space of the bat individual is enlarged, the optimizing capability of the bat algorithm in a high-dimensional space is improved, and the bat individual is prevented from falling into local optimization.
Wherein, the Levy process is a markov process, the function value satisfies a Heavy-tailed distribution, and the formula and the calculation formula of the occurrence probability P (μ) of the function value μ are respectively as follows:
u＝rand(n,1)×σ
v＝rand(n,1)
where n is the dimension of the solution and the gamma function is
α is the step factor.
And step 306, loading the 18 solved current values into corresponding coils in the setting device respectively to control the magnetic field intensity at the position of the tumor to be zero or close to zero.
Specifically, the currents of the 18 peripheral coils are calculated according to an algorithm, and the current values are input into the coils to determine the range of the tumor region so that the static magnetic field of the region where the tumor is located is zero as much as possible.
In the above embodiment, the centers of circles of all coils in the setting device are all located on the same plane, and on the plane, the magnetic field strengths in the vertical direction can be mutually offset, and in the subsequent process, when the tumor position is located on the plane, the calculated amount is greatly reduced, and in addition, 18 current values are solved by modeling the position information of the tumor and bat algorithm, so that the error in the calculation process can be effectively reduced, the calculation efficiency is improved, and the accuracy of the calculation result is ensured.
Simulation experiments were also performed with comsol simulation software. Specifically, assume that the tumor center position is located in the XY plane, which is (0,0,0), (0,100,0), (200, 0), respectively. The number of turns of the coil is 1000, the wLevyOLBA algorithm carries out solving calculation on a preset solving model according to the tumor position, 18 solved current values are substituted into a simulation model of comsol simulation software, and the solved current values are shown in tables 1-3, wherein the table 1 is the current of 18 coils at the tumor position (0,0, 0); table 2 is the current for the 18 coils at the tumor location (0,100, 0); table 3 currents of 18 coils at tumor location (200, 0):
TABLE 1
TABLE 2
TABLE 3
The 18 current values solved in tables 1 to 3 are respectively substituted into the simulation model of the comsol simulation software to obtain the anti-true effect graphs shown in fig. 5A, fig. 5B and fig. 5C, where fig. 5A is the simulation effect graph with the tumor center at (0,0,0), fig. 5B is the simulation effect graph with the tumor center at (0,100,0), and fig. 5C is the simulation effect graph with the tumor center at (200, 0), and from the simulation effect graph, the tumor position area has a smaller magnetic field intensity compared with the surrounding area, a zero magnetic field or a near zero magnetic field area is formed at the tumor position center, the tumor is heated at the tumor position, and the tissues around the tumor are protected, thereby verifying the effectiveness of the setting device and the setting method in the above embodiments.
The technical scheme of the invention is explained in detail by combining the drawings, the technical scheme of the invention provides a new setting scheme of the zero magnetic field point position, the whole device is reasonable in structure, the magnetic field at the position of the tumor in the device can be accurately controlled to be zero or close to zero, other organs and cells in a non-focus area are effectively ensured, and the device has higher practical and research values.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (3)
1. A setting device for a zero magnetic field point position, comprising: the device comprises: 6 groups of coil assemblies, wherein 6 groups of coil assemblies enclose a center circle, each group of coil assemblies comprises 3 coils, the circle center of each coil is uniformly distributed on the center circle according to every 120 degrees, a connecting line between the circle center of the center circle and the circle center of each coil divides the center circle into 18 areas in a clock clockwise direction, wherein each coil is loaded with a corresponding current value, and the current value of each coil is determined according to the position information of a tumor in the setting device so as to control the magnetic field intensity at the position of the tumor to be zero or close to zero;
the setting method comprises the following steps: acquiring position information of a tumor in the setting device;
bringing the acquired position information into a preset solving model, and carrying out solving analysis based on an improved bat algorithm to obtain 18 current values;
in the established predetermined solving model, the predetermined solving model is a mathematical model established based on the geometric shape of the setting device, the coil is approximated to be a circular current, and the expressions of magnetic field components excited by the circular current at a space point are respectively obtained according to the Biot-Saval law under a rectangular coordinate system established by taking the circle center of the circular current as the origin and the central axis as the z axis:
wherein, mu 0 The magnetic induction intensity distribution of a circular current magnetic field in the whole three-dimensional space can be obtained by calculating by utilizing a Simpson integration method through mathematical software, when a tumor is positioned on an XY plane, the magnetic induction intensities of a Z axis are mutually offset to be zero, the magnetic field intensity generated by each coil at the position of the tumor can be calculated according to the above formula, and the mode of the magnetic field intensity at the position of the tumor can be obtained through vector addition, wherein the vector addition formula is as follows:
wherein, BX i Magnetic field generated for coil i in X-axis direction at tumor positionInduction intensity, BY i For the magnetic induction produced by coil i in the direction of the Y-axis at the tumor site, BZ i The magnetic induction intensity generated by the coil i in the Z-axis direction at the position of the tumor is adopted, B is the magnetic induction intensity at the position of the tumor, and the mode of B is required to be zero or close to zero in order to meet the requirement that the magnetic field intensity at the position of the tumor is zero and close to zero;
respectively loading the 18 current values obtained by the solution into corresponding coils in the setting device to control the magnetic field intensity at the position of the tumor to be zero or close to zero;
the improved bat algorithm is as follows:
wherein x is i t+1 Is the spatial position of the bat individual i in the t +1 th search; x is the number of * Is the current global optimal position; levy (λ) is a random perturbation vector whose step size μ obeys Levy distribution, λ -3/2 is a scale parameter;
representing a vector operation; t is max Is the maximum number of iterations and t is the current number of iterations.
2. The setting device of the zero magnetic field point position according to claim 1, wherein each of the coils has a diameter of 80mm, and the center circle has a diameter of 800 mm.
3. The setting device of the zero magnetic field point position according to claim 1, wherein the predetermined solution model is a mathematical model constructed based on the geometry of the setting device.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811170387.XA CN109276811B (en) | 2018-10-09 | 2018-10-09 | Setting device and setting method for zero magnetic field point position |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811170387.XA CN109276811B (en) | 2018-10-09 | 2018-10-09 | Setting device and setting method for zero magnetic field point position |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109276811A CN109276811A (en) | 2019-01-29 |
| CN109276811B true CN109276811B (en) | 2022-08-05 |
Family
ID=65177254
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811170387.XA Active CN109276811B (en) | 2018-10-09 | 2018-10-09 | Setting device and setting method for zero magnetic field point position |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109276811B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1135360A (en) * | 1995-04-25 | 1996-11-13 | 王兆伦 | Electronic flux knife mode and method thereof |
| CN1545985A (en) * | 2003-12-03 | 2004-11-17 | 宏 朱 | Magnetic heat treatment therapeutic apparatus for tumour |
| CN105169560A (en) * | 2015-06-12 | 2015-12-23 | 郑州轻工业学院 | Device and method used for controlling magnetic nano particle heating area |
| CN107968613A (en) * | 2017-11-20 | 2018-04-27 | 无锡开放大学 | A kind of permanent magnet synchronous motor rotational speed governor based on Recurrent Fuzzy Neural Network |
| CN108416421A (en) * | 2018-03-09 | 2018-08-17 | 大连大学 | The dynamic Algorithm of Firepower Allocation of bat algorithm is improved based on DDE |
-
2018
- 2018-10-09 CN CN201811170387.XA patent/CN109276811B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1135360A (en) * | 1995-04-25 | 1996-11-13 | 王兆伦 | Electronic flux knife mode and method thereof |
| CN1545985A (en) * | 2003-12-03 | 2004-11-17 | 宏 朱 | Magnetic heat treatment therapeutic apparatus for tumour |
| CN105169560A (en) * | 2015-06-12 | 2015-12-23 | 郑州轻工业学院 | Device and method used for controlling magnetic nano particle heating area |
| CN107968613A (en) * | 2017-11-20 | 2018-04-27 | 无锡开放大学 | A kind of permanent magnet synchronous motor rotational speed governor based on Recurrent Fuzzy Neural Network |
| CN108416421A (en) * | 2018-03-09 | 2018-08-17 | 大连大学 | The dynamic Algorithm of Firepower Allocation of bat algorithm is improved based on DDE |
Non-Patent Citations (2)
| Title |
|---|
| "Parameters identification of Krasnosel"skii-Pokrovskii model for piezo-actuated stages using a modified bat optimization algorithm based on levy fights trajectory";Rui Xu 等;《2017 Chinese Automation Congress (CAC)》;20180101;第5304-5309页 * |
| "磁性纳米粒子在肿瘤热疗中的应用";张向梅 马吉明 苏日建 等;《生物技术世界》;20140630(第6期);第85-86页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109276811A (en) | 2019-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Van Heijst et al. | MR-guided breast radiotherapy: feasibility and magnetic-field impact on skin dose | |
| EP2514482B1 (en) | Treatment planning apparatus and particle therapy system | |
| US20210339049A1 (en) | Radiotherapy treatment plans using differentiable dose functions | |
| US7809107B2 (en) | Method for controlling modulation strength in radiation therapy | |
| CN110691627A (en) | Distance compensator for radiation therapy | |
| Bedford et al. | Constrained segment shapes in direct‐aperture optimization for step‐and‐shoot IMRT | |
| JP2006142043A (en) | Automated optimization method for radiotherapy dose | |
| US7196603B2 (en) | Magnetic coil design using optimization of sinusoidal coefficients | |
| CN105169560B (en) | A kind of device and method controlling magnetic nanometer heating region | |
| US20200346033A1 (en) | Radiation treatment planning with multiple target subset optimization | |
| Locke et al. | Trajectory optimization in radiotherapy using sectioning (TORUS) | |
| Sheng et al. | Atlas-guided prostate intensity modulated radiation therapy (IMRT) planning | |
| US20210244970A1 (en) | Systems and methods for planning, controlling and/or delivering radiotherapy and radiosurgery using combined optimization of dynamic axes (coda) | |
| Modiri et al. | Inverse 4D conformal planning for lung SBRT using particle swarm optimization | |
| US9694203B2 (en) | Method for preparing an irradiation plan for a moving target volume without a movement compensation | |
| Van Kesteren et al. | The dosimetric impact of leaf interdigitation and leaf width on VMAT treatment planning in Pinnacle: comparing Pareto fronts | |
| Li et al. | Extension and validation of a GPU‐Monte Carlo dose engine gDPM for 1.5 T MR‐LINAC online independent dose verification | |
| CN109276811B (en) | Setting device and setting method for zero magnetic field point position | |
| Zhang et al. | A surrogate-based metaheuristic global search method for beam angle selection in radiation treatment planning | |
| Smith et al. | Optimization of spot and trimmer position during dynamically collimated proton therapy | |
| Cevik et al. | Modeling and comparison of alternative approaches for sector duration optimization in a dedicated radiosurgery system | |
| CN106902480B (en) | A kind of parallel Quantum annealing target spot distribution calculation method | |
| de Battisti et al. | An automated optimization tool for high-dose-rate (HDR) prostate brachytherapy with divergent needle pattern | |
| JP2015150184A (en) | Treatment planning device and storage medium storing arithmetic method for treatment plan | |
| JP6823299B2 (en) | Radiation therapy planning device and radiation therapy system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |