中子剂量快速确定方法技术

本发明专利技术公开了中子剂量快速确定方法，解决了现有技术中确定硼中子俘获的三维辐照剂量的时间较长、准确率不高的问题。本发明专利技术包括1)在CPU平台上形成中子辐照的体模数据，进行体模数据的离散化，形成计算体元，并确定三维空间上每个计算体元的物质标号；2)选定用于确定中子束剂量的GPU平台，通过CPU平台向GPU平台传送计算数据；3)在GPU上利用蒙特卡罗程序并行计算三维辐照剂量矩阵D的每个单元的剂量，GPU平台向CPU平台返回辐照剂量矩阵数据D。本发明专利技术实现了将原有在CPU中运行的方式修改为在GPU中运行，极大地提高了运行的速率，并且有效保证了计算的精度。

Fast determination method of neutron dose

The invention discloses a fast determination method of neutron dose, which solves the problems of long time and low accuracy of three dimensional radiation dose determination in boron neutron capture in the prior art. The present invention includes 1) phantom data formation of neutron irradiation on the CPU platform, discrete phantom data, the formation of calculated body element, and calculate the three-dimensional space calculated body element material label; 2) selected for determining neutron beam dose calculation of the GPU platform, GPU platform to transfer data through the CPU platform 3); each unit using parallel computation of 3D dose matrix D Monte Carlo program on the GPU GPU platform to CPU dose, radiation dose return matrix data D. The invention realizes that the original operation mode in CPU is modified to run in GPU, which greatly improves the speed of operation and ensures the accuracy of calculation effectively.

