• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    放射性廢物處置研究進(jìn)展

    2011-04-03 04:52:21易樹平馬海毅鄭春苗
    地球?qū)W報(bào) 2011年5期
    關(guān)鍵詞:核素運(yùn)移廢物

    易樹平, 馬海毅, 鄭春苗

    1)廣東省電力設(shè)計(jì)研究院, 廣東廣州 510663;

    2)北京大學(xué)水資源研究中心, 北京 100871

    放射性廢物處置研究進(jìn)展

    易樹平1,2), 馬海毅1), 鄭春苗2)

    1)廣東省電力設(shè)計(jì)研究院, 廣東廣州 510663;

    2)北京大學(xué)水資源研究中心, 北京 100871

    放射性廢物的處置是制約核能可持續(xù)發(fā)展的關(guān)鍵因素, 目前已成為國際社會(huì)關(guān)注的熱點(diǎn)問題之一。針對(duì)處置場地核素運(yùn)移污染的風(fēng)險(xiǎn)問題, 對(duì)放射性廢物的處置及其選址、核素運(yùn)移試驗(yàn)和核素運(yùn)移模型進(jìn)行了回顧和論述。指出采用多重屏障系統(tǒng)進(jìn)行放射性廢物的處置, 其安全性是可以得到保障的; 處置場的選址應(yīng)遵循就近原則, 并應(yīng)從環(huán)境水文地質(zhì)的角度來構(gòu)建場地的選址及評(píng)價(jià)體系; 核素運(yùn)移試驗(yàn)應(yīng)注重室內(nèi)試驗(yàn)和原位試驗(yàn)相結(jié)合, 以保證所得結(jié)果和參數(shù)的可靠性; 核素運(yùn)移模型研究經(jīng)歷了由對(duì)流-彌散模型→Kd模型→多組分反應(yīng)性模型的發(fā)展階段, 并取得了重要進(jìn)展, 但在處理非均質(zhì)介質(zhì)和非均質(zhì)反應(yīng)、礦物溶解-沉淀與水動(dòng)力行為的相互影響以及穩(wěn)定性和實(shí)際應(yīng)用能力等方面還亟待提高。我國放射性廢物的處置與研究較為滯后, 需在處置場選址規(guī)范體系、核素運(yùn)移試驗(yàn)以及模擬技術(shù)等方面取得突破和提高。

    放射性廢物處置; 水文地質(zhì); 環(huán)境影響評(píng)價(jià); 核素運(yùn)移試驗(yàn); 核素運(yùn)移模型

    20世紀(jì)40年代以來, 核科學(xué)的飛速發(fā)展和應(yīng)用 產(chǎn)生了大量的放射性廢物。由于放射性核廢物對(duì)人體甚至整個(gè)生態(tài)系統(tǒng)具有極大的危害性, 其影響可達(dá)幾百年到數(shù)萬年甚至更長的時(shí)間, 因此世界各國對(duì)放射性廢物處置的研究都非常重視(Bradbury et al., 1985; Kittel, 1989)。

    放射性廢物處置的任務(wù)是在廢物可能對(duì)人類造成不可接受的危險(xiǎn)的時(shí)間內(nèi), 將廢物中的放射性核素限制在處置場范圍內(nèi), 防止核素以不可接受的濃度或數(shù)量向環(huán)境釋放而影響人類的健康與安全(Han et al., 1997; Kickmaier et al., 1997)。然而, 盡管半個(gè)多世紀(jì)以來放射性廢物的處置和研究取得了較大的進(jìn)展, 但在場地選址、參數(shù)獲取及污染風(fēng)險(xiǎn)評(píng)價(jià)等方面仍存在不足, 增加了放射性廢物安全處置的不確定性。

    目前, 世界各國正積極推行低碳經(jīng)濟(jì), 清潔能源的呼聲越來越高, 中國正“積極發(fā)展核電”。因此,全面和正確地認(rèn)識(shí)放射性廢物的處置就尤為重要。但是, 我國的放射性廢物處置還處于起步階段, 因此有必要對(duì)國際上放射性廢物的處置概況進(jìn)行回顧,總結(jié)處置經(jīng)驗(yàn), 探討處置過程中涉及的核心課題及其研究現(xiàn)狀就具有重要理論和實(shí)際意義。

    1 放射性廢物的處置

    1.1 歷史回顧

    人類歷史上第一次放射性廢物的處置發(fā)生于1944年美國田納西州, 采用的處置方法為簡單的壕溝掩埋(Han et al., 1997)。在核廢處置的早期階段,其它核能大國如英國、前蘇聯(lián)和法國等也有類似核廢處置的報(bào)道。因此早期放射性廢物的處置并未引起重視。

    20世紀(jì)50年代后(NAS-NRC, 1957), 放射性廢物的處置開始受到廣泛重視。世界各國對(duì)放射性廢物的處置方式展開了廣泛的研究, 通過比較地表處置、海底處置、地質(zhì)處置、太空處置及其它一些處置方式, 考慮到放射性泄露的風(fēng)險(xiǎn)性、技術(shù)可行性、回收利用性及工程造價(jià)等因素, 國際上普遍認(rèn)為采用多重屏障系統(tǒng)的地質(zhì)處置是較為理想的處置方式。其中, 由于高放廢物危害大且衰變周期長, 故采用深部地質(zhì)處置的形式; 而低中放廢物由于衰變至人類可接受的時(shí)間較短, 則采用近地表或地質(zhì)處置的方式。有的國家如芬蘭則要求所有放射性廢物都采用深部地質(zhì)處置。

    低中放廢物的近地表處置是指在地表或近地表設(shè)施內(nèi)對(duì)核廢物進(jìn)行終端處置, 處置設(shè)施最大深度一般是30 m或者更淺。由于低中放廢物數(shù)量大, 處置周期短, 處置技術(shù)難度相對(duì)較小, 因此半個(gè)世紀(jì)以來發(fā)展速度較快。目前, 國際上已成功建成并運(yùn)行上百個(gè)低中放廢物處置場(Han et al., 1997; 劉平輝等, 2000)。這些處置場根據(jù)處置設(shè)施的不同, 可劃分為多種處置方式(Kittel, 1989), 并采用壓縮包裝層(水泥)、導(dǎo)水層(砂礫)、阻水層(粘土)以及生物層(植被)等作為多重屏障系統(tǒng)。

    與低中放廢物處置不同的是, 盡管關(guān)于高放廢物處置的研究較為深入, 但是真正能用于接納高放廢物的處置場則未見報(bào)道。20世紀(jì)60年代后, 歐美和一些亞洲國家建設(shè)了多座深部地下試驗(yàn)場來研究高放廢物深部地質(zhì)處置的一些關(guān)鍵問題(Kickmaier et al., 1997; Yllera et al., 2004; Wersin et al., 2004),這些地下試驗(yàn)場按功能可分為兩種, 即一般性試驗(yàn)場和特定場址試驗(yàn)場, 前者與高放廢物處置庫場址沒有直接聯(lián)系, 僅開展試驗(yàn)研究; 而后者則是在預(yù)選場址上建設(shè)的地下設(shè)施, 具有方法學(xué)研究和場址評(píng)價(jià)雙重作用。

    我國目前已建成甘肅玉門和廣東北龍兩個(gè)近地表低中放廢物處置場(浦永寧等, 2005), 并計(jì)劃在西北、西南、華南、華東、東北五個(gè)地區(qū)各建一個(gè)低中放廢物處置基地; 另一方面, 關(guān)于高放廢物的深部地質(zhì)處置, 我國也正在甘肅北山開展相關(guān)研究(Wang, 2009)。

    1.2 經(jīng)驗(yàn)認(rèn)識(shí)

    從國際上半個(gè)多世紀(jì)以來放射性廢物的處置和研究可以獲得如下四個(gè)基本認(rèn)識(shí): 1)處置過于簡單或不當(dāng)會(huì)造成嚴(yán)重后果, 早期采用的簡陋的地表淺埋處置方式已經(jīng)被淘汰, 如美國 West Vally, Sheffield和Maxey Flats處置場由于放射性泄露而關(guān)閉; 2)處置場應(yīng)盡量選在靠近產(chǎn)生廢物的地區(qū), 因?yàn)閺U物運(yùn)輸是個(gè)值得重視的事故風(fēng)險(xiǎn)問題。由于運(yùn)送貨物的卡車發(fā)生事故, 美國 Beatty和 Hanford處置場都因此曾遭到關(guān)閉; 3)采用多重屏障系統(tǒng)進(jìn)行低中放廢物的處置, 通過正確的選址、設(shè)計(jì)、建造和管理, 其安全性是可以得到保證的。如法國Manche處置場, 自1969年開始運(yùn)行到1994年關(guān)閉退役以來,至今尚未發(fā)現(xiàn)對(duì)環(huán)境和公眾構(gòu)成危險(xiǎn); 4)公眾的接受和認(rèn)可是放射性廢物處置場選址評(píng)價(jià)過程中需要特別重視的一個(gè)問題, 如我國華東低中放廢物處置場就是由于這一原因而一直沒能夠確定下來。

    核素在地質(zhì)介質(zhì)中的運(yùn)移規(guī)律是放射性廢物處置研究和環(huán)境評(píng)價(jià)的核心問題。因此, 從地質(zhì)、水文地質(zhì)及環(huán)境評(píng)價(jià)的角度看, 放射性廢物處置研究包括如下三個(gè)方面的核心課題: 1)處置場的水文地質(zhì)選址。放射性廢物的處置對(duì)場址的水文地質(zhì)條件要求較高。一般而言, 處置場址多選在人跡罕至、研究程度較低的地方, 地質(zhì)、水文地質(zhì)資料缺乏。因此, 選址過程作為獲取水文地質(zhì)資料的途徑, 在整個(gè)場址的論證、選擇、設(shè)計(jì)及后期的環(huán)境影響評(píng)價(jià)中均具有無可替代的作用; 2)核素運(yùn)移試驗(yàn)。地下水科學(xué)的進(jìn)步在很大程度上是由參數(shù)獲取的新技術(shù)和新手段帶動(dòng)起來的, 核素運(yùn)移試驗(yàn)作為研究核素運(yùn)移最直接、最有效的手段, 關(guān)系著獲取基礎(chǔ)數(shù)據(jù)的可靠程度, 同時(shí)也是場地性能評(píng)價(jià)的參數(shù)來源和重要依據(jù); 3)核素運(yùn)移模型。由于處置場地性能評(píng)價(jià)需要考慮至少數(shù)百年以上的時(shí)間尺度, 因此必須通過建立模型來量化和預(yù)測。同時(shí), 核素運(yùn)移模型也是后期環(huán)境評(píng)價(jià)的重要依據(jù)。

    2 研究現(xiàn)狀

    2.1 水文地質(zhì)選址

    國際上放射性廢物處置場的選址和建設(shè)起步較早, 并有相當(dāng)完備的法律體系和相應(yīng)的工作流程。然而, 就水文地質(zhì)選址的技術(shù)要求而言, 則僅有些原則性的規(guī)定。如美國聯(lián)邦法規(guī)10CFR61規(guī)定處置場必須在構(gòu)造相對(duì)穩(wěn)定的地區(qū), 距離地下水位要有一定距離等; 加拿大、日本、法國也只有一些與美國相似的原則性規(guī)定, 均沒有專門針對(duì)放射性廢物處置選址的操作性強(qiáng)的技術(shù)標(biāo)準(zhǔn)。同樣, 我國現(xiàn)有的技術(shù)標(biāo)準(zhǔn)(如 GB 9132-88)和核安全導(dǎo)則(如 HJ/T 23-1998)等也只有一些原則性的規(guī)定, 而對(duì)處置場選址具體的勘察、水文地質(zhì)及巖土特性參數(shù)的獲取則關(guān)注較少, 缺乏操作性, 不能有效地指導(dǎo)低中放廢物處置場地的選址。

    放射性廢物處置選址的不足集中表現(xiàn)在選址工作與環(huán)評(píng)應(yīng)用脫節(jié)方面。選址側(cè)重于場址勘察、試驗(yàn)及參數(shù)獲取; 而環(huán)評(píng)則致力于數(shù)據(jù)應(yīng)用及安全性評(píng)價(jià)。由于專業(yè)差異, 往往導(dǎo)致上述兩個(gè)方面交流較少甚至脫節(jié), 從而出現(xiàn)如下兩種不利局面: 1)選址所得數(shù)據(jù)不是環(huán)評(píng)所需要的; 2)環(huán)評(píng)所需的某些數(shù)據(jù)在野外和室內(nèi)試驗(yàn)中沒有獲取到。鑒于此, 放射性廢物處置場地的選址工作必須與環(huán)境評(píng)價(jià)結(jié)合起來, 在不同階段互相指導(dǎo), 以達(dá)到參數(shù)的有效獲取和環(huán)境影響的正確評(píng)價(jià)。

    2.2 核素運(yùn)移試驗(yàn)

    國際上核素運(yùn)移試驗(yàn)始于20世紀(jì)50年代末60年代初, 而我國在這方面的研究則起步于 80年代初。核素運(yùn)移試驗(yàn)的目的是為核廢處置場址的可行性論證、建造、運(yùn)行、封閉及長期安全性評(píng)價(jià)提供基礎(chǔ)參數(shù)。半個(gè)世紀(jì)以來, 核素運(yùn)移試驗(yàn)研究取得了顯著成果, 積累了大量的試驗(yàn)數(shù)據(jù)和技術(shù)資料(Kickmaier et al., 1997)。根據(jù)試驗(yàn)方式和規(guī)模的不同,核素運(yùn)移試驗(yàn)可分為室內(nèi)試驗(yàn)和原位試驗(yàn)兩類。兩類試驗(yàn)互有優(yōu)、缺點(diǎn), 目前都為研究者所廣泛使用,有些大型項(xiàng)目甚至兩種試驗(yàn)方法都采用, 相互印證以獲取更為可靠的試驗(yàn)結(jié)果, 或以此來研究試驗(yàn)的尺度效應(yīng)(Hadermann et al., 1996)。

    2.2.1 室內(nèi)試驗(yàn)

    在核素運(yùn)移試驗(yàn)中, 室內(nèi)試驗(yàn)是指在實(shí)驗(yàn)室內(nèi)以處置場址的主體材料作介質(zhì), 模擬核素在介質(zhì)中的吸附、滯留和擴(kuò)散等作用機(jī)制, 來研究核素在介質(zhì)中的運(yùn)移規(guī)律。室內(nèi)試驗(yàn)由于相對(duì)簡單、經(jīng)濟(jì)、易于操作, 可在較短的時(shí)間內(nèi)掌握核素運(yùn)移的一些特性常數(shù)而被大量應(yīng)用。如García-Gutiérrez等(2006)采用室內(nèi)試驗(yàn)對(duì)中性、陰離子和陽離子核素在FEBEX膨潤土中的擴(kuò)散系數(shù)和有效孔隙度等進(jìn)行了研究; Van Loon等(2009)采用室內(nèi)試驗(yàn)分別對(duì)完整和破碎的Opalinus粘土對(duì)Cs的吸附特性進(jìn)行了對(duì)比研究, 結(jié)果表明Cs在完整和破碎Opalinus粘土上的吸附作用本質(zhì)上是相同的。我國也有一些學(xué)者開展了關(guān)于核素運(yùn)移的室內(nèi)研究, 如葉明呂等(1994)用土柱試驗(yàn)法研究了90Sr在粉碎沸石柱中的吸附與運(yùn)移行為, 結(jié)果表明沸石對(duì)90Sr的吸附能力很強(qiáng);蘇銳、李春江等(2000)對(duì)核素在花崗巖基質(zhì)中的運(yùn)移進(jìn)行了研究, 獲得了134Cs在花崗巖中的擴(kuò)散系數(shù)、阻滯系數(shù)和孔隙度等; 李書紳等(2004)對(duì)237Np、238Pu、和241Am在黃土包氣帶中的運(yùn)移進(jìn)行了模擬試驗(yàn), 得到了三種核素在黃土中的垂向擴(kuò)散系數(shù)和延遲系數(shù); 王青海等(2004)用批示法對(duì)鍶在砂巖和花崗巖中的分配系數(shù)及吸附機(jī)制進(jìn)行了對(duì)比研究。另外, 我國在找礦領(lǐng)域也做了放射性核素在沉積物等礦床內(nèi)的采樣測試分析研究(高柏等, 2006; 劉廣山等, 2008)。

    室內(nèi)試驗(yàn)的另一個(gè)優(yōu)點(diǎn)是對(duì)同一試驗(yàn)可采用多種方法進(jìn)行互相印證。如在測定某種核素的擴(kuò)散參數(shù)時(shí), 可利用的方法有透入擴(kuò)散試驗(yàn)、穿過擴(kuò)散試驗(yàn)以及其它一些試驗(yàn)方法(García-Gutiérrez et al., 2006)。然而, 由于室內(nèi)試驗(yàn)所采用的樣品通常較小,而且在準(zhǔn)備過程中會(huì)受到不同程度的擾動(dòng), 因此所得的結(jié)果就會(huì)產(chǎn)生差異, 同時(shí)尺度效應(yīng)也是學(xué)術(shù)界普遍關(guān)注的一個(gè)問題, 因此室內(nèi)試驗(yàn)測試結(jié)果的適宜性受到限制和質(zhì)疑(Bradbury et al., 1985; Kickmaier et al., 1997)。

    2.2.2 原位試驗(yàn)

    為克服室內(nèi)試驗(yàn)的缺點(diǎn)和探討可能存在的尺度效應(yīng), 核素運(yùn)移試驗(yàn)逐漸被移到室外, 有針對(duì)性地進(jìn)行野外現(xiàn)場示蹤試驗(yàn)??紤]到核素的危害性和更真實(shí)地模擬核素運(yùn)移環(huán)境, 原位試驗(yàn)一般在地下研究實(shí)驗(yàn)室(URL)進(jìn)行。第一個(gè) URL建立于美國 20世紀(jì)60年代初, 歐洲第一個(gè)URL則建成于1966年,隨后如比利時(shí)的HADES、英國的RCF、瑞典的Asp?、瑞士的GTS及Mont Terri以及芬蘭、法國、德國、西班牙等國也在進(jìn)行或計(jì)劃進(jìn)行相應(yīng)的研究。根據(jù)功能的不同, 目前國際上的URLs可分為兩種, 即一般性實(shí)驗(yàn)室和特定場址實(shí)驗(yàn)室, 前者與處置場沒有直接聯(lián)系, 僅開展方法學(xué)研究; 后者則是在預(yù)選處置場上建設(shè)的地下設(shè)施, 具有方法學(xué)研究和場址評(píng)價(jià)雙重作用, 如美國的Yucca試驗(yàn)場(Karasaki et al., 2009)。

    在原位試驗(yàn)研究方面, 我國與國外的研究差距較大, 僅開展了少量野外試驗(yàn)工作。如金月如等(1994)開展了野外試驗(yàn)場包氣帶土壤中示蹤核素活度分布的直接測定研究; 中國輻射防護(hù)研究院于1997年建成了一座用于淺部含水層核素運(yùn)移規(guī)律的野外地下研究設(shè)施(URF), 該設(shè)施屬于一般性實(shí)驗(yàn)室, 僅作為方法學(xué)試驗(yàn)研究之用(李書紳等, 2002)。隨后郭擇德等(2000), 王志明等(2003)相繼在該URF內(nèi)開展了一系列核素在黃土中的現(xiàn)場示蹤試驗(yàn), 結(jié)果表明黃土沉積層對(duì)3H、60Co、85,90Sr、134Cs、237Np、238Pu及241Am等核素具有強(qiáng)烈的吸附滯留作用。

    原位試驗(yàn)通常測試兩種核素運(yùn)移濃度數(shù)據(jù): 一種是投源鉆孔中示蹤劑濃度隨時(shí)間序列的稀釋數(shù)據(jù);另一種則是在試驗(yàn)結(jié)束后, 通過切割鉆孔附近巖土樣并進(jìn)行核素濃度測試所獲得的巖芯數(shù)據(jù)。兩種數(shù)據(jù)可分別或同時(shí)運(yùn)用, 相互印證以求取較為準(zhǔn)確的核素運(yùn)移參數(shù)??梢? 原位試驗(yàn)很好地克服了室內(nèi)試驗(yàn)的不足, 所測得數(shù)據(jù)更具有代表性, 因此受到各國核素運(yùn)移研究的重視。然而, 原位試驗(yàn)也存在耗資大, 歷時(shí)長、操作復(fù)雜等問題, 同時(shí)也具有一些無法克服的不理想因素, 如過濾器、鉆孔與圍巖間的空隙及干擾破碎帶等, 均會(huì)影響核素運(yùn)移結(jié)果,為數(shù)據(jù)解釋帶來較大的不確定性, 因此要求研究者必須具有良好的判斷能力和豐富的研究手段(Samper et al., 2010)。

    核素運(yùn)移試驗(yàn)示蹤劑的應(yīng)用經(jīng)歷了一個(gè)由簡單到復(fù)雜的過程。由早期的非吸附性示蹤劑(中性及陰離子如 HTO、HDO、82Br-、123I-、3He等)逐步到各種弱吸附性示蹤劑(如22,24Na+、85Sr2+、133Ba2+、86Rb等)以及需要較長試驗(yàn)時(shí)間的強(qiáng)吸附示蹤劑(如134,137Cs+等)演化。近兩年, 化學(xué)成分復(fù)雜的示蹤劑也被逐漸利用起來, 如234,235,238U+,60Co+,237Np,152Eu,99Tc,113Sn,75Se等(Frick et al., 1992; Smith et al., 2001)。

    需要說明的是, 一種試驗(yàn)方法的優(yōu)缺點(diǎn)是在試驗(yàn)過程中逐步發(fā)現(xiàn)的, 因此, 各國學(xué)者均在針對(duì)目前試驗(yàn)方法的不足, 進(jìn)一步完善試驗(yàn)技術(shù), 目前已成為核素試驗(yàn)研究的重要內(nèi)容。如 García-Gutiérrez等(2006)提出了一種大尺度固源擴(kuò)散試驗(yàn)方法, 同時(shí)考慮核素運(yùn)移的原位性和克服目前原位試驗(yàn)中存在的過濾器、空隙和破碎帶等引起的不確定性, 以獲取更準(zhǔn)確的Cox粘土的擴(kuò)散和阻滯參數(shù)。另外, 天然類比研究(Mckinley, 1987)以及直接針對(duì)核事故的核素運(yùn)移調(diào)查研究也是重要的研究方法(Rudenko et al., 2003)。

    2.3 核素運(yùn)移模型

    核廢處置場地安全性評(píng)價(jià)所涉及的時(shí)間尺度少則數(shù)百年, 多則上萬甚至百萬年, 遠(yuǎn)遠(yuǎn)超出直接測試的試驗(yàn)范圍, 因此需要基于核素在地下介質(zhì)中的運(yùn)移機(jī)理和過程, 建立相應(yīng)的模型來定量預(yù)測核素在處置場內(nèi)的運(yùn)移。核素運(yùn)移模型研究的基本意義包括三個(gè)方面: 1)增強(qiáng)對(duì)核素運(yùn)移的理解; 2)印證研究方法的適宜性及評(píng)價(jià)試驗(yàn)數(shù)據(jù)的有效性; 3)提高模擬技術(shù)以進(jìn)行核素地質(zhì)處置的安全性評(píng)價(jià)(Hadermann et al., 1996)。

    核素在地下介質(zhì)中的運(yùn)移包括水動(dòng)力運(yùn)移過程和地球化學(xué)反應(yīng)過程兩個(gè)方面。水動(dòng)力運(yùn)移過程主要包括對(duì)流、機(jī)械彌散及分子擴(kuò)散等作用, 可用對(duì)流-彌散模型(Bear, 1961)來描述; 而地球化學(xué)反應(yīng)過程則包括水相絡(luò)合、表面絡(luò)合、氧化-還原作用、酸堿反應(yīng)、離子交換、吸附作用、沉淀-溶解作用等(Gillham et al., 1982), 可采用地球化學(xué)模型來刻畫。20世紀(jì)80年代以前, 研究者們通常把水動(dòng)力運(yùn)移過程和地球化學(xué)反應(yīng)過程作為獨(dú)立的系統(tǒng)進(jìn)行模擬。

    然而, 地球化學(xué)過程和水動(dòng)力過程在實(shí)際情況下是相互影響的。由于水動(dòng)力過程在數(shù)學(xué)原理上能簡化為代數(shù)系統(tǒng), 20世紀(jì)80年代后Kd模型通常被用來模擬核素的地球化學(xué)吸附行為。這種方法將復(fù)雜的化學(xué)系統(tǒng)簡化為由溶質(zhì)在固、液相中的分布系數(shù)(Kd)來控制。盡管這種方法不屬于確定性動(dòng)力模型,無法揭示發(fā)生在系統(tǒng)內(nèi)部的化學(xué)反應(yīng)機(jī)制, 但其作為與化學(xué)反應(yīng)有關(guān)的運(yùn)移模型, 在缺乏詳盡基礎(chǔ)數(shù)據(jù)的時(shí)候卻非常有效, 因此目前Kd模型仍被廣泛研究和應(yīng)用(Van Der Lee et al., 2001)。但是, Kd模型畢竟是一種對(duì)化學(xué)過程的概化。在真實(shí)地下水系統(tǒng)中,一方面由于溶質(zhì)與圍巖礦物的動(dòng)態(tài)反應(yīng), 水化學(xué)組分總是隨時(shí)間和空間的不同而改變, 因此也就受水動(dòng)力系統(tǒng)的影響; 另一方面, 介質(zhì)的孔隙度與滲透性也受礦物沉淀和溶解的影響。因此, Kd不是恒定的而且在有的情況下變化比較大。已有研究表明在溶液成分濃度受礦物溶解-沉淀影響強(qiáng)烈(Reardon, 1981; Liu et al., 1989)和水化學(xué)作用強(qiáng)烈的地區(qū)(Yeh et al., 1991), 運(yùn)用Kd法是不合適的。因此, 就需要一種耦合技術(shù)來模擬同時(shí)產(chǎn)生的地球化學(xué)和水動(dòng)力運(yùn)移過程, 溶質(zhì)反應(yīng)性運(yùn)移模型就是基于這個(gè)目的而發(fā)展起來的(Grindrod et al., 1996; Darban et al., 2008)。

    直到90年代初, 反應(yīng)性運(yùn)移模型的研究進(jìn)展仍然很小, 通常只是對(duì)其進(jìn)行介紹并呼吁學(xué)術(shù)界進(jìn)行研究。然而20年來, 溶質(zhì)反應(yīng)性運(yùn)移模型的研究取得了巨大的進(jìn)步, 這可從近年來逐步增加的文獻(xiàn)發(fā)表量上顯示出來, 在此不對(duì)文獻(xiàn)做詳細(xì)評(píng)論, 可參閱如下論文及相關(guān)文獻(xiàn)(Wersin et al., 2004; García-Gutiérrez et al., 2006; Van Der Lee et al., 2001; Xu et al., 1999; Zheng et al., 1999; Prommer et al., 2003; Palut et al., 2003; Samper et al., 2006; Appelo et al., 2007; Ma et al., 2010)。

    結(jié)合核素運(yùn)移試驗(yàn)和核廢處置場地實(shí)際案例,國際上研究開發(fā)了為數(shù)眾多的數(shù)值模擬工具。目前,反應(yīng)性運(yùn)移模型已經(jīng)能夠正確用于模擬動(dòng)力試驗(yàn)及野外觀測, 并可以同時(shí)考慮多種地球化學(xué)過程及多種化學(xué)組分。常用的核素運(yùn)移計(jì)算程序包括HYTEC (Van Der Lee, 1997)、3DHYDROGEOCHEM(Cheng et al., 1998)、TOUGHREACT(Xu et al., 1999)、MT3D系列程序(Zheng et al., 1999)、PHT3D(Prommer et al., 2003)及其余一系列程序, 具體可參考 Yim et al.(2000)、Van Der Lee et al.(2001)、Zheng et al.(2002)等文獻(xiàn), 里面統(tǒng)計(jì)了當(dāng)前常用的數(shù)十種計(jì)算程序及其應(yīng)用條件和功能特征等。

    我國學(xué)者在核素運(yùn)移模型方面也開展了部分工作。如王榕樹等(1994)研究了表觀擴(kuò)散系數(shù)與地下水流速等因素對(duì)核素?cái)U(kuò)散的影響; 陳家軍等(1994)、王金生等(1996)研究了水動(dòng)力彌散方程與核素在孔隙介質(zhì)中非平衡吸附模式的耦合模型, 結(jié)果表明二維非平衡吸附數(shù)值模式能較好地描述核素運(yùn)移實(shí)際;李春江等(1998)、蘇銳等(2000)、李金軒(2001)對(duì)花崗巖體單裂隙中的核素運(yùn)移模型進(jìn)行了研究, 分析了水流變化對(duì)核素運(yùn)移的影響及控制核素運(yùn)移的主要參數(shù); 另外, 其他學(xué)者如錢天偉等(2002)、王青海等(2004)也開展了部分核素在黃土及基巖裂隙中的模擬研究工作。然而, 應(yīng)該認(rèn)識(shí)到我國關(guān)于核素運(yùn)移模型的研究與國際上存在較大差距, 具體表現(xiàn)在三個(gè)方面: 1)對(duì)核素運(yùn)移計(jì)算程序的開發(fā)差距明顯; 2)實(shí)際應(yīng)用實(shí)例的數(shù)量和質(zhì)量遠(yuǎn)遠(yuǎn)落后于歐美發(fā)達(dá)國家; 3)高水平研究工作較少, 高水平和高質(zhì)量的文章與發(fā)達(dá)國家差距較大。上述差距的存在一方面是由于我國在核素運(yùn)移模型方面的研究起步較晚, 一開始就沒有跟上國際發(fā)展步伐; 另一方面則主要在于重視程度不夠, 投入力度不大等因素。

    綜合國內(nèi)外研究狀況, 模型研究極大地強(qiáng)化了對(duì)核素在地質(zhì)介質(zhì)中的運(yùn)移和阻滯過程的理解和掌握。然而, 對(duì)概念模型的理解不足以及數(shù)據(jù)的缺乏,限制了模型的應(yīng)用(Van Der Lee et al., 2001; Berner, 1998; Zheng, 2009)。核素運(yùn)移模型研究存在如下不足: 1)由于基本數(shù)據(jù)的欠缺, Kd法作為與化學(xué)反應(yīng)有關(guān)的模型之一, 仍在實(shí)際模擬中廣為運(yùn)用, 但其穩(wěn)定性和有效性仍然需要進(jìn)一步的研究和論證; 2)耦合水動(dòng)力運(yùn)移和地球化學(xué)反應(yīng)的多組分反應(yīng)性運(yùn)移模型取得了巨大發(fā)展, 然而由于缺乏對(duì)熱力學(xué)動(dòng)態(tài)或平衡常數(shù)的認(rèn)識(shí), 常常對(duì)模型的應(yīng)用形成約束。特別是我國, 關(guān)于核素運(yùn)移的反應(yīng)性模擬實(shí)例還未見報(bào)道; 3)模型反向參數(shù)估算取得了巨大進(jìn)展,目前可對(duì)溶液組分、核素運(yùn)移固、液相參數(shù)等進(jìn)行反推預(yù)算, 解決了由于技術(shù)局限而無法測試組分濃度的困難, 也為某些條件下無法采用解析法來解釋試驗(yàn)數(shù)據(jù)時(shí)提供了備選方法。但是, 由于模型參數(shù)反向估算計(jì)算量大, 同時(shí)存在計(jì)算不收斂等問題,計(jì)算效率和精度等還需要進(jìn)一步提高, 計(jì)算方法也有待于進(jìn)一步改進(jìn); 4)含水介質(zhì)的地球化學(xué)組成在不同尺度上均存在非均質(zhì)性, 這就要求模型應(yīng)當(dāng)具有處理非均質(zhì)屬性的能力, 因此模型研究應(yīng)在刻畫介質(zhì)和反應(yīng)的非均質(zhì)性方面作出進(jìn)一步的突破; 5)礦物沉淀和溶解對(duì)水動(dòng)力行為的影響仍然很難予以量化, 而這一作用又對(duì)核素在含水介質(zhì)中的運(yùn)移具有重要影響, 因此模型研究還有待在這個(gè)方面予以加強(qiáng)。

    3 結(jié)論及討論

    放射性廢物的處置從地質(zhì)、水文地質(zhì)及環(huán)境評(píng)價(jià)的角度看, 主要包括水文地質(zhì)選址、核素運(yùn)移試驗(yàn)和核素運(yùn)移模型三個(gè)方面的核心課題。場址的選擇應(yīng)本著就近原則并采用多重屏障系統(tǒng)進(jìn)行廢物處置, 同時(shí)應(yīng)基于環(huán)境水文地質(zhì)的角度來構(gòu)建場選址評(píng)價(jià)體系。核素運(yùn)移試驗(yàn)主要包括室內(nèi)試驗(yàn)和原位試驗(yàn)兩種, 試驗(yàn)示蹤劑的運(yùn)用呈現(xiàn)出多元化、復(fù)雜化的特征, 并催生出試驗(yàn)方法學(xué)這一研究熱點(diǎn)。核素運(yùn)移模型研究經(jīng)歷了由對(duì)流-彌散模型→Kd模型→多組分反應(yīng)性模型的發(fā)展過程, 并取得了巨大進(jìn)展。但是由于對(duì)核素運(yùn)移機(jī)理、概念認(rèn)識(shí)不足以及熱、動(dòng)力學(xué)關(guān)鍵參數(shù)的缺乏限制了核素運(yùn)移模型的應(yīng)用。模擬模型還應(yīng)在處理非均質(zhì)介質(zhì)和非均質(zhì)反應(yīng)、描述礦物溶解-沉淀與水動(dòng)力行為的相互影響以及穩(wěn)定性和實(shí)際應(yīng)用能力等方面做出突破和提高。

    我國放射性廢物處置及研究還處于初步階段,處置場址勘測體系和標(biāo)準(zhǔn)仍待建立, 核素運(yùn)移試驗(yàn)技術(shù)需要完善, 核素運(yùn)移數(shù)據(jù)儲(chǔ)備急待加強(qiáng); 在核素運(yùn)移模型方面, 我國在計(jì)算程序的開發(fā)及模擬實(shí)例上與國際上的差距則更為明顯, 滿足不了目前因大力發(fā)展核電而導(dǎo)致的核廢處置需求。因此, 不論從理論研究意義還是從放射性廢物處置的實(shí)際需要來看, 我國均需要在上述幾個(gè)方面加強(qiáng)研究。

    陳家軍, 戴杰, 郭擇德, 趙英杰.1994.非平衡吸附放射性核素遷移數(shù)值模型研究[J].水文地質(zhì)工程地質(zhì), 21(2): 37-40.

    高柏, 孫占學(xué), 王廣才, 張文, 劉金輝, 吳文珍.2006.砂巖鈾礦體定位條件的鈾、釷含量及其放射性同位素研究[J].地球?qū)W報(bào), 27(6): 590-594.

    郭擇德, 王志明, 李書紳, 姜鋼, 黃慶春, 楊端節(jié), 田中忠夫, 向井雅之, 神山秀雄.2000.黃土包氣帶中放射性核素遷移的現(xiàn)場試驗(yàn)[J].輻射防護(hù), 20(1-2): 21-30.

    金月如, 張秀珍, 馬如維, 樊耀國, 毛永, 夏曉彬, 韓月琴, 郭彩萍.1994.土壤中核素遷移實(shí)驗(yàn)示蹤核素活度分布的直接測定[J].輻射防護(hù)通訊, 14(4): 87, 166.

    李春江, 楊天行.1998.花崗巖體單裂隙中核素遷移數(shù)學(xué)模型:Ⅰ.連續(xù)注入實(shí)驗(yàn)解析模型[J].核化學(xué)與放射化學(xué), 20(4): 221-227.

    李金軒.2001.非恒定流條件下單裂隙核素遷移模型[J].勘察科學(xué)技術(shù), (6): 7-10.

    李書紳, 倪世偉, 郭擇德, 程萍, 王志明, 韓新生, 姚來根, 蔡勻.2004.237Np、238Pu和241Am在包氣帶黃土中遷移的模擬實(shí)驗(yàn)[J].輻射防護(hù), 24(6): 347-355.

    李書紳, 王志明, 趙英杰, 李禎堂, 郭擇德, 郭亮天, 許兆義, 楊成永.2002.建造用于含水層中核素遷移試驗(yàn)的地下研究設(shè)施的技術(shù)要求[J].輻射防護(hù), 22(4): 200-206.

    劉廣山, 李冬梅, 易勇, 劉素美, 白潔, 張經(jīng).2008.膠州灣沉積物的放射性核素含量分布與沉積速率[J].地球?qū)W報(bào), 29(6): 769-777.

    劉平輝, 管太陽, 王勇, 黃國夫, 余運(yùn)祥.2000.陸地中低放核廢物地質(zhì)處置的發(fā)展與現(xiàn)狀[J].華東地質(zhì)學(xué)院學(xué)報(bào), 23(3): 229-234.

    浦永寧, 范選林, 韓國勝, 劉艷, 歐紅.2005.對(duì)我國低中放固體廢物處置的若干問題的探討[C].北京: 中國學(xué)術(shù)期刊電子雜志出版社.

    錢天偉, 李書紳, 武貴賓.2002.地下水多組分反應(yīng)溶質(zhì)遷移模型的研究進(jìn)展[J].水科學(xué)進(jìn)展, 13(1): 116-121.

    蘇銳, 李春江, 王駒, 高宏成.2000.花崗巖體單裂隙中核素遷移數(shù)學(xué)模型 III.擴(kuò)散模型及其有限單元法解[J].核化學(xué)與放射化學(xué), 22(2): 80-86.

    王金生, 李書紳, 王志明, 郭擇德, 陳家軍.1996.核素遷移的二維非平衡吸附數(shù)值模式研究[J].應(yīng)用生態(tài)學(xué)報(bào), 7(1): 89-93.

    王青海, 王蘭生, 李曉紅.2004.基巖裂隙水中 90Sr遷移的數(shù)值模擬[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 34(3): 405-409.

    王榕樹, 馮為.1994.放射性核素在地質(zhì)介質(zhì)中的遷移研究[J].核化學(xué)與放射化學(xué), 16(2): 117-121.

    王志明, 李書紳, 劉春立, 姜凌, 楊月娥, 江洪, 李冰, 李丹.2003.核素在弱含水層三維遷移的現(xiàn)場示蹤試驗(yàn)[J].原子能科學(xué)技術(shù), 37(2): 149-153.

    葉明呂, 陸誓俊, 秦春扣, 許麗紅, 何阿弟, 唐志華, 徐國慶,范選林, 顧綺芳, 杜志超.1994.放射性鍶在沸石中的吸附與遷移研究[J].核化學(xué)與放射化學(xué), 16(4): 199-204.

    References:

    APPELO C A J, WERSIN P.2007.Multicomponent diffusion modeling in clay systems with application to the diffusion of Tritium, Iodide, and Sodium in Opalinus Clay[J].Environment Science Technology, 41(14): 5002-5007.

    BEAR J.1961.On the tensor form of dispersion in porous media[J].J.Geophys.Res., 66(4): 1185-1197.

    BERNER U.1998.Geochemical modelling of repository systems: limitations of the thermodynamic approach[J].Radiochim.Acta, 82: 423-428.

    BRADBURY M H, STEPHEN I G.1985.Scientific Basis for Nuclear Waste Management IX[Z].Pittsburg: Material Research Society: 81.

    CHEN Jia-jun, DAI Jie, GUO Ze-de, ZHAO Ying-jie.1994.Study on the numerical models on the nonequilibrium sorption of radionuclide[J].Hydrogeology and Engineering Geology, 21(2): 37-40(in Chinese with English abstract).

    CHENG H P, YEH G T.1998.Development and demonstrative application of a 3-D numerical model of subsurface flow, heat transfer, and reactive chemical transport: 3D hydrogeochem[J].J.Contam.Hydrol., 34(1-2): 47-83.

    DARBAN A K, YONG R N, RAVAJ S.2008.Coupled chemical speciation-solute transport model for prediction of solute transport in clay buffers[J].Applied Clay Science,doi:10.1016/j.clay.2008.11.002.

    FRICK U, ALEXANDER W R, BAEYENS B, BOSSART P, BRADBURY M H, BüHLER C, EIKENBERG J, FIERZ T, HEER W, H?HN E, Mckinley I G, SMITH P A.1992.Grimsel Test Site: The radionuclide migration experiment—overview of investigations 1985-1990[R].Nagra Technical Report Series NTB 91-04.Nagra, Wettingen, Switzerland.

    GAO Bai, SUN Zhan-xue, WANG Guang-cai, ZHANG Wen, LIU Jin-hui, WU Wen-zhen.2006.The Delineation of Ore Bodies in Sandstone Type Uranium Deposits by U and Th Concentrations and Their Radioisotopes[J].Acta Geoscientica Sinica, 27(6): 590-594(in Chinese with English abstract).

    GARCíA-GUTIéRREZ M, CORMENZANA J L, MISSANA T, MINGARRO M, MARTíNA P L.2006.Large-scale laboratory diffusion experiments in clay rocks[J].Physics and Chemistry of the Earth, 31(10-14): 523-530.

    GILLHAM R W, CHERRY J A.1982.Contaminant Migration in Saturated Unconsolidated Geologic Deposits[Z].Geological Society of American Special Paper, 189: 31-61.

    GRINDROD P, TAKASE H.1996.Reactive chemical transport within engineered barriers[J].Journal of Contaminant Hydrology, 21(1-4): 283-296.

    GUO Ze-de, WANG Zhi-ming, LI Shu-shen, JIANG Gang, HUANG Qing-chun, YANG Duan-jie, TANAKA T, MUKAI M, KAMIYAMA H.2000.Field test of radionuclide migration in loess aerated zone[J].Radiation Protection, 20(1-2): 21-30( in Chinese with English abstract).

    HADERMANN J, HEER W.1996.The Grimsel (Switzerland) migration experiment: integrating field experiments, laboratory investigations and modelling[J].Journal of Contaminant Hydrology, 21(1-4): 87-100.

    HAN K W, HEINONEN J, BONNE A.1997.Radioactive Waste Disposal: Global Experience and Challenges[R].IAEA Bulletin, 39(1): 33.

    JIN Yue-ru, ZHANG Xiu-zhen, MA Ru-wei, FAN Yao-guo, MAO Yong, XIA Xiao-bin, HAN Yue-qin, GUO Cai-ping.1994.Direct measurement of radionuclide activities in migration experiments in soils[J].Radiation Protection Bullitin, 14(4): 87, 166(in Chinese with English abstract).

    KARASAKI K, APPS J, DOUGHTY C.2009.Feature Detection, Characterization and Confirmation Methodology: Final Report[R].Lawrence Berkeley National Laboratory: LBNL Paper LBNL-1358E.Retrieved from: http://escholarship.org/uc/item/7wp4f7jg.

    KICKMAIER W, MCKINLEY I.1997.A review of research carried out in European rock laboratories[J].Nuclear Engineeering and Design, 176(1-2): 75-81.

    KITTEL J H.1989.NEAR-SRUFACE LAND DISPOSAL, Radioactive Waste Management Handbook (volume 1)[M].Chur·London·Paris·New York·Melbourne: Harwood Academic Publishers.

    LI Chun-jiang, YANG Tian-xing.1998.A Mathematical Model for Radionuclide Migration in Single Fractured Granite: Ⅰ.Analytical Solution for A Continuous Injection Experiment[J].Journal of Nuclear and Radiochemistry, 20(4): 221-227(in Chinese with English abstract).

    LI Jin-xuan.2001.Single Fracture Nuclide Transport Model Under Time-Varying Groundwater Conditions[J].Site Investigation Science and Technology, (6): 7-10(in Chinese with English abstract).

    LI Shu-shen, NI Shi-wei, GUO Ze-de, CHENG Ping, WANG Zhi-ming, HAN Xin-sheng, YAO Lai-gen, CAI Yun.2004.Simulation experiments for migration of237Np、238Pu and241Am in aerate zone loess[J].Radiation Protection, 24(6): 347-355(in Chinese with English abstract).

    LI Shu-shen, WANG Zhi-ming, ZHAO Ying-jie, LI Zhen-tang, GUO Ze-de, GUO Liang-tian, XU Zhao-yi, YANG Cheng-yong.2002.Technical requirements for construction of underground research facility used for nuclide migration test in aqiufer[J].Radiation Protection, 22(4): 200-206(in Chinese with English abstract).

    LIU C W, NARASIMHAN T N.1989.Redox-controlled multiplespecies reactive chemical transport: 2.Verification and application[J].Water Resour.Res., 25(5): 883-910.

    LIU Guang-shan, LI Dong-mei, YI Yong, LIU Su-mei, BAI Jie,ZHANG Jin.2008.Radionuclide distribution in sediments and and sedimentary rates in the Jiaozhou Bay[J].Acta Geoscientica Sinica, 29(6): 769-777(in Chinese with English abstract).

    LIU Ping-hui, GUAN Tai-yang, WANG Yong, HUANG Guo-fu, YU Yun-xiang.2000.Development and Current Situation on Geological Disposal of Low and Intermediate-Level Waste on Land[J].Journal of East China Geological Institute, 23(3): 229-234(in Chinese with English abstract).

    MA R, ZHENG C, PROMMER H, GRESKOWIAK J, LIU C, ZACHARA J, ROCKHOLD M.2010.A field-scale reactive transport model for U(VI) migration influenced by coupled multirate mass transfer and surface complexation reactions[J].Water Resour.Res., 46: W05509,doi:10.1029/2009WR 008168.

    MCKINLEY I G.1987.Int Conf Chem and Migration Behavior of Ac and Fish[Z].Munich: Prod in Geosphere, Abst No: 129/1.

    NAS-NRC.1957.The Disposal of Radioactive Waste on Land[R].Report of the Committee on Waste Disposal of the Division of Earth Sciences, Apr.1957.

    PALUT J M, MONTARNAL P, GAUTSCHI A, TEVISSEN E, MOUCHE E.2003.Characterisation of HTO diffusion properties by an in situ tracer experiment in Opalinus clay at Mont Terri[J].Journal of Contaminant Hydrology, 61(1-4): 203-218.

    PROMMER H, BARRY D A, ZHENG C.2003.MODFLOW/ MT3DMS-Based Reactive Multicomponent Transport Modeling[J].Ground Water, 41(2): 247-257.

    PU Yong-Ning, FAN Xuan-lin, HAN Guo-sheng, LIU Yan, OU Hong.2005.Discussion on the disposal of LILW in China[C].Beijing: China Academic Journal Electronic Publishing House(in Chinese).

    QIAN Tian-wei, LI Shu-shen, WU Gui-bin.2002.Advances in study on transport models of multicomponent reacting solutes in groundwater[J].Advances in Water Science, 13(1): 116-121(in Chinese with English abstract).

    REARDON E J.1981.Kd’s – Can they be used to describe reversible ion sorption reaction in contaminant migration?[J].Ground Water, 19(3): 279-286.

    RUDENKO L I, KHAN V E, PANASYUK N I.2003.Physicochemical Study of the Mechanism of Radionuclide Migration from the Shelter Objecter and Its Service Area to Groundwater[J].Radiochemistry, 45(3): 293-297.

    SAMPER J, YANG C, NAVES A, YLLERA A, HERNANDEZ A, MOLINERO J, SOLER J M, HERNAN P, MAYOR J C, ASTUDILLO J.2006.A fully 3-D anisotropic numerical model of the DI-B in situ diffusion experiment in the Opalinus clay formation[J].Physics and Chemistry of the Earth, Parts A/B/C, 31(10-14): 531-540.

    SAMPER J, YI S, NAVES A.2010.Analysis of the parameter identifiability of the in situ diffusion and retention (DR) experiments[J].Physics and Chemistry of the Earth, Parts A/B/C, 35(6-8): 207-216.

    SMITH P A, ALEXANDER W R, KICKMAIER W, OTA K, FRIEG B, MCKINLEY I G.2001.Development and testing of radionuclide transport models for fractured rock: examples from the Nagra/JNC Radionuclide Migration Programme in the Grimsel Test Site, Switzerland[J].Journal of Contaminant Hydrology, 47(2-4): 335-348.

    SU Rui, LI Chun-jiang, WANG Ju, GAO Hong-cheng.2000.A Mathematical Model for Radionuclide Transport in Single Fractured Granite III.Diffusion Model and Solution by Finite Element Method[J].Journal of Nuclear and Radiochemistry, 22(2): 80-86(in Chinese with English abstract).

    VAN DER LEE J, DE WINDT L.2001.Present state and future directions of modelling of geochemistry in hydrogeological systems[J].Journal of Contaminant Hydrology, 47(2-4): 265-282.

    VAN DER LEE J.1997.HYTEC, un modèle couplé hydro-géochimique de migration de polluants et de collo?des[R].Technical Report LHMrRDr97r02.CIG, école des Mines de Paris, Fontainebleau, France.

    VAN LOON L R, BAEYENS B, BRADBURY M H.2009.The sorption behaviour of caesium on Opalinus Clay: A comparison between intact and crushed material[J].Applied Geochemistry,doi:10.1016/j.apgeochem.2009.03.003.

    WANG Jin-sheng, LI Shu-shen, WANG Zhi-ming, GUO Ze-de, CHEN Jia-jun.1996.Study on the 2D numerical model for nonequilibrium sorption[J].Chinese Journal of Applied Ecology, 7(1): 89-93(in Chinese with English abstract).

    WANG Ju.2009.The disposal of HLW in China[Z].Beijing: Beijing Research Institute of Uranium Geology.

    WANG Qing-hai, WANG Lan-sheng, LI Xiao-hong.2004.Numerical simulation of the 90Sr migration in fissure water[J].Journal of Jilin University(Earth Science Edition), 34(3): 405-409(in Chinese).

    WANG Rong-shu, FENG Wei.1994.The Study on Radionuclide Migration in Geologic Media[J].Journal of Nuclear and Radiochemistry, 16(2): 117-121(in Chinese with English abstract).

    WANG Zhi-ming, LI Shu-shen, LIU Chun-li, JIANG Ling, YANG Yue-e, JIANG Hong, LI Bing, LI Dan.2003.Three-dimension Migration Test of Radio-nuclides in Weak Phreatic Aquifer[J].Atomic Energy Scicence and Technology, 37(2): 149-153(in Chinese with English abstract).

    WERSIN P, VAN LOON L R, SOLER J M, YLLERA A, EIKENBERG J, GIMMI TH, HERNáN P, BOISSON J Y.2004.Long-term Diffusion Experiment at Mont Terri: First Results from Field and Laboratory Data[J].Appl.Clay Sci., 26(1-4): 123-135.

    XU T, PRUESS K, BRIMHALL G.1999.An improved equilibrium-kinetics speciation algorithm for redox reactions in variably saturated subsurface flow systems[J].Comput.Geosci., 25(6): 655-666.

    YE Ming-lü, LU Shi-jun, QIN Chun-kou, XU Li-hong, HE A-di, TANG Zhi-hua, XU Guo-qing, FAN Xuan-lin, GU Qi-fang, DU Zhi-chao.1994.Investigation of Sorption and Migration of Radioactive Strontium on Clinoptilolite and Mordenite[J].Chinese Journal of Nuclear Science and Engineering, 16(4): 199-204(in Chinese with English abstract).

    YEH G T, TRIPATHI V S.1991.A model for simulating transport of reactive multispecies components: model development and demonstration[J].Water Resour.Res., 27(12): 3075-3094.

    YIM M S, SIMONSON S A.2000.Performance Assessment Models for Low Level Radioactive Waste Disposal Facilities: A Review[J].Progress in Nuclear Energy, 36(1): 1-38.

    YLLERA A, HERNáNDEZ A, MINGARRO M, QUEJIDO A, SEDANO L A, SOLER J M, SAMPER J, MOLINERO J, BARCALA J M, MARTíN P L, FERNáNDEZ M, WERSIN P, RIVAS P, HERNáN P.2004.DI-B experiment: planning, design and performance of an in situ diffusion experiment in the Opalinus Clay formation[J].Applied Clay Sci., 26(1-4): 181-196.

    ZHENG C, BENNETT G D.2002.Applied Contaminant Transport Modeling (Second Edition)[M].New York: John Wiley & Sons, Inc.

    ZHENG C, WANG P.1999.PMT3DMS: A Modular Three-Dimensional Multi-Species Model for Simulation of Advection.Dispersion and Chemical Reactions of Contaminants in Groundwater Systems: Documentation and User’s Guide[Z].SERDP-99-1, U.S.Army Engineer Research and Development Center, Vicksburg, MS.

    ZHENG C.2009.Recent developments and future directions for MT3DMS and Related Transport Codes[J].Ground Water, 47(5): 620-625.

    Advances in Research on Disposal of Radioactive Waste

    YI Shu-ping1,2), MA Hai-yi1), ZHENG Chun-miao2)
    1) Guangdong Electric Power Design Institute, Guangzhou, Guangdong 510600;
    2) PKU Center for Water Research, Peking University, Beijing 100871

    Disposal of radioactive wastes (DRW) will eventually hinder the sustainable development of nuclear energy.Safety DRW has become an issue that receives worldwide attentions.Aimed at tackling the exposure risk of radioactive waste transport in the groundwater system of disposal sites, this paper presents a review on DRW along with site selection, radionuclide transport experiments and numerical models.Safety of DRW can be achieved by disposal of radioactive wastes with a multi-barrier system.A disposal site should be located at an area near the place where the wastes are produced.The site investigation and assessment system should focus on aspects of environmental hydrogeology.Both lab and in-situ test should be considered for obtaining correct results and parameters for the transport of radionuclides in a groundwater system.Numerical models for radionuclides transport have made great progress through advection-dispersion and linear sorption (Kd) to multi-component reactive transport models.However, the numerical methods need improvement in dealing with heterogeneous media and chemical reactions, interactions between mineral dissolution/precipitation and hydrodynamic process, model validity as well as the capability of application to real cases.DRW in China has been relatively poorly studied and work needs to be improved in such aspects as site selection standard series, experimental techniques and modeling approaches.

    disposal of radioactive waste; hydrogeology; assessment of environment affects; radionuclide transport experiments; radionuclide transport models

    X753; X820.3; P641

    A

    10.3975/cagsb.2011.05.09

    本文由國家自然科學(xué)基金項(xiàng)目(編號(hào): 40911130505)資助。

    2011-07-03; 改回日期: 2011-08-06。責(zé)任編輯: 閆立娟。

    易樹平, 男, 1979年生。博士。主要從事地下水循環(huán)和溶質(zhì)運(yùn)移研究。E-mail: ysp3145@163.com。

    猜你喜歡
    核素運(yùn)移廢物
    核素分類開始部分的6種7核素小片分布
    深度睡眠或有助于有效排出廢物
    中老年保健(2021年4期)2021-12-01 11:19:40
    廢物不廢的西瓜皮
    曲流河復(fù)合點(diǎn)壩砂體構(gòu)型表征及流體運(yùn)移機(jī)理
    核素分類的4量子數(shù)
    東營凹陷北帶中淺層油氣運(yùn)移通道組合類型及成藏作用
    電子廢物何處去
    開采過程中上覆急傾斜巖層運(yùn)移規(guī)律模擬與研究
    川西坳陷孝泉-新場地區(qū)陸相天然氣地球化學(xué)及運(yùn)移特征
    廢物巧利用
    熟女电影av网| 亚洲欧美日韩高清专用| 黄色丝袜av网址大全| 亚洲国产欧美网| 在线天堂最新版资源| 最新中文字幕久久久久| 国产精品久久久久久人妻精品电影| 无人区码免费观看不卡| 看片在线看免费视频| 国产极品精品免费视频能看的| 美女免费视频网站| 亚洲国产精品成人综合色| 最新在线观看一区二区三区| 在线观看舔阴道视频| 日韩人妻高清精品专区| 亚洲无线在线观看| 在线观看免费午夜福利视频| 午夜视频国产福利| 婷婷精品国产亚洲av在线| 变态另类丝袜制服| 国产成人aa在线观看| 久久久久久久亚洲中文字幕 | 在线观看免费视频日本深夜| 99久久精品热视频| 桃红色精品国产亚洲av| 亚洲人成网站在线播| 老司机深夜福利视频在线观看| 国产精品影院久久| 亚洲av电影在线进入| 久99久视频精品免费| 亚洲在线观看片| 亚洲国产高清在线一区二区三| 99久久精品热视频| 欧美xxxx黑人xx丫x性爽| 欧美另类亚洲清纯唯美| av福利片在线观看| 欧美黄色片欧美黄色片| 三级毛片av免费| 美女被艹到高潮喷水动态| 国内精品美女久久久久久| 美女黄网站色视频| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 国产精品综合久久久久久久免费| 精品久久久久久久毛片微露脸| 最后的刺客免费高清国语| 国产日本99.免费观看| 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 岛国视频午夜一区免费看| 嫁个100分男人电影在线观看| 尤物成人国产欧美一区二区三区| 精品一区二区三区视频在线观看免费| 女生性感内裤真人,穿戴方法视频| 久久久国产成人免费| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 他把我摸到了高潮在线观看| 久久亚洲真实| 亚洲人成网站高清观看| 亚洲专区中文字幕在线| 国产一区二区在线av高清观看| 欧美色视频一区免费| 精品久久久久久久久久免费视频| 亚洲成人免费电影在线观看| АⅤ资源中文在线天堂| 亚洲成人久久爱视频| 国产av麻豆久久久久久久| 久久久久久人人人人人| 亚洲av一区综合| 成年免费大片在线观看| 成年版毛片免费区| 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| 久久久精品大字幕| 亚洲精品在线观看二区| 日本免费a在线| 国产单亲对白刺激| 在线观看66精品国产| 午夜a级毛片| 夜夜爽天天搞| 成年免费大片在线观看| 99久久精品一区二区三区| 一级黄色大片毛片| 99国产极品粉嫩在线观看| 最近视频中文字幕2019在线8| 在线免费观看不下载黄p国产 | 国产在线精品亚洲第一网站| 又粗又爽又猛毛片免费看| 久久久久久大精品| 色综合欧美亚洲国产小说| 成年版毛片免费区| 五月玫瑰六月丁香| 老司机午夜十八禁免费视频| 亚洲中文日韩欧美视频| 在线视频色国产色| 亚洲av熟女| 欧洲精品卡2卡3卡4卡5卡区| 中文亚洲av片在线观看爽| 黄片小视频在线播放| or卡值多少钱| 国产视频内射| 久久久久国产精品人妻aⅴ院| 国产一区二区在线观看日韩 | 麻豆一二三区av精品| 午夜福利在线在线| 亚洲av熟女| 亚洲第一电影网av| 一区二区三区高清视频在线| av在线天堂中文字幕| 婷婷丁香在线五月| 国产精品99久久99久久久不卡| www.色视频.com| 男人的好看免费观看在线视频| 夜夜躁狠狠躁天天躁| 成年女人看的毛片在线观看| 综合色av麻豆| 精品久久久久久成人av| 九九在线视频观看精品| 亚洲熟妇中文字幕五十中出| 欧美中文日本在线观看视频| 欧美区成人在线视频| 深夜精品福利| 美女高潮的动态| 婷婷亚洲欧美| 亚洲熟妇中文字幕五十中出| 性欧美人与动物交配| 亚洲国产精品合色在线| 精品乱码久久久久久99久播| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 天堂√8在线中文| 午夜激情欧美在线| 国产黄片美女视频| 国产成年人精品一区二区| 国产美女午夜福利| 亚洲欧美日韩东京热| 国内揄拍国产精品人妻在线| 给我免费播放毛片高清在线观看| 香蕉久久夜色| 欧美黄色片欧美黄色片| 99久久精品国产亚洲精品| 国产毛片a区久久久久| 久久久久久久久大av| 日韩亚洲欧美综合| 搡老熟女国产l中国老女人| 午夜a级毛片| 青草久久国产| 天天一区二区日本电影三级| 老汉色av国产亚洲站长工具| 久久精品国产99精品国产亚洲性色| 成人国产综合亚洲| 国产成年人精品一区二区| x7x7x7水蜜桃| 九九在线视频观看精品| 免费看美女性在线毛片视频| 黑人欧美特级aaaaaa片| 国产乱人视频| 欧美成人一区二区免费高清观看| 啦啦啦韩国在线观看视频| 村上凉子中文字幕在线| 男人的好看免费观看在线视频| 午夜视频国产福利| 国产真人三级小视频在线观看| 欧美成人一区二区免费高清观看| 中文资源天堂在线| 亚洲av免费在线观看| 三级毛片av免费| 又黄又爽又免费观看的视频| 久久精品综合一区二区三区| 欧美一级毛片孕妇| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 午夜福利成人在线免费观看| 日本免费一区二区三区高清不卡| 免费搜索国产男女视频| 中文字幕高清在线视频| 欧美在线黄色| bbb黄色大片| 亚洲国产精品合色在线| 九九热线精品视视频播放| 亚洲欧美精品综合久久99| 午夜两性在线视频| 国产精品精品国产色婷婷| 亚洲久久久久久中文字幕| 国产一级毛片七仙女欲春2| 露出奶头的视频| 色av中文字幕| 国产精品久久久久久亚洲av鲁大| 男插女下体视频免费在线播放| 久久久久久久久中文| avwww免费| 好男人在线观看高清免费视频| 久久精品国产综合久久久| 男女做爰动态图高潮gif福利片| 国产精品嫩草影院av在线观看 | 首页视频小说图片口味搜索| 日本黄色视频三级网站网址| 91久久精品国产一区二区成人 | 国产精品一及| 老司机午夜福利在线观看视频| 亚洲人成网站在线播放欧美日韩| a级毛片a级免费在线| 内地一区二区视频在线| 国产精品影院久久| 欧美xxxx黑人xx丫x性爽| 亚洲中文日韩欧美视频| 日本精品一区二区三区蜜桃| 女人被狂操c到高潮| 国产精品99久久久久久久久| 国产午夜精品久久久久久一区二区三区 | 亚洲 国产 在线| 亚洲av电影在线进入| 欧美丝袜亚洲另类 | 国产真实乱freesex| 午夜免费男女啪啪视频观看 | 一个人观看的视频www高清免费观看| 精品人妻一区二区三区麻豆 | 亚洲人成网站高清观看| 婷婷亚洲欧美| 啦啦啦免费观看视频1| 中文字幕熟女人妻在线| 熟女人妻精品中文字幕| 99热这里只有是精品50| 一级黄色大片毛片| 在线观看一区二区三区| 欧美zozozo另类| 亚洲av电影不卡..在线观看| 国产探花极品一区二区| 日本五十路高清| 亚洲av电影在线进入| 免费在线观看成人毛片| 成人精品一区二区免费| 宅男免费午夜| 亚洲欧美精品综合久久99| 两个人的视频大全免费| 18美女黄网站色大片免费观看| 色综合婷婷激情| 香蕉av资源在线| 两人在一起打扑克的视频| 午夜两性在线视频| 国产黄色小视频在线观看| 欧美最黄视频在线播放免费| 国产高清激情床上av| 很黄的视频免费| 成人特级黄色片久久久久久久| 精品免费久久久久久久清纯| 99久久99久久久精品蜜桃| 宅男免费午夜| 欧美国产日韩亚洲一区| 精华霜和精华液先用哪个| 精品熟女少妇八av免费久了| 久久精品91蜜桃| 性色av乱码一区二区三区2| 精品国内亚洲2022精品成人| eeuss影院久久| 三级毛片av免费| 亚洲 欧美 日韩 在线 免费| 亚洲精品影视一区二区三区av| 韩国av一区二区三区四区| 岛国视频午夜一区免费看| 欧美极品一区二区三区四区| 在线视频色国产色| 欧美色欧美亚洲另类二区| 日本黄大片高清| 久久亚洲真实| 日韩精品中文字幕看吧| 午夜福利免费观看在线| 中文在线观看免费www的网站| 久久国产乱子伦精品免费另类| 久9热在线精品视频| 婷婷精品国产亚洲av| 亚洲av一区综合| 国产av麻豆久久久久久久| 亚洲激情在线av| 久久草成人影院| 岛国在线观看网站| 久久精品影院6| 在线播放国产精品三级| 国内少妇人妻偷人精品xxx网站| 久久久久久久午夜电影| 美女被艹到高潮喷水动态| 18禁国产床啪视频网站| 亚洲久久久久久中文字幕| 国产三级在线视频| 一级黄片播放器| 精品一区二区三区视频在线观看免费| 99riav亚洲国产免费| 日韩国内少妇激情av| АⅤ资源中文在线天堂| 欧美一级a爱片免费观看看| 热99在线观看视频| 精品久久久久久成人av| 一a级毛片在线观看| 日本黄色片子视频| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线av高清观看| 免费看日本二区| 亚洲av五月六月丁香网| 天堂动漫精品| 美女高潮的动态| 久久精品亚洲精品国产色婷小说| 桃色一区二区三区在线观看| 老司机深夜福利视频在线观看| 男女下面进入的视频免费午夜| 亚洲精华国产精华精| 美女高潮的动态| 女同久久另类99精品国产91| 哪里可以看免费的av片| 国产中年淑女户外野战色| 久久精品国产综合久久久| 不卡一级毛片| 国产成人系列免费观看| 夜夜爽天天搞| 国产欧美日韩精品一区二区| 久久久久亚洲av毛片大全| 母亲3免费完整高清在线观看| 真实男女啪啪啪动态图| 欧美绝顶高潮抽搐喷水| 99久久精品热视频| 久久午夜亚洲精品久久| 波多野结衣巨乳人妻| 人妻夜夜爽99麻豆av| 欧美性感艳星| 国产淫片久久久久久久久 | 天天躁日日操中文字幕| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看| 深夜精品福利| 女警被强在线播放| 国产成+人综合+亚洲专区| 热99在线观看视频| 亚洲最大成人手机在线| 亚洲国产日韩欧美精品在线观看 | 国产色婷婷99| 国产av不卡久久| 99热这里只有是精品50| 十八禁网站免费在线| 久久久久久久亚洲中文字幕 | 午夜福利视频1000在线观看| 精品一区二区三区视频在线观看免费| 久久九九热精品免费| 一级黄色大片毛片| 一本久久中文字幕| 亚洲国产高清在线一区二区三| 久久久久国内视频| 亚洲电影在线观看av| 性色avwww在线观看| 麻豆成人av在线观看| 免费人成视频x8x8入口观看| 黄色成人免费大全| 久久久久九九精品影院| 日本黄色视频三级网站网址| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 毛片女人毛片| 91久久精品电影网| 在线观看免费午夜福利视频| 欧美又色又爽又黄视频| av片东京热男人的天堂| 18禁黄网站禁片免费观看直播| 又爽又黄无遮挡网站| 色在线成人网| 午夜福利成人在线免费观看| 黄色女人牲交| 九色国产91popny在线| 在线免费观看的www视频| 日本三级黄在线观看| 日本 av在线| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品电影| 99久久成人亚洲精品观看| 99久久精品热视频| 两个人看的免费小视频| www.999成人在线观看| 两个人看的免费小视频| 变态另类丝袜制服| 中文字幕av在线有码专区| 国产亚洲欧美在线一区二区| 国产亚洲精品综合一区在线观看| a级毛片a级免费在线| 91在线观看av| 91在线精品国自产拍蜜月 | 好男人在线观看高清免费视频| 身体一侧抽搐| 国产成人aa在线观看| 欧美+亚洲+日韩+国产| 午夜激情福利司机影院| 精品久久久久久,| 首页视频小说图片口味搜索| 在线观看一区二区三区| 午夜久久久久精精品| 国产不卡一卡二| 国产欧美日韩一区二区精品| 免费看光身美女| 久久久国产精品麻豆| 亚洲午夜理论影院| 一夜夜www| www日本在线高清视频| 欧美一区二区亚洲| 国产免费男女视频| 亚洲国产精品成人综合色| 一区二区三区国产精品乱码| 日韩欧美免费精品| 老司机在亚洲福利影院| 日本一二三区视频观看| 日本免费一区二区三区高清不卡| 成年免费大片在线观看| 国产精品香港三级国产av潘金莲| 日韩欧美国产在线观看| 少妇的逼好多水| av专区在线播放| 成人特级av手机在线观看| 欧美乱码精品一区二区三区| 亚洲av电影在线进入| 无限看片的www在线观看| 欧美激情在线99| 三级男女做爰猛烈吃奶摸视频| 99riav亚洲国产免费| 精品一区二区三区人妻视频| 热99re8久久精品国产| 欧美日本亚洲视频在线播放| 精品久久久久久成人av| 亚洲av免费在线观看| 日日夜夜操网爽| 亚洲狠狠婷婷综合久久图片| 久久精品亚洲精品国产色婷小说| 成人午夜高清在线视频| 变态另类成人亚洲欧美熟女| 国产亚洲av嫩草精品影院| 国产高清视频在线播放一区| 免费在线观看成人毛片| 一级a爱片免费观看的视频| 亚洲精品乱码久久久v下载方式 | 色综合婷婷激情| 日韩欧美免费精品| 激情在线观看视频在线高清| 中国美女看黄片| 色哟哟哟哟哟哟| 熟女人妻精品中文字幕| 免费看a级黄色片| 精品日产1卡2卡| 亚洲激情在线av| 精品国内亚洲2022精品成人| 中文字幕熟女人妻在线| 欧美精品啪啪一区二区三区| 久久精品91无色码中文字幕| 一区二区三区激情视频| 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线观看二区| 91麻豆av在线| 久久99热这里只有精品18| 国产成人a区在线观看| 国产成年人精品一区二区| 成人精品一区二区免费| 男女之事视频高清在线观看| 波野结衣二区三区在线 | 国产精品99久久99久久久不卡| 五月玫瑰六月丁香| 高潮久久久久久久久久久不卡| 成年女人永久免费观看视频| 国产免费男女视频| 精品久久久久久久末码| 老熟妇仑乱视频hdxx| 欧美3d第一页| 免费av毛片视频| 日日夜夜操网爽| 超碰av人人做人人爽久久 | 舔av片在线| 亚洲最大成人中文| 禁无遮挡网站| 日韩人妻高清精品专区| 色av中文字幕| 亚洲精品美女久久久久99蜜臀| 9191精品国产免费久久| 亚洲精品影视一区二区三区av| 精品国产超薄肉色丝袜足j| 精品一区二区三区视频在线 | 精品久久久久久久久久久久久| 十八禁网站免费在线| 欧美日本视频| 亚洲精华国产精华精| 免费在线观看影片大全网站| 成人鲁丝片一二三区免费| 他把我摸到了高潮在线观看| 天堂av国产一区二区熟女人妻| 超碰av人人做人人爽久久 | 性色av乱码一区二区三区2| 国语自产精品视频在线第100页| 欧美乱妇无乱码| 亚洲av成人不卡在线观看播放网| 亚洲av不卡在线观看| 丰满的人妻完整版| 啦啦啦观看免费观看视频高清| 午夜福利免费观看在线| 婷婷精品国产亚洲av| 精品欧美国产一区二区三| 99久久久亚洲精品蜜臀av| 欧美3d第一页| 亚洲,欧美精品.| netflix在线观看网站| 啪啪无遮挡十八禁网站| ponron亚洲| 色综合欧美亚洲国产小说| 日本精品一区二区三区蜜桃| 国产免费男女视频| 亚洲最大成人手机在线| 亚洲成人精品中文字幕电影| xxx96com| 观看美女的网站| 亚洲国产高清在线一区二区三| aaaaa片日本免费| av欧美777| 精品国内亚洲2022精品成人| 最近最新中文字幕大全免费视频| 不卡一级毛片| 99久久综合精品五月天人人| 国产高清视频在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看| 俺也久久电影网| 久久久久亚洲av毛片大全| 久久伊人香网站| 欧美一区二区国产精品久久精品| 少妇高潮的动态图| 别揉我奶头~嗯~啊~动态视频| 五月伊人婷婷丁香| 怎么达到女性高潮| 亚洲欧美一区二区三区黑人| 精品国产亚洲在线| 国产单亲对白刺激| 波多野结衣高清作品| 两性午夜刺激爽爽歪歪视频在线观看| 欧美成狂野欧美在线观看| 精品久久久久久,| 亚洲av成人精品一区久久| 午夜亚洲福利在线播放| 身体一侧抽搐| 亚洲av电影不卡..在线观看| 久久香蕉国产精品| 老熟妇乱子伦视频在线观看| 无人区码免费观看不卡| 日韩欧美精品v在线| 亚洲国产中文字幕在线视频| 日本 av在线| a在线观看视频网站| 国产精品 欧美亚洲| 国产蜜桃级精品一区二区三区| 中文在线观看免费www的网站| 久久香蕉国产精品| 国产亚洲精品av在线| 琪琪午夜伦伦电影理论片6080| 免费人成在线观看视频色| 天堂动漫精品| 日日摸夜夜添夜夜添小说| 极品教师在线免费播放| 国语自产精品视频在线第100页| 亚洲av免费高清在线观看| 亚洲 国产 在线| 国产精品久久久久久久久免 | 成人精品一区二区免费| 欧美性猛交╳xxx乱大交人| 亚洲国产日韩欧美精品在线观看 | 欧美精品啪啪一区二区三区| 国产黄色小视频在线观看| 国内精品一区二区在线观看| 丰满乱子伦码专区| 中文字幕高清在线视频| 国产一区二区在线av高清观看| 欧美一级a爱片免费观看看| 免费高清视频大片| 日日摸夜夜添夜夜添小说| 精品久久久久久久毛片微露脸| 国产午夜精品论理片| 最近最新中文字幕大全免费视频| 黄色日韩在线| 国产亚洲av嫩草精品影院| 日韩欧美 国产精品| 一级黄片播放器| 村上凉子中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩一区二区三| 亚洲在线观看片| 国产精品日韩av在线免费观看| 韩国av一区二区三区四区| 日日干狠狠操夜夜爽| 三级男女做爰猛烈吃奶摸视频| 国产亚洲av嫩草精品影院| 日韩人妻高清精品专区| 色视频www国产| 少妇丰满av| 欧美大码av| 亚洲国产欧美人成| 日韩国内少妇激情av| 午夜老司机福利剧场| 日韩免费av在线播放| 免费看日本二区| 免费无遮挡裸体视频| 精品久久久久久久毛片微露脸| 亚洲精品456在线播放app | 俺也久久电影网| 久久久久久九九精品二区国产| 99久久精品国产亚洲精品| 国产一区二区在线观看日韩 | 高清毛片免费观看视频网站| 不卡一级毛片| 国产亚洲欧美98| 国产精品久久久人人做人人爽| 国产伦精品一区二区三区四那| 亚洲av第一区精品v没综合| 三级男女做爰猛烈吃奶摸视频| 亚洲人与动物交配视频| 亚洲片人在线观看| 亚洲在线观看片| 亚洲av电影在线进入| 99久久九九国产精品国产免费|