CN106934082A

CN106934082A - Hlw Geological Repository site suitability evaluation methods based on Delphi method

Info

Publication number: CN106934082A
Application number: CN201511020305.XA
Authority: CN
Inventors: 苏锐; 周志超
Original assignee: Beijing Research Institute of Uranium Geology
Current assignee: Beijing Research Institute of Uranium Geology
Priority date: 2015-12-30
Filing date: 2015-12-30
Publication date: 2017-07-07

Abstract

Field is administered the invention belongs to the elimination of nuclear facilities and radwaste, and in particular to a kind of Hlw Geological Repository site suitability evaluation methods based on Delphi method.The method of the present invention is comprised the following steps：Step 1. determines the index of Hlw Geological Repository site suitability quantitative assessment；Step 2. builds the Codes and Standards of indices quantification；Step 3. weight of Delphi method agriculture products；Step 4. calculates the Hlw Geological Repository site final score value of suitability.The present invention is to set up based on evaluation index, carry out agriculture products weight by means of Delphi method, it is core with the Codes and Standards for building indices quantification, solve the problem of Hlw Geological Repository site suitability quantitative assessment, it is adaptable to which Hlw Geological Repository etc. is large-scale, complicated and site suitability quantitative evaluation of safety evaluation cycle overlength.

Description

Hlw Geological Repository site suitability evaluation methods based on Delphi method

Technical field

Field is administered the invention belongs to the elimination of nuclear facilities and radwaste, and in particular to a kind of height based on Delphi method Put waste geological disposal storehouse site suitability evaluation methods.

Background technology

High radioactivity nuclear waste, abbreviation high-level waste, its radioactivity is strong, toxicity big, heating, and nucleic half-life period is extremely long, Must take reliable means that it is effectively isolated with environment for human survival.From in terms of current world development situation, using deep geology The geological disposal mode buried (in the deep about 500-1000 meters geologic body of earth's surface) and isolate with biosphere is considered as a kind of Effectively, safe disposal means, it has also become international consensus.The safe disposal of high-level waste is in science, technology, humanity and sociology Aspect faces a series of great, complicated problems, and its difficult point is that how to realize that high-level waste is carried out with biosphere abundant, thorough Bottom, reliably isolate；How to realize isolation time more than 10,000 years, even 100,000 years；How the reliability of isolation is effectively predicted Deng.Wherein, how to select qualified site, how to evaluate site suitability be solve high-level waste safe disposal section Learn one of technical barrier.

China《The Law on the Prevention and Control of Radioactive Pollution》" high-level radioactive solid waste carries out the deep geology concentrated to clear stipulaties Disposal ".According to State Council's announcement in 2011《Radwaste safety management regulations》, China mainland area is low, reclaimed water keeps flat Penetrating property solid waste disposal facility should meet the security isolation requirement of more than 300 years after closing；High-level radioactive solid waste The security isolation requirement of more than 10,000 years should be met after being closed with alpha activity solid waste depth geological disposal facility.Thus, Gao Fang Waste geological disposal storehouse engineering is also referred to as 10000 years engineerings, and the screening of its site and suitability evaluation be both different from general civilian work Journey, also different from existing all nuclear engineerings.

China's Hlw Geological Repository sieving disposal site and appraisal are started from 1985, have tentatively screened 6 in the whole nation Pre-selected zone and 21 key sections and 7 candidate sites, but sieving disposal site and appraisal greatly are still within whole height and put useless The initial stage of thing geological disposal research and development.Either pre-selected zone screening or the screening of key sections and candidate site, adopt Method is expert's qualitative subjective evaluation.

Quantitative evaluation is the information age, particularly the important development direction of big data epoch management, extensive use In natural science and Social Science Research, for science decision provides support.In view of Hlw Geological Repository engineering has The features such as construction condition complexity, safe class high, period of service limit for length (more than 10,000 years), disposal storehouse site suitability quantification is commented Valency is the inevitable requirement of site decision-making during high-level waste geology treatment research and development.

Delphi method was initiated by Hull silent (Helmer) and Gordon (Gordon) in the forties in 20th century, nineteen forty-six, beautiful Rand Corporation of state be avoid brainstorming exist submit to authority or obey blindly majority defects, first in this way Qualitative forecasting is carried out, later the method is by rapid widely used.As it was previously stated, Delphi method is only in itself a Weight prediction side Method, can not be carried out site suitability and quantified in itself by expert according to the weight for specifying program investigation index, therefore Delphi method Evaluate, only suitability quantitative assessment provides indispensable weight parameter and is input into.In other words, it is fixed for site suitability For amount is evaluated, escheat of only having the right is inadequate, it is necessary to which combining with quantitative evaluation index can just carry out quantitative evaluation.

How for Hlw Geological Repository site suitability carries out quantitative evaluation to be nuclear facilities using Delphi method The retired technical problem with radwaste Treatment process field urgent need to resolve.

The content of the invention

The technical problem to be solved in the invention is：Hlw Geological Repository site suitability is difficult to quantitatively comment Valency.

Technical scheme is as described below：

A kind of Hlw Geological Repository site suitability evaluation methods based on Delphi method, comprise the following steps： Step 1. determines the index of Hlw Geological Repository site suitability quantitative assessment；Step 2. builds the rule of indices quantification Then and standard；Step 3. weight of Delphi method agriculture products；It is suitable that step 4. calculates Hlw Geological Repository site The final score value of property.

Index includes described in step 1：Site Selection stage index, Regional survey stage index and site are evaluated and true Recognize stage index.

In step 2, Site Selection stage index suitability scale, Regional survey rank are represented respectively using PQV, QQV and SQV Duan Zhibiao suitabilities scale and site are evaluated and the stage of recognition index suitability scale；

PQV_i1 ^j1Represent the suitability scale of 1 two-level index of jth in the i-th 1 first class index of Site Selection stage； QQV_i2 ^j2Represent the suitability scale of 2 two-level index of jth in the i-th 2 first class index of Regional survey stage；SQV_i3 ^j3Represent field Location evaluate and the i-th 3 first class index of the stage of recognition in 3 two-level index of jth suitability scale；

Step 2.1. Site Selection stages each index suitability scale value

If the suitability scale of i-th first class index of Site Selection stage is PQV_i, wherein i=1~5；

I, geological conditions first class index PQV₁

(1) candidate's country rock thickness PRT

Value rules and methods：If argillaceous rocks and the continuously distributed thickness of rock salt geologic body are less than 100m, other rock classes are thick Degree is less than 500m, then suitability scale PQV₁ ¹=0；If argillaceous rocks and the continuously distributed thickness of rock salt geologic body are less than 100m, its Its rock class thickness is less than 500m, suitability scale PQV₁ ¹=10；

(2) candidate's country rock buried depth PRD

Value rules and methods：If candidate's country rock buried depth meets underground 500m~1000m depth boundses, PQV₁ ²=10； Otherwise, PQV₁ ²=0；

(3) candidate's country rock formation dip PFA

Value rules and methods：For non-laminar class candidate's Rock Mass, PQV₁ ³=10；For stratiform class geologic body, when PFA<At 5 °, PQV₁ ³=10, as 5 °≤PFA<At 10 °, PQV₁ ³=5, other situations PQV₁ ³=0；

II, following natural trend condition first class index PQV₂

(1) the dynamic peak accelerator PSL of regional earthquake

Value rules and methods：If PSL>0.2g, PQV₂ ¹=0；If 0.15g<PSL≤0.2g, PQV₂ ¹=1；If 0.1g<PSL≤0.15g, PQV₂ ¹=4；If 0.05g<PSL≤0.1g, PQV₂ ¹=7；If PSL≤0.05g, PQV₂ ¹=10；

(2) to Regional Geotectonic secondary unit border apart from PGZ

Value rules and methods：If pre-selected zone border is to Regional Geotectonic secondary unit border apart from PGZ< 25km, PQV₂ ²=0, i.e., should reduce the scope of pre-selected zone；If PGZ>100km, PQV₂ ²=10；If 25km≤PGZ≤ 100km, PQV₂ ²Determined according to equipartition principle and linear interpolation method；

(3) to the strong zone of influence of extreme geologic process apart from PGA

Value rules and methods：If PGA<25km, PQV₂ ³=0；If PGA>100km, PQV₂ ³=10；If 25km ≤ PGA≤100km, PQV₂ ³Determined according to equipartition principle and linear interpolation method；

(4) natural geologic hazard development degree PND

Value rules and methods：Natural geologic hazard development degree is divided into strong, medium, weak and four grades of agensis, The corresponding suitability scale PQV of this four grades₂ ⁴It is followed successively by 0,2,7,10；

III, geochemical conditions first class index PQV₃

(1) the main surrounding rock type PTR in pre-selected zone

Value rules and methods：When PTR is this four rocks of CRYSTALLINE ROCKS, argillaceous rocks, tufa stone and rock salt, PQV₃ ¹= 10；Otherwise, PQV₃ ¹=0；

IV, environmental protection condition first class index PQV₄

(1) to Head streams protection zone apart from PSS

Value rules and methods：If PSS=0, PQV₄ ¹=0；If PSS>80km, then PQV₄ ¹=10；If 0< PSS≤80km, according to equipartition principle and linear interpolation principle value；

(2) the quantity PEP of environmental sensitive area above the provincial level

Value rules and methods：Count with preselect district center in round dot 80km radius laws and regulations specify it is provincial And the quantity of national environmental sensitive area, if PEP<20, then PQV₄ ²=10；If PEP>500, then PQV₄ ²=1；Such as Really 20≤PEP≤500, are determined by equipartition principle and linear interpolation method；(3) the quantity PWP of important water body

Value rules and methods：Statistics with preselect district center as in round dot 80km radius in downstream direction laws and regulations The quantity of the provincial and national important water body of regulation, if PWP<10, then PQV₄ ³=10；If PWP>100, then PQV₄ ³=1；If 10≤PWP≤100, determined by equipartition principle and linear interpolation method；(4) pre-selected zone absolute altitude is high Degree PHE

Value rules and methods：If the height above sea level of pre-selected zone highest point is more than 3500m, PQV₄ ⁴=0；Otherwise, PQV₄ ⁴ =10；

V, social economy and humane condition first class index PQV₅

(1) the absolute density of population PAP in region

Value rules and methods：The absolute density of population in region refers to and resided in pre-selected zone population and floating population's summation and preselect The area ratio in area, if PAP >=400 people/km², PQV₅ ¹=0；If 200 people/km²≤PAP<400 people/km², PQV₅ ¹=1； If 100 people/km²≤PAP<200 people/km², PQV₅ ¹=3；If 20 people/km²≤PAP<100 people/km², PQV₅ ¹=6；If PAP<20 people/km², PQV₅ ¹=10；

(2) Main Towns quantity PTW

Value rules and methods：Main Towns quantity refers to medium and above scale cities and towns the quantity included in pre-selected zone, If PTW≤5, PQV₅ ²=10；If 5<PTW≤50, PQV₅ ²=7；If 50<PTW≤100, PQV₅ ²=3； If PTW>100, PQV₅ ²=1；

(3) region total output value PDP per capita

Value rules and methods：Provincial administrative area gives birth to per capita where this phase zone per capita total output value pre-selected zone Produce total value to characterize, if pre-selected zone is averagely added up afterwards across multiple provincial administrative areas according to area；Assuming that the whole nation produces always per capita Value PDP_nRepresent, Provincial administrative area per capita total output value come top 10 total output value per capita average value PDP₁₀ Represent, if PDP≤PDP_n, PQV₅ ³=10；If PDP >=PDP₁₀, PQV₅ ³=1；In the case of other, PQV₅ ³=5；

Step 2.2. Regional survey stages each index suitability scale value

If the suitability scale QQV of i-th first class index of Regional survey stage_iRepresent, wherein i=1~5；

1st, geological conditions first class index QQV₁

(1) candidate's country rock area QGA

Value rules and methods：If QGA<50km², QQV₁ ¹=1；If QGA>300km², then QQV₁ ¹=10；If 50km²≤QGA≤300km², QQV₁ ¹According to equipartition principle and linear interpolation principle value；(2) candidate's country rock thickness QRT

Value rules and methods：If QRT<100m, QQV₁ ²=0；If QRT>500m, QQV₁ ²=10；If 100m≤ QRT≤500m, QQV₁ ²According to equipartition principle and linear principle value；

(3) candidate's country rock buried depth QRD

Value rules and methods：If candidate's country rock buried depth meets underground 500m~1000m depth boundses, QQV₁ ³=10； Otherwise, QQV₁ ³=0；

(4) rift structure distribution surface density QSD

Value rules and methods：Rift structure refers to the rift structure for cutting and passing through pre-selection location, rift structure distribution face Density refers to the rift structure quantity in unit area under the conditions of uniform or more uniform distribution；If QSD<1/km², QQV₁ ⁴= 10；If QSD>10/km², QQV₁ ⁴=1；If 1/km²≤ QSD≤100/km², QQV₁ ⁴According to equipartition principle and line Property principle value；II, hydrogeologic condition first class index QQV₂

(1) groundwater level depth QWB

Value rules and methods：Assuming that the buried depth in proposed disposal storehouse is MZ, if QWB≤MZ, QQV₂ ¹=5；Otherwise, QQV₂ ¹=10；

(2) underground water hydraulic gradient QHI

Value rules and methods：If QHI<1 ‰, QQV₂ ²=10；If QHI>100 ‰, QQV₂ ²=0；If 1 ‰≤ QHI≤100 ‰, QQV₂ ²According to equipartition principle and linear principle value；(3) rock mass ground infiltration coefficient accounting QK

Value rules and methods：Using mathematical statistics method, with rock mass hyposmosis coefficient≤1 × 10^-8M/s proportions are big Small determination QQV₂ ³Value；When hyposmosis coefficient≤1 × 10 of rock mass in certain location^-8Accounting P { QK │ QK≤1 × 10 of m/s^-8m/ S }=1 when, QQV₂ ³=10；When 0.01<P{QK│QK≤1×10^-8M/s } ＜ 1 when, QQV₂ ³=10 × P；When P QK │ QK≤1 × 10^-8M/s }≤0.01 when, QQV₂ ³=0.1；

(4) to region discharge area apart from QPD

Value rules and methods：If QPD<50km, QQV₂ ⁴=1；If QPD>200km, QQV₂ ⁴=10；If 50km ≤ QPD≤200km, QQV₂ ⁴By equipartition principle and linear principle value；

(5) the transmissibility QTF of rift structure

Value rules and methods：Important Faults construction refers to the rift structure for cutting and passing through rock mass；In this stage, according to disconnected Split the property and charges feature qualitative discrimination of construction, be divided into water transmitting ability it is poor, in, strong three kinds of situations, its corresponding suitability Scale QQV₂ ⁵Value is followed successively by 10,1,0.1；