【技术实现步骤摘要】
中子剂量快速确定方法
本专利技术涉及一种中子剂量的确定方法，具体涉及中子剂量快速确定方法。
技术介绍
为了实现硼中子俘获治疗，除了必要的热中子照射设备外，快速准确评估热中子辐照后在人体中的剂量分布也是一个关键环节。迄今为止，硼中子俘获治疗肿瘤的方法已有70多年历史，尽管方法先进、有效，但一直未能取得临床应用的突破，主要技术瓶颈为三个方面：制造适用于医院大产额中子源的困难；研制富集于人体含硼药物的困难；热中子辐照剂量的准确快速确定的困难。硼中子俘获治疗在对一些肿瘤的临床试治中都出现过因辐照过量导致正常组织坏死的情况，也出现了在患部辐照剂量不足形成日后肿瘤的复发的情况。这些情况反映了实施辐照剂量与处方辐照剂量的差距。目前国际上辐照剂量计算软件有:美国的SERA，2000年，由SeppalT，KankaanrantaL，Se.T等人在日本大阪举行的第九届中子俘获治疗癌症国际研讨会上发表；美国的MacNCTPlan，2002年，由SantaCruzSJ，KigerWIII等人在德国埃森市举行的第十届中子俘获治疗癌症国际研讨会上发表；日本的JCDS，2004年，由KumadaH，YamamotoK，MatsumuraA在Phys.Med.Biol.杂志第49卷上发表。上述软件主体计算均采用蒙特卡罗程序MCNP模拟，硼中子俘获治疗的单一网格模型结构简单，由完全相同的立方体网格(或长方体网格)组成，排列有序，随着网格尺寸的变小，网格数多达几十万，蒙特卡罗程序MCNP在模拟硼中子俘获治疗模型时计算粒子射线与网格面的交点要花费大量时间，除此以外，材料的质量密度、原子密度以及网格的六个面号等都是网格编号的函数，这些数组很庞大，搜索时间相当长。
技术实现思路
本专利技术的目的在于解决现有技术中确定硼中子俘获的三维辐照剂量的时间较长、准确率不高的问题，提供一种解决上述问题的中子剂量快速确定方法。本专利技术通过下述技术方案实现：中子剂量快速确定方法，包括：1)在CPU平台上形成中子辐照的体模数据，进行体模数据的离散化，形成计算体元，并确定三维空间上每个计算体元的体元标号；所述计算体元包括大小成正倍数的关系的辐照区体元和非辐照区体元，辐照区体元小于非辐照区体元，计算体元均为立方体；2)选定用于确定中子束剂量的GPU平台，通过CPU平台向GPU平台传送计算数据；3)在GPU上利用蒙特卡罗程序并行计算三维辐照剂量矩阵D的每个单元的剂量，GPU平台向CPU平台返回辐照剂量矩阵数据D。虽然在GPU平台上并行计算的方式在其它领域有类似的应用，但是在中子辐射治疗中尚未见报道。在现有技术中，通常是仅仅在CPU平台上采用蒙特卡罗程序进行中子辐照剂量的计算，该方式由于计算的数据量庞大，导致计算速率较低。本专利技术将现有的计算过程搬入GPU中进行，采用GPU平台对蒙特卡罗程序并行运算，由于GPU本身的特性，可以使计算速率得到极大地提高。由于在采用现有技术所述的蒙特卡罗程序进行计算时，因为计算规则的限制，导致现有技术中的蒙特卡罗程序并不能完全适用于GPU与CPU中并行计算，因而，本专利技术优化了辐照的几何位置对应于计算体元中体元点的算法，通过该算法的优化设置，进而使蒙特卡罗程序能在GPU中运行时极大提高计算速率，并且有效保证计算的精度。本专利技术中采用的蒙特卡罗程序中除辐照的几何位置与计算体元中体元点的对应计算方式以外，其余均为现有技术，因此下面着重阐述辐照的几何位置与计算体元中体元点的对应计算方式，其他计算过程不再赘述。蒙特卡罗程序中，该辐照的几何位置对应于计算体元中体元点的算法如下：假设粒子所在空间点坐标为P(X,Y,Z)，计算I＝X/dX，J＝Y/dY，K＝Z/dZ，然后再利用IX、IY和IZ所确定的对应关系得到点P在三维辐照剂量矩阵D中对应的体元标号。其中，I、J、K分别为P点在数组IX、IY和IZ中的体元标号，dX为计算体元在X方向上的大小，dY为计算体元在Y方向上的大小，dZ为计算体元在Z方向上的大小；IX是在X方向上按dX大小进行离散化所形成的体元与计算体元在X方向的对应关系；IY是在Y方向上按dY大小进行离散化所形成的体元与计算体元在Y方向的对应关系；IZ是在Z方向上按dZ大小进行离散化所形成的体元与计算体元在Z方向的对应关系。进一步，本专利技术的具体实现方法如下：1)在CPU平台上形成中子辐照的体模数据，所述体模数据包括人体CT数据、体表轮廓、肿瘤区域和危及器官区域；然后进行体模数据的离散化，形成计算体元，并获得三维空间上的人体CT值矩阵H[NX*NY*NZ]，器官标志矩阵M[NX*NY*NZ]、X方向体元的长度PX[NX]、Y方向体元的长度PY[NY]、Z方向体元的长度PZ[NZ]、IX[2*NX]、IY[2*NY]和IZ[2*NZ]。根据人体CT值矩阵的值可以确定每个计算体元对应的人体CT值。根据器官标志矩阵的值可以确定每个计算体元所属的器官，进而判断该计算体元是否是肿瘤、或者是正常组织、危及器官等。为了提高计算速率，同时保证计算的精度，本专利技术中使用了不同大小的计算体元，其中肿瘤区域的计算体元较小，而非肿瘤区域的计算体元较大。即，所述计算体元包括大小成正倍数的关系的辐照区体元和非辐照区体元，辐照区体元小于非辐照区体元，计算体元均为立方体；所述辐照区体元中：dY＝dY＝dZ＝0.25cm，所述非辐照区体元的大小为非辐照区体元大小的两倍。在本专利技术中，非肿瘤区域的计算体元大小是肿瘤区域的计算体元大小的8倍，即在X、Y、Z方向上均为2倍。当采用非均匀体元时，不能简单利用体元空间点的坐标(X，Y，Z)除以体元大小(dX，dY，dZ)来获得体元的体元编号。为了加快计算速度，本专利技术中，预先按照小体元的大小来计算空间点的坐标和实际使用的体元的编号的对应关系，并保存在数组IX、IY、IZ中，并设置数组IZ、IY、IZ的长度分别为2*NX、2*NY和2*NZ。根据人体CT值矩阵H[NX*NY*NZ]确定每个计算体元的物质标号和元素组成，得到体元元素矩阵VOX[NX*NY*NZ]，其中，NX、NY、NZ分别为X、Y、Z方向上的体元数目；在CPU平台上设置中子束辐照参数，该中子束辐照参数包括：入射角度、射野形状、辐照时间、中子源强度；在CPU平台上根据体元所在的区域设置每个器官的硼含量B[NORG]，其中NORG为器官数目。2)选定用于确定中子束剂量的GPU平台，通过CPU平台向GPU平台传送计算数据；计算数据包括NX、NY、NZ、dX、dY、dZ、M、PX、PY、PZ、IX、IY、IZ、VOX、NORG、B、D、入射角度、射野形状、辐照时间、中子源强度。3)在GPU平台上利用蒙特卡罗程序通过多数据流SIMD的处理方式计算三维辐照剂量矩阵D的每个计算体元的剂量分布，GPU平台向CPU平台返回辐照剂量矩阵数据；蒙特卡罗程序中，辐照的几何位置对应于计算体元中体元点的算法如下：假设粒子所在空间点坐标为P(X,Y,Z)，计算I＝X/dX，J＝Y/dY，K＝Z/dZ，I、J、K分别为P点在数组IX、IY和IZ中的标号，然后利用IX、IY和IZ所确定的对应关系得到点P在三维辐照剂量矩阵D中对应的体元标号。CPU，英文CentralProcessingUnit的缩写，意指中央处理机。本文档来自技高网...

