High-energy electron accelerator safety interlocking method, device, equipment and medium

Abstract

The invention relates to the field of safety interlocking control, in particular to a high-energy electron accelerator safety interlocking method, device and equipment and a medium. The method comprises the following steps: acquiring and analyzing historical electronic information (such as radiation dose and beam parameters) and historical equipment load corresponding to the electron accelerator, deeply knowing the behavior mode of the accelerator in the past working state, calculating the historical equipment stability degree, and determining critical electronic information by combining the historical data, so as to obtain the electronic information of the accelerator. The method is beneficial for taking measures in advance to avoid potential dangers, monitoring current electronic information in real time and comparing the current electronic information with critical information, quickly and accurately judging whether the accelerator is about to enter an unstable or dangerous state or not, and immediately generating a safety trigger instruction once the working state is determined to be adjusted to a trigger state. Personnel safety and stable operation of equipment are effectively guaranteed, and the accuracy of safety interlocking triggering of the high-energy electron accelerator is remarkably improved.

Description

Technical Field

[0001] The present application relates to the field of safety interlock control, and particularly to a safety interlock method, device, equipment and medium for a high-energy electron accelerator. Background Art

[0002] As an important scientific research and industrial equipment, high-energy electron accelerators play an important role in particle physics research, medical treatment and other fields. With the progress of technology, high-energy electron accelerators are more and more widely used, and higher requirements are also put forward for their performance and safety. In order to ensure the safety of operators and the validity of experimental data, effective safety measures must be taken to prevent accidents.

[0003] At present, the safety control of high-energy electron accelerators mainly relies on traditional hardware protection mechanisms, such as physical shielding and interlock switches. These methods can reduce risks to a certain extent, but in the face of complex dynamic change environments, especially when precise triggering of safety interlocks is required, traditional methods seem powerless.

[0004] Specifically, when the existing safety control system detects and responds to various abnormal situations that may occur during the operation of the electron accelerator, the single threshold setting method may not comprehensively cover all potential risk points, resulting in the inability to trigger the safety interlock in a timely manner in some cases, thus increasing the risk of safety accidents. Therefore, how to improve the accuracy of triggering the safety interlock of high-energy electron accelerators has become an urgent technical problem to be solved. Summary of the Invention

[0005] In order to improve the accuracy of triggering the safety interlock of high-energy electron accelerators, the present application provides a safety interlock method, device, equipment and medium for a high-energy electron accelerator.

[0006] In a first aspect, the present application provides a safety interlock method for a high-energy electron accelerator, adopting the following technical solutions: A safety interlock method for a high-energy electron accelerator includes: Obtaining historical electron information and historical equipment load corresponding to the electron accelerator, and real-time monitoring of the current electron information and current equipment information of the electron accelerator, where the historical electron information includes electron information corresponding to a historical period, and the historical period is the time period between the working state of the electron accelerator changing from a normal state to a triggered state, and the electron information includes radiation dose and beam parameters; Calculating the historical equipment stability degree corresponding to the electron accelerator; Analyzing the historical electron information, the historical equipment load and the historical equipment stability degree to determine the critical electron information corresponding to the electron accelerator; Based on the critical electronic information and the current electronic information, determine whether the working state of the electron accelerator is adjusted from the normal state to the trigger state; If it is determined that the working state of the electron accelerator is adjusted from the normal state to the trigger state, a safety trigger instruction is generated.

[0007] By adopting the above technical solution, obtaining and analyzing the historical electronic information (such as radiation dose, beam current parameters) and historical equipment load corresponding to the electron accelerator can deeply understand the behavior pattern of the accelerator in the past working state, especially the key features before the conversion between normal operation and the trigger state. Secondly, calculating the historical equipment stability provides an important basis for evaluating the overall performance and potential risks of the accelerator. Then, combining these historical data to determine the critical electronic information helps to take measures in advance to avoid potential dangers. Real-time monitoring of the current electronic information and comparing it with the critical information can quickly and accurately determine whether the accelerator is about to enter an unstable or dangerous state. Once it is determined that the working state needs to be adjusted to the trigger state, a safety trigger instruction is immediately generated, effectively ensuring the safety of personnel and the stable operation of the equipment. In summary, through the comprehensive application of historical data analysis, real-time monitoring and warning mechanisms, this method significantly improves the accuracy of the safety interlock trigger of the high-energy electron accelerator, thereby further enhancing the safety performance and operation efficiency of the electron accelerator.

[0008] In a possible implementation manner, calculate the historical equipment stability corresponding to the electron accelerator. The calculation process of the historical equipment stability corresponding to each historical period includes any one of the following: Obtain each historical position and the historical preset orbit of the electron accelerator, and determine whether each historical position is located on the historical preset orbit; if there is a historical position not on the historical preset orbit, determine the historical position not on the historical preset orbit as the target historical position, and calculate the historical deviation between each target historical position and the historical preset orbit. Obtain the deviation allowable range and the historical confidence level corresponding to the electron accelerator, calculate the confidence interval corresponding to the average value of the historical deviation based on the historical confidence level, and calculate the historical equipment stability corresponding to the electron accelerator based on the formula, where the formula is , S is the historical equipment stability, is the maximum value of the confidence interval corresponding to the average value of the historical deviation, is the minimum value of the deviation allowable range, is the maximum value of the deviation allowable range; Obtain the historical beam current intensity and historical beam energy corresponding to each historical moment of the electron accelerator during the historical period, determine the beam current intensity fluctuation and beam energy fluctuation, perform frequency-domain analysis on the beam current intensity fluctuation and beam energy fluctuation respectively, establish the beam current intensity frequency-domain relationship and beam energy fluctuation frequency relationship, and establish the beam current intensity amplitude matrix and beam energy amplitude matrix. Based on the beam current intensity frequency-domain relationship, the beam energy fluctuation frequency relationship, the beam current intensity amplitude matrix, and the beam energy amplitude matrix, determine the historical equipment stability degree corresponding to the electron accelerator.

[0009] By adopting the above technical solutions, through the comprehensive application of two methods, namely position deviation analysis and beam parameter frequency-domain analysis, the historical equipment stability degree of the electron accelerator is comprehensively and deeply evaluated not only from the perspective of physical position accuracy but also from the perspective of beam control efficiency. This multi-dimensional stability evaluation strategy provides a solid foundation for timely detecting potential faults, optimizing equipment operation strategies, and ensuring the long-term safe operation of the accelerator.

[0010] In a possible implementation manner, analyze the historical electronic information, the historical equipment load, and the historical equipment stability degree to determine the critical electronic information corresponding to the electron accelerator, including: Obtain the equipment load threshold corresponding to the electron accelerator, and based on the equipment load threshold, screen out at least one first historical period from multiple historical periods, and determine the historical periods other than the first historical period as the second historical period. The historical equipment load corresponding to the first historical period does not exceed the equipment load threshold, and the historical equipment stability degree corresponding to the first historical period is not less than the equipment stability degree threshold; Based on the beam parameters corresponding to each first historical period, establish the fault time-domain relationship and beam time-domain relationship corresponding to the electron accelerator; Screen out the historical moments when the historical equipment load of the electron accelerator exceeds the equipment load threshold from each second historical period as the second historical moments, and obtain the second historical electronic information corresponding to each second historical moment from the historical electronic information; Perform frequency-domain analysis on the beam parameters and radiation doses corresponding to the second historical electronic information respectively to obtain the beam frequency-domain relationship and radiation frequency-domain relationship corresponding to the electron accelerator; Based on the fault time-domain relationship, beam time-domain relationship, beam frequency-domain relationship, and radiation frequency-domain relationship corresponding to the electron accelerator, determine the critical electronic information corresponding to the electron accelerator.

[0011] By adopting the above technical solutions, by setting the equipment load threshold and the equipment stability threshold, the first historical period representing the operation state of the equipment under a lower load and relatively stable is screened out from the historical periods, and based on the beam parameters of the first historical period, the fault time-domain relationship and the beam time-domain relationship are established, which helps to deeply understand the operation law and potential fault characteristics of the equipment in the normal state. For the second historical period, especially those historical moments when the equipment load exceeds the threshold, by obtaining the corresponding second historical electronic information at these moments and performing frequency-domain analysis on the beam parameters and radiation dose, the beam frequency-domain relationship and the radiation frequency-domain relationship are obtained, thereby revealing the variation law of the beam and radiation parameters in the frequency domain under the high-load state of the equipment, and providing key information in the frequency domain for identifying the critical state. Finally, by synthesizing the fault time-domain relationship, the beam time-domain relationship, the beam frequency-domain relationship, and the radiation frequency-domain relationship, the critical information that the equipment is about to enter an unstable or dangerous state can be identified more accurately.

