CN-121583366-B - Radiation environment monitoring data dosage calculation method and system

CN121583366BCN 121583366 BCN121583366 BCN 121583366BCN-121583366-B

Abstract

The application relates to the technical field of radiation environment monitoring and discloses a radiation environment monitoring data dosage calculation method and system. The method comprises the steps of obtaining radiation source data, calculating initial radiation field distribution based on the data, obtaining preliminary dose deposition values of all spatial points of a monitoring area, analyzing interaction among particles, determining energy redistribution coefficients of all the spatial points, utilizing the coefficients to adjust the preliminary dose deposition values point by point to obtain corrected three-dimensional dose concentration field data, calculating gradient vectors among the spatial points to generate an optimized gradient distribution map, extracting boundary point sets with gradient change exceeding a threshold value, constructing an equivalent surface grid to form a three-dimensional equivalent surface structure, obtaining iterative update coefficients of time sequence radiation source data to obtain a concentration field time sequence evolution sequence, identifying a high concentration area based on the sequence, fusing the high concentration area and the three-dimensional equivalent surface structure, and generating a visualized distribution model reflecting dose concentration.

Inventors

SHI DONGFENG
WANG JIAYU
Mao xiyin
LIANG SHIMIN
XU ZHIDAN
ZhuGe Wenting
CHEN QIE
CHEN XUE
ZHANG QINYI
ZHONG DAN

Assignees

杭州旭辐检测技术有限公司

Dates

Publication Date: 20260508
Application Date: 20260127

Claims (7)