【技术保护点】
中子剂量快速确定方法，其特征在于，包括：1)在CPU平台上形成中子辐照的体模数据，进行体模数据的离散化，形成计算体元，并确定三维空间上每个计算体元的体元标号；所述计算体元包括大小成正倍数的关系的辐照区体元和非辐照区体元，辐照区体元小于非辐照区体元，计算体元均为立方体；2)选定用于确定中子束剂量的GPU平台，通过CPU平台向GPU平台传送计算数据；3)在GPU平台上利用蒙特卡罗程序计算三维辐照剂量矩阵D的每个计算体元的剂量，GPU平台向CPU平台返回辐照剂量矩阵数据D。

【技术特征摘要】
1.中子剂量快速确定方法，其特征在于，包括：1)在CPU平台上形成中子辐照的体模数据，进行体模数据的离散化，形成计算体元，并确定三维空间上每个计算体元的体元标号；所述计算体元包括大小成正倍数的关系的辐照区体元和非辐照区体元，辐照区体元小于非辐照区体元，计算体元均为立方体；2)选定用于确定中子束剂量的GPU平台，通过CPU平台向GPU平台传送计算数据；3)在GPU平台上利用蒙特卡罗程序计算三维辐照剂量矩阵D的每个计算体元的剂量，GPU平台向CPU平台返回辐照剂量矩阵数据D。2.根据权利要求1所述的中子剂量快速确定方法，其特征在于，所述蒙特卡罗程序中，辐照的几何位置对应于计算体元中体元点的算法如下：假设粒子所在空间点坐标为P(X,Y,Z)，计算I＝X/dX，J＝Y/dY，K＝Z/dZ，然后再利用IX、IY和IZ所确定的对应关系得到点P在三维辐照剂量矩阵D中对应的体元标号；其中，I、J、K分别为P点在数组IX、IY和IZ中的体元标号，dX为计算体元在X方向上的大小，dY为计算体元在Y方向上的大小，dZ为计算体元在Z方向上的大小；IX、IY和IZ所确定的对应关系如下：IX是在X方向上按dX大小进行离散化所形成的体元与计算体元在X方向的对应关系；IY是在Y方向上按dY大小进行离散化所形成的体元与计算体元在Y方向的对应关系；IZ是在Z方向上按dZ大小进行离散化所形成的体元与计算体元在Z方向的对应关系。3.根据权利要求2所述的中子剂量快速确定方法，其特征在于，步骤1)的过程为：在CPU平台上...

【专利技术属性】
技术研发人员：王慧芳，周烨，
申请(专利权)人：四川瑶天纳米科技有限责任公司，
类型：发明
国别省市：四川,51