[0012] In a possible implementation manner, the beam parameters include the particle passing speed and the particle passing quantity. The establishment of the fault time-domain relationship and the beam time-domain relationship corresponding to the electron accelerator based on the beam parameters corresponding to each of the first historical periods includes: Determine the historical duration corresponding to each of the first historical periods, and establish the fault time-domain relationship corresponding to the electron accelerator based on each historical duration; Based on the historical electronic information, obtain the particle passing quantity and the particle passing speed corresponding to each of the first historical periods, and establish the beam time-domain relationship corresponding to the electron accelerator based on the particle passing quantity and the particle passing speed corresponding to each of the first historical periods.

[0013] By adopting the above technical solutions, the historical duration corresponding to each first historical period is determined. By clarifying the time range of each period, the change of the equipment operation state and time can be more precisely correlated, so as to capture the timing characteristics of potential faults. And based on the historical electronic information, the particle passing quantity and the particle passing speed corresponding to each first historical period are extracted. Based on the particle passing quantity and the particle passing speed corresponding to each first historical period, the beam time-domain relationship corresponding to the electron accelerator is established, forming a time-beam change curve, which not only shows the fluctuation of the beam parameters with time, but also reveals the stability and periodic characteristics of the beam under the normal operation state of the equipment.

[0014] In a possible implementation manner, based on the fault time-domain relationship, the beam time-domain relationship, the beam frequency-domain relationship, and the radiation frequency-domain relationship corresponding to the electron accelerator, determining the critical electronic information corresponding to the electron accelerator includes: Establish a decision tree model based on the fault time domain relationship and beam current time domain relationship corresponding to the electron accelerator; Based on the decision tree model, predict the abnormal state of the electron accelerator in the current cycle, where the abnormal state of the equipment is abnormal or non-abnormal, and the abnormal state includes equipment load abnormality and / or equipment stability degree abnormality; When the abnormal state of the electron accelerator in the current cycle is non-abnormal, determine the first critical value based on the beam current frequency domain relationship and the radiation frequency domain relationship, and determine the second critical value based on the fault time domain relationship and the beam current time domain relationship. Based on the first critical value and the second critical value, determine the critical electron information corresponding to the electron accelerator; When the abnormal state of the electron accelerator in the current cycle is abnormal, determine the correlation degree between the equipment load and the equipment stability degree. Based on the correlation degree, the fault time domain relationship and the beam current time domain relationship, determine the first target critical value, and based on the beam current frequency domain relationship and the radiation frequency domain relationship, determine the second target critical value.

[0015] By adopting the above technical solution, a decision tree model is established based on the fault time domain relationship and beam current time domain relationship of the electron accelerator, and the decision tree model is used to predict the abnormal state of the electron accelerator in the current cycle, thereby improving the prediction accuracy. When the prediction result is non-abnormal, the solution further determines the first critical value based on the beam current frequency domain relationship and the radiation frequency domain relationship, and at the same time determines the second critical value by combining the fault time domain relationship and the beam current time domain relationship. The setting of these two critical values provides a quantitative standard for judging whether the equipment is close to the critical state. When the prediction result is abnormal, analyze the correlation degree between the equipment load and the equipment stability degree, determine the first target critical value, and at the same time determine the second target critical value by combining the information of the beam current frequency domain and the radiation frequency domain.

[0016] In a possible implementation manner, determining the correlation degree between the equipment load and the equipment stability degree includes: Sort the historical equipment load and the historical equipment stability degree respectively according to the historical start time corresponding to each historical cycle to obtain a load sequence and a stability sequence; Based on the load sequence and the stability sequence, establish a first regression equation of the equipment stability degree with respect to its own lag term, and a second regression equation of the equipment stability degree with respect to its own lag term and the equipment load lag term; Based on the first regression equation and the second regression equation, determine the causal relationship direction between the equipment load and the equipment stability degree; Calculate the correlation coefficient between the equipment load and the equipment stability degree, and combine the causal relationship direction and the correlation coefficient to determine the correlation degree between the equipment load and the equipment stability degree.

[0017] By adopting the above technical solution, the historical equipment loads and the historical equipment stability degrees are sorted according to the historical start time corresponding to each historical period, and a load sequence and a stability sequence are obtained respectively. Based on the load sequence and the stability sequence, two regression equations are established: the first regression equation of the equipment stability degree with respect to its own lag term, and the second regression equation of the equipment stability degree with respect to its own lag term and the lag term of the equipment load. By comparing indexes such as the goodness of fit or significance test of the first regression equation and the second regression equation, the causal relationship direction between the equipment load and the equipment stability degree is determined, which helps to clarify who affects whom, and the correlation coefficient between the equipment load and the equipment stability degree is calculated. Combining the causal relationship direction, the association degree between the two is comprehensively evaluated. The correlation coefficient quantifies the linear correlation degree between the two, while the causal relationship direction provides information about the action direction. The combination of the two makes the evaluation result more comprehensive and accurate.

[0018] In a possible implementation manner, determining whether the working state of the electron accelerator is adjusted from the normal state to the trigger state based on the critical electronic information and the current electronic information includes: Comparing the current radiation dose with the critical radiation dose, and comparing the current beam parameter with the critical beam parameter to determine whether the current radiation dose exceeds the critical radiation dose and whether the current beam parameter exceeds the critical beam parameter; If the current radiation dose exceeds the critical radiation dose, or / and, the current beam parameter exceeds the critical beam parameter, it is determined that the working state of the electron accelerator is adjusted from the normal state to the trigger state; If the current radiation dose does not exceed the critical radiation dose and the current beam parameter does not exceed the critical beam parameter, it is determined that the working state of the electron accelerator remains the normal state.

[0019] By adopting the above technical solution, the current radiation dose is compared with the critical radiation dose to evaluate whether the current radiation level has reached or exceeded the threshold that may cause safety risks. At the same time, the current beam parameters are compared with the critical beam parameters to determine whether the beam state is abnormal. According to the comparison results, if the current radiation dose exceeds the critical radiation dose, or the current beam parameters exceed the critical beam parameters (or both exceed), this means that the operating state of the electron accelerator has exceeded the safe range. Therefore, the solution determines that the working state of the electron accelerator should be adjusted from the normal state to the trigger state. If the current radiation dose does not exceed the critical radiation dose and the current beam parameters also do not exceed the critical beam parameters, this indicates that the electron accelerator is still within the safe and stable operating state. Therefore, it is determined that its working state remains the normal state. This judgment helps to maintain the continuous and stable operation of the equipment, while avoiding unnecessary intervention and shutdown.

[0020] In a second aspect, the present application provides a safety interlock device for a high-energy electron accelerator, adopting the following technical solution: A safety interlock device for a high-energy electron accelerator, comprising: An acquisition module, configured to acquire the historical electronic information and historical equipment load corresponding to the electron accelerator, and to monitor in real time the current electronic information and current equipment information of the electron accelerator. The historical electronic information includes the electronic information corresponding to a historical period, and the historical period is the time period between the working state of the electron accelerator changing from the normal state to the trigger state. The electronic information includes the radiation dose and the beam parameters; A calculation module, configured to calculate the historical equipment stability degree corresponding to the electron accelerator; An analysis module, configured to analyze the historical electronic information, the historical equipment load, and the historical equipment stability degree to determine the critical electronic information corresponding to the electron accelerator; A determination module, configured to determine whether the working state of the electron accelerator is adjusted from the normal state to the trigger state based on the critical electronic information and the current electronic information; A generation module, configured to generate a safety trigger instruction if it is determined that the working state of the electron accelerator is adjusted from the normal state to the trigger state.

[0021] In a third aspect, the present application provides an electronic device, adopting the following technical solution: An electronic device, which includes: At least one processor; A memory; At least one application program, where the at least one application program is stored in a memory and configured to be executed by at least one processor, and the at least one application program is configured to: execute the high-energy electron accelerator safety interlock method described in the first aspect above.

[0022] In a fourth aspect, the present application provides a computer-readable storage medium, adopting the following technical solution: A computer-readable storage medium, including: a computer program stored with the ability to be loaded and executed by a processor for the high-energy electron accelerator safety interlock method described in the first aspect above.

[0023] In summary, the present application includes the following beneficial technical effects: By acquiring and analyzing the historical electronic information (such as radiation dose, beam current parameters) corresponding to the electron accelerator and the historical equipment load, it is possible to deeply understand the behavior pattern of the accelerator in the past working state, calculate the historical equipment stability degree, which provides an important basis for evaluating the overall performance and potential risks of the accelerator. Combining these historical data to determine the critical electronic information helps to take measures in advance to avoid potential dangers. Real-time monitoring of the current electronic information and comparing it with the critical information can quickly and accurately determine whether the accelerator is about to enter an unstable or dangerous state. Once it is determined that the working state needs to be adjusted to the trigger state, a safety trigger instruction is immediately generated, effectively ensuring the safety of personnel and the stable operation of the equipment, and significantly improving the accuracy of the high-energy electron accelerator safety interlock trigger. Description of the Drawings

