Radionuclide diffusion path simulation method and device, electronic equipment and storage medium

Abstract

The embodiment of the invention provides a radionuclide diffusion path simulation method and device, electronic equipment and a storage medium, and relates to the technical field of nuclide diffusion analysis. The method comprises the following steps: establishing a geometric discrete model of a matrix according to the geometric shape of the matrix where the radionuclide is located and internal pore distribution, and determining the temperature and the nuclide concentration of each discrete point in the geometric discrete model according to the temperature distribution of the matrix, the initial concentration distribution of the radionuclide and the geometric discrete model; and according to the temperature of each matrix discrete point and the nuclide concentration, multiple diffusion simulation is carried out, and a diffusion path of the radionuclide in the matrix is constructed. By introducing key parameters such as matrix temperature distribution, pore structure information and initial concentration distribution, high-precision prediction simulation of a nuclide diffusion path is realized.

Description

Technical Field

[0001] The present invention relates to the technical field of nuclide diffusion analysis, and in particular to a method, device, electronic equipment and storage medium for simulating the diffusion path of radioactive nuclides. Background Art

[0002] The study of the diffusion and migration behavior of nuclides in media is an important research topic in fields such as nuclear safety assessment, geological disposal of radioactive waste, and environmental remediation. Traditional diffusion models are primarily based on Fick's law, establishing a diffusion equation driven by concentration gradients. However, under complex environmental conditions, the regulation of diffusion behavior by temperature gradients is becoming a research hotspot. In particular, in geological repositories for high-level radioactive waste, the non-uniform distribution of the temperature field caused by nuclide decay heat can significantly alter the anisotropic characteristics of the diffusion coefficient tensor, thus significantly affecting the accuracy of the prediction of nuclide migration paths. Summary of the Invention

[0003] In view of this, an object of embodiments of the present invention is to provide a method, apparatus, electronic device, and storage medium for simulating the diffusion path of radionuclides to at least partially alleviate the above-mentioned problem.

[0004] In order to achieve the above objectives, the technical solutions adopted in the embodiments of the present invention are as follows: In a first aspect, an embodiment of the present invention provides a method for simulating the diffusion path of radionuclides, comprising: Acquiring a matrix in which the radionuclide is located, and determining the temperature distribution of the matrix, the geometric shape of the matrix, the internal pore distribution of the matrix, the medium type of the internal pores, and the initial concentration distribution of the radionuclide; Establishing a geometric discrete model of the substrate based on the geometric shape and the internal pore distribution; wherein the geometric discrete model includes a plurality of discrete points, and the discrete points are divided into substrate discrete points located on the substrate and pore discrete points located in the internal pores; Determining the temperature and nuclide concentration of each discrete point according to the temperature distribution, the initial concentration distribution, and the geometric discrete model; Multiple diffusion simulations are performed based on the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix.

[0005] Optionally, performing multiple diffusion simulations based on the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix includes: For each diffusion simulation, the diffusion vector and diffusion flux of each discrete point in the matrix are calculated according to the temperature and nuclide concentration of each discrete point in the matrix; updating the nuclide concentration of each discrete point of the matrix after the current diffusion simulation according to each diffusion flux; The diffusion path of the radioactive nuclide in the matrix is ​​constructed according to the diffusion vectors of each diffusion simulation.

[0006] Optionally, calculating the diffusion vector and diffusion flux of each discrete point in the matrix according to the temperature and nuclide concentration of each discrete point in the matrix includes: For each of the matrix discrete points, taking the matrix discrete point as a target matrix discrete point, finding a first discrete point with the highest temperature, a second discrete point with the lowest temperature, a third discrete point with the highest nuclide concentration, and a fourth discrete point with the lowest nuclide concentration among the matrix discrete points within a circle with the target matrix discrete point as the center and a preset radius; Calculating a temperature gradient vector of the target substrate discrete point according to the temperature of the target substrate discrete point, the temperature of the first discrete point, and the temperature of the second discrete point; Calculating a nuclide concentration gradient vector of the target matrix discrete point based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the third discrete point, and the nuclide concentration of the fourth discrete point; Calculating the cross product of the temperature gradient vector and the nuclide concentration gradient vector to obtain a diffusion vector of the target matrix discrete point in the matrix; Taking the discrete point pointed to by the diffusion vector as the diffusion discrete point of the target matrix discrete point; The diffusion flux from the target matrix discrete point to the diffusion discrete point is calculated based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the diffusion discrete point, the diffusion coefficient of the radioactive nuclide in the matrix, and the adsorption coefficient of the matrix to the radioactive nuclide.

