A method and system for detecting a radio source device transaction based on posture analysis

Abstract

The application discloses a kind of based on gesture analysis's radioactive source equipment abnormal transaction detection method and system, specifically related to the measurement field of radiation beam current, including obtaining the attitude data of radioactive source equipment, by analyzing attitude data and constructing beam principal axis vector to express the radiation propagation direction of radioactive source equipment under current attitude;By constructing radiation action region and judging its space inclusion relation to identify whether the irradiation area of radioactive source under current attitude is out of bounds;By analyzing the detection data to construct the measured intensity principal axis vector and compare with beam principal axis vector to identify whether the radiation propagation direction is shifted. By constructing the mapping relationship between device attitude and beam principal axis, and jointly judging the out-of-bounds state of space action region and the shift state of radiation propagation direction, the accurate identification of posture-induced radiation directionality abnormality is realized, so as to effectively solve the problem that the traditional system cannot perceive beam out-of-control path.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of radiation beam current measurement, and more particularly to a radioactive source equipment abnormality detection method and system based on attitude analysis. BACKGROUND

[0002] In the existing radiation monitoring technology system, the operation safety of the radioactive source equipment usually depends on the real-time monitoring of the output radiation intensity and the dose threshold control in the preset area. In the radiation intensity measurement device, fixed point or area array detectors are widely used to perform dose coverage sensing on the environment around the radiation source. The judgment is mainly focused on whether the radiation intensity exceeds the set upper limit.

[0003] However, in actual application, the radioactive source equipment may be slightly detached, tilted or turned over due to non-human factors such as loose base, mechanical impact, vibration transmission or structural fatigue. Although such abnormality does not cause significant changes in the source power or total dose, it may still maintain the appearance of "output dose compliance", but it substantially changes the spatial direction of the beam, causing the radiation energy to be directly directed to the area that is not originally protected or the space where personnel activities exist.

[0004] Since such "radiation direction drift type abnormality" does not change the dose value at the measurement point, the traditional dose monitoring system will completely fail to identify it, thereby forming a high-risk, invisible and non-alarm mechanism radiation safety blind area.

[0005] Therefore, the problem in the current radiation measurement system is that the existing radiation intensity-based measurement system cannot identify the spatial radiation direction deviation caused by the abnormal attitude of the radioactive source equipment, and thus cannot effectively detect the direction out-of-control type abnormality. SUMMARY

[0006] In order to overcome the above-mentioned defects of the prior art, the embodiments of the present application provide a radioactive source equipment abnormality detection method and system based on attitude analysis, which constructs the mapping relationship between the device attitude and the beam main axis, and jointly judges the out-of-bound state of the space action area and the direction deviation state of the radiation propagation, to accurately identify the attitude-induced radiation direction abnormality, thereby effectively solving the problem that the traditional system cannot sense the beam out-of-control path.

[0007] To achieve the above-mentioned purpose, the present application provides the following technical scheme: a radioactive source equipment abnormality detection method based on attitude analysis, comprising:

[0008] S1, acquiring attitude data of a radioactive source equipment, constructing a beam main axis vector to express the radiation propagation direction of the radioactive source equipment under the current attitude by analyzing the attitude data;

[0009] S2, identifying whether the irradiation area of the radiation source in the current posture exceeds the boundary by constructing the radiation action area and judging the spatial inclusion relationship thereof;

[0010] S3, identifying whether the radiation propagation direction is deviated by constructing the measured intensity principal axis vector from the detection data and comparing the same with the beam principal axis vector;

[0011] S4, generating the abnormal event label by combining the spatial action area exceeding boundary label and the radiation direction deviation state to represent the radiation pointing abnormality caused by the posture change.

[0012] In a preferred embodiment, in S1, further comprising: acquiring posture data by a posture measurement component of the radiation source device, the posture data including three direction angle parameters of pitch angle, yaw angle and roll angle, the posture data being used to represent the spatial orientation state of the device;

[0013] performing three-axis rotation transformation on the posture data to solve the exit surface normal vector of the device at the current time, and taking the normal vector as the initial construction reference of the beam propagation direction; constructing the beam principal axis vector of the radiation source at the current time by combining the normal vector and the preset source structure parameters, the beam principal axis vector being used to represent the propagation trend of the radiation direction in the three-dimensional space; the preset source structure parameters including the exit cone angle and the symmetry calibration value.

[0014] In a preferred embodiment, in S2, further comprising: jointly solving the geometric envelope structure of the radiation action area by the beam principal axis vector and the exit cone angle of the radiation, the geometric envelope structure including the central axis, the conical opening angle and the boundary extension range;