III, geochemical conditions first class index QQV₃

(1) underground water TDS, i.e. QTW

Value rules and methods：If QTW>100g/l, QQV₃ ¹=1；If QTW<1g/l, QQV₃ ¹=10；If 1g/l ≤ QTW≤100g/l, QQV₃ ¹Determined according to equipartition principle and linear interpolation method；(2) underground water Eh values, i.e. QEH

Value rules and methods：If QEH=0, QQV₃ ²=6；If QEH<0, QQV₃ ²=10；If QEH>0, QQV₃ ² =1；

(3) groundwater age QAW

Value rules and methods：If QAW≤30 year, QQV₃ ³=0；If 30 years<QAW≤1000 year, QQV₃ ³=3；Such as Really 1000<QAW≤10000 year, QQV₃ ³=6；If QAW>, QQV in 10000₃ ³=10；

(4) sulfate concentration QSC in underground water

Value rules and methods：If QSC<500mg/L, QQV₃ ⁴=10；If QSC>5000mg/L, QQV₃ ⁴=1；Such as Fruit 500mg/L≤QSC≤5000mg/L, QQV₃ ⁴Determine according to equipartition principle and linear interpolation method；(5) chlorion is dense in underground water Degree QCC

Value rules and methods：If QCC<1000mg/l, QQV₃ ⁵=10；If QCC>5000mg/l, QQV₃ ⁵=1；Such as Fruit 1000mg/l≤QCC≤5000mg/l, QQV₃ ⁵Determine according to equipartition principle and linear interpolation method；IV, environmental protection condition one Level index QQV₄

(1) the quantity QEP of environmental sensitive area above the provincial level

Value rules and methods：Count with Section Center in round dot 25km radius laws and regulations specify it is provincial and The quantity of national environmental sensitive area；If QEP<5, QQV₄ ¹=10；If QEP>50, QQV₄ ¹=1；5≤QEP≤ 50, QQV₄ ¹Determine by equipartition principle and linear interpolation method；(2) the quantity QWP of important water body

Value rules and methods：Statistics is to preselect the province that Section Center specifies as laws and regulations in round dot 25km radius The quantity of level and national important water body；If QWP<5, QQV₄ ²=10；If QWP>20, QQV₄ ²=1；If 5≤ QWP≤20, QQV₄ ²Determine by equipartition principle and linear interpolation method；V, social economy and humane condition first class index QQV₅

(1) the absolute density of population QAP in region

Value rules and methods：The absolute density of population QAP in region refers in the range of with Section Center as center of circle radius 15km Resident population and floating population's summation and area ratio, QAP is more big more is not suitable for addressing；If QAP<5 people/km², QQV₅ ¹= 10；If QAP>200 people/km², QQV₅ ¹=0；If 5 people/km²≤ QAP≤200 people/km², QQV₅ ¹By equipartition principle and linearly Interpolation method determines；(2) Main Towns quantity QTW

Value rules and methods：Main Towns quantity QTW refers to medium in the range of with Section Center as center of circle radius 15km And the quantity in above scale cities and towns；If QTW>50, QQV₅ ²=1；If 10<QTW≤50, QQV₅ ²=3；If 3<QTW≤ 10, then QQV₅ ²=6；If QTW≤3, QQV₅ ²=10；

(3) total output value QDP per capita

Value rules and methods：This stage, total output value QDP referred to giving birth to per capita for local administrative area where pre-selection location per capita Total value is produced, if pre-selection location is averagely added up afterwards across multiple local administrative areas according to area；Assuming that national local administrative area The average value of total output value is QDP per capita_x, total output value comes giving birth to per capita for inverse the 50th per capita in national local administrative area Produce total value QDP_-50Represent；If QDP>2QDP_x, QQV₅ ³=1；If QDP_x<QDP≤2QDP_x, QQV₅ ³=3；If QDP_-50 <QDP≤QDP_x, QQV₅ ³=6；If QDP≤QDP_-50, QQV₅ ³=10；

Evaluate and each index suitability scale value of the stage of recognition step 2.3. sites

If the score value SQV with i-th first class index of the stage of recognition is evaluated in site_iRepresent, wherein i=1~10；

1st, geological conditions first class index SQV₁

(1) candidate's country rock area SGA

Value rules and methods：If candidate's country rock area SGA<10km², SQV₁ ¹=1；If 10km²≤SGA<20km², SQV₁ ¹=3；If 20km²≤SGA<50km², SQV₁ ¹=6；If SGA >=50km², SQV₁ ¹=10；

(2) candidate's country rock thickness SRT

Value rules and methods：If candidate's country rock thickness SRT<100m, SQV₁ ²=0；If 100m≤SRT<200m, SQV₁ ²=3；If 200m≤SRT<400m, SQV₁ ²=6；If SRT >=400m, SQV₁ ²=10；

(3) rock thermal conductivity factor SCT

Value rules and methods：The SQV of SCT the maximum in candidate site₁ ³=10, SCT recklings in candidate site SQV₁ ³=1, the SQV of other candidate sites₁ ³Determine by equipartition principle and linear interpolation method；

(4) rift structure distribution surface density SSD

Value rules and methods：This stage only counts the rift structure for cutting and passing through candidate site, if rift structure Distribution surface density SSD>5/km², SQV₁ ⁴=1；If 5/km²≥SSD>1/km², SQV₁ ⁴=3；If 1/km²≥SSD>0.2/ km², SQV₁ ⁴=6；If SSD≤0.2/km², SQV₁ ⁴=10；

(5) to the Weighted distance SFD of Important Faults construction

Value rules and methods：So-called Important Faults construction refers to the rift structure for cutting and having passed through pre-selection location, herein Stage scope of statistics is 25km scopes around candidate site；Assuming that having Important Faults in the range of 25km around i-th candidate site Construction m, its distance DS_jRepresent；The weight that every Important Faults are constructed from the close-by examples to those far off is given according to apart from site, k is used_jTable Show, apart from DS_jWith weight k_jRelation be shown below：

If apart from k_jInterior Important Faults construction quantity Y_kRepresent, wherein k is the interval quantity of distance, then i-th candidate The Weighted distance index S FD that location constructs to Important Faults_iIt is expressed as follows：

For all candidate sites, the SQV of SFD the maximum₁ ⁵The SQV of=10, SFD reckling₁ ⁵=1, other sites SQV₁ ⁵Determine by equipartition principle and linear interpolation method；

II, following natural trend condition first class index SQV₂

(1) strong geologic process probability of happening SPG

Value rules and methods：Assuming that the probability that a certain geologic process of j-th candidates site occurs is P_i, i=1,2, 3 ..., K, then select these geologic process probability of happening the maximum to represent the probability that strong geologic process occurs in the candidate site, That is SPG_j=Max (P_i, i=1,2,3 ..., K)；It is required that SPG_jLess than a ten thousandth, work as SPG_jDuring less than a ten thousandth, j-th The SQV of candidate site₂ ¹Value take 10；The SQV of other sites₂ ¹Determine by equipartition principle and linear interpolation method.

(2) extreme meteorology and climate disaster probability of happening SPT

Value rules and methods：Assuming that the probability that a certain extreme meteorology and climate disaster of j-th candidates site occurs is P_i, I=1,2,3 ..., K, then select these meteorology and climate disaster probability of happening the maximum to represent the candidate site and extreme meteorology occur The probability of climate damage, i.e. SPT_j=Max (P_i, i=1,2,3 ..., K)；It is required that SPT_jThe probability of generation is less than 150 points One of, work as SPT_jWhen the probability of generation is less than 1/150th, the SQV of j-th candidates site₂ ²Value is 10, works as SPT_j When the probability of generation is 100 percent, the SQV of j-th candidates site₂ ²Value is 0, the SQV of other sites₂ ²By equipartition principle Determine with linear interpolation method.

(3) site to active fault apart from SMP

Value rules and methods：If SMP≤5km, SQV₂ ³=0；If 5km<SMP≤10km, SQV₂ ³=1；If 10km<SMP≤50km, SQV₂ ³=3；If 50km<SMP≤100km, SQV₂ ³=6；If SMP>100km, SQV₂ ³=10；

(4) place SL-2 grades of acceleration peak SSL

Value rules and methods：The SQV of SSL the maximum in candidate site₂ ⁴=1, SSL recklings in candidate site SQV₂ ⁴=10, the SQV of other candidate sites₂ ⁴Value presses equipartition principle and linear interpolation method determines；

(5) near field periphery in site points to site SZX without earthquake magnitude seismogenic fault band high

Value rules and methods：If being site near field in the range of the 25km of candidate site border periphery, if candidate site Site near field periphery there are the seismic structure bands that occur not less than 7 grades of earthquakes and the candidate site is pointed to, then from SZX Angle considers that the candidate site is not suitable for as site；Otherwise, consider that the candidate site is suitable as site from SZX angles； Location near field periphery points to site without earthquake magnitude seismogenic fault band high, then SQV₂ ⁵Value 10, site near field periphery has Earthquake magnitude seismogenic fault band high points to site, then value is 0.

(6) site SZJ is pointed to without seismogenic fault in the near field of site

Value rules and methods：If being site near field in the range of the 25km of candidate site border periphery, if candidate site Site near field in have seismogenic fault and point to candidate site, then consider that the candidate site is not suitable for as field from SZJ angles Location；Otherwise, consider that the candidate site is suitable as site from SZJ angles；Site is pointed to without seismogenic fault in the near field of site, SQV₂ ⁶Value is 10, has seismogenic fault to point to site, SQV in the near field of site₂ ⁶Value is 0.

III, hydrogeologic condition first class index SQV₃

(1) underground water hydraulic gradient SHI

Value rules and methods：If SHI<1 ‰, SQV₃ ¹=10；If SHI>100 ‰, SQV₃ ¹=0；If 1 ‰≤ SHI≤100 ‰, SQV₃ ¹According to equipartition principle and linear principle value；(2) rock mass ground infiltration coefficient accounting SK

Value rules and methods：Using mathematical statistics method, with Rockmass Permeabuity Coefficient≤1 × 10^-9The proportion of m/s is big Small determination SQV₃ ²Value；When infiltration coefficient SK≤1 × 10 of rock mass in certain site^-9Accounting P { SK │ SK≤1 × 10 of m/s^-9m/ S }=1 when, SQV₃ ²=10；When accounting meets 0.01<P{SK│SK≤1×10^-9M/s } ＜ 1 when, SQV₃ ²=10 × P；Work as accounting Meet P { SK │ SK≤1 × 10^-9M/s }≤0.01 when, SQV₃ ²=0.1；

(3) to region discharge area apart from SPD

Value rules and methods：If SPD<50km, SQV₃ ³=1；If SPD>200km, SQV₃ ³=110；If 50km ≤ SPD≤200km, SQV₃ ³According to equipartition principle and linear principle value；

(4) the coefficient of transmissibility STF of rift structure

Value rules and methods：In this stage, emphasis considers the transmissibility of the rift structure for cutting and passing through candidate site, Obtain quantitative data；Using mathematical statistics method, with coefficient of transmissibility≤1 × 10 of rift structure^-5m²The probability size of/s determines The suitability scale SQV of certain candidate site₃ ⁴Value；When STF≤1 × 10 in certain candidate site^-5m²/ s probability P STF │ STF≤ 1×10^-5m²/ s }=1 when, SQV₃ ⁴=10；When probability 0.01<P{STF│STF≤1×10^-5m²/ s } ＜ 1 when, SQV₃ ⁴=10 × P；When probability P { STF │ STF≤1 × 10^-5m²/ s }≤0.01 when, SQV₃ ⁴=0.1；

(5) the hydraulic connection index S HR in storehouse and important water body is disposed

Value rules and methods：The important water body in the range of the downstream 80km of candidate site is considered, if disposal storehouse candidate Ground water regime has hydraulic connection with downstream water where location, then site is not suitable for, SQV₃ ⁵Value is 0；Otherwise, site is suitable SQV₃ ⁵Value is 10；

IV, geochemical conditions first class index SQV₄

(1) adsorption coefficient SKD of the Characteristic radionuclides in rock forming mineral

Value rules and methods：The SQV of SKD values reckling in candidate site₄ ¹=1, SKD the maximum in candidate site SQV₄ ¹=10, the SQV of other candidate sites₄ ¹Determine by equipartition principle and linear interpolation method；

(2) adsorption coefficient SKF of the Characteristic radionuclides in fracture filling

Value rules and methods：The SQV of SKF values reckling in candidate site₄ ²=1, SKF the maximum in candidate site SQV₄ ²=10, the SQV of other candidate sites₄ ²Determine by equipartition principle and linear interpolation method；

(3) underground water TDS, i.e. STW

Value rules and methods：If STW>100g/l, SQV₄ ³=1；If STW<1g/l, SQV₄ ³=10；If 1g/l ≤ STW≤100g/l, SQV₄ ³Determined according to equipartition principle and linear interpolation method；(4) groundwater age SAW

Value rules and methods：If 30 years<SAW≤1000 year, SQV₄ ⁴=3；If 1000<SAW≤10000 year, SQV₄ ⁴=6；If SAW>, SQV in 10000₄ ⁴=10；

(5) sulfate concentration SSC in underground water

Value rules and methods：If SSC<500mg/L, SQV₄ ⁵=10；If SSC>5000mg/L, SQV₄ ⁵=1；Such as Fruit 500mg/L≤SSC≤5000mg/L, SQV₄ ⁵Determine according to equipartition principle and linear interpolation method；(6) chlorion is dense in underground water Degree SCC

Value rules and methods：If SCC<1000mg/l, SQV₄ ⁶=10；If SCC>5000mg/l, SQV₄ ⁶=1；Such as Fruit 1000mg/l≤SCC≤5000mg/l, SQV₄ ⁶Determine according to equipartition principle and linear interpolation method；(7) radioactive substance is in rock Effective diffusion cofficient SDE in stone

Value rules and methods：SDE maximum SQV in candidate site₄ ⁷=10, SDE minimum SQV in candidate site₄ ⁷= 1, the SQV of other candidate sites₄ ⁷Determine by equipartition principle and linear interpolation method；

(8) solubility SSD of the radioactive substance in ground water regime

Value rules and methods：SSD maximum SQV in candidate site₄ ⁸=1, SSD minimum SQV in candidate site₄ ⁸= 10, the SQV of other candidate sites₄ ⁸Determine by equipartition principle and linear interpolation method；

V, mankind's activity first class index

(1) site to available resources ground apart from SUD