[0024] Figure 1 is a schematic flowchart of a high-energy electron accelerator safety interlock method provided by an embodiment of the present application; Figure 2 is a schematic block diagram of a high-energy electron accelerator safety interlock device provided by an embodiment of the present application; Figure 3 is a schematic diagram of an electronic device provided by an embodiment of the present application. Detailed Embodiments

[0025] The following is a further detailed description of the present application in conjunction with the attached Figure 1 - attached Figure 3 to further illustrate the present application.

[0026] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are some, but not all, of the embodiments of the present application. All other embodiments obtained by those of ordinary skill in the art without creative efforts based on the embodiments in the present application belong to the scope of protection of the present application.

[0027] To facilitate the understanding of the technical solution proposed in this application, several elements introduced in the description of this application are introduced here first. It should be understood that the following introduction only facilitates the understanding of these elements in order to understand the content of the embodiments of this application, and does not necessarily cover all possible situations.

[0028] A high-energy electron accelerator uses the forces of electric and magnetic fields to propel electrons and achieves the required energy and speed through multiple accelerations and focusing. In the accelerator, electrons are injected into the acceleration structure and obtain an accelerating force by applying an electric or magnetic field, thereby reaching a high-energy state.

[0029] A safety interlock is an integrated system that ensures the correct operation and interlocking of equipment through interlocking devices installed between different devices or systems to guarantee the safety of operators and equipment. It uses mechanical or electrical control to sequentially connect each relevant instrument and equipment to make them interact with each other, achieving the purpose of protecting equipment and ensuring safe production.

[0030] The safety interlock plays a crucial role in a high-energy electron accelerator. It mainly monitors and controls the operating state of the accelerator in real time through a series of interconnected devices and systems. When an abnormal situation or potential danger occurs in the accelerator, the safety interlock can respond quickly and take necessary protective measures, such as cutting off the power supply, stopping the beam current, etc., to prevent accidents and ensure the safety of personnel and equipment.

[0031] The embodiments of this application provide a safety interlock method for a high-energy electron accelerator. As Figure 1 shown, in the method provided in the embodiments of this application, it is executed by an electronic device, which can be a server or a terminal device. Among them, the server can be an independent physical server, a server cluster or a distributed system composed of multiple physical servers, or a cloud server providing cloud computing services. The terminal device can be a smart phone, a tablet computer, a notebook computer, a desktop computer, etc., but is not limited thereto. The terminal device and the server can be directly or indirectly connected through wired or wireless communication methods. This application does not limit this here. The method includes steps S101 - step S105, where: Step S101: Obtain the historical electronic information and historical equipment load corresponding to the electron accelerator, and monitor the current electronic information and current equipment information of the electron accelerator in real time.

[0032] Among them, the historical electronic information includes the electronic information corresponding to the historical period, where the historical period is the time period between the working state of the electron accelerator changing from the normal state to the trigger state. The electronic information includes radiation dose and beam parameters. Specifically, the radiation dose is a physical quantity measuring the radiation intensity around the electron accelerator, usually expressed in specific units (such as gray, sievert, etc.). Excessively high radiation dose may cause harm to personnel and the environment. The beam parameters are parameters describing the characteristics of the electron beam, including beam intensity, beam energy, beam position, etc. These parameters affect the performance and application effect of the electron accelerator.

[0033] Specifically, the electronic device communicates with the data storage system of the electron accelerator and extracts historical electronic information and historical device load data from the storage database using a data interface protocol (such as HTTP, TCP / IP, etc.). For real-time monitoring, the electronic device is connected to various sensors installed on the electron accelerator (such as radiation dose sensors, beam parameter sensors, device load sensors, etc.) and reads the data output by the sensors at a certain time interval (such as every second, every minute, etc.) to obtain the current electronic information and current device information.

[0034] Step S102: Calculate the historical device stability degree corresponding to the electron accelerator.

[0035] Among them, the historical device stability degree is a quantitative index comprehensively reflecting the stability of the device state of the electron accelerator during the historical operation process. The higher the value, the more stable the device operation and the less affected by various factors.

[0036] In this embodiment, the calculation process of the historical device stability degree corresponding to each historical period can be implemented by Method 1 or Method 2, where: Method 1: Obtain each historical position and the historical preset orbit of the electron accelerator, and determine whether each historical position is on the historical preset orbit; if there is a historical position not on the historical preset orbit, then determine the historical position not on the historical preset orbit as the target historical position, and calculate the historical deviation between each target historical position and the historical preset orbit, obtain the deviation allowable range and historical confidence level corresponding to the electron accelerator, calculate the confidence interval corresponding to the average value of the historical deviation based on the historical confidence level, and calculate the historical device stability degree corresponding to the electron accelerator based on the formula, where the formula is , S is the historical device stability degree, is the maximum value of the confidence interval corresponding to the average value of the historical deviation, is the minimum value of the deviation allowable range, is the maximum value of the deviation allowable range.

[0037] Specifically, by interacting with the positioning system or related monitoring devices of the electron accelerator, each historical position information corresponding to the beam current of the electron accelerator in each historical cycle is obtained from the database of the electron accelerator. At the same time, the historical preset orbit information is read from the database of the electron accelerator. Further, each obtained historical position is compared with the historical preset orbit. For each historical position, it is judged whether it is on the historical preset orbit. If there are historical positions that are not on the historical preset orbit, these positions are marked as target historical positions.

[0038] Furthermore, for each target historical position, the distance between it and the historical preset orbit is calculated to obtain the historical deviation corresponding to each target historical position, and the deviation allowable range and historical confidence level (such as a 95% confidence level) corresponding to the electron accelerator are obtained from the database of the electron accelerator. This range stipulates the acceptable upper and lower limits of the position deviation of the electron accelerator.

[0039] Based on the obtained historical confidence level, a suitable statistical method (such as a confidence interval calculation method based on normal distribution) is used to calculate the confidence interval of the average historical deviation, and the maximum value of the confidence interval corresponding to the average value of the above-obtained historical deviation, the minimum value of the deviation allowable range, and the maximum value of the deviation allowable range are substituted into the formula , and the historical equipment stability degree corresponding to the electron accelerator is calculated. Where S is the historical equipment stability degree, is the maximum value of the confidence interval corresponding to the average value of the historical deviation, is the minimum value of the deviation allowable range, is the maximum value of the deviation allowable range.

[0040] Method 2: Obtain the historical beam current intensity and historical beam energy corresponding to each historical moment of the electron accelerator in the historical cycle, determine the beam current intensity fluctuation situation and the beam energy fluctuation situation, perform frequency domain analysis on the beam current intensity fluctuation situation and the beam energy fluctuation situation respectively to establish the beam current intensity frequency domain relationship and the beam energy fluctuation frequency relationship, and establish the beam current intensity amplitude matrix and the beam energy amplitude matrix. Based on the beam current intensity frequency domain relationship, the beam energy fluctuation frequency relationship, the beam current intensity amplitude matrix, and the beam energy amplitude matrix, determine the historical equipment stability degree corresponding to the electron accelerator.

[0041] Obtain the historical beam current intensity and historical beam energy data corresponding to each historical moment of the electron accelerator in each historical cycle from the database of the electron accelerator, and form a curve that changes with time, that is, the beam current intensity fluctuation situation and the beam energy fluctuation situation. Further, frequency domain analysis can be performed on the obtained curves of the beam current intensity fluctuation situation and the beam energy fluctuation situation. Specifically, the Fourier transform method can be used to convert the fluctuation curve in the time domain into a frequency domain representation, so as to obtain the beam current intensity frequency domain relationship and the beam energy fluctuation frequency relationship, and reveal the distribution of beam parameters at different frequencies.

[0042] Furthermore, according to the results of the frequency domain analysis, a beam current intensity amplitude matrix and a beam energy amplitude matrix are established respectively, and the stability of the beam parameters is evaluated by calculating the concentration degree of the amplitude distribution of each frequency component and the proportion of a specific frequency component, and then the stability degree of the corresponding historical equipment of the electron accelerator is determined. Among them, the elements in the matrix represent the amplitude sizes corresponding to different frequency components, reflecting the contribution degree of each frequency component to the beam parameter fluctuation.

[0043] Specifically, after performing frequency domain analysis on beam parameters (such as beam current intensity, particle passing speed), the amplitude values corresponding to each frequency component are obtained. In order to calculate the concentration degree of the amplitude distribution, first calculate the mean value of these amplitude values, and then, for each amplitude value, calculate the square of its deviation from the mean value, and sum all the squared deviations and divide by the number of amplitude values to obtain the variance. The smaller the variance, the more concentrated the amplitude distribution, and the more stable the beam parameters in the frequency domain.

