Method for determining the quantity of plutonium in the presence of curium by passive neutron measurement

Abstract

A method for determining the quantity of plutonium in a radioactive sample in the presence of curium, comprising the steps of: Measuring (301) a train of electrical pulses using a radiation detection system, Determining (302), respectively in a first and a second time window, a first and a second number of second-order pulse coincidences between the electrical pulses supplied by the detectors, The second time window being formed by the union of the first time window and a third time window, The first and third time windows being chosen such that the number of second-order pulse coincidences related to plutonium is greater than that related to curium in one of the two windows between the first and third time windows and less than that related to curium in the other window,Determine (304) the quantity of plutonium by applying a nonlinear regression model taking as parameters: the total number of measured pulses, the first number of second-order pulse coincidences, and the second number of second-order pulse coincidences, the model having been previously trained on a set of values ​​simulated from a numerical model. Figure 3,

Description

Title of the invention: Method for determining the quantity of plutonium in the presence of curium by passive neutron measurement

[0001] The invention relates to the field of passive neutron measurement for the characterization of radioactive waste or nuclear materials (solid or in solution). More specifically, it relates to the measurement of the coincidences of neutrons and gamma radiation from spontaneous fissions for the determination of the mass of plutonium in packages of radioactive waste, samples of nuclear materials (ingots, pots of powder, pellets or fuel rods, machining scrap, or any other form), or even solutions of nuclear materials (spent fuel reprocessing processes).