Value rules and methods：The candidate site maximum for SUD, SQV₅ ¹=10, the candidate site minimum for SUD, SQV₅ ¹=1, the SQV of other candidate sites₅ ¹Determine by equipartition principle and linear interpolation method；

(2) probability SFH is flooded in the site that mankind's activity causes

Value rules and methods：The candidate site maximum for SFH, SQV₅ ²=1, the candidate site minimum for SFH, SQV₅ ²=10, the SQV of other candidate sites₅ ²Determine by equipartition principle and linear interpolation method；

VI, construction and engineering specifications first class index

(1) rock strength and stress ratio SRP

Value rules and methods：Assuming that rock uniaxiality strength R_CWith the maximum principal stress σ for measuring₁Ratio represents with RP, That is RP=R_C/σ₁；When the probability of RP >=5 in certain site is 1, i.e. SRP=P { RP │ RP_i>=5 }=1, SQV₆ ¹=10, wherein i= 1,2 ..., K, K are the experiment quantity completed in certain site；As 0.1≤SRP ＜ 1, SQV₆ ¹=10 × P；As SRP ＜ 0.1, SQV₆ ¹=1；

(2) rock mass integrity index index S KV

Value rules and methods：Rock mass completeness be divided into it is complete, more complete, relatively broken, broken, smashed to pieces five etc. Level, its corresponding rock mass integrity index Kv scope is respectively：Kv>0.75,0.55<Kv≤0.75,0.35<Kv≤0.55, 0.15≤Kv≤0.35, Kv<0.15, its corresponding suitability scale SQV₆ ²Value is respectively 10,4,1,0.1,0；

(3) under temperature and effect of irradiation rock damaging stress value STR

Value rules and methods：The candidate site maximum for STR, SQV₆ ³=10, the candidate site minimum for STR, SQV₆ ³=1, the SQV of other candidate sites₆ ³Determine by equipartition principle and linear interpolation method；

(4) under temperature and effect of irradiation rock creep rate SRR

Value rules and methods：The candidate site maximum for SRR, SQV₆ ⁴=1, the candidate site minimum for SRR, SQV₆ ⁴=10, the SQV of other candidate sites₆ ⁴Determine by equipartition principle and linear interpolation method；

(5) place is than SBG high

Value rules and methods：The candidate site maximum for SBG, SQV₆ ⁵=1, the candidate site minimum for SBG, SQV₆ ⁵=10, the SQV of other candidate sites₆ ⁵Determine by equipartition principle and linear interpolation method；

VII, waste traffic condition first class index SQV₇

(1) to spentnuclear fuel commercialization reprocessing plant apart from SDM

Value rules and methods：The candidate site maximum for SDM, SQV₇ ¹=1, the candidate site minimum for SDM, SQV₇ ¹=10, the SQV of other candidate sites₇ ¹Determine by equipartition principle and linear interpolation method；

(2) the distance SGD of reconstruction transport routes is increased newly

Value rules and methods：The candidate site maximum for SGD, SQV₇ ²=1, the candidate site minimum for SGD, SQV₇ ²=10, the SQV of other candidate sites₇ ²Determine by equipartition principle and linear interpolation method；

(3) waste freight STD

Value rules and methods：The candidate site maximum for STD, SQV₇ ³=1, the candidate site minimum for STD, SQV₇ ³=10, the SQV of other candidate sites₇ ³Determine by equipartition principle and linear interpolation method；

(4) the exposure dose of radiation SPR produced in transporting

Value rules and methods：The candidate site maximum for SPR, SQV₇ ⁴=1, the candidate site minimum for SPR, SQV₇ ⁴=10, the SQV of other candidate sites₇ ⁴Determine by equipartition principle and linear interpolation method；

VIII, environmental protection condition first class index SQV₈

(1) to the Weighted distance index S DS of national nature reserve

Value rules and methods：Scope of statistics is 80km around candidate site, if 80km scopes around i-th candidate site Inside there are national nature reserve m, its distance DS_jRepresent；According to the weight for from the close-by examples to those far off giving apart from site each protection zone, Use k_jRepresent, apart from DS_jWith weight k_jRelation be shown below：

If apart from k_jInterior protection zone quantity Y_kRepresent, wherein k is the interval quantity of distance, then i-th candidate site is to state The Weighted distance index S DS of family's level protection zone_iIt is expressed as follows：

During SDS=0, SQV₈ ¹=0；During SDS ≠ 0, the candidate site maximum for SDS, SQV₈ ¹=10, it is minimum for SDS Candidate site, SQV₈ ¹=1, other situations press equipartition principle and linear interpolation method determines；

(2) to the Weighted distance index S DA of provincial protection zone

Value rules and methods：Scope of statistics is 80km around candidate site, if 80km scopes around i-th candidate site Inside there are provincial protection zone m, its distance DA_jRepresent；According to the weight for from the close-by examples to those far off giving apart from site each protection zone, use k_jRepresent, distance is shown below with the relation of weight：

If apart from k_jInterior protection zone quantity Y_kRepresent, wherein k is the interval quantity of distance, then the candidate site is to provincial guarantor Protect the Weighted distance index S DA in area_iIt is expressed as follows：

During SDA=0, SQV₈ ²=0；During SDA ≠ 0, the candidate site maximum for SDA, SQV₈ ²=10, it is minimum for SDA Candidate site, SQV₈ ²=1, other situations press equipartition principle and linear interpolation method determines；(3) to the weighting of the important water body in downstream Range index SDW

Value rules and methods：Scope of statistics is 80km around candidate site, it is assumed that i-th candidate site downstream peripheral There are important water body m, its distance DW in the range of 80km_jRepresent；According to apart from site downstream is from the close-by examples to those far off given, each is important The weight of water body, uses k_jRepresent, distance is shown below with the relation of weight：

If apart from k_jInterior protection zone quantity Y_kRepresent, wherein k is the interval quantity of distance, then the candidate site is to arriving downstream The Weighted distance index S DW of important water body_iIt is expressed as follows：

During SDW=0, SQV₈ ³=0；During SDW ≠ 0, the candidate site maximum for SDW, SQV₈ ³=10；It is minimum for SDW Candidate site, SQV₈ ³=1, other situations press equipartition principle and linear interpolation method determines；

IX, Land_use change condition first class index SQV₉

(1) land use pattern SLC

Value rules and methods：If land type is arable land, field, forest land, the meadow in farming land, SQV₉ ¹=5；Such as Fruit land use pattern is that the business in construction land takes ground, industrial and mineral warehouse land, residential land, public administration and public service Land used, special use area, SQV₉ ¹=0；If land type is other soils in unused land, SQV₉ ¹=10；If soil is sharp It is the type in addition to the above-mentioned type enumerated, SQV with type₉ ¹=7；

(2) Land_use change value SVL

Value rules and methods：It is high, medium and low and without four grades of value that Land_use change value is divided into value, Its suitability scale SQV₉ ²Value is followed successively by 1,3,6,10；

X, social economy and humane condition first class index SQV₁₀

(1) population of ethnic minorities ratio SMP

Value rules and methods：Assuming that being to evaluate area, SMP parameter values and suitability with site surrounding area radius 80km Scale value SQV₁₀ ¹Between relation it is as follows:

(2) whether local government supports SCS

Value rules and methods：Investigation statisticses scope is counties and cities' level for being related in the range of 80km around site and with Shangdi Fang Zhengfu；If local government actively supports, SQV₁₀ ²=10；If local government does not oppose, SQV₁₀ ²=6；If local political affairs Mansion opposition, SQV₁₀ ²=0；

(3) the local public supports the ratio SPS of site

Value rules and methods：Permanent resident population to candidate site region carries out survey, and scope of statistics is field The counties and cities level administrative region being related in the range of 80km around location；If the local public supports the ratio SPS of site<50%, SQV₁₀ ³=0；If the local public supports ratio SPS >=50% of site, SQV₁₀ ³Equipartition principle and line are pressed according to different SPS values Property interpolation method determine.

Beneficial effects of the present invention are：

(1) the Hlw Geological Repository site suitability evaluation methods based on Delphi method of the invention, to set up Based on evaluation index, agriculture products weight is carried out by means of Delphi method, be with the Codes and Standards for building indices quantification Core, solves the problem of Hlw Geological Repository site suitability quantitative assessment, it is adaptable to high-level waste geology treatment Storehouse etc. is large-scale, complicated and site suitability quantitative evaluation of safety evaluation cycle overlength (more than 10,000 years)；

(2) the Hlw Geological Repository site suitability evaluation methods based on Delphi method of the invention are used, it is complete Into in high-level waste geology treatment Xinjiang pre-selected zone Jamin Soviet Union and Tian Hu locations candidate's sieving disposal site and evaluation study project The suitability evaluation of key sections and candidate site are elected, and are full marks, Jamin Soviet Union location suitability scale with 10 points by evaluating It is that 8.26, Tian Hu locations suitability scale is 8.08, therefore, it is determined that the suitability in Jamin Soviet Union location is better than Tian Hu locations, main table Present Jamin Soviet Union location has that rock mass area is larger, complete, permeability is relatively low, Fracture Networks agensis, rift structure or alteration The advantages of band water transmitting ability is extremely weak.

Specific embodiment

The Hlw Geological Repository site suitability based on Delphi method of the invention is commented with reference to embodiment Valency method is described in detail.

Hlw Geological Repository site suitability evaluation methods based on Delphi method of the invention, including following 4 Individual step：

Step 1. determines the index of Hlw Geological Repository site suitability quantitative assessment

Determine the basis of index this method of Hlw Geological Repository site suitability quantitative assessment.Present invention reference 《World's core geological sciences》Journal article is " at the beginning of the index system of Hlw Geological Repository site suitability quantitative evaluation The principle determined in spy " and the index system for proposing, determine the finger of Hlw Geological Repository site suitability quantitative assessment Mark.

Index is divided into Site Selection stage, Regional survey stage, site evaluation and the stage of recognition three phases altogether, by 10 First class index and 59 two-level index are combined.Wherein, the site suitability quantitative evaluation index in Site Selection stage by Geological conditions, following natural trend condition, geochemical conditions, environmental protection condition and social economy and humane condition etc. 5 Individual first class index and 15 two-level index are constituted；The site suitability quantitative evaluation index in Regional survey stage is by geology bar 5 first class index such as part, hydrogeologic condition, geochemical conditions, environmental protection condition and social economy and humane condition Constituted with 19 two-level index；Site evaluate and the stage of recognition site suitability quantitative evaluation index by geological conditions, not Coming natural trend, hydrogeology, geochemistry, construction and engineering specifications, mankind's activity, waste transport, environmental protection, soil makes Constituted with 10 first class index such as, social economy and humane condition and 43 two-level index.The first class index in each stage and two grades Indicator combination relation is as shown in table 1.

The Hlw Geological Repository site suitability quantitative assessing index list of table 1

Step 2. builds the Codes and Standards of indices quantification

The Codes and Standards for building indices quantification are the cores of this method, it is therefore an objective to setting up index and site suitability Quantitative relationship, i.e., establish each index parameter for Hlw Geological Repository site suitability quantitative assessing index With the functional relation between suitability scale.

So-called index parameter, refers to data, the information of the description disposal a certain feature in storehouse site or attribute.So-called suitability Scale, refers to the quantification numerical value for measuring first class index or two-level index to site suitability influence degree.

Represent Site Selection stage, Regional survey stage, site evaluation and the stage of recognition respectively using PQV, QQV and SQV Index suitability scale.Specifically, PQV_i1 ^j1Represent 1 two-level index of jth in the i-th 1 first class index of Site Selection stage Suitability scale；QQV_i2 ^j2Represent the suitability mark of 2 two-level index of jth in the i-th 2 first class index of Regional survey stage Degree；SQV_i3 ^j3Represent site evaluate and the i-th 3 first class index of the stage of recognition in 3 two-level index of jth suitability scale.

Step 2.1. Site Selection stages each index suitability scale value

In this stage, site region screening index system is by geological conditions, following natural trend condition, geochemistry bar 5 first class index such as part, environmental protection condition, social economy and humane condition and 15 two-level index are constituted, as shown in table 1. Assuming that i-th suitability scale of first class index is PQV_i, wherein i=1~5.

1st, geological conditions first class index PQV₁

(1) candidate's country rock thickness PRT

Value rules and methods：Candidate's country rock thickness is the bigger the better.If argillaceous rocks and the continuously distributed thickness of rock salt geologic body Degree is less than 100m, and other rock class thickness are less than 500m, then suitability scale PQV₁ ¹=0；If argillaceous rocks and rock salt geologic body connect Continuous spreading depth is less than 100m, and other rock class thickness are less than 500m, suitability scale PQV₁ ¹=10.

(2) candidate's country rock buried depth PRD

Value rules and methods：It is general to require that disposal storehouse is imbedded in underground 500m~1000m depth boundses, the depth mark Standard is international consensus, if meeting the depth bounds, PQV₁ ²=10；Otherwise, PQV₁ ²=0.

(3) candidate's country rock formation dip PFA

Value rules and methods：For non-laminar class candidate's Rock Mass, PQV₁ ³=10；For stratiform class geologic body, when PFA<At 5 °, PQV₁ ³=10, as 5 °≤PFA<At 10 °, PQV₁ ³=5, other situations PQV₁ ³=0.

2nd, following natural trend condition first class index PQV₂

(1) the dynamic peak accelerator PSL of regional earthquake

Value rules and methods：With the earthquake that dynamic parameter zoning map (GB18306-2001) of national standard Earthquake In China is given Dynamic peak accelerator partition value is foundation, and the region in low value is more preferable than being in the region suitability of high level.If PSL> 0.2g, PQV₂ ¹=0；If 0.15g<PSL≤0.2g, PQV₂ ¹=1；If 0.1g<PSL≤0.15g, PQV₂ ¹=4；If 0.05g<PSL≤0.1g, PQV₂ ¹=7；If PSL≤0.05g, PQV₂ ¹=10.

(2) to Regional Geotectonic secondary unit border apart from PGZ

Value rules and methods：The regional geotectonic classification write with Zhang Peizhen academician and classification defined two Level structural unit be foundation, if pre-selected zone border to Regional Geotectonic secondary unit border apart from PGZ<25km, then PQV₂ ²=0, i.e., should reduce the scope of pre-selected zone；If PGZ>100km, then PQV₂ ²=10；If 25km≤PGZ≤ 100km, PQV₂ ²According to equipartition principle and linear interpolation method determination, PGZ index parameters value and suitability scale value PQV₂ ²Between Relation is as shown in table 2.

Table 2PGZ parameter values and suitability scale PQV₂ ²Value relation table

PGZ(km)	<25	25~33	33~42	42~50	50~58	58~67	67~75	75~83	83~92	92~100	>100
												PQV₂ ²	0	1	2	3	4	5	6	7	8	9	10

