PERFORMANCE ASSESSMENT SYSTEM FOR DEEP GEOLOGIC REPOSITORY FOR RADIOACTIVE WASTE DISPOSAL

Abstract

A performance assessment system for deep geologic repository for radioactive waste disposal is introduced to integrate a number of independent sub-system to perform the repository assessments in a systematic way under computer-based environment. Basically, the sub-system includes the input data file preparation sub-system for near-field/far-field multiple running, the near-field/far-field multiple running sub-system and the uncertainty and sensitivity analysis sub-system. With the system, the assessment for the deep geologic repository for radioactive waste disposal in many aspects can be achieved more completely and precisely.

Description

TECHNICAL FILED

The present invention generally relates to a performance assessment system of a deep geologic repository for the radioactive waste disposal, more particularly, to a system of assessing the long-term resistance function of the nuclide transportation before the actual construction and operation of the final deep geologic repository for the radioactive waste disposal.

TECHNICAL BACKGROUND

As FIG. 1 showed, the concept of the radioactive waste deep geologic repository equipping with the multiple barriers that has been considered as the most feasible and reliable final disposal method for the radioactive waste globally. Radioactive waste shall be long-term and permanently isolated from biosphere, thus it sets the multiple barrier system to dispose the radioactive waste. Basically speaking, the multiple barrier system is composed of the engineered barrier and natural barrier systems, which can be used to retard the release and transportation of the radioactive nuclides in order to ensure the safety and reliability of the final repository; therefore, before actual construction and operation of the final repository, the long-term retarding function of the nuclide transportation shall be precisely and completely assessed in advance.

In addition, the construction of the multiple barrier system will expend considerable resources (time and money), and the isolation effect between the radioactive waste and biosphere after completing the construction will acutely affect human living and life in the future; therefore, the assessment process has become extremely important before actual construction of the repository.

Presently, in the field of disposing radioactive waste, there is not yet a professional and complete assessment system which can precisely and completely to assess the isolation effect between the buried and disposed radioactive waste and the biosphere in order to be the basis of constructing the final repository for the radioactive waste disposal.

This said invention of the radioactive waste deep geologic repository performance assessment system provides a precise and complete assessment direction for assuring the safety and reliability of the radioactive waste final repository, which can precisely assess the long-tem retarding function of nuclide transportation and the isolation effect between the radioactive waste and the biosphere before actually constructing and operating the radioactive waste final repository; in addition, it is undoubtedly an optimal solution in the field of assessing the radioactive waste disposal.

TECHNICAL SUMMARY

Main purpose of this said invention is to provide a performance assessment system for supplying complete assessment information on the long-term retarding effect of the radioactive waste nuclide transportation before actually constructing and operating the radioactive waste final repository for the radioactive waste deep geologic repository.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is the conceptual diagram of the radioactive waste deep geologic repository which is equipped with the multiple barriers.

FIG. 2 is the main framework of this said invention for the radioactive waste deep geologic repository performance assessment system (hereinafter referred to as this said system).

FIG. 3A is the diagram of the analytic and applied scope for this said system's near-field release assessment sub-system, far-field release assessment sub-system and biosphere dose assessment sub-system.

FIG. 3B is the conceptual diagram of transportation for this said system's near-field release assessment sub-system near-field.

FIG. 3C is the conceptual diagram of transportation for this said system's far-field release assessment sub-system far-field.

FIG. 4A is the diagram of near-field release assessment data input sub-system for this said system's basic data input sub-system.

FIG. 4B is the diagram of the far-field release assessment data input sub-system for this said system's basic data input sub-system.

FIG. 4C is the diagram of the biosphere dose assessment data input sub-system for this said system's basic data input sub-system.

FIG. 5 is the diagram of the sampled result from the implemented nuclide parameter data sampling process for this said system's parameter sampling sub-system.

FIG. 6 is the diagram of the preparation function of the data input file for this said system's multiple running of the said near-field release assessment sub-system.

FIG. 7 is the diagram of the preparation function of the data input file for this said system's multiple running of the said near-field release assessment sub-system of the said mear-field release assessment sub-system after the execution of the preparation function is complete.

FIG. 8 is the diagram of implementing the multiple running function for this said system's near-field release assessment sub-system.

FIG. 9 is the diagram of parameter sensitivity for implementing the near-field release assessment sub-system of this said system.

FIG. 10A is the diagram parameter sensitivity for implementing the near-field release assessment sub-system of another said system.

FIG. 10B is the diagram parameter sensitivity for implementing the near-field release assessment sub-system of another said system.

FIG. 10C is the diagram parameter sensitivity for implementing the near-field release assessment sub-system of another said system.

FIG. 11 is the diagram of the sampled result from the implemented near-field release assessment sub-system for this said system by using the Monte Carlo Random Sampling process.

FIG. 12 is the diagram of the sampled result from the implemented near-field release assessment sub-system for this said system by using the Latin Hypercube Sampling process.

FIG. 13 is the diagram of the sampled result from the implemented near-field release assessment sub-system for another said system by using the Latin Hypercube Sampling process.

FIG. 14 is the diagram after implementing the multiple running of the said near-field release assessment sub-system (multiple running) for this said system.

FIG. 15 is the diagram after implementing the multiple running of the said near-field release assessment sub-system (multiple running) for another said system.

FIG. 16 is the diagram of the file control sub-system for implementing the far-field release assessment sub-system of this said system.

FIG. 17 is the diagram of the nuclide decay chain, half-life, sorption coefficient data for implementing the file control sub-system of far-field release assessment sub-system in this said system.

FIG. 18 is the diagram of the data of natural barrier system (NBS) property for implementing the far-field release assessment sub-system of this said system.

FIG. 19 is the diagram of establishing the new data file in the file control sub-system of implementing the far-field release assessment sub-system for another said system.

FIG. 20 is the diagram of the result after implementing the far-field release assessment sub-system for this said system.

FIG. 21A is the diagram of implementing the file control sub-system of the far-field release assessment sub-system for another said system.

FIG. 21B is the diagram of implementing the file control sub-system of the far-field release assessment sub-system for another said system.

FIG. 22A is the diagram of implementing the Review function sub-system of the far-field release assessment sub-system for this said system.

FIG. 22B is the diagram of implementing the Review function sub-system of the far-field release assessment sub-system for another said system.

FIG. 22C is the diagram of implementing the Review function sub-system of the far-field release assessment sub-system for another said system.

FIG. 22D is the diagram of implementing the Review function sub-system of the far-field release assessment sub-system for another said system.

FIG. 22E is the diagram of implementing the Review function sub-system of the far-field release assessment sub-system for another said system.

FIG. 23 is the diagram of implementing the Drawing function sub-system of the far-field release assessment sub-system for this said system.

FIG. 24 is the diagram of implementing the preparation function of data input file for the multiple running of the said far-field release assessment sub-system in this said system.

FIG. 25 is the diagram of implementing the data input of the far-field release assessment sub-system for this said system.

FIG. 26 is the diagram of implementing multiple running of the said far-field release assessment sub-system for this said system.

FIG. 27 is the diagram of implementing the correlative function for the variable sensitivity of the far-field release assessment sub-system in this said system.

FIG. 28 is the diagram of implementing the multiple running of random sampling for the far-field release assessment sub-system in this said system.

FIG. 29 is the diagram of implementing the multiple running of Latin hypercube sampling for the far-field release assessment sub-system in this said system.

FIG. 30 is the diagram of implementing the multiple running of Latin hypercube sampling for the far-field release assessment sub-system in another said system.

FIG. 31 is the diagram of implementing the multiple running of the said far-field release assessment sub-system for another said system.

FIG. 32 is the diagram of implementing the multiple running for this said system.

FIG. 33 is the diagram of implementing the file mergence function for the multiple running function sub-system in this said system.

FIG. 34 is the diagram of the result after implemented the multiple running function sub-system of this said system.

FIG. 35 is the diagram of the result after implemented the multiple running function sub-system of another said system.

FIG. 36 is the diagram of implementing the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 37 is the diagram of implementing the multiple running of the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 38 is the diagram of saved data that after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 39 is the diagram of the saved fixed time release rate CCDF data after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 40 is the diagram of saved percentage total release rate curve after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 41 is the diagram of probability analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 42 is the diagram of the fixed time release rate in the probability analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 43 is the diagram of the release rate peak in the probability analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 44 is the diagram of the peak occurrence time in the probability analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 45 is the diagram of the sensitivity analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 46 is the diagram of the fixed time release rate in the sensitivity analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 47 is the diagram of the time analysis in the fixed time release rate of the sensitivity analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 48 is the diagram of the time analysis in the fixed time release rate of the sensitivity analysis after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 49 is the diagram of the result that obtained from the sensitivity analysis and Rank transformation after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 50 is the diagram of the result that obtained from the sensitivity analysis and data Log transformation after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 51 is the diagram of the drawing scatter plot function that obtained from the sensitivity analysis and Drawing function after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 52 is the diagram of the figure magnification function that obtained from the Drawing function after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 53 is the diagram of the scatter plot parameter name tag display function item that obtained from the Drawing function after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 54 is the diagram of the magnified CCDF function that obtained from the Drawing function after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 55 is the diagram of adding the basic assessed result after implemented the uncertainty and parameter sensitivity analysis function for this said system.

FIG. 56 is the diagram of the program verification function item after implemented the uncertainty and parameter sensitivity analysis function for this said system.

DETAILED DESCRIPTION OF EXEMPLIFICATION

This said system, in order to assure the safety and reliability of the radioactive waste final repository, will thus propose a precise and complete assessment solution for the long-term retarding function of the nuclide transportation before actually constructing and operating the final repository.