[0015] performing spatial inclusion judgment on the radiation action area constructed at the current time and the preset target irradiation area, judging whether the current radiation action area is completely covered within the irradiation permitted boundary; if the judgment result is not included, generating the spatial action area exceeding boundary label, otherwise, judging the current radiation action area as covering effective and maintaining the current irradiation state.

[0016] In a preferred embodiment, in S3, further comprising: collecting detection data by an intensity detection point array arranged at a spatial fixed position, the detection data including the spatial coordinates and the corresponding measured radiation dose values of each detection point, the detection data being used to describe the spatial distribution state of the current radiation intensity;

[0017] performing spatial interpolation on the detection data and performing directional principal component extraction to construct the measured intensity principal axis vector, the measured intensity principal axis vector representing the central direction feature corresponding to the upper limit direction of the energy density in the current radiation intensity distribution;

[0018] The direction included angle calculation is performed between the measured intensity principal axis vector and the beam principal axis vector, if the included angle is greater than a preset deviation threshold, the current radiation direction deviation state is identified.

[0019] In a preferred embodiment, in S4, further comprising: judging whether a space action area boundary label has been generated at the current time, and whether a radiation direction deviation state has been identified, if both conditions are true, it is determined that the current is a posture change induced space radiation pointing abnormal state; otherwise, it is determined to be a non-abnormal state.

[0020] The posture data associated with the posture change induced space radiation pointing abnormal state, the beam principal axis vector and the measured intensity principal axis vector are constructed into an abnormal event element set; the abnormal event element set is converted into an abnormal event label, the abnormal event label is used to represent the radiation propagation direction deviation and its space action influence induced by the radiation source device under certain posture change conditions.

[0021] In a preferred embodiment, in S1, further comprising: obtaining the dynamic change of the attitude angle and the angular velocity, taking the disturbance function thereof as the path integral kernel, and solving the three-dimensional beam principal axis vector under the exit cone angle modulation, so as to build a posture disturbance beam principal axis construction model:

[0022]

[0023] Where θ p (τ), θ y (τ) and θ r (τ) are time functions of pitch angle θ p , yaw angle θ y and roll angle θ r ; ω p (τ) represents the time derivative of pitch angle θ p at time point τ; ω y (τ) represents the time derivative of yaw angle θ y at time point τ; ω r (τ) represents the time derivative of roll angle θ r at time point τ; β p is the attitude symmetry bias value in the pitch direction; β y is the attitude symmetry bias value in the yaw direction; β r is the attitude symmetry bias value in the roll direction; α is the exit cone angle of the radiation source; n x , n y , n z represent the initial calibration normal components of the device exit surface along the X, Y and Z axis directions in the three-dimensional coordinate system; t is the current calculation time; represents the three-dimensional beam principal axis vector formed under the action of posture disturbance.

[0024] In a preferred embodiment, in S2, further comprising: constructing a geometric drift ratio and a maximum out-of-bound amplitude of the irradiation region volume according to the offset of the beam principal axis in the three axes of space;

[0025]

[0026] X c =C z +V x

[0027] Y c =C y +V y

[0028] Z c =C z +V z

[0029] Δ 最大越界 =max{(|X c -C x |-R x ),(|Y c -C y |-R y ),(|Z c -C z |

[0030] -R z})

[0031] wherein C x , C y , C z represent the center point coordinates of the permitted irradiation region in three-dimensional space; R x , R y , R z represent the permitted boundary radius of the irradiation region in each spatial axis direction; V x , V y , V z represent the components in X, Y, Z three axis directions obtained by ; X c , Y c , Z c represent the spatial center point coordinates of the irradiation region in X, Y, Z three axis directions at the current time; the offset volume ratio represents the volume ratio of the irradiation region spatial offset in the permitted region; Δ 最大越界 is the upper limit out-of-bound distance of the irradiation region in any direction at the current time.

[0032] In a preferred embodiment, in S3, further comprising: extracting the principal gradient direction from the measured radiation intensity distribution function, and calculating the unit directional angle with the theoretical principal axis vector, for determining the directional drift; defining φ to represent the angle between the unit vectors, expressed as:

[0033]

[0034] wherein the angle φ between the unit vectors is used to measure the degree of deviation of the radiation propagation direction; arccos(·) is the inverse cosine function, in the above formula, the inverse cosine function is used to calculate the angle between two vectors according to the ratio of the vector dot product and the modulus; represents the three-dimensional intensity variation rate vector calculated according to the spatial coordinate difference of the detection points in X, Y, Z directions based on the measured dose values of the multiple radiation detection points in space; represents the dot product between the radiation intensity gradient vector and the beam principal axis vector; represents the Euclidean norm of the radiation intensity gradient vector; represents the Euclidean norm of the beam principal axis vector.