[0007] Optionally, calculating the temperature gradient vector of the target substrate discrete point according to the temperature of the target substrate discrete point, the temperature of the first discrete point, and the temperature of the second discrete point includes: Determining whether the temperature of the target substrate discrete point is greater than or equal to the temperature of the first discrete point; If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the target substrate discrete point minus the temperature of the second discrete point; If not, determining whether the temperature of the target substrate discrete point is equal to the temperature of the second discrete point; If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the target substrate discrete point; If not, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the second discrete point.

[0008] Optionally, calculating the nuclide concentration gradient vector of the target matrix discrete point according to the nuclide concentration of the target matrix discrete point, the nuclide concentration of the third discrete point, and the nuclide concentration of the fourth discrete point includes: Determining whether the nuclide concentration of the target matrix discrete point is greater than or equal to the nuclide concentration of the third discrete point; If so, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the target matrix discrete point minus the nuclide concentration of the fourth discrete point; If not, determining whether the nuclide concentration of the target matrix discrete point is equal to the nuclide concentration of the fourth discrete point; If so, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the third discrete point minus the nuclide concentration of the target matrix discrete point; If not, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the third discrete point minus the nuclide concentration of the fourth discrete point.

[0009] Optionally, after the step of using the discrete point pointed to by the diffusion vector as the diffusion discrete point of the target matrix discrete point, the method further includes: Determining whether the diffusion discrete point belongs to the pore discrete point; If so, searching outwards from the diffusion discrete point along the diffusion vector until a first discrete point belonging to the matrix discrete point is found, and changing the diffusion discrete point to the first discrete point; The diffusion vector is changed so that the target matrix discrete point points to the first discrete point.

[0010] Optionally, the calculating the diffusion flux from the target matrix discrete point to the diffusion discrete point based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the diffusion discrete point, the diffusion coefficient of the radioactive nuclide in the matrix, and the adsorption coefficient of the matrix to the radioactive nuclide includes: Determining whether the nuclide concentration of the diffusion discrete point is 0; If so, calculating the estimated concentration of the diffusion discrete point by using the diffusion coefficient, and calculating the diffusion flux from the target matrix discrete point to the diffusion discrete point by using the adsorption coefficient; If not, the diffusion flux from the target matrix discrete point to the diffusion discrete point is calculated using the adsorption coefficient.

[0011] In a second aspect, an embodiment of the present invention provides a radionuclide diffusion path simulation device, comprising: a matrix confirmation unit, configured to obtain a matrix in which the radionuclide is located, and determine a temperature distribution of the matrix, a geometric shape of the matrix, an internal pore distribution of the matrix, a medium type of the internal pores, and an initial concentration distribution of the radionuclide; a geometric discrete model establishing unit, configured to establish a geometric discrete model of the substrate according to the geometric shape and the internal pore distribution; wherein the geometric discrete model includes a plurality of discrete points, the discrete points being divided into substrate discrete points located on the substrate and pore discrete points located in the internal pores; A discrete point temperature and nuclide concentration confirmation unit, configured to determine the temperature and nuclide concentration of each discrete point according to the temperature distribution, the initial concentration distribution, and the geometric discrete model; The diffusion simulation unit is used to perform multiple diffusion simulations according to the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix.

[0012] In a third aspect, an embodiment of the present invention provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and running on the processor, wherein the processor implements any of the above-described methods when executing the computer program.

[0013] In a fourth aspect, an embodiment of the present invention provides a storage medium having a computer program stored thereon, wherein the computer program implements any of the above-described methods when executed by a processor.

[0014] The present invention provides a method, device, electronic device, and storage medium for simulating radionuclide diffusion paths. By incorporating key parameters such as substrate temperature distribution, pore structure information, and initial concentration distribution, these methods achieve high-precision prediction and simulation of radionuclide diffusion paths. This not only improves the accuracy and applicability of simulation results, but also plays an important role in fields such as nuclear waste disposal and environmental remediation.

