CN112326708A - 基于瞬发伽马光谱的人体组织密度和元素组成重建算法 - Google Patents
基于瞬发伽马光谱的人体组织密度和元素组成重建算法 Download PDFInfo
- Publication number
- CN112326708A CN112326708A CN202011222491.6A CN202011222491A CN112326708A CN 112326708 A CN112326708 A CN 112326708A CN 202011222491 A CN202011222491 A CN 202011222491A CN 112326708 A CN112326708 A CN 112326708A
- Authority
- CN
- China
- Prior art keywords
- hydrogen
- equation
- spectrum
- gamma
- gamma spectrum
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 44
- 239000000203 mixture Substances 0.000 title claims abstract description 8
- 229910052739 hydrogen Inorganic materials 0.000 claims abstract description 29
- 239000001257 hydrogen Substances 0.000 claims abstract description 29
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 7
- 229910052799 carbon Inorganic materials 0.000 claims description 14
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 13
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 13
- 229910052760 oxygen Inorganic materials 0.000 claims description 13
- 239000001301 oxygen Substances 0.000 claims description 13
- 230000005855 radiation Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 3
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 claims description 2
- 229910052791 calcium Inorganic materials 0.000 claims description 2
- 239000011575 calcium Substances 0.000 claims description 2
- 239000003623 enhancer Substances 0.000 claims description 2
- 230000014509 gene expression Effects 0.000 claims description 2
- 150000002431 hydrogen Chemical class 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 claims 1
- LYCAIKOWRPUZTN-UHFFFAOYSA-N Ethylene glycol Chemical compound OCCO LYCAIKOWRPUZTN-UHFFFAOYSA-N 0.000 abstract description 12
- 239000000463 material Substances 0.000 abstract description 12
- 229920003229 poly(methyl methacrylate) Polymers 0.000 abstract description 8
- 239000004926 polymethyl methacrylate Substances 0.000 abstract description 8
- UWJJYHHHVWZFEP-UHFFFAOYSA-N pentane-1,1-diol Chemical compound CCCCC(O)O UWJJYHHHVWZFEP-UHFFFAOYSA-N 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 239000000470 constituent Substances 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 9
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- CSJDCSCTVDEHRN-UHFFFAOYSA-N methane;molecular oxygen Chemical compound C.O=O CSJDCSCTVDEHRN-UHFFFAOYSA-N 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 238000005025 nuclear technology Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000002727 particle therapy Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/22—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
- G01N23/225—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion
- G01N23/2255—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion using incident ion beams, e.g. proton beams
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/07—Investigating materials by wave or particle radiation secondary emission
- G01N2223/081—Investigating materials by wave or particle radiation secondary emission incident ion beam, e.g. proton
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/601—Specific applications or type of materials density profile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/612—Specific applications or type of materials biological material
- G01N2223/6126—Specific applications or type of materials biological material tissue
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明系统研究了不同元素的光谱组成,从理论上推出了组织密度和元素重建的新算法,解决了氢元素的非线性问题。并验证了三种常见碳氢氧不同组合材料的元素组成和密度,它们是聚甲基丙烯酸甲酯,戊二醇和乙二醇。分子的伽马光谱结构可以分解成其组成元素的伽马光谱的线性和非线性的组合,其组合参数可以通过伽马计数的峰值通过最小二乘法求解,进而确定元素组成以及其密度。
Description
技术领域
粒子治疗,核技术,医学成像,蒙特卡洛技术。
背景技术
由于瞬发伽马是射程验证最有前景的方法,它得以广泛的研究。世界各地的研究组研发了不同的伽马相机,有的可以测量一维的束流射程,有的也可以用来三维成像。有些还进行了临床验证测试。
同时蒙特卡洛方法是研究粒子相互作用的有力工具,新的理论方法或者结构设计都可以用蒙特卡洛程序进行模拟验证。
瞬发伽马光谱还可以用来测量人体组织的元素,但这方面的研究还很少,对碳元素和氧元素的初步研究发现具有很好的线性关系。由于氢元素的中子俘获辐射造成的非线性难题,鲜有对氢元素的研究。
本发明系统研究了不同元素的光谱组成,从理论上推出了组织密度和元素重建的新算法,解决了氢元素的非线性问题。并验证了三种常见碳氢氧不同组合材料的元素组成和密度,它们是聚甲基丙烯酸甲酯,戊二醇和乙二醇,分子结构如图1所示。
发明内容
聚甲基丙烯酸甲酯的缩写是PMMA,它的分子式是C5H8O2,当具有一定能量的质子或离子轰击一个PMMA分子,得到伽马光谱同时如果用质子轰击单个元素的原子核,得到光谱和那么分子的伽马光谱和元素的光谱有方程(2)所示的关系。
氢原子核的中子俘获辐射有部分的中子是碳元素产生的,有部分是氧元素产生的,因此方程(2)进一步分解成方程(3)的形式。其中的系数表达式如方程(4)和(5)所示。
其中和代表碳和氧的中子产生贡献率,是氢的增强因子。由于氢的中子俘获只产生2.23MeV的伽马光子,因此可以单独列出。为数学上方便处理,把方程(3)分裂成一个方程组(6)。其中表示样品x的伽马光谱,表示该光谱在2.23MeV处的计数。
当样品中有M摩尔的分子,则元素的摩尔含量为:MC=M*NC, MO=M*NO和MH=M*NH。方程(6)两边都乘以摩尔数,得到方程组(7a) (7b)。
由于光谱是连续性的,为了得到数值解,我们取其中的一些特征峰来求解方程(7),特征峰如图2所示。因此需要对(7a)进行离散化,得到矩阵方程 (10),优化目标是min(Am-B)。
进一步通过测量束流的射程R,得到体积,还可以求解样品的密度。
V=S*R (14)。
如果样品中还有别的元素比如钙,则重复碳和氧的分析,如果增加的元素可以发生中子俘获反应,则重复氢的非线性分析。
算法验证
这里将通过通过一系列蒙特卡洛仿真实验来验证上述算法。用150MeV 的质子轰击一个细长的靶体,长度要比高能量的射程大,一般18cm,半径为2cm。然后用虚拟探测器测量从靶体中出射的所有伽马光子,能量从0-10MeV,能量分辨设置为0.01MeV。如图4所示。
为得到单质的摩尔谱,先用质子轰击一系列碳氢材料组合,得到伽马光谱如图5所示,再除以射程内分子的摩尔数,得到碳元素的摩尔谱,如图6所示。
再用质子轰击一系列氢氧材料组合,得到伽马光子谱如图7所示,再除以射程内分子的摩尔数,得到氧元素的摩尔谱,如图8所示。
从图6和图8的摩尔谱中提取氢的中子俘获辐射伽马峰2.23MeV的摩尔计数,如图9所示,它随着氢含量的增加而非线性增加,然后当氢超过10-12又开始减少。因此可以认为氢含量是2.23MeV光子的增强剂。体现在方程(4)和 (5)的增强因子。
最后用质子轰击样品材料聚甲基丙烯酸甲酯,戊二醇和乙二醇,得到相应的伽马光谱和特征峰,如图10所示。这里把特征峰按能量分成几组:全能段包含1.5~7MeV的特征峰;低能段包含1.5~4MeV的特征峰;高能段包含4~7MeV 的特征峰;中高能段包含2~7MeV的特征峰。通过求解方程(10),每一组都可以求得元素相应的摩尔含量。
三种样品求解的结果汇总到表1-表3中。可以看到中高能段求得的碳氧比和真实值偏差较小。因此建议用中高能段的特征峰来求解。
表1聚甲基丙烯酸甲酯的碳氧含量
表2戊二醇的碳氧含量
表3乙二醇的碳氧含量
通过求解方程(13)和(15)可以得到氢含量以及样品密度。如表4所示。氢含量由于其非线性难题,最后的偏差还是稍大,但对于氢含量较低的分子,结果符合的比较好。对于三个样品材料,求解的密度都符合的比较好。
表4氢含量求解
用重建算法重新绘制伽马光谱,如图11和图12所示。在各个特征峰处符合的都比较好。因此验证本发明的重建算法的有效性。
说明书附图
图1是分子结构图
图2是原理图
图3是PMMA峰值点图
图4是实验几何设置图
图5是碳氢材料组合的伽马光谱图
图6是碳氢材料组合的每摩尔伽马光谱图
图7是氢氧材料组合的伽马光谱图
图8是氢氧材料组合的每摩尔伽马光谱图
图9是氢2.23MeV计数的非线性图
图10是样品材料的伽马光谱以及特征峰图
图11是重建算法得到的伽马光谱0-4MeV图
图12是重建算法得到的伽马光谱4-7MeV图。
Claims (6)
6.根据权利要求1所述的重建算法,其特征在于,如果样品中还有别的元素比如钙,不发生中子俘获,则重复权利要求2到权利要求4中碳和氧的分析求解过程,如果增加的元素可以发生中子俘获反应,则重复则重复权利要求2到权利要求5中氢的非线性分析求解过程。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011127332 | 2020-10-20 | ||
CN2020111273328 | 2020-10-20 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112326708A true CN112326708A (zh) | 2021-02-05 |
CN112326708B CN112326708B (zh) | 2023-10-24 |
Family
ID=74315981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011222491.6A Active CN112326708B (zh) | 2020-10-20 | 2020-11-05 | 基于瞬发伽马光谱的人体组织密度和元素组成重建算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112326708B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DK498979A (da) * | 1978-11-27 | 1980-05-28 | Philips Nv | Fremgangsmaade og apparat til tomografisk undersoegelse af stler rukturer ved skandering med roentgenstraaler eller gammastraa |
AU8343182A (en) * | 1981-05-11 | 1982-11-18 | Australian Atomic Energy Commission | Neutron method of elemental analysis |
FR2618225A1 (fr) * | 1987-07-15 | 1989-01-20 | France Etat Ponts Chaussees | Procede et installation d'analyse par bombardement neutronique d'un flux de materiaux en vrac. |
US20120046867A1 (en) * | 2010-08-20 | 2012-02-23 | Daniel Faber | Portable System for Analyzing and Determining Elemental Composition of Rock Samples |
CN103069268A (zh) * | 2010-07-22 | 2013-04-24 | 于利奇研究中心有限公司 | 利用标准化样本容器的用于确定中子流的中子激发分析 |
US20160097716A1 (en) * | 2014-09-29 | 2016-04-07 | Zyomed Corp. | Systems and methods for blood glucose and other analyte detection and measurement using collision computing |
CN106291656A (zh) * | 2016-08-01 | 2017-01-04 | 汪金龙 | 瞬发伽马成像系统 |
CN108645880A (zh) * | 2018-05-11 | 2018-10-12 | 南京航空航天大学 | 一种大体积样品的能谱解析方法 |
CN109740279A (zh) * | 2019-01-14 | 2019-05-10 | 东华理工大学 | 一种特征频率谱库的电可控中子伽马能谱解析方法 |
-
2020
- 2020-11-05 CN CN202011222491.6A patent/CN112326708B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DK498979A (da) * | 1978-11-27 | 1980-05-28 | Philips Nv | Fremgangsmaade og apparat til tomografisk undersoegelse af stler rukturer ved skandering med roentgenstraaler eller gammastraa |
AU8343182A (en) * | 1981-05-11 | 1982-11-18 | Australian Atomic Energy Commission | Neutron method of elemental analysis |
FR2618225A1 (fr) * | 1987-07-15 | 1989-01-20 | France Etat Ponts Chaussees | Procede et installation d'analyse par bombardement neutronique d'un flux de materiaux en vrac. |
CN103069268A (zh) * | 2010-07-22 | 2013-04-24 | 于利奇研究中心有限公司 | 利用标准化样本容器的用于确定中子流的中子激发分析 |
US20120046867A1 (en) * | 2010-08-20 | 2012-02-23 | Daniel Faber | Portable System for Analyzing and Determining Elemental Composition of Rock Samples |
US20160097716A1 (en) * | 2014-09-29 | 2016-04-07 | Zyomed Corp. | Systems and methods for blood glucose and other analyte detection and measurement using collision computing |
CN106291656A (zh) * | 2016-08-01 | 2017-01-04 | 汪金龙 | 瞬发伽马成像系统 |
CN108645880A (zh) * | 2018-05-11 | 2018-10-12 | 南京航空航天大学 | 一种大体积样品的能谱解析方法 |
CN109740279A (zh) * | 2019-01-14 | 2019-05-10 | 东华理工大学 | 一种特征频率谱库的电可控中子伽马能谱解析方法 |
Non-Patent Citations (3)
Title |
---|
JIN-LONG WANG等: "Pixelatedpromptgammaimagingdetectorforonlinemeasurementof protonbeam:MonteCarlofeasibilitystudybyFLUKA", 《RADIATION DETECTION TECHNOLOGYAND METHODS》, vol. 2, no. 1, pages 4 * |
李刚;吕伟国;崔志文;马俊;王克协;谢荣华;刘继生;吕秀梅;: "套管井偶极弯曲波频散向高频偏移的特性", 地球物理学报, no. 01, pages 358 - 394 * |
李聪;王志朋;马玉刚;龚亚林;张伟;: "MC法模拟质子治疗中靶区定位", 核电子学与探测技术, no. 12, pages 58 - 62 * |
Also Published As
Publication number | Publication date |
---|---|
CN112326708B (zh) | 2023-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Anderson et al. | Inclusive μ-pair production at 150 GeV by π+ mesons and protons | |
Balewski et al. | Total cross section of the reaction pp→ pK+ Λ close to threshold | |
Agakishiev et al. | Study of dielectron production in C+ C collisions at 1AGeV | |
Molnár et al. | Prompt-gamma activation analysis using the k 0 approach | |
Kang et al. | Squared form factors of valence-shell excitations of atomic argon studied by high-resolution inelastic x-ray scattering | |
Liamsuwan et al. | A Monte Carlo track structure simulation code for the full-slowing-down carbon projectiles of energies 1 keV u–1–10 MeV u–1 in water | |
Belkic | Theory of heavy ion collision physics in hadron therapy | |
Kusakabe et al. | Charge-transfer processes in collisions of H+ ions with H 2, D 2, CO, and CO 2 molecules in the energy range 0.2–4.0 keV | |
Chicilo et al. | Dose profiles and x-ray energy optimization for microbeam radiation therapy by high-dose, high resolution dosimetry using Sm-doped fluoroaluminate glass plates and Monte Carlo transport simulation | |
Spyrou et al. | Cross section measurements of (p, γ) reactions on Pd isotopes relevant to the p process | |
CN112326708B (zh) | 基于瞬发伽马光谱的人体组织密度和元素组成重建算法 | |
Carlson et al. | Measurement of neutral cascade production from negative-kaon-hydrogen at 1.8 GeV/c, K− p→ K 0 Ξ 0 | |
Raisbeck et al. | Production Cross Sections of Be Isotopes in C and O Targets Bombarded by 2.8-GeV α Particles: Implications for Factorization | |
Griffin et al. | Double photoionization of Be and Mg atoms using the R-matrix-with-pseudostates method | |
Peterson et al. | Quasifree pion single charge exchange at 500 MeV | |
Nicolanti et al. | Calculation of electron interaction models in N2 and O2 | |
Yamaki et al. | Systematic study of individual charge-changing cross sections of intermediate-energy secondary beams | |
De Paolis et al. | Investigating the E2 nuclear resonance effect in kaonic atoms | |
Cata-Danil et al. | Measurement of cross sections and thick target yields for (α, γ) process on 63Cu | |
Schatz et al. | 17O (n, α) 14C—Bottle-neck for primordial nucleosynthesis? | |
Inokuti et al. | Physics as an element of Radiation Research | |
Guess et al. | The Nd 150 (He 3, t) and Sm 150 (t, He 3) reactions with applications to β β decay of Nd 150 | |
Lee et al. | Development of graphene ion-chamber for radiation dosimetry | |
Hayward | Study of the giant resonances with virtual and real photons | |
Bucurescu et al. | Excited states of the odd-odd nucleus Eu 158 from the (d, α) reaction |
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 |