(3) to the strong zone of influence of extreme geologic process apart from PGA

Value rules and methods：If PGA<25km, then PQV₂ ³=0；If PGA>100km, then PQV₂ ³=10；If 25km≤PGA≤100km, PQV₂ ³Equipartition principle value, PGA index parameters value and suitability scale value are pressed according to distance PQV₂ ³Between relation it is as shown in table 3.

Table 3PGA parameter values and suitability scale PQV₂ ³Value relation table

(4) natural geologic hazard development degree PND

Value rules and methods：This report with《Western part of China geological disaster figure (1:250 ten thousand)》The standard and subregion of definition It is foundation, natural geologic hazard development degree can be divided into strong, medium, weak and four grades of agensis, this four grade correspondences Suitability scale PQV₂ ⁴It is followed successively by 0,2,7,10.

3rd, geochemical conditions first class index PQV₃

(1) the main surrounding rock type PTR in pre-selected zone

Value rules and methods：When PTR is four rocks such as CRYSTALLINE ROCKS, argillaceous rocks, tufa stone and rock salt, PQV₃ ¹= 10；Otherwise, PQV₃ ¹=0.

4th, environmental protection condition first class index PQV₄

(1) to Head streams protection zone apart from PSS

Value rules and methods：If PSS=0, PQV₄ ¹=0；If PSS>80km, then PQV₄ ¹=10；If 0< PSS≤80km, according to equipartition principle and linear interpolation principle value, PSS parameter value and suitability scale PQV₄ ¹Between relation As shown in table 4.

Table 4PSS parameter values and suitability scale PQV₄ ¹Value relation table

PSS(km)	0	0~9	9~18	18~27	27~36	36~45	45~54	54~63	63~72	72~80	>80
												PQV₄ ¹	0	1	2	3	4	5	6	7	8	9	10

(2) the quantity PEP of environmental sensitive area above the provincial level

Value rules and methods：Principal statistical is to preselect district center as laws and regulations regulation in round dot 80km radius The quantity of provincial and national environmental sensitive area.If PEP<20, then PQV₄ ²=10；If PEP>500, then PQV₄ ²= 1；If 20≤PEP≤500, determined by equipartition principle and linear interpolation method.PEP parameter values and PQV₄ ²Between relation As shown in table 5.

Table 5PEP parameter values and suitability scale PQV₄ ²Value relation table

(3) the quantity PWP of important water body

Value rules and methods：Principal statistical with preselect district center as in round dot 80km radius in downstream direction law The quantity of the provincial and national important water body of rules and regulations.If PWP<10, then PQV₄ ³=10；If PWP>100, then PQV₄ ³=1；If 10≤PWP≤100, determined by equipartition principle and linear interpolation method.PWP parameter values and suitability mark Relation between angle value is as shown in table 6.

Table 6PWP parameter values and suitability scale PQV₄ ³Value relation table

(4) pre-selected zone absolute height above sea level PHE

Value rules and methods：If the height above sea level of pre-selected zone highest point is more than 3500m, PQV₄ ⁴=0；Otherwise, PQV₄ ⁴ =10.

5th, social economy and humane condition first class index PQV₅

(1) the absolute density of population PAP in region

Value rules and methods：The absolute density of population in region refers to and resided in pre-selected zone population and floating population's summation and preselect The area ratio in area.If PAP >=400 people/kilometer², PQV₅ ¹=0；If 200 people/kilometer²≤PAP<400 people/kilometer², PQV₅ ¹=1；If 100 people/kilometer²≤PAP<200 people/kilometer², PQV₅ ¹=3；If 20 people/kilometer²≤PAP<100 people/public affairs In², PQV₅ ¹=6；If PAP<20 people/kilometer², PQV₅ ¹=10.

(2) Main Towns quantity PTW

Value rules and methods：Main Towns quantity refers to medium and above scale cities and towns the quantity included in pre-selected zone. If PTW≤5, PQV₅ ²=10；If 5<PTW≤50, PQV₅ ²=7；If 50<PTW≤100, PQV₅ ²=3； If PTW>100, PQV₅ ²=1.

(3) region total output value PDP per capita

Value rules and methods：Provincial administrative area gives birth to per capita where this phase zone per capita total output value pre-selected zone Total value is produced to characterize.If pre-selected zone is averagely added up afterwards across multiple provincial administrative areas according to area.Assuming that the whole nation produces always per capita Value PDP_nRepresent, Provincial administrative area per capita total output value come top 10 total output value per capita average value PDP₁₀ Represent.If PDP≤PDP_n, PQV₅ ³=10；If PDP >=PDP₁₀, PQV₅ ³=1；In the case of other, PQV₅ ³=5.

Step 2.2. Regional survey stages each index suitability scale value

In this stage, site region screening index system is by geological conditions, hydrogeologic condition, geochemical conditions, ring 5 first class index such as border protective condition, social economy and humane condition and 19 two-level index are constituted, as shown in table 1.Assuming that the The i suitability scale QQV of first class index_iRepresent, wherein i=1~5.

1st, geological conditions first class index QQV₁

(1) candidate's country rock area QGA

Value rules and methods：If QGA<50km², QQV₁ ¹=1；If QGA>300km², then QQV₁ ¹=10；If 50km²≤QGA≤300km², QQV₁ ¹According to equipartition principle and linear interpolation principle value.Candidate's country rock area QGA parameter values with Relation between suitability scale value is as shown in table 7.

Table 7QGA parameter values and suitability scale QQV₁ ¹Value relation table

(2) candidate's country rock thickness QRT

Value rules and methods：If QRT<100m, QQV₁ ²=0；If QRT>500m, QQV₁ ²=10；If 100m≤ QRT≤500m, QQV₁ ²According to equipartition principle and linear principle value.Relation between QRT parameter values and suitability scale value is such as Shown in table 8.

Table 8QRT parameter values and suitability scale QQV₁ ²Value relation table

(3) candidate's country rock buried depth QRD

Value rules and methods：Value rule with the Site Selection stage is identical, general to require that disposal storehouse is imbedded in underground In 500m~1000m depth boundses, if meeting the depth bounds, QQV₁ ³=10；Otherwise, QQV₁ ³=0.

(4) rift structure distribution surface density QSD

Value rules and methods：Rift structure refers to the rift structure for cutting and passing through pre-selection location, rift structure distribution face Density refers to the rift structure quantity in unit area under the conditions of uniform or more uniform distribution.If QSD<1/km², QQV₁ ⁴= 10；If QSD>10/km², QQV₁ ⁴=1；If 1/km²≤ QSD≤100/km², QQV₁ ⁴According to equipartition principle and line Property principle value.Relation between QSD parameter values and suitability scale value is as shown in table 9.

Table 9QSD parameter values and suitability scale QQV₁ ⁴Value relation table

QSD (bars/km²)	<1	1~2	2~3	3~4	4~5	5~6	6~7	7~8	8~10	>10
											QQV₁ ⁴	10	9	8	7	6	5	4	3	2	1

2nd, hydrogeologic condition first class index QQV₂

(1) groundwater level depth QWB

Value rules and methods：Assuming that the buried depth in proposed disposal storehouse is MZ, if QWB≤MZ, QQV₂ ¹=5；Otherwise, QQV₂ ¹=10.

(2) underground water hydraulic gradient QHI

Value rules and methods：If QHI<1 ‰, QQV₂ ²=10；If QHI>100 ‰, QQV₂ ²=0；If 1 ‰≤ QHI≤100 ‰, QQV₂ ²According to equipartition principle and linear principle value.Relation between QHI parameter values and suitability scale value is such as Shown in table 10.

Table 10QHI parameter values and suitability scale QQV₂ ²Value relation table

QHI(‰)	<1	1~12	12~23	23~34	34~45	45~56	56~67	67~78	78~89	89~100	>100
												QQV₂ ²	10	9	8	7	6	5	4	3	2	1	0

(3) rock mass ground infiltration coefficient accounting QK

Value rules and methods：Using mathematical statistics method, with rock mass hyposmosis coefficient≤1 × 10^-8M/s proportions are big Small determination QQV₂ ³Value.When hyposmosis coefficient≤1 × 10 of rock mass in certain location^-8Accounting P { QK │ QK≤1 × 10 of m/s^-8m/ S }=1 when, QQV₂ ³=10；When 0.01<P{QK│QK≤1×10^-8M/s } ＜ 1 when, QQV₂ ³=10 × P；When P QK │ QK≤1 × 10^-8M/s }≤0.01 when, QQV₂ ³=0.1.

(4) to region discharge area apart from QPD

Value rules and methods：If QPD<50km, QQV₂ ⁴=1；If QPD>200km, QQV₂ ⁴=10；If 50km ≤ QPD≤200km, QQV₂ ⁴By equipartition principle and linear principle value.To region discharge area apart from QPD index parameters value with it is suitable Relation between suitable property scale value is as shown in table 11.

Table 11QPD parameter values and suitability scale QQV₂ ⁴Value relation table

(5) the transmissibility QTF of rift structure

Value rules and methods：Important Faults construction refers to the rift structure for cutting and passing through rock mass.In this stage, general root According to the property and charges feature qualitative discrimination of rift structure, be divided into water transmitting ability it is poor, in, strong three kinds of situations, its is corresponding suitable Suitable property scale QQV₂ ⁵Value is followed successively by 10,1,0.1.

3rd, geochemical conditions first class index QQV₃

(1) underground water TDS, i.e. QTW

Value rules and methods：If QTW>100g/l, QQV₃ ¹=1；If QTW<1g/l, QQV₃ ¹=10；If 1g/l ≤ QTW≤100g/l, QQV₃ ¹Determined according to equipartition principle and linear interpolation method.Between QTW parameter values and suitability scale value Relation is as shown in table 12.

Table 12QTW parameter values and suitability scale QQV₃ ¹Value relation table

QTW(g/l)	<1	1~13	13~25	25~37	37~49	49~61	61~73	73~85	85~100	>100
											QQV₃ ¹	10	9	8	7	6	5	4	3	2	1

(2) underground water Eh values, i.e. QEH

Value rules and methods：If QEH=0, QQV₃ ²=6；If QEH<0, QQV₃ ²=10；If QEH>0, QQV₃ ² =1.

(3) groundwater age QAW

Value rules and methods：If QAW≤30 year, QQV₃ ³=0；If 30 years<QAW≤1000 year, QQV₃ ³=3；Such as Really 1000<QAW≤10000 year, QQV₃ ³=6；If QAW>, QQV in 10000₃ ³=10.

(4) sulfate concentration QSC in underground water

Value rules and methods：If QSC<500mg/L, QQV₃ ⁴=10；If QSC>5000mg/L, QQV₃ ⁴=1；Such as Fruit 500mg/L≤QSC≤5000mg/L, QQV₃ ⁴Determine according to equisection method and linear interpolation method.QSC parameter values and suitability scale Relation between value is as shown in table 13.

Table 13QSC parameter values and suitability scale QQV₃ ⁴Value relation table

(5) chlorine ion concentration QCC in underground water

Value rules and methods：If QCC<1000mg/l, QQV₃ ⁵=10；If QCC>5000mg/l, QQV₃ ⁵=1；Such as Fruit 1000mg/l≤QCC≤5000mg/l, QQV₃ ⁵Determine according to equisection method and linear interpolation method.QCC parameter values and suitability mark Relation between angle value is as shown in table 14.

Table 14QCC parameter values and suitability scale QQV₃ ⁵Value relation table

4th, environmental protection condition first class index QQV₄

(1) the quantity QEP of environmental sensitive area above the provincial level

Value rules and methods：The province that principal statistical is specified with Section Center as laws and regulations in round dot 25km radius Level and the quantity of national environmental sensitive area.If QEP<5, QQV₄ ¹=10；If QEP>50, QQV₄ ¹=1；5≤ QEP≤50, QQV₄ ¹Determine by equipartition principle and linear interpolation method.Relation between QEP parameter values and suitability scale value As shown in Table 15.

Table 15QEP parameter values and suitability scale QQV₄ ¹Value relation table

QEP (individual)	<5	5~11	11~17	17~23	23~29	29~35	35~40	40~45	45~50	>50
											QQV₄ ¹	10	9	8	7	6	5	4	3	2	1

(2) the quantity QWP of important water body

Value rules and methods：Principal statistical is specified with preselecting Section Center as laws and regulations in round dot 25km radius Provincial and national important water body quantity.If QWP<5, QQV₄ ²=10；If QWP>20, QQV₄ ²=1；If 5 Individual≤QWP≤20, QQV₄ ²Determine by equipartition principle and linear interpolation method.Between QWP parameter values and suitability scale value Relation is as shown in table 16.

Table 16QWP parameter values and suitability scale QQV42 value relation tables

QWP (individual)	<5	5~7	7~9	9~11	11~13	13~15	15~17	17~19	19~20	>20
											QQV₄ ²	10	9	8	7	6	5	4	3	2	1

5th, social economy and humane condition first class index QQV₅

(1) the absolute density of population QAP in region

Value rules and methods：The absolute density of population QAP in region refers in the range of with Section Center as center of circle radius 15km Resident population and floating population's summation and area ratio, QAP is more big more is not suitable for addressing.If QAP<5 people/km², QQV₅ ¹= 10；If QAP>200 people/km², QQV₅ ¹=0；If 5 people/km²≤ QAP≤200 people/km², QQV₅ ¹By equipartition principle and linearly Interpolation method determines.Relation between QAP parameter values and suitability scale value is as shown in table 17.

Table 17QWP parameter values and suitability scale QQV₅ ¹Value relation table

(2) Main Towns quantity QTW

Value rules and methods：Main Towns quantity QTW refers to medium in the range of with Section Center as center of circle radius 15km And the quantity in above scale cities and towns.If QTW>50, QQV₅ ²=1；If 10<QTW≤50, QQV₅ ²=3；If 3<QTW≤ 10, then QQV₅ ²=6；If QTW≤3, QQV₅ ²=10.

(3) total output value QDP per capita

Value rules and methods：This stage, total output value QDP referred to giving birth to per capita for local administrative area where pre-selection location per capita Produce total value.If pre-selection location is averagely added up afterwards across multiple local administrative areas according to area.Assuming that national local administrative area The average value of total output value is QDP per capita_x；Total output value comes giving birth to per capita for inverse the 50th per capita in national local administrative area Produce total value QDP_-50Represent.If QDP>2QDP_x, QQV₅ ³=1；If QDP_x<QDP≤2QDP_x, QQV₅ ³=3；If QDP_-50 <QDP≤QDP_x, QQV₅ ³=6；If QDP≤QDP_-50, QQV₅ ³=10.