[0015] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments. It should be understood that the following drawings only illustrate certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without paying any creative work.

[0017] Figure 1 A schematic structural block diagram of an electronic device provided by an embodiment of the present invention; Figure 2 A schematic flow chart of a method for simulating the diffusion path of radionuclides provided in an embodiment of the present invention; Figure 3 A schematic diagram of a substrate and internal pore distribution provided by an embodiment of the present invention; Figure 4 A schematic diagram of a geometric discrete model of a substrate provided by an embodiment of the present invention; Figure 5 Another schematic flow chart of a method for simulating the diffusion path of radionuclides provided in an embodiment of the present invention; Figure 6 A schematic diagram of a diffusion vector in a geometric discrete model provided by an embodiment of the present invention; Figure 7 A schematic structural block diagram of a radionuclide diffusion path simulation device provided in an embodiment of the present invention.

[0018] Icons: 100-electronic device; 101-memory; 102-communication interface; 103-processor; 104-communication bus; 300-radioactive nuclide diffusion path simulation device; 310-matrix confirmation unit; 320-geometric discrete model establishment unit; 330-discrete point temperature and nuclide concentration confirmation unit; 340-diffusion simulation unit. DETAILED DESCRIPTION

[0019] 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.

[0020] 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.

[0021] It should be noted that similar reference numerals and letters represent similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined or explained in subsequent drawings. At the same time, in the description of the present invention, the terms "first", "second", etc. are used only to distinguish the description and should not be understood as indicating or implying relative importance.

[0022] It should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply the existence of any such actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or device comprising the element.

[0023] As described in the background art, traditional diffusion models mainly establish a concentration gradient-driven diffusion equation based on Fick's law, without considering the effect of temperature gradient on diffusion behavior.

[0024] Based on the above situation, the embodiments of the present invention provide a method, device, electronic device and storage medium for simulating the diffusion path of radioactive nuclides. By introducing key parameters such as matrix temperature distribution, pore structure information and initial concentration distribution, high-precision prediction and simulation of the nuclide diffusion path can be achieved.

[0025] To implement the process steps and functions of each example of the present invention, please refer to Figure 1 , Figure 1 This is a schematic block diagram of an electronic device provided in an embodiment of the present invention. The electronic device 100 includes a memory 101 and a processor 103. The memory 101 and processor 103 are electrically connected to each other, directly or indirectly, to enable data transmission or interaction. For example, these components can be electrically connected to each other via one or more communication buses 104 or signal lines. The memory 101 can be used to store software programs and modules, and the processor 103 executes the software programs and modules stored in the memory 101, thereby performing various functional applications and data processing.

[0026] The electronic device 100 may be, but is not limited to, a personal computer (PC), a server, a distributed computer, or the like. It is understood that the electronic device 100 is not limited to a physical server and may also be a virtual machine on a physical server, a virtual machine built on a cloud platform, or other computer that provides the same functionality as the server or virtual machine. The operating system of the electronic device 100 may be, but is not limited to, Windows, Linux, or the like.

[0027] The memory 101 may be, but is not limited to, a random access memory (RAM), a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), etc.

[0028] The communication connection between the electronic device 100 and an external device is achieved through at least one communication interface 102 (which can be wired or wireless).

[0029] Processor 103 may be an integrated circuit chip with signal processing capabilities. During implementation, each step of the embodiments of the present invention may be completed by hardware integrated logic circuits in processor 103 or by software instructions. Processor 103 may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0030] I understand. Figure 1 The structure shown is for illustration only. The electronic device 100 may further include Figure 1 More or fewer components than shown, or with Figure 1 Different configurations shown. Figure 1 Each component shown in the figure can be implemented by hardware, software or a combination thereof.

[0031] The following is an exemplary description of the radionuclide diffusion path simulation method provided by the present invention. Figure 2 For, see Figure 2 , the execution subject of this method can be the above Figure 1 The electronic device 100 shown in FIG. 1 includes the following steps: Figure 2 The following steps are described: S210: Obtain a matrix in which the radioactive nuclides are located, and determine the temperature distribution of the matrix, the geometric shape of the matrix, the internal pore distribution of the matrix, the medium type of the internal pores, and the initial concentration distribution of the radioactive nuclides.

