Dose evaluation method for radioactive waste disposal site based on biosphere evolution
By employing a dose assessment method for radioactive waste disposal sites based on biosphere evolution, the problem of dose assessment instability caused by dynamic environmental evolution has been solved, enabling more accurate long-term safety assessments and providing a reliable risk spectrum basis.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-12-22
- Publication Date
- 2026-04-07
AI Technical Summary
Existing technologies for dose assessment at radioactive waste disposal sites neglect the dynamic evolution of the environment, leading to unromantic and unreliable dose assessment results, especially the lack of correlation between nuclide migration calculations and biosphere identification on long-term timescales.
A dose assessment method for radioactive waste disposal sites based on biosphere evolution is adopted. By simulating environmental parameters through geographic information systems and climate evolution data, the location of key groups and exposure pathways are dynamically identified. Combined with a radionuclide migration model, radionuclide concentrations and public doses are calculated, based on dynamic environmental changes.
It improves the authenticity and credibility of dose assessment, provides a more reliable basis for safety decisions, and transforms qualitative environmental uncertainty into quantitative dose result range by simulating the biosphere evolution path under different climate scenarios.
Smart Images
Figure CN121806084A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of radionuclide migration assessment technology, and in particular to a dose assessment method for radioactive waste disposal sites based on biosphere evolution. Background Technology
[0002] The safe disposal of radioactive waste, especially high-level waste and long-lived intermediate-level waste, is a critical issue concerning environmental safety and public health. Disposal systems must ensure effective isolation of radionuclides from the biosphere over timescales of tens of thousands of years or even longer. Therefore, assessing the long-term safety performance of disposal sites, particularly predicting potential future radiation doses to the public, is a core basis for site selection, design, safety approval, and post-closure oversight of disposal facilities.
[0003] Currently, long-term safety assessments of radioactive waste disposal both domestically and internationally typically employ the following dose assessment methods: First, based on the current climate, hydrology, geology, and land use conditions of the disposal site, one or more representative, static biospheres are identified, and corresponding public exposure pathways (such as drinking water and ingestion of agricultural products) are determined. Second, simplified nuclide migration models (such as convection-diffusion models or reservoir models using fixed parameters) are used to simulate the process of nuclides migrating from the disposal engineering barrier to the biosphere via the geological medium. Finally, by combining biosphere parameters and dose conversion factors, the potential public exposure dose is calculated.
[0004] However, this dose assessment method ignores the dynamic evolution that inevitably occurs in the natural environment, neglects the time scale of thousands of years, and the fact that nuclide migration calculations and biosphere identification are relatively independent and lack correlation. As a result, the dose assessment and the final calculation results may be based on an incorrect "spatiotemporal matching" relationship, thus affecting the robustness and credibility of the safety assessment conclusions. Summary of the Invention
[0005] To address the shortcomings of existing technologies, this application aims to provide a dose assessment method for radioactive waste disposal sites based on biosphere evolution, which links dose assessment with dynamic environmental changes and significantly improves the authenticity and reliability of assessment results.
[0006] To achieve the above objectives, this application provides a dose assessment method for radioactive waste disposal sites based on biosphere evolution, comprising: Using geographic information systems and regional climate evolution data, environmental parameters of the environment surrounding the treatment site at future time scales are simulated and generated. Based on environmental parameters of the surrounding environment of the treatment site at future time scales, watershed tools are used to dynamically identify the locations of key clusters and sets of irradiated pathways around the treatment site at different future times. Based on the set of irradiation pathways around the treatment site at different future times, the source term data of the treatment site, and the address medium parameters, the nuclide concentration in the near-field aqueous phase is calculated; The corresponding public dose is calculated based on the environmental parameters corresponding to the selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in the near-field water phase.
[0007] Furthermore, the specific steps for calculating and determining the corresponding nuclide concentration in the near-field aqueous phase based on source term data of the disposal site at different future times include: Based on the source term data of the treatment site at different future times around the treatment site, we filter from the FEPs database and construct the interaction matrix of the treatment site at different future times according to the corresponding set of irradiated pathways. Based on the interaction matrix, the set of irradiated pathways, and the address medium parameters, a corresponding nuclide migration model is established using the cell method to calculate the nuclide concentration in the near-field aqueous phase.
[0008] Furthermore, the environmental parameters include at least: precipitation, evaporation, water flow velocity in water fissures, dispersion coefficient, flow wetting surface area, matrix diffusion coefficient, effective surface density of soil, river irrigation rate, vegetation cover type, and surface runoff path.
[0009] Furthermore, the set of irradiated paths includes at least: water source type, average irrigation rate, irrigation water source type, main crop type, and food intake.
[0010] Furthermore, the nuclide migration model uses the following formula: ; ; ; ; ; ; ; in, The contact area between the chamber and the rock. The flow porosity in the rock. The diffusion coefficient of the nuclide in water. is the residence time of the pore water in the rock in contact with the reservoir; Length is the length of the contact between the flowing water and the reservoir; Darcy Velocity is the Darcy velocity of the water in the rock. Let t be the concentration of nuclides in the near-field aqueous phase of the treatment reservoir; in, The concentration of radionuclide n in the fissure. The concentration of the parent radionuclide p in the fissure. The concentration of radionuclide n in the matrix pore water. This represents the concentration of the parent radionuclide p in the matrix pore water. Let be the flux of the nuclide entering the far-field fracture at time t. For time, The velocity of the water flow in the water-bearing fracture. The dispersion coefficient is... For the flow wetting surface area, The matrix diffusion coefficient is... The distance along the flow direction. To the depth of penetration into the matrix, Let n be the decay rate of nuclide n. It is the retardation factor for radionuclide n in the matrix. It is the retardation factor of the parent radionuclide p in the matrix.
[0011] Furthermore, the specific steps for calculating the corresponding public dose based on the environmental parameters corresponding to the selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in the near-field aqueous phase adopt the following formula: ; in, This is the ingested dose; This refers to the concentration of radionuclides in ingested products / drinking water; Ingested dose conversion factor; This refers to the amount of product ingested / water consumed.
[0012] Furthermore, the specific steps for calculating the corresponding public dose based on the environmental parameters corresponding to the selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in the near-field aqueous phase also employ the following formula: ; ; ; in, It refers to the concentration of sediment deposited on the soil surface due to irrigation water; It is a concentration factor for radionuclides in plant products; It refers to the concentration of radionuclides in the river water; It is the average irrigation rate during the growing season; P is the effective surface density of the soil. It is the accumulation time of radioactive nuclides on the soil surface; It is the share of river water for irrigation. It is the effective rate constant for the removal of radionuclides from the soil surface. It is the rate constant for the elimination of other nuclides besides decay; It refers to the half-life of a nuclide; t h It refers to the time from harvesting to consumption of agricultural products.
[0013] Furthermore, the specific steps for calculating the corresponding public dose based on the environmental parameters corresponding to the selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in the near-field aqueous phase also employ the following formula: ; in, It refers to the proportion of radioactive nuclides in animal products; It refers to the time from the slaughter of animals to the consumption of their products by humans. This refers to the daily feed intake of animals; It is the specific activity of radionuclides in animal feed.
[0014] Furthermore, the method also includes: At each time step on the time scale, environmental parameters, key group locations, and exposure pathway sets are automatically updated to calculate public doses.
[0015] The dose assessment method for radioactive waste disposal sites based on biosphere evolution proposed in this application uses a "variable interaction matrix" to dynamically link the state of the biosphere at different future periods with the migration paths of radionuclides, thus avoiding the disconnect between migration calculations and exposure scenarios in traditional methods.
[0016] This application presents a dose assessment method for radioactive waste disposal sites based on biosphere evolution. By coupling climate evolution prediction with biosphere dynamics identification, dose assessment on a millennium-scale is established on the scientific basis of environmental dynamics, significantly improving the accuracy and reliability of the prediction results.
[0017] The dose assessment method for radioactive waste disposal sites based on biosphere evolution proposed in this application can simulate evolution paths under different climate scenarios, transforming long-term, qualitative environmental uncertainties into quantitative dose result ranges, and providing a more reliable risk spectrum basis for safety decisions.
[0018] Other features and advantages of this application will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing this application. Attached Figure Description
[0019] The accompanying drawings are provided to further illustrate the present application and form part of the specification. Together with the embodiments of the present application, they serve to explain the present application but do not constitute a limitation thereof. In the drawings: Figure 1This is a schematic flowchart of the dose assessment method for radioactive waste disposal sites based on biosphere evolution, as described in Embodiment 1 of this application. Figure 2 This is a schematic diagram for identifying a watershed. Figure 3 This is a schematic diagram of the interaction matrix; Figure 4 This is a schematic diagram of the effective dose curve for the public over the next 100,000 years for a treatment site. Detailed Implementation
[0020] Embodiments of this application will now be described in more detail with reference to the accompanying drawings. While some embodiments of this application are shown in the drawings, it should be understood that this application can be implemented in various forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided to provide a more thorough and complete understanding of this application. It should be understood that the drawings and embodiments of this application are for illustrative purposes only and are not intended to limit the scope of protection of this application.
[0021] It should be understood that the steps described in the method embodiments of this application may be performed in different orders and / or in parallel. Furthermore, the method embodiments may include additional steps and / or omit the steps shown. The scope of this application is not limited in this respect.
[0022] The term "comprising" and its variations as used herein are open-ended inclusions, meaning "including but not limited to". The term "based on" means "at least partially based on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment"; the term "some embodiments" means "at least some embodiments". Definitions of other terms will be given in the description below.
[0023] It should be noted that the terms "one" and "multiple" used in this application are illustrative rather than restrictive, and those skilled in the art should understand that, unless explicitly stated otherwise in the context, they should be understood as "one or more". "Multiple" should be understood as two or more.
[0024] The embodiments of this application will now be described in detail with reference to the accompanying drawings.
[0025] Example 1 One embodiment of this application provides a dose assessment method for radioactive waste disposal sites based on biosphere evolution, which will be described below. Figure 1-4 The dose assessment method for radioactive waste disposal sites based on biosphere evolution, as described in this application, is as follows: Step S101: Using geographic information systems and regional climate evolution data, simulate and generate environmental parameters of the environment surrounding the treatment site on future time scales; In this embodiment of the application, an overlay analysis is performed using a geographic information system. Based on regional climate evolution data (temperature and precipitation prediction data) under different carbon emission scenarios and a digital elevation model (DEM) around the disposal site, environmental parameters that change over time are simulated and generated. In the embodiments of this application, regional climate evolution data can be obtained from regional climate evolution models such as the GCM model.
[0026] In the embodiments of this application, regional climate evolution data include, but are not limited to: precipitation, evaporation, river irrigation rate, vegetation cover type, and surface runoff path.
[0027] In the embodiments of this application, environmental parameters include, but are not limited to: precipitation, evaporation, water flow velocity in water fissures, dispersion coefficient, flow wetting surface area, matrix diffusion coefficient, effective surface density of soil, river irrigation rate, vegetation cover type, and surface runoff path.
[0028] Understandably, this step involves a quantitative analysis of the long-term evolution of the biosphere, which refers to the near-surface environment directly related to radionuclide migration and public exposure to radiation, including the near-surface atmospheric layer, surface hydrosphere, soil sphere, and biological communities.
[0029] For example, climate simulations predict that 1,000 years after treatment, the region will evolve into a "warm and humid forest" state, with output parameters of 800 mm annual precipitation and tree-dominated vegetation; and 10,000 years after treatment, the region will evolve into a "dry and cold grassland" state, with output parameters of 300 mm annual precipitation and shrub-dominated vegetation. These environmental parameters will directly replace the relevant variables in the dosage formula.
[0030] Step S102: Using watershed tools, dynamically identify the locations of key clusters and sets of irradiated pathways around the treatment site at different future time scales based on environmental parameters of the environment surrounding the treatment site. In the embodiments of this application, based on the environmental parameters of the surrounding environment of the treatment site at future time scales and DEM (digital elevation model) data, the watershed tool is used to dynamically identify the locations of catchment areas, lakes and groundwater outcrops at different future times, determine habitable areas with agricultural conditions, and finally determine the location coordinates of key groups and the set of irradiation pathways that change over time.
[0031] In the embodiments of this application, the key group location coordinates refer to the possible water source areas around the future treatment site identified through evolution, and then determined to be the locations where people may be located.
[0032] In the embodiments of this application, the set of irradiated pathways includes, but is not limited to: drinking water source type, average irrigation rate, irrigation water source type, main crop type, food intake, etc. For example, Figure 2 This is a watershed identification diagram. It uses a digital elevation model as a background and uses contour lines or color gradients to present the topographic undulations of the treatment area (such as mountains, valleys, and basins), mark the water flow path, and divide the confluence areas of different water systems.
[0033] In this embodiment, the watershed tool in GIS is used to calculate the confluence path, and the surface runoff is calculated by combining the precipitation parameter from the environmental parameters. The DEM provides the watershed geometric boundary, and the precipitation parameter provides the input flux. Based on both, the dynamic value of river / lake water volume can be calculated. Based on water volume Determine whether human habitation is met (e.g., by setting a minimum water availability threshold) to determine public residence locations; simultaneously, calculate... It is directly used as a dilution factor for the concentration of radionuclides in water bodies in subsequent dose calculation formulas. For example, at T=5000, a DEM shows a potential lake basin in a low-lying area. If the climate data predicts abundant rainfall at this time, GIS determines that a "lake-type biosphere" has formed in this area, and the public living near the lake will be exposed to radiation through "consuming lake fish"; if the climate data predicts drought, then this area is only a "dry valley", and the public may migrate to downstream rivers, and the exposure pathway becomes "drinking river water".
[0034] Step S103: Calculate the nuclide concentration in the near-field aqueous phase based on the set of irradiation pathways around the treatment site at different future times, the source term data of the treatment site, and the address medium parameters; First, based on the source term data of the treatment site around the treatment site at different future times, we filter from the FEPs database and construct the interaction matrix of the treatment site around the treatment site at different future times according to the corresponding set of irradiated pathways. The FEPs database is a resource library that filters specific elements for a site from the Feature, Event, and Process (FEPs) database using source item data from the disposal site. The interaction matrix uses these elements to construct the system topology (i.e., the path of nuclides from medium A to medium B).
[0035] Based on the set of irradiated pathways, different interaction matrices are constructed for different time windows. For example, during the wet season, the interaction matrix contains a strong interaction pathway of "aquifer - surface water - sediment"; during the dry season, this pathway may break and transform into a pathway of "deep groundwater - soil (through capillary action)".
[0036] Then, based on the interaction matrix, the set of irradiated pathways, and the address medium parameters, the corresponding nuclide migration model is established using the cell method to calculate the nuclide concentration in the near-field aqueous phase.
[0037] In this embodiment of the application, the nuclide migration model adopts the following formula: ; ; ; ; ; ; in, The contact area (m2) between the chamber and the rock. The flow porosity in the rock. denoted as the diffusion coefficient of the nuclide in water (m² / a). , where is the residence time (a) of the pore water in the rock in contact with the reservoir; Length is the length (m) of the contact between the flowing water and the reservoir; and Darcy Velocity is the Darcy velocity of the water in the rock (m / a). in, The concentration of radionuclide n in the fissure. The concentration of the parent radionuclide p in the fissure. The concentration of radionuclide n in the matrix pore water. This represents the concentration of the parent radionuclide p in the matrix pore water. Let be the flux of the nuclide entering the far-field fracture at time t. For time, The velocity of the water flow in the water-bearing fracture. The dispersion coefficient is... For the flow wetting surface area, The matrix diffusion coefficient is... The distance along the flow direction. To the depth of penetration into the matrix, Let n be the decay rate of nuclide n. It is the retardation factor for radionuclide n in the matrix. It is the retardation factor of the parent radionuclide p in the matrix.
[0038] In the embodiments of this application, as shown in the above formula, when considering the process of radionuclides diffusing into bedrock, the equivalent water flow rate is used to represent the diffusion transport of radionuclides to the far field, that is, the outflow from the reservoir in direct contact with the far field. This flow transports dissolved radionuclides with a concentration equal to the concentration released into the far field reservoir. The value of the equivalent water flow rate depends on the contact area size, water flux (Darcy velocity), flow porosity, and diffusion coefficient.
[0039] In obtaining Then, it was applied to the inlet boundary condition (z=0) of the far-field nuclide migration model to calculate the flux of nuclides entering the far-field fracture. The boundary conditions are as follows: ; In the embodiments of this application, at the boundary between the crack and the matrix, i.e. At that time, the concentration in the fissure is the same as the concentration in the matrix, that is: By setting the sum of advection and diffusion flow on the fracture cross section to equal the near-field release at the inflow boundary, a model for the release of radionuclides from the near field to the far field is established. The boundary conditions inside the transfer matrix limit the radionuclide flow driven solely by diffusion to zero, while at the interface with the biosphere, the radionuclide flow exiting the fracture is limited to advection.
[0040] It should be noted that the specific values of the above parameters can be obtained through field measurements or laboratory measurements.
[0041] Step S104: Calculate the corresponding public dose based on the environmental parameters corresponding to the selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in the near-field water phase; In this application's implementation, the dosage resulting from public ingestion of agricultural products irrigated by the watershed's water and from drinking water is considered: in, It is the ingested dose (Sv); It is the concentration of radionuclides in ingested products / drinking water (Bq / kg or Bq / L). Ingested dose conversion factor (Sv / Bq); Product intake / water intake (kg or L).
[0042] Among them, plant product concentration : in, It is the concentration of soil surface deposits caused by irrigation water (Bq / m³). 2 ); Concentration factor (Bqkg) of radionuclides in plant products -1 (Fresh) / Bqkg -1 (Dry) Select the corresponding biological parameter library based on the soil type and vegetation type identified in the steps; The concentration of nuclides in the river water (Bq / m³) 3 The value is obtained by dividing the flux released by the geological sphere by the river flow, where the river flow comes from watershed and climate analysis. It is the average irrigation rate during the growing season (m 3 / (m 2 a)), derived from climate evolution data. When the evolution analysis determines an arid climate, a higher setting is used. Value; in humid climates, set a lower value. The value may be zero; P is the effective surface density of soil (kg (dry) / m³). 2 ); It is the cumulative time of radionuclides on the soil surface (a). This represents the share of irrigation water from the river. If GIS analysis determines that there are no surface rivers flowing through the settlement area but there are groundwater wells, then the value is 0 (using the groundwater model instead); if there are rivers, then the value is 1. It is the effective rate constant for the removal of radionuclides from the soil surface (a -1 ), It is the rate constant of other nuclide elimination processes besides decay (a -1 ), It is the half-life of a nuclide (a -1 ); t h It is the time from harvest to consumption of crops (a).
[0043] Among them, animal product concentration : in, It is the proportion of radionuclides in animal products (d / kg). It is the time from the slaughter of an animal to the consumption of its products by humans (a). It is the daily feed intake of animals (kg (dry weight) / d). It is the specific activity of radionuclides in animal feed (Bq / kg (dry weight)).
[0044] It should be noted that the specific values of the above parameters can be obtained through field measurements or laboratory measurements.
[0045] In the embodiments of this application, during actual simulation calculations, environmental parameters, key group locations, and radiation pathway sets can be automatically updated at each time step on the time scale to calculate public dose, thereby comprehensively generating a public annual effective dose curve for the long-term evolution process. For example, ... Figure 4 As shown, Figure 4 The effective dose curve for the public at a treatment plant during a 10-week annual inspection is shown.
[0046] The above description is merely a partial embodiment of this application and an explanation of the technical principles employed. Those skilled in the art should understand that the scope of disclosure in this application is not limited to technical solutions formed by specific combinations of the above-described technical features, but should also cover other technical solutions formed by arbitrary combinations of the above-described technical features or their equivalents without departing from the above-described concept. For example, technical solutions formed by substituting the above features with (but not limited to) technical features with similar functions disclosed in this application.
[0047] Furthermore, while the operations are described in a specific order, this should not be construed as requiring these operations to be performed in the specific order shown or in sequential order. Multitasking and parallel processing may be advantageous in certain environments. Similarly, while several specific implementation details are included in the above discussion, these should not be construed as limiting the scope of this application. Certain features described in the context of individual embodiments may also be implemented in combination in a single embodiment. Conversely, various features described in the context of a single embodiment may also be implemented individually or in any suitable sub-combination in multiple embodiments.
[0048] Although the subject matter has been described using language specific to structural features and / or methodological logic, it should be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or actions described above. Rather, the specific features and actions described above are merely illustrative examples of implementing the claims.
Claims
1. A dose assessment method for radioactive waste disposal sites based on biosphere evolution, characterized in that, include: Using geographic information systems and regional climate evolution data, environmental parameters of the environment surrounding the treatment site at future time scales are simulated and generated. Based on environmental parameters of the surrounding environment of the treatment site at future time scales, watershed tools are used to dynamically identify the locations of key clusters and sets of irradiated pathways around the treatment site at different future times. Based on the set of irradiation pathways around the treatment site at different future times, the source term data of the treatment site, and the address medium parameters, the nuclide concentration in the near-field aqueous phase is calculated; The corresponding public dose is calculated based on the environmental parameters corresponding to the selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in the near-field water phase.
2. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 1, characterized in that, The specific steps for calculating and determining the corresponding nuclide concentration in the near-field aqueous phase based on source term data of the treatment site at different future times include: Based on the source term data of the treatment site at different future times around the treatment site, we filter from the FEPs database and construct the interaction matrix of the treatment site at different future times according to the corresponding set of irradiated pathways. Based on the interaction matrix, the set of irradiated pathways, and the address medium parameters, a corresponding nuclide migration model is established using the cell method to calculate the nuclide concentration in the near-field aqueous phase.
3. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 1, characterized in that, The environmental parameters include at least: precipitation, evaporation, water flow velocity in water fissures, dispersion coefficient, flow wetting surface area, matrix diffusion coefficient, effective surface density of soil, river irrigation rate, vegetation cover type, and surface runoff path.
4. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 2, characterized in that, The set of irradiated pathways includes at least: water source type, average irrigation rate, irrigation water source type, main crop type, and food intake.
5. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 2, characterized in that, The nuclide migration model uses the following formula: ； ； ； ； ； ； ； in, The contact area between the chamber and the rock. The flow porosity in the rock. Let be the diffusion coefficient of the nuclide in water. is the residence time of the pore water in the rock in contact with the reservoir; Length is the length of the contact between the flowing water and the reservoir; Darcy Velocity is the Darcy velocity of the water in the rock. Let t be the concentration of nuclides in the near-field aqueous phase of the treatment reservoir; in, The concentration of radionuclide n in the fissure. This represents the concentration of the parent radionuclide p in the fracture. The concentration of radionuclide n in the matrix pore water. This represents the concentration of the parent radionuclide p in the matrix pore water. Let be the flux of the nuclide entering the far-field fracture at time t. For time, The velocity of water flow in the water-bearing fracture. The dispersion coefficient is... For the flow wetting surface area, The matrix diffusion coefficient is... The distance along the flow direction. To the depth of penetration into the matrix, Let n be the decay rate of nuclide n. It is the retardation factor for radionuclide n in the matrix. It is the retardation factor of the parent radionuclide p in the matrix.
6. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 1, characterized in that, The specific steps for calculating the corresponding public dose based on environmental parameters corresponding to a selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in near-field water phase are as follows: ； in, This is the ingested dose; This refers to the concentration of radionuclides in ingested products / drinking water; Ingested dose conversion factor; This refers to the amount of product ingested / water consumed.
7. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 1, characterized in that, The specific steps for calculating the corresponding public dose based on environmental parameters corresponding to a selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in near-field water also employ the following formula: ； ； ； in, It refers to the concentration of sediment deposited on the soil surface due to irrigation water; It is a concentration factor for radionuclides in plant products; It refers to the concentration of radionuclides in the river water; It is the average irrigation rate during the growing season; P is the effective surface density of the soil. It is the accumulation time of radioactive nuclides on the soil surface; It is the share of river water for irrigation. It is the effective rate constant for the removal of radionuclides from the soil surface. It is the rate constant for the elimination of other nuclides besides decay; It refers to the half-life of a nuclide; t h It refers to the time from harvesting to consumption of agricultural products.
8. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 1, characterized in that, The specific steps for calculating the corresponding public dose based on environmental parameters corresponding to a selected time, the location of key groups and the set of irradiation pathways, and the concentration of radionuclides in near-field water also employ the following formula: ； in, It refers to the proportion of radioactive nuclides in animal products; It refers to the time from the slaughter of animals to the consumption of their products by humans. This refers to the daily feed intake of animals; It is the specific activity of radionuclides in animal feed.
9. The dose assessment method for radioactive waste disposal sites based on biosphere evolution according to claim 1, characterized in that, The method further includes: At each time step on the time scale, environmental parameters, key group locations, and exposure pathway sets are automatically updated to calculate public doses.