Site is evaluated and the stage of recognition site suitability quantitative evaluation index system is by geological conditions, following natural trend Condition, hydrogeologic condition, geochemical conditions, mankind's activity, construction and engineering specifications, waste traffic condition, environmental protection 10 first class index such as condition, Land_use change condition, social economy and humane condition and 43 two-level index are constituted, such as the institute of table 1 Show.Assuming that i-th score value SQV of first class index_iRepresent, wherein i=1~10.

1st, geological conditions first class index SQV₁

(1) candidate's country rock area SGA

Value rules and methods：If candidate's country rock area SGA<10km², SQV₁ ¹=1；If 10km²≤SGA<20km², SQV₁ ¹=3；If 20km²≤SGA<50km², SQV₁ ¹=6；If SGA >=50km², SQV₁ ¹=10.

(2) candidate's country rock thickness SRT

Value rules and methods：If candidate's country rock thickness SRT<100m, SQV₁ ²=0；If 100m≤SRT<200m, SQV₁ ²=3；If 200m≤SRT<400m, SQV₁ ²=6；If SRT >=400m, SQV₁ ²=10.

(3) rock thermal conductivity factor SCT

Value rules and methods：The SQV of SCT the maximum in candidate site₁ ³=10, SCT recklings in candidate site SQV₁ ³=1, the SQV of other candidate sites₁ ³Determine by equipartition principle and linear interpolation method.

(4) rift structure distribution surface density SSD

Value rules and methods：This stage only counts the rift structure for cutting and passing through candidate site.If fracture structure Make distribution surface density SSD>5/km², SQV₁ ⁴=1；If 5/km²≥SSD>1/km², SQV₁ ⁴=3；If 1/km²≥SSD>0.2/ km², SQV₁ ⁴=6；If SSD≤0.2/km², SQV₁ ⁴=10.

(5) to the Weighted distance SFD of Important Faults construction

Value rules and methods：In this stage, scope of statistics is 25km scopes around candidate site.So-called Important Faults structure Make the rift structure for referring to and cutting and passed through pre-selection location.Assuming that having Important Faults in the range of 25km around i-th candidate site Construction m, its distance DS_jRepresent；The weight that every Important Faults are constructed from the close-by examples to those far off is given according to apart from site, k is used_jTable Show, distance is as shown in table 18 with the relation of weight；Assuming that Important Faults construction quantity Y in the corresponding distance of table 18 is interval_kTable Show, wherein k is the interval quantity of distance.The Weighted distance index S FD that then i-th candidate site constructs to Important Faults_iRepresent It is as follows：

Distance and weight mapping table that the site of table 18 constructs to Important Faults

For all candidate sites, the SQV of SFD the maximum₁ ⁵The SQV of=10, SFD reckling₁ ⁵=1, other sites SQV₁ ⁵Determine by equipartition principle and linear interpolation method.

2nd, following natural trend condition first class index SQV₂

(1) strong geologic process probability of happening SPG

Value rules and methods：Assuming that the probability that a certain geologic process of j-th candidates site occurs is P_i(i=1,2, 3 ..., K), then select these geologic process probability of happening the maximum to represent the probability that strong geologic process occurs in the candidate site, That is SPG_j=Max (P_i, i=1,2,3 ..., K).According to the relevant administrative provisions of China, the disposal storehouse safety evaluation phase is at least 10,000 Year, therefore the general maximum probability for requiring the strong geologic process of generation is to meet for 10000 years one, i.e. SPG_jShould be less than a ten thousandth.When SPG_jLess than a ten thousandth, the SQV of j-th candidates site₂ ¹Value is 10, the SQV of other sites₂ ¹It is by equipartition principle and linearly slotting The determination of value method,

(2) extreme meteorology and climate disaster probability of happening SPT

Value rules and methods：Assuming that the probability that a certain extreme meteorology and climate disaster of j-th candidates site occurs is P_i (i=1,2,3 ..., K), then select these meteorology and climate disaster probability of happening the maximum to represent the candidate site and extreme gas occur As the probability of climate damage, i.e. SPT_j=Max (P_i, i=1,2,3 ..., K).According to domestic and foreign experience, the storehouse safe operation phase is disposed It is 150 years, therefore typically requires SPT_jThe probability of generation should be less than 1/150th.Work as SPT_jThe probability of generation is less than one When 1/150th, the SQV of j-th candidates site₂ ²Value is 10, works as SPT_jWhen the probability of generation is 100 percent, jth The SQV of individual candidate site₂ ²Value is 0, the SQV of other sites₂ ²Determine by equipartition principle and linear interpolation method.

(3) site to active fault apart from SMP

Value rules and methods：If during SMP≤5km, SQV₂ ³=0；If 5km<SMP≤10km, SQV₂ ³=1；If 10km<SMP≤50km, SQV₂ ³=3；If 50km<SMP≤100km, SQV₂ ³=6；If SMP>100km, then SQV₂ ³=10.

(4) place SL-2 grades of acceleration peak SSL

Value rules and methods：According to relevant standard, by the place of field investigation candidate site different with acquisition is calculated SL-2 grades of acceleration peak.The SQV of SSL the maximum in candidate site₂ ⁴=1；The SQV of SSL recklings in candidate site₂ ⁴ =10, the SQV of other candidate sites₂ ⁴Value presses equipartition principle and linear interpolation method determines.

Value rules and methods：Assuming that be site near field in the range of the 25km of candidate site border periphery, if candidate The site near field periphery of location has the seismic structure bands that occur not less than 7 grades of earthquakes and points to the candidate site, then from SZX angles consider that the candidate site is not suitable for as site；Otherwise, consider that the candidate site is suitable as site from SZX angles. Site near field periphery points to site without earthquake magnitude seismogenic fault band high, then SQV₂ ⁵Value is 10, and site near field periphery is attached Closely there is earthquake magnitude seismogenic fault band high to point to site, then value is 0.

(6) site SZJ is pointed to without seismogenic fault in the near field of site

Value rules and methods：Assuming that be site near field in the range of the 25km of candidate site border periphery, if candidate There is seismogenic fault in the site near field of location and point to candidate site, then consider that the candidate site is not suitable for conduct from SZJ angles Site；Otherwise, consider that the candidate site is suitable as site from SZJ angles.Site is pointed to without seismogenic fault in the near field of site, SQV₂ ⁶Value is 10, has seismogenic fault to point to site, SQV in the near field of site₂ ⁶Value is 0；The SQV of other candidate sites₂ ⁶Value Determine by equipartition principle and linear interpolation method.

3rd, hydrogeologic condition first class index SQV₃

(1) underground water hydraulic gradient SHI

Value rules and methods：If SHI<1 ‰, SQV₃ ¹=10；If SHI>100 ‰, SQV₃ ¹=0；If 1 ‰≤ SHI≤100 ‰, SQV₃ ¹According to equipartition principle and linear principle value.Relation between SHI parameter values and suitability scale value is such as Shown in table 19.

Table 19SHI parameter values and suitability scale SQV₃ ¹Value relation table

SHI(‰)	<1	1~12	12~23	23~34	34~45	45~56	56~67	67~78	78~89	89~100	>100
												SQV₃ ¹	10	9	8	7	6	5	4	3	2	1	0

(2) rock mass ground infiltration coefficient accounting SK

Value rules and methods：Using mathematical statistics method, with Rockmass Permeabuity Coefficient≤1 × 10^-9The proportion of m/s is big Small determination SQV₃ ²Value.When infiltration coefficient SK≤1 × 10 of rock mass in certain site^-9Accounting P { SK │ SK≤1 × 10 of m/s^-9m/ S }=1 when, SQV₃ ²=10；When accounting meets 0.01<P{SK│SK≤1×10^-9M/s } ＜ 1 when, SQV₃ ²=10 × P；Work as accounting Meet P { SK │ SK≤1 × 10^-9M/s }≤0.01 when, SQV₃ ²=0.1.

(3) to region discharge area apart from SPD

Value rules and methods：If SPD<50km, SQV₃ ³=1；If SPD>200km, SQV₃ ³=110；If 50km ≤ SPD≤200km, SQV₃ ³According to equipartition principle and linear principle value.Relation such as table 20 between SPD and suitability scale value It is shown.

Table 20SPD parameter values and suitability scale SQV₃ ³Value relation table

(4) the coefficient of transmissibility STF of rift structure

Value rules and methods：In this stage, emphasis considers the transmissibility of the rift structure for cutting and passing through candidate site, Quantitative data is obtained by methods such as hydrogeological tests.Using mathematical statistics method, with coefficient of transmissibility≤1 of rift structure × 10^-5m²The probability size of/s determines the suitability scale SQV of certain candidate site₃ ⁴Value.When STF≤1 × 10 in certain candidate site^-5m²Probability P { STF │ STF≤1 × 10 of/s^-5m²/ s }=1 when, SQV₃ ⁴=10；When probability 0.01<P{STF│STF≤1×10^- ⁵m²/ s } ＜ 1 when, SQV₃ ⁴=10 × P；When probability P { STF │ STF≤1 × 10^-5m²/ s }≤0.01 when, SQV₃ ⁴=0.1.

Value rules and methods：It is main to consider in the range of the downstream 80km of candidate site, if having especially important beyond 80km Water body, also need investigation.Site is not if ground water regime and downstream water where disposal storehouse candidate site have hydraulic connection Suitably, SQV₃ ⁵Value is 0；Otherwise, the suitable SQV in site₃ ⁵Value is 10, the SQV of other candidate sites₃ ⁵By equipartition principle and line Property interpolation method determine.

4th, geochemical conditions first class index SQV₄

Value rules and methods：Determined according to different SKD values sizes, the SQV of SKD values reckling in candidate site₄ ¹=1, The SQV of SKD the maximum in candidate site₄ ¹=10, the SQV of other candidate sites₄ ¹It is true by equipartition principle and linear interpolation method It is fixed.

Value rules and methods：Determined according to different SKF values sizes, the SQV of SKF values reckling in candidate site₄ ²=1, The SQV of SKF the maximum in candidate site₄ ²=10, the SQV of other candidate sites₄ ²It is true by equipartition principle and linear interpolation method It is fixed.

(3) underground water TDS, i.e. STW

Value rules and methods：If STW>100g/l, SQV₄ ³=1；If STW<1g/l, SQV₄ ³=10；If 1g/l ≤ STW≤100g/l, SQV₄ ³Determined according to equipartition principle and linear interpolation method.Between STW parameter values and suitability scale value Relation is as shown in table 21.

Table 21STW parameter values and suitability scale SQV₄ ³Value relation table

STW(g/l)	<1	1~13	13~25	25~37	37~49	49~61	61~73	73~85	85~100	>100
											SQV₄ ³	10	9	8	7	6	5	4	3	2	1

(4) groundwater age SAW

Value rules and methods：If 30 years<SAW≤1000 year, SQV₄ ⁴=3；If 1000<SAW≤10000 year, SQV₄ ⁴=6；If SAW>, SQV in 10000₄ ⁴=10.

(5) sulfate concentration SSC in underground water

Value rules and methods：If SSC<500mg/L, SQV₄ ⁵=10；If SSC>5000mg/L, SQV₄ ⁵=1；Such as Fruit 500mg/L≤SSC≤5000mg/L, SQV₄ ⁵Determine according to equisection method and linear interpolation method.SSC parameter values and suitability scale Relation between value is as shown in table 22.

Table 22SSC parameter values and suitability scale SQV₄ ⁵Value relation table

(6) chlorine ion concentration SCC in underground water

Value rules and methods：If SCC<1000mg/l, SQV₄ ⁶=10；If SCC>5000mg/l, SQV₄ ⁶=1；Such as Fruit 1000mg/l≤SCC≤5000mg/l, SQV₄ ⁶Determine according to equisection method and linear interpolation method.SCC parameter values and suitability mark Relation between angle value is as shown in table 23.

Table 23SCC parameter values and suitability scale SQV₄ ⁶Value relation table

(7) effective diffusion cofficient SDE of the radioactive substance in rock

Value rules and methods：SDE maximum SQV in candidate site₄ ⁷=10, SDE minimum SQV in candidate site₄ ⁷= 1, the SQV of other candidate sites₄ ⁷Determine by equipartition principle and linear interpolation method.

(8) solubility SSD of the radioactive substance in ground water regime

Value rules and methods：SSD maximum SQV in candidate site₄ ⁸=1, SSD minimum SQV in candidate site₄ ⁸= 10, the SQV of other candidate sites₄ ⁸Determine by equipartition principle and linear interpolation method.

5th, mankind's activity first class index

(1) site to available resources ground apart from SUD

Value rules and methods：It is more remote more suitable apart from these resources.If SUD=0, show that site has covered available Resource, the site is not suitable for.The candidate site maximum for SUD, SQV₅ ¹=10, the candidate site minimum for SUD, SQV₅ ¹ =1, the SQV of other candidate sites₅ ¹Determine by equipartition principle and linear interpolation method.

(2) probability SFH is flooded in the site that mankind's activity causes

Value rules and methods：The candidate site maximum for SFH, SQV₅ ²=1, the candidate site minimum for SFH, SQV₅ ²=10, the SQV of other candidate sites₅ ²Determine by equipartition principle and linear interpolation method.

6th, build and engineering specifications first class index

(1) rock strength and stress ratio SRP

Value rules and methods：Assuming that rock uniaxiality strength R_CWith the maximum principal stress σ for measuring₁Ratio represents with RP, That is RP=R_C/σ₁.When RP >=10, stress problem will not produce influence to engineering-built；As 5≤RP ＜ 10, can produce slight Palisades cracking phenomena；As 3≤RP ＜ 5, obvious palisades cracking phenomena can be produced；As RP ＜ 3, have Rock-burst and go out It is existing.The assignment of SRP indexs is determined on the basis of the statistical probability distribution of the RP values of certain candidate site Nei Ge testing site of this research use, And scale value is determined with the probability size of RP >=5.When the probability of RP >=5 in certain site is 1, i.e. SRP=P { RP │ RP_i≥ 5 }=1, SQV₆ ¹=10, wherein i=1,2 ..., K, K are the experiment quantity completed in certain site；As 0.1≤SRP ＜ 1, SQV₆ ¹=10 × P；As SRP ＜ 0.1, SQV₆ ¹=1.

(2) rock mass integrity index index S KV

Value rules and methods：Rock mass completeness be divided into it is complete, more complete, relatively broken, broken, smashed to pieces five etc. Level, its corresponding rock mass integrity index Kv scope is respectively：Kv>0.75,0.55<Kv≤0.75,0.35<Kv≤0.55, 0.15≤Kv≤0.35, Kv<0.15, its corresponding suitability scale SQV₆ ²Value is respectively 10,4,1,0.1,0.

