CN113126142B - 高能粒子探测器性能评估方法及系统 - Google Patents
高能粒子探测器性能评估方法及系统 Download PDFInfo
- Publication number
- CN113126142B CN113126142B CN202110411873.1A CN202110411873A CN113126142B CN 113126142 B CN113126142 B CN 113126142B CN 202110411873 A CN202110411873 A CN 202110411873A CN 113126142 B CN113126142 B CN 113126142B
- Authority
- CN
- China
- Prior art keywords
- energy
- flux
- cut
- particle
- particle detector
- 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.)
- Expired - Fee Related
Links
- 239000002245 particle Substances 0.000 title claims abstract description 77
- 238000011156 evaluation Methods 0.000 title claims abstract description 7
- 230000004907 flux Effects 0.000 claims abstract description 51
- 230000008859 change Effects 0.000 claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 230000005358 geomagnetic field Effects 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 238000007619 statistical method Methods 0.000 claims abstract description 5
- 238000012937 correction Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 description 10
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 230000005855 radiation Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种高能粒子探测器性能评估方法及系统。该方法包括：通过AP8计算高能质子的全向通量；使用数值模拟方法得到地磁场长周期变化导致的高能质子垂直截止刚度的长周期变化趋势；拟合得到截止刚度的长周期相对变化与时间的函数曲线关系；通过上述函数曲线关系，对AP8计算得到的粒子通量进行修正；通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征；根据高能质子对应截止刚度的投掷角分布特征，结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量；获取电磁卫星高能粒子探测器探测到的粒子通量；使用上述获得的卫星观测和数值模型通量，对两者进行统计分析，得到包括其整体趋势、线性相关性、绝对偏差、标准差以及均方差大小的用于评价探测器性能的参量。本发明提供的高能粒子探测器性能评估方法及系统可以实现对粒子探测器在轨性能的有效评估。
Description
技术领域
本发明涉及高能粒子探测技术领域，特别是涉及一种高能粒子探测器性能评估方法及系统。
背景技术
近地空间高能粒子辐射作为地球物理场卫星探测的主要参量，其主要通过卫星搭载的高能粒子探测器实现对高能粒子的能量、通量、投掷角等开展测量。卫星发射之前，将对高能粒子探测器的性能进行地面标定试验，主要通过地面束线的方法对其能量线性、探测效率、几何因子等进行有效的验证。卫星入轨之后，需要对探测器性能进行定期评估，并进行必要的参数调整，以获得最佳的探测结果。但由于一般探测器不会搭载星上标定源，因此无法获得探测器探测位置的物理量真值。目前主要方法是使用数值模型与探测器探测结果比对、星地交叉比对以及不同卫星的不同探测器交叉比对等三种方法。
其中星地交叉通常使用地面中子台站的观测数据，通过反演得到高能粒子通量，此方法对地面台站的布局以及数量要求较高。由于地面台站和卫星处于不同的高度，同时其反演高能粒子通量的方法和过程影响最终结果，可能导致比测校验结果不理想。
对于不同卫星的高能粒子探测器的交叉对比。由于不同卫星的轨道高度、轨道倾角、探测器探测原理、探测能段范围、分辨率以及探测器指向等原因，要找到完全符合两者交叉的探测数据，相对困难。
鉴于上述两种方法存在的难题，目前主要方法为使用数值模型和卫星观测结果进行比对分析。数值模型由于是统计模型，相对星地交叉比测校验以及不同卫星的比测校验具有明显的优势。
但目前使用范围最广泛AP8模型，是美国于上世纪60年代开发的模型。众所周知，近地空间高能粒子通量与地球磁场结构有着重要的关系。由于地磁场具有长周期变化特征，尤其是近年来地磁变化明显加快(IGRF13模型提前发布)，使得地球磁场对高能粒子的屏蔽作用将发生明显变化，导致太阳/银河宇宙线进入磁层导致的背景高能粒子环境改变。因此需要对磁场变化导致的高能粒子截止刚度进行重新衡量，并对AP8模型进行一定的修正，这也是美国开发新模型的原因。
目前见于国外网站和期刊组织的AP9模型是美国最新开发的质子辐射环境模型，但由于美国技术封锁等原因，我们国内一直无法获取其程序。因此需要通过方法解决此问题。
地球内辐射带<100MeV的质子主要来源于太阳质子，>100MeV的质子主要来自于宇宙线反照中子产生，这种机制称为宇宙线反照中子源(CRAND)，包括太阳质子和银河宇宙线都可以称为宇宙线。因此宇宙线进入磁层对于内辐射带的动态变化具有重要作用。由于地球磁场的屏蔽作用，宇宙线粒子不能随意的进入地球磁层空间，而是具有能量和方向依赖性。
地磁截止刚度是定量衡量地球磁场对高能粒子屏蔽效应的参数，描述了高能粒子穿越地磁屏蔽效应到达指定观测点的带电粒子刚度阈值。高能粒子刚度作为描述高能粒子能量大小的物理量，其定义如下：R＝mvc/Ze，其中R表示粒子刚度，m为高能粒子质量，v为粒子速度，c为光速，Z为电荷数，e为元电荷电量。刚度在数值上等于单位单核的能量，单位为伏特。截止刚度表征到达给定观测点的最小刚度值。
发明内容
本发明要解决的技术问题是提供一种高能粒子探测器性能评估方法及系统，可以实现对粒子探测器在轨性能的有效评估。
为解决上述技术问题，本发明提供了一种高能粒子探测器性能评估方法，所述方法包括：通过AP8计算高能质子的全向通量；使用数值模拟方法得到地磁场长周期变化导致的高能质子垂直截止刚度的长周期变化趋势；拟合得到截止刚度的长周期相对变化与时间的函数曲线关系；通过上述函数曲线关系，对AP8计算得到的粒子通量进行修正；通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征；根据高能质子对应截止刚度的投掷角分布特征，结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量；获取电磁卫星高能粒子探测器探测到的粒子通量；使用上述获得的卫星观测和数值模型通量，对两者进行统计分析，得到包括其整体趋势、线性相关性、绝对偏差、标准差以及均方差大小的用于评价探测器性能的参量。
在一些实施方式中，由高能粒子垂直方向截止刚度的长周期变化得到其时间关系函数为：
其中，R为高能粒子刚度，
为位置矢量，t为时间。
在一些实施方式中，通过上述函数曲线关系，对AP8计算得到的粒子通量进行修正，包括：
其中，F为使用AP8模型给定的高能质子的全向通量，F′为修正后的全向通量，g为函数曲线关系。
在一些实施方式中，通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征，包括：通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征函数f(α)，其中，f(α)为对应投掷角α的高能质子通量。
在一些实施方式中，根据高能质子对应截止刚度的投掷角分布特征，结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量，包括：
其中，α为投掷角，f(α)为对应投掷角α的高能质子通量。
此外，本发明还提供了一种高能粒子探测器性能评估系统，所述系统包括：一个或多个处理器；存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现根据前文所述的高能粒子探测器性能评估方法。
采用这样的设计后，本发明至少具有以下优点：
可以实现对粒子探测器在轨性能的有效评估；有效地修正由于地磁场长周期变化导致的AP8模型在当前的不准确造成的粒子探测器和数值模型的不吻合现象；明显提高观测和数值模型的线性相关性；使观测和数值模型的均方差明显变小。
附图说明
上述仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，以下结合附图与具体实施方式对本发明作进一步的详细说明。
图1是高能粒子探测器性能评估方法的流程示意图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
参见图1，本申请方法的具体步骤如下：
通过AP8计算高能质子的全向通量。
使用数值模拟方法得到地磁场长周期变化导致的高能质子垂直截止刚度的长周期变化趋势。
拟合得到截止刚度的长周期相对变化与时间的函数曲线关系。
通过上述函数曲线关系，对AP8计算得到的粒子通量进行修正。
通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征。
根据高能质子对应截止刚度的投掷角分布特征，结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量。
获取电磁卫星高能粒子探测器探测到的粒子通量。
使用上述获得的卫星观测和数值模型通量，对两者进行统计分析，得到包括其整体趋势、线性相关性、绝对偏差、标准差以及均方差大小的用于评价探测器性能的参量。
高能粒子垂直方向截止刚度的长周期变化。
由于截止刚度对纬度的依赖性明显大于经度依赖性，但地磁场的变化的经度依赖性明显，因此为了对全球进行高精度的空间变化特征分析，我们使用纬度方向间隔1度，经度方向5(10)度对全球进行网格的划分。计算得到全球垂直方向截止刚度的时空演化特征。考虑长周期地磁场的变化，我们仅使用国际地磁参考场(IGRF)内源场作为背景场进行数值模拟计算。为了使研究结果具有最大的适用性，选取地磁平静期间进行研究。
由高能粒子垂直方向截止刚度的长周期变化得到其时间关系函数为
其中R为高能粒子刚度，
为位置矢量，t为时间。
利用垂直方向截止刚度长时间变化的函数曲线关系，对AP8模型进行修正。
假设使用AP8模型给定的高能质子的全向通量为F,进行时间修正之后关系为
高能粒子截止刚度的投掷角依赖性计算方法。
通常情况下，粒子投掷角分布满足关于90°对称的高斯分布形态，但初步研究表明，投掷角分布关于90°度存在一定偏差，因此为了使结果真实可靠，通过计算电磁卫星高能粒子探测器位置的高能质子截止刚度对应的投掷角分布特征，可以得到f(α)的函数形态。其中α为粒子投掷角，f(α)为对应投掷角α的高能质子通量。
结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量。
AP8模型可以计算不同能量不同位置的高能质子通量，但其无法给出不同投掷角对应的高能质子通量。全向通量为投掷角的积分函数，满足如下关系
通过上述关系式子，可以得到对应能段和投掷角的高能质子通量。
通过选取地磁平静期间的高能粒子探测器探测数据，结合修正之后的AP8模型的结果。对两者进行统计分析，得到包括其整体趋势、线性相关性、绝对偏差、标准差以及均方差大小等用于评价探测器性能的参量。
本发明的目的在于，为解决现有使用AP8模型进行张衡一号卫星高能粒子探测器评估导致的偏差较大的现状。通过使用上述引进的技术创新方法：可以实现对粒子探测器在轨性能的有效评估；有效地修正由于地磁场长周期变化导致的AP8模型在当前的不准确造成的粒子探测器和数值模型的不吻合现象；明显提高观测和数值模型的线性相关性；使观测和数值模型的均方差明显变小。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。
Claims (3)
1.一种高能粒子探测器性能评估方法，其特征在于，包括：
通过AP8计算高能质子的全向通量；
使用数值模拟方法得到地磁场长周期变化导致的高能质子垂直截止刚度的长周期变化趋势；
拟合得到截止刚度的长周期相对变化与时间的函数曲线关系；
通过上述函数曲线关系，对AP8计算得到的粒子通量进行修正；
通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征；
根据高能质子对应截止刚度的投掷角分布特征，结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量；
获取电磁卫星高能粒子探测器探测到的粒子通量；
使用上述获得的卫星观测和数值模型通量，对两者进行统计分析，得到包括其整体趋势、线性相关性、绝对偏差、标准差以及均方差大小的用于评价探测器性能的参量；
由高能粒子垂直方向截止刚度的长周期变化得到其时间关系函数为：
其中，R为高能粒子刚度，
为位置矢量，t为时间；
通过上述函数曲线关系，对AP8计算得到的粒子通量进行修正，包括：
其中，F为使用AP8模型给定的高能质子的全向通量，F′为修正后的全向通量，g为函数曲线关系；
通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征，包括：
通过数值计算得到卫星观测点的高能质子对应截止刚度的投掷角分布特征函数f(α)，其中，f(α)为对应投掷角α的高能质子通量。
2.根据权利要求1所述的高能粒子探测器性能评估方法，其特征在于，根据高能质子对应截止刚度的投掷角分布特征，结合AP8模型得到的修正之后的全向通量，得到不同投掷角的高能粒子通量，包括：
其中，α为投掷角，f(α)为对应投掷角α的高能质子通量。
3.一种高能粒子探测器性能评估方法系统，其特征在于，包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序，
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现根据权利要求1至2任意一项所述的高能粒子探测器性能评估方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110411873.1A CN113126142B (zh) | 2021-04-16 | 2021-04-16 | 高能粒子探测器性能评估方法及系统 |
| AU2021103770A AU2021103770A4 (en) | 2021-04-16 | 2021-06-30 | Method and system for evaluating performance of energetic particle detector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110411873.1A CN113126142B (zh) | 2021-04-16 | 2021-04-16 | 高能粒子探测器性能评估方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113126142A CN113126142A (zh) | 2021-07-16 |
| CN113126142B true CN113126142B (zh) | 2022-04-01 |
Family
ID=76776808
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110411873.1A Expired - Fee Related CN113126142B (zh) | 2021-04-16 | 2021-04-16 | 高能粒子探测器性能评估方法及系统 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN113126142B (zh) |
| AU (1) | AU2021103770A4 (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459372A (zh) * | 2014-11-10 | 2015-03-25 | 中国科学院新疆理化技术研究所 | 基于p-i-n结构的位移损伤剂量探测方法 |
| GB201708270D0 (en) * | 2017-05-23 | 2017-07-05 | Cheyney Design & Dev Ltd | Improvements in or relating to detectors |
| CN108426856A (zh) * | 2018-03-26 | 2018-08-21 | 武汉大学 | 评估激光雷达卫星测量大气co2浓度性能的综合分析系统 |
| CN111722267A (zh) * | 2020-06-23 | 2020-09-29 | 中国科学院国家空间科学中心 | 一种磁层能量粒子事件爆发源的探测追踪方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2164361C1 (ru) * | 1999-10-18 | 2001-03-20 | Кумахов Мурадин Абубекирович | Линза для управления излучением в виде потока нейтральных или заряженных частиц, способ изготовления таких линз и содержащее такие линзы аналитическое устройство, устройство для лучевой терапии и устройства для контактной и проекционной литографии |
| US7138629B2 (en) * | 2003-04-22 | 2006-11-21 | Ebara Corporation | Testing apparatus using charged particles and device manufacturing method using the testing apparatus |
| EP2195684A2 (en) * | 2007-09-21 | 2010-06-16 | Purdue Research Foundation | Detection of solar events |
| CN104297585A (zh) * | 2014-09-22 | 2015-01-21 | 中国空间技术研究院 | 一种用高能质子进行空间位移损伤效应评估试验的方法 |
| CN105676255B (zh) * | 2014-11-21 | 2018-12-04 | 中国科学院空间科学与应用研究中心 | 一种分布式空间高能质子的多方向探测装置 |
| CN106291649B (zh) * | 2015-05-25 | 2018-12-04 | 中国科学院空间科学与应用研究中心 | 一种利用地磁刚度进行空间离子探测器定标的方法 |
| CN105974905B (zh) * | 2016-05-10 | 2018-08-17 | 中国民航大学 | 面向航空数据总线单粒子翻转故障的仿真测试系统及方法 |
| CN112363201B (zh) * | 2020-11-10 | 2023-03-14 | 国家卫星气象中心(国家空间天气监测预警中心) | 一种中轨道卫星高能电子能谱数据的获取方法 |
-
2021
- 2021-04-16 CN CN202110411873.1A patent/CN113126142B/zh not_active Expired - Fee Related
- 2021-06-30 AU AU2021103770A patent/AU2021103770A4/en not_active Ceased
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459372A (zh) * | 2014-11-10 | 2015-03-25 | 中国科学院新疆理化技术研究所 | 基于p-i-n结构的位移损伤剂量探测方法 |
| GB201708270D0 (en) * | 2017-05-23 | 2017-07-05 | Cheyney Design & Dev Ltd | Improvements in or relating to detectors |
| CN108426856A (zh) * | 2018-03-26 | 2018-08-21 | 武汉大学 | 评估激光雷达卫星测量大气co2浓度性能的综合分析系统 |
| CN111722267A (zh) * | 2020-06-23 | 2020-09-29 | 中国科学院国家空间科学中心 | 一种磁层能量粒子事件爆发源的探测追踪方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| AU2021103770A4 (en) | 2021-08-19 |
| CN113126142A (zh) | 2021-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lipari et al. | Diffuse Galactic gamma-ray flux at very high energy | |
| Litvak et al. | LEND neutron data processing for the mapping of the Moon | |
| Bruel et al. | Gamma rays, electrons and positrons up to 3 TeV with the Fermi Gamma-ray Space Telescope | |
| Panda et al. | Ionospheric TEC variations at low latitude Indian region | |
| Mitra et al. | Reconstructing air shower parameters with LOFAR using event specific GDAS atmosphere | |
| Bruno et al. | Geomagnetically trapped, albedo and solar energetic particles: Trajectory analysis and flux reconstruction with PAMELA | |
| CN114384564B (zh) | 一种基于多源数据驱动的电离层层析成像方法 | |
| CN113126142B (zh) | 高能粒子探测器性能评估方法及系统 | |
| Vashenyuk et al. | Relativistic solar protons in the event of January 20, 2005: Model studies | |
| Anghel et al. | Kalman filter-based algorithms for monitoring the ionosphere and plasmasphere with GPS in near-real time | |
| Pak et al. | A numerical method to analyze geometric factors of a space particle detector relative to omnidirectional proton and electron fluxes | |
| Shprits et al. | Data assimilation in the near-earth electron radiation environment | |
| Badavi et al. | Low Earth orbit assessment of proton anisotropy using AP8 and AP9 trapped proton models | |
| Magkos et al. | Deflections of ultra-high energy cosmic rays by the Milky Way magnetic field: how well can they be corrected? | |
| Pu et al. | Drift shell tracing and secular variation of inner zone high energy proton environment in the SAA | |
| Sarkar et al. | Monte Carlo simulation of CRAND protons trapped at low Earth orbits | |
| Christe et al. | The RHESSI microflare height distribution | |
| Alpatov et al. | Analyzing existing applied models of the ionosphere to calculate radio wave propagation and a possibility of their use for radar-tracking systems. II. Domestic models | |
| Astapov et al. | Investigation of local deformations of muon flux angular distribution during CME with GSE-mapping technique | |
| CN117008154B (zh) | 一种基于松弛因子逆时衰减函数的快速电离层层析方法 | |
| Gebauer et al. | Measurement of anisotropies in cosmic ray arrival directions with the Alpha Magnetic Spectrometer on the ISS | |
| US20230142341A1 (en) | Method and system for monitoring the precipitation of particles in the magnetosphere | |
| Zanimonskiy et al. | Modeling of TEC variations based on signals from near zenith GNSS satellite observed by dense regional network | |
| Gargano et al. | The DAMPE experiment: 2 year in orbit | |
| Mikhailov et al. | Observations of Trapped Electrons and Positrons with E> 50 MeV in the Inner Radiation Belt by the PAMELA Magnetic Spectrometer |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220401 |