Electronic cyclotron radiation meter calibration method and system based on plasma displacement control

Abstract

The application discloses a method and system for calibrating an electron cyclotron radiation gauge based on plasma displacement control, and relates to the technical field of electron cyclotron radiation gauge calibration. The method comprises the following steps: obtaining a recursive relationship between adjacent channels based on the condition that the measurement temperature of the same plasma in each measurement space position is consistent; obtaining the calibration coefficients of each channel in the electron cyclotron radiation gauge by using the time difference of adjacent channels, the recursive relationship and the preset value of the first channel in the electron cyclotron radiation gauge; and calibrating each channel in the electron cyclotron radiation gauge by using the calibration coefficients of each channel. The method and system can realize the measurement space overlap between adjacent channels or even multiple channels, and can accurately calibrate the electron cyclotron radiation gauge based on accurate plasma control and mature displacement measurement means, thereby greatly improving the calibration accuracy.

Description

Technical Field

[0001] The present invention relates to the technical field of electron cyclotron radiometer calibration, and more particularly to an electron cyclotron radiometer calibration method and system based on plasma displacement control. Background Art

[0002] Electron cyclotron radiometers (ECRs) are an important microwave diagnostic tool for electron temperature profiles and temperature disturbances in tokamak plasma physics research. Due to their high temporal and spatial resolution, as well as their ability to measure small, localized regions, ECRs are crucial diagnostic tools in tokamak plasma physics. At tokamaks worldwide, ECRs measure electron temperature profiles and temperature disturbances, enabling a wide range of studies, including magnetohydrodynamic instabilities, microscopic turbulent transport, rupture prediction, identification and localization of neoclassical tearing modes, electron thermal transport, and high-energy particle instabilities.

[0003] A crucial prerequisite for using diagnostic data from electron cyclotron radiometers for research is accurate system calibration and the acquisition of calibration coefficients. Two common calibration methods exist: one is the dual-temperature method, which is complex and makes it difficult to accurately measure the emissivity of a blackbody source. Furthermore, since microwave radiation in plasma is simulated through mirror reflection, there are numerous uncertainties, resulting in considerable errors in the calibration results. The other is the magnetic field difference method, which varies the strength of the toroidal magnetic field between the two discharge guns to adjust the relative calibration coefficient between the front and rear channels. This method is simple, but the measurement intermediate frequency interval of traditional multi-channel electron cyclotron radiometers is constant. In a tokamak, the toroidal magnetic field decreases inversely along the radial direction. Therefore, when using the magnetic field difference method to calibrate a traditional multi-channel electron cyclotron radiometer, the measurement positions of the front and rear channels cannot completely overlap between the calibration guns, resulting in errors in the final calibration coefficients. This method cannot effectively meet the requirements for accurate electron temperature and disturbance distributions in various studies, impacting the accuracy of some important research results. Summary of the Invention

[0004] The purpose of the present invention is to provide a calibration method and system for an electron cyclotron radiometer based on plasma displacement control, which solves the problems in traditional calibration methods such as the complexity of the measurement process, the uncertainty of the blackbody source radiation coefficient, non-ideal in-situ calibration, and difficulty in error analysis; and the magnetic field difference method is limited by the fact that the positions of adjacent measurement channels between calibration guns cannot be completely overlapped. Although the calibration accuracy can be improved through mathematical smoothing, there are considerable errors due to factors such as plasma displacement and discharge instability. This method and system can achieve measurement space overlap between adjacent channels or even multiple channels. Based on accurate plasma control and mature displacement measurement means, the electron cyclotron radiometer can be accurately calibrated, greatly improving the calibration accuracy.

[0005] The above technical objectives of the present invention are achieved through the following technical solutions:

[0006] In a first aspect, the present application provides a calibration method for an electron cyclotron radiometer based on plasma displacement control, comprising the following specific steps:

[0007] Based on the toroidal magnetic field distribution in the tokamak device, the measurement space position corresponding to each channel of the electron cyclotron radiometer under the toroidal magnetic field distribution is calculated, and the measurement position interval between different channels is calculated using the measurement space position of each channel;