[0035] In a preferred embodiment, in S4, further comprising: based on the directional deviation angle and the deviation volume ratio, constructing a joint event scoring function:

[0036]

[0037] wherein δ 阈值 is the directional deviation angle threshold; ∈ is the response adjustment parameter; ε is a very small positive value, which is used to prevent the logarithmic term from being 0; E 事件 represents the score value of the posture-induced radiation pointing anomaly event;

[0038] The response adjustment parameter ∈ is expressed as:

[0039]

[0040] wherein θ i (τ) represents the time function of the pitch angle θ p , the yaw angle θ y and the roll angle θ r , reflecting the dynamic change state of the direction posture angle of the radiation source device at time τ; p, y, r correspond to the pitch direction, the yaw direction and the roll direction, respectively; represents the instantaneous change rate of the posture angle at time τ; β i represents the symmetry calibration value in the source structure parameter; t is the current calculation time.

[0041] A radiation source device abnormality detection system based on posture analysis, comprising a posture decomposition module, a region boundary judgment module, an offset identification module, and an abnormality calibration module.

[0042] The posture decomposition module is configured to obtain posture data of the radiation source device, parse the posture data, and construct a beam principal axis vector to represent the radiation propagation direction of the radiation source device in the current posture.

[0043] The region boundary judgment module is configured to construct a radiation action region and judge the spatial inclusion relationship thereof to identify whether the irradiation region of the radiation source in the current posture is out of boundary.

[0044] The offset identification module is configured to parse detection data, construct a measured intensity principal axis vector, and compare the measured intensity principal axis vector with the beam principal axis vector to identify whether the radiation propagation direction is offset.

[0045] The abnormality calibration module is configured to generate an abnormality event label by combining the spatial action region out-of-boundary label and the radiation direction offset state to represent the radiation pointing abnormality caused by the posture change.

[0046] Technical effects and advantages of the present application:

[0047] 1. The beam principal axis vector is constructed based on the pitch angle, the yaw angle, and the roll angle to identify the radiation pointing drift problem caused by the device posture change, and the technical defect that the traditional dose monitoring cannot detect directional abnormality is solved.

[0048] 2. The geometric envelope structure of the current irradiation region is constructed by the principal axis vector and the exit cone angle, and the spatial inclusion judgment is performed to identify whether the current radiation is out of boundary to the non-permitted region, and the structural identification of the out-of-space risk is realized.

[0049] 3. The intensity principal axis vector is reconstructed by analyzing the dose data of the spatial detection point, and the angle comparison is performed with the theoretical principal axis to identify whether the actual energy propagation direction is out of control, and the identification ability of the non-dose type abnormality is improved.

[0050] 4. The abnormality event label is generated only when the spatial out-of-boundary state and the direction offset state are both established by joint judgment, false positives are excluded, and the abnormality identification result with structural consistency and high confidence is ensured. BRIEF DESCRIPTION OF DRAWINGS

[0051] Fig. 1 The figure is a method step framework diagram of the present application.

[0052] Fig. 2 The figure is a system module schematic diagram of the present application. DETAILED DESCRIPTION

[0053] With reference to the drawings of the embodiments of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described, obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present application.

[0054] With reference to the drawings of the embodiments of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described, obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present application. Figs. 1-2 According to the method, the posture analysis-based radioactive source equipment abnormality detection method comprises the following steps.

[0055] S1, acquiring posture data of a radioactive source equipment, constructing a beam principal axis vector by analyzing the posture data to express a radiation propagation direction of the radioactive source equipment in a current posture;

[0056] S2, judging whether the irradiation area of the radioactive source in the current posture is out of boundary by constructing a radiation action area and judging a spatial inclusion relationship thereof;

[0057] S3, constructing a measured intensity principal axis vector by analyzing detection data and comparing the measured intensity principal axis vector with the beam principal axis vector to identify whether the radiation propagation direction is deviated;

[0058] S4, generating an abnormality event label by combining a spatial action area out-of-boundary label and a radiation direction deviation state to represent a radiation pointing abnormality caused by a posture change.