[0032] For example, there may be multiple radionuclides, such as Cs and Se. The diffusion of each radionuclide is independent. Therefore, in the embodiment of the present invention, only the diffusion of one radionuclide is described.

[0033] S220: Establishing a geometric discrete model of the matrix based on the geometric shape and internal pore distribution, wherein the geometric discrete model includes a plurality of discrete points, which are divided into matrix discrete points located on the matrix and pore discrete points located in the internal pores.

[0034] S230: Determine the temperature and nuclide concentration of each discrete point based on the temperature distribution, the initial concentration distribution, and the geometric discrete model.

[0035] S240: Based on the temperature and nuclide concentration of each discrete point in the matrix, multiple diffusion simulations are performed to construct the diffusion path of the radioactive nuclide in the matrix.

[0036] First, the matrix in which the radionuclides reside is determined. For example, if the matrix material is granite, this method simulates the diffusion path within the matrix. The temperature distribution of the matrix, the matrix geometry, the internal pore distribution, the type of medium in the internal pores, and the initial concentration distribution of the radionuclides in the matrix are then determined.

[0037] For example, see Figure 3 , a cross section of the matrix is ​​rectangular, in which there are some pores. The geometric discrete model of the matrix is ​​established, see Figure 4 The geometric discrete model is represented by discrete points, wherein the dark discrete points indicate that the discrete points are located on the matrix material, and the light discrete points indicate that the discrete points are located on the internal pores.

[0038] After establishing the geometric discrete model, the temperature of each discrete point can be determined based on the temperature distribution of the matrix, and the nuclide concentration at each discrete point can be determined based on the initial nuclide concentration distribution. Determining the temperature and nuclide concentration at discrete points can be achieved using numerical methods such as the finite element method, the finite volume method, or the finite difference method.

[0039] Finally, after establishing the geometric discrete model and determining the temperature and nuclide concentration of each discrete point, multiple diffusion simulations can be performed based on the temperature and nuclide concentration of each discrete point in the matrix to obtain the diffusion path of each discrete point in the matrix for each diffusion. Based on the diffusion paths of each diffusion, the complete diffusion path of the radioactive nuclide in the matrix can be constructed.

[0040] This method establishes a geometric discrete model for the matrix in which the radioactive nuclides are located, and determines the temperature and nuclide concentration of each discrete point in the geometric discrete model based on the temperature distribution of the matrix and the initial concentration distribution of the nuclides. Based on the temperature and nuclide concentration of each discrete point, a high-precision prediction simulation of the nuclide diffusion path is achieved.

[0041] In one possible implementation, see Figure 5 , step S240 may include the following steps: S241: For each diffusion simulation, the diffusion vector and diffusion flux of each discrete point in the matrix are calculated based on the temperature and nuclide concentration of each discrete point in the matrix.

[0042] S242: Based on each diffusion flux, the nuclide concentration of each discrete point of the matrix after the current diffusion simulation is updated.

[0043] S243: Constructing a diffusion path of the radioactive nuclide in the matrix based on each diffusion vector of each diffusion simulation.

[0044] For each diffusion simulation, a study is conducted. In a diffusion simulation, each discrete point is calculated to diffuse to other discrete points, although the diffusion flux may be 0, because each point may also receive the nuclide concentration of other discrete points. For example, if the initial nuclide concentration of a discrete point is 0, it will not diffuse outward for the first time and receive the diffusion flux of other points. In the next diffusion simulation, it has a nuclide concentration and may diffuse outward. Therefore, we calculate each discrete point of the matrix in each diffusion mode.

[0045] After a diffusion simulation, each matrix discrete point has a diffusion vector and diffusion flux. Of course, the diffusion vector and diffusion flux can be 0. Figure 6 The diffusion vector and diffusion flux represent how much nuclide concentration diffuses from one discrete point to another. Therefore, after obtaining a diffusion vector and diffusion flux for each matrix discrete point, each matrix discrete point may have both outward-diffused nuclide concentrations and inward-diffused nuclide concentrations from other discrete points. The updated nuclide concentration at each matrix discrete point is equal to the original nuclide concentration minus the outward-diffused nuclide concentration and the added inward-diffused nuclide concentration.