The function assessment of the radioactive waste final repository is a difficult process which has involved numerous influential factors, and the geologic heterogeneity of rock and the geographic environment change in a long period of time; in addition, as for many factors (or can be called as the parameter or variable) of affecting the isolation function of the repository, their values certainly are unable to be ascertained; thus, when assessing the function of a repository, it is usually common to set each factor's value as certain reasonable distribution pattern and scope, and by means of the parameter sampling to conduct multiple computer calculations, and to use the uncertainty and parameter sensitivity analysis, such as Monte Carlo assessment process and technique, to dispose then.

According to the path of transporting the radioactive nuclide, the total system function assessment of the repository can be divided into various sub-systems, such as the near-field transportation, far-field transportation (or entitled as the geological migration) and the biosphere transportation to carry out the assessment. The near-field transportation assessment includes that how the nuclide can be transported by passing through those barriers: the waste form, waste canister, buffered material layer, backfilled layer, and the excavation disturbed zone (EDZ), etc., (collectively entitled as the engineered barrier system, EBS); as for the far-field transportation assessment, it shall be assessed for that how the nuclide can be penetrated through the host rock to transport to human's living environments; and as for the biosphere transportation assessment, it shall be assessed for that how the nuclide can reach to reach human body by means of drinking well water and food-chain of the biosphere.

Please refer to FIG. 2, which is a main framework figure of this said system. From FIG. 2, we know that this said system is included 9 major sub-systems: the basic data input sub-system21, fixed parameter setting sub-system22, distributive parameter sampling sub-system23, Latin hypercube or random sampling sub-system24, near-field release assessment sub-system25, far-field release assessment sub-system26, biosphere dose assessment sub-system27, uncertainty analysis sub-system28 and sensitivity analysis sub-system29.

Among which, the scope of analytic application for the near-field release assessment sub-system, far-field release assessment sub-system and the biosphere dose assessment sub-system will be shown as in FIG. 3A; the concept of near-field transportation will be shown as in FIG. 3, and the idea of the far-field transportation will be shown in FIG. 3 C.

In addition, the uncertainty and sensitivity analysis sub-systems include two sets of sampling techniques: Latin Hypercube Sampling and Random Sampling, for further application, and adopt the Stepwise Regression Analysis technique to conduct the analysis of the parameter sensitivity.

Herein described the functions of these 9 major sub-systems that contained in this said system as follows; first of all, please refer to FIG. 4A, FIG. 4B and FIG. 4C. FIG. 4A is the figure of the near-field release assessment data input sub-system in this said system's basic data input sub-system, and from it we can know that the input data of the near-field release assessment data input sub-system included:

(1). Waste property data 42 includes the data of the starting time of inventory and the time of complete dissolution of the waste form;

(2). Waste canister property data 43 includes the data of life-span, corrosion product density, inner radius, outer radius, length, corrosion product porosity, and diffusion coefficient in the corrosion product;

(3). Buffered layer property data 44 includes the data of density, porosity, outer radius, and diffusion coefficient;

(4). Excavation disturbed zone (EDZ) property data 45 includes the data of rock density, outer radius, porosity and diffusion coefficient;

(5). Host rock property data 46 includes the data of Darcy flow rate, fracture diffusion coefficient, fracture spacing, fracture opening;

(6). Title of released nuclide, half-life and decay phase data 47;

(7). Host nuclide, half-life and sub-nuclide data 48;

(8). Nuclide title, inventory and Instant Release Fraction (IRF) data 49; and

(9). Chemical element solubility and sorption coefficient data 410.

Please also refer to FIG. 4B, which is the figure of the far-field release assessment data input sub-system in this said system's basic data input sub-system, from FIG. 4B we can know that the input data by the far-field release assessment data input sub-system included:

• • (1). Geometric property data 401 includes the data of the geological transportation distance, waste pit spacing, nuclide sorption coefficient, fracture surface sorption depth, waste canister length, fracture spacing, fracture opening, fracture transportation division number, rock mass diffusion division number;
  • (2). Host rock property data 402 includes the data of the density and porosity data;
  • (3). Transportation property data 403 includes the data of Darcy flow rate, fracture diffusion coefficient, rock mass diffusion coefficient and the dispersivity;
  • (4). Nuclide flux input data 404;
  • (5). Nuclide concentration output time 405;
  • (6). Release rate assessment of the nuclide decay chain data 406, which includes the data of host nuclide, half-life and sub-nuclide; and
  • (7). Chemical element's sorption coefficient data in the host rock 407.

Next, please refer to FIG. 4C, which is the figure of biosphere dose assessment data input sub-system in this said system's basic data input sub-system. Form FIG. 4C we know that the input data of this biosphere dose assessment data input sub-system included: far-field nuclide release rate, the well entering percentage of nuclide 431, the annual water output volume of the well 432, annual water consumption for individual 433 and the annual dose rate.

After completed the input of basic data, next the function and operation of this said invention's near-field release assessment sub-system will be explained. Before implementing this said system's near-field release assessment sub-system, the distributive parameter sampling sub-system of this system shall be implemented in advance, which included 2 sampling methods: Latin Hypercube Sampling and Monte Carlo Random Sampling.

Please refer to FIG. 5, which is the figure of their result after implemented the nuclide parameter data sampling by this said system's parameter sampling sub-system.

After completed the sampling process, the near-field release assessment of the radioactive waste near-field release can then be implemented; and before implementing multiple running of the near-field release assessment sub-system, the preparation function of the data input file shall be implemented in advance as showed in FIG. 6. FIG. 6 showed that the preparation function of the data input file for multiple running of this said invention's near-field release assessment sub-system, and its main data related input area is: (1) the disposal facility design and the geologic property correlative setting zone 61 which contained 3 data setting zones—the well-obtained data of the parameter list 611; the correlative setting list 612 which contained the well-obtained parameter data and used the near-field release assessment sub-system to assess the variable sensitivity; and the random parameter list 613 which is contained in the near-field release assessment sub-system; (2) Chemical element solubility correlative setting zone 62; (3) Chemical element's correlative setting zone 63 of sorption coefficient for buffered materials; (4) Chemical element's correlative setting zone 64 of sorption coefficient for host rock; (5) Chemical element's correlative setting zone 65 sorption coefficient for the erosion object in waste tanks, etc.

As implied in the title, the correlative setting zone 61 of disposal facility design and geologic property is used to connect with uncertain parameters that are related the disposal facility design and the geologic property; the correlative setting zone 62 of chemical element solubility is used to connect with uncertain solubility of chemical element; the correlative setting zone 63 is used to connect with uncertain sorption coefficient of buffered materials for chemical element; the correlative setting zone 64 is used to connect with the uncertain sorption coefficient of host rock for chemical element; and the correlative setting zone 65 is used to connect with the uncertain sorption coefficient of the erosion object in waste tanks for chemical element. Each property correlative zone contained 3 data display zones—that is, the left-side well-sampled parameter title listing zone, the right-side random parameter listing zone or element title listing zone of the near-field release assessment sub-system; and correlative listing zone in the middle.

After implemented the distributive parameter sampling sub-system, this said invention system can then be implemented the “Single Running” or “Multiple Running” function for near-field release assessment sub-system; at this moment, users shall input proper data in advance to provide the near-field release assessment sub-system for carrying out the assessment calculation. The screen of the result after inputted the data input file is showed as in FIG. 4A, at this time, if users selected the “Single Running” function 41, then the system will directly implement single calculation of the near-field release assessment sub-system in accordance with the data that displayed on the screen.

After inputted a complete data input file of the near-field release assessment sub-system, the “multiple running” function of this said invention can be implemented then. The implementation figure of the “multiple running” function for the near-field release assessment sub-system is showed as in FIG. 8, users can select the parameter, which they wanted to explore its variable sensitivity, from the list that titled as the “parameter list of the well-obtained data” 81 in the system, and such parameter will be automatically added into the list of the “correlative setting list of the well-obtained parameter and the exploration of variable sensitivity by using the near-field release assessment sub-system” 82 as showed in FIG. 9.

After selected the required parameter from the “correlative setting list of the well-obtained parameter and the exploration of variable sensitivity by using the near-field release assessment sub-system”, the system will display a data input zone for “the random parameter list of near-field release assessment sub-system” 101 as showed in FIG. 10A. In the list of “the random parameter list of near-field release assessment sub-system” 101, users can select the correlative parameter from the near-field release assessment sub-system which they wanted to assess its variable sensitivity, and it will yield the result of zone 102 as showed in FIG. 10A. Now, it needs to be carefully concerned that such sensitivity exploring parameter shall not be correlated with the parameter of the near-field release assessment sub-system; however, the correlating parameter has to be in the selected status, and the correlation can then be established.

Related parameter that selected from the correlative near-field release assessment sub-system will be showed the selected status in the parameter value column of the near-field release assessment sub-system as showed in 103 and 104 of FIG. 10B, and which indicated the parameter of assessing variable sensitivity has been selected already.

Similarly speaking, users can select the parameter of assessing variable sensitivity from those data zones: “chemical element solubility”, “chemical element's sorption coefficient in host rock”, “chemical element's sorption coefficient in buffered materials” and “chemical element's sorption coefficient in corrosion object in waste tanks”.

After completed the correlation setting of the assessing sensitivity parameter, the “multiple running” function 105 of FIG. 10A shall be implemented then; if only used the Monte Carlo Random Sampling to conduct the sampling process, the result will be showed as the same as in FIG. 11; if only adopted the Latin Hypercube Sampling to implement the sampling process, then only the item of “random arrangement” 1101 function for the “data arrangement” in the parameter sampling system can be implemented then, and the result will be as same as the result of using the Monte Carlo Random Sampling; in addition, the result will be showed as the same as in FIG. 11; if increasingly implemented the item of “non-correlative arrangement” 1201 function for the “data arrangement”, and then the result will be showed as the same as in FIG. 12; in addition, if increasingly implemented the item of “specific correlative arrangement” 1301 function for the “data arrangement”, and then the result will be showed as the same as in FIG. 13.

