JP2024502075A - 放射線照射システム及びその制御方法 - Google Patents
放射線照射システム及びその制御方法 Download PDFInfo
- Publication number
- JP2024502075A JP2024502075A JP2023540520A JP2023540520A JP2024502075A JP 2024502075 A JP2024502075 A JP 2024502075A JP 2023540520 A JP2023540520 A JP 2023540520A JP 2023540520 A JP2023540520 A JP 2023540520A JP 2024502075 A JP2024502075 A JP 2024502075A
- Authority
- JP
- Japan
- Prior art keywords
- irradiation system
- radiation irradiation
- particles
- neutron
- particle
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 213
- 230000005855 radiation Effects 0.000 title claims abstract description 186
- 238000004364 calculation method Methods 0.000 claims abstract description 175
- 238000004088 simulation Methods 0.000 claims abstract description 139
- 230000008569 process Effects 0.000 claims abstract description 115
- 238000006243 chemical reaction Methods 0.000 claims abstract description 20
- 238000012216 screening Methods 0.000 claims abstract description 5
- 239000002245 particle Substances 0.000 claims description 286
- 238000011282 treatment Methods 0.000 claims description 216
- 238000007493 shaping process Methods 0.000 claims description 52
- 238000010521 absorption reaction Methods 0.000 claims description 51
- 229910052796 boron Inorganic materials 0.000 claims description 39
- 230000001225 therapeutic effect Effects 0.000 claims description 39
- 238000002560 therapeutic procedure Methods 0.000 claims description 38
- 239000000463 material Substances 0.000 claims description 31
- 230000001133 acceleration Effects 0.000 claims description 19
- 230000009467 reduction Effects 0.000 claims description 18
- 206010028980 Neoplasm Diseases 0.000 claims description 15
- 239000006185 dispersion Substances 0.000 claims description 14
- 229910052744 lithium Inorganic materials 0.000 claims description 13
- 229910052757 nitrogen Inorganic materials 0.000 claims description 13
- 229910052791 calcium Inorganic materials 0.000 claims description 12
- 229910052799 carbon Inorganic materials 0.000 claims description 12
- 229910052801 chlorine Inorganic materials 0.000 claims description 12
- 229910052739 hydrogen Inorganic materials 0.000 claims description 12
- 230000033001 locomotion Effects 0.000 claims description 12
- 229910052749 magnesium Inorganic materials 0.000 claims description 12
- 229910052760 oxygen Inorganic materials 0.000 claims description 12
- 229910052698 phosphorus Inorganic materials 0.000 claims description 12
- 229910052708 sodium Inorganic materials 0.000 claims description 12
- 229910052790 beryllium Inorganic materials 0.000 claims description 11
- 229910052717 sulfur Inorganic materials 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 229910052782 aluminium Inorganic materials 0.000 claims description 9
- 229910052731 fluorine Inorganic materials 0.000 claims description 9
- 229910052734 helium Inorganic materials 0.000 claims description 9
- 229910052710 silicon Inorganic materials 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 230000008033 biological extinction Effects 0.000 claims description 6
- 230000007423 decrease Effects 0.000 claims description 6
- 239000008280 blood Substances 0.000 claims description 5
- 210000004369 blood Anatomy 0.000 claims description 5
- 210000003205 muscle Anatomy 0.000 claims description 5
- 210000000988 bone and bone Anatomy 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 description 31
- 210000004027 cell Anatomy 0.000 description 20
- 210000001519 tissue Anatomy 0.000 description 19
- 210000004881 tumor cell Anatomy 0.000 description 14
- 238000009826 distribution Methods 0.000 description 12
- 239000003814 drug Substances 0.000 description 11
- 229940079593 drug Drugs 0.000 description 11
- 238000001959 radiotherapy Methods 0.000 description 10
- 239000011575 calcium Substances 0.000 description 8
- 239000000460 chlorine Substances 0.000 description 8
- 239000011777 magnesium Substances 0.000 description 8
- 239000011734 sodium Substances 0.000 description 8
- 239000006096 absorbing agent Substances 0.000 description 7
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- 238000000342 Monte Carlo simulation Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 201000011510 cancer Diseases 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 229910052700 potassium Inorganic materials 0.000 description 5
- 239000004698 Polyethylene Substances 0.000 description 4
- 230000001413 cellular effect Effects 0.000 description 4
- 230000008021 deposition Effects 0.000 description 4
- 230000004992 fission Effects 0.000 description 4
- 229920000573 polyethylene Polymers 0.000 description 4
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- ZAMOUSCENKQFHK-UHFFFAOYSA-N Chlorine atom Chemical compound [Cl] ZAMOUSCENKQFHK-UHFFFAOYSA-N 0.000 description 3
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 3
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 3
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 3
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 150000002431 hydrogen Chemical class 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 239000011574 phosphorus Substances 0.000 description 3
- 239000011591 potassium Substances 0.000 description 3
- 229910018072 Al 2 O 3 Inorganic materials 0.000 description 2
- KLZUFWVZNOTSEM-UHFFFAOYSA-K Aluminium flouride Chemical compound F[Al](F)F KLZUFWVZNOTSEM-UHFFFAOYSA-K 0.000 description 2
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 2
- 229910052786 argon Inorganic materials 0.000 description 2
- ATBAMAFKBVZNFJ-UHFFFAOYSA-N beryllium atom Chemical compound [Be] ATBAMAFKBVZNFJ-UHFFFAOYSA-N 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 229910052745 lead Inorganic materials 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000002139 neutron reflectometry Methods 0.000 description 2
- 229910052759 nickel Inorganic materials 0.000 description 2
- -1 polyethylene Polymers 0.000 description 2
- 230000002105 relative biological effectiveness Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000011593 sulfur Substances 0.000 description 2
- 235000013619 trace mineral Nutrition 0.000 description 2
- 239000011573 trace mineral Substances 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 229910016569 AlF 3 Inorganic materials 0.000 description 1
- 229910004261 CaF 2 Inorganic materials 0.000 description 1
- GUTLYIVDDKVIGB-OUBTZVSYSA-N Cobalt-60 Chemical compound [60Co] GUTLYIVDDKVIGB-OUBTZVSYSA-N 0.000 description 1
- YZCKVEUIGOORGS-OUBTZVSYSA-N Deuterium Chemical compound [2H] YZCKVEUIGOORGS-OUBTZVSYSA-N 0.000 description 1
- 201000010915 Glioblastoma multiforme Diseases 0.000 description 1
- 238000012879 PET imaging Methods 0.000 description 1
- WUKWITHWXAAZEY-UHFFFAOYSA-L calcium difluoride Chemical compound [F-].[F-].[Ca+2] WUKWITHWXAAZEY-UHFFFAOYSA-L 0.000 description 1
- 229910001634 calcium fluoride Inorganic materials 0.000 description 1
- 230000005880 cancer cell killing Effects 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 229910052805 deuterium Inorganic materials 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000009201 electron therapy Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 208000005017 glioblastoma Diseases 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000011344 liquid material Substances 0.000 description 1
- 229910001338 liquidmetal Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 201000001441 melanoma Diseases 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 238000002727 particle therapy Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000012636 positron electron tomography Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 229910052715 tantalum Inorganic materials 0.000 description 1
- GUVRBAGPIYLISA-UHFFFAOYSA-N tantalum atom Chemical compound [Ta] GUVRBAGPIYLISA-UHFFFAOYSA-N 0.000 description 1
- 238000002626 targeted therapy Methods 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
- A61N2005/1035—Simulated annealing
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N2005/1085—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
- A61N2005/109—Neutrons
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiation-Therapy Devices (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Radiation (AREA)
- Heating, Cooling, Or Curing Plastics Or The Like In General (AREA)
Abstract
本発明の実施例に係る放射線照射システム及びその制御方法は、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションすることにより、計算速度を大幅に向上させ、計算時間を減少させることができる。
Description
本発明の一態様は、放射線照射システムに関し、本発明の別の態様は、放射線照射システムの制御方法に関する。
原子科学の発展に従って、コバルト60、線形加速器、電子ビームなどの放射線療法は、既にがん治療の主な手段の1つとなった。しかしながら、従来の光子又は電子療法は、放射線そのものの物理的条件の制限で腫瘍細胞を殺すとともに、ビーム経路上の数多くの正常組織に損傷を与える。また、腫瘍細胞の放射線に対する感受性の度合いが異なるため、従来の放射線療法では、放射線耐性の高い悪性腫瘍(例えば、多形神経膠芽腫(glioblastoma multiforme)、黒色腫(melanoma))に対する治療効果が高くない。
腫瘍の周囲の正常組織への放射線損傷を低減するために、化学療法(chemotherapy)における標的療法の構想が、放射線療法に適用されている。また、放射線耐性の高い腫瘍細胞に対し、現在では、陽子線治療、重粒子治療、中性子捕捉療法などの、生物学的効果比(relative biological effectiveness、RBE)の高い放射線源が積極的に開発されている。このうち、中性子捕捉療法は、上記2つの構想を組み合わせたものであり、例えば、ホウ素中性子捕捉療法(Boron Neutron Capture Therapy、BNCT)では、ホウ素含有薬物が腫瘍細胞に特異的に集まり、高精度なビームの制御と合わせることで、従来の放射線に比べて、より良いがん治療オプションを提供する。
ホウ素中性子捕捉療法は、ホウ素(10B)含有薬物が熱中性子に対して大きい捕捉断面積を持つという特性を利用し、10B(n,α)7Li中性子捕捉と核分裂反応により、4Heと7Liという2種の重荷電粒子を生成する。2種の粒子の総飛程が約1つの細胞のサイズに相当するため、生体への放射線損傷が細胞レベルに抑えられる。ホウ素含有薬物が腫瘍細胞に選択的に集まり、適切な中性子ビーム源と合わせることで、正常組織に大きな損傷を与えない前提で、腫瘍細胞を部分的に殺すという目的を達成することができる。
放射粒子ががん細胞をできる限り多く殺すとともに正常細胞に対する損傷を低減するために、通常、患者を治療する前にCT又はPETの画像診断を行い、画像診断結果に基づいて人体の組織材料情報を取得し、材料情報及び放射線源に基づいて計算モデルを確立し、放射粒子の人体内の輸送過程をシミュレーションし、最終的に放射粒子の人体内の線量分布を取得し、そして患者に対する線量分布の最適な計画を患者の治療スキームとして選択する。
現在、放射線治療計画システムにおける線量計算モジュールは、主に、モンテカルロ法で放射粒子をシミュレーションして取得されたものである。従来の放射線治療は、光子及び電子の運動過程をシミュレーションする必要があるのに対して、放射線照射治療は、中性子及び光子の運動過程をシミュレーションする必要がある。現在、モンテカルロ法は、線量計算の最正確な方法であるが、計算は、時間がかかり、メモリ消費量が大きい。
現在、放射線治療の計算には、MCNP、Geant4のような汎用性モンテカルロプログラムが使用されることが多く、MCNPは、最初に原子炉設計計算に使用され、Geant4は、最初に高エネルギー物理計算に使用され、設計初期に放射線治療の物理シーンが考慮されないため、放射線治療の計算に対して特殊な最適化を行わない。放射線治療計画の作成は、通常、所定の時間内に完了する必要があり、例えば、放射線照射治療は、放射線治療計画システムが1時間内に治療スキームを作成できることが要求され、線量計算過程は、放射線治療計画の作成時間の大部分を占めるため、放射線治療計画の作成時間を減少させるために線量計算方法を最適化する必要がある。
従来技術の欠陥を解消するために、本発明の第1態様に係る放射線照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成し、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションする治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含む。
他の実施例において、治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaから選択される1つ又は複数の元素である。
より具体的には、治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaである。
別の実施例において、治療計画モジュールは、各元素の対応する中性子断面積データベースから、294K及び0K以外の温度に対応するデータベースを除去する。
放射線照射システムは、好ましくは、中性子捕捉療法システムであり、ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、中性子発生装置は、加速器及びターゲットを含み、加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつターゲットと作用して中性子線を発生させ、ビーム整形体は、中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、中性子発生装置により発生した中性子線は、ビーム整形体を通過して治療台上の被照射体に照射される。
本発明の第2態様に係る放射線照射システムの制御方法において、ビーム照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含み、放射線照射システムの制御方法は、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションするステップを含む。
他の実施例において、治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaから選択される1つ又は複数の元素である。
より具体的には、治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaである。
別の実施例において、制御方法は、治療計画モジュールが、各元素の対応する中性子断面積データベースから、294K及び0K以外の温度に対応するデータベースを除去するステップをさらに含む。
放射線照射システムは、好ましくは、中性子捕捉療法システムであり、ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、中性子発生装置は、加速器及びターゲットを含み、加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつターゲットと作用して中性子線を発生させ、ビーム整形体は、中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、中性子発生装置により発生した中性子線は、ビーム整形体を通過して治療台上の被照射体に照射される。
本発明の実施例に係る放射線照射システム及びその制御方法は、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションすることにより、治療計画モジュールの計算速度を大幅に向上させ、計算時間を減少させることができる。
本発明の第3態様に係る放射線照射システムは、治療用ビームを発生させ、かつ前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含む。
放射線照射システムは、好ましくは、中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される。
一実施例において、CPUのプロセス並列、スレッド並列及びGPU加速により粒子をシミュレーションし、プロセス並列及びスレッド並列は、いずれもマルチコアCPUを使用して並列計算を実現し、CPUのプロセス並列及びスレッド並列の計算過程は、まず、システムがプロセス数又はスレッド数を取得して、値nを得るステップと、次に、システムがシミュレーションすべき粒子をn個に均一に分けるステップと、そして、各スレッド又は各プロセスにおいてそれぞれ1個の粒子を独立してシミュレーションし、かつカウントするステップと、最後に、システムが各プロセス又は各スレッドにおいて取得されたカウント数を統計して、最終的な線量を取得するステップと、を含む。
一実施例において、GPU加速は、GPUのマルチコアプロセッサの並列計算により実現され、GPU加速によるシミュレーション計算の過程は、まず、システムが乱数、断面積データなどをCPUメモリからGPUグラフィックメモリに伝送し、次に、GPU中の各プロセッサが単一の粒子をシミュレーションし、計算し、カウントし、かつカウント結果をグローバルカウントに統計するステップと、そして、システムが、シミュレーションされていない粒子があるか否かを判断し、シミュレーションされていない粒子がない場合、カウント結果をGPUメモリからCPUメモリに伝送し、シミュレーションされていない粒子がある場合、前のステップに戻り、全ての粒子のシミュレーションが完了するまで、シミュレーションされていない粒子をシミュレーションし続け、カウントするステップと、を含む。
本発明の第4態様に係る放射線照射システムの制御方法において、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、前記放射線照射システムの制御方法は、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得するステップを含む。
一実施例において、CPUのプロセス並列、スレッド並列及びGPU加速により粒子をシミュレーションする。
好ましくは、プロセス並列及びスレッド並列はいずれも、マルチコアCPUを使用して並列計算を実現する。
一実施例において、CPUのプロセス並列及びスレッド並列の計算過程は、まず、システムがプロセス数又はスレッド数を取得して、値nを取得するステップと、次に、システムが、シミュレーションすべき粒子をn個に均一に分けるステップと、そして、各スレッド又は各プロセスにおいてそれぞれ1個の粒子を独立してシミュレーションし、かつカウントするステップと、最後に、システムが各プロセス又は各スレッドにおいて取得されたカウント数を統計して、最終的な線量を取得するステップと、を含む。
好ましくは、GPU加速は、GPUのマルチコアプロセッサの並列計算により実現される。
GPU加速によるシミュレーション計算の過程は、まず、システムが乱数、断面積データなどをCPUメモリからGPUグラフィックメモリに伝送し、次に、GPU中の各プロセッサが単一の粒子をシミュレーションし、計算し、カウントし、かつカウント結果をグローバルカウントに統計するステップと、そして、システムが、シミュレーションされていない粒子があるか否かを判断し、シミュレーションされていない粒子がない場合、カウント結果をGPUメモリからCPUメモリに伝送し、シミュレーションされていない粒子がある場合、前のステップに戻り、全ての粒子のシミュレーションが完了するまで、シミュレーションされていない粒子をシミュレーションし続け、カウントするステップと、を含む。
本発明の実施例に係る放射線照射システム及びその制御方法は、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得することにより、計算速度を大幅に向上させ、計算時間を減少させることができる。
本発明の第5態様に係る放射線照射システムは、治療用ビームを発生させ、かつ前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、光子の半吸収厚さが第1所定の値以下である場合、光子に対するシミュレーションを停止する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含む。
放射線照射システムは、好ましくは、中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される。
一実施例において、前記治療計画モジュールは、式(1-1)で光子の半吸収厚さtを計算し、
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される。
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される。
第1所定の値は、1つの細胞のサイズであり、好ましくは、第1所定の値は、0.2mmであり、光子エネルギーが第2所定の値以下である場合、光子エネルギーに対応する光子の半吸収厚さは、第1所定の値以下であり、第2所定の値は、10KeVである。
本発明の第6態様に係る放射線照射システムの制御方法において、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、前記放射線照射システムの制御方法は、光子の半吸収厚さが第1所定の値以下である場合、光子に対するシミュレーションを停止するステップを含む。
一実施例において、前記治療計画モジュールは、式(1-1)で光子の半吸収厚さtを計算し、
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される。
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される。
第1所定の値は、1つの細胞のサイズである。
好ましくは、第1所定の値は、0.2mmである。
光子エネルギーが第2所定の値以下である場合、光子エネルギーに対応する光子の半吸収厚さは、第1所定の値以下である。
好ましくは、第2所定の値は、10KeVである。
本発明の実施例に係る放射線照射システム及びその制御方法は、光子の半吸収厚さが第1所定の値以下である場合、光子に対するシミュレーションを停止することにより、計算精度を保証する前提で、計算時間を減少させることができる。
本発明の第7態様に係る放射線照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、分散減少により粒子をシミュレーションする治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含む。
さらに、分散減少は、暗黙吸収と、ウェイトウインドウゲームと、ロシアンルーレット及びスプリッティングとを含み、放射線照射システムは、中性子捕捉療法システムであり、ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、中性子発生装置は、加速器及びターゲットを含み、加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつターゲットと作用して中性子線を発生させ、ビーム整形体は、中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、中性子発生装置により発生した中性子線は、ビーム整形体を通過して治療台上の被照射体に照射される。
本発明の第8態様に係る放射線照射システムの制御方法において、ビーム照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含み、放射線照射システムの制御方法は、治療計画モジュールが分散減少により粒子をシミュレーションするステップを含む。
さらに、分散減少は、暗黙吸収と、ウェイトウインドウゲームと、ロシアンルーレット及びスプリッティングとを含み、分散減少により粒子をシミュレーションするステップは、1つの線源粒子を取得するステップS1と、粒子が格子セル内で衝突するか否かを判断し、衝突すると判断した場合、ステップS3、S4を順に実行し、衝突しないと判断した場合、ステップS5、S6を順に実行するステップS2と、暗黙吸収処理を行うステップS3と、重みがウェイトウインドウより低いか否かを判断し、低いと判断した場合、ステップS7を実行し、低くないと判断した場合、ステップS2に戻る、ステップS4と、ロシアンルーレット及びスプリッティング処理を行うステップS5と、ロシアンルーレット処理を行うか否かを判断し、ロシアンルーレット処理を行うと判断した場合、ステップS7を実行し、ロシアンルーレット処理を行わないと判断した場合、ステップS2に戻る、ステップS6と、消滅させるか否かを判断し、消滅させると判断した場合、ステップS8を実行し、消滅させないと判断した場合、ステップS9を実行した後にステップS2に戻る、ステップS7と、粒子の処理が完了するか否かを判断し、完了すると判断した場合、フローを終了し、完了しないと判断した場合、ステップS1に戻る、ステップS8と、重みを消滅確率で割るステップS9と、を含む。
さらに、ロシアンルーレット及びスプリッティングは、各格子の格子空間インポータンスを計算し、かつ記録するステップS1と、粒子を取得するステップS2と、粒子に対してウェイトウインドウ検出及び操作を行うステップS3と、粒子が格子境界を跨ぐ前後の格子空間インポータンスIn及びIn+1を計算するステップS4と、InがIn+1より大きいか否かを比較し、大きくない場合、ステップS6を実行した後にステップS3に戻り、大きい場合、ステップS7を実行する、ステップS5と、粒子をスプリッティングし、かつ粒子の重みを低下させるステップS6と、粒子を消滅させるか否かを判断し、粒子を消滅させると判断した場合、ステップS9を実行し、粒子を消滅させないと判断した場合、ステップS8を実行した後にステップS3に戻る、ステップS7と、粒子の重みを増加させるステップS8と、粒子のシミュレーションが完了するか否かを判断し、シミュレーションが完了していないと判断した場合、ステップS2に戻り、シミュレーションが完了したと判断した場合、終了するステップS9と、を含む。
さらに、暗黙吸収は、粒子を取得するステップS1と、粒子が格子セル内で衝突するか否かを判断し、衝突すると判断した場合、ステップS3を実行し、衝突しないと判断した場合、ステップS1に戻る、ステップS2と、重みに散乱発生確率をかけるステップS3と、粒子の重みが最も低い重みより小さいか否かを判断し、小さいと判断した場合、ステップS5を実行し、小さくないと判断した場合、ステップS2に戻る、ステップS4と、粒子を消滅させるか否かを判断し、粒子を消滅させると判断した場合、ステップS7を実行し、粒子を消滅させないと判断した場合、ステップS6を実行した後にステップS2に戻る、ステップS5と、重みを消滅確率で割るステップS6と、粒子のシミュレーションが完了するか否かを判断し、シミュレーションが完了していないと判断した場合、ステップS1に戻り、シミュレーションが完了したと判断した場合、終了するステップS7と、を含む。
さらに、ウェイトウインドウゲームは、粒子運動をシミュレーションするステップS1と、粒子の重みがウェイトウインドウ範囲内にあるか否かを判断し、範囲内にあると判断した場合、ステップS1に戻り、範囲内にないと判断した場合、ステップS3を実行するステップS2と、重みがウェイトウインドウより大きいか否かを判断し、大きいと判断した場合、ステップS4を実行した後にステップS1に戻り、大きくないと判断した場合、ステップS5を実行するステップS3と、粒子をスプリッティングし、重みを低下させるステップS4と、消滅させるか否かを判断し、消滅させると判断した場合、終了し、消滅させないと判断した場合、ステップS6を実行した後にステップS1に戻る、ステップS5と、粒子の重みを増加させるステップS6と、を含む。
一実施例において、格子空間インポータンスは、下式の随伴輸送方程式を解くことにより取得され、
式中、Φ*は、随伴束であり、S*は、随伴線源であり、vは、粒子運動の速度であり、Ωは、粒子運動の方向であり、Σtは、粒子と物質とが衝突する反応断面積であり、Σsは、散乱断面積であり、rは、粒子の所在する位置であり、Eは、粒子のエネルギーであり、tは、時間である。
式中、Φ*は、随伴束であり、S*は、随伴線源であり、vは、粒子運動の速度であり、Ωは、粒子運動の方向であり、Σtは、粒子と物質とが衝突する反応断面積であり、Σsは、散乱断面積であり、rは、粒子の所在する位置であり、Eは、粒子のエネルギーであり、tは、時間である。
より具体的には、粒子は、格子空間インポータンスがInの格子から格子空間インポータンスがIn+1の格子に運動し、In+1>Inであれば、m=In+1/Inにし、粒子をm個の粒子にスプリッティングし、各粒子の重みを元の重みの1/mに低下させ、In+1<Inであれば、粒子に対してロシアンルーレット手法を行い、P=In+1/Inにし、0~1から1つの乱数xをサンプリングし、xがPより小さければ、粒子が生き残り、重みに1/Pをかけ、そうでなければ、粒子を消滅させ、粒子のシミュレーションを終了する。
他の実施例において、粒子の重みは、上限が10で、下限が0.25であり、粒子の重みwが10より大きい場合、wの整数部分をw1にし、小数部分をw2にし、0~1から1つの乱数xをサンプリングし、xがw2より小さい場合、w=w1+1にし、xがw2より大きい場合、w=w1にし、次に、1つの粒子をw個の粒子にスプリッティングして、各粒子の重みが1に低下するまでシミュレーションを行う。
本発明の実施例に係る放射線照射システム及びその制御方法は、分散減少により粒子をシミュレーションすることにより、計算精度を保証する前提で、計算時間を減少させることができる。
本発明の第9態様に係る放射線照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、非一様格子を使用してシミュレーション計算を行う治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含む。
さらに、治療計画モジュールは、重要領域の格子の寸法を非重要領域の格子の寸法より小さく設定する。
さらに、治療計画モジュールは、腫瘍の所在する領域の格子の寸法を0.4mm以下に設定する。
さらに、治療計画モジュールは、血液、空気、骨の所在する領域の格子の寸法を1.6mm以上に設定する。
さらに、治療計画モジュールは、正常な筋肉の所在する領域の格子の寸法を0.8mmより大きく1.6mmより小さく設定する。
放射線照射システムは、好ましくは、中性子捕捉療法システムであり、ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、中性子発生装置は、加速器及びターゲットを含み、加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつターゲットと作用して中性子線を発生させ、ビーム整形体は、中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、中性子発生装置により発生した中性子線は、ビーム整形体を通過して治療台上の被照射体に照射される。
本発明の第10態様に係る放射線照射システムの制御方法において、ビーム照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含み、放射線照射システムの制御方法は、治療計画モジュールが非一様格子を使用してシミュレーション計算を行うステップを含む。
該制御方法は、治療計画モジュールが重要領域の格子の寸法を非重要領域の格子の寸法より小さく設定するステップをさらに含む。
該制御方法は、治療計画モジュールが腫瘍の所在する領域の格子の寸法を0.4mm以下に設定するステップをさらに含む。
該制御方法は、治療計画モジュールが、血液、空気、骨の所在する領域の格子の寸法を1.6mm以上に設定し、正常な筋肉の所在する領域の格子の寸法を0.8mmより大きく1.6mmより小さく設定するステップをさらに含む。
本発明の実施例に係る放射線照射システム及びその制御方法は、非一様格子を使用してシミュレーション計算を行うことにより、計算時間を顕著に増加させない前提で、重要領域の計算精度を向上させ、かつ非重要領域の計算正確度を満たすとともに計算時間を減少させることができる。
以下、図面を参照しながら本発明の実施例をさらに詳細に説明することにより、当業者であれば、明細書の文字を参照して実施することができる。
中性子捕捉療法は、効果的ながん治療の手段として、近年ではその適用が増加しており、そのうち、ホウ素中性子捕捉療法が最も一般的なものとなった。ホウ素中性子捕捉療法に使用される中性子は、原子炉又は加速器によって供給することができる。ホウ素中性子捕捉療法(Boron Neutron Capture Therapy、BNCT)は、ホウ素(10B)含有薬物が熱中性子に対して大きい捕捉断面積を有するという特性を利用し、10B(n,α)7Li中性子捕捉と核分裂反応により、4He及び7Liという2種の重荷電粒子を生成し、2種の重荷電粒子は、平均エネルギーが約2.33MeVであり、高い線エネルギー付与(Linear Energy Transfer、LET)及び短い飛程という特徴を持ち、α粒子の線エネルギーと飛程は、それぞれ150KeV/μm、8μmであり、7Li重荷電粒子の場合は、それぞれ175KeV/μm、5μmであり、2種の粒子の総飛程が約1つの細胞のサイズに相当するため、生体への放射線損傷が細胞レベルに抑えられる。ホウ素含有薬物が腫瘍細胞に選択的に集まり、適切な中性子ビーム源と合わせることで、正常組織に大きな損傷を与えない前提で、腫瘍細胞を正確に殺すという目的を達成することができる。
ホウ素中性子捕捉療法の中性子源が原子炉又は荷電粒子とターゲットとの原子核反応によるものに関わらず、生成するのはすべて混合放射線場であり、即ち、ビームは、低エネルギーから高エネルギーまでの中性子及び光子を含み、深部腫瘍のホウ素中性子捕捉療法について、熱外中性子を除き、他の放射線の含有量が多ければ多いほど、正常組織での非選択的線量沈着を引き起こす割合も大きくなるため、これらの不必要な線量沈着を引き起こす放射線をできる限り低減する必要がある。エアビームの品質要素の他、中性子による人体における線量分布をさらに理解するために、本発明の実施例は、人間の頭部組織の人工器官を使用して線量分布の計算を行い、そして人工器官におけるビームの品質要素を中性子ビーム設計の参考とする。
国際原子力機構(IAEA)は、臨床ホウ素中性子捕捉療法に使用される中性子源に対して、エアビームの品質要素に関する5提案を出しており、この5提案は異なる中性子源の利点と欠点を比較するために利用できる他、中性子生成経路の選定及びビーム整形体の設計をする時の参考として利用できる。この5提案はそれぞれ以下のとおりである。
熱外中性子束Epithermal neutron flux>1x109n/cm2s
高速中性子汚染Fast neutron contamination<2×10-13Gy-cm2/n
光子汚染Photon contamination<2×10-13Gy-cm2/n
熱中性子と熱外中性子束との比thermal to epithermal neutron flux ratio<0.05
中性子流とフラックスとの比epithermal neutron current to flux ratio>0.7
注:熱外中性子エネルギー領域は、0.5eV~40KeVであり、熱中性子エネルギー領域は、0.5eVより小さく、高速中性子エネルギー領域は、40KeVより大きい。
高速中性子汚染Fast neutron contamination<2×10-13Gy-cm2/n
光子汚染Photon contamination<2×10-13Gy-cm2/n
熱中性子と熱外中性子束との比thermal to epithermal neutron flux ratio<0.05
中性子流とフラックスとの比epithermal neutron current to flux ratio>0.7
注:熱外中性子エネルギー領域は、0.5eV~40KeVであり、熱中性子エネルギー領域は、0.5eVより小さく、高速中性子エネルギー領域は、40KeVより大きい。
図1に示すように、本実施例における放射線照射システムは、ホウ素中性子捕捉療法システム100であり、中性子ビーム照射装置10、治療計画モジュール20及び制御モジュール30を含む。中性子ビーム照射装置10は、治療用中性子ビームNを発生させ、治療用中性子ビームNを被照射体200に照射して被照射部位を形成する。治療計画モジュール20は、中性子ビーム照射装置10により発生した治療用中性子ビームNのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成し、治療計画に基づいて照射治療時の中性子ビーム照射装置10に対する被照射部位の位置及び対応する照射時間を決定する。治療計画に基づいて決定された位置に被照射体200を位置決めした後、治療を開始することができ、制御モジュール30は、治療計画モジュール20から現在の被照射体200に対応する治療計画を呼び出し、かつ治療計画に基づいて中性子ビーム照射装置10の照射を制御する。制御モジュール30は、他のデータ情報、例えば、中性子ビーム照射装置10のデータ、被照射体200のデータなどを受信してもよい。
図2に示すように、該実施例において、中性子ビーム照射装置10は、中性子発生装置11、ビーム整形体12、コリメータ13及び治療台14を含み、中性子発生装置11は、加速器111及びターゲットTを含み、加速器111は、荷電粒子(例えば、陽子、デューテリウム原子核など)を加速し、陽子線のような荷電粒子線Pを発生させ、荷電粒子線Pは、ターゲットTに照射されて、ターゲットTと作用して中性子線(中性子ビーム)Nを発生させ、ターゲットTは、好ましくは、金属ターゲットである。必要な中性子収率及びエネルギー、提供可能な加速荷電粒子のエネルギー及び電流、並びに金属ターゲットの物理的・化学的特性などにより、適切な原子核反応が選択され、一般的に検討されている原子核反応は、7Li(p,n)7Be及び9Be(p,n)9Bであり、これらの2種の反応は、いずれも吸熱反応である。2種の原子核反応は、エネルギー閾値がそれぞれ1.881MeVと2.055MeVであり、ホウ素中性子捕捉療法の理想的な中性子源がkeVエネルギーレベルの熱外中性子であるため、理論的には、エネルギーが閾値よりわずかに高い陽子を金属リチウムターゲットに衝撃させることで、比較的低いエネルギーの中性子を生成することができ、あまり多くの減速処理を必要とせずに臨床的に使用することができるが、金属リチウム(Li)及び金属ベリリウム(Be)の2種のターゲットは、閾値エネルギーの陽子と作用する断面が高くなく、十分に大きな中性子束を発生させるために、一般的に比較的高いエネルギーを持つ陽子を選択して原子核反応を引き起こす。理想的なターゲットは、中性子収率が高く、生成した中性子のエネルギー分布が熱外中性子エネルギー領域(以下、詳細に説明する)に近く、強い透過性を有する放射線があまり多く発生せず、安全で、安価で、操作しやすく、かつ耐高温などの特性を有するが、実際には全ての要件を満たす原子核反応を見つけることは不可能であり、本発明の実施例においてリチウム金属で製造されたターゲットを使用する。当業者に周知のように、ターゲットTは、リチウム、ベリリウム以外の金属材料で製造されてもよく、例えば、タンタル(Ta)又はタングステン(W)などで形成され、ターゲットTは、円板状であってもよく、他の固体形状であってもよく、液体材料(液体金属)を使用してもよい。加速器111は、線形加速器、サイクロトロン、シンクロトロン、シンクロサイクロトロンであってもよく、中性子発生装置11は、加速器及びターゲットを使用せず、原子炉であってもよい。ホウ素中性子捕捉療法の中性子源が原子炉又は加速器による荷電粒子とターゲットとの原子核反応によるものに関わらず、生成するのは実際にすべて混合放射線場であり、即ち、ビームは、低エネルギーから高エネルギーまでの中性子及び光子を含む。深部腫瘍のホウ素中性子捕捉療法について、熱外中性子を除き、他の放射線の含有量が多ければ多いほど、正常組織での非選択的線量沈着を引き起こす割合も大きくなるため、これらの不必要な線量を引き起こす放射線をできる限り低減する必要がある。また、被照射体の正常組織について、様々な放射線が多すぎて、同様に不必要な線量沈着を引き起こすことを回避すべきである。
中性子発生装置11により発生した中性子ビームNは、順にビーム整形体12とコリメータ13を通過して治療台14上の被照射体200に照射される。ビーム整形体12は、中性子発生装置11により発生した中性子ビームNのビーム品質を調整可能であり、コリメータ13は、中性子ビームNを集め、治療中に中性子ビームNに高い標的性を備えさせる。理解できるように、本発明は、コリメータを有さずに、ビームがビーム整形体12から出た後に治療台14上の被照射体200に直接的に照射されてもよい。
ビーム整形体12は、反射体121、減速体122、熱中性子吸収体123、放射線遮蔽体124及びビーム出口125をさらに含み、中性子発生装置11により発生した中性子のエネルギースペクトルが広く、治療需要を満たす熱外中性子の以外、他の種類の中性子及び光子の含有量をできる限り減少させて、操作者又は被照射体に損傷を引き起こすことを回避する必要があるため、中性子発生装置10から出た中性子は、減速体22を通して高速中性子のエネルギー(>40keV)が熱外中性子エネルギー領域(0.5eV~40keV)に調整され、かつ熱外中性子(<0.5eV)をできる限り減少させる必要があり、減速体22は、高速中性子と作用する断面が大きく、熱外中性子と作用する断面が小さい材料で製造され、好ましい実施例として、減速体122は、D2O、AlF3、Fluental(登録商標)、CaF2、Li2CO3、MgF2及びAl2O3のうちの少なくとも1種で製造され、反射体121は、減速体122を囲み、かつ減速体122を通過して周囲に拡散した中性子を中性子ビームNに反射して中性子の利用率を向上させ、中性子反射能力が高い材料で製造され、好ましい実施例として、反射体121は、Pb又はNiのうちの少なくとも1種で製造され、減速体122は、後部に熱中性子吸収体123があり、熱中性子と作用する断面が大きい材料で製造され、好ましい実施例として、熱中性子吸収体123は、Li-6で製造され、減速体122を通過した熱中性子を吸収して中性子ビームNにおける熱中性子の含有量を減少させ、治療時に浅層正常組織に対して多すぎる線量を与えることを回避し、理解できるように、熱中性子吸収体は、減速体と一体であってもよく、減速体の材料にLi-6が含まれ、放射線遮蔽体124は、ビーム出口125以外の部分からしみ出る中性子及び光子を遮蔽し、放射線遮蔽体124の材料は、光子遮蔽材料と中性子遮蔽材料のうちの少なくとも1種を含み、好ましい実施例として、放射線遮蔽体124の材料は、光子遮蔽材料の鉛(Pb)と中性子遮蔽材料のポリエチレン(PE)を含む。コリメータ13は、ビーム出口125の後部に設置され、コリメータ13から出た熱外中性子ビームは、被照射体200に照射され、浅層正常組織を通過した後に熱中性子に減速されて腫瘍細胞Mに到着する。理解できるように、ビーム整形体20は、さらに他の構造であってもよく、治療に必要な熱外中性子ビームを取得すればよく、説明しやすくするために、コリメータ13が設置された場合、コリメータ13の出口は、後述するビーム出口125と見なすこともできる。本実施例において、被照射体200とビーム出口125との間に放射線遮蔽装置15がさらに設置され、ビーム出口125から出たビームの被照射体の正常組織に対する放射線を遮蔽し、理解できるように、放射線遮蔽装置15が設置されなくてもよい。
被照射体200は、ホウ素(B-10)含有薬物を服用するか又は注射された後、ホウ素含有薬物が腫瘍細胞Mに選択的に集まり、次にホウ素(B-10)含有薬物が熱中性子に対して高い捕捉断面積を有するという特性を利用して、10B(n,α)7Li中性子捕捉及び核分裂反応により、4He及び7Liという2種の重荷電粒子を生成する。2種の荷電粒子は、平均エネルギーが約2.33MeVであり、高い線エネルギー付与(Linear Energy Transfer、LET)及び短い飛程という特徴を持ち、α粒子の線エネルギー付与と飛程は、それぞれ150keV/μm、8μmであり、7Li重荷電粒子の線エネルギー付与と飛程は、175keV/μm、5μmであり、2種の粒子の総飛程が約1つの細胞のサイズに相当するため、生体への放射線損傷が細胞レベルに抑えられ、正常組織に大きな損傷を与えない前提で、腫瘍細胞を部分的に殺すという目的を達成することができる。
ホウ素中性子捕捉療法システム100は、コンクリート構造の建築物内に全体的に収容され、具体的には、ホウ素中性子捕捉療法システム100は、照射室101及び荷電粒子ビーム発生室102をさらに含み、治療台14上の被照射体200は、照射室101において中性子ビームNの照射による治療を受け、荷電粒子ビーム発生室102は、加速器111を少なくとも部分的に収容し、ビーム整形体12は、照射室101と荷電粒子ビーム発生室102との仕切り壁103内に少なくとも部分的に収容される。理解できるように、仕切り壁103は、照射室101と荷電粒子ビーム発生室102とを完全に隔てるものであってもよく、照射室101と荷電粒子ビーム発生室102とを部分的に隔てるものであってもよく、照射室101と荷電粒子ビーム発生室102とが連通する。ターゲットTは、1つ又は複数であってもよく、荷電粒子線Pは、選択的にそのうちの1つ又はいくつかのターゲットTと作用するか又は同時に複数のターゲットTと作用して、1つ又は複数の治療用中性子ビームNを発生させる。ターゲットTの数に応じて、ビーム整形体12、コリメータ13及び治療台14は、1つ又は複数であってもよく、複数の治療台は、同一の照射室内に設置されてもよく、各治療台に個別の照射室が設置されてもよい。照射室101と荷電粒子ビーム発生室102は、コンクリート壁W(仕切り壁103を含む)で囲まれた空間であり、コンクリート構造は、ホウ素中性子捕捉療法システム100の動作過程において漏れた中性子及び他の放射線を遮蔽することができる。
がん細胞を最大限に殺すとともに正常組織への放射線損傷を低減できるために、治療計画モジュール20の設定において、熱外中性子及び光子の線量分布の精度が非常に重要になる。放射線照射システムの適用シーンにおいて、線量計算は、光子及び中性子断面積データベースに読み込まれ、中性子の断面積データベースが非常に大きいため、放射線照射システムの治療計画を作成するソフトウェアのインストールは、大きい記憶空間を占め、かつ該ソフトウェアを使用して線量計算シミュレーションを行う場合に占められたメモリ空間が大きく、かつ計算時間が長くなる。本発明の以下の実施例は、一連の治療計画モジュールの最適化方法を提供し、実行メモリを低下させるとともに線量計算時間を減少させて、放射線治療計画スキームを迅速に作成する要求を満たす。以下、図面を参照しながら各最適化方法を説明する。
治療計画モジュール20においてデータベースを最適化する。
人体内に合計60種あまりの元素が含有され、そのうち、20種あまりの元素が必須要素であり、人体の正常な生理学的機能を維持することに重要な意義を有する。人体内の含有量が高い元素は、炭素、水素、酸素、窒素、リン、塩素、ナトリウム、マグネシウム、カリウム、カルシウムなどであり、炭素、水素、酸素、窒素は、人体の有機物を構成する主要な元素であり、人体総重量の96%を占め、人体内の0.01%より大きい残りの大量元素は、それぞれカルシウム、カリウム、リン、硫黄、塩素、マグネシウム、ナトリウムであり、それらが人体総重量を占める百分率は、それぞれ1.5%、0.35%、1%、0.25%、0.15%、0.05%、0.15%である。以上より、炭素、水素、酸素、窒素、カルシウム、カリウム、リン、硫黄、塩素、マグネシウム、ナトリウムの含有量は、既に人体総重量の99.45%に達し、モンテカルロに100種あまりの元素のシミュレーション及び計算システムが記憶され、毎回、モンテカルロを使用してシミュレーション及び計算を行う場合、いずれも各元素をシミュレーション及び計算するため、かかる時間が長い。また、多くの元素に対応する計算及びシミュレーションシステムが記憶されているため、モンテカルロのデータベースが非常に大きい。しかしながら、放射線照射システムの適用シーンにおいて、モンテカルロにより人体内の光子、中性子などの粒子に対してシミュレーション計算を行う場合、人体総重量を占める百分率が大きい(0.01%より大きい)人体内の元素に対してシミュレーション計算を行うだけで正確な計算結果を取得することができ、いくつかの微量元素に対するシミュレーション計算は、長い時間がかかり、該計算結果は、人体総重量の万分の一以上を占める元素に対してシミュレーション計算を行って取得された結果に比べて、差値が0.001%以内にある。
また、ある元素、例えば、Feは、人体内の占める重量比が大きいが、中性子及び光子と反応せず、又は中性子及び光子との反応強度が大きくないため、このような元素のシミュレーション計算結果による総体的な計算結果への影響を無視してもよい。ある元素、例えば、ホウ素は、人体内の占める重量比が小さいが、放射線照射システムの適用シーンにおいて人体に大量のホウ素薬物を注射し、かつ中性子との反応が激しく、最終的な計算結果への影響が大きい。
元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおける、中性子及び光子のシミュレーション計算結果に対して影響力を有する元素(H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaから選択される1つ又は複数の元素)をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションすることにより、計算精度を保証する前提でシミュレーション計算時間を減少させて、モンテカルロの実行効率を向上させることができ、また、スクリーニングされたこれらの元素に対応するデータベース容量は、元のデータベース容量の2%~3%のみであり、磁気ディスクのメモリ要求が大幅に低下し、製造コストが低減される。
一方では、各元素に対応する中性子断面積データベースにおいて、モンテカルロシステムに、選択のために、多くの温度、例えば、0K、1200K、2500K、250K、294K、600K、900Kに対応する断面積データが記憶されているが、放射線照射システムの適用シーンにおいて、人体モデルのみに対してシミュレーション計算を行うだけでよく、人体の正常温度が310K~315Kである。モンテカルロの実行原理に基づいて、温度が入力された後、システムは、記憶されている多くの温度において、入力された温度に最も近い温度を自動的にマッチングし、かつ該温度をシミュレーションパラメータとしてシミュレーション計算を行うため、人体の温度が入力された後、システムは、データベースに記憶されている温度294Kをマッチングし、かつ294Kをパラメータとしてシミュレーション計算を行う。また、ドップラー効果の計算時に、使用される温度は、0Kである。つまり、放射線照射システムの適用シーンにおいて、データベースに記憶されている温度のうち、294K及び0Kの2つのみが使用され、他の温度は、放射線照射システムの適用シーンに対して余分であるため、294K及び0K以外の他の温度に対応するデータベースを除去し、シミュレーション計算構造に影響を与えない前提で、データベースのサイズを約元の三分の一に減少させ、データベースの実行コストが低減される。
表1は、それぞれ294K、310K及び330Kの温度で、モンテカルロシステムが同一の人体モデルに対して中性子線量シミュレーション計算を行った結果である。
表1のデータから分かるように、温度の小範囲の変動による中性子線量の最終的な計算結果への影響を無視してもよく、放射線照射システムの適用シーンにおいて、システムに記憶されている294Kに対応するデータベースを保留するだけで、使用要求を満たすことができ、かつ最終的な計算結果の誤差が許容可能な範囲内にある。
粒子のシミュレーション過程においてプロセスを遮断して粒子のシミュレーション時間を減少させる。
放射線照射システムの適用シーンにおいて、線量計算は、中性子線量計算及び光子線量計算を含む。光子は、被照射体の体内の運動過程において、光子エネルギー及び光子の半吸収厚さが徐々に低下するが、光子の吸収断面積が光子エネルギー及び光子の半吸収厚さの低下につれて迅速に増大する。光子の半吸収厚さが1つの細胞のサイズより小さい場合、光子がその所在する細胞内に吸収される確率は、非常に大きく、この場合、光子がその所在する細胞内に全てのエネルギーが沈着することを直接的に設定し、そのシミュレーションを継続せず、このような処理後に算出した線量分布と光子のシミュレーションを継続して算出した線量分布との誤差が0.1%以内にある。
人体の細胞のサイズが2~200μmであり、通常、放射線照射システムの適用において、モデルの格子の最小寸法が0.8mmであると、光子の半吸収厚さが格子の最小寸法の四分の一より小さい場合、光子のシミュレーションを停止することができる。
具体的には、式(1-1)で光子の半吸収厚さtを計算し、
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される。
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される。
以下、人体の骨格を例として、異なる光子エネルギー値の骨格での半吸収厚さを計算する。
表2から分かるように、骨格において、光子エネルギーが10KeVより小さい場合、光子の半吸収厚さが0.2mmより小さく、つまり、光子エネルギーが10KeVより小さい場合、光子がその現在所在する格子内に吸収される可能性が非常に大きいため、この場合、光子のシミュレーションを停止して取得された線量と、光子をさらにシミュレーションして取得された線量との差が許容可能な範囲内にある。
光子の細胞組織での半吸収厚さが、光子が通過する材料及び光子エネルギーにより決定されるため、各種の材料に対応するカットオフエネルギーを設定することができ、光子エネルギーがこのカットオフエネルギーより低い場合、光子のシミュレーションを停止することができる。以下、人体モデルにおいてシミュレーション計算を行うことにより、設定された異なる光子カットオフエネルギーに対応する光子線量を取得する。
表2から分かるように、光子エネルギーが10KeV以下である場合、光子のシミュレーション計算を停止し、最終的に算出した光子線量と、光子のシミュレーション計算を継続して取得された線量分布との誤差が0.1%以内にある。
以上より、光子のシミュレーション計算の過程において、光子エネルギーが所定のカットオフエネルギーより小さい場合、光子の半吸収厚さが1つの細胞のサイズより小さく、この場合、光子のシミュレーションを停止し、線量分布の計算結果に影響を与えない場合に光子線量の計算時間を減少させることができる。
他の実施例において、実際のシミュレーション精度の要求に応じて第1所定の値及び第2所定の値を設定し、光子の半吸収厚さが第1所定の値以下である場合、或いは、光子エネルギーが第2所定の値以下である場合、光子のシミュレーションを停止し、第1所定の値は、1つの細胞のサイズより大きくてもよく、1つの細胞のサイズより小さくてもよく、第2所定の値は、12KeV、16KeVなどであってもよい。
並列技術を使用して線量計算モジュールを加速する。
放射線照射システムの治療計画モジュールは、モンテカルロ法で、異なる線源粒子のシミュレーション過程が完全に独立するものであり、つまり、異なる粒子のシミュレーション順序が計算結果に影響を与えず、このような計算特性に基づいて、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得する。コンピュータは、CPUのプロセス並列、スレッド並列という2種の並列及びGPU加速を提供することができる。
図3に示すように、いわゆるCPUのプロセス並列及びスレッド並列の計算過程は、まず、システムがプロセス数又はスレッド数を取得して、値nを取得するステップと、次に、システムがシミュレーションすべき粒子をn個に均一に分けるステップと、そして、各スレッド又は各プロセスにおいてそれぞれ1個の粒子を独立してシミュレーションし、かつカウントするステップと、最後に、システムが各プロセス又は各スレッドにおいて取得されたカウント数を統計して、最終的な線量を取得するステップと、を含む。粒子を異なるプロセス又はスレッドに割り当てるとともにシミュレーションし計算し、シミュレーション及び計算の時間を、1つのプロセス又はスレッドを使用して全ての粒子をシミュレーションし計算することに必要な時間のn分の1(プロセス)又は2n分の1に減少させる。
プロセス並列及びスレッド並列はいずれも、マルチコアCPUを使用して並列計算を実現し、GPU加速は、GPUのマルチコアプロセッサの並列計算により実現される。プロセス及びスレッド並列効果は、CPUのコス数によって限定され、一般的なスタンドアローンの場合、コア数が4つであり、このように、プロセス並列及びスレッド並列は、速度を最大限4倍(プロセス)又は8倍(スレッド)に向上させ、1つのGPUは、複数のプロセッサが集積され、論理的に計算速度は、複数倍向上することができる。
図4に示すように、いわゆるGPU加速によるシミュレーション計算の過程は、まず、システムが乱数、断面積データなどをCPUメモリからGPUグラフィックメモリに伝送し、次に、GPU中の各プロセッサが単一の粒子をシミュレーションし、計算し、カウントし、かつカウント結果をグローバルカウントに統計するステップと、そして、システムが、シミュレーションされていない粒子があるか否かを判断し、シミュレーションされていない粒子がない場合、カウント結果をGPUメモリからCPUメモリに伝送し、シミュレーションされていない粒子がある場合、前のステップに戻り、全ての粒子のシミュレーションが完了するまで、シミュレーションされていない粒子をシミュレーションし続け、カウントするステップと、を含む。
GPUのプロセッサが多いため、GPU加速によるシミュレーション計算は、計算速度を大幅に向上させ、計算時間を減少させることができる。また、CPUは、価格が約6000元のIntel Xeon processor with 2.27GHzを使用し、GPUは、価格が約9000元のNVIDIA Tesla C2050を使用し、各GPUは、合計448個のプロセッサを有し、該448個のプロセッサは、448個のCPUに相当する。CPUの計算時間は、GPU計算時間の50~70倍である。これにより、CPUに比べて、GPUを使用してシミュレーション計算を行うと、価格及びかかる時間の面で大きい優位性を有する。
モンテカルロの分散減少法を使用して加速及び収束を行う。
モンテカルロ法は、プログラムの計算速度を向上させるために豊かな分散減少法を有し、治療計画システムに使用可能な分散減少法は、暗黙吸収、ウェイトウインドウゲーム、ロシアンルーレット及びスプリッティングなどを含み、ロシアンルーレット及びスプリッティングは、ロシアンルーレット手法及びスプリッティング手法で構成される。以下、分散減少法のロシアンルーレット手法、スプリッティング手法、格子空間インポータンス、暗黙吸収及びウェイトウインドウゲームを解釈する。
ロシアンルーレット手法について、通常、粒子の重みは、上限が10に設定され、下限が0.25に設定される。粒子の重みwがある所定の値より小さく(例えば、0.25)に低下する場合、0~1から1つの乱数xをサンプリングし、xがwより小さい場合、粒子が生き残り、粒子の重みが1に戻り、xがw以上であれば、粒子を消滅させ、粒子のシミュレーションを終了する。
スプリッティング手法について、粒子の重みwがある値(例えば、10)より大きい場合、wの整数部分をw1にし、小数部分をw2にし、0~1から1つの乱数xをサンプリングし、xがw2より小さい場合、w=w1+1にし、xがw2より大きい場合、w=w1にし、次に、1つの粒子をw個の粒子にスプリッティングして、各粒子の重みが1に低下するまでシミュレーションを行う。
格子空間インポータンスについて、モデルの異なる領域にある粒子による線量の寄与は、差異があり、空間インポータンスで粒子による線量の寄与を表し、シミュレーションを行う前に、腫瘍の位置、モデル特性などのパラメータに基づいて、各格子の格子空間インポータンスを計算する。粒子は、格子空間インポータンスがInの格子から格子空間インポータンスがIn+1の格子に移動し、In+1>Inであれば、m=In+1/Inにさせ、粒子をm個の粒子にスプリッティングし、各粒子の重みを元の1/mに低下させ、In+1<Inであれば、粒子に対してロシアンルーレット手法を行い、P=In+1/Inにし、0~1から1つの乱数xをサンプリングし、xがPより小さければ、粒子が生き残り、重みに1/Pをかけ、そうでなければ、粒子を消滅させ、粒子のシミュレーションを終了する。
格子空間インポータンスの計算方法について、格子空間インポータンスは、随伴束を計算することにより取得され、随伴束は、下式の随伴輸送方程式を解くことにより取得され、
式中、Φ*は、随伴束であり、S*は、随伴線源であり、vは、粒子運動の速度であり、Ωは、粒子運動の方向であり、Σtは、粒子と物質とが衝突する反応断面積であり、Σsは、散乱断面積であり、rは、粒子の所在する位置であり、Eは、粒子のエネルギーであり、tは、時間である。
式中、Φ*は、随伴束であり、S*は、随伴線源であり、vは、粒子運動の速度であり、Ωは、粒子運動の方向であり、Σtは、粒子と物質とが衝突する反応断面積であり、Σsは、散乱断面積であり、rは、粒子の所在する位置であり、Eは、粒子のエネルギーであり、tは、時間である。
暗黙吸収について、粒子と物質とが相互作用する場合、粒子は、散乱するしかなく、吸収されず、衝突するたびに粒子の重みにP散乱/P総をかける。粒子の重みがある値(例えば、0.25)に低下する場合、ロシアンルーレット手法で粒子を処理する。
ウェイトウインドウゲームについて、粒子の重みがある値、例えば、10より大きい場合、スプリッティング手法を行い、粒子の重みがある値、例えば、0.25より小さい場合、ロシアンルーレット手法を行う。ここでの2つの値は、ソフトウェア性能に基づいて設定される。
まとめて言えば、粒子は、格子空間インポータンスが小さい箇所から格子空間インポータンスが大きい箇所に運動する場合、粒子の重みが向上し、粒子が消滅せず、重みが低下する場合について、粒子は、格子空間インポータンスが大きい箇所から格子空間インポータンスが小さい箇所に運動する場合、粒子の重みが低下し、粒子に対してスプリッティング又は暗黙吸収処理を行う。熱中性子の人体内での吸収断面積(主に、N、Bなどの元素)が大きいため、中性子源から離れた箇所の格子の線量計算収束速度を加速するために、暗黙吸収により粒子がこれらの格子に運動するべきである。粒子が治療領域から離れる場合、粒子のシミュレーションを継続することは、治療領域の線量計算に対する意義が大きくなく、この場合、計算を継続すると、計算リソースを浪費するため、格子空間インポータンスを設定することにより、このような状況の発生確率を低減することができる。粒子の重みが小さすぎる場合、粒子によるカウントの寄与が非常に小さく、シミュレーションを継続すると、計算リソースを浪費し、粒子の重みが大きすぎる場合、単一のカウントが大きすぎることを引き起こし、間違った結果が生じるリスクがあるため、ウェイトウインドウを設定して粒子の重みを適切な範囲にするべきである。したがって、暗黙吸収と組み合わせて、格子空間インポータンス及びウェイトウインドウゲームを合理的に設定することによりモンテカルロプログラムの収束速度を加速する必要がある。
図5に示すように、分散減少法で粒子をシミュレーションするフローは、1つの線源粒子を取得するステップS1と、粒子が格子セル内で衝突するか否かを判断し、衝突すると判断した場合、ステップS3、S4を順に実行し、衝突しないと判断した場合、ステップS5、S6を順に実行するステップS2と、暗黙吸収処理を行うステップS3と、重みがウェイトウインドウより低いか否かを判断し、低いと判断した場合、ステップS7を実行し、低くないと判断した場合、ステップS2に戻る、ステップS4と、ロシアンルーレット及びスプリッティング処理を行うステップS5と、ロシアンルーレット処理を行うか否かを判断し、ロシアンルーレット処理を行うと判断した場合、ステップS7を実行し、ロシアンルーレット処理を行わないと判断した場合、ステップS2に戻る、ステップS6と、消滅させるか否かを判断し、消滅させると判断した場合、ステップS8を実行し、消滅させないと判断した場合、ステップS9を実行した後にステップS2に戻る、ステップS7と、粒子の処理が完了するか否かを判断し、完了すると判断した場合、フローを終了し、完了しないと判断した場合、ステップS1に戻る、ステップS8と、重みを消滅確率で割るステップS9と、を含む。
図6に示すように、ロシアンルーレット及びスプリッティングのフローは、各格子の格子空間インポータンスを計算し、かつ記録するステップS1と、粒子を取得するステップS2と、粒子に対してウェイトウインドウ検出及び操作を行うステップS3と、粒子が格子境界を跨ぐ前後の格子空間インポータンスIn及びIn+1を計算するステップS4と、InがIn+1より大きいか否かを比較し、大きくない場合、ステップS6を実行した後にステップS3に戻り、大きい場合、ステップS7を実行する、ステップS5と、粒子をスプリッティングし、かつ粒子の重みを低下させるステップS6と、粒子を消滅させるか否かを判断し、粒子を消滅させると判断した場合、ステップS9を実行し、粒子を消滅させないと判断した場合、ステップS8を実行した後にステップS3に戻る、ステップS7と、粒子の重みを増加させるステップS8と、粒子のシミュレーションが完了するか否かを判断し、シミュレーションが完了していないと判断した場合、ステップS2に戻り、シミュレーションが完了したと判断した場合、終了するステップS9と、を含む。
図7に示すように、暗黙吸収のフローは、粒子を取得するステップS1と、粒子が格子セル内で衝突するか否かを判断し、衝突すると判断した場合、ステップS3を実行し、衝突しないと判断した場合、ステップS1に戻る、ステップS2と、重みに散乱発生確率をかけるステップS3と、粒子の重みが最も低い重みより小さいか否かを判断し、小さいと判断した場合、ステップS5を実行し、小さくないと判断した場合、ステップS2に戻る、ステップS4と、粒子を消滅させるか否かを判断し、粒子を消滅させると判断した場合、ステップS7を実行し、粒子を消滅させないと判断した場合、ステップS6を実行した後にステップS2に戻る、ステップS5と、重みを消滅確率で割るステップS6と、粒子のシミュレーションが完了するか否かを判断し、シミュレーションが完了していないと判断した場合、ステップS1に戻り、シミュレーションが完了したと判断した場合、終了するステップS7と、を含む。
図8に示すように、ウェイトウインドウゲームのフローは、粒子運動をシミュレーションするステップS1と、粒子の重みがウェイトウインドウ範囲内にあるか否かを判断し、範囲内にあると判断した場合、ステップS1に戻り、範囲内にないと判断した場合、ステップS3を実行するステップS2と、重みがウェイトウインドウより大きいか否かを判断し、大きいと判断した場合、ステップS4を実行した後にステップS1に戻り、大きくないと判断した場合、ステップS5を実行するステップS3と、粒子をスプリッティングし、重みを低下させるステップS4と、消滅させるか否かを判断し、消滅させると判断した場合、終了し、消滅させないと判断した場合、ステップS6を実行した後にステップS1に戻る、ステップS5と、粒子の重みを増加させるステップS6と、を含む。
分散減少法による線量計算への加速効果を比較するために、暗黙吸収を使用する場合と暗黙吸収を使用しない場合との計算時間及び結果の標準偏差を比較する。同一の人体モデルに対して、1千万個の粒子を計算し、暗黙吸収手法を使用する場合の計算時間が1584sであり、線量計算の標準偏差が11%であり、暗黙吸収手法を使用しない場合の計算時間が1258sであり、線量計算の標準偏差が14.5%であり、標準偏差を暗黙吸収手法を使用する場合の標準偏差と同じになるように低下させるために、シミュレーションされる粒子数を1600万個に増加させ、シミュレーション計算時間を2045sに増加させる必要があり、暗黙吸収手法を使用する場合の計算時間より461s増加し、すなわち、暗黙吸収法を使用すると、計算時間が約20%低下する。
非一様矩形格子を使用してシミュレーション計算を行う。
従来の人体モデルの計算は、一様格子を使用することが多く、ある場合に部分領域の精密な計算のために、格子の数を増加させなければならないため、計算時間が長くなり、実行メモリが指数関数的に増加することになる。CT又はPETから読み取られたモデルにより取得されたものは、一様格子であることが多く、多くの接続領域の材料は同じであり、例えば、空気、血液などであり、これらの領域のうちのいくつかに対して、細かい線量分布を正確に計算する必要がないため、粗い格子で細い格子を代替することができ、腫瘍などの重要な位置に対して、より精密な計算を実現するように、より細い格子を分ける必要がある。非一様格子を使用して線量を計算し、計算精度を向上させるとともに、計算時間及び実行メモリが顕著に増加せず、計算時間を顕著に増加させない場合に重要領域の計算精度を向上させ、かつ非重要領域の計算正確度を満たす場合に計算時間を減少させることができる。
以下、異なる寸法の格子を設定してシミュレーション計算を行い、異なる格子寸法によるシミュレーション時間及び計算精度への影響を比較する。
本実験において、使用された混合格子は、中性子の入射方向にあり、前の5cmの部分は、0.4mmの格子で構成され、中間の5cmの部分は、0.8mmの格子で構成され、残りの部分は、1.6mmの格子で構成され、表3から分かるように、非一様格子を使用することにより、格子の総数が元の0.4mmの格子の数の2/5に低下し、格子に必要なメモリも元の1/3に低下し、計算時間が元の34.7%に減少することができ、0.4mmの格子を使用して計算して取得された結果を基準として、混合格子を使用して計算して取得された結果における中性子の誤差が0.1%より小さく、光子の誤差が0.2%より小さい。
混合格子の構成パターンは、以上に挙げられたパターンに限定されず、被照射体の具体的な状況に応じて設定することができる。一般的な場合、格子の寸法は、該領域の重要度により決定され、例えば、腫瘍の所在する領域の格子は、寸法が小さく、4mm以下の格子に設定され、血液、空気及び骨の所在する領域の格子は、寸法が大きく、1.6mm以上に設定され、他の正常筋肉などの組織の所在する領域の格子の寸法は、0.8mmより大きく1.6mmより小さく設定される。
線量計算に使用されるアルゴリズムは、モンテカルロアルゴリズムであり、モンテカルロアルゴリズムは、利点として計算精度が高く、欠点として収束速度が遅く、計算時間が長いため、計算効率の最適化は、最も重要な一部である。本願の実施例1~実施例5の最適化方法は、異なる程度で計算効率を向上させ、演算時間を減少させることができる。
以上に本発明の例示的な具体的な実施形態を説明して、当業者が本発明を理解することを容易にするが、明らかに、本発明は、具体的な実施形態の範囲に限定されず、当業者にとって、様々な変化が添付の特許請求の範囲で限定及び決定される本発明の精神及び範囲内にあれば、これらの変化が明らかで、いずれも本発明の特許請求の範囲内にある。
ホウ素中性子捕捉療法は、ホウ素(10B)含有薬物が熱中性子に対して大きい捕捉断面積を持つという特性を利用し、 10
B(n,α)
7
Li中性子捕捉と核分裂反応により、4Heと7Liという2種の重荷電粒子を生成する。2種の粒子の総飛程が約1つの細胞のサイズに相当するため、生体への放射線損傷が細胞レベルに抑えられる。ホウ素含有薬物が腫瘍細胞に選択的に集まり、適切な中性子ビーム源と合わせることで、正常組織に大きな損傷を与えない前提で、腫瘍細胞を部分的に殺すという目的を達成することができる。
本発明の第2態様に係る放射線照射システムの制御方法において、放射線照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含み、放射線照射システムの制御方法は、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションするステップを含む。
本発明の第4態様に係る放射線照射システムの制御方法において、前記放射線照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、前記放射線照射システムの制御方法は、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得するステップを含む。
本発明の第6態様に係る放射線照射システムの制御方法において、前記放射線照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、前記放射線照射システムの制御方法は、光子の半吸収厚さが第1所定の値以下である場合、光子に対するシミュレーションを停止するステップを含む。
本発明の第8態様に係る放射線照射システムの制御方法において、放射線照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含み、放射線照射システムの制御方法は、治療計画モジュールが分散減少により粒子をシミュレーションするステップを含む。
本発明の第10態様に係る放射線照射システムの制御方法において、放射線照射システムは、治療用ビームを発生させ、治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、治療用ビームのパラメータ及び被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、治療計画に基づいてビーム照射装置の照射を制御する制御モジュールと、を含み、放射線照射システムの制御方法は、治療計画モジュールが非一様格子を使用してシミュレーション計算を行うステップを含む。
ビーム整形体12は、反射体121、減速体122、熱中性子吸収体123、放射線遮蔽体124及びビーム出口125をさらに含み、中性子発生装置11により発生した中性子のエネルギースペクトルが広く、治療需要を満たす熱外中性子の以外、他の種類の中性子及び光子の含有量をできる限り減少させて、操作者又は被照射体に損傷を引き起こすことを回避する必要があるため、中性子発生装置11から出た中性子は、減速体122を通して高速中性子のエネルギー(>40keV)が熱外中性子エネルギー領域(0.5eV~40keV)に調整され、かつ熱外中性子(<0.5eV)をできる限り減少させる必要があり、減速体122は、高速中性子と作用する断面が大きく、熱外中性子と作用する断面が小さい材料で製造され、好ましい実施例として、減速体122は、D2O、AlF3、Fluental、CaF2、Li2CO3、MgF2及びAl2O3のうちの少なくとも1種で製造され、反射体121は、減速体122を囲み、かつ減速体122を通過して周囲に拡散した中性子を中性子ビームNに反射して中性子の利用率を向上させ、中性子反射能力が高い材料で製造され、好ましい実施例として、反射体121は、Pb又はNiのうちの少なくとも1種で製造され、減速体122は、後部に熱中性子吸収体123があり、熱中性子と作用する断面が大きい材料で製造され、好ましい実施例として、熱中性子吸収体123は、Li-6で製造され、減速体122を通過した熱中性子を吸収して中性子ビームNにおける熱中性子の含有量を減少させ、治療時に浅層正常組織に対して多すぎる線量を与えることを回避し、理解できるように、熱中性子吸収体は、減速体と一体であってもよく、減速体の材料にLi-6が含まれ、放射線遮蔽体124は、ビーム出口125以外の部分からしみ出る中性子及び光子を遮蔽し、放射線遮蔽体124の材料は、光子遮蔽材料と中性子遮蔽材料のうちの少なくとも1種を含み、好ましい実施例として、放射線遮蔽体124の材料は、光子遮蔽材料の鉛(Pb)と中性子遮蔽材料のポリエチレン(PE)を含む。コリメータ13は、ビーム出口125の後部に設置され、コリメータ13から出た熱外中性子ビームは、被照射体200に照射され、浅層正常組織を通過した後に熱中性子に減速されて腫瘍細胞Mに到着する。理解できるように、ビーム整形体12は、さらに他の構造であってもよく、治療に必要な熱外中性子ビームを取得すればよく、説明しやすくするために、コリメータ13が設置された場合、コリメータ13の出口は、後述するビーム出口125と見なすこともできる。本実施例において、被照射体200とビーム出口125との間に放射線遮蔽装置15がさらに設置され、ビーム出口125から出たビームの被照射体の正常組織に対する放射線を遮蔽し、理解できるように、放射線遮蔽装置15が設置されなくてもよい。
表3から分かるように、光子エネルギーが10KeV以下である場合、光子のシミュレーション計算を停止し、最終的に算出した光子線量と、光子のシミュレーション計算を継続して取得された線量分布との誤差が0.1%以内にある。
Claims (50)
- 放射線照射システムであって、
治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、
前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成し、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、前記放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションする治療計画モジュールと、
前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含むことを特徴とする放射線照射システム。 - 前記治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaから選択される1つ又は複数の元素である、ことを特徴とする請求項1に記載の放射線照射システム。
- 前記治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaである、ことを特徴とする請求項2に記載の放射線照射システム。
- 前記治療計画モジュールは、各元素の対応する中性子断面積データベースから、294K及び0K以外の温度に対応するデータベースを除去する、ことを特徴とする請求項1~3のいずれか一項に記載の放射線照射システム。
- 中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される、ことを特徴とする請求項1~3のいずれか一項に記載の放射線照射システム。
- 放射線照射システムの制御方法であって、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、
前記放射線照射システムの制御方法は、元素の、人体内の重量割合と、中性子及び光子との反応強度とを総合して、放射線照射システムの適用シーンにおいて中性子及び光子のシミュレーション計算結果に対して影響力を有する元素をスクリーニングし、シミュレーション過程において、スクリーニングされた元素のみをシミュレーションするステップを含む、ことを特徴とする放射線照射システムの制御方法。 - 前記治療計画モジュールによってスクリーニングした元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaから選択される1つ又は複数の元素である、ことを特徴とする請求項6に記載の放射線照射システムの制御方法。
- 前記治療計画モジュールによってスクリーニングされた元素は、H、He、Li、Be、B、C、N、O、F、Na、Mg、Al、Si、P、S、Cl、Ar、K及びCaである、ことを特徴とする請求項7に記載の放射線照射システムの制御方法。
- 治療計画モジュールが、各元素の対応する中性子断面積データベースから、294K及び0K以外の温度に対応するデータベースを除去するステップをさらに含む、ことを特徴とする請求項6~8のいずれか一項に記載の放射線照射システムの制御方法。
- 前記放射線照射システムは、中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される、ことを特徴とする請求項6~8のいずれか一項に記載の放射線照射システムの制御方法。
- 治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、
前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得する治療計画モジュールと、
前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含むことを特徴とする放射線照射システム。 - 中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される、ことを特徴とする請求項11に記載の放射線照射システム。
- CPUのプロセス並列、スレッド並列及びGPU加速により粒子をシミュレーションし、プロセス並列及びスレッド並列は、いずれもマルチコアCPUを使用して並列計算を実現し、CPUのプロセス並列及びスレッド並列の計算過程は、まず、システムがプロセス数又はスレッド数を取得して、値nを取得するステップと、次に、システムが、シミュレーションすべき粒子をn個に均一に分けるステップと、そして、各スレッド又は各プロセスにおいてそれぞれ1個の粒子を独立してシミュレーションし、かつカウントするステップと、最後に、システムが各プロセス又は各スレッドにおいて取得されたカウント数を統計して、最終的な線量を取得するステップと、を含む、ことを特徴とする請求項11又は12に記載の放射線照射システム。
- GPU加速は、GPUのマルチコアプロセッサの並列計算により実現され、GPU加速によるシミュレーション計算の過程は、まず、システムが乱数、断面積データなどをCPUメモリからGPUグラフィックメモリに伝送し、次に、GPU中の各プロセッサが単一の粒子をシミュレーションし、計算し、カウントし、かつカウント結果をグローバルカウントに統計するステップと、そして、システムが、シミュレーションされていない粒子があるか否かを判断し、シミュレーションされていない粒子がない場合、カウント結果をGPUメモリからCPUメモリに伝送し、シミュレーションされていない粒子がある場合、前のステップに戻り、全ての粒子のシミュレーションが完了するまで、シミュレーションされていない粒子をシミュレーションし続け、カウントするステップと、を含む、ことを特徴とする請求項13に記載の放射線照射システム。
- 放射線照射システムの制御方法であって、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、
前記放射線照射システムの制御方法は、異なる線源粒子のシミュレーションタスクを異なるプロセス又はスレッドに割り当て、各プロセス又は各スレッドの計算タスクを完了した後にまとめて、最終的な計算結果を取得するステップを含む、ことを特徴とする放射線照射システムの制御方法。 - CPUのプロセス並列、スレッド並列及びGPU加速により粒子をシミュレーションする、ことを特徴とする請求項15に記載の放射線照射システムの制御方法。
- プロセス並列及びスレッド並列はいずれも、マルチコアCPUを使用して並列計算を実現する、ことを特徴とする請求項16に記載の放射線照射システムの制御方法。
- CPUのプロセス並列及びスレッド並列の計算過程は、まず、システムがプロセス数又はスレッド数を取得して、値nを取得するステップと、次に、システムが、シミュレーションすべき粒子をn個に均一に分けるステップと、そして、各スレッド又は各プロセスにおいてそれぞれ1個の粒子を独立してシミュレーションし、かつカウントするステップと、最後に、システムが各プロセス又は各スレッドにおいて取得されたカウント数を統計して、最終的な線量を取得するステップと、を含む、ことを特徴とする請求項17に記載の放射線照射システムの制御方法。
- GPU加速は、GPUのマルチコアプロセッサの並列計算により実現される、ことを特徴とする請求項16に記載の放射線照射システムの制御方法。
- GPU加速によるシミュレーション計算の過程は、まず、システムが乱数、断面積データなどをCPUメモリからGPUグラフィックメモリに伝送し、次に、GPU中の各プロセッサが単一の粒子をシミュレーションし、計算し、カウントし、かつカウント結果をグローバルカウントに統計するステップと、そして、システムが、シミュレーションされていない粒子があるか否かを判断し、シミュレーションされていない粒子がない場合、カウント結果をGPUメモリからCPUメモリに伝送し、シミュレーションされていない粒子がある場合、前のステップに戻り、全ての粒子のシミュレーションが完了するまで、シミュレーションされていない粒子をシミュレーションし続け、カウントするステップと、を含む、ことを特徴とする請求項19に記載の放射線照射システムの制御方法。
- 治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、
前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、光子の半吸収厚さが第1所定の値以下である場合、光子に対するシミュレーションを停止する治療計画モジュールと、
前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含むことを特徴とする放射線照射システム。 - 中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される、ことを特徴とする請求項21に記載の放射線照射システム。
- 前記治療計画モジュールは、式(1-1)で光子の半吸収厚さtを計算し、
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される、ことを特徴とする請求項21又は22に記載の放射線照射システム。 - 第1所定の値は、1つの細胞のサイズであり、好ましくは、第1所定の値は、0.2mmであり、光子エネルギーが第2所定の値以下である場合、光子エネルギーに対応する光子の半吸収厚さは、第1所定の値以下であり、第2所定の値は、10KeVである、ことを特徴とする請求項21又は22に記載の放射線照射システム。
- 放射線照射システムの制御方法であって、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、
前記放射線照射システムの制御方法は、光子の半吸収厚さが第1所定の値以下である場合、光子に対するシミュレーションを停止するステップを含む、ことを特徴とする放射線照射システムの制御方法。 - 前記治療計画モジュールは、式(1-1)で光子の半吸収厚さtを計算し、
式中、μは、光子の線減衰係数であり、光子が通過する材料及び光子エネルギーにより決定される、ことを特徴とする請求項25に記載の放射線照射システムの制御方法。 - 第1所定の値は、1つの細胞のサイズである、ことを特徴とする請求項26に記載の放射線照射システムの制御方法。
- 第1所定の値は、0.2mmである、ことを特徴とする請求項26に記載の放射線照射システムの制御方法。
- 光子エネルギーが第2所定の値以下である場合、光子エネルギーに対応する光子の半吸収厚さは、第1所定の値以下である、ことを特徴とする請求項25に記載の放射線照射システムの制御方法。
- 第2所定の値は、10KeVである、ことを特徴とする請求項29に記載の放射線照射システムの制御方法。
- 治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、
前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、分散減少により粒子をシミュレーションする治療計画モジュールと、
前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含むことを特徴とする放射線照射システム。 - 前記分散減少は、暗黙吸収と、ウェイトウインドウゲームと、ロシアンルーレット及びスプリッティングとを含み、前記放射線照射システムは、中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される、ことを特徴とする請求項31に記載の放射線照射システム。
- 放射線照射システムの制御方法であって、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、
前記放射線照射システムの制御方法は、前記治療計画モジュールが分散減少により粒子をシミュレーションするステップを含む、ことを特徴とする放射線照射システムの制御方法。 - 前記分散減少は、暗黙吸収と、ウェイトウインドウゲームと、ロシアンルーレット及びスプリッティングとを含み、分散減少により粒子をシミュレーションするステップは、1つの線源粒子を取得するステップS1と、粒子が格子セル内で衝突するか否かを判断し、衝突すると判断した場合、ステップS3、S4を順に実行し、衝突しないと判断した場合、ステップS5、S6を順に実行するステップS2と、暗黙吸収処理を行うステップS3と、重みがウェイトウインドウより低いか否かを判断し、低いと判断した場合、ステップS7を実行し、低くないと判断した場合、ステップS2に戻る、ステップS4と、ロシアンルーレット及びスプリッティング処理を行うステップS5と、ロシアンルーレット処理を行うか否かを判断し、ロシアンルーレット処理を行うと判断した場合、ステップS7を実行し、ロシアンルーレット処理を行わないと判断した場合、ステップS2に戻る、ステップS6と、消滅させるか否かを判断し、消滅させると判断した場合、ステップS8を実行し、消滅させないと判断した場合、ステップS9を実行した後にステップS2に戻る、ステップS7と、粒子の処理が完了するか否かを判断し、完了すると判断した場合、フローを終了し、完了しないと判断した場合、ステップS1に戻る、ステップS8と、重みを消滅確率で割るステップS9と、を含む、ことを特徴とする請求項33に記載の放射線照射システムの制御方法。
- 前記ロシアンルーレット及びスプリッティングは、各格子の格子空間インポータンスを計算し、かつ記録するステップS1と、粒子を取得するステップS2と、粒子に対してウェイトウインドウ検出及び操作を行うステップS3と、粒子が格子境界を跨ぐ前後の格子空間インポータンスIn及びIn+1を計算するステップS4と、InがIn+1より大きいか否かを比較し、大きくない場合、ステップS6を実行した後にステップS3に戻り、大きい場合、ステップS7を実行する、ステップS5と、粒子をスプリッティングし、かつ粒子の重みを低下させるステップS6と、粒子を消滅させるか否かを判断し、粒子を消滅させると判断した場合、ステップS9を実行し、粒子を消滅させないと判断した場合、ステップS8を実行した後にステップS3に戻る、ステップS7と、粒子の重みを増加させるステップS8と、粒子のシミュレーションが完了するか否かを判断し、シミュレーションが完了していないと判断した場合、ステップS2に戻り、シミュレーションが完了したと判断した場合、終了するステップS9と、を含む、ことを特徴とする請求項34に記載の放射線照射システムの制御方法。
- 前記暗黙吸収は、粒子を取得するステップS1と、粒子が格子セル内で衝突するか否かを判断し、衝突すると判断した場合、ステップS3を実行し、衝突しないと判断した場合、ステップS1に戻る、ステップS2と、重みに散乱発生確率をかけるステップS3と、粒子の重みが最も低い重みより小さいか否かを判断し、小さいと判断した場合、ステップS5を実行し、小さくないと判断した場合、ステップS2に戻る、ステップS4と、粒子を消滅させるか否かを判断し、粒子を消滅させると判断した場合、ステップS7を実行し、粒子を消滅させないと判断した場合、ステップS6を実行した後にステップS2に戻る、ステップS5と、重みを消滅確率で割るステップS6と、粒子のシミュレーションが完了するか否かを判断し、シミュレーションが完了していないと判断した場合、ステップS1に戻り、シミュレーションが完了したと判断した場合、終了するステップS7と、を含む、ことを特徴とする請求項34に記載の放射線照射システムの制御方法。
- ウェイトウインドウゲームは、粒子運動をシミュレーションするステップS1と、粒子の重みがウェイトウインドウ範囲内にあるか否かを判断し、範囲内にあると判断した場合、ステップS1に戻り、範囲内にないと判断した場合、ステップS3を実行するステップS2と、重みがウェイトウインドウより大きいか否かを判断し、大きいと判断した場合、ステップS4を実行した後にステップS1に戻り、大きくないと判断した場合、ステップS5を実行するステップS3と、粒子をスプリッティングし、重みを低下させるステップS4と、消滅させるか否かを判断し、消滅させると判断した場合、終了し、消滅させないと判断した場合、ステップS6を実行した後にステップS1に戻る、ステップS5と、粒子の重みを増加させるステップS6と、を含む、ことを特徴とする請求項34に記載の放射線照射システムの制御方法。
- 格子空間インポータンスは、下式の随伴輸送方程式を解くことにより取得され、
式中、Φ*は、随伴束であり、S*は、随伴線源であり、vは、粒子運動の速度であり、Ωは、粒子運動の方向であり、Σtは、粒子と物質とが衝突する反応断面積であり、Σsは、散乱断面積であり、rは、粒子の所在する位置であり、Eは、粒子のエネルギーであり、tは、時間である、ことを特徴とする請求項35に記載の放射線照射システムの制御方法。 - 粒子は、格子空間インポータンスがInの格子から格子空間インポータンスがIn+1の格子に運動し、In+1>Inであれば、m=In+1/Inにし、粒子をm個の粒子にスプリッティングし、各粒子の重みを元の重みの1/mに低下させ、In+1<Inであれば、粒子に対してロシアンルーレット手法を行い、P=In+1/Inにし、0~1から1つの乱数xをサンプリングし、xがPより小さければ、粒子が生き残り、重みに1/Pをかけ、そうでなければ、粒子を消滅させ、粒子のシミュレーションを終了する、ことを特徴とする請求項38に記載の放射線照射システムの制御方法。
- 粒子の重みは、上限が10で、下限が0.25であり、粒子の重みwが10より大きい場合、wの整数部分をw1にし、小数部分をw2にし、0~1から1つの乱数xをサンプリングし、xがw2より小さい場合、w=w1+1にし、xがw2より大きい場合、w=w1にし、次に、1つの粒子をw個の粒子にスプリッティングして、各粒子の重みが1に低下するまでシミュレーションを行う、ことを特徴とする請求項34~39のいずれか一項に記載の放射線照射システムの制御方法。
- 治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、
前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールであって、非一様格子を使用してシミュレーション計算を行う治療計画モジュールと、
前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含むことを特徴とする放射線照射システム。 - 前記治療計画モジュールは、重要領域の格子の寸法を非重要領域の格子の寸法より小さく設定する、ことを特徴とする請求項41に記載の放射線照射システム。
- 前記治療計画モジュールは、腫瘍の所在する領域の格子の寸法を0.4mm以下に設定する、ことを特徴とする請求項41に記載の放射線照射システム。
- 前記治療計画モジュールは、血液、空気、骨の所在する領域の格子の寸法を1.6mm以上に設定する、ことを特徴とする請求項41に記載の放射線照射システム。
- 前記治療計画モジュールは、正常な筋肉の所在する領域の格子の寸法を0.8mmより大きく1.6mmより小さく設定する、ことを特徴とする請求項41に記載の放射線照射システム。
- 中性子捕捉療法システムであり、前記ビーム照射装置は、中性子発生装置、ビーム整形体及び治療台を含み、前記中性子発生装置は、加速器及びターゲットを含み、前記加速器は、荷電粒子を加速して荷電粒子線を発生させ、かつ前記ターゲットと作用して中性子線を発生させ、前記ビーム整形体は、前記中性子発生装置により発生した中性子線を所定のビーム品質に調整可能であり、前記中性子発生装置により発生した中性子線は、前記ビーム整形体を通過して前記治療台上の前記被照射体に照射される、ことを特徴とする請求項41~45のいずれか一項に記載の放射線照射システム。
- 放射線照射システムの制御方法であって、前記ビーム照射システムは、治療用ビームを発生させ、前記治療用ビームを被照射体に照射して被照射部位を形成するビーム照射装置と、前記治療用ビームのパラメータ及び前記被照射部位の医用画像データに基づいて線量シミュレーション計算を行い、かつ治療計画を作成する治療計画モジュールと、前記治療計画に基づいて前記ビーム照射装置の照射を制御する制御モジュールと、を含み、
前記放射線照射システムの制御方法は、前記治療計画モジュールが非一様格子を使用してシミュレーション計算を行うステップを含む、ことを特徴とする放射線照射システムの制御方法。 - 前記治療計画モジュールが重要領域の格子の寸法を非重要領域の格子の寸法より小さく設定するステップをさらに含む、ことを特徴とする請求項47に記載の放射線照射システムの制御方法。
- 前記治療計画モジュールが腫瘍の所在する領域の格子の寸法を0.4mm以下に設定するステップをさらに含む、ことを特徴とする請求項47に記載の放射線照射システムの制御方法。
- 前記治療計画モジュールが、血液、空気、骨の所在する領域の格子の寸法を1.6mm以上に設定し、正常な筋肉の所在する領域の格子の寸法を0.8mmより大きく1.6mmより小さく設定するステップをさらに含む、ことを特徴とする請求項47に記載の放射線照射システムの制御方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011637729 | 2020-12-31 | ||
CN202011637729.1 | 2020-12-31 | ||
PCT/CN2021/140429 WO2022143333A1 (zh) | 2020-12-31 | 2021-12-22 | 放射线照射系统及其控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2024502075A true JP2024502075A (ja) | 2024-01-17 |
Family
ID=82135377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023540520A Pending JP2024502075A (ja) | 2020-12-31 | 2021-12-22 | 放射線照射システム及びその制御方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20230330433A1 (ja) |
EP (1) | EP4272813A1 (ja) |
JP (1) | JP2024502075A (ja) |
CN (5) | CN114681815A (ja) |
TW (3) | TWI830626B (ja) |
WO (1) | WO2022143333A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118098503A (zh) * | 2022-11-28 | 2024-05-28 | 中硼(厦门)医疗器械有限公司 | 蒙特卡罗计算的优化方法及中子捕获治疗系统 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5844241A (en) * | 1996-07-19 | 1998-12-01 | City Of Hope | System and method for determining internal radioactivity and absorbed dose estimates |
WO2002024278A1 (en) * | 2000-09-22 | 2002-03-28 | Numerix Llc | Improved radiation therapy treatment method |
WO2012066498A1 (en) * | 2010-11-18 | 2012-05-24 | Northern Oncology (Pty) Ltd | Brachytherapy seed, methodology and calculating dose of brachytherapy and method of treatment |
US10737121B2 (en) * | 2011-06-27 | 2020-08-11 | Adelphi Technology, Inc. | Neutron source for neutron capture therapy |
CN103955567B (zh) * | 2014-04-10 | 2017-02-08 | 中国科学院近代物理研究所 | 基于gpu的多粒子输运仿真方法 |
CN104376217B (zh) * | 2014-11-20 | 2016-05-18 | 中国科学院合肥物质科学研究院 | 一种基于蒙特卡罗自适应降低方差的辐射屏蔽计算方法 |
WO2016162999A1 (ja) * | 2015-04-09 | 2016-10-13 | 三菱電機株式会社 | ボーラス装置、治療計画装置及び粒子線治療装置 |
JP6849966B2 (ja) * | 2016-11-21 | 2021-03-31 | 東芝エネルギーシステムズ株式会社 | 医用画像処理装置、医用画像処理方法、医用画像処理プログラム、動体追跡装置および放射線治療システム |
CN108295384B (zh) * | 2017-01-11 | 2020-02-28 | 南京中硼联康医疗科技有限公司 | 基于医学影像的组织元素质量比例解构方法及几何模型建立方法 |
CN109985316B (zh) * | 2017-12-29 | 2021-08-20 | 北京连心医疗科技有限公司 | 一种复杂射野的放疗剂量快速计算设备和存储介质 |
CN110310720B (zh) * | 2018-03-20 | 2022-05-10 | 北京连心医疗科技有限公司 | 一种蒙特卡罗剂量计算方法、设备和储存介质 |
CN110310743B (zh) * | 2018-03-20 | 2022-05-06 | 北京连心医疗科技有限公司 | 一种蒙特卡罗剂量计算方法、设备和存储介质 |
CN109011221A (zh) * | 2018-09-04 | 2018-12-18 | 东莞东阳光高能医疗设备有限公司 | 一种剂量引导的中子俘获治疗系统及其操作方法 |
CA3135517C (en) * | 2019-04-17 | 2023-12-19 | Wei-Lin Chen | Neutron capture therapy system |
TW202042864A (zh) * | 2019-05-21 | 2020-12-01 | 睿健生技股份有限公司 | 放射治療系統、控制醫療儀器的控制方法及控制裝置 |
-
2021
- 2021-12-22 JP JP2023540520A patent/JP2024502075A/ja active Pending
- 2021-12-22 WO PCT/CN2021/140429 patent/WO2022143333A1/zh active Application Filing
- 2021-12-22 EP EP21914071.2A patent/EP4272813A1/en active Pending
- 2021-12-23 TW TW112109935A patent/TWI830626B/zh active
- 2021-12-23 TW TW112109932A patent/TWI823800B/zh active
- 2021-12-23 CN CN202111587192.7A patent/CN114681815A/zh active Pending
- 2021-12-23 CN CN202111587219.2A patent/CN114681816A/zh active Pending
- 2021-12-23 CN CN202111587191.2A patent/CN114681814A/zh active Pending
- 2021-12-23 CN CN202111587225.8A patent/CN114681818A/zh active Pending
- 2021-12-23 CN CN202111587222.4A patent/CN114681817A/zh active Pending
- 2021-12-23 TW TW110148435A patent/TWI800177B/zh active
-
2023
- 2023-06-29 US US18/215,900 patent/US20230330433A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
CN114681814A (zh) | 2022-07-01 |
TWI800177B (zh) | 2023-04-21 |
TW202326748A (zh) | 2023-07-01 |
TW202326749A (zh) | 2023-07-01 |
CN114681816A (zh) | 2022-07-01 |
EP4272813A1 (en) | 2023-11-08 |
WO2022143333A1 (zh) | 2022-07-07 |
CN114681817A (zh) | 2022-07-01 |
CN114681818A (zh) | 2022-07-01 |
TWI830626B (zh) | 2024-01-21 |
CN114681815A (zh) | 2022-07-01 |
TW202228147A (zh) | 2022-07-16 |
US20230330433A1 (en) | 2023-10-19 |
TW202325359A (zh) | 2023-07-01 |
TWI823800B (zh) | 2023-11-21 |
TW202326750A (zh) | 2023-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rasouli et al. | Design and optimization of a beam shaping assembly for BNCT based on D–T neutron generator and dose evaluation using a simulated head phantom | |
Muraro et al. | Challenges in Monte Carlo simulations as clinical and research tool in particle therapy: a review | |
Qin et al. | Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy | |
US20230330433A1 (en) | Radiation irradiation system and control method therefor | |
Franciosini et al. | GPU-accelerated Monte Carlo simulation of electron and photon interactions for radiotherapy applications | |
Türkmen et al. | A new method in beam shaping: Multi-Objective Genetic Algorithm method coupled with a Monte-Carlo based reactor physics code | |
Alcocer-Ávila et al. | Cell survival prediction in hadrontherapy with the NanOx biophysical model | |
TWI842437B (zh) | 放射線照射系統及其控制方法 | |
TWI842438B (zh) | 放射線照射系統及其控制方法 | |
RU2820986C1 (ru) | Система радиационного облучения и способ управления ею | |
Çeçen | Monte Carlo radiation transport in external beam radiotherapy | |
Rahmani et al. | Dose calculation and in-phantom measurement in BNCT using response matrix method | |
Li et al. | GPU-based cross-platform Monte Carlo proton dose calculation engine in the framework of Taichi | |
Amato et al. | Applications of the monte carlo method in medical physics | |
Mirzapour et al. | Fast Monte-Carlo photon transport employing GPU-based parallel computation | |
Holmes et al. | Fast Monte Carlo Dose Calculation in Proton Therapy | |
Chang et al. | A Deterministic Dose Calculation Algorithm for Boron Neutron Capture Therapy | |
Douglass | Development of an integrated stochastic radiobiological model for electromagnetic particle interactions in a 4D cellular geometry. | |
Laakkonen | Monte Carlo simulations of spectral properties of clinical photon beams with Geant4 | |
Aarts | Optimisation of Filter Compositions in Boron Neutron Capture Therapy | |
Goorley | A comparison of three gadolinium based approaches to cancer therapy | |
Ivanchenko et al. | 1 The Monte Carlo | |
Satoh et al. | Monte Carlo Simulation Using Japanese Voxel Phantoms to Analyze the Contribution ofParticle Types and Their Energy Distributions to Organ Doses upon External Neutron Exposure | |
Wemple et al. | Modifications to rtt_MC for Fast Neutron Therapy Treatment Planning | |
Nahum et al. | Monte-Carlo and Grid-Based-Deterministic Models for Patient Dose Computation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230810 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230810 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240516 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240528 |