[0046] Finally, after each diffusion simulation, the nuclide concentration of each discrete point of the matrix is ​​updated, and the diffusion simulation is performed again. This is repeated a preset number of times, and the diffusion path of the radioactive nuclide in the matrix can be constructed based on the diffusion vectors of each diffusion simulation.

[0047] Each diffusion simulation may select discrete points within a preset range to calculate the diffusion direction of each discrete point of the substrate. The above step S241 may include the following steps: S2411: For each matrix discrete point, with the matrix discrete point as the target matrix discrete point, find the first discrete point with the highest temperature, the second discrete point with the lowest temperature, the third discrete point with the highest nuclide concentration, and the fourth discrete point with the lowest nuclide concentration among the matrix discrete points within a circle with the target matrix discrete point as the center and a preset radius.

[0048] S2412: Calculate the temperature gradient vector of the target substrate discrete point according to the temperature of the target substrate discrete point, the temperature of the first discrete point, and the temperature of the second discrete point.

[0049] S2413: Calculate the nuclide concentration gradient vector of the target matrix discrete point based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the third discrete point, and the nuclide concentration of the fourth discrete point.

[0050] S2414: Calculate the cross product of the temperature gradient vector and the nuclide concentration gradient vector to obtain the diffusion vector of the target matrix discrete point in the matrix.

[0051] S2415: The discrete point pointed to by the diffusion vector of the target matrix discrete point is used as the diffusion discrete point of the target matrix discrete point.

[0052] S2416: Calculate the diffusion flux from the target matrix discrete point to the diffusion discrete point based on the nuclide concentration at the target matrix discrete point, the nuclide concentration at the diffusion discrete point, the diffusion coefficient of the radioactive nuclide in the matrix, and the adsorption coefficient of the matrix to the radioactive nuclide.

[0053] For each matrix discrete point, first identify the discrete point with the highest temperature, the lowest temperature, the highest nuclide concentration, and the lowest nuclide concentration within a certain range of the matrix discrete point. Based on the highest and lowest temperature discrete points, a temperature gradient vector is calculated. This temperature gradient vector points from the highest temperature discrete point to the lowest temperature discrete point, and its value is the value of the highest temperature discrete point minus the lowest temperature discrete point. Similarly, the nuclide concentration gradient vector can be calculated. The temperature gradient vector is then multiplied by the nuclide concentration gradient vector to obtain a new vector, which is the diffusion vector of the target matrix discrete point in the matrix.

[0054] The point pointed by the diffusion vector from the target matrix discrete point is then used as the diffusion discrete point, that is, the target matrix discrete point will diffuse towards this diffusion discrete point. The diffusion flux from the target matrix discrete point to the diffusion discrete point is then calculated based on the nuclide concentrations at these two points, the diffusion coefficient of the radionuclide in the matrix, and the matrix's adsorption coefficient of the radionuclide.

[0055] Step S2412 may include: S24121: Determine whether the temperature of the target substrate discrete point is greater than or equal to the temperature of the first discrete point.

[0056] S24122: If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the target substrate discrete point minus the temperature of the second discrete point.

[0057] S24123: If not, determine whether the temperature of the target substrate discrete point is equal to the temperature of the second discrete point.

[0058] S24124: If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the target substrate discrete point.

[0059] S24125: If not, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the second discrete point.

[0060] For example, the temperature of the target matrix discrete point is T i , the temperature of the first discrete point is T max , the temperature of the second discrete point is T min , first determine T i Is it greater than or equal to T max If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the target substrate discrete point minus the temperature of the second discrete point, that is, the direction of the temperature gradient vector of the target substrate discrete point is from the target substrate discrete point to the second discrete point, and the value is T i -T min If not, then judge T i Is it equal to T min If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the target substrate discrete point; if not, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the second discrete point.

[0061] S2413 may include: S24131: Determine whether the nuclide concentration of the target matrix discrete point is greater than or equal to the nuclide concentration of the third discrete point.

[0062] S24132: If so, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the target matrix discrete point minus the nuclide concentration of the fourth discrete point.