1. A radiation environment monitoring data based dose calculation method, the method comprising: Step S1, acquiring radiation source data, calculating initial radiation field distribution based on the radiation source data, obtaining preliminary dose deposition values of all spatial points in a monitored area, analyzing inter-particle interactions reflected by the preliminary dose deposition values, and determining energy redistribution coefficients corresponding to all the spatial points; The method comprises the steps of S1, determining energy redistribution coefficients corresponding to all space points, wherein the step comprises the steps of carrying out numerical analysis on the preliminary dose deposition values by adopting a finite difference method, simulating collision events of radiation particles in space, recording the occurrence positions and the intensities of the collision events, calculating energy transfer amounts among the particles, determining initial energy redistribution coefficients of all the space points based on the energy transfer amounts and total energy of the corresponding space points, constructing an energy redistribution model according to the initial energy redistribution coefficients and combining the space distribution characteristics and the intensity distribution rules of the collision events, adjusting the energy distribution relation among the space points in a monitoring area according to the energy redistribution model, carrying out iterative optimization on the initial energy coefficients by using the adjusted space point energy distribution relation as a constraint condition by a gradient descent method, and outputting final energy redistribution coefficients; The method comprises the steps of setting a local acting domain for each space point, calculating energy transfer probability among each point in the local acting domain based on collision event distribution characteristics in an energy redistribution model, constructing a Lagrangian dynamics system by taking a preliminary dose deposition value of each space point as a generalized coordinate and taking an intensity distribution rule in the energy redistribution model as a potential energy field, obtaining a dose distribution evolution track by solving a Euler-Lagrangian equation with constraint, extracting a generalized coordinate convergence value when the system reaches a steady state from the evolution track, taking the generalized coordinate convergence value as a dose deposition value after adjustment of each space point, and recalculating energy gradient distribution among adjacent space points to obtain an adjusted energy distribution relation; Step S2, utilizing the energy redistribution coefficient to carry out point-by-point adjustment on the preliminary dose deposition value to obtain corrected three-dimensional dose concentration field data, and calculating gradient vectors among all spatial points based on the corrected three-dimensional dose concentration field data to obtain an optimized gradient distribution diagram; The method comprises the steps of S2, obtaining corrected three-dimensional dose concentration field data, namely, carrying out point-by-point adjustment on the preliminary dose deposition value according to the energy redistribution coefficient, calculating the dose value adjusted by each space point, marking the corresponding space point as a potential high-concentration point if the adjusted dose value exceeds a preset threshold value, recording three-dimensional coordinate positions and corresponding dose value information aiming at the potential high-concentration point, and generating initial three-dimensional dose concentration field data by integrating the adjusted dose values of all the space points in a monitoring area; S3, extracting a boundary point set with gradient change exceeding a preset threshold value from the optimized gradient distribution diagram, constructing an isosurface grid based on the boundary point set to form a three-dimensional isosurface structure, acquiring time sequence radiation source data, and iteratively updating the energy redistribution coefficient to obtain a time sequence evolution sequence of a concentration field; And S4, identifying and marking a high-concentration region based on the time sequence evolution sequence of the concentration field, fusing the high-concentration region and the three-dimensional isosurface structure, and generating a visualized distribution model reflecting the dose concentration according to a fusion result.
2. The radiation environment monitoring data based dose calculation method according to claim 1, wherein obtaining preliminary dose deposition values at each spatial point in the monitored area in step S1 comprises: The method comprises the steps of extracting radiation source data from a preset database, wherein the radiation source data comprise positions, energy and types of radiation sources, processing the radiation source data by adopting a Monte Carlo simulation method, simulating a propagation path and an energy deposition process of particles in a monitoring area space, calculating a dose deposition value on each space point according to a simulation result, generating distribution data of preliminary dose deposition values, carrying out gridding processing on the distribution data of the preliminary dose deposition values, constructing a mapping relation between the space points and the dose values, and outputting the gridded preliminary dose deposition values.
3. The radiation environment monitoring data based dose calculation method according to claim 1, wherein the obtaining of the optimized gradient profile in step S2 comprises: calculating concentration differences between adjacent space points according to the corrected three-dimensional dose concentration field data, determining gradient vectors between the adjacent space points by taking the concentration differences as components and combining the intervals between the adjacent space points, calculating the direction offset angle of each gradient vector and the regional average gradient vector, judging whether the direction offset angle exceeds a preset range, if so, performing smoothing treatment on the gradient vectors, constructing an initial gradient distribution map according to the gradient vectors after the smoothing treatment, and verifying the spatial continuity and the direction consistency of the gradient vectors for the initial gradient distribution map to generate an optimized gradient distribution map.
4. The radiation environment monitoring data dose calculation method according to claim 1, wherein the forming of the three-dimensional isosurface structure in step S3 comprises: identifying and extracting space points with gradient amplitude exceeding a preset threshold value as a boundary point set according to the optimized gradient distribution diagram, adopting an interpolation algorithm to calculate connection relations among boundary points according to the boundary point set, constructing a triangular grid surface according to the connection relations to generate an isosurface grid, obtaining a preliminary three-dimensional isosurface structure according to the isosurface grid, verifying the topological integrity of the grid according to the preliminary three-dimensional isosurface structure, and optimizing and adjusting the isosurface grid to generate final three-dimensional isosurface structure data.
5. The radiation environment monitoring data dose calculation method according to claim 4, wherein the time-series evolution sequence of the concentration field is obtained in step S3, comprising: Obtaining time sequence radiation source data corresponding to different time stamps from a preset database, simulating the dynamic change process of a radiation field in a monitoring area along with time according to the time sequence radiation source data and the three-dimensional isosurface structure, iteratively updating the energy redistribution coefficient through change simulation, adjusting concentration field distribution according to the updated energy redistribution coefficient to generate an initial concentration field time sequence evolution sequence, verifying the time continuity of the data according to the initial concentration field time sequence evolution sequence, and performing smoothing treatment on the initial concentration field time sequence evolution sequence to generate a final concentration field time sequence evolution sequence.
6. The radiation environment monitoring data based dose calculation method according to claim 1, wherein generating a visualized distribution model reflecting dose concentration according to the fusion result in step S4 comprises: Extracting dynamic change information of a high-concentration region according to a time sequence evolution sequence of a concentration field, marking a high-concentration key region in the three-dimensional equivalent surface structure by combining the high-concentration region mark, carrying out data fusion on the high-concentration key region and the time sequence evolution sequence through a weighted average method to construct a preliminary visual distribution model, generating a three-dimensional dynamic dose distribution map according to the preliminary visual distribution model, calculating rendering coverage rate aiming at the three-dimensional dynamic dose distribution map, verifying visual integrity of data, and outputting a final visual distribution model reflecting dose concentration.
7. A radiation environment monitoring based data dose calculation system for implementing a radiation environment monitoring based data dose calculation method as claimed in any one of claims 1-6, the system comprising: The acquisition module is used for acquiring radiation source data, calculating initial radiation field distribution based on the radiation source data, obtaining preliminary dose deposition values on all spatial points in a monitoring area, analyzing particle interactions reflected by the preliminary dose deposition values, and determining energy redistribution coefficients corresponding to all the spatial points; The optimizing module is used for utilizing the energy redistribution coefficient to carry out point-by-point adjustment on the preliminary dose deposition value to obtain corrected three-dimensional dose concentration field data, and calculating gradient vectors among all space points based on the corrected three-dimensional dose concentration field data to obtain an optimized gradient distribution diagram; The iteration updating module is used for extracting a boundary point set with gradient change exceeding a preset threshold value from the optimized gradient distribution diagram, constructing an isosurface grid based on the boundary point set to form a three-dimensional isosurface structure, acquiring time sequence radiation source data, and iteratively updating the energy redistribution coefficient to obtain a time sequence evolution sequence of a concentration field; And the visualization module is used for identifying and marking a high-concentration region based on the time sequence evolution sequence of the concentration field, fusing the high-concentration region and the three-dimensional isosurface structure, and generating a visualized distribution model reflecting the dose concentration according to a fusion result.