[0008] Establishing a first functional relationship between plasma displacement and time based on a preset waveform for controlling plasma displacement in the tokamak device;

[0009] Using the measurement position intervals and the first functional relationship between different channels in the electron cyclotron radiometer, the time difference when adjacent channels measure the same plasma position under the condition of plasma horizontal displacement fluctuation is calculated;

[0010] Based on the condition that the measured temperature of the same plasma in each measurement space position is consistent, a recursive relationship between adjacent channels is obtained, and the calibration coefficients of each channel in the electron cyclotron radiometer are obtained through the time difference and recursive relationship of adjacent channels, as well as the preset value of the first channel in the electron cyclotron radiometer. Each channel in the electron cyclotron radiometer is calibrated using the calibration coefficients of each channel.

[0011] On the basis of the above technical solution, the present invention can also be improved as follows.

[0012] Furthermore, the measurement position intervals between the above different channels are specifically as follows:

[0013]

[0014] Where, I Nrepresents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer, f N+1 Indicates the intermediate frequency of the N+1th channel of the electron cyclotron radiometer.

[0015] Furthermore, the above time difference is obtained by:

[0016] I N =D(t N +t0+n×T0), where D=f(t);

[0017] Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, D represents the first functional relationship, t N It represents the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, n represents a positive integer, and T0 represents the period of horizontal displacement of the plasma.

[0018] Furthermore, the above recursive relationship is specifically:

[0019]

[0020] Where E1 represents the radiation signal intensity of the toroidal magnetic field distribution received by the first channel of the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, t1 represents the time difference between the first and second channels of the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, n represents a positive integer, E2 represents the radiation signal intensity of the toroidal magnetic field distribution received by the second channel of the electron cyclotron radiometer, C1 represents the calibration coefficient of the first channel of the electron cyclotron radiometer, and C2 represents the calibration coefficient of the second channel of the electron cyclotron radiometer.

[0021] Furthermore, the calibration coefficients of the above channels are specifically:

[0022]

[0023] Where, E N It represents the radiation signal intensity of the toroidal magnetic field distribution received by the Nth channel of the electron cyclotron radiometer, t N Indicates the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, C NIndicates the calibration coefficient of the Nth channel in the electron cyclotron radiometer, C N+1 Indicates the calibration coefficient of the N+1th channel in the electron cyclotron radiometer, E N+1 It represents the radiation signal intensity of the toroidal magnetic field distribution received by the N+1th channel of the electron cyclotron radiometer, t N+1 It represents the time difference between the N+1th channel and the N+2th channel in the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, and n represents a positive integer.

[0024] Furthermore, the above-mentioned toroidal magnetic field distribution is specifically as follows:

[0025]

[0026] Where B represents the toroidal magnetic field distribution, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, and R represents the measurement position.

[0027] Furthermore, the above measurement positions are specifically:

[0028]

[0029] Where R N represents the measurement position of the Nth channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer.

[0030] In a second aspect, the present application provides an electron cyclotron radiometer calibration system based on plasma displacement control, which is applied to any electron cyclotron radiometer calibration method based on plasma displacement control in the first aspect, comprising:

[0031] The first module is used to calculate the measurement space position corresponding to each channel of the electron cyclotron radiometer under the toroidal magnetic field distribution in the tokamak device, and use the measurement space position of each channel to calculate the measurement position interval between different channels;

[0032] The second module is configured to establish a first functional relationship between plasma displacement and time based on a preset waveform for controlling plasma displacement in the tokamak device;

[0033] The third module is used to calculate the time difference when adjacent channels measure the same plasma position under the horizontal displacement fluctuation of each channel by using the measurement position interval between different channels in the electron cyclotron radiometer and the first functional relationship;

[0034] The fourth module is used to obtain the recursive relationship between adjacent channels based on the condition that the measured temperature of the same plasma in each measurement space position is consistent, and to obtain the calibration coefficients of each channel in the electron cyclotron radiometer through the time difference and recursive relationship between adjacent channels and the preset value of the first channel in the electron cyclotron radiometer.