[0059] In S1, the posture data is acquired by a posture measurement component of the radioactive source equipment, the posture data comprises three direction angle parameters of a pitch angle, a yaw angle and a roll angle, and the posture data is used to represent a spatial orientation state of the equipment; in addition, in actual application, the posture measurement component of the radioactive source equipment comprises but is not limited to an inertial measurement assembly composed of an accelerometer, a gyroscope and a magnetometer.

[0060] The posture data is subjected to three-axis rotation transformation to solve an outgoing surface normal vector of the equipment at a current time, and the normal vector is taken as an initial construction reference of the beam propagation direction; a beam principal axis vector of the radioactive source at the current time is constructed by combining the normal vector with preset source structure parameters, and is used to represent a propagation trend of the radiation direction in a three-dimensional space; the preset source structure parameters comprise an outgoing cone angle and a symmetry calibration value.

[0061] In S2, the geometric envelope structure of the radiation action area is solved by the beam principal axis vector and the outgoing cone angle of the radiation, and the geometric envelope structure comprises a central axis, a conical opening angle and a boundary extension range.

[0062] The current time is constructed and the preset target irradiation area is executed to perform a spatial inclusion judgment, and whether the current radiation effect area is completely covered in the irradiation permission boundary is judged; if the judgment result is not included, a spatial effect area out-of-bound label is generated, otherwise the current radiation effect area is determined to be covered effectively, and the current irradiation state is maintained.

[0063] In S3, it also includes: collecting detection data through the intensity detection point array arranged in the fixed position of the space, the detection data including the spatial coordinates of each detection point and the corresponding measured radiation dose value, the detection data being used to describe the spatial distribution state of the current radiation intensity;

[0064] The detection data is subjected to spatial interpolation and directional principal component extraction to construct a measured intensity principal axis vector, the measured intensity principal axis vector representing the central directional characteristics corresponding to the energy density upper limit direction in the current radiation intensity distribution;

[0065] The measured intensity principal axis vector and the beam principal axis vector are executed to perform a directional angle calculation, and if the angle is greater than a preset deviation threshold, it is identified as a current radiation direction deviation state.

[0066] In S4, it also includes: judging whether the spatial effect area out-of-bound label has been generated at the current time, and whether the radiation direction deviation state has been identified, if both conditions are established, it is determined that the current is a spatial radiation pointing abnormal state caused by the posture change; otherwise, it is determined to be a non-abnormal state;

[0067] The posture data associated with the posture change induced spatial radiation pointing abnormal state, the beam principal axis vector and the measured intensity principal axis vector are constructed into an abnormal event element set; the abnormal event element set is converted into an abnormal event label, the abnormal event label being used to represent the radiation propagation direction deviation and its spatial effect caused by the radiation source device under the specific posture change condition.

[0068] It should be noted that for the formula structure involved in the present scheme, the dimensionless term can be used as a proportional or structural adjustment factor, and when combined with a quantity with a unit, it only plays a numerical scaling role and does not introduce a new physical dimension, so it will not change or confuse the unit system of the whole expression; such combination of "dimensionless term and quantity unit term" can be understood as a composite structure expression form commonly used in mathematical and physical modeling, which conforms to the principle of dimensional consistency and has a clear physical interpretation basis;

[0069] Secondly, in the formula structure of the scheme, if multiple variable terms with different physical units are involved, including but not limited to time, mass or energy variables, their joint occurrence is to express the cooperative modeling relationship of multiple physical mechanisms. Each variable can be mapped by a function, combined by a ratio or adjusted by a normalization to form a unified structure with clear units and explicit meaning, and the overall expression conforms to the principle of dimensional consistency and the general norm of engineering modeling.

[0070] In the scheme, if constants, weights, adjustment factors, threshold parameters, and proportion coefficients are designed, they are adjustable control parameters for different application environments, and their values depend on the target device configuration, data input characteristics, and performance optimization goals. In the implementation stage, they are set within a reasonable range through model verification, performance constraints, or engineering calibration. Although these parameters do not have a unique value, they have clear adjustment logic and calculation paths, and belong to the deterministic setting process in engineering implementation. The purpose of such setting is to ensure that the scheme has both general adaptability and reproducibility and operability, without affecting its technical clarity and implementability.

[0071] In S1, the dynamic changes of the attitude angle and angular velocity are also included. The disturbance function of the path integral kernel is solved to obtain the three-dimensional beam principal axis vector under the modulation of the exit cone angle, and the attitude disturbance beam principal axis construction model is built:

[0072]