(3) under temperature and effect of irradiation rock damaging stress value STR

Value rules and methods：The damaging stress value of rock characterizes the long-term strength of rock, pacifies for a long time for disposal storehouse Full property is significant.The candidate site maximum for STR, SQV₆ ³=10, the candidate site minimum for STR, SQV₆ ³=1, The SQV of other candidate sites₆ ³Determine by equipartition principle and linear interpolation method.

(4) under temperature and effect of irradiation rock creep rate SRR

Value rules and methods：The candidate site maximum for SRR, SQV₆ ⁴=1, the candidate site minimum for SRR, SQV₆ ⁴=10, the SQV of other candidate sites₆ ⁴Determine by equipartition principle and linear interpolation method.

(5) place is than SBG high

Value rules and methods：The candidate site maximum for SBG, SQV₆ ⁵=1, the candidate site minimum for SBG, SQV₆ ⁵=10, the SQV of other candidate sites₆ ⁵Determine by equipartition principle and linear interpolation method.

7th, waste traffic condition first class index SQV₇

(1) to spentnuclear fuel commercialization reprocessing plant apart from SDM

Value rules and methods：The candidate site maximum for SDM, SQV₇ ¹=1, the candidate site minimum for SDM, SQV₇ ¹=10, the SQV of other candidate sites₇ ¹Determine by equipartition principle and linear interpolation method.

(2) the distance SGD of reconstruction transport routes is increased newly

Value rules and methods：The candidate site maximum for SGD, SQV₇ ²=1, the candidate site minimum for SGD, SQV₇ ²=10, the SQV of other candidate sites₇ ²Determine by equipartition principle and linear interpolation method.

(3) waste freight STD

Value rules and methods：The candidate site maximum for STD, SQV₇ ³=1, the candidate site minimum for STD, SQV₇ ³=10, the SQV of other candidate sites₇ ³Determine by equipartition principle and linear interpolation method.

(4) the exposure dose of radiation SPR produced in transporting

Value rules and methods：The candidate site maximum for SPR, SQV₇ ⁴=1, the candidate site minimum for SPR, SQV₇ ⁴=10, the SQV of other candidate sites₇ ⁴Determine by equipartition principle and linear interpolation method.

8th, environmental protection condition first class index SQV₈

(1) to the Weighted distance index S DS of national nature reserve

Value rules and methods：Scope of statistics is 80km around candidate site.Assuming that 80km models around i-th candidate site There are national nature reserve m, its distance DS in enclosing_jRepresent；According to the power for from the close-by examples to those far off giving apart from site each protection zone Weight, uses k_jRepresent, distance is as shown in table 24 with the relation of weight；Assuming that quantity Y in protection zone in the corresponding distance of table 24 is interval_k Represent, wherein k is the interval quantity of distance, then Weighted distance index S DS of the i-th candidate site to national nature reserve_iTable Show as follows.SDS values are bigger, and candidate site is more suitable, its suitability scale value SQV₈It is bigger；If to national nature reserve away from From being zero, then show that candidate site in protection zone, is not suitable for.

Distance and weight mapping table of the site of table 24 to national nature reserve

During SDS=0, SQV₈ ¹=0；During SDS ≠ 0, the candidate site maximum for SDS, SQV₈ ¹=10, it is minimum for SDS Candidate site, SQV₈ ¹=1, other situations press equipartition principle and linear interpolation method determines.

(2) to the Weighted distance index S DA of provincial protection zone

Value rules and methods：Scope of statistics is 80km around candidate site.Assuming that 80km models around i-th candidate site There are provincial protection zone m, its distance DA in enclosing_jRepresent；According to the weight for from the close-by examples to those far off giving apart from site each protection zone, Use k_jRepresent, distance is as shown in Table 25 with the relation of weight；Assuming that quantity Y in protection zone in the corresponding distance of table 25 is interval_kTable Show, wherein k is the interval quantity of distance, then Weighted distance index S DA of the candidate site to provincial protection zone_iIt is expressed as follows.SDA Value is bigger, and candidate site is more suitable, and its suitability scale value is bigger；If the distance to provincial protection zone is zero, show candidate Site is not suitable in protection zone.

During SDA=0, SQV₈ ²=0；During SDA ≠ 0, the candidate site maximum for SDA, SQV₈ ²=10, it is minimum for SDA Candidate site, SQV₈ ²=1, other situations press equipartition principle and linear interpolation method determines.

Distance and weight mapping table of the site of table 25 to provincial protection zone

(3) to the Weighted distance index S DW of the important water body in downstream

Value rules and methods：Scope of statistics is 80km around candidate site.Assuming that i-th candidate site downstream peripheral There are important water body m, its distance DW in the range of 80km_jRepresent；According to apart from site downstream is from the close-by examples to those far off given, each is important The weight of water body, uses k_jRepresent, distance is as shown in table 26 with the relation of weight；Assuming that protection zone in the corresponding distance of table 26 is interval Quantity Y_kRepresent, wherein k is the interval quantity of distance, then Weighted distance index S DW of the candidate site to the important water body in downstream_i It is expressed as follows.SDW values are bigger, and candidate site is more suitable, and its scale value is bigger；If the Weighted distance to the important water body in downstream is Zero, then show that candidate site and water body are in same position, it is not suitable for.

During SDW=0, SQV₈ ³=0；During SDW ≠ 0, the candidate site maximum for SDW, SQV₈ ³=10；It is minimum for SDW Candidate site, SQV₈ ³=1, other situations press equipartition principle and linear interpolation method determines.

Weighted distance and weight mapping table of the site of table 26 to the important water body in downstream

9th, Land_use change condition first class index SQV₉

(1) land use pattern SLC

Value rules and methods：Land use pattern refers to《Law of Land Administration of the People's Republic of China》With《Land_use change shows Shape is classified》The land type of (GB/T 21010-2007) definition.Different land use pattern indicates Land_use change valency indirectly Value.If land type is arable land (01), field (02), forest land (03), meadow (043), SQV in farming land₉ ¹Value compared with It is small, such as SQV₉ ¹=5；If land use pattern be construction land in business take (05), industrial and mineral warehouse land (06), house Land used (07), public administration and public service land used (08), special use area (09), then SQV₉ ¹Value is 0；If land type is Other soils (124,125,126,127) in unused land, then SQV₉ ¹Value is 10；If land use pattern is except above-mentioned Type outside the type enumerated, then SQV₉ ¹Value is larger, such as SQV₉ ¹=7.

(2) Land_use change value SVL

Value rules and methods：It is high, medium and low and without value four etc. that Land_use change value can be divided into value Level, its suitability scale SQV₉ ²Value is followed successively by 1,3,6,10.

10th, social economy and humane condition first class index SQV₁₀

(1) population of ethnic minorities ratio SMP

Value rules and methods：Assuming that with site surrounding area radius 80km as area is evaluated, usually, evaluating in area The ratio SMP that the population of ethnic minorities account for the region total population is smaller more suitable.SMP parameter values and suitability scale value SQV₁₀ ¹It Between relation it is as follows.

(2) whether local government supports SCS

Value rules and methods：Investigation statisticses scope is counties and cities' level for being related in the range of 80km around site and with Shangdi Fang Zhengfu.If local government actively supports, SQV₁₀ ²=10；If local government does not oppose, SQV₁₀ ²=6；If local political affairs Mansion opposition, SQV₁₀ ²=0.

(3) the local public supports the ratio SPS of site

Value rules and methods：Permanent resident population to candidate site region carries out survey, and scope of statistics is field The counties and cities level administrative region being related in the range of 80km around location.If the local public supports the ratio SPS of site<50%, SQV₁₀ ³=0；If the local public supports ratio SPS >=50% of site, SQV₁₀ ³Equipartition principle and line are pressed according to different SPS values Property interpolation method determine.

Step 3. weight of Delphi method agriculture products

Using the weight of Delphi method agriculture products, weight assignment is realized by expert investigation, weight is suitable site The property indispensable |input paramete of quantitative evaluation.

Step 4. calculates the Hlw Geological Repository site final score value of suitability

The Hlw Geological Repository site final score value of suitability is calculated, this is the final output result of this method, is The integrated ordered foundation of Hlw Geological Repository site suitability.Present invention reference《World's core geological sciences》Journal article The computational methods and step proposed in " Preliminary study on index system of Hlw Geological Repository site suitability quantitative evaluation ", Calculate the Hlw Geological Repository site final score value of suitability.

Claims

1. a kind of Hlw Geological Repository site suitability evaluation methods based on Delphi method, it is characterised in that：Including Following steps：

Step 1. determines the index of Hlw Geological Repository site suitability quantitative assessment；

Step 2. builds the Codes and Standards of indices quantification；

Step 3. weight of Delphi method agriculture products；

Step 4. calculates the Hlw Geological Repository site final score value of suitability.

2. a kind of Hlw Geological Repository site suitability evaluation side based on Delphi method according to claim 1 Method, it is characterised in that：Index includes described in step 1：Site Selection stage index, Regional survey stage index and site are commented Valency and the stage of recognition index.

3. a kind of Hlw Geological Repository site suitability evaluation side based on Delphi method according to claim 2 Method, it is characterised in that：In step 2, Site Selection stage index suitability scale, region are represented respectively using PQV, QQV and SQV Investigation phase index suitability scale and site are evaluated and the stage of recognition index suitability scale；

PQV_i1 ^j1Represent the suitability scale of 1 two-level index of jth in the i-th 1 first class index of Site Selection stage；QQV_i2 ^j2Table Show the suitability scale of 2 two-level index of jth in the i-th 2 first class index of Regional survey stage；SQV_i3 ^j3Represent site evaluate and The suitability scale of 3 two-level index of jth in the i-th 3 first class index of the stage of recognition；

Step 2.1. Site Selection stages each index suitability scale value

I, geological conditions first class index PQV₁

(1) candidate's country rock thickness PRT

Value rules and methods：If argillaceous rocks and the continuously distributed thickness of rock salt geologic body are less than 100m, other rock class thickness are small In 500m, then suitability scale PQV₁ ¹=0；If argillaceous rocks and the continuously distributed thickness of rock salt geologic body are less than 100m, other rocks Class thickness is less than 500m, suitability scale PQV₁ ¹=10；

(2) candidate's country rock buried depth PRD

Value rules and methods：If candidate's country rock buried depth meets underground 500m~1000m depth boundses, PQV₁ ²=10；Otherwise, PQV₁ ²=0；

(3) candidate's country rock formation dip PFA

Value rules and methods：For non-laminar class candidate's Rock Mass, PQV₁ ³=10；For stratiform class geologic body, work as PFA< At 5 °, PQV₁ ³=10, as 5 °≤PFA<At 10 °, PQV₁ ³=5, other situations PQV₁ ³=0；

II, following natural trend condition first class index PQV₂

(1) the dynamic peak accelerator PSL of regional earthquake

Value rules and methods：If PSL>0.2g, PQV₂ ¹=0；If 0.15g<PSL≤0.2g, PQV₂ ¹=1；If 0.1g< PSL≤0.15g, PQV₂ ¹=4；If 0.05g<PSL≤0.1g, PQV₂ ¹=7；If PSL≤0.05g, PQV₂ ¹=10；

(2) to Regional Geotectonic secondary unit border apart from PGZ

Value rules and methods：If pre-selected zone border is to Regional Geotectonic secondary unit border apart from PGZ<25km, PQV₂ ²=0, i.e., should reduce the scope of pre-selected zone；If PGZ>100km, PQV₂ ²=10；If 25km≤PGZ≤100km, PQV₂ ²Determined according to equipartition principle and linear interpolation method；

(3) to the strong zone of influence of extreme geologic process apart from PGA

Value rules and methods：If PGA<25km, PQV₂ ³=0；If PGA>100km, PQV₂ ³=10；If 25km≤PGA ≤ 100km, PQV₂ ³Determined according to equipartition principle and linear interpolation method；

(4) natural geologic hazard development degree PND

Value rules and methods：Natural geologic hazard development degree is divided into strong, medium, weak and four grades of agensis, and this four The corresponding suitability scale PQV of individual grade₂ ⁴It is followed successively by 0,2,7,10；

III, geochemical conditions first class index PQV₃

(1) the main surrounding rock type PTR in pre-selected zone

Value rules and methods：When PTR is this four rocks of CRYSTALLINE ROCKS, argillaceous rocks, tufa stone and rock salt, PQV₃ ¹=10；It is no Then, PQV₃ ¹=0；

IV, environmental protection condition first class index PQV₄

(1) to Head streams protection zone apart from PSS

Value rules and methods：If PSS=0, PQV₄ ¹=0；If PSS>80km, then PQV₄ ¹=10；If 0<PSS≤ 80km, according to equipartition principle and linear interpolation principle value；

(2) the quantity PEP of environmental sensitive area above the provincial level

Value rules and methods：Count to preselect district center as the provincial and state of laws and regulations regulation in round dot 80km radius The quantity of family's level environmental sensitive area, if PEP<20, then PQV₄ ²=10；If PEP>500, then PQV₄ ²=1；If 20 Individual≤PEP≤500, are determined by equipartition principle and linear interpolation method；(3) the quantity PWP of important water body

Value rules and methods：Statistics is to preselect district center to specify in downstream direction laws and regulations in round dot 80km radius Provincial and national important water body quantity, if PWP<10, then PQV₄ ³=10；If PWP>100, then PQV₄ ³= 1；If 10≤PWP≤100, determined by equipartition principle and linear interpolation method；(4) pre-selected zone absolute height above sea level PHE

Value rules and methods：If the height above sea level of pre-selected zone highest point is more than 3500m, PQV₄ ⁴=0；Otherwise, PQV₄ ⁴=10；

V, social economy and humane condition first class index PQV₅

(1) the absolute density of population PAP in region

Value rules and methods：The absolute density of population in region refers to resident population and floating population's summation and pre-selected zone in pre-selected zone Area ratio, if PAP >=400 people/km², PQV₅ ¹=0；If 200 people/km²≤PAP<400 people/km², PQV₅ ¹=1；If 100 people/km²≤PAP<200 people/km², PQV₅ ¹=3；If 20 people/km²≤PAP<100 people/km², PQV₅ ¹=6；If PAP< 20 people/km², PQV₅ ¹=10；

(2) Main Towns quantity PTW

Value rules and methods：Main Towns quantity refers to medium and above scale cities and towns the quantity included in pre-selected zone, if PTW≤5, PQV₅ ²=10；If 5<PTW≤50, PQV₅ ²=7；If 50<PTW≤100, PQV₅ ²=3；If PTW>100, PQV₅ ²=1；

(3) region total output value PDP per capita

