CN109192273B - 硼中子俘获治疗的束流评估方法、装置、设备及存储介质 - Google Patents

硼中子俘获治疗的束流评估方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN109192273B
CN109192273B CN201811051944.6A CN201811051944A CN109192273B CN 109192273 B CN109192273 B CN 109192273B CN 201811051944 A CN201811051944 A CN 201811051944A CN 109192273 B CN109192273 B CN 109192273B
Authority
CN
China
Prior art keywords
depth
spectrum
neutron
dose
sub
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
Application number
CN201811051944.6A
Other languages
English (en)
Other versions
CN109192273A (zh
Inventor
黄红
冷志涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guoke Neutron Medical Technology Co ltd
Original Assignee
Dongguan Dongyangguang Gaoneng Medical Equipment Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dongguan Dongyangguang Gaoneng Medical Equipment Co ltd filed Critical Dongguan Dongyangguang Gaoneng Medical Equipment Co ltd
Priority to CN201811051944.6A priority Critical patent/CN109192273B/zh
Publication of CN109192273A publication Critical patent/CN109192273A/zh
Application granted granted Critical
Publication of CN109192273B publication Critical patent/CN109192273B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture

Landscapes

  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Surgery (AREA)
  • Urology & Nephrology (AREA)
  • Engineering & Computer Science (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

公开了一种硼中子俘获治疗的束流评估方法、装置、设备及存储介质，所述方法包括：分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；获取待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取待评估光子能谱内各个子能谱区间的注量，作为第二注量；分别将第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，将第二注量代入第三深度剂量公式，以获得第三曲线数据；通过第一曲线数据、第二曲线数据与第三曲线数据，确定优势深度与治疗增益，以评估超热中子束流的治疗效果，有效提高了超热中子束流治疗效果的评估效率，实现了快速调整中子整形装置的有益效果。

Description

硼中子俘获治疗的束流评估方法、装置、设备及存储介质
技术领域
本发明实施例涉及放疗技术领域，尤其涉及一种硼中子俘获治疗的束流评估方法、装置、设备及存储介质。
背景技术
随着现代化医疗技术的快速发展，硼中子俘获治疗在临床医学得到了广泛应用。
BNCT(Boron Neutron Capture Therapy，硼中子俘获治疗)是一种具有“固有”安全性的生物靶向放射治疗模式，通过将具有亲肿瘤组织的无毒的含硼药物注入人体血液，待硼药物富集在肿瘤组织后，采用超热中子照射肿瘤部位。超热中子进入人体组织内，经过慢化成热中子能够与癌组织中的10B核素发生辐射俘获反应，释放出α粒子(4He)和锂粒子(7Li)。由于这些粒子的射程很短(与细胞核的直径差不多)以及LET(Linear EnergyTransfer，传能线密度)值很高，导致绝大多数的α粒子(4He)和锂粒子(7Li)能够沉积在肿瘤细胞中，从而达到在癌细胞内部破坏癌细胞的效果。又由于α粒子(4He)和锂粒子(7Li)的射程很短，因此只能杀死癌细胞而不损伤周围组织，因此其治疗效果优于常规光子放疗(X射线机、医用直线加速器等)和质子放疗。
在BNCT治疗过程中，为了保证病人正常组织所受剂量在安全范围的前提下，尽可能杀死肿瘤组织，病人在治疗前通常会采用TPS(Treatment Planning System，治疗计划系统)治疗系统以实现对超热中子束流的评估，并制定相应的治疗计划。在现有技术中，在质子打靶产生中子后，需要将中子放入中子整形装置以产生一定范围能量的超热中子束流。但是，由于中子整形装置的高复杂度，每次调整中子整形装置均需花费很长时间，导致人力与物力资源严重浪费，效率较低。
发明内容
本发明实施例提供一种硼中子俘获治疗的束流评估方法、装置、设备及存储介质，能够快速评估超热中子束流的治疗效果。
第一方面，本发明实施例提供了一种硼中子俘获治疗的束流评估方法，包括：
分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；
获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；
分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；其中，所述第一曲线数据为所述待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据，所述第二曲线数据为所述待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据，所述第三曲线数据为所述待评估光子能谱内各个子能谱区间的光子在组织各个深度下对应的剂量数据；
通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果。
第二方面，本发明实施例还提供了一种硼中子俘获治疗的束流评估装置，包括：
区间划分模块，用于分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；
注量获取模块，用于获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；
数据获取模块，用于分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；其中，所述第一曲线数据为所述待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据，所述第二曲线数据为所述待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据，所述第三曲线数据为所述待评估光子能谱内各个子能谱区间的光子在人体组织各个深度下对应的剂量数据；
治疗效果评估模块，用于通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果。
第三方面，本发明实施例还提供了一种设备，所述设备包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现本发明任意实施例提供的硼中子俘获治疗的束流评估方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明任意实施例提供的硼中子俘获治疗的束流评估方法。
本发明实施例提供了一种硼中子俘获治疗的束流评估技术方案，分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果，在调整中子整形装置过程中，避免反复通过蒙特卡洛程序评估超热中子束流时间过长的现象，在保证设备输出精准度的同时，实现了快速调整中子整形装置，以使设备尽早使用的有益效果。
附图说明
图1是本发明实施例一中的一种硼中子俘获治疗的束流评估方法的流程图；
图2A是本发明实施例二中的一种硼中子俘获治疗的束流评估方法的流程图；
图2B是本发明实施例二中的一种中子能谱内的子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型后不同成分的深度剂量曲线；
图3A是本发明实施例三中的一种硼中子俘获治疗的束流评估方法的流程图；
图3B是本发明实施例三中的一种超热中子束流在正常组织和肿瘤组织内产生的深度剂量曲线示意图；
图4是本发明实施例四中的一种硼中子俘获治疗的束流评估装置的结构图；
图5是本发明实施例五中的一种设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
在相关技术中，硼中子俘获治疗的方法作为当前最新的一种放射线疗法，已经在杀死癌细胞的过程中表现出较好的治疗效果。但是，为了保证治疗的效果和治疗过程中病人的安全，在治疗前必须确定用于治疗的射线束的特性，即测量超热中子能谱以及获取其在体中产生的吸收剂量，以制定合适的治疗计划。但是，现有技术对于评估超热中子束流的时间较长，影响设备的使用，,因此，本发明实施例提供了一种硼中子俘获治疗的剂量快速计算以及超热中子束流评估的方法，用于提高中子整形装置的调整效率，尽早用于医用以帮助癌症患者，具体方法如下：
实施例一
图1为本发明实施例一提供的一种硼中子俘获治疗的束流评估方法的流程图，本实施例可适用于各种硼中子俘获治疗的束流评估的情况，该方法可以由本发明实施例提供的硼中子俘获治疗的束流评估装置来执行，该装置可采用软件和/或硬件的方式实现，该装置可集成在任何提供硼中子俘获治疗的束流评估功能的设备中，例如可以是电脑。如图1所示，具体包括：
S110、分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间。
由于质子打靶后产生的中子束能量范围较宽，包括快中子、热中子、超热中子等，而在放射治疗过程中需要使用超热中子进行照射，超热中子在体内慢化成热中子，热中子与硼元素进行辐射俘获反应，因此将中子束通过中子整形装置(Beam Shaping Assembly，BSA)整合为一定能量范围的超热中子，用以照射肿瘤组织。其中，中子整形装置可以将快中子慢化，低能中子过滤以及光子过滤等，因此，经过中子整形装置处理过的中子束大部分为超热中子，成为超热中子束，并能够控制在一定能量范围内。虽然中子整形装置能够滤除大部分的光子，但是仍然会有部分残留，剩余的光子照射人体组织，同样能够对人体组织产生影响。因此，在硼中子俘获治疗的过程中，需要控制超热中子束流的中子和携带的光子在人体组织中的总剂量不超过人体能接受剂量的最大值，以保证治疗安全。
在本实施例中，待评估中子能谱与待评估光子能谱可以是超热中子束流包含的不同能量的中子和光子。其中，待评估中子能谱可以是对超热中子束流中不同能量中子的统计，例如可以是以横坐标为中子能量，纵坐标为中子个数的曲线；待评估光子能谱可以是对超热中子束流中不同能量光子的统计，例如可以是以横坐标为光子能量，纵坐标为光子个数的曲线。
具体的，待评估中子能谱与待评估光子能谱的获取可以通过放置在中子整形装置出口的中子探测仪和光子探测仪而获得。将获取的待评估中子能谱分成多个子能谱区间，如将获取的待评估中子能谱分成13个子能谱区间；将获取的待评估光子能谱分成多个子能谱区间，如将获取的待评估光子能谱分成12个子能谱区间。
具体例如，待评估中子能谱段可分为13段，分别为1e-9MeV到1e-8MeV，1e-8MeV到1e- 7MeV，1e-7MeV到5e-7MeV，5e-7MeV到1e-6MeV，1e-6MeV到5e-6MeV，5e-6MeV到1e-4MeV，1e-4MeV到5e-4MeV，5e-4MeV到1e-3MeV，1e-3MeV到1e-2MeV，1e-2MeV到1e-1MeV，1e-1MeV到1MeV，1MeV到10MeV，10MeV到20MeV，由于BNCT主要使用的是超热中子，因此在超热中子段划分比较细。待评估光子能谱可分为12段，分别为0.01MeV到0.1MeV，0.1MeV到1MeV，1MeV到3MeV，3MeV到5MeV，5MeV到7MeV，7MeV到9MeV，9MeV到11MeV，11MeV到13MeV，13MeV到15MeV，15MeV到17MeV，17MeV到19MeV，19MeV到21MeV。
需要说明的是，本实施例对待评估中子能谱与待评估光子能谱的获取、以及待评估中子能谱与待评估光子能谱的子能谱区间个数不做限定。
S120、获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量。
其中，注量可为在单位面积下，中子整形装置出口通过的中子数或光子数。第一注量可为在单位面积内，中子整形装置出口通过的各个子能谱区间的中子数；第二注量可为在单位面积内，中子整形装置出口通过的各个子能谱区间的光子数。通过获取第一注量，用于记录待评估中子能谱的各个子能谱区间在单位面积内通过的中子个数；通过获取第二注量，用于记录待评估光子能谱的各个子能谱区间在单位面积内通过的光子个数。
S130、分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据。
其中，所述第一曲线数据为所述待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据，所述第二曲线数据为所述待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据，所述第三曲线数据为所述待评估光子能谱内各个子能谱区间的光子在人体组织各个深度下对应的剂量数据。
由于正常组织与肿瘤组织中的硼浓度不同，导致中子进入正常组织与肿瘤组织中产生的深度剂量数据有所不同。第一深度剂量公式为根据第一响应矩阵与第一注量的关系所得到的待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据。其中，第一响应矩阵为预设中子能谱中各个子能谱区间的单位注量在人体正常组织各个深度造成的剂量率。第二深度剂量公式为根据第二响应矩阵与第二注量的关系所得到的待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据。其中，第二响应矩阵为预设中子能谱中各个子能谱区间的单位注量在人体肿瘤组织各个深度造成的剂量率。第三深度剂量公式为根据第三响应矩阵与第三注量的关系所得到的待评估光子能谱内各个子能谱区间的光子在组织各个深度下对应的剂量数据。其中，第三响应矩阵为预设光子能谱中各个子能谱区间的单位注量在人体组织各个深度造成的剂量率。
需要说明的是，人体组织包含肿瘤组织与正常组织，由于光子在肿瘤组织与正常组织中的运动相同，因此研究光子在人体组织内的运动时，无需区分肿瘤细胞与正常细胞。
通过第一深度剂量公式、第二深度剂量公式以及第三深度剂量公式能够获取人体组织中各个深度所含有的剂量率。
S140、通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果。
其中，优势深度AD指的是在肿瘤组织内束流剂量率等于正常组织最大剂量率时所对应的深度，AD可以用来进行束流评估，AD数值越大，可以表明治疗肿瘤的深度越深。其中，正常组织最大剂量率时ADDR，可以表示正常组织能够接受的最大剂量率。治疗增益TG表示在肿瘤组织内的有效深度前，肿瘤的平均剂量率与正常组织最大剂量率ADDR的比值。治疗增益TG的数据越大，表明超热中子束流的治疗效果越好。
具体的，通过将第一曲线数据、第二曲线数据与第三曲线数据进行处理，可以确定肿瘤组织的深度剂量曲线以及正常组织的深度剂量曲线。在正常组织的深度剂量曲线中，确定正常组织最大剂量率ADDR。在肿瘤组织的深度剂量曲线中，确定与正常组织最大剂量率ADDR数值相同的剂量率所对应的深度，作为优势深度AD。在肿瘤组织的深度剂量曲线中，确定优势深度AD之前的肿瘤的平均剂量率与正常组织最大剂量率ADDR的比值，作为治疗增益TG。结合优势深度AD与治疗增益TG，可以评估超热中子束流的治疗效果。
本发明实施例分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果，在调整中子整形装置过程中，避免反复通过蒙特卡洛程序评估超热中子束流时间过长的现象，在保证设备输出精准度的同时，实现了快速调整中子整形装置，以使设备尽早使用的有益效果。
示例性的，在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度之后，还包括：通过所述正常组织最大剂量率、所述正常组织内的平均硼浓度以及所述正常组织的最大剂量率，确定所述超热中子束流的照射时间。
其中，正常组织内的平均硼浓度可为照射超热中子束流前，检测的人体血液中的硼浓度。正常组织的最大剂量率可为人体正常组织接收的最大剂量。
在本实施例中，还可以通过正常组织最大剂量率ADDR，确定超热中子束流在患者肿瘤组织上的照射时间。计算照射时间的公式如下：
ADDR*B*T＝Dm
其中，B表示正常组织的平均硼浓度，可以从人体血液中检测而得；T表示照射时间；Dm表示处方剂量，即医生给定的正常组织能够接受的最大剂量。由上述公式可知，正常组织最大剂量率越大，患者能接受的照射时间越短。
本发明实施例通过正常组织最大剂量率、正常组织内的平均硼浓度以及正常组织的最大剂量率，确定超热中子束流的照射时间，有助于控制治疗时间。
实施例二
图2A为本发明实施例二提供的一种硼中子俘获治疗的束流评估方法的流程图，本实施例在上述实施例的基础上还包括第一深度剂量公式、第二深度剂量公式以及第三深度剂量公式的构建。
相应的，本实施例的方法包括：
S210、构建第一深度剂量公式、第二深度剂量公式以及第三深度剂量公式。
其中，所述第一深度剂量公式包括：
Figure BDA0001794796750000111
其中，Di表示待评估中子能谱在正常组织深度i处引起的总剂量率；Rim表示在预设中子能谱中，单位注量第m个子能谱区间的中子在正常组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第二深度剂量公式包括：
Figure BDA0001794796750000112
其中，Ei表示待评估中子能谱在肿瘤组织深度i处引起的总剂量率；Oim表示在预定中子能谱中，单位注量第m个子能谱区间的中子在肿瘤组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第三深度剂量公式包括：
Figure BDA0001794796750000121
其中，Fj表示待评估光子能谱在人体组织深度j处引起的总剂量率；Pjm表示在预定光子能谱中，单位注量第m个子能谱区间的光子在人体组织深度j处引起的剂量率；φm表示待评估光子能谱中，第m个子能谱区间内光子的注量。
S220、分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间。
S230、获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量。
S240、分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据。
示例性的，在预设中子能谱中，获取单位注量各个子能谱区间的中子在正常组织或肿瘤组织不同深度处引起的剂量率，包括：将所述预设中子能谱分成多个子能谱区间；依次将所述预设中子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取所述各个子能谱区间在预设深度相对应的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据；基于所述正常组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据，确定所述预设中子能谱中单位注量各个子能谱区间的中子在正常组织的不同深度引起的剂量率；基于所述肿瘤组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据，确定所述预设中子能谱中单位注量各个子能谱区间的中子在肿瘤组织的不同深度引起的剂量率。
其中，预设中子能谱为单位注量能量值范围为10e-9MeV到20MeV的中子能谱。由于在中子整形装置出口检测的待评估中子能谱的能量范围可以处于预设中子能谱的任意范围中，因此使用预设中子能谱构建第一深度剂量公式与第二深度剂量公式，能够涵盖各种能量范围的待评估中子能谱的情况。
其中，蒙特卡洛程序是基于蒙特卡罗方法的用于计算三维复杂几何结构中的中子、光子、电子或者耦合中子/光子/电子输运问题的通用软件包，可以包括MCNP(MonteCarlo N Particle Transport Code)、Geant4(Geometry And Tracking)，FLUKA(FLUktuierende KAsade)，EGS(Electron-Gamma Shower)等程序。在本实施例中，可以采用上述任意一种蒙特卡洛程序模拟不同单能中子或单能光子在人体部位模型中深度剂量分布情况，本实施例对此不作限定。、
其中，人体部位模型可以是在蒙特卡洛程序中构建的关于人体各个部位的模型，例如可以是，人体肝模型、人体脑模型、人体胃模型以及人体肺模型等。
其中，硼深度剂量数据指的是在人体组织中预设深度的硼剂量率分布，热中子深度剂量数据指的是在热中子与人体组织中氮元素反应产生的预设深度的剂量率分布；快中子深度剂量数据指的是在人体组织中预设深度的快中子剂量率分布；光子深度剂量数据指的是在人体组织中预设深度的光子剂量率分布。在本实施例中，硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据来源于10B(n,α)7Li、14N(n,p)14C、1H(n,n’)1H、1H(n,γ)2H等多个反应。硼深度剂量数据可以是以横坐标组织深度，纵坐标为硼深度剂量率的数据；热中子深度剂量数据可以是以横坐标组织深度，纵坐标为热中子剂量率的数据；快中子深度剂量数据可以是以横坐标组织深度，纵坐标为快中子深度剂量率的数据；光子深度剂量数据可以是以横坐标组织深度，纵坐标为光子深度剂量率的数据。
具体的，本实施例可以将预设中子能谱分成多个子能谱区间(例如分为13个子能谱区间)，则依次将预设中子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取各个子能谱区间在预设深度相对应的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据。
需要说明的是，本实施例可以利用MATLAB、Excel等能够实现数据处理的制图工具对蒙特卡洛程序(如MCNP)输出的各种成分深度剂量数据的文件进行数据处理，且本实施例不对制图工具的软件类型进行限定。。如图2B所示，图2B为中子能谱内的一个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型后不同成分的深度剂量曲线。其中，曲线20为1ppm浓度的硼深度-剂量曲线，曲线21为热中子深度-剂量曲线，曲线22为快中子深度-剂量曲线，曲线23为光子深度-剂量曲线。
具体的，若在预设中子能谱分为13个子能谱区间，作为13组，则将每组子能谱区间在正常组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据进行累加。其中，累加过程可以是将同一个深度值条件下的剂量率的数据进行累加，以得到预设中子能谱中单位注量各个子能谱区间的中子在正常组织的预设深度引起的剂量率。其中，剂量率的数据可以直接在蒙特卡洛程序(如MCNP)中获取。
具体例如，预设中子能谱中单位注量各个子能谱区间的中子在正常组织的不同深度引起的剂量率的计算公式可以为：
Dw＝WgDg+WBDB+WNDN+WfnDfn
其中，Dw单位注量各个子能谱区间的中子在正常组织不同深度处引起的剂量率；Dg表示γ剂量，也即光子深度剂量数据中的剂量率；DB为硼元素在正常组织(或肿瘤组织)内中引起的吸收剂量率，也即硼深度剂量数据中的剂量率；DN为热中子与人体组织氮元素反应引起的剂量率，也即热中子深度剂量数据中的剂量率；Dfn为快中子深度剂量数据中的剂量率；Wg、WB、WN以及Wfn为加权因子，其中，Wg为在靶目标肿瘤和正常组织中被认为是1.0，WN和Wfn的数值为3.2，WB的数值与正常组织和肿瘤组织有关，正常组织中的WB数值为1.3。
需要说明的是，通过预设中子能谱中单位注量各个子能谱区间的中子在正常组织的不同深度引起的剂量率的计算公式，可以得到Rim，也即在在预设中子能谱中，单位注量第m个子能谱区间的中子在正常组织深度i处引起的剂量率。
同理，若将每组子能谱区间在肿瘤组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据进行累加，可以得到预设中子能谱中单位注量各个子能谱区间的中子在肿瘤组织的预设深度引起的剂量率。
具体例如，预设中子能谱中单位注量各个子能谱区间的中子在肿瘤组织的不同深度引起的剂量率的计算公式可以为：
Dr＝WgDg+WBDB+WNDN+WfnDfn
其中，Dr为单位注量各个子能谱区间的中子在肿瘤组织不同深度处引起的剂量率；Dg表示γ剂量，也即光子深度剂量数据中的剂量率；WB的数值与正常组织和肿瘤组织有关，肿瘤组织中的WB数值为3.8。
需要说明的是，由于肿瘤组织与正常组织中的硼元素含量不同，因此，加权因子WB的数值在正常组织和肿瘤组织中的有所不同。通过预设中子能谱中单位注量各个子能谱区间的中子在肿瘤组织的不同深度引起的剂量率的计算公式，可以得到Oim，也即在预定中子能谱中，单位注量第m个子能谱区间的中子在肿瘤组织深度i处引起的剂量率。
示例性的，在预设光子能谱中，获取单位注量各个子能谱区间的光子在组织不同深度处引起的剂量率，包括：将所述预设光子能谱分成多个子能谱区间；依次将所述预设光子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取所述预设光子能谱中单位注量各个子能谱区间在预设深度相对应的光子深度剂量数据；基于所述人体组织中的光子深度剂量数据，确定所述预设光子能谱中单位注量各个子能谱区间的光子在人体组织的不同深度引起的剂量率。
其中，预设光子能谱为单位注量能量值范围为0.01MeV到20MeV的光子能谱。由于在光子整形装置出口检测的待评估光子能谱的能量范围可以处于预设光子能谱的任意范围中，因此使用预设光子能谱构建第三深度剂量公式，能够涵盖各种能量范围的待评估光子能谱的情况。
具体的，本实施例可以将预设光子能谱分成多个子能谱区间(例如分为12个子能谱区间)，则依次将预设光子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取各个子能谱区间在预设深度相对应的光子深度剂量数据。基于人体组织中的光子深度剂量数据，确定预设光子能谱中单位注量各个子能谱区间的光子在人体组织的不同深度引起的剂量率，也即Pjm。
本发明实施例通过构建第一深度剂量公式、第二深度剂量公式以及第三深度剂量公式，以评估所述超热中子束流的治疗效果，有效提高了超热中子束流治疗效果的评估效率。
实施例三
图3A为本发明实施例三提供的一种硼中子俘获治疗的束流评估方法的流程图，本实施例在上述实施例的基础上进行优化，提供了优化的通过所述多组第一曲线数据、所述多组第二曲线数据与所述多组第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果的处理方法。
相应的，本实施例的方法包括：
S310、分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间。
S320、获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量。
S330、分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据。
S340、将所述第一曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在正常组织内的深度剂量曲线。
其中，第一曲线总数据为待评估中子能谱中不同子能谱区间的中子在正常组织中产生的剂量率；第三曲线总数据为待评估光子能谱中不同子能谱区间的光子在组织中产生的剂量率。
由于超热中子束流中包含待评估中子能谱和待评估光子能谱，因此，若确定超热中子束流在正常组织内的深度剂量数据，则将第一曲线数据与第三曲线数据进行累加，并按照累加后的数据处理为超热中子束流在正常组织内的深度剂量曲线。其中，第一曲线数据与第三曲线数据进行累加为将同一深度值的条件下，剂量率值进行累加的计算过程。
S350、将所述第二曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在肿瘤组织内的深度剂量曲线。
由于超热中子束流中包含待评估中子能谱和待评估光子能谱，因此，若确定超热中子束流在肿瘤组织内的深度剂量数据，则将第二曲线总数据与第三曲线总数据进行累加，并按照累加后的数据处理为超热中子束流在正常组织内的深度剂量曲线。
S360、在所述正常组织内的深度剂量曲线中，将纵坐标的最大值作为所述正常组织最大剂量率。
具体例如，图3B为超热中子束流在正常组织和肿瘤组织内产生的深度剂量曲线示意图。其中，曲线301为超热中子束流在正常组织内产生的深度剂量曲线；曲线302是超热中子束流在肿瘤组织内产生的深度剂量曲线。因此，正常组织最大剂量率为在曲线301内的纵坐标最大值A点，其对应的纵坐标的数值为ADDR值。
S370、在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度，以及将所述优势深度以内的平均剂量率与所述正常组织最大剂量率的比值，作为所述治疗增益。
如图3B所示，优势深度为在肿瘤组织内的深度剂量曲线302中，与正常组织最大剂量率A点相等的剂量率所对应的B点的深度，其在横坐标对应的数值为AD值。治疗增益为优势深度AD以内的平均剂量率与正常组织最大剂量率ADDR的比值。
需要说明的是，本实施例中曲线301的硼浓度可为10ppm；曲线302的硼浓度可为30ppm。
S380、根据所述优势深度与所述治疗增益，以评估所述超热中子束流的治疗效果。
本发明实施例通过获取超热中子束流在正常组织和肿瘤组织内的深度剂量曲线，在正常组织和肿瘤组织内的深度剂量曲线中分析得到正常组织最大剂量率、优势深度，以提高评估所述超热中子束流的治疗效果的准确性。
实施例四
图4为本发明实施例四提供的一种硼中子俘获治疗的束流评估装置的结构示意图，本实施例可适用于各种硼中子俘获治疗的束流评估的情况，该装置可采用软件和/或硬件的方式实现，该装置可集成在任何提供硼中子俘获治疗的束流评估功能的设备中，例如可以是电脑。如图4所示，具体包括：区间划分模块41、注量获取模块42、数据获取模块43和治疗效果评估模块44。
区间划分模块41，用于分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；
注量获取模块42，用于获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；
数据获取模块43，用于分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；其中，所述第一曲线数据为所述待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据，所述第二曲线数据为所述待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据，所述第三曲线数据为所述待评估光子能谱内各个子能谱区间的光子在组织各个深度下对应的剂量数据；
治疗效果评估模块44，用于通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果。
本发明实施例提供了一种硼中子俘获治疗的装置，通过分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果，在调整中子整形装置过程中，避免反复通过蒙特卡洛程序评估超热中子束流时间过长的现象，在保证设备输出精准度的同时，实现了快速调整中子整形装置，以使设备尽早使用的有益效果。
在上述实施例基础上，所述第一深度剂量公式包括：
Figure BDA0001794796750000211
其中，Di表示待评估中子能谱在正常组织深度i处引起的总剂量率；Rim表示在预设中子能谱中，单位注量第m个子能谱区间的中子在正常组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第二深度剂量公式包括：
Figure BDA0001794796750000212
其中，Ei表示待评估中子能谱在肿瘤组织深度i处引起的总剂量率；Oim表示在预定中子能谱中，单位注量第m个子能谱区间的中子在肿瘤组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第三深度剂量公式包括：
Figure BDA0001794796750000213
其中，Fj表示待评估光子能谱在人体组织深度j处引起的总剂量率；Pjm表示在预定光子能谱中，单位注量第m个子能谱区间的光子在人体组织深度j处引起的剂量率；φm表示待评估光子能谱中，第m个子能谱区间内光子的注量。
在上述实施例基础上，还包括：中子剂量率确定模块。
中子剂量率确定模块，用于将所述预设中子能谱分成多个子能谱区间；依次将所述预设中子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取所述各个子能谱区间在预设深度相对应的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据；基于所述正常组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据，确定所述预设中子能谱中单位注量各个子能谱区间的中子在正常组织的不同深度引起的剂量率；基于所述肿瘤组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据，确定所述预设中子能谱中单位注量各个子能谱区间的中子在肿瘤组织的不同深度引起的剂量率。
在上述实施例基础上，还包括：光子剂量率确定模块。
光子剂量率确定模块，用于将所述预设光子能谱分成多个子能谱区间；依次将所述预设光子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取所述预设光子能谱中单位注量各个子能谱区间在预设深度相对应的光子深度剂量数据；基于所述人体组织中的光子深度剂量数据，确定所述预设光子能谱中单位注量各个子能谱区间的光子在人体组织的不同深度引起的剂量率。
在上述实施例基础上，所述治疗效果评估模块44具体用于：
将所述第一曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在正常组织内的深度剂量曲线；将所述第二曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在肿瘤组织内的深度剂量曲线；在所述正常组织内的深度剂量曲线中，将纵坐标的最大值作为所述正常组织最大剂量率；在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度，以及将所述优势深度以内的平均剂量率与所述正常组织最大剂量率的比值，作为所述治疗增益；根据所述优势深度与所述治疗增益，以评估所述超热中子束流的治疗效果。
在上述实施例基础上，还包括：照射时间确定模块。
照射时间确定模块，用于在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度之后，通过所述正常组织最大剂量率、所述正常组织内的平均硼浓度以及所述正常组织的最大剂量率，确定所述超热中子束流的照射时间。
实施例五
本实施五例提供了一种设备。图5是本发明实施例四提供的一种设备的结构示意图。参见图5，该设备包括：
一个或多个处理器510；
存储器520，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器510执行，使得所述一个或多个处理器510实现如本发明实施例提出的硼中子俘获治疗的束流评估方法。
图5中以一个处理器510为例；处理器510和存储器520可以通过总线或其他方式连接，图5中以通过总线连接为例。
存储器520作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的硼中子俘获治疗的束流评估方法对应的程序指令/模块。处理器510通过运行存储在存储器520中的软件程序、指令以及模块，从而实现上述的硼中子俘获治疗的束流评估方法。
存储器520可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据机顶盒的使用所创建的数据等。此外，存储器520可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器520可进一步包括相对于处理器510远程设置的存储器，这些远程存储器可以通过网络连接至机顶盒。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
本实施例提出的设备与本发明实施例提出的硼中子俘获治疗的束流评估方法属于同一发明构思，未在本实施例中详尽描述的技术细节可参见上述实施例，并且本实施例与上述实施例具有相同的有益效果。
实施例六
本实施例提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明实施例提出的硼中子俘获治疗的束流评估方法。
本实施例提出的计算机可读存储介质与本发明实施例提出的硼中子俘获治疗的束流评估方法属于同一发明构思，未在本实施例中详尽描述的技术细节可参见上述实施例，并且本实施例与上述实施例具有相同的有益效果。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(RandomAccess Memory，RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台机顶盒(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (9)

1.一种硼中子俘获治疗的束流评估方法，其特征在于，包括：
分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；
获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；
分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；其中，所述第一曲线数据为所述待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据，所述第二曲线数据为所述待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据，所述第三曲线数据为所述待评估光子能谱内各个子能谱区间的光子在人体组织各个深度下对应的剂量数据；
通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果；
其中，通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果，包括：
将所述第一曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在正常组织内的深度剂量曲线；将所述第二曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在肿瘤组织内的深度剂量曲线；在所述正常组织内的深度剂量曲线中，将纵坐标的最大值作为所述正常组织最大剂量率；在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度，以及将所述优势深度以内的平均剂量率与所述正常组织最大剂量率的比值，作为所述治疗增益；根据所述优势深度与所述治疗增益，以评估所述超热中子束流的治疗效果。
2.根据权利要求1所述的方法，其特征在于，所述第一深度剂量公式包括：
Figure FDA0003109548560000021
其中，Di表示待评估中子能谱在正常组织深度i处引起的总剂量率；Rim表示在预设中子能谱中，单位注量第m个子能谱区间的中子在正常组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第二深度剂量公式包括：
Figure FDA0003109548560000022
其中，Ei表示待评估中子能谱在肿瘤组织深度i处引起的总剂量率；Oim表示在预定中子能谱中，单位注量第m个子能谱区间的中子在肿瘤组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第三深度剂量公式包括：
Figure FDA0003109548560000023
其中，Fj表示待评估光子能谱在人体组织深度j处引起的总剂量率；Pjm表示在预定光子能谱中，单位注量第m个子能谱区间的光子在人体组织深度j处引起的剂量率；φm表示待评估光子能谱中，第m个子能谱区间内光子的注量。
3.根据权利要求2所述的方法，其特征在于，在预设中子能谱中，获取单位注量各个子能谱区间的中子在正常组织或肿瘤组织不同深度处引起的剂量率，包括：
将所述预设中子能谱分成多个子能谱区间；
依次将所述预设中子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取所述各个子能谱区间在预设深度相对应的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据；
基于所述正常组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据，确定所述预设中子能谱中单位注量各个子能谱区间的中子在正常组织的不同深度引起的剂量率；
基于所述肿瘤组织中的硼深度剂量数据、热中子深度剂量数据、快中子深度剂量数据以及光子深度剂量数据，确定所述预设中子能谱中单位注量各个子能谱区间的中子在肿瘤组织的不同深度引起的剂量率。
4.根据权利要求2所述的方法，其特征在于，在预设光子能谱中，获取单位注量各个子能谱区间的光子在人体组织不同深度处引起的剂量率，包括：
将所述预设光子能谱分成多个子能谱区间；
依次将所述预设光子能谱中的各个子能谱区间的能量值输入蒙特卡洛程序构建的人体部位模型，以获取所述预设光子能谱中单位注量各个子能谱区间在预设深度相对应的光子深度剂量数据；
基于所述人体组织中的光子深度剂量数据，确定所述预设光子能谱中单位注量各个子能谱区间的光子在人体组织的不同深度引起的剂量率。
5.根据权利要求1所述的方法，其特征在于，在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度之后，还包括：
通过所述正常组织最大剂量率、所述正常组织内的平均硼浓度以及所述正常组织的最大剂量率，确定所述超热中子束流的照射时间。
6.一种硼中子俘获治疗的束流评估装置，其特征在于，包括：
区间划分模块，用于分别将超热中子束流中的待评估中子能谱与待评估光子能谱分成多个子能谱区间；
注量获取模块，用于获取所述待评估中子能谱内各个子能谱区间的注量，作为第一注量，以及获取所述待评估光子能谱内各个子能谱区间的注量，作为第二注量；
数据获取模块，用于分别将所述第一注量代入第一深度剂量公式与第二深度剂量公式，以获得第一曲线数据与第二曲线数据，以及，将所述第二注量代入第三深度剂量公式，以获得第三曲线数据；其中，所述第一曲线数据为所述待评估中子能谱内各个子能谱区间的中子在正常组织各个深度下对应的剂量数据，所述第二曲线数据为所述待评估中子能谱内各个子能谱区间的中子在肿瘤组织各个深度下对应的剂量数据，所述第三曲线数据为所述待评估光子能谱内各个子能谱区间的光子在人体组织各个深度下对应的剂量数据；
治疗效果评估模块，用于通过所述第一曲线数据、所述第二曲线数据与所述第三曲线数据，确定优势深度与治疗增益，以评估所述超热中子束流的治疗效果；
其中，所述治疗效果评估模块具体用于：将所述第一曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在正常组织内的深度剂量曲线；将所述第二曲线数据与所述第三曲线数据进行累加，以确定所述超热中子束流在肿瘤组织内的深度剂量曲线；在所述正常组织内的深度剂量曲线中，将纵坐标的最大值作为所述正常组织最大剂量率；在所述肿瘤组织内的深度剂量曲线中，将与所述正常组织最大剂量率相等的剂量率所对应的深度作为所述优势深度，以及将所述优势深度以内的平均剂量率与所述正常组织最大剂量率的比值，作为所述治疗增益；根据所述优势深度与所述治疗增益，以评估所述超热中子束流的治疗效果。
7.根据权利要求6所述的装置，其特征在于，所述第一深度剂量公式包括：
Figure FDA0003109548560000051
其中，Di表示待评估中子能谱在正常组织深度i处引起的总剂量率；Rim表示在预设中子能谱中，单位注量第m个子能谱区间的中子在正常组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第二深度剂量公式包括：
Figure FDA0003109548560000052
其中，Ei表示待评估中子能谱在肿瘤组织深度i处引起的总剂量率；Oim表示在预定中子能谱中，单位注量第m个子能谱区间的中子在肿瘤组织深度i处引起的剂量率；φm表示待评估中子能谱中，第m个子能谱区间内中子的注量；
所述第三深度剂量公式包括：
Figure FDA0003109548560000061
其中，Fj表示待评估光子能谱在人体组织深度j处引起的总剂量率；Pjm表示在预定光子能谱中，单位注量第m个子能谱区间的光子在人体组织深度j处引起的剂量率；φm表示待评估光子能谱中，第m个子能谱区间内光子的注量。
8.一种电子设备，其特征在于，所述设备包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-5中任一所述的硼中子俘获治疗的束流评估方法。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-5中任一所述的硼中子俘获治疗的束流评估方法。
CN201811051944.6A 2018-09-10 2018-09-10 硼中子俘获治疗的束流评估方法、装置、设备及存储介质 Active CN109192273B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811051944.6A CN109192273B (zh) 2018-09-10 2018-09-10 硼中子俘获治疗的束流评估方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811051944.6A CN109192273B (zh) 2018-09-10 2018-09-10 硼中子俘获治疗的束流评估方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN109192273A CN109192273A (zh) 2019-01-11
CN109192273B true CN109192273B (zh) 2021-09-28

Family

ID=64915915

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811051944.6A Active CN109192273B (zh) 2018-09-10 2018-09-10 硼中子俘获治疗的束流评估方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN109192273B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020111085A1 (ja) * 2018-11-26 2020-06-04 国立大学法人筑波大学 ホウ素中性子捕捉療法（ｂｎｃｔ）用治療計画システム、ｂｎｃｔシステム及びｂｎｃｔ
JP7378729B2 (ja) * 2019-12-25 2023-11-14 住友重機械工業株式会社 中性子線測定装置、及び中性子線測定方法
CN113096724B (zh) * 2021-04-09 2022-08-02 中国科学院近代物理研究所 一种四维剂量计算系统及存储介质
WO2024020874A1 (zh) * 2022-07-27 2024-02-01 上海联影医疗科技股份有限公司 一种能谱获取方法、系统和存储介质
CN118486473A (zh) * 2024-07-16 2024-08-13 华硼中子科技(杭州)有限公司 用于bnct的剂量分布位移装置设计方法、存储介质、终端和计算机程序产品

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006071588A (ja) * 2004-09-06 2006-03-16 Toshiba Corp 中性子検出器の寿命診断システム
CN105445777A (zh) * 2015-12-01 2016-03-30 南京航空航天大学 一种硼中子俘获治疗剂量测量方法
CN105930636A (zh) * 2016-03-29 2016-09-07 中北大学 一种自动确定目标函数权重的放射治疗方案优化系统
JP2016214760A (ja) * 2015-05-25 2016-12-22 株式会社東芝 ホウ素中性子捕捉療法用治療装置及びその制御方法
EP3195901A1 (en) * 2016-01-20 2017-07-26 Ion Beam Applications S.A. Method and device for determining an interest of applying a qa procedure to a treatment plan in radiation therapy
CN107224675A (zh) * 2016-03-25 2017-10-03 南京中硼联康医疗科技有限公司 硼中子捕获治疗系统
CN107290774A (zh) * 2017-07-21 2017-10-24 四川瑶天纳米科技有限责任公司 中子剂量快速确定方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006071588A (ja) * 2004-09-06 2006-03-16 Toshiba Corp 中性子検出器の寿命診断システム
JP2016214760A (ja) * 2015-05-25 2016-12-22 株式会社東芝 ホウ素中性子捕捉療法用治療装置及びその制御方法
CN105445777A (zh) * 2015-12-01 2016-03-30 南京航空航天大学 一种硼中子俘获治疗剂量测量方法
EP3195901A1 (en) * 2016-01-20 2017-07-26 Ion Beam Applications S.A. Method and device for determining an interest of applying a qa procedure to a treatment plan in radiation therapy
CN107224675A (zh) * 2016-03-25 2017-10-03 南京中硼联康医疗科技有限公司 硼中子捕获治疗系统
CN105930636A (zh) * 2016-03-29 2016-09-07 中北大学 一种自动确定目标函数权重的放射治疗方案优化系统
CN107290774A (zh) * 2017-07-21 2017-10-24 四川瑶天纳米科技有限责任公司 中子剂量快速确定方法

Also Published As

Publication number Publication date
CN109192273A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109192273B (zh) 硼中子俘获治疗的束流评估方法、装置、设备及存储介质
Gao et al. First demonstration of the FLASH effect with ultrahigh dose rate high-energy X-rays
Patriarca et al. Experimental set-up for FLASH proton irradiation of small animals using a clinical system
Alaei et al. Imaging dose from cone beam computed tomography in radiation therapy
Zhu et al. Proton therapy verification with PET imaging
JP6695981B2 (ja) ビーム照射角度の評価方法
CN109173081B (zh) 硼中子俘获治疗的束流评估方法、装置、设备及存储介质
Mesbahi et al. Effect of wedge filter and field size on photoneutron dose equivalent for an 18 MV photon beam of a medical linear accelerator
JP2020503959A (ja) 医用画像に基づく組織元素質量比率の脱構築方法及び幾何的モデルの構築方法
Olivera et al. Maximum likelihood as a common computational framework in tomotherapy
Allahverdi et al. Monte Carlo estimation of electron contamination in a 18 MV clinical photon beam
US20230347176A1 (en) Radiotherapy system and therapy plan generation method therefor
CN107292075A (zh) 增进放射治疗系统计算效益的方法
AU2024200068A1 (en) Irradiation parameter selection apparatus and usage method thereof and control system comprising said apparatus and usage method thereof
Yang et al. Thoracic Organ Doses and Cancer Risk from Low Pitch Helical 4‐Dimensional Computed Tomography Scans
Chow et al. Dose-volume consistency and radiobiological characterization between prostate IMRT and VMAT plans
Mamballikalam et al. Time motion study to evaluate the impact of flattening filter free beam on overall treatment time for frameless intracranial radiosurgery using Varian TrueBeam® linear accelerator
Penyu et al. Effect of cerrobend thickness on electron transmission in megaelectronvolt (MeV) produced by clinical linear accelerator
Aslian et al. Neutron and photon out-of-field doses at cardiac implantable electronic device (CIED) depths
Arndt et al. Dosimetric calibration and characterization for experimental mouse thoracic irradiation using orthovoltage X rays
US20230111230A1 (en) Radiotherapy system and treatment plan generation method therefor
Alem-Bezoubiria et al. Sensitivity Analysis to Neutron Dose at Patient's Organs According to Some Parameters in a LINAC in the Case of a Bladder Treatment
Blombäck Monte Carlo evaluation of static and dynamic 6FFF treatments-Evaluation of dose distributions calculated with AAA, Acuros XB, and Collapsed Cone (RayStation and DoseCheck)
Sudarsana Luqyana Latifa, Ni Nyoman Ratini, I. Wayan Balik Sudarsana b, Ida Bagus Made Suryatika, I. Gusti Ayu Ratnawati, Anak Agung Ngurah Gunawan and Rozi Irhas b
Yu et al. Sensitivity analysis of biological washout and depth selection for a machine learning based dose verification framework in proton therapy

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230704

Address after: Room 1001, Unit 1, Building 1, No. 1 Yiran Road, Songshanhu Park, Dongguan City, Guangdong Province, 523808

Patentee after: Guoke Neutron Medical Technology Co.,Ltd.

Address before: 523808 No. 1318 Building A1, Guangda West Valley, Songshan Lake High-tech Industrial Development Zone, Dongguan City, Guangdong Province

Patentee before: DONGGUAN DONGYANGGUANG GAONENG MEDICAL EQUIPMENT Co.,Ltd.