(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211090683.5 (22)申请日 2022.09.07 (71)申请人 西北核技 术研究所 地址 710024 陕西省西安市灞桥区平峪路 28号 (72)发明人 商鹏　左应红　朱金辉　牛胜利　 伏琰军　陈统　 (74)专利代理 机构 西安智邦专利商标代理有限 公司 6121 1 专利代理师 徐秦中 (51)Int.Cl. G06F 30/20(2020.01) G06F 111/10(2020.01) G06F 119/02(2020.01) (54)发明名称 一种基于插值的瞬发核辐射环境参数的快 速获取方法 (57)摘要 本发明为解决现有瞬发核辐射环境参数的 获取方法存在无法兼顾计算结果精度与计算求 解时间的问题， 而提供了一种基于插值的瞬发核 辐射环境参数快速获取方法。 本发 明基于蒙特卡 罗法预计算二维数据表进行插值， 从距离和源高 两个方向分别插值得到了目标点处的瞬发核辐 射环境参数， 取二者平均作为最终待测目标处的 瞬发核辐射环 境参数结果， 保持较高计算精度的 同时， 有效减少计算时间， 从而快速给出任意距 离和任意源高条件下地面上方一定高度处的瞬 发核辐射环境 参数。 权利要求书2页 说明书6页 附图3页 CN 115481532 A 2022.12.16 CN 115481532 A 1.一种基于插值的瞬发核辐射环境 参数快速获取 方法， 其特 征在于， 包括以下步骤： 步骤1)、 建立瞬发核辐射输运模型； 步骤2)、 利用蒙特卡罗法， 在二维矩形域上通过调整辅助探测器的距离d和辐射源高度 h， 计算给出n×m个结点(di,hj)(i＝0,1, …,n‑1； j＝0,1, …,m‑1)的瞬发核辐射环境参数值 P(hi,dj)， 并生成二维数据表， 其中， m≥8,n≥8,且均为 正整数； 步骤3)、 确定待测目标场景参数： 目标距离插值点dtar(d0≤dtar≤dn‑1)与目标源高插值 点htar(h0≤htar≤hm‑1)； 步骤4)、 查找目标距离插值点dtar所在二维数据表di数列中的左边界dleft和右边界 dright， 并根据已知的左、 右边界在二维矩形域中分别获取其对应的一维源高数据 结点序列 值P(dleft,hj)和P(dright,hj)； 步骤5)、 分别以P(dleft,hj)和P(dright,hj)作为插值结点序列， 插值计算左、 右边界序列 上目标源高插值 点htar处的核辐射环境 参数近似值P(dleft,htar)和P(dright,htar)； 步骤6)、 以左、 右边界插值结果P(dleft,htar)和P(dright,htar)为数据点， 插值计算目标距 离插值点dtar处的瞬发核辐射环境 参数值P1(dtar,htar)； 步骤7)、 查找目标源高插值点htar在二维数据表hj数列中的下边界hdown和上边界hup， 并 根据已知的上、 下边界在二维矩形域中分别获取其对应的一维距离数据结点序列P(di, hdown)和P(di,hup)； 步骤8)、 分别以P(di,hdown)和P(di,hup)作为插值结点序列， 插值计算上、 下边界序列上 目标距离插值 点dtar处的核辐射环境 参数近似值P(dtar,hdown)和P(dtar,hup)； 步骤9)、 以上、 下边界插值结果P(dtar,hdown)和P(dtar,hup)为数据点， 插值计算目标源高 插值点htar处的瞬发核辐射环境 参数值P2(dtar,htar)； 步骤10)、 取步骤6)的插值结果P1(dtar,htar)与步骤9)的插值结果P2(dtar,htar)的平均值 作为目标距离插值点dtar和目标源高插值点htar的瞬发核辐射环 境参数值， 即待测目标场景 的核辐射环境 参数值为P(dtar,htar)＝(P1+P2)/2。 2.根据权利要求1所述的基于插值的瞬发核辐射环境 参数快速获取 方法， 其特 征在于： 步骤4)和步骤7)中， 均采用二 分法查找目标距 离插值点dtar所在二维数据表di数列中的 左边界dleft和右边界dright， 以及目标源高插值点htar在二维数据表hj数列中的下边界hdown和 上边界hup。 3.根据权利要求2所述的基于插值的瞬发核辐射环境 参数快速获取 方法， 其特 征在于： 步骤6)和步骤9)中， 均采用线性插值方法进行插值计算。 4.根据权利要求3所述的基于插值的瞬发核辐射环境 参数快速获取 方法， 其特 征在于： 步骤5)和步骤8)中， 均采用拉格朗日插值法进行插值计算。 5.根据权利要求1 ‑4任一所述的基于插值的瞬发核辐射环境参数快速获取方法， 其特 征在于： 步骤5)、 步骤6)、 步骤8)和步骤9)中， 所述插值计算过程 为： ①将作为插值数据点的核辐射环境 参数值取以10为底的对数； ②进行插值计算； ③将得到的插值结果作为指数进行一次以10为底的指数运算转换为对应插值点的瞬 发核辐射环境 参数值。权　利　要　求　书 1/2 页 2 CN 115481532 A 26.根据权利要求 4所述的基于插值的瞬发核辐射环境 参数快速获取 方法， 其特 征在于： 步骤5)和步骤8)中采用拉格朗日插值法插值计算时， 选取的插值结点至少为4个， 其确 定规则为： 当插值点位于结点区间中且其前后结点数目均大于4个时， 选取的插值结点为8个； 当 插值点靠近结点区间任一端且其一侧结点数目小于4个时， 选取的插值结点为4～7个； 当插 值结点位于结点区间外时， 选取的插值结点 为结点区间中靠 近插值点的4个结点。 7.根据权利要求5所述的基于插值的瞬发核辐射环境 参数快速获取 方法， 其特 征在于： 步骤2)中， 在二维矩形域上， n ×m个结点的选取是非均匀的， 在辐射源附近结点的选取 密集， 距离辐射源较 远处结点的选取稀疏。权　利　要　求　书 2/2 页 3 CN 115481532 A 3