[0035] In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements any one of the methods in the first aspect when executing the computer program.

[0036] In a fourth aspect, the present application provides a non-transitory computer-readable storage medium, which stores computer instructions, and the computer instructions enable a computer to execute any one of the methods in the first aspect.

[0037] Compared with the prior art, the present invention has at least the following beneficial effects:

[0038] In this application, based on accurate plasma displacement control technology, by calculating the relationship between the overlapping time of the displacement plasma space position between adjacent channels, and comparing the radiation signal intensity of the adjacent channels measuring the overlapping moment in multiple cycles, the calibration coefficients of all channels are recursively obtained; since the circumferential magnetic field intensity of each tokamak device in the low-pressure plasma is approximately an ideal inverse proportional function, the measurement position between each channel of the electron cyclotron radiometer can be accurately calculated based on the measurement frequency. Then the waveform and period of the plasma horizontal displacement are pre-designed. Since the electron temperature distribution is an ideal magnetic surface function, the time for adjacent channels to measure the same plasma space position can be calculated by establishing a functional relationship between the plasma horizontal displacement and the measurement interval between adjacent channels. Therefore, according to the radiation intensity timing signals received by different channels, the calibration coefficients of each channel can be compared and recursively deduced, which greatly improves the calibration accuracy of the electron cyclotron radiometer.

[0039] In the present application, based on the plasma displacement control method, adjacent channels are used to measure the same plasma space at different times simultaneously, and finally the calibration coefficients of all channels of the electron cyclotron radiometer are obtained through the recursive relationship; wherein, the key to realize the calibration of the electron cyclotron radiometer is plasma displacement control. This method can not only use adjacent single channels to measure the same plasma space, but also use adjacent multiple channels to measure the same plasma space simultaneously; in addition to being used in the calibration of the electron cyclotron radiometer of the tokamak device, this method can also be applied to any magnetic confinement device in which the electron cyclotron radiometer measures the local electron temperature distribution, such as a stellarator. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, constitute a part of this application, and do not constitute a limitation of the embodiments of the present invention. In the drawings:

[0041] Figure 1 A flow chart of a calibration method according to an embodiment of the present invention;

[0042] Figure 2 Schematic diagram of the measurement position interval between the Nth channel and the N+1th channel in an embodiment of the present invention;

[0043] Figure 3 Schematic diagram of the time difference between adjacent channels in an embodiment of the present invention;

[0044] Figure 4 Schematic diagram of the connection of the calibration system in an embodiment of the present invention. DETAILED DESCRIPTION

[0045] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings herein can be arranged and designed in various different configurations.

[0046] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention as claimed, but rather merely represents selected embodiments of the present invention. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without creative effort shall fall within the scope of protection of the present invention.

[0047] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, it does not need to be further defined or explained in subsequent drawings.

[0048] Example 1: In order to solve the problems in the traditional calibration method, that the absolute in-situ calibration of the blackbody source is limited by the complexity of the measurement process, the uncertainty of the blackbody source radiation coefficient, the non-ideal in-situ calibration, the difficulty of error analysis, and the magnetic field difference method is limited by the position of the adjacent measurement channels between the calibration guns cannot be completely overlapped, although the calibration accuracy can be improved through mathematical smoothing, there are considerable errors due to factors such as plasma displacement and discharge instability. This embodiment provides an electron cyclotron radiometer calibration method based on plasma displacement control, which can achieve measurement space overlap between adjacent channels or even multiple channels. Based on accurate plasma control and mature displacement measurement methods, the electron cyclotron radiometer can be accurately calibrated, greatly improving the calibration accuracy, such as Figure 1 As shown, the method includes the following specific steps:

[0049] S1. Based on the toroidal magnetic field distribution in the tokamak device, the measurement space position corresponding to each channel of the electron cyclotron radiometer under the toroidal magnetic field distribution is calculated, and the measurement position interval between different channels is calculated using the measurement space position of each channel.

[0050] Optionally, the measurement position intervals between the above different channels are specifically:

[0051]