When implementing the type of selecting data arrangement from the “selecting data arrangement type, this said system will display the result as showed in FIG. 14, and this said system will confirm that whether users will save/store related data of multiple running into the certain category or not. At this moment, the implementation of calculation for this said system's near-field release assessment sub-system can then be considered as a completion, select the block of “assessed result” from the “near-field release assessment sub-system”, this said system will display the near-field release assessed result as showed in FIG. 15 then.

As showed in FIG. 4A, this said system's near-field release assessment sub-system includes functions of “File” 411, “Save As” 412, “Insert” 413, “Clear” 414, “Review” 415, “Drawing” 416 and “Work Directory” 417; in addition, those functions and follow-up explanations are almost identical to those functions in this said system's far-field release assessment sub-system, such as “File”, “Save As”, “Insert”, “Clear”, “Review”, “Drawing” and “Work Directory”. Except the “File” function in this said system's near-field release assessment sub-system, it is not included those functions of file mergence , file name change and file delete , other functions will be identical to each other, thus it will not explain herein, and it will be explained in this said system's far-field release assessment sub-system then.

Next, it will explain the function and its operation method for this said system's far-field release assessment sub-system, such as this said system's “far-field release assessment sub-system” function as showed in FIG. 4B. After entered the said far-field release assessment sub-system, from FIG. 4B, we know the said sub-system is included those file control and Drawing control functions, such as “File” 421, “Save As” 422, “Insert” 423, “Clear” 424, “Review” 425, “Drawing”, and “Work Directory” 427.

As showed in FIG. 4B, after completed the basic data input for this said invention's far-field release assessment sub-system, its data will be included “Geometry property” 401, “Host Rock property” 402, “Transport property” 40, “program setting for the said far-field release assessment sub-system” 408, “Nuclide Flux Input File” 404, “Nuclide Concentration Output Time” 405, “the nuclide decay chain that needs to be conducted the release rate assessment” 406 and “element's sorption coefficient in the host rock” 407 nuclide property data setting and the assessment implementing functions.

“Geometry property” 401, “Host Rock property” 402 and “Transport property” 403 are jointly titled as the data of natural barrier system (NBS) property. These 2 parts, “the nuclide decay chain that needs to be conducted the release rate assessment” 406 and “element's sorption coefficient in the host rock” 407, can be jointly named as the data of the nuclide decay chain, half-life and sorption coefficient. As for these 3 parts, the “program setting for the said far-field release assessment sub-system”408, “Nuclide Flux Input File” 404 and “Nuclide Concentration Output Time”405, can be jointly named as the system setting data of calculation implementation.

With particular attention, in the “File name” column of “Nuclide Flux Input File” 404, you have to key in the correct operation is performed near-field release assessment sub-system after the output of the nuclide flux output data, as the implementation of far-field release assessment sub-system required for the nuclide data flux input file

The function and its operation method of this said far-field release assessment sub-system will be described as follows, basically speaking, the nuclide transportation data used by this said far-field release assessment sub-system that will be the result data after assessed this said near-field release assessment sub-system, that is, the input data used by this said far-field release assessment sub-system is the assessed output data for this said near-field release assessment sub-system; thus, the “File” function will be explained firstly for this said far-field release assessment sub-system to understand how to make the assessed result of the near-field release assessment sub-system to be the assessment data for this said far-field release assessment sub-system. First of all, as showed in FIG. 16, there are 4 sub-function items will be displayed in the scroll menu of the “File” function items, such as Open Old File , Establish New Data File , Implement Previous Data File and File Processing , and the operation and function of these 4 sub-function items for this said function item that will be explained.

As showed in FIG. 16, 4 sub-functions in the Open Old File , such as the files of nuclide decay chain, half-life, sorption coefficient data 161; data of natural barrier system (NBS) property 162, setting data of calculation implementation 163 and above-mentioned 3 data (one complete implementation case) 164. If the data existed, then the established old file can be selected from Open Old File function item. Thus, operation and function of these 4 sub-function items of this said function item that will be introduced respectively.

When selecting the data of the nuclide decay chain, half-life and sorption coefficient 161 function item as showed in FIG. 17, users can be selected the Open Old File and key the nuclide-related data in the “data of the nuclide decay chain, half-life and sorption coefficient”161.

It needs to put particular attention that users have to select the output data of the “title of released nuclide, half-life and decay chain” that yielded from implementing the calculation of the near-field release assessment sub-system which will then be able to consistent with the nuclide transportation types of the near-filed assessment.

If the data input is good, then the data will display as showed in FIG. 4B, now, if users implemented “Single Running” function 4201, this said system will directly implement single calculation of this said far-field release assessment sub-system according to the data on the screen.

After selected certain nuclide data from the “the nuclide decay chain that needs to be conducted the release rate assessment” function, when selected the “Delete” function, such nuclide data will be deleted; when certain element's isotope nuclide has been completely deleted, then such element's sorption coefficient in the zone of “element's sorption coefficient in the host rock” will be automatically deleted as well. Users can make modification and revision of parameter input in the parameter and noted input zone.

When implementing the data of natural barrier system (NBS) property 162 function, as showed in FIG. 18, this said system will display a file list block, and the listed File name is the old File name for the data of natural barrier system (NBS) property that established by users before. After selected such file, the previous data file can be opened to input the data of natural barrier system (NBS) property.

When selecting the program implementing setting data 163 function, such system will display a file list, and the listed file is the File name of this said system program implementing setting data that established by users previously. After selected the file, the previous data file can be opened then.

After selected the abovementioned 3 data (one complete implementation case) 164 function item, such system will display a file list, the listed file is the File name that established by users previously. After selected the file, the previous data file can be opened then. When implementing this function item, it can be concurrently read the aforesaid set file data of files of nuclide decay chain, half-life, sorption coefficient data 161, data of natural barrier system (NBS) property 162, and program implementing setting data 1633.

The main purpose of these abovementioned 4 sub-function item in the Open Old File function is to increase the freedom for users to select different data files freely to compose of a new implementing parameter content for this said far-field release assessment sub-system, or directly click on the abovementioned 3 data (one complete implementation case) function item to read a complete data to implement this said far-field release assessment sub-system.

If users are the first-time users for this said system, then there's no previous data file available; however, now the function item of Establish New Data File can be made use of establishing the New Data File, as showed in FIG. 19. After selected the Establish New Data File function, users can establish new data in such function to conduct the assessment.

It needs to put extra attention, the analyzing nuclide shall be identical to the nuclide that analyzed by implementing the near-field release assessment sub-system; thus, when implementing the near-field release assessment sub-system, this said system will be automatically yielded the output data file of the “title of released nuclide, half-life and decay chain” and to be used for implementing the far-field release assessment sub-system; as a result, it can be identical to the transportation nuclide type for this said near-filed assessment, in the FIG. 19, “the nuclide decay chain that needs to be conducted the release rate assessment” 191 will display the nuclide data that needs to conduct the assessment.

After all nuclide data have completely inputted, users can select their required element from the list of “element's sorption coefficient in the host rock” 192, and then they can modify and revise each element's sorption coefficient.

After completely inputted the related data, users can click on “Save As” function to save the inputted related data.

Within these 3 data zones, such as the “data of natural barrier system (NBS) property”, “sub-system implementing setting data” and “data of the nuclide decay chain, half-life and sorption coefficient”, data has to be available and integral in these zones; otherwise, the Single Running function or multiple running function of this said far-field release assessment sub-system cannot be implemented.

After selected the function item of Implement Previous Data File , such system will display a file list, listed File name is the File name that established by users previously; in addition, after users selected their required data file, then implemented the FIG. 19's “Single Running” 193 to implement the calculation and assessment. After completely implemented this said far-field release assessment sub-system, a nuclide release flux file and an implementation file will be generated then.

After completely implemented the far-field release assessment sub-system to conduct the calculation and assessment, its diagram is as showed in FIG. 20. In FIG. 20, these calculated results of this said far-field release assessment sub-system can be drew as a time-changed figure of the nuclide release flux 201, and the data of calculated result that displayed by words 202.

file function in the far-field release assessment sub-system, its File Processing function item can be divided into 3 sub-function items, such as the file mergence , filename changing and file delete , and their function property and operation methods will be introduced as follows.

After selected the file mergence , as showed in FIG. 21A, the upper zone is the merging file list 211, when selecting the merging file, then the File name will be duplicated to the bottom menu 212, after implemented is as showed in FIG. 21B. The emerging file data that has already merged and displayed by words in 214, as well as drew the figure of merged data 213.

When implementing the filename changing or file delete , users can change the file name or delete the name change or delete the file, and it will not explain herein.

After users opened Old File or newly added the data file, the can use the “Save As” function to save file. The “Save As” function has 4 sub-function items, such as the data of the nuclide decay chain, half-life and sorption coefficient , data of natural barrier system (NBS) property , program implementing setting data , and abovementioned 3 data (one complete implementation case) . After additionally increased and modified the data, users can select different sub-function items to save different file data, and it will not explain herein.

After opened Old File or newly added the data file, the “Insert” function can then be applied. After selected the “Insert” function item, users can use such function item to insert other nuclide items to connect and form a new nuclide data content.

After opened Old File, the “Clear” function item can then be applied; in addition, after used such function, users can clear nuclide decay chain, half-life, sorption coefficient, decay chain, and the element's sorption coefficient in the host rock.

After implemented “Review” function, as showed in FIG. 22A, the “Review” function item will then contained following 3 sub-function items, such as the program implementing data input file 2201, program implementing output file 2202 and program implementing output explanatory file 2203, and they will be introduced as follows.

After selected the program implementing data input file 2201 function item, the program implementing data input file function item is also included 2 sub-function items, such as the latest saved file 22011 and previously established file 22012, as showed in FIG. 22B. If users have not yet implemented the function of Implement Previous Data File or Save As above mentioned 3 data (one complete implementation case) , and then the sub-function item of latest saved file for the program implementing data input file is unable to work then.