[0073] where θ p (τ), θ y (τ), and θ r (τ) are the time functions of the pitch angle θ p , the yaw angle θ y , and the roll angle θ r , respectively; ω p (τ) represents the time derivative of the pitch angle θ p at time point τ, i.e., the pitch angle change rate, which is used to represent the instantaneous attitude dynamic strength of the device in that direction; ω y (τ) represents the time derivative of the yaw angle θ y at time point τ, i.e., the yaw angle change rate, which is used to reflect the change of the rotation rate of the device around the vertical axis with time; ω r (τ) represents the time derivative of the roll angle θ r at time point τ, i.e., the roll angle change rate, which is used to describe the attitude change sensitivity of the device in the roll direction; in summary, ω p (τ), ω y (τ), and ω r (τ) are the time derivatives (angular velocities) of the corresponding angles; β pis a yaw direction attitude symmetry bias value, the yaw direction attitude symmetry bias value is used to represent a tendency of the main axis direction to drift caused by structural asymmetry of the device around the vertical axis; β y is a yaw direction attitude symmetry bias value, the yaw direction attitude symmetry bias value is used to represent a tendency of the main axis direction to drift caused by structural asymmetry of the device around the vertical axis; β r is a roll direction attitude symmetry bias value, the roll direction attitude symmetry bias value is used to reflect a deviation in the angle offset response caused by uneven structural distribution when the device is laterally rotated; α is an emission cone angle of the radiation source; n x , n y , n z represents initial calibration normal components of the device emission surface in the X, Y, and Z axis directions in the three-dimensional coordinate system, and the initial calibration normal components of the three axis directions are used to define a spatial reference of the beam propagation direction in the undisturbed state; t is a current calculation time; represents a three-dimensional beam main axis vector formed under the action of the attitude disturbance.

[0074] In S2, the geometric drift ratio and the maximum out-of-bound amplitude of the radiation action region volume are also constructed according to the offset of the beam main axis in the spatial three-axis.

[0075]

[0076] X c =C z +V x

[0077] Y c =C y +V y

[0078] Z c =C z +V z