[0052] Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer, f N+1 Indicates the intermediate frequency of the N+1th channel of the electron cyclotron radiometer.

[0053] Optionally, the above-mentioned toroidal magnetic field distribution is specifically as follows:

[0054]

[0055] Where B represents the toroidal magnetic field distribution, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, and R represents the measurement position.

[0056] Optionally, the above measurement positions are specifically:

[0057]

[0058] Where R N represents the measurement position of the Nth channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer.

[0059] S2, establishing a first functional relationship between plasma displacement and time based on a preset waveform for controlling plasma displacement in the tokamak device.

[0060] S3, using the measurement position intervals between different channels in the electron cyclotron radiometer and the first functional relationship, calculate the time difference when adjacent channels measure the same plasma position under the horizontal displacement fluctuation of the plasma.

[0061] Optionally, the time difference is obtained by:

[0062] I N =D(t N +t0+n×T0), where D=f(t);

[0063] Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, D represents the first functional relationship, t N It represents the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, n represents a positive integer, and T0 represents the period of horizontal displacement of the plasma.

[0064] S4. Based on the condition that the measured temperature of the same plasma in each measurement space position is consistent, a recursive relationship between adjacent channels is obtained, and the calibration coefficients of each channel in the electron cyclotron radiometer are obtained through the time difference and recursive relationship of adjacent channels, as well as the preset value of the first channel in the electron cyclotron radiometer, and each channel in the electron cyclotron radiometer is calibrated using the calibration coefficients of each channel.

[0065] Optionally, the above recursive relationship is specifically:

[0066]

[0067] Where E1 represents the radiation signal intensity of the toroidal magnetic field distribution received by the first channel of the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, t1 represents the time difference between the first and second channels of the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, n represents a positive integer, E2 represents the radiation signal intensity of the toroidal magnetic field distribution received by the second channel of the electron cyclotron radiometer, C1 represents the calibration coefficient of the first channel of the electron cyclotron radiometer, and C2 represents the calibration coefficient of the second channel of the electron cyclotron radiometer.

[0068] Optionally, the calibration coefficients of the above channels are specifically:

[0069]

[0070] Where, E N It represents the radiation signal intensity of the toroidal magnetic field distribution received by the Nth channel of the electron cyclotron radiometer, tN Indicates the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, C N Indicates the calibration coefficient of the Nth channel in the electron cyclotron radiometer, C N+1 Indicates the calibration coefficient of the N+1th channel in the electron cyclotron radiometer, E N+1 It represents the radiation signal intensity of the toroidal magnetic field distribution received by the N+1th channel of the electron cyclotron radiometer, t N+1 It represents the time difference between the N+1th channel and the N+2th channel in the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, and n represents a positive integer.

[0071] Example 2: The electron cyclotron radiometer calibration method proposed in this example is based on plasma displacement control. Since plasma displacement can be accurately measured by magnetic diagnosis under low specific pressure conditions, the measurable functional relationship between the time sequence and the measurement space interval between different channels can be used to find the different times when different channels are measured in the same plasma space. Then, the calibration coefficients of all channels can be obtained through recursive relationships. Depending on the different tokamaks, the electron cyclotron radiometer calibration based on plasma displacement control has the following four steps:

[0072] S1, first of all, we need to calculate the measurement position intervals of different channels of the electron cyclotron radiometer according to the toroidal magnetic field distribution of different tokamaks; the measurement positions of different channels of the electron cyclotron radiometer depend on its local oscillator frequency F0, intermediate frequency f N , the intermediate frequency is the intermediate frequency of the Nth channel of the electron cyclotron radiometer, and it also depends on the toroidal magnetic field intensity distribution (toroidal magnetic field distribution). In the tokamak device, the toroidal magnetic field intensity distribution is a good inverse proportional function in the radial direction, and its toroidal magnetic field distribution can be expressed by the following formula:

[0073]

[0074] Where B represents the toroidal magnetic field distribution, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, and R represents the measurement position. The secondary X-mode frequency measured by each channel of the electron cyclotron radiometer can be expressed as follows:

[0075]