Value rules and methods：The production per capita in provincial administrative area is total where this phase zone per capita total output value pre-selected zone Value is characterized, if pre-selected zone is averagely added up afterwards across multiple provincial administrative areas according to area；Assuming that the whole nation per capita total output value use PDP_nRepresent, Provincial administrative area per capita total output value come top 10 total output value per capita average value PDP₁₀Table Show, if PDP≤PDP_n, PQV₅ ³=10；If PDP >=PDP₁₀, PQV₅ ³=1；In the case of other, PQV₅ ³=5；

Step 2.2. Regional survey stages each index suitability scale value

1st, geological conditions first class index QQV₁

(1) candidate's country rock area QGA

Value rules and methods：If QGA<50km², QQV₁ ¹=1；If QGA>300km², then QQV₁ ¹=10；If 50km²≤ QGA≤300km², QQV₁ ¹According to equipartition principle and linear interpolation principle value；(2) candidate's country rock thickness QRT

Value rules and methods：If QRT<100m, QQV₁ ²=0；If QRT>500m, QQV₁ ²=10；If 100m≤QRT≤ 500m, QQV₁ ²According to equipartition principle and linear principle value；(3) candidate's country rock buried depth QRD

Value rules and methods：If candidate's country rock buried depth meets underground 500m~1000m depth boundses, QQV₁ ³=10；Otherwise, QQV₁ ³=0；

(4) rift structure distribution surface density QSD

Value rules and methods：Rift structure refers to the rift structure for cutting and passing through pre-selection location, rift structure distribution surface density Refer to the rift structure quantity in unit area under the conditions of uniform or more uniform distribution；If QSD<1/km², QQV₁ ⁴=10；Such as Fruit QSD>10/km², QQV₁ ⁴=1；If 1/km²≤ QSD≤100/km², QQV₁ ⁴According to equipartition principle and linear principle Value；II, hydrogeologic condition first class index QQV₂

(1) groundwater level depth QWB

Value rules and methods：Assuming that the buried depth in proposed disposal storehouse is MZ, if QWB≤MZ, QQV₂ ¹=5；Otherwise, QQV₂ ¹ =10；

(2) underground water hydraulic gradient QHI

Value rules and methods：If QHI<1 ‰, QQV₂ ²=10；If QHI>100 ‰, QQV₂ ²=0；If 1 ‰≤QHI≤ 100 ‰, QQV₂ ²According to equipartition principle and linear principle value；(3) rock mass ground infiltration coefficient accounting QK

Value rules and methods：Using mathematical statistics method, with rock mass hyposmosis coefficient≤1 × 10^-8M/s proportion sizes are true Determine QQV₂ ³Value；When hyposmosis coefficient≤1 × 10 of rock mass in certain location^-8Accounting P { QK │ QK≤1 × 10 of m/s^-8M/s }=1 When, QQV₂ ³=10；When 0.01<P{QK│QK≤1×10^-8M/s } ＜ 1 when, QQV₂ ³=10 × P；When P { QK │ QK≤1 × 10^-8m/ S }≤0.01 when, QQV₂ ³=0.1；

(4) to region discharge area apart from QPD

Value rules and methods：If QPD<50km, QQV₂ ⁴=1；If QPD>200km, QQV₂ ⁴=10；If 50km≤QPD ≤ 200km, QQV₂ ⁴By equipartition principle and linear principle value；(5) the transmissibility QTF of rift structure

Value rules and methods：Important Faults construction refers to the rift structure for cutting and passing through rock mass；In this stage, according to fracture structure The property and charges feature qualitative discrimination made, be divided into water transmitting ability it is poor, in, strong three kinds of situations, its corresponding suitability scale QQV₂ ⁵Value is followed successively by 10,1,0.1；

III, geochemical conditions first class index QQV₃

(1) underground water TDS, i.e. QTW

Value rules and methods：If QTW>100g/l, QQV₃ ¹=1；If QTW<1g/l, QQV₃ ¹=10；If 1g/l≤QTW ≤ 100g/l, QQV₃ ¹Determined according to equipartition principle and linear interpolation method；

(2) underground water Eh values, i.e. QEH

Value rules and methods：If QEH=0, QQV₃ ²=6；If QEH<0, QQV₃ ²=10；If QEH>0, QQV₃ ²=1；

(3) groundwater age QAW

Value rules and methods：If QAW≤30 year, QQV₃ ³=0；If 30 years<QAW≤1000 year, QQV₃ ³=3；If 1000<QAW≤10000 year, QQV₃ ³=6；If QAW>, QQV in 10000₃ ³=10；

(4) sulfate concentration QSC in underground water

Value rules and methods：If QSC<500mg/L, QQV₃ ⁴=10；If QSC>5000mg/L, QQV₃ ⁴=1；If 500mg/L≤QSC≤5000mg/L, QQV₃ ⁴Determine according to equipartition principle and linear interpolation method；(5) chlorine ion concentration in underground water QCC

Value rules and methods：If QCC<1000mg/l, QQV₃ ⁵=10；If QCC>5000mg/l, QQV₃ ⁵=1；If 1000mg/l≤QCC≤5000mg/l, QQV₃ ⁵Determine according to equipartition principle and linear interpolation method；IV, environmental protection condition one-level Index QQV₄

(1) the quantity QEP of environmental sensitive area above the provincial level

Value rules and methods：Count the provincial and country with Section Center as laws and regulations regulation in round dot 25km radius The quantity of level environmental sensitive area；If QEP<5, QQV₄ ¹=10；If QEP>50, QQV₄ ¹=1；5≤QEP≤50, QQV₄ ¹Determine by equipartition principle and linear interpolation method；(2) the quantity QWP of important water body

Value rules and methods：Count with preselect Section Center in round dot 25km radius laws and regulations specify it is provincial and The quantity of national important water body；If QWP<5, QQV₄ ²=10；If QWP>20, QQV₄ ²=1；If 5≤QWP ≤ 20, QQV₄ ²Determine by equipartition principle and linear interpolation method；V, social economy and humane condition first class index QQV₅

(1) the absolute density of population QAP in region

Value rules and methods：The absolute density of population QAP in region refers to resident in the range of with Section Center as center of circle radius 15km Population and floating population's summation and area ratio, QAP is more big more is not suitable for addressing；If QAP<5 people/km², QQV₅ ¹=10；Such as Fruit QAP>200 people/km², QQV₅ ¹=0；If 5 people/km²≤ QAP≤200 people/km², QQV₅ ¹By equipartition principle and linear interpolation Method determines；(2) Main Towns quantity QTW

Value rules and methods：Main Towns quantity QTW refer to medium in the range of with Section Center as center of circle radius 15km and with The quantity in upper scale cities and towns；If QTW>50, QQV₅ ²=1；If 10<QTW≤50, QQV₅ ²=3；If 3<QTW≤10, then QQV₅ ²=6；If QTW≤3, QQV₅ ²=10；

(3) total output value QDP per capita

Value rules and methods：This stage, total output value QDP referred to that the production per capita of local administrative area where pre-selection location is total per capita Value, if pre-selection location is averagely added up afterwards across multiple local administrative areas according to area；Assuming that national local administrative area is per capita The average value of total output value is QDP_x, in national local administrative area per capita total output value to come the production per capita of 50th reciprocal total Value QDP_-50Represent；If QDP>2QDP_x, QQV₅ ³=1；If QDP_x<QDP≤2QDP_x, QQV₅ ³=3；If QDP_-50<QDP ≤QDP_x, QQV₅ ³=6；If QDP≤QDP_-50, QQV₅ ³=10；

1st, geological conditions first class index SQV₁

(1) candidate's country rock area SGA

Value rules and methods：If candidate's country rock area SGA<10km², SQV₁ ¹=1；If 10km²≤SGA<20km², SQV₁ ¹ =3；If 20km²≤SGA<50km², SQV₁ ¹=6；If SGA >=50km², SQV₁ ¹=10；

(2) candidate's country rock thickness SRT

Value rules and methods：If candidate's country rock thickness SRT<100m, SQV₁ ²=0；If 100m≤SRT<200m, SQV₁ ²= 3；If 200m≤SRT<400m, SQV₁ ²=6；If SRT >=400m, SQV₁ ²=10；

(3) rock thermal conductivity factor SCT

Value rules and methods：The SQV of SCT the maximum in candidate site₁ ³=10, the SQV of SCT recklings in candidate site₁ ³= 1, the SQV of other candidate sites₁ ³Determine by equipartition principle and linear interpolation method；

(4) rift structure distribution surface density SSD

Value rules and methods：This stage only counts the rift structure for cutting and passing through candidate site, if rift structure is distributed Surface density SSD>5/km², SQV₁ ⁴=1；If 5/km²≥SSD>1/km², SQV₁ ⁴=3；If 1/km²≥SSD>0.2/km², SQV₁ ⁴=6；If SSD≤0.2/km², SQV₁ ⁴=10；

(5) to the Weighted distance SFD of Important Faults construction

Value rules and methods：So-called Important Faults construction refers to the rift structure for cutting and having passed through pre-selection location, in this stage Scope of statistics is 25km scopes around candidate site；Assuming that there is Important Faults to construct m around i-th candidate site in the range of 25km It is individual, its distance DS_jRepresent；The weight that every Important Faults are constructed from the close-by examples to those far off is given according to apart from site, k is used_jRepresent, away from From DS_jWith weight k_jRelation be shown below：

k_{j} = \{\begin{matrix} 0.01 & {DS}_{j} \leq 1 \\ 0.03 & 1 < {DS}_{j} \leq 2 \\ 0.04 & 2 < {DS}_{j} \leq 3 \\ 0.06 & 3 < {DS}_{j} \leq 4 \\ 0.1 & 4 < {DS}_{j} \leq 5 \\ 0.13 & 5 < {DS}_{j} \leq 7 \\ 0.17 & 7 < {DS}_{j} \leq 10 \\ 0.21 & 10 < {DS}_{j} \leq 15 \\ 0.26 & 15 < {DS}_{j} \leq 25 \end{matrix}

If apart from k_jInterior Important Faults construction quantity Y_kRepresent, wherein k is the interval quantity of distance, then i-th candidate site is arrived The Weighted distance index S FD of Important Faults construction_iIt is expressed as follows：

{SFD}_{i} = Σ_{k = 1}^{9} (Σ_{i = 1}^{Y_{k}} {DS}_{j} \cdot k_{j}) / Y_{k};

For all candidate sites, the SQV of SFD the maximum₁ ⁵The SQV of=10, SFD reckling₁ ⁵=1, the SQV of other sites₁ ⁵Press Equipartition principle and linear interpolation method determine；

II, following natural trend condition first class index SQV₂

(1) strong geologic process probability of happening SPG

Value rules and methods：Assuming that the probability that a certain geologic process of j-th candidates site occurs is P_i, i=1,2,3 ..., K, then select these geologic process probability of happening the maximum to represent the probability that strong geologic process occurs in the candidate site, i.e. SPG_j =Max (P_i, i=1,2,3 ..., K)；It is required that SPG_jLess than a ten thousandth, when meeting less than a ten thousandth, j-th candidates The SQV of location₂ ¹Value take 10；The SQV of other sites₂ ¹Determine by equipartition principle and linear interpolation method；

(2) extreme meteorology and climate disaster probability of happening SPT

Value rules and methods：Assuming that the probability that a certain extreme meteorology and climate disaster of j-th candidates site occurs is P_i, i=1, 2,3 ..., K, then select these meteorology and climate disaster probability of happening the maximum to represent the candidate site and extreme meteorology and climate calamity occur Harmful probability, i.e. SPT_j=Max (P_i, i=1,2,3 ..., K)；It is required that SPT_jThe probability of generation is less than 1/150th, when SPT_jWhen the probability of generation is less than 1/150th, the SQV of j-th candidates site₂ ²Value is 10, works as SPT_jWhat is occurred is general When rate is 100 percent, SQV₂ ²Value is 0, the SQV of other sites₂ ²Determine by equipartition principle and linear interpolation method；

(3) site to active fault apart from SMP

Value rules and methods：If SMP≤5km, SQV₂ ³=0；If 5km<SMP≤10km, SQV₂ ³=1；If 10km< SMP≤50km, SQV₂ ³=3；If 50km<SMP≤100km, SQV₂ ³=6；If SMP>100km, SQV₂ ³=10；

(4) place SL-2 grades of acceleration peak SSL

Value rules and methods：The SQV of SSL the maximum in candidate site₂ ⁴=1, the SQV of SSL recklings in candidate site₂ ⁴= 10, the SQV of other candidate sites₂ ⁴Value presses equipartition principle and linear interpolation method determines；

Value rules and methods：If being site near field in the range of the 25km of candidate site border periphery, if the field of candidate site Location near field periphery has and occurs to be not less than 7 grades of seismic structure bands of earthquake and point to the candidate site, then from SZX angles Consider that the candidate site is not suitable for as site；Otherwise, consider that the candidate site is suitable as site from SZX angles；Site is near Area periphery nearby points to site without earthquake magnitude seismogenic fault band high, then SQV₂ ⁵Value is 10, and site near field periphery has height Earthquake magnitude seismogenic fault band points to site, then value is 0, the SQV of other candidate sites₂ ⁵Value presses equipartition principle and linear interpolation method is true It is fixed；

(6) site SZJ is pointed to without seismogenic fault in the near field of site

Value rules and methods：If being site near field in the range of the 25km of candidate site border periphery, if the field of candidate site There is seismogenic fault in the near field of location and point to candidate site, then consider that the candidate site is not suitable for as site from SZJ angles；It is no Then, consider that the candidate site is suitable as site from SZJ angles；In the near field of site site, SQV are pointed to without seismogenic fault₂ ⁶Take It is 10 to be worth, and has seismogenic fault to point to site, SQV in the near field of site₂ ⁶Value is 0, the SQV of other candidate sites₂ ⁶Value presses decile Principle and linear interpolation method determine；

III, hydrogeologic condition first class index SQV₃

(1) underground water hydraulic gradient SHI

Value rules and methods：If SHI<1 ‰, SQV₃ ¹=10；If SHI>100 ‰, SQV₃ ¹=0；If 1 ‰≤SHI≤ 100 ‰, SQV₃ ¹According to equipartition principle and linear principle value；(2) rock mass ground infiltration coefficient accounting SK

Value rules and methods：Using mathematical statistics method, with Rockmass Permeabuity Coefficient≤1 × 10^-9The proportion size of m/s is true Determine SQV₃ ²Value；When infiltration coefficient SK≤1 × 10 of rock mass in certain site^-9Accounting P { SK │ SK≤1 × 10 of m/s^-9M/s }=1 When, SQV₃ ²=10；When accounting meets 0.01<P{SK│SK≤1×10^-9M/s } ＜ 1 when, SQV₃ ²=10 × P；When accounting meets P {SK│SK≤1×10^-9M/s }≤0.01 when, SQV₃ ²=0.1；

(3) to region discharge area apart from SPD

Value rules and methods：If SPD<50km, SQV₃ ³=1；If SPD>200km, SQV₃ ³=110；If 50km≤SPD ≤ 200km, SQV₃ ³According to equipartition principle and linear principle value；

(4) the coefficient of transmissibility STF of rift structure

Value rules and methods：In this stage, emphasis considers the transmissibility of the rift structure for cutting and passing through candidate site, obtains Quantitative data；Using mathematical statistics method, with coefficient of transmissibility≤1 × 10 of rift structure^-5m²The probability size of/s determines that certain is waited The suitability scale SQV of selected scenes location₃ ⁴Value；When STF≤1 × 10 in certain candidate site^-5m²/ s probability P STF │ STF≤1 × 10^-5m²/ s }=1 when, SQV₃ ⁴=10；When probability 0.01<P{STF│STF≤1×10^-5m²/ s } ＜ 1 when, SQV₃ ⁴=10 × P；When Probability P { STF │ STF≤1 × 10^-5m²/ s }≤0.01 when, SQV₃ ⁴=0.1；