[0063] S24133: If not, determine whether the nuclide concentration of the target matrix discrete point is equal to the nuclide concentration of the fourth discrete point.

[0064] S24134: If so, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the third discrete point minus the nuclide concentration of the target matrix discrete point.

[0065] S24135: If not, the nuclide concentration gradient vector of the target matrix discrete point is equal to the vector of the nuclide concentration at the third discrete point minus the nuclide concentration at the fourth discrete point.

[0066] Steps S24131 to S24135 are logically identical to the above-mentioned steps S24121 to S24125 and will not be repeated here.

[0067] After determining the diffusion discrete point in step S2415, it is also possible to determine whether the discrete point is a pore discrete point. If so, the diffusion discrete point can be corrected. In a possible implementation, after step S2415, the method can further include: S24151: Determine whether the diffusion discrete point belongs to the pore discrete point.

[0068] S24152: If yes, start searching outward along the diffusion vector from the diffusion discrete point until the first discrete point belonging to the matrix discrete point is found, and change the diffusion discrete point to the first discrete point.

[0069] S24153: Change the diffusion vector so that the target matrix discrete point points to the first discrete point.

[0070] Determine whether the diffusion discrete point belongs to a pore discrete point. If not, no processing is required. If so, start searching outward along the diffusion vector from the diffusion discrete point and use the first matrix discrete point found as the corrected diffusion discrete point.

[0071] There are many ways to calculate the diffusion flux. In one possible implementation, step S2416 may include: S24161: Determine whether the nuclide concentration at the diffusion discrete point is 0.

[0072] S24162: If so, calculate the estimated concentration of the diffusion discrete point using the diffusion coefficient, and calculate the diffusion flux from the target matrix discrete point to the diffusion discrete point using the adsorption coefficient.

[0073] S24163: If not, calculate the diffusion flux from the target matrix discrete point to the diffusion discrete point using the adsorption coefficient.

[0074] First, determine whether the nuclide concentration at the diffusion discrete point is 0. If it is 0, first calculate the estimated concentration of the diffusion discrete point using the diffusion coefficient using the following formula: Cb = Da×Ca / Db, where Cb is the estimated concentration at the diffusion discrete point, Ca is the nuclide concentration at the target matrix discrete point, Da and Db are diffusion coefficients, and D=D0exp(-Q / RT), where D0 is the frequency factor, Q is the diffusion activation energy (for example, the diffusion coefficient of Cs in granite is D0=8.48×10-13m2 / s, and the diffusion activation energy is 189kJ / mol; the diffusion coefficient of Se in granite is D0=9.75×10-13m2 / s, and the diffusion activation energy is 303kJ / mol), R is the gas constant, and T is the absolute temperature of the corresponding discrete point (i.e., the temperature of the discrete point plus 273.15).

[0075] The diffusion flux from the target matrix discrete point to the diffusion discrete point is then calculated using the following formula: (1-Kd) × (Ca-Cb), where Kd is the distribution coefficient. For example, the adsorption distribution coefficient of Cs in granite is 562 mg / L, and the adsorption distribution coefficient of Se in granite is 11.88 mg / L.

[0076] In addition, if the diffusion discrete points are corrected in step S24152 and the medium type of the internal pores is liquid, the diffusion flux needs to be multiplied by a correction coefficient, for example, 0.9.

[0077] Furthermore, the present invention also provides a radionuclide diffusion path simulation device, see Figure 7 The radionuclide diffusion path simulation device 300 includes: The matrix confirmation unit 310 is used to obtain the matrix in which the radioactive nuclides are located, and determine the temperature distribution of the matrix, the geometric shape of the matrix, the internal pore distribution of the matrix, the medium type of the internal pores, and the initial concentration distribution of the radioactive nuclides.

[0078] The geometric discrete model building unit 320 is used to build a geometric discrete model of the matrix according to the geometric shape and internal pore distribution; wherein the geometric discrete model includes multiple discrete points, which are divided into matrix discrete points located on the matrix and pore discrete points located in the internal pores.

[0079] The discrete point temperature and nuclide concentration confirmation unit 330 is used to determine the temperature and nuclide concentration of each discrete point according to the temperature distribution, the initial concentration distribution and the geometric discrete model.