[0076] Where e represents the electron charge, m e Representing the electron mass, by substituting equation (1) into equation (2), we can obtain the measurement space position (i.e., measurement position) of different channels:

[0077]

[0078] Specifically, let the measurement position difference between the Nth channel and the N+1th channel be I N ,but:

[0079]

[0080] Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer, f N+1 Indicates the intermediate frequency of the N+1th channel of the electron cyclotron radiometer.

[0081] S2, the second step is to measure the displacement of the plasma. The change of its displacement with time is usually a trigonometric function or a sawtooth change. Then, a functional relationship between the plasma displacement and time is established. From the plasma core with R0 as the time reference, a functional relationship between the plasma displacement and time is established. Since the electron temperature is a good magnetic surface function, when the plasma is horizontally displaced, the radial distribution of the electron temperature also fluctuates horizontally. Therefore, the radial position change of each spatial point is consistent. The plasma displacement function is expressed as: D = f(t), (5); where t represents time and D is the plasma displacement function. Its specific form is determined by the preset waveform for plasma displacement control.

[0082] S3, the third step is to calculate the time when each channel of the electron cyclotron radiometer overlaps with the adjacent channel in the plasma spatial position under the fluctuation of the plasma horizontal displacement. By equating equation (4) and equation (5), we can get the time difference t when the Nth channel and the N+1th channel of the electron cyclotron radiometer measure the same plasma position. N , the period of plasma horizontal displacement can be set to T0, and the time when each channel overlaps with the adjacent channel in measuring the plasma spatial position can be solved by the following formula:

[0083] I N =D(t N +t0+n×T0), (6);

[0084] Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, D represents the first functional relationship, t NIt represents the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, n represents a positive integer, and T0 represents the period of horizontal displacement of the plasma.

[0085] S4. Finally, the calibration coefficients of each channel of the electron cyclotron radiometer are obtained by averaging and comparing the values ​​of each channel with its adjacent channels at the moment when the corresponding plasma space overlaps. The calibration coefficients of each channel are then obtained through a recursive relationship. Since the electron temperature measured in the same plasma space is consistent, the radiation signal intensity measured in adjacent channels multiplied by the calibration coefficient is equal. For the first and second channels, the following recursive relationship holds:

[0086]

[0087] Wherein, E1 represents the radiation signal intensity of the toroidal magnetic field distribution received by the first channel of the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, t1 represents the time difference between the first channel and the second channel of the electron cyclotron radiometer, T0 represents the period of horizontal displacement of the plasma, n represents a positive integer, E2 represents the radiation signal intensity of the toroidal magnetic field distribution received by the second channel of the electron cyclotron radiometer, C1 represents the calibration coefficient of the first channel of the electron cyclotron radiometer, and C2 represents the calibration coefficient of the second channel of the electron cyclotron radiometer.

[0088] Specifically, the calibration coefficient C1 of the first channel can be set to 1, so the calibration coefficient of the second channel is:

[0089] Since the calibration coefficient C1 is 1,

[0090] Furthermore, the calibration coefficients of the above channels can be expressed as:

[0091]

[0092] Where, E N It represents the radiation signal intensity of the toroidal magnetic field distribution received by the Nth channel of the electron cyclotron radiometer, t N Indicates the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, C N Indicates the calibration coefficient of the Nth channel in the electron cyclotron radiometer, C N+1 Indicates the calibration coefficient of the N+1th channel in the electron cyclotron radiometer, E N+1 It represents the radiation signal intensity of the toroidal magnetic field distribution received by the N+1th channel of the electron cyclotron radiometer, t N+1It represents the time difference between the N+1th channel and the N+2th channel in the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, and n represents a positive integer.

[0093] Specifically, it is assumed that the maximum radius R0 of the tokamak device is 178 cm, the central toroidal magnetic field strength B0 is 1.7 T, the electron cyclotron radiometer has 16 channels, the local oscillator frequency is 84 GHz, the intermediate frequency frequency is 1-16 GHz, and the intermediate frequency interval is 1 GHz; the measurement position difference between the Nth channel and the N+1th channel is I N It can be calculated by formula (3), as Figure 2 As shown, the absolute measurement position of the electron cyclotron radiometer remains unchanged. As the plasma moves horizontally, the electron temperature at the same plasma spatial position will be measured at different times between adjacent channels.