[0079] Δ 最大越界 =max{(|X c -C x |-R x ),(|Y c -C y |-R y ),(|Z c -C z |-R z )

[0080] wherein C x , C y , C z represent the center point coordinates of the permitted irradiation region in the three-dimensional space; R x , Ry , R z denotes the allowable boundary radius of the irradiation region in each spatial axis direction; V x , V y , V z respectively denote the components in X, Y, Z three axis directions obtained by , which are used to describe the spatial propagation direction of the beam in the three-dimensional coordinate system; X c , Y c , Z c respectively denote the spatial center point coordinates of the irradiation region in X, Y, Z three axis directions at the current moment; the offset volume ratio represents the volume ratio of the spatial offset of the irradiation region to the permitted region; Δ 最大越界 is the upper limit out-of-bound distance of the irradiation region in any direction at the current moment.

[0081] In S3, it further includes: extracting the main gradient direction from the measured radiation intensity distribution function, and calculating the unit direction angle with the theoretical main axis vector, which is used to determine the direction drift; define φ to represent the angle between the unit vectors, which is represented as:

[0082]

[0083] Wherein the angle φ between the unit vectors is used to measure the degree of offset of the radiation propagation direction; arccos(·) is the inverse cosine function, which is used to calculate the angle between two vectors according to the ratio of vector dot product and modulus in the above formula; denotes the three-dimensional intensity change rate vector calculated based on the measured dose values of multiple radiation detection points in space, according to the spatial coordinate differences of the detection points in X, Y, Z three directions, which is used to represent the main rising direction of the current radiation energy in space; denotes the dot product between the radiation intensity gradient vector and the beam main axis vector, which is used to reflect the consistency degree of the two in the direction; denotes the Euclidean norm of the radiation intensity gradient vector, that is, the modulus of the gradient vector, which is used to quantify the spatial intensity of the direction gradient; denotes the Euclidean norm of the beam main axis vector, that is, the modulus of the main axis vector, which is used to quantify the spatial amplitude of the theoretical propagation direction.

[0084] In S4, it further includes: based on the direction offset angle and the offset volume ratio, a joint event scoring function is constructed:

[0085]

[0086] Wherein δ 阈值is a direction deviation angle determination threshold; ∈ is a response adjustment parameter; ε is a minimum positive value for preventing the logarithmic term from being 0; E 事件 represents a score value of a posture-induced radiation direction change event;

[0087] The response adjustment parameter ∈ is represented as:

[0088]

[0089] where θ i (τ) represents the time function of the pitch angle θ p , the yaw angle θ y and the roll angle θ r , reflecting the dynamic change state of the direction posture angle of the radiation source device at time τ; p, y and r correspond to the pitch direction, the yaw direction and the roll direction, respectively; represents the instantaneous change rate of the posture angle at time τ, that is, the posture angle velocity; θ i includes the pitch angle θ p , the yaw angle θ y and the roll angle θ r ; β i represents a symmetry calibration value in the source structure parameter, which is used to describe the static angular deviation of each posture direction of the radiation source device relative to the ideal exit axis in the structural design, and is used to correct the structural alignment difference between the posture angle and the exit direction; t is the current calculation time; the summation symbol ∑ i∈{p,y,r} represents the sum of the integral contributions of the pitch direction p, the yaw direction y and the roll direction r; in addition, all the integrals in the formula are path integrals, representing the disturbance response energy or the joint gradient influence accumulated from the initial state (τ=0) to the current time (τ=t).

[0090] The radiation source device anomaly detection system based on posture analysis includes a posture deconstruction module, a region judgment module, an offset identification module and an anomaly calibration module.

[0091] The posture deconstruction module is used to obtain the posture data of the radiation source device, and the beam principal axis vector is constructed by analyzing the posture data to express the radiation propagation direction of the radiation source device under the current posture.

[0092] The region judgment module is used to identify whether the irradiation region of the radiation source under the current posture is out of boundary by constructing the radiation action region and judging the spatial inclusion relationship thereof.

[0093] The offset identification module is used to identify whether the radiation propagation direction is deviated by constructing the measured intensity principal axis vector by analyzing the detection data and comparing it with the beam principal axis vector.

[0094] The abnormality calibration module is configured to generate an abnormality event label by combining a space action area boundary crossing label and a radiation direction deviation state to represent a radiation direction abnormality caused by a change in the posture.

[0095] Overall, the present scheme is used for systematic reflection on the significant blind spot of traditional radiation monitoring technology in identifying spatial directionality abnormalities; in the operation of existing radioactive source equipment, mainstream dose monitoring systems mainly rely on fixed points or area array detectors to monitor whether the output intensity exceeds the threshold; however, such systems cannot effectively identify the beam propagation direction drift problem caused by changes in the device posture; under non-human intervention, the radioactive source equipment may change its geometric posture in the pitch, yaw, or roll direction due to reasons such as loose base, impact vibration, or structural fatigue, although the dose output does not exceed the standard and the total power is not abnormal, the beam may be directed to an area that was not originally protected, forming a blind area irradiation risk with extremely high hidden dangers;

[0096] To solve this structural technical gap, the present scheme constructs a set of posture analysis-based radioactive source equipment abnormality detection method, breaking through the traditional "dose over-limit alarm" logic, taking "whether the beam propagation direction is out of control" as the core evaluation standard, and realizing real-time abnormality identification starting from the modeling of the device state;

[0097] The first step is to deconstruct the posture to obtain the pitch, yaw, and roll three-axis angle parameters and their angular velocity changes of the radioactive source equipment, analyze the orientation of the equipment in three-dimensional space through a posture measurement component (such as an inertial measurement unit), and construct a beam principal axis vector based on the exit cone angle and structural symmetry calibration value; this vector not only reflects the theoretical direction of the current radiation direction, but also serves as the main reference axis for subsequent space action area and intensity distribution comparison;

[0098] The second step is to determine the boundary of the area, that is, based on the principal axis vector and the exit cone angle, to construct the geometric envelope structure of the radiation action area at the current time, calculate the relationship between the space center point and the permitted irradiation boundary, and output whether there is a space boundary crossing conclusion; this judgment not only focuses on whether it "crosses the boundary", but also quantifies the spatial scale changes such as the offset volume ratio as input factors for subsequent abnormality joint scoring;

[0099] The third step is to identify the deviation by obtaining real-time radiation intensity data from the detection points arranged in the fixed position in the scene, performing spatial interpolation and gradient analysis, extracting the energy density direction principal component to construct the measured intensity principal axis vector, and then calculating the unit angle between the beam principal axis; if the direction angle exceeds the deviation threshold, it is determined that the current radiation direction has deviated out of control; the calculation in this part not only integrates the physical gradient field information, but also captures the substantial propagation deviation hidden under the normal surface dose through principal axis alignment calculation;

[0100] The fourth step is the transaction calibration. The transaction calibration is the fusion determination of the logical intersection of the second step and the third step: only when the spatial action area is out of bounds and the radiation propagation direction has been shifted, a real radiation direction change caused by attitude change is constituted; this determination avoids false positives (for example, short-term and slight attitude fluctuations do not cause actual shifts) and enhances the recognition confidence; once the transaction determination condition is met, the attitude data, the theoretical main axis vector and the measured direction vector corresponding to the time are constructed as an event element set, a transaction event label is generated, and the abnormal example is tracked, analyzed and alarmed.

[0101] The purpose of the present scheme is to fundamentally break the existing dose threshold-centered discrimination mode, to pre-position the abnormality recognition to the device physical attitude disturbance level, to actively model the beam runaway path, to improve the identification ability of the "invisible beam shift", and to ultimately guarantee the safety of personnel and the controllability of the radiation of the target area; the method has engineering adaptability and lays a data structure foundation for subsequent beam stability modeling and intelligent diagnosis of the radiation source device.

[0102] The above only describes the preferred embodiments of the present application and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for detecting anomalies in radioactive source equipment based on attitude analysis, characterized in that, The method comprises the following steps: S1, acquiring attitude data of the radioactive source device, and constructing a beam principal axis vector by analyzing the attitude data to represent the radiation propagation direction of the radioactive source device in the current attitude; S2, constructing a radiation action region and judging its spatial inclusion relationship to identify whether the irradiation region of the radioactive source in the current attitude is out of bounds; S3, constructing a measured intensity principal axis vector by analyzing the detection data and comparing it with the beam principal axis vector to identify whether the radiation propagation direction is deviated; S4, generating an abnormal event label by combining the spatial action region out-of-bounds label and the radiation direction deviation state to represent the radiation pointing abnormality caused by the attitude change; In S2, the beam principal axis vector and the exit cone angle of the radiation are used to jointly solve the geometric envelope structure of the radiation action region, which includes a central axis, a conical opening angle, and a boundary extension range; The current time radiation action region is compared with the preset target irradiation region to perform spatial inclusion judgment, and it is judged whether the current radiation action region is completely covered within the irradiation permission boundary; if the judgment result is not included, a spatial action region out-of-bounds label is generated, otherwise the current radiation action region is determined to be covered effectively, and the current irradiation state is maintained; In S4, it is further judged whether the spatial action region out-of-bounds label has been generated at the current time, and whether the radiation direction deviation state has been identified, if both conditions are true, it is determined that the current is a spatial radiation pointing abnormality state caused by the attitude change; otherwise, it is determined to be a non-abnormal state; The attitude data, the beam principal axis vector, and the measured intensity principal axis vector associated with the spatial radiation pointing abnormality state caused by the attitude change are constructed into an abnormal event element set; the abnormal event element set is converted into an abnormal event label, which is used to represent the radiation propagation direction deviation caused by the radioactive source device under specific attitude change conditions and its spatial action influence; In S3, the main gradient direction is extracted from the measured radiation intensity distribution function, and the unit direction angle between the theoretical principal axis vector is calculated to determine the direction drift; define Ф as the angle between the unit vectors, which is represented as: wherein the included angle Ф between the unit vectors is used to measure the degree of deviation of the radiation propagation direction; arccos(·) is the inverse cosine function, and in the above formula, the inverse cosine function is used to calculate the included angle between two vectors according to the ratio of the vector dot product to the modulus; represents a three-dimensional intensity change rate vector calculated according to the spatial coordinate difference of the detection points in the X, Y and Z directions on the basis of the measured dose values of the plurality of radiation detection points in the space; represents the dot product between the radiation intensity gradient vector and the beam principal axis vector; represents the Euclidean norm of the radiation intensity gradient vector; represents the Euclidean norm of the beam principal axis vector; The measured intensity principal axis vector represents the central direction characteristic corresponding to the upper limit direction of the energy density in the current radiation intensity distribution.

2. The radioactive source device abnormality detection method based on attitude analysis according to claim 1, wherein: In S1, the attitude data is acquired by the attitude measurement component of the radioactive source device, and the attitude data includes three direction angle parameters, i.e., pitch angle, yaw angle, and roll angle, which represent the spatial orientation state of the device; The three-axis rotation transformation is performed on the attitude data to solve the exit surface normal vector of the device at the current time, and the normal vector is used as the initial construction reference of the beam propagation direction; the beam principal axis vector of the radioactive source at the current time is constructed by combining the normal vector and the preset source structure parameters, which represents the propagation trend of the radiation direction in the three-dimensional space; the preset source structure parameters include the exit cone angle and the symmetry calibration value.

3. The radioactive source device abnormality detection method based on attitude analysis according to claim 2, wherein: In S3, further comprising: collecting detection data by an array of intensity detection points arranged at fixed positions in space, the detection data including spatial coordinates of each detection point and corresponding measured radiation dose values, the detection data being used to describe the spatial distribution state of the current radiation intensity; Performing spatial interpolation on the detection data and performing directional principal component extraction to construct a measured intensity principal axis vector; Performing directional angle calculation on the measured intensity principal axis vector and the beam principal axis vector, and if the angle is greater than a preset deviation threshold, identifying a current radiation direction deviation state.

4. The radiation source device abnormality detection method based on attitude analysis according to claim 3, further comprising: In S1, further comprising: obtaining the dynamic changes of the attitude angle and the angular velocity, taking the disturbance function thereof as a path integral kernel, and solving the three-dimensional beam principal axis vector under the emission cone angle modulation to build an attitude disturbance beam principal axis construction model: where θ p , )θ y , )θ r , are time functions of the pitch angle θ p , the yaw angle θ y and the roll angle θ r , respectively; ω p ( ) represents the time derivative of the pitch angle θ p at the time point ; ω y ( ) represents the time derivative of the yaw angle θ y at the time point ; ω r ( ) represents the time derivative of the roll angle θ r at the time point ; β p is the attitude symmetry bias in the pitch direction; β y is the attitude symmetry bias in the yaw direction; β r is the attitude symmetry bias in the roll direction; α is the emission cone angle of the radiation source; n x , n y , n z represent the initial calibration normal components of the device emission surface along the X, Y, Z three-axis directions in the three-dimensional coordinate system; t is the current calculation time; represents the three-dimensional beam principal axis vector formed under the action of the attitude disturbance.