[0080] The diffusion simulation unit 340 is used to perform multiple diffusion simulations based on the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix.

[0081] In summary, the radionuclide diffusion path simulation method, device, electronic device, and storage medium provided by the present invention achieve high-precision prediction and simulation of radionuclide diffusion paths by incorporating key parameters such as substrate temperature distribution, pore structure information, and initial concentration distribution. This not only improves the accuracy and applicability of simulation results but also plays an important role in fields such as nuclear waste disposal and environmental remediation.

[0082] In the embodiments provided by the present invention, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are merely illustrative. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architectures, functions, and operations of the devices, methods, and computer program products according to multiple embodiments of the present invention. In this regard, each box in the flowchart or block diagram can represent a module, program segment, or part of the code, which contains one or more executable instructions for implementing the specified logical functions. It should also be noted that in some alternative implementations, the functions marked in the boxes can also occur in an order different from that marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, as well as the combination of boxes in the block diagram and / or flowchart, can be implemented using a dedicated hardware-based system that performs the specified functions or actions, or can be implemented using a combination of dedicated hardware and computer instructions.

[0083] In addition, the functional modules in the various embodiments of the present invention may be integrated together to form an independent part, or each module may exist independently, or two or more modules may be integrated to form an independent part.

[0084] If a function is implemented as a software module and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or the portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a computer-readable storage medium and includes several instructions for causing a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the various embodiments of the method of the present invention. The aforementioned storage medium includes various media that can store program code, such as USB flash drives, mobile hard drives, read-only memories (ROMs), random access memories (RAMs), magnetic disks, or optical disks.

[0085] The above are merely preferred embodiments of the present invention and are not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention shall be included within the scope of protection of the present invention.

[0086] It will be apparent to those skilled in the art that the present invention is not limited to the details of the exemplary embodiments described above and that the invention can be embodied in other specific forms without departing from the spirit or essential characteristics of the invention. Therefore, the embodiments should be considered in all respects as illustrative and non-restrictive, and the scope of the invention is defined by the appended claims, not the foregoing description, and all variations within the meaning and range of equivalents of the claims are intended to be included therein. Any reference sign in a claim should not be construed as limiting the claim to which it relates.

Claims

1. A method for simulating the diffusion path of radionuclides, characterized in that: include: Acquiring a matrix in which the radionuclide is located, and determining the temperature distribution of the matrix, the geometric shape of the matrix, the internal pore distribution of the matrix, the medium type of the internal pores, and the initial concentration distribution of the radionuclide; Establishing a geometric discrete model of the substrate based on the geometric shape and the internal pore distribution; wherein the geometric discrete model includes a plurality of discrete points, and the discrete points are divided into substrate discrete points located on the substrate and pore discrete points located in the internal pores; Determining the temperature and nuclide concentration of each discrete point according to the temperature distribution, the initial concentration distribution, and the geometric discrete model; Multiple diffusion simulations are performed based on the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix.

2. The method according to claim 1, characterized in that The method of performing multiple diffusion simulations based on the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix includes: For each diffusion simulation, the diffusion vector and diffusion flux of each discrete point in the matrix are calculated according to the temperature and nuclide concentration of each discrete point in the matrix; updating the nuclide concentration of each discrete point of the matrix after the current diffusion simulation according to each diffusion flux; The diffusion path of the radioactive nuclide in the matrix is ​​constructed according to the diffusion vectors of each diffusion simulation.

3. The method according to claim 2, characterized in that Calculating the diffusion vector and diffusion flux of each discrete point in the matrix according to the temperature and nuclide concentration of each discrete point in the matrix includes: For each of the matrix discrete points, taking the matrix discrete point as a target matrix discrete point, finding a first discrete point with the highest temperature, a second discrete point with the lowest temperature, a third discrete point with the highest nuclide concentration, and a fourth discrete point with the lowest nuclide concentration among the matrix discrete points within a circle with the target matrix discrete point as the center and a preset radius; Calculating a temperature gradient vector of the target substrate discrete point according to the temperature of the target substrate discrete point, the temperature of the first discrete point, and the temperature of the second discrete point; Calculating a nuclide concentration gradient vector of the target matrix discrete point based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the third discrete point, and the nuclide concentration of the fourth discrete point; Calculating the cross product of the temperature gradient vector and the nuclide concentration gradient vector to obtain a diffusion vector of the target matrix discrete point in the matrix; Taking the discrete point pointed to by the diffusion vector as the diffusion discrete point of the target matrix discrete point; The diffusion flux from the target matrix discrete point to the diffusion discrete point is calculated based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the diffusion discrete point, the diffusion coefficient of the radioactive nuclide in the matrix, and the adsorption coefficient of the matrix to the radioactive nuclide.