[0094] Furthermore, the plasma horizontal displacement control is designed. Assuming that the horizontal displacement control waveform is a sine function with an amplitude of 5 cm and a period T0 of 0.1 second, the plasma displacement function can be expressed as: D = 5sin(20×π×t); when the plasma displacement and the measurement position interval between adjacent channels are consistent, the time difference when adjacent channels measure the same plasma space can be obtained as follows: Figure 3 As shown, according to the displacement function of the plasma, then: Therefore, we can know that for the Nth and N+1th channels, the electron temperature measured at the time t0+n×T0 when the plasma core position and the center position coincide with each other in the Nth channel, and the electron temperature measured at the time t0+n×T0 when the Nth channel and the N+1th channel coincide with each other are the same. N The electron temperature measured at time +n×T0 is the same.

[0095] Furthermore, by finally bringing the time-series radiation signal intensity measured by each channel and equation (11) into equation (9), the calibration coefficient of the first channel can be set to 1; thus, the relative calibration coefficients of all channels can be accurately obtained, and then the core electron temperature is given by the Thomson scattering signal, or the kinetic pressure is calculated to obtain the absolute calibration coefficients of all channels of the electron cyclotron radiometer.

[0096] Example 3: The present application provides an electron cyclotron radiometer calibration system based on plasma displacement control, which is applied to the electron cyclotron radiometer calibration method based on plasma displacement control in any one of Example 1 or Example 2, such as Figure 4 As shown, this may include:

[0097] The first module is used to calculate the measurement space position corresponding to each channel of the electron cyclotron radiometer under the toroidal magnetic field distribution in the tokamak device, and use the measurement space position of each channel to calculate the measurement position interval between different channels.

[0098] The second module is used to establish a first functional relationship between plasma displacement and time based on a preset waveform for controlling plasma displacement in the tokamak device.

[0099] The third module is used to calculate the time difference when adjacent channels measure the same plasma position under the horizontal displacement fluctuation of each channel by using the measurement position interval between different channels in the electron cyclotron radiometer and the first functional relationship.

[0100] The fourth module is used to obtain the recursive relationship between adjacent channels based on the condition that the measured temperature of the same plasma in each measurement space position is consistent, and to obtain the calibration coefficients of each channel in the electron cyclotron radiometer through the time difference and recursive relationship between adjacent channels and the preset value of the first channel in the electron cyclotron radiometer.

[0101] Example 4: An embodiment of the present application provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method of any one of Example 1 or Example 2 is implemented.

[0102] Example 5: The embodiment of the present application provides a non-transitory computer-readable storage medium, which stores computer instructions, and the computer instructions enable a computer to execute the method of any one of Example 1 or Example 2.