[0044] Obtain the frequencies that do not belong to the frequency domain threshold from all frequencies as specific frequencies, obtain the amplitude values corresponding to each specific frequency, sum the amplitude values of all specific frequencies to obtain the sum of the amplitude values of specific frequencies, and at the same time, calculate the sum of the amplitude values of all frequencies to obtain the sum of the amplitude values of all frequencies. Calculate the proportion of the specific frequency component, that is, the ratio of the sum of the amplitude values of specific frequencies to the sum of the amplitude values of all frequencies. The lower the proportion, the smaller the contribution of the beam parameter fluctuation at these specific frequencies, and the more stable the beam parameters are reflected from the perspective of specific frequencies.

[0045] Step S103: Analyze the historical electronic information, historical equipment load, and historical equipment stability degree to determine the critical electronic information corresponding to the electron accelerator.

[0046] Among them, the critical electronic information is the critical value range of electronic information during the operation of the electron accelerator. When the actual electronic information exceeds this range, it may mean that the working state of the accelerator is abnormal and safety measures need to be triggered.

[0047] Specifically, data analysis methods can be used to determine the critical electron information of the electron accelerator. More specifically, time-domain analysis can be first performed on historical electron information, historical equipment load, and historical equipment stability to observe the changing trends and fluctuation rules of these data over time; then frequency-domain analysis can be carried out to find out the periodic components and characteristic frequencies in the data. At the same time, correlation analysis (such as calculating Pearson correlation coefficient, Spearman correlation coefficient, etc.) is used to determine the correlation and causal relationship between various factors. According to the analysis results, combined with the characteristics of electron information in the normal state and abnormal state within the historical period, critical values are set for electron information such as radiation dose and beam parameters, and these critical values constitute the critical electron information.

[0048] Step S104: Based on the critical electron information and the current electron information, determine whether the working state of the electron accelerator is adjusted from the normal state to the trigger state.

[0049] Among them, the working state is a description of the state of the electron accelerator during operation, which is divided into the normal state and the trigger state. The normal state means that the accelerator operates stably according to the design requirements, and all parameters are within the safe range; the trigger state means that the accelerator has an abnormal situation, which may cause harm to the equipment, personnel or environment, and corresponding safety measures need to be taken.

[0050] Specifically, the currently monitored current electron information is compared with the critical electron information. For each electron information parameter (such as radiation dose, beam intensity, etc.), it is judged whether it exceeds the corresponding critical value range. If one or more parameters exceed the critical value and the exceeding situation lasts for a preset duration (to avoid misjudgment caused by instantaneous fluctuations, the preset duration can be set according to the actual situation or historical experience, and the embodiments of the present application do not limit this), it is considered that the working state of the electron accelerator is adjusted from the normal state to the trigger state.

[0051] Step S105: If it is determined that the working state of the electron accelerator is adjusted from the normal state to the trigger state, generate a safety trigger instruction.

[0052] After the electronic device determines that the working state of the electron accelerator enters the trigger state, a safety trigger instruction is generated according to the preset instruction format. Among them, the safety trigger instruction contains relevant information about the trigger state (such as the triggered electron information parameters, trigger time, etc.), and then the instruction is sent to the control system of the electron accelerator through the communication interface, and the control system executes corresponding safety measures according to the instruction (such as cutting off the power supply, closing the protection door, etc.).

[0053] An embodiment of the present application provides a safety interlock method for a high-energy electron accelerator. By obtaining and analyzing the historical electron information (such as radiation dose, beam current parameters) and historical equipment load corresponding to the electron accelerator, it is possible to deeply understand the behavior pattern of the accelerator in the past working state. Calculating the historical equipment stability provides an important basis for evaluating the overall performance and potential risks of the accelerator. Combining these historical data to determine the critical electron information helps to take measures in advance to avoid potential hazards. Real-time monitoring of the current electron information and comparing it with the critical information can quickly and accurately determine whether the accelerator is about to enter an unstable or dangerous state. Once it is determined that the working state needs to be adjusted to the trigger state, a safety trigger instruction is immediately generated, effectively ensuring the safety of personnel and the stable operation of the equipment, and significantly improving the accuracy of the safety interlock trigger of the high-energy electron accelerator.

[0054] A possible implementation manner of the embodiment of the present application. In the above step S103, analyzing the historical electron information, historical equipment load, and historical equipment stability to determine the critical electron information corresponding to the electron accelerator includes: Obtain the equipment load threshold corresponding to the electron accelerator, and based on the equipment load threshold, screen out at least one first historical period from multiple historical periods, and determine the historical periods other than the first historical period as the second historical period. The historical equipment load corresponding to the first historical period does not exceed the equipment load threshold, and the historical equipment stability corresponding to the first historical period is not less than the equipment stability threshold; Based on the beam current parameters corresponding to each first historical period, establish the fault time domain relationship and beam current time domain relationship corresponding to the electron accelerator; Screen out the historical moments when the historical equipment load of the electron accelerator exceeds the equipment load threshold from each second historical period as the second historical moments, and obtain the second historical electron information corresponding to each second historical moment from the historical electron information; Perform frequency domain analysis on the beam current parameters and radiation dose corresponding to the second historical electron information respectively to obtain the beam current frequency domain relationship and radiation frequency domain relationship corresponding to the electron accelerator; Based on the fault time domain relationship, beam current time domain relationship, beam current frequency domain relationship, and radiation frequency domain relationship corresponding to the electron accelerator, determine the critical electron information corresponding to the electron accelerator.

[0055] Specifically, read the equipment load threshold and equipment stability threshold corresponding to the electron accelerator from the database corresponding to the electron accelerator. Then, the electronic device traverses the data of multiple historical periods. For each historical period, check whether the corresponding historical equipment load does not exceed the equipment load threshold and whether the historical equipment stability is not less than the equipment stability threshold. If these two conditions are met, mark this historical period as the first historical period; the historical periods that do not meet the conditions are marked as the second historical period.

[0056] Further, extract the beam parameter data within each first historical period, including data such as beam intensity and beam energy that vary with time. For the fault time-domain relationship, analyze the correlation between possible fault conditions (such as beam interruption and abnormal parameter fluctuations) and time within the first historical period. For example, statistically analyze information such as the time points and durations when faults occur, and summarize the occurrence patterns of faults in the time dimension. For the beam time-domain relationship, plot the curves of beam parameters changing with time, and analyze the time-changing trends of beam parameters in the normal and stable operating states, such as the periodic fluctuations of beam intensity and the steady changes of beam energy, so as to establish the corresponding relationship between beam parameters and time.

[0057] For each second historical period, check the historical equipment load data at each historical moment one by one. When it is found that the historical equipment load at a certain historical moment exceeds the equipment load threshold, mark this historical moment as the second historical moment. Then, based on these second historical moments, extract the corresponding second historical electronic information at the corresponding moments from the stored historical electronic information, including beam parameters and radiation dose. Further, perform frequency-domain analysis on the beam parameter and radiation dose data in the second historical electronic information. For beam parameters, the frequency-domain analysis method of Fourier transform can be used to convert the time-domain signal of beam parameters changing with time into a frequency-domain signal, and analyze the distribution of beam parameters at different frequencies. For radiation dose data, also perform frequency-domain conversion and analysis to determine the frequency characteristics of radiation dose fluctuations and establish the radiation frequency-domain relationship.

[0058] Even further, comprehensively analyze the fault time-domain relationship, beam time-domain relationship, beam frequency-domain relationship, and radiation frequency-domain relationship. For beam parameters, combine the changing trends in the normal state in the beam time-domain relationship and the frequency characteristics of abnormal fluctuations in the beam frequency-domain relationship to determine the critical change ranges of beam parameters in the time domain and frequency domain, such as the maximum fluctuation amplitude of beam intensity and the maximum amplitude at specific frequencies. For radiation dose, determine the critical value and fluctuation range of radiation dose based on the radiation frequency-domain relationship and the possible accompanying radiation anomalies in the fault time-domain relationship. These determined critical ranges and values of beam parameters and radiation dose constitute the corresponding critical electronic information of the electron accelerator.

[0059] In a possible implementation manner of the embodiment of the present application, in the above embodiment, the beam parameters include the particle passing speed and the number of particles passing through. Based on the beam parameters corresponding to each first historical period, establish the fault time-domain relationship and beam time-domain relationship corresponding to the electron accelerator, including: Determine the historical duration corresponding to each first historical period, and establish the fault time-domain relationship corresponding to the electron accelerator based on each historical duration; Based on historical electronic information, obtain the number of particles passing through and the particle passing speed corresponding to each first historical period, and establish the beam current time-domain relationship corresponding to the electron accelerator based on the number of particles passing through and the particle passing speed corresponding to each first historical period.