5. The radiation source device abnormality detection method based on attitude analysis according to claim 4, further comprising: In S2, further comprising: constructing a geometric drift ratio of the radiation affected area volume and a maximum out-of-bound amplitude according to the offset of the beam principal axis in the spatial three-axis; X c = C x + V x Y c = C y + V y Z c = C z + V z Δ 最大越界 = max{(|X c -C x |-R x ), (|Y c -C y |-R y ), (|Z c -C z |-R z )} wherein C x , C y , C z represents the coordinate of the center point of the permitted irradiation region in three-dimensional space; R x , R y , R z represents the permitted boundary radius of the irradiation region in each spatial axis direction; V x , V y , V z respectively represent the components in the X, Y, Z three axis directions obtained by ; X c , Y c , Z c respectively represent the spatial center point coordinates of the radiation action region in the X, Y, Z three axis directions at the current moment; the offset volume ratio represents the volume ratio of the spatial offset of the radiation action region to the permitted region; Δ 最大越界 is the upper limit out-of-bound distance of the radiation action region in any direction at the current moment.

6. The radiation source device abnormality detection method based on attitude analysis according to claim 5, further comprising: In S4, further comprising: constructing a joint event scoring function based on the directional deviation angle and the offset volume ratio: where δ 阈值 is a direction offset angle decision threshold; ∈ is a response adjustment parameter; ε is a small positive value to prevent the logarithm term from being zero; E 事件 represents a score value for a posture-induced radiation pointing anomaly event; The response adjustment parameter ∈ is expressed as: where θ i ( ) denotes the time function of the pitch angle θ p , the yaw angle θ y and the roll angle θ r , reflecting the dynamic change state of the directional attitude angle of the radiation source device at time ; p, y, r correspond to the pitch direction, the yaw direction, and the roll direction, respectively; denotes the instantaneous change rate of the attitude angle at time ; β i denotes the symmetry calibration value in the source structure parameters; t is the current calculation time.

7. A system for detecting a change of a radioactive source device based on posture analysis, the system performing the method for detecting a change of a radioactive source device based on posture analysis according to claim 6. The system comprises an attitude deconstruction module, a region boundary judgment module, an offset identification module, and an abnormality calibration module. The attitude deconstruction module is used to obtain attitude data of the radiation source device, and to express the radiation propagation direction of the radiation source device under the current attitude by analyzing the attitude data and constructing a beam principal axis vector. The region boundary judgment module is used to identify whether the irradiation region of the radiation source under the current attitude is out of bound by constructing a radiation affected region and judging the spatial inclusion relationship thereof. The offset identification module is used to identify whether the radiation propagation direction is deviated by constructing a measured intensity principal axis vector from the analysis of the detection data and comparing it with the beam principal axis vector. The abnormality calibration module is used to generate an abnormality event label by combining the out-of-bound label of the spatial affected region and the radiation direction deviation state to represent the radiation pointing abnormality caused by the attitude change.

Images