Description

Radiation environment monitoring data dosage calculation method and system Technical Field The application relates to the technical field of radiation environment monitoring, in particular to a radiation environment monitoring data-based dosage calculation method and system. Background The radiation environment monitoring is a core technical support for guaranteeing the safe operation of nuclear facilities, the prevention and control of radiation pollution and emergency response, and is directly related to public health, ecological environment safety and social stability, and the accurate calculation and dynamic tracking of radiation dose are key technical requirements in the field. In actual scenes such as radiation leakage emergency treatment, medical radiation therapy planning and the like, the spatial distribution, intensity change and time evolution rule of radiation particles in a monitoring area are determined through dose calculation, and scientific basis is provided for risk assessment, protective measure formulation and emergency decision making. The current mainstream dose calculation method mainly comprises Monte Carlo simulation, but has the obvious defects that firstly, only an initial dose deposition value in a statistical sense can be output, the energy redistribution effect caused by particle collision and scattering is difficult to accurately capture, the boundary of a three-dimensional dose concentration field is fuzzy, the identification precision of a high-concentration risk area is insufficient, the requirement of fine monitoring cannot be met, secondly, a collaborative iteration mechanism of an energy redistribution coefficient and a spatial structure is lacking in a dynamic scene, after time sequence radiation source data is updated, a correction process and boundary extraction are disjointed, a concentration field time sequence evolution sequence is not consistent, a radiation diffusion track and intensity change cannot be truly restored, and emergency response efficiency is easily influenced by data distortion in the emergency scene. In addition, the fusion depth of the monitoring data in the prior art is insufficient, environmental data optimization models such as medium properties and topography features are not fully combined, the visualization degree of the dose distribution result is low, dynamic changes of a high-concentration area are difficult to intuitively present, risk situations cannot be quickly mastered by monitoring staff, and the technology is restricted to land in an actual scene. Therefore, it is necessary to develop a dose calculation method based on radiation environment monitoring data, which has both accuracy, dynamics and visualization characteristics, so as to solve the pain point in the prior art, and become an urgent need in the radiation environment monitoring field. Disclosure of Invention In order to solve the technical problems, the application provides a radiation environment monitoring data-based dose calculation method and system, which are used for describing particle interaction and radiation field space-time evolution rules, greatly improving the accuracy and dynamic analysis capability of radiation dose calculation and realizing visual presentation of dose distribution. In a first aspect, the present application provides a radiation environment monitoring data based dose calculation method, the method comprising: Step S1, acquiring radiation source data, calculating initial radiation field distribution based on the radiation source data, obtaining preliminary dose deposition values of all spatial points in a monitored area, analyzing inter-particle interactions reflected by the preliminary dose deposition values, and determining energy redistribution coefficients corresponding to all the spatial points; Step S2, utilizing the energy redistribution coefficient to carry out point-by-point adjustment on the preliminary dose deposition value to obtain corrected three-dimensional dose concentration field data, and calculating gradient vectors among all spatial points based on the corrected three-dimensional dose concentration field data to obtain an optimized gradient distribution diagram; S3, extracting a boundary point set with gradient change exceeding a preset threshold value from the optimized gradient distribution diagram, constructing an isosurface grid based on the boundary point set to form a three-dimensional isosurface structure, acquiring time sequence radiation source data, and iteratively updating the energy redistribution coefficient to obtain a time sequence evolution sequence of a concentration field; And S4, identifying and marking a high-concentration region based on the time sequence evolution sequence of the concentration field, fusing the high-concentration region and the three-dimensional isosurface structure, and generating a visualized distribution model reflecting the dose concentration according to a fusion result. In a sec