[0060] Among them, the fault time-domain relationship is used to describe the time-domain distribution of the number of faults occurring in the electron accelerator within a certain period of time. The beam current time-domain relationship is used to describe the time-domain distribution of the number of particles passing through the electron accelerator within a certain period of time, and shows the variation law of the beam current intensity with time.

[0061] After obtaining the first historical period, extract the start time and end time of each first historical period from the stored first historical period data, and calculate the difference between the end time and the start time to obtain the historical duration corresponding to each first historical period. For example, if the start time of a first historical period is 9:00 and the end time is 10:00, then the historical duration of this period is 1 hour. Further, within each first historical period, count the number of faults occurring according to the fault record information. And with the historical duration as the horizontal axis and the number of faults occurring as the vertical axis, mark the historical duration and the number of faults occurring corresponding to each first historical period, so as to establish the time-domain distribution of the number of faults occurring within a certain period of time, that is, the fault time-domain relationship.

[0062] Screen out the data on the number of particles passing through and the particle passing speed within each first historical period from the historical electronic information. For each first historical period, divide its historical duration into several equal time intervals. For example, divide a 1-hour historical period into 60 1-minute time intervals. Within each divided time interval, the representative value of the number of particles passing through in this interval can be obtained by summing or averaging the number of particles passing through data recorded in this time interval, so as to obtain the number of particles passing through.

[0063] Further, with time as the horizontal axis and the number of particles passing through as the vertical axis, mark the number of particles passing through corresponding to each time interval. Analyze the distribution of these data points, such as the change trend of the number of particles passing through with time, whether there are periodic fluctuations, etc., so as to establish the time-domain distribution of the number of particles passing through within a certain period of time, that is, the beam current time-domain relationship.

[0064] In a possible implementation manner of the embodiments of the present application, in the above embodiments, based on the fault time-domain relationship, beam current time-domain relationship, beam current frequency-domain relationship, and radiation frequency-domain relationship corresponding to the electron accelerator, determine the critical electron information corresponding to the electron accelerator, including: Based on the fault time-domain relationship and beam current time-domain relationship corresponding to the electron accelerator, establish a decision tree model; Based on the decision tree model, predict the abnormal state of the electron accelerator in the current cycle. The abnormal state of the equipment is abnormal or non-abnormal, and the abnormal state includes equipment load abnormality and / or equipment stability abnormality. When the abnormal state of the electron accelerator in the current cycle is non-abnormal, determine the first critical value based on the beam current frequency domain relationship and the radiation frequency domain relationship, and determine the second critical value based on the fault time domain relationship and the beam current time domain relationship. Based on the first critical value and the second critical value, determine the critical electron information corresponding to the electron accelerator. When the abnormal state of the electron accelerator in the current cycle is abnormal, determine the correlation degree between the equipment load and the equipment stability degree. Based on the correlation degree, the fault time domain relationship, and the beam current time domain relationship, determine the first target critical value, and determine the second target critical value based on the beam current frequency domain relationship and the radiation frequency domain relationship.

[0065] After obtaining the fault time domain relationship and the beam current time domain relationship, the number of fault occurrences and the fault occurrence time in the fault time domain relationship, as well as the number of particles passing through and the time in the beam current time domain relationship, etc., can be integrated. These data are organized into data samples suitable for constructing a decision tree according to a preset format. Each data sample contains multiple features (the number of faults at different time points, the number of particles passing through) and corresponding class labels (whether the equipment is abnormal). Use a suitable decision tree algorithm (such as ID3, C4.5, or CART algorithm) to construct a decision tree model based on the data samples and the features contained in the data samples. The decision tree model will recursively divide the data set according to the features until certain stopping conditions are met (such as all samples in the node belong to the same class, the depth of the tree reaches a preset value, etc.).

[0066] Furthermore, obtain the current duration from the start time to the current time in the current cycle and the beam current-related information (such as the number of particles passing through, the particle passing speed), and organize these information in the data format when constructing the decision tree model, and provide it as input data to the decision tree model. The decision tree model makes inferences according to the structure of the tree and the judgment conditions of the nodes. Starting from the root node, select the appropriate branch according to the feature values of the input data, and gradually traverse the decision tree downward until reaching the leaf node. The class label corresponding to the leaf node is the predicted abnormal state of the equipment (abnormal or non-abnormal).

[0067] When the equipment abnormal state of the electron accelerator in the current cycle is non-abnormal, the beam current frequency-domain relationship and the radiation frequency-domain relationship can be analyzed. For the beam current frequency-domain relationship, observe the distribution and variation of beam current parameters at different frequencies, and determine the normal fluctuation range of beam current parameters in the frequency domain, such as the maximum amplitude of the beam current intensity at a specific frequency. For the radiation frequency-domain relationship, also find the normal fluctuation range of the radiation dose in the frequency domain. Determine the boundary values of these normal fluctuation ranges in the frequency domain as the first critical value. According to the fault time-domain relationship and the beam current time-domain relationship, analyze the variation range of the number of faults occurring and the number of particles passing through under normal conditions. For example, count the maximum number of faults occurring, the maximum and minimum values of the number of particles passing through in the historical cycle. Determine the boundary values of these normal variation ranges in the time domain as the second critical value. Combine the first critical value and the second critical value, and combine the critical values of the beam current parameters and the radiation dose in the time domain and the frequency domain to form the critical electron information corresponding to the electron accelerator. Among them, the first critical value reflects the normal fluctuation boundary of the beam current parameters and the radiation dose in the frequency domain, and is used to judge whether there is an abnormality in the frequency domain; the second critical value reflects the boundary of the normal variation range of the number of faults occurring and the number of particles passing through in the time domain, and is used to judge whether there is an abnormality in the time domain.

[0068] When the equipment abnormal state of the electron accelerator in the current cycle is abnormal, relevant analysis methods (such as Pearson correlation coefficient, Spearman correlation coefficient) can be used to determine the degree of association between the equipment load and the equipment stability. Specifically, collect the corresponding values of the equipment load and the equipment stability in the historical data, and then calculate the correlation coefficient between the two according to the selected relevant analysis method. The closer the absolute value of the correlation coefficient is to 1, the stronger the degree of association; the closer it is to 0, the weaker the degree of association. Further, analyze by combining the degree of association between the equipment load and the equipment stability, as well as the fault time-domain relationship and the beam current time-domain relationship. Consider the influence of the degree of association on the occurrence of faults and the beam current situation under the equipment abnormal state. For example, if the degree of association is strong, it means that the abnormality of the equipment load may greatly affect the equipment stability and the beam current parameters.

[0069] Specifically, based on the historical start time and the historical end time of each historical cycle, calculate the historical duration corresponding to each historical cycle, and calculate the average value and the standard deviation of the historical duration corresponding to the historical time period. Based on the above-mentioned degree of association, use the first formula: , calculate the historical duration adjustment factor, where is the historical duration adjustment factor, is the duration influence coefficient, which can be determined according to the actual situation and is used to measure the influence weight of the degree of association on the number of faults occurring, is the degree of association. After obtaining the historical duration adjustment factor, multiply the historical duration adjustment factor by the historical duration mean and the historical duration standard deviation respectively to obtain the updated historical duration mean and the updated historical duration standard deviation. Obtain the duration confidence level corresponding to the electron accelerator, and calculate the first target critical value corresponding to the historical duration based on the duration confidence level, the updated historical duration mean, and the updated historical duration standard deviation. More specifically, based on the second formula: , calculate the first target critical value corresponding to the historical duration, where is the first target critical value corresponding to the historical duration, is the updated historical duration mean, is the duration confidence level, is the updated historical duration standard deviation.

[0070] Similarly, calculate the number of particles passing through corresponding to each historical period, and calculate the mean value and standard deviation of the number of particles passing through corresponding to the historical period. Based on the above degree of association and the third formula: , calculate the particle passing number adjustment factor, where is the particle passing number adjustment factor, is the quantity influence coefficient, which can be determined according to the actual situation and is used to measure the influence weight of the degree of association on the number of faults, is the degree of association. After obtaining the particle passing number adjustment factor, multiply the particle passing number adjustment factor by the particle passing number mean and the particle passing number standard deviation respectively to obtain the updated particle passing number mean and the updated particle passing number standard deviation. Obtain the particle passing number confidence level corresponding to the electron accelerator, and calculate the first target critical value corresponding to the particle passing number based on the particle passing number confidence level, the updated particle passing number mean, and the updated particle passing number standard deviation. More specifically, based on the fourth formula: , calculate the first target critical value corresponding to the particle passing number, where is the first target critical value corresponding to the particle passing number, is the updated particle passing number mean, is the particle passing number confidence level, is the updated particle passing number duration standard deviation.

[0071] Based on the first target critical value corresponding to the particle passing number and the first target critical value corresponding to the historical duration, obtain the first target critical value.