When selected the function item of previously established file and the proper file data, as well as implemented, the implemented result is as showed in FIG. 22C. File inputted by such function is the input file that established for the data of this said far-field release assessment sub-system.

After selected the previously established file of the program implementing output file function 2202, and after selected the proper file data and implemented, the implemented result is showed in FIG. 22D, and the input data of this said function is the nuclide release flux output data when implemented the said far-field release assessment sub-system.

After selected the previously established file of the program implementing output explanatory file function 2203, and after selected the proper file data and implemented, the implemented result is showed in FIG. 22E, and the input data of this said function is the output explanatory file when implemented the said far-field release assessment sub-system.

In FIG. 4B, after selected the “Drawing” 426 function item, “Drawing” function item is contained 5 sub-function items, such as the recently implemented case , previously implemented case , modified Y-axis , modified X-axis and adding figure .

If users have not yet calculated and implemented the far-field release assessment sub-system, then only the sub-function item of previously implemented case is effective, and other 4 sub-function items will be temporarily ineffective. When users after selected the sub-function items of recently implemented case or previously implemented case , after properly selected the file, and then those sub-function items of modified Y-axis , modified X-axis and adding figure can then be effective. Operation and function of these 5 sub-function items under this said function item will be explained as follows.

After selected the recently implemented case function item, this said system will display the output result on the screen of data and figure that recently implemented the far-field release assessment sub-system, as showed in FIG. 20; at this moment, users can conduct other sub-function items of the “Drawing” function item, such as modified Y-axis , modified X-axis and adding figure , to clearly observe the changing situation of release flux for each nuclide.

After selected the previously implemented case function item, and selected the file, this said system will display the screen of the output result and figure as showed in FIG. 20.

After selected the modified Y-axis function item, the modified Y-axis function item is included 2 sub-function items, such as the maximum value and minimum value . If selected the maximum value function item, as showed in FIG. 23, users can input Y-axis's maximum value 2301, and then change Y-axis's maximum value in the figure. Similarly, if adopted the minimum value function item, and the minimum value can be changed in the figure.

As for the effect of modified X-axis function item, as the introduction of the aforesaid modified Y-axis function item, X-axis's maximum/minimum value can be changed in the figure, and it will not explain herein.

After selected the adding figure function item, this said system will display a file list, and after selected the designated file, the figure of such file can be stacked onto the original figure; as a result, users can then be displayed different output results on a same screen, and it will also not explain herein.

After selected the “Work Directory” function item, “Work Directory” contained 3 sub-function items, such as Display the Current Work Directory , Change Work Directory and Establish New Work Directory , as implied by the names, their functions are respectively notifying users of the current directory path of system and data, and users can then be selected the designated path to understand the directory path for their system and data, and they can also be established the path of a new directory to add a new directory as well.

Next, it will be explained the preparation function of the input file for multiple running data in the said far-field release assessment sub-system, the implementation of such preparation function of the input file for multiple running data in the said far-field release assessment sub-system is similar to the implementation of preparation function of the input file for multiple running data in the said near-field release assessment sub-system, except to implement the parameter sampling system in advance, the type of sampling parameter arrangement that adopted by this said far-field release assessment sub-system shall be identical to the type of sampling parameter arrangement that used by the said near-field release assessment sub-system. Since the output data file of nuclide release rate that generated from the said near-field release assessment sub-system shall be used when implementing the said far-field release assessment sub-system; in addition, as considering the consistence of the parameter for further analysis process, users have been recommended to continuously implement the preparation function of the input file for multiple running data in the said far-field release assessment sub-system after completely implemented the preparation function of the input file for multiple running data in the said near-field release assessment sub-system to facilitate the further analysis of the near-field and far-field nuclide release uncertainty and parameter sensitivity.

The preparation function of the input file for the multiple running of the said far-field release assessment sub-system data, as showed in FIG. 24, which can be divided into 2 major nuclide property data correlative zones: (1) disposal facility design and geologic property correlative setting zone 241 (the said parameter list of the well-obtained data 2411, the parameter of well-obtained data and the correlative setting list of the said far-field release assessment sub-system that needs to be explored the variable sensitivity 2412, the random parameter list in the said far-field release assessment sub-system 2413). (2) the correlative setting zone of chemical element's sorption coefficient in the host rock 242.

As implied in the name, the property correlative zone 241 is used to connect with uncertain parameters that related to the disposal facility design and geologic property, etc.; the property correlative zone 242 is used to connect with the uncertain sorption coefficient for related chemical element in host rock. Each property correlative zone is contained 3 blocks, such as the well-sampled parameter name list block 2411 on the left; the random parameter list block 2413 in the said far-field release assessment sub-system on the right; or the element name list block 2423; and the correlative list block in the middle. Functions of these 2 property correlative zones will be explained as follows.

After implemented the parameter sampling system and completed the multiple running of the said near-field release assessment sub-system, this said system will be automatically accessed to the said far-field release assessment sub-system, users can then select the proper file data to input a complete data input file of the said far-field release assessment sub-system to the far-field release assessment sub-system, as showed in FIG. 25. If, at this moment, users selected the function of “Single Running” 251, and then a single assessing function can then be directly implemented for the said far-field release assessment sub-system.

Next, implemented the function of “multiple running” 261 in FIG. 26, in the list 262 of “parameter list of the well-obtained data”, selected the parameter of assessing variable sensitivity, and then such parameter will be automatically added into the list 263 of “parameter table of the well-obtained data and using the far-field release assessment sub-system to explore the correlative setting for variable sensitivity”, users can also be deleted the parameter from the list of the “parameter table of the well-obtained data and using the far-field release assessment sub-system to explore the correlative setting for variable sensitivity”.

As showed in FIG. 27, in the list 271 of “parameter table of the well-obtained data and using the far-field release assessment sub-system to explore the correlative setting for variable sensitivity”, after selected the a parameter, this said system will display a “the random parameter list in the said far-field release assessment sub-system” 272, and in such list, users can then be selected the parameter of the said far-field release assessment sub-system that related to the parameter needs to be assessed its variable sensitivity.

Similarly, selected the parameter from the list of “chemical element's sorption coefficient in the host rock” 273 that needs to be assessed the variable sensitivity, then the parameter will be added into the attached list, and users can also delete the parameter that they want to delete in the list, and it will also not explain herein and it will also not explain herein.

After completed the correlation setting of the sensitivity parameter that needs to be assessed, then selected the “multiple running” function 274 from the FIG. 27; in addition, adopted the Monte Carlo Random Sampling method to sample in the multiple running of the said near-field release assessment sub-system. Since Monte Carlo Random Sampling is only adopted the random sampling method, thus the current multiple running of the said far-field release assessment sub-system will only display the option of “using random arrangement data”, as showed in FIG. 28; in the multiple running of the said near-field release assessment sub-system the Latin Hypercube Sampling is adopted to sample and implemented the option of “using random arrangement data”, and then the result for the multiple running of the said far-field release assessment sub-system will be consistent to the sampling operation by using Monte Carlo Random Sampling, as showed in FIG. 28; if adopted Latin Hypercube Sampling to sample the multiple running of the said near-field release assessment sub-system, and after implemented the option of “using non-correlative arrangement data”, then this multiple running of the said far-field release assessment sub-system will display the result as showed in FIG. 29; if applied Latin Hypercube Sampling to sample the multiple running of the said near-field release assessment sub-system, and after implemented the option of “using specific correlative arrangement data”, then this multiple running of the said far-field release assessment sub-system will display the result as showed in FIG. 30. Since the consistence for the calculation and analysis under the multiple running of the said far-field release assessment sub-system system, the previously selected parameter sampling method for the multiple running of the said near-field release assessment sub-system that can only be displayed without any change or modification.

In FIG. 30, after selected the “selecting data arrangement type”, the calculated result is as showed in FIG. 31. This said system will make sure of whether users will save the related data of this multiple running into the current directory or not, or users can change the sub-directory or can self establish a new sub-directory. In addition, this said system will arrange the parameter that needs to be analyzed the sensitivity according to the selected type of parameter arrangement to orderly write the parameter into the file they named in order to maintain the consistence in the type of parameter sampling arrangement for the near-field release assessment sub-system and multiple running of the said far-field release assessment sub-system system.

Next, this said system will request users to input the nuclide flux output file name after completed the implementation of the said far-field release assessment sub-system that required for implementing the near-field release assessment sub-system.

This said system according to the selected data arrangement type to automatically complete the number of data input file that required for the said far-field release assessment sub-system, then this said system will automatically switch to the near-field release assessment sub-system and the multiple running of the said far-field release assessment sub-system; thus, the pre-operation process is completed for the multiple running of far-field release assessment sub-system.

Next, this said system's multiple running function will be introduced as follows, this said invention's multiple running function is designed by focusing on the multiple running near-field release assessment sub-system and/or far-field release assessment sub-system; therefore, before implemented this said multiple running function, the preparation function of the data input file shall be implemented in advance for the near-field release assessment sub-system and/or multiple running of the said far-field release assessment sub-system.

This said system's multiple running function is as showed in FIG. 32, and from FIG. 6.1.1, this said system multiple running function is included those sub-functions, such as “File”321, “Drawing”322 and “Work Directory” 323. From FIG. 32, this said system's multiple running function is also included those functions, such as “multiple running of the said near-field release assessment sub-system”324 and “multiple running of the said far-field release assessment sub-system”325. File names displayed in the function of “single click/double click on these following files” 326 that can be implemented the multiple running process, and these files have been established after respectively implemented the “preparation system of the input file for multiple running data of the said near-field release assessment sub-system data” and “preparation system of the input file for multiple running data of the said far-field release assessment sub-system data”. After selected certain file from these files, the bottom values of “final implementation round” 328 and “nuclide number”329 will be automatically set then. After accessed to this said system's multiple running function, the function of “Work Directory” can be selected to change the Work Directory to those sub-directory items which have already saved the multiple running files for the near-field release assessment sub-system and the said far-field release assessment sub-system.