Value rules and methods：The important water body in the range of the downstream 80km of candidate site is considered, if disposal storehouse candidate site institute There is hydraulic connection in ground water regime and downstream water, then site is not suitable for, SQV₃ ⁵Value is 0；Otherwise, the suitable SQV in site₃ ⁵ Value is 10, the SQV of other candidate sites₃ ⁵Determine by equipartition principle and linear interpolation method；

IV, geochemical conditions first class index SQV₄

Value rules and methods：The SQV of SKD values reckling in candidate site₄ ¹=1, the SQV of SKD the maximum in candidate site₄ ¹= 10, the SQV of other candidate sites₄ ¹Determine by equipartition principle and linear interpolation method；

Value rules and methods：The SQV of SKF values reckling in candidate site₄ ²=1, the SQV of SKF the maximum in candidate site₄ ²= 10, the SQV of other candidate sites₄ ²Determine by equipartition principle and linear interpolation method；

(3) underground water TDS, i.e. STW

Value rules and methods：If STW>100g/l, SQV₄ ³=1；If STW<1g/l, SQV₄ ³=10；If 1g/l≤STW ≤ 100g/l, SQV₄ ³Determined according to equipartition principle and linear interpolation method；(4) groundwater age SAW

Value rules and methods：If 30 years<SAW≤1000 year, SQV₄ ⁴=3；If 1000<SAW≤10000 year, SQV₄ ⁴ =6；If SAW>, SQV in 10000₄ ⁴=10；

(5) sulfate concentration SSC in underground water

Value rules and methods：If SSC<500mg/L, SQV₄ ⁵=10；If SSC>5000mg/L, SQV₄ ⁵=1；If 500mg/L≤SSC≤5000mg/L, SQV₄ ⁵Determine according to equipartition principle and linear interpolation method；(6) chlorine ion concentration in underground water SCC

Value rules and methods：If SCC<1000mg/l, SQV₄ ⁶=10；If SCC>5000mg/l, SQV₄ ⁶=1；If 1000mg/l≤SCC≤5000mg/l, SQV₄ ⁶Determine according to equipartition principle and linear interpolation method；(7) radioactive substance is in rock In effective diffusion cofficient SDE

Value rules and methods：SDE maximum SQV in candidate site₄ ⁷=10, SDE minimum SQV in candidate site₄ ⁷=1, its The SQV of his candidate site₄ ⁷Determine by equipartition principle and linear interpolation method；

(8) solubility SSD of the radioactive substance in ground water regime

Value rules and methods：SSD maximum SQV in candidate site₄ ⁸=1, SSD minimum SQV in candidate site₄ ⁸=10, its The SQV of his candidate site₄ ⁸Determine by equipartition principle and linear interpolation method；

V, mankind's activity first class index

(1) site to available resources ground apart from SUD

Value rules and methods：The candidate site maximum for SUD, SQV₅ ¹=10, the candidate site minimum for SUD, SQV₅ ¹ =1, the SQV of other candidate sites₅ ¹Determine by equipartition principle and linear interpolation method；

(2) probability SFH is flooded in the site that mankind's activity causes

Value rules and methods：The candidate site maximum for SFH, SQV₅ ²=1, the candidate site minimum for SFH, SQV₅ ²= 10, the SQV of other candidate sites₅ ²Determine by equipartition principle and linear interpolation method；

VI, construction and engineering specifications first class index

(1) rock strength and stress ratio SRP

Value rules and methods：Assuming that rock uniaxiality strength R_CWith the maximum principal stress σ for measuring₁Ratio represents with RP, i.e. RP =R_C/σ₁；When the probability of RP >=5 in certain site is 1, i.e. SRP=P { RP │ RP_i>=5 }=1, SQV₆ ¹=10, wherein i=1, 2 ..., K, K are the experiment quantity completed in certain site；As 0.1≤SRP ＜ 1, SQV₆ ¹=10 × P；As SRP ＜ 0.1, SQV₆ ¹=1；

(2) rock mass integrity index index S KV

Value rules and methods：Rock mass completeness is divided into complete, more complete, relatively broken, broken, smashed to pieces five grades, its Corresponding rock mass integrity index Kv scopes are respectively：Kv>0.75,0.55<Kv≤0.75,0.35<Kv≤0.55,0.15≤Kv ≤ 0.35, Kv<0.15, its corresponding suitability scale SQV₆ ²Value is respectively 10,4,1,0.1,0；

(3) under temperature and effect of irradiation rock damaging stress value STR

Value rules and methods：The candidate site maximum for STR, SQV₆ ³=10, the candidate site minimum for STR, SQV₆ ³ =1, the SQV of other candidate sites₆ ³Determine by equipartition principle and linear interpolation method；

(4) under temperature and effect of irradiation rock creep rate SRR

Value rules and methods：The candidate site maximum for SRR, SQV₆ ⁴=1, the candidate site minimum for SRR, SQV₆ ⁴= 10, the SQV of other candidate sites₆ ⁴Determine by equipartition principle and linear interpolation method；

(5) place is than SBG high

Value rules and methods：The candidate site maximum for SBG, SQV₆ ⁵=1, the candidate site minimum for SBG, SQV₆ ⁵= 10, the SQV of other candidate sites₆ ⁵Determine by equipartition principle and linear interpolation method；

VII, waste traffic condition first class index SQV₇

(1) to spentnuclear fuel commercialization reprocessing plant apart from SDM

Value rules and methods：The candidate site maximum for SDM, SQV₇ ¹=1, the candidate site minimum for SDM, SQV₇ ¹= 10, the SQV of other candidate sites₇ ¹Determine by equipartition principle and linear interpolation method；

(2) the distance SGD of reconstruction transport routes is increased newly

Value rules and methods：The candidate site maximum for SGD, SQV₇ ²=1, the candidate site minimum for SGD, SQV₇ ²= 10, the SQV of other candidate sites₇ ²Determine by equipartition principle and linear interpolation method；

(3) waste freight STD

Value rules and methods：The candidate site maximum for STD, SQV₇ ³=1, the candidate site minimum for STD, SQV₇ ³= 10, the SQV of other candidate sites₇ ³Determine by equipartition principle and linear interpolation method；

(4) the exposure dose of radiation SPR produced in transporting

Value rules and methods：The candidate site maximum for SPR, SQV₇ ⁴=1, the candidate site minimum for SPR, SQV₇ ⁴= 10, the SQV of other candidate sites₇ ⁴Determine by equipartition principle and linear interpolation method；

VIII, environmental protection condition first class index SQV₈

(1) to the Weighted distance index S DS of national nature reserve

Value rules and methods：Scope of statistics is 80km around candidate site, if having in the range of 80km around i-th candidate site National nature reserve m, its distance DS_jRepresent；According to the weight for from the close-by examples to those far off giving apart from site each protection zone, k is used_j Represent, apart from DS_jWith weight k_jRelation be shown below：

k_{j} = \{\begin{matrix} 0.01 & {DS}_{j} \leq 5 \\ 0.03 & 5 < {DS}_{j} \leq 10 \\ 0.04 & 10 < {DS}_{j} \leq 20 \\ 0.06 & 20 < {DS}_{j} \leq 30 \\ 0.1 & 30 < {DS}_{j} \leq 40 \\ 0.13 & 40 < {DS}_{j} \leq 50 \\ 0.17 & 50 < {DS}_{j} \leq 60 \\ 0.21 & 60 < {DS}_{j} \leq 70 \\ 0.26 & 70 < {DS}_{j} \leq 80 \end{matrix}

If apart from k_jInterior protection zone quantity Y_kRepresent, wherein k is the interval quantity of distance, then i-th candidate site is to national level The Weighted distance index S DS of protection zone_iIt is expressed as follows：

{SDS}_{i} = Σ_{k = 1}^{9} (Σ_{i = 1}^{Y_{k}} {DS}_{j} \cdot k_{j}) / Y_{k};

During SDS=0, SQV₈ ¹=0；During SDS ≠ 0, the candidate site maximum for SDS, SQV₈ ¹=10, the time minimum for SDS Selected scenes location, SQV₈ ¹=1, other situations press equipartition principle and linear interpolation method determines；

(2) to the Weighted distance index S DA of provincial protection zone

Value rules and methods：Scope of statistics is 80km around candidate site, if having in the range of 80km around i-th candidate site Provincial protection zone m, its distance DA_jRepresent；According to the weight for from the close-by examples to those far off giving apart from site each protection zone, k is used_jTable Show, distance is shown below with the relation of weight：

k_{j} = \{\begin{matrix} 0.01 & {DS}_{j} \leq 5 \\ 0.03 & 5 < {DS}_{j} \leq 10 \\ 0.04 & 10 < {DS}_{j} \leq 20 \\ 0.06 & 20 < {DS}_{j} \leq 30 \\ 0.1 & 30 < {DS}_{j} \leq 40 \\ 0.13 & 40 < {DS}_{j} \leq 50 \\ 0.17 & 50 < {DS}_{j} \leq 60 \\ 0.21 & 60 < {DS}_{j} \leq 70 \\ 0.26 & 70 < {DS}_{j} \leq 80 \end{matrix}

If apart from k_jInterior protection zone quantity Y_kRepresent, wherein k is the interval quantity of distance, then the candidate site to provincial protection zone Weighted distance index S DA_iIt is expressed as follows：

{SDS}_{i} = Σ_{k = 1}^{9} (Σ_{i = 1}^{Y_{k}} {DA}_{j} \cdot k_{j}) / Y_{k}

During SDA=0, SQV₈ ²=0；During SDA ≠ 0, the candidate site maximum for SDA, SQV₈ ²=10, the time minimum for SDA Selected scenes location, SQV₈ ²=1, other situations press equipartition principle and linear interpolation method determines；(3) to the Weighted distance of the important water body in downstream Index S DW

Value rules and methods：Scope of statistics is 80km around candidate site, it is assumed that i-th candidate site downstream peripheral 80km model There are important water body m, its distance DW in enclosing_jRepresent；Each important water body of downstream is from the close-by examples to those far off given according to apart from site Weight, uses k_jRepresent, distance is shown below with the relation of weight：

k_{j} = \{\begin{matrix} 0.01 & {DS}_{j} \leq 5 \\ 0.03 & 5 < {DS}_{j} \leq 10 \\ 0.04 & 10 < {DS}_{j} \leq 20 \\ 0.06 & 20 < {DS}_{j} \leq 30 \\ 0.1 & 30 < {DS}_{j} \leq 40 \\ 0.13 & 40 < {DS}_{j} \leq 50 \\ 0.17 & 50 < {DS}_{j} \leq 60 \\ 0.21 & 60 < {DS}_{j} \leq 70 \\ 0.26 & 70 < {DS}_{j} \leq 80 \end{matrix}

If apart from k_jInterior protection zone quantity Y_kRepresent, wherein k is the interval quantity of distance, then the candidate site is important to downstream is arrived The Weighted distance index S DW of water body_iIt is expressed as follows：

{SDS}_{i} = Σ_{k = 1}^{9} (Σ_{i = 1}^{Y_{k}} {DW}_{j} \cdot k_{j}) / Y_{k}

During SDW=0, SQV₈ ³=0；During SDW ≠ 0, the candidate site maximum for SDW, SQV₈ ³=10；The time minimum for SDW Selected scenes location, SQV₈ ³=1, other situations press equipartition principle and linear interpolation method determines；

IX, Land_use change condition first class index SQV₉

(1) land use pattern SLC

Value rules and methods：If land type is arable land, field, forest land, the meadow in farming land, SQV₉ ¹=5；If native Ground use pattern is that the business in construction land takes ground, industrial and mineral warehouse land, residential land, public administration and public service use Ground, special use area, SQV₉ ¹=0；If land type is other soils in unused land, SQV₉ ¹=10；If Land_use change Type is the type in addition to the above-mentioned type enumerated, SQV₉ ¹=7；

(2) Land_use change value SVL

Value rules and methods：Land_use change value is divided into that value is high, medium and low and without four grades of value, and it is fitted Suitable property scale SQV₉ ²Value is followed successively by 1,3,6,10；

X, social economy and humane condition first class index SQV₁₀

(1) population of ethnic minorities ratio SMP

Value rules and methods：Assuming that being to evaluate area, SMP parameter values and suitability scale with site surrounding area radius 80km Value SQV₁₀ ¹Between relation it is as follows:

{SQV}_{9}^{1} = \{\begin{matrix} 10 & S M P < 8.5 % \\ 7 & 8.5 % \leq S M P < 23.3 % \\ 5 & 23.3 % \leq S M P < 50 % \\ 3 & 50 % \leq S M P < 66.7 % \\ 1 & 66.75 % \leq S M P < 85 % \\ 0.1 & S M P \leq 85 % \end{matrix}

(2) whether local government supports SCS

Value rules and methods：Investigation statisticses scope is the counties and cities' level and local political affairs above being related in the range of 80km around site Mansion；If local government actively supports, SQV₁₀ ²=10；If local government does not oppose, SQV₁₀ ²=6；If local government is anti- It is right, SQV₁₀ ²=0；

(3) the local public supports the ratio SPS of site

Value rules and methods：Permanent resident population to candidate site region carries out survey, and scope of statistics is site week Enclose the counties and cities level administrative region being related in the range of 80km；If the local public supports the ratio SPS of site<50%, SQV₁₀ ³= 0；If the local public supports ratio SPS >=50% of site, SQV₁₀ ³Equipartition principle and linear interpolation are pressed according to different SPS values Method determines.