[0072] Further, analyze the beam current frequency-domain data, calculate its spectral characteristics (such as power spectral density), and use the fast Fourier transform (FFT) to convert the beam current data in the time domain into frequency-domain data to obtain the changes in the beam current spectrum under abnormal device states. Similarly, perform similar processing on the radiation frequency-domain data, calculate the power spectral density of the radiation, analyze the energy distribution of the radiation in the frequency domain, and obtain the changes in the radiation intensity at specific frequencies in the radiation spectrum when the device is abnormal. Conduct statistical analysis on the beam current frequency-domain and radiation frequency-domain data under abnormal device states. For example, calculate the mean and standard deviation statistics for each frequency component. For the beam current frequency domain, according to the results of the statistical analysis, determine the amplitude value of the frequency component that exceeds the mean plus a certain multiple of the standard deviation (such as 2 or 3 times the standard deviation) as the second target critical value of the beam current frequency domain. For example, if the amplitude of a certain frequency component often exceeds the mean plus 2 times the standard deviation when the device is abnormal, then this value is used as the second target critical value of the beam current frequency domain at that frequency. For the radiation frequency domain, use a similar method to determine the second target critical value of the radiation frequency domain based on the statistical analysis of the radiation frequency-domain data. For example, use the radiation intensity value that exceeds the mean plus 3 times the standard deviation as the second target critical value corresponding to that frequency in the radiation frequency domain.

[0073] A possible implementation manner of the embodiment of the present application. In the above embodiment, determining the degree of association between the device load and the device stability includes: Sort the historical device load and the historical device stability respectively according to the historical start time corresponding to each historical period to obtain a load sequence and a stability sequence; Based on the load sequence and the stability sequence, establish a first regression equation of the device stability with respect to its own lag term, and a second regression equation of the device stability with respect to its own lag term and the lag term of the device load; Based on the first regression equation and the second regression equation, determine the causal relationship direction between the device load and the device stability; Calculate the correlation coefficient between the device load and the device stability, and combine the causal relationship direction and the correlation coefficient to determine the degree of association between the device load and the device stability.

[0074] Extract the historical equipment load and historical equipment stability data corresponding to each historical period from the stored historical data, and at the same time obtain the historical start time of each historical period. Then, sort the historical equipment load data according to the chronological order of the historical start times to form a load sequence arranged in increasing order of time. Similarly, sort the historical equipment stability data according to the historical start times to obtain a stability sequence. For example, if there are three historical periods with start times t1, t2, and t3 (t1 < t2 < t3) respectively, the equipment load and stability data corresponding to time t1 will be ranked first, those corresponding to t2 will be ranked second, and those corresponding to t3 will be ranked last.

[0075] Furthermore, select multiple-order lags, comprehensively consider the data of multiple historical periods, and establish the first regression equation and the second regression equation. Specifically, based on the stability sequence, construct a regression equation for the equipment stability with respect to its own lag terms. Assume that the equipment stability is represented by , the first-order lag term is , then the first regression equation can be expressed as , where is the constant term, , are the regression coefficients, n is the number of historical periods, is the error term. Among them, methods such as the least squares method can be used to estimate the coefficients , , in the equation using the data of the load sequence and the stability sequence.

[0076] Based on the first regression equation, add the lag term of the equipment load. The equipment load is represented by , the first-order lag term is , then the second regression equation is , where is the constant term, , and are the regression coefficients, is the error term. Similarly, use methods such as the least squares method to estimate the coefficients , , and in the equation.

[0077] Propose the null hypothesis : The equipment load is not the Granger cause of the equipment stability, that is, the coefficients of the equipment load lag terms in the second regression equation are all 0 (such as = 0) in the above equation. By calculating the F statistic, compare the residual sum of squares of the two regression equations. The calculation formula of the F statistic is , where is the sum of squared residuals of the first regression equation, is the sum of squared residuals of the second regression equation, q is the number of lag terms of the equipment load, n is the number of historical periods, and p is the number of lag terms of the equipment stability itself.

[0078] Furthermore, according to the calculated F statistic and the preset significance level (such as 0.05), look up the critical value in the F distribution table. If the F statistic is greater than the critical value and the corresponding p value is less than the significance level, reject the original hypothesis and determine that the equipment load is the Granger cause of the equipment stability, that is, the change in the equipment load may cause a change in the equipment stability; otherwise, determine that the equipment load is not the Granger cause of the equipment stability. At the same time, the roles of the equipment load and the equipment stability can be exchanged, and the above process can be repeated to determine whether the equipment stability is the Granger cause of the equipment load, so as to determine the direction of the causal relationship between the two.

[0079] Even further, the Pearson correlation coefficient can be selected to calculate the correlation between the equipment load and the equipment stability. Specifically, calculate the mean and standard deviation of the equipment load sequence and the stability sequence, and then calculate according to the Pearson correlation coefficient formula where and are the i-th values of the equipment load sequence and the stability sequence respectively, and are the means of the two sequences respectively, and n is the number of historical periods.

[0080] Combine the calculated correlation coefficient and the previously determined causal relationship direction to comprehensively judge the degree of association between the equipment load and the equipment stability. Specifically, when the causal relationship is clear, determine the coefficient corresponding to the causal relationship as 1, when the causal relationship is not clear, determine the coefficient corresponding to the causal relationship as 0, and calculate the degree of association based on the fifth formula and the relevant washing and causal relationship, where the fifth formula is , A is the degree of association, r is the correlation coefficient, and c is the coefficient corresponding to the causal relationship.

[0081] In a possible implementation manner of the embodiment of the present application, in the above step S104, based on the critical electronic information and the current electronic information, determining whether the working state of the electron accelerator is adjusted from the normal state to the trigger state includes: Compare the current radiation dose with the critical radiation dose, and compare the current beam parameters with the critical beam parameters to determine whether the current radiation dose exceeds the critical radiation dose and whether the current beam parameters exceed the critical beam parameters; If the current radiation dose exceeds the critical radiation dose, and / or the current beam parameters exceed the critical beam parameters, it is determined that the operating state of the electron accelerator is adjusted from the normal state to the trigger state; If the current radiation dose does not exceed the critical radiation dose and the current beam parameters do not exceed the critical beam parameters, it is determined that the operating state of the electron accelerator remains the normal state.

[0082] Obtain the current radiation dose data in real time from the radiation dose monitor installed around the electron accelerator, and at the same time obtain the current beam parameters (such as beam intensity, particle passing speed, number of particles passing through, etc.) from the beam monitoring device. In addition, extract the critical radiation dose and critical beam parameters from the pre-calculated and stored critical electron information.

[0083] Further, numerically compare the current radiation dose with the critical radiation dose. There may be upper and lower limit values for the critical radiation dose (for example, the radiation dose should be within a specific range during normal operation), and check whether the current radiation dose exceeds this range. Similarly, for each item in the current beam parameters, such as beam intensity, particle passing speed, etc., compare it with the corresponding critical beam parameter respectively to determine whether it exceeds the critical range.

[0084] When the current radiation dose exceeds the range of the critical radiation dose, or any one of the current beam parameters exceeds the corresponding critical beam parameter range, or both situations occur simultaneously, mark that the operating state of the electron accelerator has changed, that is, update the operating state from the previously recorded normal state to the trigger state, and make corresponding state records and identifications in the internal system.

[0085] If it is found that the current radiation dose is within the normal range specified by the critical radiation dose, and each item of the current beam parameters is also within the range allowed by the corresponding critical beam parameters, then the operating state of the electron accelerator will be maintained as the normal state. At the same time, continue to monitor the current radiation dose and current beam parameters in real time to detect possible abnormal situations in a timely manner.

[0086] See Figure 2 ., the high-energy electron accelerator safety interlock device 20 may specifically include: an acquisition module 201, a calculation module 202, an analysis module 203, a determination module 204, and a generation module 205, where: A high-energy electron accelerator safety interlock device 20, including: An acquisition module 201, configured to acquire historical electronic information and historical equipment load corresponding to an electron accelerator, and to monitor in real time the current electronic information and current equipment information of the electron accelerator. The historical electronic information includes electronic information corresponding to a historical period, and the historical period is a time period between when the working state of the electron accelerator changes from a normal state to a trigger state. The electronic information includes radiation dose and beam parameters; A calculation module 202, configured to calculate the historical equipment stability degree corresponding to the electron accelerator; An analysis module 203, configured to analyze the historical electronic information, historical equipment load, and historical equipment stability degree to determine the critical electronic information corresponding to the electron accelerator; A determination module 204, configured to determine whether the working state of the electron accelerator changes from a normal state to a trigger state based on the critical electronic information and the current electronic information; A generation module 205, configured to generate a safety trigger instruction if it is determined that the working state of the electron accelerator changes from a normal state to a trigger state.