The “File” function of this said system's multiple running function item is included 3 sub-function items, such as file mergence , filename changing and file delete , after selected the file mergence , as showed in FIG. 33, users can select those files they want to merge, and then the File name will be duplicated to the bottom menu. After selected the merging files, implemented the mergence function, this said system will conduct the file mergence function. This said system has file mergence function since the number of nuclide is too many, and it will be consumed a lot of time when implemented this said far-field release assessment sub-system; therefore, firstly to divided nuclide into several files to be individually implemented (only the decay chain related can be divided) to save the calculation time, wait to completely implement all divided files, and then merged the result of each implemented file into a complete output file.

Users can use the functions of filename changing and file delete to change the file names and delete the files that they want to change and delete.

The operation of “Drawing” function of this said system's multiple running function is totally identical to the “Drawing” 426 function in FIG. 4B, and it will not explain herein.

After selected the “Work Directory” function item, the “Work Directory” function item is included 3 sub-function items, such as Display the Current Work Directory , Change Work Directory and Establish New Work Directory , and the introduction is as follows.

This said system's multiple running functions contained those functions, such as Display the Current Work Directory , Change Work Directory and Establish New Work Directory , and make users to understand the current directory path for their system and data, and they can also be selected the designated directory path and established the path of a new directory and save the data into such new directory as well, and it will not explain herein.

“Suspension” function item can only be used for carrying out the program of multiple running near-field release assessment sub-system or the far-field release assessment sub-system, its function is to terminate the currently operating multiple running procedures. Click on the “suspension” function item that can suspend the implementing multiple running system, the current time consumption for this said system will be no longer to increase; however, the current starting of this said system for the near-field release assessment sub-system or the far-field release assessment sub-system that will not suspend for implementation, and it needs to be manually shut down or automatically suspended after completely the implementation.

If only needed to implement the multiple running of the said near-field release assessment sub-system function, wait for completing the parameter sampling and establishing the multiple running of those data input files, and said near-field release assessment sub-system system. After selected the “proceeding multiple running” function item (as showed in FIG. 326.1.1) of the “multiple running” function, it will access into the multiple running system; then, after selected the files from the list of “multiple running of the said near-field release assessment sub-system”, this said system will state the pre-set sampling number and nuclide number. Users shall be notified, when conducting the parameter sampling, it's better to complete the multiple running of the said far-field release assessment sub-system function, and complete the consistence for parameter arrangement type in order to facilitate to the use for further assessment and analysis of the said far-field nuclide release. If the data input file has been established previously, then the multiple running of the said near-field release assessment sub-system function, and it can be directly selected the “proceeding multiple running” function item of the “multiple running” system function.

FIG. 32, as for the “initiate implementation round” 327 column value in the “multiple running of the said near-field release assessment sub-system”324 list, its default value is 1, “final implementation round” 328 column value, and the default value is the parameter sampling number, both values can be changed and modified, and it means that users would like to start calculating by selecting the number of sampling data, the scope will be 1˜ parameter sampling number. Also, users can directly change the number, but the “initiate implementation round” 327 column value is unable to be less than 1, if its is less than 1, this said system is considered as 1, and it is unable to exceed “final implementation round” 328 column value, if it is exceeded the “final implementation round” 328 column value, the system will be considered as the “final implementation round” 328 column value. The “final implementation round” 328 column value cannot be less than the column value of “initiate implementation round”327; if so, the “initiate implementation round”327 column value, such system will be considered as the column value of the “initiate implementation round”327, and it cannot be exceeded the parameter sampling number; if so, this said system will be considered as the parameter sampling number.

In FIG. 32, the column of “nuclide number” 329 in the “multiple running of the said near-field release assessment sub-system” list is indicated that the set analyzing nuclide number in the data input file is unable to be modified or changed. After selected the “delete the data input file of the said near-field release assessment sub-system after implemented” 330 function, all input file after completely implemented the multiple running of the said near-field release assessment sub-system data that will be deleted then.

When each parameter has been completely inputted, and after completely implemented the first multiple running for the near-field release assessment sub-system, then the system will draw the first-round nuclide release quantity figure and will be automatically initiate the 2^nd-round analysis of implementation, as showed in FIG. 34, till completed all currently selected sampling number as showed in FIG. 35.

FIG. 34 showed the Drawing condition that the result after implemented the multiple running. From FIG. 34, each complete implementation for the time relationship between the total nuclide release flux and time that can be drew and displayed as in FIG. 341. If the number of rounds for previously implementing the multiple running, this said system will not be implemented again, and will directly be drew the Figure for the result which can save a great deal of time of implementation, such figure has also equipped a word block 342, which can input notes and words, and such FIG. 341 can be also shrunk and magnified to facilitate users' viewing.

From FIG. 34, the system function of “multiple running” is contained 4 display items, such as: “set time consumption” 343, “implementation round” 344, “current time consumption” 345, and “time increment” 346, and they will be introduced as follows: “set time consumption” 343 means the multiple running system that after the set time (unit is second), it will start to check and inspect all calling programs have completely implemented or not, and the default value is 1(second), which can also be changed. “time increment” 346 means that since the time setting, what is the time interval for this said system will check and inspect whether the calling program has been completely implemented or not, and the default value is 1 second, and which can be changed as well. “Implementation round” 344 means that the number of rounds for the multiple running system is implemented currently. “Current time consumption” 345 means the consumed time (seconds) for those rounds of current implementation for this said multiple running system.

Operation methods of implementing the multiple running process for the said far-field release assessment sub-system function are similar to those methods of implementing the multiple running of the said near-field release assessment sub-system function, and it will not be explained herein.

Next, it will be explained that this said invention's parameter sensitivity and uncertainty analysis function, and completely implemented the multiple running of the said near-field release assessment sub-system function and the multiple running of the said far-field release assessment sub-system function, and it can then be conducted the sensitivity and uncertainty analysis for the parameter of this said invention.

Uncertainty and parameter sensitivity analysis is only focused to conduct the analysis on the result that obtained from implemented the multiple running system. When users completed the near-field release assessment sub-system or the multiple running of the said far-field release assessment sub-system process, users can then be used this said system to conduct the uncertainty and parameter sensitivity analysis on abovementioned near-field release assessment sub-system or the multiple running result.