4. The method according to claim 3, characterized in that The step of calculating the temperature gradient vector of the target substrate discrete point according to the temperature of the target substrate discrete point, the temperature of the first discrete point, and the temperature of the second discrete point includes: Determining whether the temperature of the target substrate discrete point is greater than or equal to the temperature of the first discrete point; If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the target substrate discrete point minus the temperature of the second discrete point; If not, determining whether the temperature of the target substrate discrete point is equal to the temperature of the second discrete point; If so, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the target substrate discrete point; If not, the temperature gradient vector of the target substrate discrete point is equal to the vector of the temperature of the first discrete point minus the temperature of the second discrete point.

5. The method according to claim 3, characterized in that Calculating the nuclide concentration gradient vector of the target matrix discrete point according to the nuclide concentration of the target matrix discrete point, the nuclide concentration of the third discrete point, and the nuclide concentration of the fourth discrete point includes: Determining whether the nuclide concentration of the target matrix discrete point is greater than or equal to the nuclide concentration of the third discrete point; If so, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the target matrix discrete point minus the nuclide concentration of the fourth discrete point; If not, determining whether the nuclide concentration of the target matrix discrete point is equal to the nuclide concentration of the fourth discrete point; If so, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the third discrete point minus the nuclide concentration of the target matrix discrete point; If not, the nuclide concentration gradient vector of the target matrix discrete point is equal to the nuclide concentration of the third discrete point minus the nuclide concentration of the fourth discrete point.

6. The method according to claim 3, characterized in that After the step of using the discrete point pointed to by the diffusion vector as the diffusion discrete point of the target matrix discrete point, the method further includes: Determining whether the diffusion discrete point belongs to the pore discrete point; If so, searching outwards from the diffusion discrete point along the diffusion vector until a first discrete point belonging to the matrix discrete point is found, and changing the diffusion discrete point to the first discrete point; The diffusion vector is changed so that the target matrix discrete point points to the first discrete point.

7. The method according to claim 3, characterized in that Calculating the diffusion flux from the target matrix discrete point to the diffusion discrete point based on the nuclide concentration of the target matrix discrete point, the nuclide concentration of the diffusion discrete point, the diffusion coefficient of the radioactive nuclide in the matrix, and the adsorption coefficient of the matrix to the radioactive nuclide includes: Determining whether the nuclide concentration of the diffusion discrete point is 0; If so, calculating the estimated concentration of the diffusion discrete point by using the diffusion coefficient, and calculating the diffusion flux from the target matrix discrete point to the diffusion discrete point by using the adsorption coefficient; If not, the diffusion flux from the target matrix discrete point to the diffusion discrete point is calculated using the adsorption coefficient.

8. A radionuclide diffusion path simulation device, characterized in that: include: a matrix confirmation unit, configured to obtain a matrix in which the radionuclide is located, and determine a temperature distribution of the matrix, a geometric shape of the matrix, an internal pore distribution of the matrix, a medium type of the internal pores, and an initial concentration distribution of the radionuclide; a geometric discrete model establishing unit, configured to establish a geometric discrete model of the substrate according to the geometric shape and the internal pore distribution; wherein the geometric discrete model includes a plurality of discrete points, the discrete points being divided into substrate discrete points located on the substrate and pore discrete points located in the internal pores; A discrete point temperature and nuclide concentration confirmation unit, configured to determine the temperature and nuclide concentration of each discrete point according to the temperature distribution, the initial concentration distribution, and the geometric discrete model; The diffusion simulation unit is used to perform multiple diffusion simulations according to the temperature and nuclide concentration of each discrete point of the matrix to construct the diffusion path of the radioactive nuclide in the matrix.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images