[0087] In a possible implementation manner of the embodiment of the present application, when the calculation module 202 calculates the historical equipment stability degree corresponding to the electron accelerator and in the calculation process of the historical equipment stability degree corresponding to each historical period, it may specifically be used for: Obtain each historical position and historical preset orbit of the electron accelerator, and determine whether each historical position is located on the historical preset orbit; if there is a historical position not on the historical preset orbit, then determine the historical position not on the historical preset orbit as the target historical position, and calculate the historical deviation between each target historical position and the historical preset orbit, obtain the deviation allowable range and historical confidence level corresponding to the electron accelerator, calculate the confidence interval corresponding to the average value of the historical deviation based on the historical confidence level, and calculate the historical equipment stability degree corresponding to the electron accelerator based on the formula, where the formula is , S is the historical equipment stability degree, is the maximum value of the confidence interval corresponding to the average value of the historical deviation, is the minimum value of the deviation allowable range, is the maximum value of the deviation allowable range; Obtain the historical beam intensity and historical beam energy corresponding to each historical moment within the historical period of the electron accelerator, and determine the beam intensity fluctuation situation and beam energy fluctuation situation, respectively perform frequency domain analysis on the beam intensity fluctuation situation and beam energy fluctuation situation to establish a beam intensity frequency domain relationship and a beam energy fluctuation frequency relationship, and establish a beam intensity amplitude matrix and a beam energy amplitude matrix, and determine the historical equipment stability degree corresponding to the electron accelerator based on the beam intensity frequency domain relationship, beam energy fluctuation frequency relationship, beam intensity amplitude matrix, and beam energy amplitude matrix.

[0088] In a possible implementation manner of the embodiment of the present application, when the analysis module 203 analyzes historical electronic information, historical equipment load, and historical equipment stability to determine the critical electronic information corresponding to the electron accelerator, it can specifically be used for: Obtain the equipment load threshold corresponding to the electron accelerator, and based on the equipment load threshold, screen out at least one first historical period from multiple historical periods, and determine the historical periods other than the first historical period as the second historical period. The historical equipment load corresponding to the first historical period does not exceed the equipment load threshold, and the historical equipment stability corresponding to the first historical period is not less than the equipment stability threshold; Based on the beam parameters corresponding to each first historical period, establish the fault time domain relationship and beam time domain relationship corresponding to the electron accelerator; Screen out the historical moments when the historical equipment load of the electron accelerator exceeds the equipment load threshold from each second historical period as the second historical moments, and obtain the second historical electronic information corresponding to each second historical moment from the historical electronic information; Perform frequency domain analysis on the beam parameters and radiation doses corresponding to the second historical electronic information respectively to obtain the beam frequency domain relationship and radiation frequency domain relationship corresponding to the electron accelerator; Based on the fault time domain relationship, beam time domain relationship, beam frequency domain relationship, and radiation frequency domain relationship corresponding to the electron accelerator, determine the critical electronic information corresponding to the electron accelerator.

[0089] In a possible implementation manner of the embodiment of the present application, the beam parameters include the particle passing speed and the number of particles passing through. When the analysis module 203 establishes the fault time domain relationship and beam time domain relationship corresponding to the electron accelerator based on the beam parameters corresponding to each first historical period, it can specifically be used for: Determine the historical duration corresponding to each first historical period, and establish the fault time domain relationship corresponding to the electron accelerator based on each historical duration; Based on the historical electronic information, obtain the number of particles passing through and the particle passing speed corresponding to each first historical period, and establish the beam time domain relationship corresponding to the electron accelerator based on the number of particles passing through and the particle passing speed corresponding to each first historical period.

[0090] In a possible implementation manner of the embodiment of the present application, when the analysis module 203 determines the critical electronic information corresponding to the electron accelerator based on the fault time domain relationship, beam time domain relationship, beam frequency domain relationship, and radiation frequency domain relationship corresponding to the electron accelerator, it can specifically be used for: Based on the fault time domain relationship and beam time domain relationship corresponding to the electron accelerator, establish a decision tree model; Based on the decision tree model, predict the abnormal state of the electron accelerator in the current cycle. The abnormal state of the equipment is abnormal or non-abnormal, and the abnormal state includes equipment load abnormality and / or equipment stability degree abnormality. When the abnormal state of the electron accelerator in the current cycle is non-abnormal, determine the first critical value based on the beam current frequency domain relationship and the radiation frequency domain relationship, and determine the second critical value based on the fault time domain relationship and the beam current time domain relationship. Determine the critical electron information corresponding to the electron accelerator based on the first critical value and the second critical value. When the abnormal state of the electron accelerator in the current cycle is abnormal, determine the correlation degree between the equipment load and the equipment stability degree. Determine the first target critical value based on the correlation degree, the fault time domain relationship, and the beam current time domain relationship, and determine the second target critical value based on the beam current frequency domain relationship and the radiation frequency domain relationship.

[0091] In a possible implementation manner of the embodiments of the present application, when the analysis module 203 determines the correlation degree between the equipment load and the equipment stability degree, it can specifically be used for: Sort the historical equipment load and the historical equipment stability degree according to the historical start time corresponding to each historical cycle to obtain a load sequence and a stability sequence. Based on the load sequence and the stability sequence, establish a first regression equation of the equipment stability degree with respect to its own lag term, and a second regression equation of the equipment stability degree with respect to its own lag term and the lag term of the equipment load. Based on the first regression equation and the second regression equation, determine the causal relationship direction between the equipment load and the equipment stability degree. Calculate the correlation coefficient between the equipment load and the equipment stability degree, and determine the correlation degree between the equipment load and the equipment stability degree in combination with the causal relationship direction and the correlation coefficient.

[0092] In a possible implementation manner of the embodiments of the present application, when the determination module 204 determines whether the working state of the electron accelerator is adjusted from the normal state to the trigger state based on the critical electron information and the current electron information, it can specifically be used for: Compare the current radiation dose with the critical radiation dose, and compare the current beam current parameters with the critical beam current parameters to determine whether the current radiation dose exceeds the critical radiation dose and whether the current beam current parameters exceed the critical beam current parameters. If the current radiation dose exceeds the critical radiation dose, or / and, the current beam current parameters exceed the critical beam current parameters, then determine that the working state of the electron accelerator is adjusted from the normal state to the trigger state. If the current radiation dose does not exceed the critical radiation dose and the current beam current parameters do not exceed the critical beam current parameters, then determine that the working state of the electron accelerator remains the normal state.

[0093] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the systems, devices, and units described above can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0094] Participate Figure 3 In addition, embodiments of the present application also introduce an electronic device from the perspective of an entity device. For example, Figure 3 as shown, Figure 3 The electronic device 300 shown in FIG. includes: a processor 301 and a memory 303. Among them, the processor 301 and the memory 303 are connected, such as through a bus 302. Optionally, the electronic device 300 may further include a transceiver 304. It should be noted that in practical applications, the transceiver 304 is not limited to one, and the structure of the electronic device 300 does not constitute a limitation to the embodiments of the present application.

[0095] The processor 301 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It can implement or execute various exemplary logical blocks, modules, and circuits described in connection with the disclosure of the present application. The processor 301 may also be a combination for implementing computing functions, such as a combination including one or more microprocessors, a combination of a DSP and a microprocessor, etc.

[0096] The bus 302 may include a path for transmitting information between the above components. The bus 302 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. The bus 302 may be divided into an address bus, a data bus, a control bus, etc. For the convenience of representation, Figure 3 only a thick line is used to represent it in the figure, but it does not mean that there is only one bus or one type of bus.

[0097] The memory 303 can be a ROM (Read Only Memory), or other types of static storage devices that can store static information and instructions, a RAM (Random Access Memory), or other types of dynamic storage devices that can store information and instructions. It can also be an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read Only Memory), or other optical disc storage, optical disc storage (including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), magnetic disk storage media, or other magnetic storage devices, or any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto.

[0098] The memory 303 is used to store the application program code for executing the solution of this application, and is controlled by the processor 301 for execution. The processor 301 is used to execute the application program code stored in the memory 303 to implement the content shown in the foregoing method embodiments.

[0099] Among them, the electronic device includes but is not limited to: mobile terminals such as mobile phones, laptop computers, digital broadcast receivers, PDAs (Personal Digital Assistants), PADs (Tablet Computers), PMPs (Portable Multimedia Players), vehicle terminals (such as vehicle navigation terminals), etc., and fixed terminals such as digital TVs, desktop computers, etc., and can also be a server, etc. Figure 3 The illustrated electronic device is only an example and should not impose any limitations on the functions and usage scope of the embodiments of this application.

[0100] The embodiments of this application provide a computer-readable storage medium, on which a computer program is stored. When it runs on a computer, it enables the computer to execute the corresponding content in the foregoing method embodiments.