This said invention's uncertainty and parameter sensitivity analysis function is as showed in FIG. 36, and it is mainly included and composed of those following sub-functions, such as “File” 361, “probability analysis” 362, “sensitivity analysis” 363, “Drawing” 364, “Work Directory” 365 and “program verification” 366, and one word display zone 367 and 4 Drawing Zones: 368, 369, 370 and 371. Word display zone 367 is mainly to display data, and the temporary result in the process of regression analysis and the regression equation that obtained from the last regression analysis. 4 Drawing zones 368, 369, 370 and 371 will be showed respectively: (1) near-field release assessment sub-system or multiple running of the said far-field release assessment sub-system result (as showed in FIG. 36's 368). (2) assessed result of the complementary cumulative distribution function (CCDF) (as showed in FIG. 36's 369). (3) assessed result's multiple scatter plot for each parameter (as showed in FIG. 36's 370), and (4) magnified figure of assessed result's scatter plot for certain parameter (as showed in FIG. 36's 371).

If previously completed the implementation of the multiple running function, and it can directly select the “uncertainty and sensitivity analysis” sub-function item of the “sensitivity analysis” system function item to access this said uncertainty and sensitivity analysis function. Operation methods of this said system will be introduced as follows.

After implemented the “File” function in the “uncertainty and sensitivity analysis” sub-function, “File” function item is also included 2 sub-function options, such as Open and Save As . Open function item is include a multiple running figure sub-function item.

After implemented the multiple running figure function item, users can select the file list of multiple running figure to select the file that needs to be analyzed, as showed in FIG. 37. Such figure is also known as the multiple running result figure, and it can be magnified/shrunk as well. Figure's upper left is the word note zone which can input word data. In such figure, it can be seen that it is 5% of the multiple running result (that is, orderly arrange all analytic groups from small to big values, the 5% values). The distributive situation and position for the 50%, 95% and average value curve figure that can be directly selected from the options of the list on figure's left-hand side (as showed in FIG. 37's Run-01), which can use to observe the distributive situation for each round curve.

After implemented the Save As function item, the Save As function item is also included 4 sub-function items, such as the fixed time release rate CCDF data 381, release rate peak CCDF data 382, peak occurrence time CCDF data 383, and percentage total release rate curve 384, as showed in FIG. 38.

After implemented the fixed time release rate CCDF data 381 function, as showed in FIG. 39, and selected the analyzing time point to complete the process of save. Similarly, users can select the release rate peak CCDF data 382 function and complete the save process of file, and after selected the peak occurrence time CCDF data 383 function to complete the file save process.

After implemented the percentage total release rate curve function as showed in FIG. 40, this said system will be automatically set the File name and save the data.

After implemented FIG. 36's “probability analysis”362 function item, as showed in FIG. 41, the function item “probability analysis” is included 3 sub-function items, such as fixed time release rate 4101, release rate peak 4102 and peak occurrence time 4103, and their further function property will be introduced as follows.

After implemented the fixed time release rate 4101 sub-function item, the timetable will be showed in FIG. 39. Users can select the time that needs to carry out the analysis and assessment, and then in the multiple running model to obtain the CCDF figure of the current annual release flow rate (Bq/year), as showed in FIG. 42.

Users implemented the release rate peak 4102 sub-function item that they can obtain the CCDF figure of the release flow rate (Bq/year) peak value for each round in the multiple running model, as showed in FIG. 43.

Users implemented the peak occurrence time 4103 sub-function item that they can obtain the CCDF figure of the release flow rate peak occurrence time (year) for each round in the multiple running model, as showed in FIG. 44.

After implemented FIG. 36's “sensitivity analysis” 363 function item, the “sensitivity analysis” function is included 3 sub-function items, such as the fixed time release rate 441, release rate peak 442 and peak occurrence time 443, as showed in FIG. 45, and their further function property will be introduced as follows.

After implemented the “fixed time release rate” 441 function, such function is also included 3 sub-function items, such as data non-transformed 4411, data Rank transformed 4412 and data Log transformed 4413, as showed in FIG. 46.

After implemented the function item of data non-transformed 4411 will display a time menu, as showed in FIG. 39. After selected the time of carrying out the analysis as showed in FIG. 47, this said system is displayed 3 columns, such as “F VALUE (>=0.01) that parameter has included by regression equation”461, “F VALUE (<=0.009) that parameter has eliminated by regression equation” 462, “tolerance (0.00001˜0.01) obtained by conducting the regression analysis” 463. In addition, the default values that set by this said system “parameter that has included by regression equation when F value (>=0.01)” and the column values is 4.0; for “parameter that has eliminated regression equation when F value (<=0.009)” and the column values is 3.9; for “tolerance (0.00001˜0.01) obtained by conducting the regression analysis” and the column values will be 0.001. When F value is set too high, then the number of parameter that has been included will become less, thus the column values of “parameter that has included by regression equation F value (>=0.01)” column values have to be slightly greater than “parameter that has eliminated by regression equation when F value (<=0.009)” column values are greater their function will be better. The result of implementation is showed in FIG. 48.

Similarly, the data Rank transformed function is identical to the data non-transformed function, and the difference between them is that each will be transformed into Rank in advance, then according to the value of each parameter data in the total data value to code it in the small to big sequence of arrangement. The minimum parameter data value is 1 (Rank=1), and the maximum parameter data value will be the number of sampling, and then transformed into the Rank value and conducted the regression analysis, as showed in FIG. 49.

data Log transformed function is mainly focused on the Log-shape distributive parameter (such as Log Uniform, Log Normal, Log Triangular, etc) to conduct the data Log transformed , and for non-Log-shape distributive parameter (such as Uniform, Normal, Triangular, Gamma, Beta, etc.) the function of data Log transformed , after implemented such function, each parameter data will be obtained the log value in advance, and then used the log values to continuously carry out the analysis on the regression equation, as showed in FIG. 50.

After selected the “release rate peak” 442 function, “release rate peak” function is also included 3 sub-function items: such as data non-transformed data Rank transformed and data Log transformed , its function is identical to abovementioned fixed time release rate 441 function, and it will also not explain herein.

After moved mouse to “peak occurrence time” 443, 3 sub-function items will be displayed: data non-transformed , data Rank transformed and data Log transformed , its function is identical abovementioned fixed time release rate 441, and it will also not explain herein.

After implemented “Drawing” 364 function item in FIG. 36, “Drawing” function is also included 6 sub-function items: modified Y-axis 3641, modified X-axis 3642, drawing scatter plot 3643, display—scatter plot parameter name tag 3644, magnified CCDF 3645 and adding figure 3646. Users shall be noted that if not yet implemented the “sensitivity analysis” function item, and for those functions, such as drawing scatter plot , display—scatter plot parameter name tag and magnified CCDF etc., are unable to use then. The following is introduced the property for such function item.

The “modified Y-axis” function and “modified X-axis” function are identical to these “modified Y-axis” function and “modified X-axis” function in previous other function sub-systems, thus it will not be explained herein.

As for the use of “drawing scatter plot” function item, it has implemented the “sensitivity analysis” function item and after implemented “drawing scatter plot” function item, as showed in FIG. 51. As showed in FIG. 51, each small figure is a scatter plot, which is a small figure that corresponded to an uncertain parameter, and these small figures from the upper left corner to bottom right corner, and from left to right, will be ranged orderly by each uncertain parameter's influential level. Parameters included into the regression line, their scatter plot can be marked by using red regression line, and each parameter's relationship diagram can be magnified in order to facilitate users to observe as showed in FIG. 52.

“Display—scatter plot parameter name tag” function item is to add the scatter plot correspondent parameter name into each small scatter plot, and after added the parameter name into small scatter plot, then the title name of this said function item will be modified to become the “hidden—scatter plot parameter name tag” as showed in FIG. 53; thus, the effect of this said function item is toggled between positive and negative which can add correspondent parameter name or hidden parameter name into each small scatter plot.

“Magnified CCDF” function item is focused on shrinking and magnifying the CCDF, and it is ineffective to other figures; thus, the said function is unable to use if it has not made the CCDF. After selected the “magnified CCDF” function item, the CCDF will be magnified to the full screen, and the title of said function item will be changed into the “shrunk CCDF”, as showed in FIG. 54. Therefore, the effect of this said function item is a continuous circulation with magnifying and shrinking the CCDF.

“Adding figure” function item is the figure of assessed result that adding base case in the multiple running figure for the purpose of comparison. After implemented the “adding figure” function item and the magnified multiple running figure as showed in FIG. 55. In such figure, the X mark curve is adopted the base case to obtain the assessed result of parameter data, and this can be more easily to assess the influence of uncertainty parameter's maximum/minimum value on the assessed result, and whether the setting scope of maximum/minimum values is biased or not.

After selected the “program verification” function item, 3 build-in examples of this said system as showed in FIG. 56. These 3 examples are extracted from the textbooks of Statistics, has standard regression analysis result which can be carried out the comparison. Users are able to use these 3 examples that provided by this said system to verify this program's validity and accuracy for the regression analysis, and its operation methods are identical to abovementioned “sensitivity analysis” function item, thus it will not be explained herein.

From abovementioned detailed introduction, the radioactive waste deep geologic repository performance assessment system disclosed by this said invention can be simplified in order to reduce the difficulty in near-field release assessment and far-field release assessment and any possible man-made error when constructing the radioactive waste deep geologic repository near-field release assessment and far-field release assessment; in addition, it will be helpful to integrate and connect each individual and independent sub-system or external program (such as FORTRAN) to facilitate conduct the safety assessment for the recycled radioactive waste deep geologic repository.

Currently, the radioactive waste deep geologic repository performance assessment system disclosed in this said invention can be calculated from analyzing the nuclide from waste tanks. The tank has broken and it will be released with following the groundwater, through the buffered material of Bentonite, excavation disturbed zone, the disposed geologic host rock, and the diffusion, advection and dispersion effects on the geologic crack to release to biosphere, and it can be analyzed the sequence of influential factors for the near-/far-field release rate.

To sum up, the structural characteristics of this said invention and each actual implementing case has been disclosed in details, and then this said invention can be significantly displayed on its the purpose and efficiency with having great originality and improvement for implementation, which really has the value of industrial usage. This said invention is a unique and exclusive operation and application that ever seen in the current market, according to the spirit of the Patent Act, this said invention case is totally conformed to the important conditions of invention patent.

However, the above mentioned is only the optimal actual case of implementation for this said invention, and cannot be the scope of limiting the implementation scope for this said invention; that is, in most cases will according to this said invention claims to conduct the equal change and modification, and all of such condition will still belong to the coverage scope of patent of this said invention.

Dear review committee member, please give your kind review and approve the application of this said invention.

Claims

1. One radioactive waste deep geologic repository performance assessment system, comprising: One basic data input sub-system is used to set the basic data of the radioactive nuclide that required for the calculation of such system;One parameter setting sub-system is used to set the parameter that required for the calculation of such system;One near-field release assessment sub-system is used to calculate and assess the result of radioactive waste near-field release;One far-field release assessment sub-system is used to calculate and assess the result of radioactive waste far-field release;One biosphere dose assessment sub-system is used to calculate and assess the radioactive waste releasing result of biosphere dose;One multiple running sub-system is used to make multiple running for the near-field release assessment sub-system, such far-field release assessment sub-system and the assessment sub-system of biosphere dose; andOne uncertainty analysis sub-system is used to calculate and assess the uncertainty that obtained from this said system multiple running; andOne sensitivity analysis sub-system is used to calculate and assess the parameter sensitivity that obtained from this said system multiple running.

2. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein said radioactive waste deep geologic repository performance assessment system; among which, such basic data input sub-system will be included one near-field release assessment data input sub-system, one far-field release assessment data input sub-system and the biosphere dose assessment data input sub-system; the input data that used to input the near-field release assessment data input sub-system included: (1). Waste property data contained the data of the inventory time and the time of complete dissolution;(2). Waste tanks property data includes the life-span, erosion density, inner radius, inner radius, length, erosion porosity, and erosion diffusion coefficient;(3). Buffered layer property data contained the data of density, porosity, outer radius and diffusion coefficient;(4). Excavation disturbed zone property data contained the data of rock density, outer radius, porosity, and diffusion coefficient;(5). Host rock property data contained the data of the Darcy flow rate, crack diffusion coefficient, crack spacing and crack opening;(6). Title of released nuclide, half-life and decay phase data;(7). Host nuclide, half-life and sub-nuclide data;(8). Nuclide title, inventory and Instant Release Fraction (IRF) data; and(9). Data that inputted by using such element solubility and sorption coefficient,

3. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein said radioactive waste deep geologic repository performance assessment system; among which, such parameter setting sub-system is included: One fixed parameter sub-system is used to set the fixed parameter; andOne distributive parameter sub-system is used to set the distributive parameter, and such distributive parameter sub-system is included and contained one Latin Hypercube Sampling sub-system that can be used to calculate the Latin Hypercube Sampling, and one Monte Carlo random Sampling sub-system that will be adopted to use to calculate the Monte Carlo Random Sampling.

4. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein the described radioactive waste deep geologic repository performance assessment system; among which such near-field release assessment sub-system is included: The preparation sub-system of data input file in the multiple running of near-field release assessment sub-system, which can be used to input the data file that required for the multiple running of the near-field release assessment sub-system; the preparation sub-system of data input file for multiple running in the near-field release assessment sub-system is included the following 5 sub-systems: (a) to (e):(a). One disposal facility design and the geologic property correlative setting sub-system is used to set the parameter of well-obtained data, the parameter of well-obtained data and adopted the near-field release assessment sub-system to assess the correlative values for the variable sensitivity and the random parameter in a near-field release assessment sub-system;(b). One chemical element solubility correlative setting sub-system is used to set the correlative value of chemical element solubility;(c). One chemical element's sorption coefficient correlative setting sub-system for buffered materials is used to set the correlative value of chemical element in the buffered materials;(d). One chemical element's sorption coefficient correlative setting sub-system in host rock is used to set the chemical element's sorption coefficient correlative value in host rock;(e). One chemical element's sorption coefficient correlative setting sub-system for the erosion in waste tanks is used to set the chemical element's sorption coefficient correlative value for the erosion in waste tanks;One Single Running sub-system is implemented one calculation according to the set parameter;One multiple running sub-system is implemented the multiple running according to the set parameter, and such multiple running sub-system is included 8 sub-systems as the following (f) to (m)=:(f). One parameter of well-obtained data sub-system is used to select the parameter of assessing the variable sensitivity;(g). One parameter of well-obtained data and the correlative setting sub-system which intended to explore the correlative setting sub-system of variable sensitivity by using the near-field release assessment sub-system are used to correlate with the parameter of well-obtained data and used the near-field release assessment sub-system to explore the sensitivity of variable;(h). One near-field release assessment sub-system in random parameter sub-system is used to set the correlative near-field release assessment sub-system of the assessing variable sensitivity parameter;(i). One chemical element solubility setting sub-system is used to set the assessing variable sensitivity for parameter;(j). One chemical element's sorption coefficient setting in sub-system host rock is used to set the parameter of assessing variable sensitivity;(k). One chemical element's sorption coefficient setting sub-system in the buffered materials is used to set the parameter of assessing variable sensitivity;(l). One chemical element's sorption coefficient setting sub-system for the erosion in waste tanks is used to set the assessing for the variable sensitivity;(m). One data arrangement sub-system is used to set the assessed result data arrangement method, and such data arrangement sub-system is also included the following sub-systems: a random arrangement sub-system is arranged by the random method of the assessed result data; a non-correlative arrangement sub-system is arranged by the non-correlative method of the assessed result data; and a specific correlative arrangement sub-system is arranged by the specific correlative method of the assessed result data.

5. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein the described radioactive waste deep geologic repository performance assessment system; among which, such near-field release assessment sub-system is included: One file function sub-system is used to control the data file that applied to such far-field release assessment sub-system, and such file function sub-system is included the following 4 sub-systems from (a) to (c):(a). One Open Old File sub-system is used to open the existing data file, and such Open Old File sub-system is included the following 4 sub-systems from (a1) to (a4): (a1). One nuclide decay chain, half-life, sorption coefficient data setting sub-system is used to select the output data file for the title of released nuclide, half-life and decay chain that generated from implementing the calculation of the near-field release assessment sub-system;(a2). One data of natural barrier system (NBS) property setting sub-system is used to select the existing data of natural barrier system (NBS) property file;(a3). One implementing calculation data setting sub-system is used to select the existed implementing calculation data file;(a4). One complete implementing sub-system is used to input such data setting sub-system, such as the nuclide decay chain, half-life, sorption coefficient, and such data of natural barrier system (NBS) property setting sub-system and the set data file of the said calculation and implementing data setting sub-system;(b). One establishing new data file sub-system is used to establish new data file;(c). One implementing previous data file sub-system is used to calculate and implement the existed data file;One save file function sub-system is used to save the data field of the said far-field release assessment sub-system, and the said save file function sub-system is included the following 4 sub-systems from (e) to (h):(e) Data save sub-system of the nuclide decay chain, half-life, sorption coefficient is used to save the data of nuclide decay chain, half-life and sorption coefficient;(f) Natural barrier system (NBS) property save sub-system is used to save the data of natural barrier system (NBS) property;(g) Program implementing the setting data save sub-system is used to save the setting data of program implementation;(h) Complete implementing case save sub-system is used to save the above mentioned (e), (f), and (g) data at a time;One insert file function sub-system is used to insert other nuclide items to connect into the new nuclide data content;One clear file function sub-system is used to clear the data of the nuclide decay chain, half-life, sorption coefficient, decay chain, element's sorption coefficient in the host rock;One review file function sub-system is used the data of the review program implementing data input file, program implementing output file, and the program implementing output explanatory file;One drawing function sub-system is used to display the result of implementing and calculating the said far-field release assessment sub-system in the form of figure, and the drawing function sub-system is included the following 5 sub-systems:One previously implemented case sub-system is used to display the data and figure of the output result from the previously implemented far-field release assessment sub-system;One previously implemented case sub-system is used to display the correlative result and figure of the selected file;One modified Y-axis sub-system is used to modify and display the maximum and minimum values for Y-axis value in the figure;One modified X-axis sub-system is used to modify and display the maximum and minimum values for X-axis value in the figure;One adding figure sub-system is used to stack the output figure of the selected file onto another figure in order to display different output results in a same figure;One Work Directory sub-system is used to display the current Work Directory for the said system, and replace the said system's Work Directory as well as establish new Work Directory for the said system.

6. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein said radioactive waste deep geologic repository performance assessment system, the far-field release assessment sub-system is included: One geometry property setting sub-system is used to set the geometry property data that required for the calculation of the said far-field release assessment sub-system;One Host Rock property setting sub-system is used to set the Host Rock property data that required for the calculation of the said far-field release assessment sub-system;One Transport property setting sub-system is used to set the transportation property data that required for the calculation of the said far-field release assessment sub-system;One Nuclide Flux Input File setting sub-system is used to set the Nuclide Flux Input File setting data that required for the calculation of the said far-field release assessment sub-system;One Nuclide Concentration Output Time setting sub-system is used to set the Nuclide Concentration Output Time setting data that required for the calculation of the said far-field release assessment sub-system;One assessing the release rate for the nuclide decay chain setting sub-system is used to set the nuclide decay chain data that required for the calculation of the said far-field release assessment sub-system; andOne element's sorption coefficient in the host rock setting sub-system is used to set the element's sorption coefficient in the host rock that required for the calculation of the said far-field release assessment sub-system.

7. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein the described radioactive waste deep geologic repository performance assessment system, the said far-field release assessment sub-system is included: One file function sub-system is used to control the data file that applied to the said far-field release assessment sub-system, and the said file function sub-system is included the following 4 sub-systems from (a) to (d):(a). One Open Old File sub-system is used to open the existed data file and the said Open Old File sub-system is included the following 4 sub-systems from (a1) to (a4): (a1). One nuclide decay chain, half-life, sorption coefficient data setting sub-system is used to select the output data of the title of released nuclide, half-life, and decay chain that generated from the calculation and implementing the said near-field release assessment sub-system;(a2). One data of natural barrier system (NBS) property setting sub-system is used to select the existed data of natural barrier system (NBS) property file;(a3). One calculation and implementing data setting sub-system is used to select the existed calculation and implementing data file;(a4). One complete implementing sub-system is used to integrally input the data setting sub-system for nuclide decay chain, half-life and sorption coefficient, the data of natural barrier system (NBS) property setting sub-system and the data file set by the calculation and implementing data setting sub-system;(b). One establishing new data file sub-system is used to establish the new data file;(c). One implementing previous data file sub-system is used to calculate and implement the existed data file;(d). One File Processing sub-system is used to control file, and the said File Processing sub-system is included one file mergence sub-system that is used to merge and merging files; one file name changing sub-system is used to change the title of file; and one file delete sub-system is used to delete the file that needs to be deleted:One save file function sub-system is used to save the data file that applied to the said far-field release assessment sub-system, and the save file function sub-system is included the following 4 sub-systems from (e) to (h):(e) Data save sub-system of the nuclide decay chain, half-life, sorption coefficient is used to save the data of nuclide decay chain, half-life and sorption coefficient;(f) Data of natural barrier system (NBS) property save sub-system is used to save the data of natural barrier system (NBS) property;(g) Program implementation setting data save sub-system is used to save the setting data of program implementation;(h) Complete implementing case save sub-system is used to save these abovementioned data of (e), (f), and (g) at a time;One insert file function sub-system is used to insert other nuclide items to connect into a new content of nuclide data;One Clear file function sub-system is used to clear the data of the nuclide decay chain, half-life, sorption coefficient, decay chain, element's sorption coefficient in the host rock;One Review file function sub-system is used to review the data of the program implementing data input file, program implementing output file and program implementing output explanatory file;One Drawing function sub-system is used to display the result of calculation and implementation of the said far-field release assessment sub-system in the form of figure, and the said Drawing function sub-system is included the following 5 sub-systems:One previously implemented case sub-system is used to display the data and figure of the output result from the previously implemented far-field release assessment sub-system;One previously implemented case sub-system is used to display the correlative result and figure of the selected file;One modified Y-axis sub-system is used to modify and display the maximum and minimum values for Y-axis value in the figure;One modified X-axis sub-system is used to modify and display the maximum and minimum values for X-axis value in the figure;One adding figure sub-system is used to stack the output figure of the selected file onto another figure in order to display different output results in a same figure;One Work Directory sub-system is used to display the current Work Directory for the said system, and replace the said system's Work Directory as well as establish new Work Directory for the said system.

8. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein the described radioactive waste deep geologic repository performance assessment system, the said far-field release assessment sub-system is included: The preparation sub-system of data input file for the multiple running of one far-field release assessment sub-system is used to input the data file that required for the multiple running of the said far-field release assessment sub-system; the preparation sub-system of data input file for the multiple running of the said far-field release assessment sub-system is included the following 2 sub-systems from (a) to (b):(a). One disposal facility design and geologic property correlative setting sub-system is used to set the parameter of well-obtained data, parameter of well-obtained data; adopting one near-field release assessment sub-system to assess the correlative value of variable sensitivity; and the random parameter of the said near-field release assessment sub-system;(b). One chemical element's sorption coefficient in the host rock setting sub-system is used to set the correlative value of chemical element's sorption coefficient in the host rock;One Single Running sub-system is used to implement a single calculation according to the set parameter;One multiple running sub-system is used to implement the multiple running according to the set parameter, and the said multiple running sub-system is included the following 5 sub-systems from (c) to (j):(c). One parameter of well-obtained data sub-system is used to select the parameter of assessing the variable sensitivity;(g). One parameter of well-obtained data and using the said far-field release assessment sub-system to explore the correlative setting sub-system for the variable sensitivity is applied the correlative parameter of well-obtained data and using the said far-field release assessment sub-system to explore the sensitivity of variable;(h). One random parameter sub-system of the said far-field release assessment sub-system is used to set the correlative near-field release assessment sub-system parameter of assessing variable sensitivity parameter;(i). One chemical element's sorption coefficient setting sub-system in the host rock is used to set the parameter of assessing the variable sensitivity;(j). One data arrangement sub-system is used to set the method of assessed result data arrangement, and the said data arrangement sub-system is included the following 3 sub-systems: a random arrangement sub-system is used to arrange the data of assessed result in a random method; a non-correlative arrangement sub-system is used to arrange the data of assessed resulting a non-correlative method; and a specific correlative arrangement sub-system is used to arrange the data of assessed result in a specific correlative method.

9. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein the described radioactive waste deep geologic repository performance assessment system, the said multiple running sub-system is include: One file control sub-system is used to provide the functions of file mergence, file name change and file delete, etc.;One Drawing sub-system is used to display the calculation result of the said near-/far-field release assessment sub-system in the form of figure, and the said Drawing function sub-system is included the following 3 sub-systems:One previously implemented case sub-system is used to display the correlative result and figure of the selected file;One modified Y-axis sub-system is used to modify and display the maximum and minimum values for Y-axis value in the figure;One modified X-axis sub-system is used to modify and display the maximum and minimum values for X-axis value in the figure;One Work Directory sub-system is used to display the current Work Directory for the said system, and replace the said system's Work Directory as well as establish new Work Directory for the said system.One multiple running sub-system of the said near-field release assessment sub-system is used to select the file data that established by the preparation system of data input file for the multiple running of the said near-field release assessment sub-system to conduct the multiple running for the said near-field release assessment sub-system, and is used to set the values of final implementation round and nuclide number;One multiple running sub-system of the said far-field release assessment sub-system is used to select the file data that established by the preparation system of data input file for the multiple running of the said far-field release assessment sub-system to conduct the multiple running for the said far-field release assessment sub-system, and is used to set the values of final implementation round and nuclide number;One suspended sub-system is used to suspend the implementation of calculation for the said near-field release assessment sub-system or far-field release assessment sub-system.One display function sub-system is used to display 4 items: set time consumption, implementation round, current time consumption and time increment; among which, the set time consumption is indicated that after the set time (unit is second), the said multiple running system will start to check whether the called program is completely implemented or not, and the default value is 1(second); the time increment is indicated that since the set time began, the said system will check in a certain time interval for whether the called program is completely implemented or not, and the default value is 1 second; the implementation round is indicated that the current number of implanting rounds for the said multiple running system; and the current time consumption is indicated that the time (seconds) consumed of the current implementing rounds for the multiple running system.

10. One radioactive waste deep geologic repository performance assessment system as recited in claim 1, wherein the described radioactive waste deep geologic repository performance assessment system, the said uncertainty analysis sub-system and sensitivity analysis sub-system will be included: One file control sub-system is used to control file, the said file control sub-system is included one open sub-system and one save sub-system. The said open sub-system is used to open a multiple running figure; and the said save sub-system is included the following 4 sub-systems; fixed time release rate CCDF data sub-system; release rate peak CCDF data sub-system; peak occurrence time CCDF data sub-system; and percentage total release rate curve sub-system. Fixed time release rate CCDF data sub-system is used to select the assessing time points to input the File name and complete the process of date save. Release rate peak CCDF data sub-system is used to input the File name and complete the process of file save. Peak occurrence time CCDF data sub-system is used to input the File name and complete the process of file save. Percentage total release rate curve sub-system is used to automatically set the File name and save the data;One probability analysis sub-system is used to conduct the calculation of the probability analysis, and the said sub-system is included the following 3 sub-systems: fixed time release rate sub-system, release rate peak sub-system and peak occurrence time sub-system, among which, fixed time release rate is used to select the time that needs to be analyzed and assessed, and can be outputted a CCDF figure of an annual release flow rate (Bq/year) in the multiple running model. Release rate peak sub-system is used to output the CCDF figure of the peak release flow rate (Bq/year) for each round in a multiple running model. Peak occurrence time sub-system is used to output the CCDF figure of the peak occurrence time (year) for each round's release flow rate in a multiple running model;One sensitivity analysis sub-system is used to conduct the calculation of the sensitivity analysis, and the said sub-system is included the following 3 sub-systems from (a) to (c): fixed time release rate sub-system, release rate peak sub-system and peak occurrence time sub-system, among which,(a) the said fixed time release rate sub-system is included the following 3 sub-systems: data non-transformed sub-system, data Rank transformed sub-system and data Log-transformed sub-system, among which, the said data non-transformed sub-system is used to select the time that needs to be analyzed, and its selectable values are the F VALUE (>=0.01) that parameter has included by regression equation; the F VALUE (<=0.009) that parameter has eliminated by regression equation; and the tolerance (0.00001˜0.01) that obtained by conducting the regression analysis. Each parameter data in the said data Rank transformed sub-system will be transformed into Rank in advance, then according to the value of each parameter data in the total data value to code it in the small to big sequence of arrangement. The minimum parameter data value is 1 (Rank=1), and the maximum parameter data value will be the number of sampling, and then transformed into the Rank value and conducted the regression analysis. Each parameter data in the data Log-transformed sub-system will be obtained the logarithmic value in advance, next to obtain the value after conducting the log-transformed process, and then continuously conduct the regression analysis;(b) The said release rate peak sub-system is included the following 3 sub-systems: data non-transformed sub-system, data Rank transformed sub-system and data Log-transformed sub-system, among which, the data non-transformed sub-system is used to select the time that needs to be analyzed, and its selectable values are the F VALUE (>=0.01) that parameter has included by regression equation; the F VALUE (<=0.009) that parameter has eliminated by regression equation; and the tolerance (0.00001˜0.01) that obtained by conducting the regression analysis. Each parameter data in the said data Rank transformed sub-system will be transformed into Rank in advance, then according to the value of each parameter data in the total data value to code it in the small to big sequence of arrangement. The minimum parameter data value is 1 (Rank=1), and the maximum parameter data value will be the number of sampling, and then transformed into the Rank value and conducted the regression analysis. Each parameter data in the data Log-transformed sub-system will be obtained the logarithmic value in advance, next to obtain the value after conducted the Log-transformed process, and then continuously conduct the regression analysis;(c) The said peak occurrence time sub-system is included the following 3 sub-systems: data non-transformed sub-system, data Rank transformed sub-system and data Logo transformed sub-system, among which, the said data non-transformed sub-system is used to select the time that needs to be analyzed, and its selectable values are the F VALUE (>=0.01) that parameter has included by regression equation; the F VALUE (<=0.009) that parameter has eliminated by regression equation; and the tolerance (0.00001˜0.01) that obtained by conducting the regression analysis. Each parameter data in the said data Rank transformed sub-system will be transformed into Rank in advance, then according to the value of each parameter data in the total data value to code it in the small to big sequence of arrangement. The minimum parameter data value is 1 (Rank=1), and the maximum parameter data value will be the number of sampling, and then transformed into the Rank value and conducted the regression analysis. Each parameter data in the data Log-transformed sub-system will be obtained the logarithmic value in advance, next to obtain the value after conducting the Log-transformed process, and then continuously conduct the regression analysis;One Work Directory control sub-system is used to control the Work Directory that required for the calculation of the said system;One program verification program verification sub-system is used to verify the accuracy of calculated result;One text display sub-system is used to display the data, the temporary result in the process of regression analysis, and the regression value that obtained from the final regression analysis;One Drawing sub-system is used to display the result of multiple running for the said near-field release assessment sub-system or far-field release assessment sub-system multiple running, the result of complementary cumulative distribution function (CCDF) for the assessed result, and the multiple scatter plot of the assessed result to individual parameter and (4) the magnified scatter plot of the assessed result to individual parameter. The said Drawing sub-system is included the following 6 sub-systems from (a) to (f): modified Y-axis sub-system; modified X-axis sub-system; Drawing scatter plot sub-system; display scatter plot parameter name tag sub-system; magnified CCDF sub-system; and adding figure sub-system, among which,(a) The said modified Y-axis sub-system is included the following 2 sub-systems: the maximum value sub-system and the minimum value sub-system, among which, the maximum value sub-system can be used to modify the maximum value of Y-axis in the multiple running analysis figure; the minimum value sub-system can be used to modify the minimum value of Y-axis in the multiple running analysis figure;(b) The said modified X-axis sub-system is included the following 2 sub-systems: the maximum value sub-system and the minimum value sub-system, among which, the maximum value sub-system can be used to modify the maximum value of X-axis in the multiple running analysis figure; the minimum value sub-system can be used to modify the minimum value of X-axis in the multiple running analysis figure;(c) The said Drawing scatter plot sub-system is used to draw and display the scatter plot;(d) The said display scatter plot parameter name tag sub-system is used to add the correspondent title of the parameter to the said scatter plot for each scatter plot;(e) The said magnified/shrunk CCDF sub-system is used to magnify or shrink CCDF figure;(f) The said adding figure sub-system is used to add the base case into the multiple running figure for the said assessed result figure on the purpose of conducting the comparative analysis.