[0103] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. Electron cyclotron radiometer calibration method based on plasma displacement control, characterized in that: The specific steps include: Based on the toroidal magnetic field distribution in the tokamak device, the measurement space position corresponding to each channel of the electron cyclotron radiometer under the toroidal magnetic field distribution is calculated, and the measurement position interval between different channels is calculated using the measurement space position of each channel; Establishing a first functional relationship between plasma displacement and time based on a preset waveform for controlling plasma displacement in the tokamak device; Using the measurement position intervals and the first functional relationship between different channels in the electron cyclotron radiometer, the time difference when adjacent channels measure the same plasma position under the condition of plasma horizontal displacement fluctuation is calculated; Based on the condition that the measured temperature of the same plasma at each measurement space position is consistent, a recursive relationship between adjacent channels is obtained, and calibration coefficients of each channel in the electron cyclotron radiometer are obtained using the time difference and the recursive relationship between the adjacent channels and a preset value of the first channel in the electron cyclotron radiometer, and each channel in the electron cyclotron radiometer is calibrated using the calibration coefficients of each channel; The measurement position intervals between the different channels are specifically: Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer, f N+1 Indicates the intermediate frequency of the N+1th channel of the electron cyclotron radiometer; The time difference is obtained by: I N =D(t N +t0+n×T0), where D=f(t); Where, I N represents the measurement position interval between the Nth channel and the N+1th channel in the electron cyclotron radiometer, D represents the first functional relationship, t N represents the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, n represents a positive integer, and T0 represents the period of horizontal displacement of the plasma; The recursive relationship is specifically: Wherein, E1 represents the radiation signal intensity of the toroidal magnetic field distribution received by the first channel of the electron cyclotron radiometer, t0 represents the moment when the core position and the center position of the plasma coincide, t1 represents the time difference between the first and second channels of the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, n represents a positive integer, E2 represents the radiation signal intensity of the toroidal magnetic field distribution received by the second channel of the electron cyclotron radiometer, C1 represents the calibration coefficient of the first channel of the electron cyclotron radiometer, and C2 represents the calibration coefficient of the second channel of the electron cyclotron radiometer; The calibration coefficients of the various channels are specifically: Where, E N It represents the radiation signal intensity of the toroidal magnetic field distribution received by the Nth channel of the electron cyclotron radiometer, t N Indicates the time difference between the Nth channel and the N+1th channel in the electron cyclotron radiometer, C N Indicates the calibration coefficient of the Nth channel in the electron cyclotron radiometer, C N+1 Indicates the calibration coefficient of the N+1th channel in the electron cyclotron radiometer, E N+1 It represents the radiation signal intensity of the toroidal magnetic field distribution received by the N+1th channel of the electron cyclotron radiometer, t N+1 It represents the time difference between the N+1th channel and the N+2th channel in the electron cyclotron radiometer, T0 represents the period of plasma horizontal displacement, and n represents a positive integer.

2. The electron cyclotron radiometer calibration method based on plasma displacement control according to claim 1, characterized in that: The toroidal magnetic field distribution is specifically as follows: Where B represents the toroidal magnetic field distribution, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, and R represents the measurement position.

3. The electron cyclotron radiometer calibration method based on plasma displacement control according to claim 2, characterized in that: The measurement positions are specifically: Where R N represents the measurement position of the Nth channel in the electron cyclotron radiometer, e represents the electron charge number, R0 represents the central maximum radius of the tokamak device, B0 represents the central toroidal magnetic field intensity of the toroidal magnetic field distribution, m e represents the electron mass, F0 represents the local oscillator frequency of the electron cyclotron radiometer, f N Indicates the intermediate frequency of the Nth channel of the electron cyclotron radiometer.

4. An electron cyclotron radiometer calibration system based on plasma displacement control, applied to the electron cyclotron radiometer calibration method based on plasma displacement control according to any one of claims 1 to 3, characterized in that: include: The first module is used to calculate the measurement space position corresponding to each channel of the electron cyclotron radiometer under the toroidal magnetic field distribution in the tokamak device, and use the measurement space position of each channel to calculate the measurement position interval between different channels; The second module is configured to establish a first functional relationship between plasma displacement and time based on a preset waveform for controlling plasma displacement in the tokamak device; The third module is used to calculate the time difference when adjacent channels measure the same plasma position under the horizontal displacement fluctuation of each channel by using the measurement position interval between different channels in the electron cyclotron radiometer and the first functional relationship; The fourth module is used to obtain the recursive relationship between adjacent channels based on the condition that the measured temperature of the same plasma in each measurement space position is consistent, and to obtain the calibration coefficient of each channel in the electron cyclotron radiometer through the time difference and recursive relationship of adjacent channels and the preset value of the first channel in the electron cyclotron radiometer, and to calibrate each channel in the electron cyclotron radiometer using the calibration coefficient of each channel.

5. An electronic device, characterized in that: The method comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the method according to any one of claims 1 to 3 is implemented when the processor executes the computer program.

6. A non-transitory computer-readable storage medium, characterized in that The non-transitory computer-readable storage medium stores computer instructions, and the computer instructions enable a computer to execute the method according to any one of claims 1 to 3.

Images