[0101] It should be understood that although the steps in the flowchart of the accompanying drawings are shown in sequence according to the indication of the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless there is a clear description in this article, the execution of these steps has no strict order limit, and they can be executed in other orders. Moreover, at least a part of the steps in the flowchart of the accompanying drawings may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times, and their execution order is not necessarily sequential, but can be executed alternately or in turn with at least a part of other steps or sub-steps or stages of other steps.

[0102] The above are only some embodiments of the present application. It should be noted that for those of ordinary skill in the art, without departing from the principle of the present application, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present application.

Claims

1. A safety interlock method for a high-energy electron accelerator, characterized in that, Including: Obtaining historical electron information and historical equipment load corresponding to an electron accelerator, and real-time monitoring current electron information and current equipment information of the electron accelerator. The historical electron information includes electron information corresponding to a historical period, and the historical period is a time period between when the working state of the electron accelerator changes from a normal state to a trigger state. The electron information includes radiation dose and beam parameters; Calculating the historical equipment stability degree corresponding to the electron accelerator; Analyzing the historical electron information, the historical equipment load, and the historical equipment stability degree to determine critical electron information corresponding to the electron accelerator; Based on the critical electron information and the current electron information, determining whether the working state of the electron accelerator changes from a normal state to a trigger state; If it is determined that the working state of the electron accelerator changes from a normal state to a trigger state, generating a safety trigger instruction.

2. The high-energy electron accelerator safety interlock method according to claim 1, wherein For the calculation of the historical equipment stability degree corresponding to the electron accelerator, the calculation process of the historical equipment stability degree corresponding to each historical period includes any one of the following: Obtain each historical position and historical preset orbit of the electron accelerator, and determine whether each historical position is on the historical preset orbit; if there is a historical position not on the historical preset orbit, determine the historical position not on the historical preset orbit as the target historical position, and calculate the historical deviation between each target historical position and the historical preset orbit, obtain the deviation allowable range and historical confidence level corresponding to the electron accelerator, calculate the confidence interval corresponding to the average value of the historical deviation based on the historical confidence level, and calculate the corresponding historical equipment stability degree of the electron accelerator based on the formula, where the formula is , S is the historical equipment stability degree, is the maximum value of the confidence interval corresponding to the average value of the historical deviation, is the minimum value of the deviation allowable range, is the maximum value of the deviation allowable range; Obtaining historical beam intensity and historical beam energy corresponding to each historical moment within the historical period of the electron accelerator, determining the beam intensity fluctuation situation and the beam energy fluctuation situation, respectively performing frequency domain analysis on the beam intensity fluctuation situation and the beam energy fluctuation situation to establish a beam intensity frequency domain relationship and a beam energy fluctuation frequency relationship, and establishing a beam intensity amplitude matrix and a beam energy amplitude matrix. Based on the beam intensity frequency domain relationship, the beam energy fluctuation frequency relationship, the beam intensity amplitude matrix, and the beam energy amplitude matrix, determining the historical equipment stability degree corresponding to the electron accelerator.

3. The high-energy electron accelerator safety interlock method according to claim 2, characterized in that, Analyzing the historical electron information, the historical equipment load, and the historical equipment stability degree to determine critical electron information corresponding to the electron accelerator, including: Obtaining an equipment load threshold corresponding to the electron accelerator, and based on the equipment load threshold, screening out at least one first historical period from multiple historical periods, and determining historical periods other than the first historical period as second historical periods. The historical equipment load corresponding to the first historical period does not exceed the equipment load threshold, and the historical equipment stability degree corresponding to the first historical period is not less than the equipment stability degree threshold; Based on the beam parameters corresponding to each first historical period, establishing a fault time domain relationship and a beam time domain relationship corresponding to the electron accelerator; Screening out historical moments when the historical equipment load of the electron accelerator in each second historical period exceeds the equipment load threshold as second historical moments, and obtaining second historical electron information corresponding to each second historical moment from the historical electron information; Respectively performing frequency domain analysis on the beam parameters and radiation dose corresponding to the second historical electron information to obtain a beam frequency domain relationship and a radiation frequency domain relationship corresponding to the electron accelerator; Based on the fault time domain relationship, the beam time domain relationship, the beam frequency domain relationship, and the radiation frequency domain relationship corresponding to the electron accelerator, determining the critical electron information corresponding to the electron accelerator.

4. The high-energy electron accelerator safety interlock method according to claim 3, wherein The beam parameters include the particle passing velocity and the particle passing quantity. Establishing the fault time domain relationship and the beam time domain relationship corresponding to the electron accelerator based on the beam parameters corresponding to each of the first historical periods includes: Determining the historical duration corresponding to each of the first historical periods, and establishing the fault time domain relationship corresponding to the electron accelerator based on each historical duration; Based on the historical electron information, obtaining the particle passing quantity and the particle passing velocity corresponding to each of the first historical periods, and establishing the beam time domain relationship corresponding to the electron accelerator based on the particle passing quantity and the particle passing velocity corresponding to each of the first historical periods.

5. The high-energy electron accelerator safety interlock method according to claim 3 or 4, characterized in that, Determining the critical electron information corresponding to the electron accelerator based on the fault time domain relationship, the beam time domain relationship, the beam frequency domain relationship, and the radiation frequency domain relationship corresponding to the electron accelerator includes: Establishing a decision tree model based on the fault time domain relationship and the beam time domain relationship corresponding to the electron accelerator; Based on the decision tree model, predicting the device abnormal state of the electron accelerator in the current period, where the device abnormal state is abnormal or non-abnormal, and the abnormal state is device load abnormality and / or device stability degree abnormality; When the device abnormal state of the electron accelerator in the current period is non-abnormal, determining a first critical value based on the beam frequency domain relationship and the radiation frequency domain relationship, and determining a second critical value based on the fault time domain relationship and the beam time domain relationship, and determining the critical electron information corresponding to the electron accelerator based on the first critical value and the second critical value; When the device abnormal state of the electron accelerator in the current period is abnormal, determining the correlation degree between the device load and the device stability degree, determining a first target critical value based on the correlation degree, the fault time domain relationship, and the beam time domain relationship, and determining a second target critical value based on the beam frequency domain relationship and the radiation frequency domain relationship.

6. The high-energy electron accelerator safety interlock method according to claim 5, characterized in that Determining the correlation degree between the device load and the device stability degree includes: Sorting the historical device load and the historical device stability degree respectively according to the historical start time corresponding to each historical period to obtain a load sequence and a stability sequence; Based on the load sequence and the stability sequence, establishing a first regression equation of the device stability degree with respect to its own lag term, and a second regression equation of the device stability degree with respect to its own lag term and the device load lag term; Based on the first regression equation and the second regression equation, determining the causal relationship direction between the device load and the device stability degree; Calculating the correlation coefficient between the device load and the device stability degree, and determining the correlation degree between the device load and the device stability degree in combination with the causal relationship direction and the correlation coefficient.

7. The high-energy electron accelerator safety interlock method according to claim 1, characterized in that, Determining whether the working state of the electron accelerator is adjusted from the normal state to the trigger state based on the critical electron information and the current electron information includes: Compare the current radiation dose with the critical radiation dose, and compare the current beam parameters with the critical beam parameters to determine whether the current radiation dose exceeds the critical radiation dose and whether the current beam parameters exceed the critical beam parameters; If the current radiation dose exceeds the critical radiation dose, or / and, the current beam parameters exceed the critical beam parameters, then determine that the working state of the electron accelerator is adjusted from the normal state to the trigger state; If the current radiation dose does not exceed the critical radiation dose and the current beam parameters do not exceed the critical beam parameters, then determine that the working state of the electron accelerator remains the normal state.

8. A safety interlock device for a high-energy electron accelerator, characterized in that, Comprising: An acquisition module, configured to acquire the historical electronic information and historical equipment load corresponding to the electron accelerator, and to monitor in real time the current electronic information and current equipment information of the electron accelerator. The historical electronic information includes the electronic information corresponding to the historical period, and the historical period is the time period between the working state of the electron accelerator changing from the normal state to the trigger state. The electronic information includes the radiation dose and the beam parameters; A calculation module, configured to calculate the historical equipment stability degree corresponding to the electron accelerator; An analysis module, configured to analyze the historical electronic information, the historical equipment load, and the historical equipment stability degree to determine the critical electronic information corresponding to the electron accelerator; A determination module, configured to determine whether the working state of the electron accelerator is adjusted from the normal state to the trigger state based on the critical electronic information and the current electronic information; A generation module, configured to generate a safety trigger instruction if it is determined that the working state of the electron accelerator is adjusted from the normal state to the trigger state.

9. An electronic device, characterized in that, This electronic device includes: At least one processor; A memory; At least one application program, where at least one application program is stored in the memory and is configured to be executed by at least one processor. The at least one application program is configured to: execute the high-energy electron accelerator safety interlock method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed on the computer, cause the computer to execute the high-energy electron accelerator safety interlock method according to any one of claims 1 to 7.

Images