[0002] The possible presence of 242Cm and 244Cm in packages of radioactive waste and nuclear materials can lead to a significant overestimation of the amount of plutonium. Indeed, to take the example of 244Cm, which is more frequently present (radioactive half-life of 18.1 years) than 242Cm, which decays rapidly (half-life of 162 days), the specific neutron emission (ns'.g1) by spontaneous fission of 244Cm is approximately 10,000 times greater than that of 240Pu, the main neutron emitter of plutonium by spontaneous fission. Consequently, the presence of curium in a waste package can potentially mask the neutron signal resulting from the spontaneous fission of even-numbered plutonium isotopes (240Pu, but we are only referring to the main emitter 240Pu for simplicity).In this case, evaluating the total plutonium mass from the equivalent mass of 240Pu (the mass that would emit as many neutrons as all the '' Pu isotopes combined) becomes highly uncertain, with a significant risk of overestimation. This poses problems at various levels, for example, for verifying a safety-criticality criterion or for materials accounting when considering fissile mass (239 and 241Pu isotopes), or for managing radioactive waste storage sites (initial acceptance criteria and long-term alpha activity, particularly for the 239 and 240Pu isotopes). In this type of case, it is not feasible to use prior information (i.e.(not measured, provided by the holder of the nuclear materials or the producer of the waste) on the relative proportions of plutonium and curium, because even a small uncertainty on this mass proportion would induce too large an error in estimating the quantity of plutonium.

[0003] In order to accurately estimate the mass of plutonium despite the presence of curium in radioactive waste or nuclear materials, there are two types of measurement methods: either an active neutron measurement, or a passive neutron measurement.

[0004] Active neutron measurement is based on the use of a pulsed neutron generator. This method, described in reference

[10] , is currently the recommended solution for obtaining a sufficiently high ratio between the signal (prompt neutrons from fission induced on plutonium, but also uranium) and the background noise (primarily passive, due to spontaneous fission of curium). However, implementing active measurement entails significantly greater technical constraints and costs than passive measurement. Indeed, the price of a neutron generator is high and depends on the required emission level; moreover, its ownership and use are subject to regulatory authorizations requiring, in particular, specific operator training, radiation protection measures (irradiation chamber) and physical protection of the device (dual-use equipment), as well as accounting for the tritium source it contains.

[0005] For these reasons, it is less expensive to search for a method based on a passive neutron measurement.

[0006] The vast majority of passive neutron measurement stations intended for the characterization of plutonium (nuclear materials, radioactive waste packages, process solutions, etc.), by detecting spontaneous fission neutrons emitted in coincidence, are equipped with 3He gaseous proportional counters because they are very sensitive to thermal neutrons and not very sensitive to gamma radiation.

[0007] In some cases, the parasitic presence of curium, a very intense emitter of neutrons by spontaneous fission, leads to a very significant overestimation of the quantity of plutonium.

[0008] Currently, there are two main passive neutron measurement methods implemented with 3He detectors for evaluating the mass of plutonium in the presence of curium, both based on the fact that curium isotopes emit more neutrons by spontaneous fission than plutonium isotopes (2.72 on average versus 2.16, with different multiplicity distributions see reference [4]).

[0009] These methods are based on the exploitation of pulse coincidence measurements from a pulse train measured by the detectors. A pulse coincidence of order N or higher, with N a strictly positive integer, is detected for each pulse when, in a time window of given size starting at the instant of arrival of the pulse, Nl additional pulses are detected in the time window.

[0010] A first method of evaluation by passive neutron measurement consists of evaluating the ratio between the count of second-order real pulse coincidences corresponding to neutron pairs and the total pulse count. This method, described in reference [5], is generally used for stations measurement systems that do not have very high detection efficiency or neutron multiplicity counting systems.

[0011] In this first method, a significant increase in the ratio between actual pulse coincidences and total count can be observed in the presence of 244Cm. However, this ratio evolves inversely with the apparent neutron detection efficiency. Consequently, any uncontrolled decrease in this apparent efficiency (such as in the case of waste packages containing elements that slow down and then absorb neutrons, or which contain only moderators that make neutrons more difficult to detect when the measurement station contains an absorbent such as cadmium in front of the detectors) could be wrongly interpreted as the presence of curium. Thus, this method generally does not allow for a precise estimation of the quantity of plutonium in the presence of curium and provides only qualitative information on the possible presence of curium and the associated risk of overestimating the plutonium mass.

[0012] A second method of evaluation by passive neutron measurement consists of measuring neutron multiplicities and in particular the ratio between the count of second-order real impulse coincidences (called doublets) and that of third-order real impulse coincidences (called triplets). This second method is described, for example, in references [1] and [8].

[0013] This second method also exploits the difference between the probability distribution of the number of neutrons emitted by spontaneous fission of 244Cm and that associated with 240Pu (the two main emitters of Cm and Pu). The ratio between the counts of triplets and doublets is higher for 244Cm than for 240Pu, which allows them to be distinguished a priori. Reference [1] indicates that it is possible to evaluate the quantity of plutonium down to a maximum fraction of 244Cm atoms of 1000 ppm, i.e., a maximum ratio between the masses of 244Cm and 240Pu of the order of 0.1%.Furthermore, implementing this method requires the use of neutron measurement stations equipped with a very large number of 3He counters (up to 5 rings of 40 3He counters surrounding a 100-liter waste drum) in order to obtain very high detection efficiencies (number of counts per emitted neutron), from 30% to 40%, and consequently statistically usable triplet count rates. However, in most cases, particularly in radioactive waste packages where matrix effects reduce this detection efficiency and where the ratio between the masses of Cm and Pu can reach several percent, this method, as before, provides only qualitative information on the risk of the presence of 244Cm, as also indicated in reference [8].

[0014] The Applicant's patent application published under number EP 3 835 831 A1 relates to a method for measuring plutonium by the passive detection of neutron coincidences with PVT (polyvinyl toluene, fast neutron sensitive) plastic scintillators as a replacement for 3He detectors, primarily due to their excessive cost and their sensitivity to accidental coincidences (much longer coincidence windows). The proposed method is based on a temporal selection of coincidence windows, taking into account, for the counting of triplets, the arrival times between the three detected particles (neutrons and prompt gamma rays) in order to better discriminate triplets due to spontaneous fissions from those due to various parasitic reactions such as (a,n) reactions, but also elastic and inelastic scattering creating crosstalk and therefore undesirable coincidences.However, this method does not take into account estimation errors related to the presence of curium 244Cm.

[0015] There is therefore a need for a method to accurately estimate the quantity of 240Pu in the presence of 244Cm in nuclear materials or in radioactive waste drums, by passive measurement of neutron coincidences, in cases where the prior art did not allow this, namely with a mass ratio of 244Cm to 240Pu greater than 0.1%. Furthermore, it is necessary to eliminate spurious gamma contributions due to fission or activation products (such as 137Cs, 60Co) potentially present in the objects to be measured, since the plastic scintillators used in detectors are sensitive to this type of radiation.

[0016] The invention proposes a new method for determining the mass of plutonium by passive neutron measurement which exploits the technique described in reference EP 3 835 831 Al but improves it in the presence of curium.

[0017] For this purpose, the invention is based on a nonlinear regression model with several explanatory variables which are real coincidence numbers of different orders counted in two different time intervals, sized so that, in one of the two intervals, the contribution of plutonium is different (higher or lower) than that of curium, and in the other interval, it is the reverse.

[0018] In the following description, the coincidence numbers of real impulses of order 2 are also called doublets. The coincidence numbers of real impulses of order 3 are also called triplets.

[0019] Taking into account two time intervals to calculate the doublets and triplets associated with the training of a non-linear regression model parameterized by explanatory variables calculated on these two distinct intervals, makes it possible to better characterize the mass of plutonium in the presence of curium unlike prior art methods.

[0020] The invention relates to a method for determining the quantity of plutonium in a radioactive sample in the presence of curium, the method comprising the steps of: - To measure a train of electrical pulses using a radiation detection system comprising several detectors arranged around the sample and capable of generating an electrical pulse in response to the detection of radiation, - Determine, respectively in a first and a second time window, a first and a second number of second-order pulse coincidences between the electrical pulses provided by the detectors, - The second time window being formed by the union of the first time window and a third time window, - The first and third time windows are chosen such that the number of second-order pulse coincidences linked to plutonium is greater than that linked to curium in one of the two windows (the first and third time windows) and less than that linked to curium in the other window, - Determine the quantity of plutonium by applying a non-linear regression model taking as parameters: the total number of measured pulses, the first number of second-order pulse coincidences and the second number of second-order pulse coincidences, the model being previously trained on a set of simulated values ​​from a numerical model.

[0021] According to a particular aspect of the invention, the second time window is chosen so as to include pairs of pulses corresponding to radiations of type (y,n) and (n,n) and to exclude pairs of pulses corresponding to radiations of type (y,y).

[0022] In one embodiment, the process according to the invention further comprises the steps of: - Determine, respectively in the first and second time windows, a first and a second number of third-order pulse coincidences between the electrical pulses provided by the detectors, - The nonlinear regression model further taking as parameters the aforementioned first and second numbers of third-order impulse coincidences.

[0023] In an embodiment, the method according to the invention further comprises the steps of: Determine a first number of third-order pulse coincidences between the electrical pulses supplied by the detectors such that The time interval between the first and second pulses is included in the first time window, and the time interval between the second and third pulses is included in the first time window. - Determine a second number of third-order pulse coincidences between the electrical pulses provided by the detectors such that the time gap between the first and second pulses is included in the second time window and the time gap between the second and third pulses is included in the second time window. - The nonlinear regression model further taking as parameters said first and second numbers of coincidence of 3rd order pulses.

[0024] According to a particular aspect of the invention, the second time window is chosen so as to include triplets of pulses corresponding to radiations of type (y,n,n) and (n,n,n) and to exclude triplets of pulses corresponding to radiations of type (y,y,y).

[0025] According to a particular aspect of the invention, the first time window is equal to [10 ns ; 20 ns], the second time window is equal to [10 ns ; 60 ns] and the third time window is equal to [20 ns ; 60 ns].

[0026] According to a particular aspect of the invention, the first time window is equal to [20 ns ; 60 ns], the second time window is equal to [10 ns ; 60 ns] and the third time window is equal to [10ns ; 20 ns].

[0027] The invention also relates to a system for determining the quantity of plutonium in a radioactive sample in the presence of curium, the system comprising several detectors arranged around the sample and capable of generating an electrical pulse in response to the detection of radiation and a computer configured to execute the steps of the process according to the invention.

[0028] According to a particular aspect of the invention, each detector comprises a plastic organic scintillator made of polyvinyl-toluene (PVT).

[0029] Other features and advantages of the present invention will become more apparent from the following description in relation to the following accompanying drawings.

[0030] [Fig-1] represents a first diagram of a passive neutron measurement system according to an embodiment of the invention,

[0031] [Fig.2] represents a second diagram of a passive neutron measurement system according to an embodiment of the invention,

[0032] [Fig.3] represents a flowchart detailing the steps for implementing a method for determining a quantity of plutonium according to an embodiment of the invention,

[0033] [Fig.4a] represents an example of coincidence determination for a pulse train,

[0034] [Fig.4b] represents a table of multiplicity order distributions for the example of [Fig.4a],

[0035] [Fig.5] represents an example of a Rossi-Alpha curve,

[0036] [Fig.6] represents an example of a neutron energy spectrum diagram of spontaneous fission of 244Cm and 240Pu,

[0037] [Fig.7] represents the respective Rossi-Alpha curves for 244Cm and 240Pu,

[0038] [Fig.8] illustrates a first example of estimating the mass of 240Pu using a first type of non-linear regression model,

[0039] [Fig.9] illustrates a first example of estimating the mass of 240Pu using a second type of non-linear regression model,

[0040] [Fig. 10] illustrates a first example of estimating the mass of 240Pu using a third type of non-linear regression model,

[0041] [Fig. 11] illustrates a second example of estimating the mass of 240Pu using the second type of nonlinear regression model,

[0042] [Fig. 12] illustrates another example of estimating the mass of 240Pu using a fourth type of nonlinear regression model,

[0043] [Fig. 13] illustrates the calculation of a third-order impulse coincidence number with time selection from a two-dimensional histogram according to a prior art principle,

[0044] [Fig. 14] illustrates another example of estimating the mass of 240Pu using a fifth type of nonlinear regression model,

[0045] Fig. 1 represents a simplified diagram of a passive neutron measurement system according to one embodiment of the invention.

[0046] In this non-limiting example, the measuring system comprises 4 detectors, each detector comprising a PVT plastic S1-S4 scintillator and a P1-P4 photomultiplier. The number of detectors can be higher, for example equal to 10 or 16.

[0047] The measurement system also includes a computer 2 connected to the photomultipliers P1-P4 via analog-to-digital converters 3. The system 1 may also include a human-machine interface module 4, for example a screen on which the measurement results are displayed.

[0048] Figure 2 represents another example of a diagram of a measurement system Passive neutron detection, this time comprising 16 detectors 203 based on PVT plastic scintillators. The system also includes a lead screen 201 arranged in the form of a tube inside the detectors and a barrel 202 in which a sample is placed. The system rests on a support pallet 204. In an alternative embodiment, the lead screen 201 can be omitted or its thickness can be adjusted according to the intensity of the gamma radiation emitted by the sample to be analyzed.

[0049] The measurement system according to the invention is equipped with PVT plastic scintillators (without neutron-gamma discrimination) as described in the patent application referenced [2].

[0050] The measurement system preferentially uses polyvinyl-toluene PVT plastic scintillators and not 3He detectors for the following reasons.

[0051] The use of 3He detectors imposes time constants on the order of several tens of microseconds, which correspond to the time required for the thermalization of neutrons, unlike plastic scintillators which directly detect fast neutrons and for which the time constants are on the order of tens of nanoseconds. This property is used in references [2] and [3] to limit the number of accidental coincidences and to differentiate the plutonium signal from that of parasitic gamma sources (presence of fission or activation products in the measured object) and neutron sources ((α,ν) reactions on light nuclei such as oxygen in plutonium oxides, in the presence of intense alpha emitters such as α'Am or Pu).

[0052] Figure 3 represents a flowchart detailing the steps for implementing a process for determining the quantity of plutonium in the presence of curium in a sample using the system described in Figures 1 and 2.

[0053] The process begins in step 301 with a coincidence measurement step performed using the system described above for a given sample. Each detector includes a photomultiplier tube capable of providing an electrical signal representative of the light signal generated in the plastic scintillator. This electrical signal is recorded for each detector. At the end of this measurement step, a pulse train is obtained, each pulse being characterized at least by its arrival time and its radiation energy.

[0054] In step 302, from the pulse train measured in step 301, a number of coincidences of order 2 is determined for two particular time windows, the choice of which will be explained later.

[0055] In step 303, an optional number of coincidences of order 3 is also determined for the same time windows as in step 302.

[0056] Finally, in step 304, a nonlinear regression model is run which is trained to determine a quantity of plutonium present in the sample from three parameters: the total number of pulses and the two second-order coincidence numbers determined in step 302 or from five parameters: the three parameters mentioned above and the two third-order coincidence numbers determined in step 303.

[0057] Fig. 4a illustrates, on an example, the principle of determining a number of N-order pulse coincidences for a measured pulse train of 400.

[0058] According to a known principle, notably described in reference [5], a coincidence of pulses of order N exists when, for a given pulse having an arrival time t0, there are N - 1 other pulses in a time window of a given size whose beginning coincides with the instant t0. The number N varies from 0 to Nd-1 where Nd is the number of detectors in the system.

[0059] In the example of [Fig.4a], the time window has a size of 100 ns and we have identified 0th order multiplicities or 1st order coincidences (no other pulse in the window), 1st order and 2nd order multiplicities (respectively 1 and 2 additional pulses in the window).

[0060] Furthermore, a distinction is made between the so-called "real and accidental" R+A windows, which are applied to all measured pulses, and the so-called "accidental" A windows, which correspond to accidental coincidences for time windows located beyond a certain time lag, taken to be approximately 10 times the detection time of a neutron in the system. In the example of [Fig. 4a], this time lag is equal to 1000 ns.

[0061] To determine the number of coincidences of order N, we first determine the number of multiplicities MN for N ranging from 0 to a predefined maximum integer value max. We obtain a distribution of the multiplicity orders given in the table in [Fig. 4b] for the particular example in [Fig. 4a]. From this distribution, we determine factorial moments of order k defined by the following relation:

[0062] (K+

[0063] (4) ^ymax^LM . 'A4 - L^N=k^kJlylN (A)

[0064] Where:

[0065] (R+A)k and (A)k are respectively the factorial moments of order k of the distributions (R+A) and (A)

[0066] Mn (R+A) and MN(A) are respectively the orders N of the multiplicities of the distributions (R+A) and (A),

[0067] “max” is the maximum order of multiplicity which is, for example, equal to the number of system detectors.

[0068] The real coincidence numbers of orders 1, 2 and 3, also called singlets, doublets and triplets (S, D and T), are then given by the following equations (see references [5] [7]):

[0069] 5= (R + A)^^

[0070] D= (R + A^À^

[0071] T=l [(r+a)2-(A)2-^[(R+A)1-(A)))

[0072] S corresponds to the total number of pulses in the pulse train.

[0073] Fig. 5 shows an example of a Rossi-Alpha curve which gives a histogram of the number of coincidences of second-order pulses as a function of the value of the time window.

[0074] In other words, this histogram is obtained by measuring, for each detection, the arrival time of the next pulse. The example in [Fig. 5] is obtained for a plutonium-240Pu source in a detection system composed of PVT plastic scintillators.

[0075] On the curve, we identified a first time window [10-60] ns which corresponds to the counting of real + accidental coincidences and a second time window [260-310] ns which corresponds to the counting of accidental coincidences.

[0076] In the time interval [0-10] ns, the coincidences of order 2 corresponding to the pairs yy are concentrated. The time interval [10-60] ns corresponds to the pairs yn and nn.

[0077] In order to be able to quantify precisely the plutonium in the presence of curium, it is proposed to exploit two distinct time intervals to count the number of coincidences of second-order pulses.

[0078] Indeed, the present invention is based on the differences between the energy spectra of spontaneous fission neutrons of 244Cm and 240Pu.

[0079] Fig. 6 shows the 601,602 energy spectra of spontaneous fission neutrons of 240Pu and 244Cm.

[0080] We can see in this figure a significant difference between the emission probability of 240Pu (curve 601) and of 244Cm (curve 602) in the energy ranges [0.2 - 2] MeV and [2-7] MeV.

[0081] Figure 7 represents a comparison of the Rossi-Alpha curves for the same doublet values ​​measured over a time window [10-60] ns between 240Pu (701) and 244Cm (702).

[0082] Figure 7 shows the average arrival time of neutrons and gamma rays. fission following a first detection which is generally that of gamma radiation due to their higher speed (30 cm.ns') than that of neutrons (a few cm.ns', depending on their energy).

[0083] A first peak is observed, corresponding to the detection of gamma-ray pairs, followed by a "bump" corresponding to that of pairs involving at least one neutron, with a shape comparable to that of the spectra in [Fig. 6] because the neutron velocity varies as the square root of their energy. A significant difference between 240Pu and 244Cm appears in [Fig. 7] in terms of the number of coincidences. of order 2 over the [20-60] ns window, and also over the [10-20] ns window. Thus, we observe a higher probability of detecting neutrons from the spontaneous fission of 244Cm in the coincidence window between 10 and 20 ns than for 240Pu, the trend being reversed between 20 and 60 ns. Indeed, the curves 701,702 intersect approximately at time t=20 ns.

[0084] More specifically, on the numerical example of [Fig.7], the values ​​of real coincidence numbers of order 2 are given in the table below for the two time intervals, respectively for a plutonium sample and for a curium sample each having a mass equivalent to 400 g of 240Pu, for the same number of coincidences in the total window [10-60] ns. 240Pu (400 g) 244Cm (400 g of 240Puéq) D [10-20] ns 286155 299059 D [20-60] ns 553343 535718

[0085] By exploiting this difference in contribution over two distinct time windows, it is possible to improve the prediction of the mass of 40Pu in a sample in the presence of 244Cm using a non-linear regression model.

[0086] Indeed, taking into account a first time window in which the contribution of 240Pu is greater than that of 244Cm and a second time window in which the contribution of 244Cm is greater than that of 240Pu as parameters of a non-linear regression model makes it possible to improve the prediction of the mass of 40Pu as will be demonstrated later.

[0087] The precise determination of the boundary between the two time windows (here equal to 20 ns) is done by experimentation by seeking the best compromise of performance, for example by testing different time windows and comparing the predictions of the regression model obtained.

[0088] For example, the two time windows can be defined such that the ratio of the respective coincidence numbers of the two elements 240Pu and 244Cm is greater than a first predefined threshold for one of the windows and less than a second predefined threshold for the other window. In the example above, the ratio between the contributions of 244Cm and 240Pu is 1.04 for the [10-20] ns window and 0.97 for the [20-60] ns window.

[0089] Step 304 of the method according to the invention consists of running a previously trained nonlinear regression model to determine a mass or quantity of 240 Pu from several coincidence number parameters.

[0090] The model is, for example, a multilinear least squares regression model or an artificial neural network as described in reference

[11] or more generally any learning model that can be trained to predict a quantity of plutonium from several measurements of pulse coincidence numbers.

[0091] Figures 8, 9 and 10 illustrate prediction results obtained with a nonlinear regression model based on a Monte Carlo method for simulated data corresponding to a source composed of plutonium and curium.

[0092] The number of configurations corresponds to a full factorial experimental design composed of three variable parameters, namely the masses of 244Cm and 240Pu, and the activity of a 60Co source (emitting two correlated gamma rays at 1173 keV and 1332 keV). A full factorial design comprising a large number of configurations of the aforementioned values ​​is used as a training basis for the prediction model. Once the model has been trained on this basis, it can be used to estimate a mass of 240Pu from new measurements performed on a new, unknown sample.

[0093] The training data provide at a minimum the masses of 244Cm and 240Pu, as well as the values ​​of the explanatory variables of the model (total number of pulses, number of coincidences in two distinct intervals). The training data can be measured using the device described in Figures 1 and 2 on several samples or can be simulated with a particle transport code based on the Monte Carlo method such as that described in reference [9].

[0094] Fig. 8 shows, on a diagram, the prediction results (on the ordinate) as a function of the actual value of the mass of 240Pu for a regression model taking as inputs only the total number of pulses (also called singlets S) and the number of coincidences of second-order pulses counted in a single time window equal to [10-60] ns (also called doublets D).

[0095] As can be seen in [Fig. 8], the prediction results are highly dispersed and do not closely match the actual values. This is because the model fails to discriminate the contribution of plutonium from that of curium in the complete [10-60] ns window in which 244Cm and 240Pu exhibit the same number of second-order coincidences.

[0096] Fig. 9 shows the prediction results when the model takes three parameters as input: the total number of pulses, a first number of real coincidences of second-order pulses counted in a first time window equal to [10-20] ns and a second number of real coincidences of second-order pulses counted in a second time window equal to [10-60] ns.

[0097] Figure 10 shows, alternatively, the prediction results when the model takes three parameters as input: the total number of pulses, a first number of real coincidences of second-order pulses counted in a first a time window equal to [20-60] ns and a second number of real coincidences of second-order pulses counted in a second time window equal to [10-60] ns.

[0098] Figures 9 and 10 show that taking into account in the regression model the three explanatory variables, namely singlets on the one hand and doublets counted in two distinct coincidence windows such as D [10-20] ns and D [10-60] ns (see [Fig. 9]), or D [20-60] ns and D [10-60] ns (see [Fig. 10]), on the other hand, makes it possible to obtain a significantly better prediction of the mass of 240Pu than by using a single window D [10-60] ns ([Fig. 8]).

[0099] The predicted values ​​are indeed much closer to the true value in Figures 9 or 10 than in [Fig. 8], regardless of the amount of curium. This result is linked to the separation of the explanatory variables into two distinct windows: a global window and a window in which the influence of plutonium or curium is greater.

[0100] The results illustrated in Figures 8 to 10 are obtained with simulated data for point sources of plutonium and curium located in the center of the 118 L drum filled with cellulose of density 0.3. In order to evaluate the impact related to the multiplication effect by induced fissions, which becomes significant for plutonium masses greater than 100 g, a second experimental design with 124 configurations is simulated for the case of a sphere of PuO2 powder of density 3, still located in the center of the same organic matrix.

[0101] The simulated variables for this second experimental plan are the volume of the PuO2 powder sphere, related to the total mass of plutonium, the mass of 240Pu and the mass of 244Cm.

[0102] Fig. 11 shows the prediction results obtained with the same regression model taking as input the total number of pulses, a first number of second-order pulse coincidences counted in a first time window equal to [10-20] ns and a second number of second-order pulse coincidences counted in a second time window equal to [10-60] ns.

[0103] Even though we observe in [Fig.1 1] a degradation of the performance of the regression model compared to the results of figures 9 and 10 for masses of 240Pu equivalent greater than 100 g, due to the multiplication effect, the performance of the model composed of the 3 explanatory variables mentioned above S, D [10-20] ns, D [10-60] ns remains nevertheless quite acceptable for practical application.

[0104] More generally, the numerical values ​​of the time intervals for counting doublets depend on the detection device considered as well as the distance between the sample and the detectors.

[0105] The intervals to be taken into account as inputs to the regression model are more generally chosen as two contiguous time intervals in which respectively one of the two elements between the 240Pu or the 244Cm contributes more than the other and vice versa.

[0106] Another embodiment of the invention is now described in which the regression model takes as input 5 parameters: the 3 parameters S, D [10-20] ns, D [10-60] ns described previously and 2 other parameters which are a first number of coincidences of 3rd order pulses counted in a first time window equal to [10-20] ns, and a second number of coincidences of 3rd order pulses counted in a second time window equal to [10-60] ns.

[0107] A regression model composed of these 5 explanatory variables makes it possible to strongly limit the influence of the multiplication effect (significant for plutonium masses greater than 100 g) on ​​the estimation of the plutonium mass.

[0108] Indeed, the use of these third-order coincidences makes it possible to better exploit the difference in neutron multiplicity between 244Cm and 240Pu. It is thus possible to observe a better prediction of the mass of 240Pu with regression models consisting of 5 explanatory variables evaluated in a time window between 10 and 20 ns and another between 10 and 60 ns.

[0109] Fig. 12 shows the prediction results obtained for a model consisting of the 5 aforementioned variables.

[0110] In one embodiment, the number of 3rd order impulse coincidences can be determined by means of a two-dimensional histogram of 3rd order coincidences.

[0111] Fig. 13 illustrates such a histogram 1300 constructed from a pulse train 1301. The histogram represents all combinations of 3 coincident pulses 1302 whose time interval between the first two pulses is shown on the A2i axis and the time interval between the 2nd and 3rd pulse is shown on the A32 axis.

[0112] The number of second-order pulse coincidences can be determined by summing the values ​​of an area of ​​interest of the ROInette histogram which is defined by the desired time windows.

[0113] According to reference [3], this area of ​​interest is given by the relation:

[0114] ROInette = (ROIbrute - Ac) -(A A21 - Ac) - (A Aî 2 - Ac)

[0115] Ac is the common accidental region of histogram 1300

[0116] A A2_i is the accidental region along the first axis of the histogram

[0117] A A3 2 is the accidental region along the second axis of the histogram

[0118] ROIbrute corresponds to the counting area in the given time window.

[0119] Thus, similarly to what has been described previously, the two explanatory variables of number of coincidences of order 3 determined for the two time windows [10-20] ns and [10-60] ns can be determined via the calculation of ROInette on this histogram 1300.

[0120] [Fig. 14] shows the prediction results obtained for a model consisting of the 5 variables by calculating the number of coincidences of order 3 using the histogram of [Fig. 13].

[0121] In summary, the nonlinear regression model executed in step 304 of the process according to the invention can consist of 3 or 5 explanatory variables.

[0122] The first three variables are: - the total number of pulses which corresponds to the number of first-order coincidences (singlets), - a first number of coincidences of order 2 (doublets) calculated in a first chosen time window such that the contribution of plutonium is greater or less than that of curium, - a first number of coincidences of order 2 (doublets) calculated in a second time window which corresponds to the concatenation of the first time window and a third time window such that the respective contributions of plutonium or curium are reversed with respect to the first time window.

[0123] The two additional variables are: - a first number of coincidences of order 3 (triplets) calculated in the first time window, - a first number of coincidences of order 3 (triplets) calculated in the second time window - the number of third-order coincidences that can be calculated directly from the pulse train (triplets) or via the construction of a two-dimensional histogram (ROInette)-

[0124] For example, the first time window is [10-20] ns, the second time window is [10-60] ns, the third time window is [20-60] ns.

[0125] Alternatively, the first time window is [20-60] ns, the second time window is [10-60] ns, the third time window is [10-20] ns. List of documents cited

[0126] [1] DH Beddingfield, AP Belian, “Detection of Cm-244 in Plutonium-Bearing Wastes at reprocessing Facilities”, Los Alamos National Laboratory, 2004

[0127] [2] V. Bottau, R. De Stefano, C. Carasco, C. Eleon, B. Perot, “Detection method of radiation and associated system implementing neutron-gamma coincidence discrimination”, patent application EP 3 835 831 Al.

[0128] [3] V. Bottau, C. Carasco, B. Pérot, C. Eleon, R. De Stefano, L. Isnel, I. Tsekhanovich, “Detection of Fission Coincidences With Plastic Scintillators for the Characterization of Radioactive Waste Drums”, TNS IEEE, volume 69, issue 4, 2022

[0129] [4] N. Ensslin, “Passive Nondestructive Assay of Nuclear Materials,” 1991

[0130] [5] JB Porcher, T. Lambert, N. Saurel, H. Schoech, L. Tondut, C. Passard, G. Granier “Recommendations for the Optimization of Passive Neutron Measurements”, CEA Report, 2013, ISSN 0429-3460,

[0131] [7] N. Ensslin, “Application Guide to Neutron Multiplicity Counting”, LA-13442-M, UC-700, November 1998.

[0132] [8] Patrick M J Chard et al, “Field Examples of Waste Assay Solutions for Curium- Contaminated Wastes”, ICEM2009-16259

[0133] [9] S. A. Pozzi, E. Padovani, M. Marseguerra, “MCNP-PoliMi, a Monte-Carlo code for corrélation measurements”, Nuclear Instruments and Methods in Physics Research, A 513 (2003) 550-558.

[0134]

[10] B. Perot et al, “The characterization of radioactive waste: a critical review of techniques implemented or under development at CEA, France”, EPJ Nuclear Sci. Technol. Volume 4, 2018

[0135]

[11] Pedregosa, et al., “Scikit-learn: Machine learning in Python.” Journal of Machine Learning Research, 2825-2830, 2011.

Claims

Demands

1. A method for determining the quantity of plutonium in a radioactive sample in the presence of curium, the method comprising the steps of: - Measuring (301) a train of electrical pulses by means of a radiation detection system comprising several detectors arranged around the sample and capable of generating an electrical pulse in response to the detection of radiation, - Determining (302), respectively in a first and a second time window, a first and a second number of second-order pulse coincidences between the electrical pulses supplied by the detectors, - The second time window being formed by the union of the first time window and a third time window,- The first and third time windows being chosen such that the number of second-order pulse coincidences related to plutonium is greater than that related to curium in one of the two windows (the first and third time windows) and less than that related to curium in the other window, - Determine (304) the quantity of plutonium by applying a nonlinear regression model taking as parameters: the total number of measured pulses, the first number of second-order pulse coincidences, and the second number of second-order pulse coincidences, the model having been previously trained on a set of simulated values ​​from a numerical model.

2. A method for determining a quantity of plutonium according to claim 1 wherein the second time window is chosen so as to include pairs of pulses corresponding to radiations of type (y,n) and (n,n) and to exclude pairs of pulses corresponding to radiations of type (y,y).

3. A method for determining a quantity of plutonium according to any one of the preceding claims, further comprising the steps of: - Determining (303), respectively in the first and second time windows, a first and a second number of third-order pulse coincidences between the electrical pulses supplied by the detectors, - The nonlinear regression model further taking as parameters the said first and second numbers of third-order pulse coincidences

4. A method for determining a quantity of plutonium according to any one of claims 1 or 2, further comprising the steps of: - Determining (303) a first number of third-order pulse coincidences between the electrical pulses supplied by the detectors such that the time interval between the first and second pulses is contained within the first time window, and the time interval between the second and third pulses is contained within the first time window, - Determining (303) a second number of third-order pulse coincidences between the electrical pulses supplied by the detectors such that the time interval between the first and second pulses is contained within the second time window, and the time interval between the second and third pulses is contained within the second time window.- The nonlinear regression model, further taking as parameters the aforementioned first and second numbers of coincidences of third-order impulses,

5. A method for determining a quantity of plutonium according to any one of claims 3 or 4 wherein the second time window is chosen so as to include triplets of pulses corresponding to radiations of type (y,n,n) and (n,n,n) and to exclude triplets of pulses corresponding to radiations of type (y,y,y).

6. A method for determining a quantity of plutonium according to any one of the preceding claims wherein the first time window is equal to [10 ns; 20 ns], the second time window is equal to [10 ns; 60 ns] and the third time window is equal to [20 ns; 60 ns].

7. A method for determining a quantity of plutonium according to any one of claims 1 to 5, wherein the first time window is [20 ns; 60 ns], the second time window is equal to [10 ns ; 60 ns] and the third time window is equal to [10ns ; 20 ns].

8. System for determining the quantity of plutonium in a radioactive sample in the presence of curium, the system comprising several detectors arranged around the sample and capable of generating an electrical pulse in response to the detection of radiation and a computer configured to execute the steps of the process according to any one of the preceding claims.

9. System for determining a quantity of plutonium according to claim 8 in which each detector comprises a plastic organic scintillator made of polyvinyl-toluene (PVT).

Images