楊 靖

(上海核工程研究設(shè)計(jì)院，上海 200233)

?

Fluidyn-Pollusol在核電廠地下水?dāng)?shù)值模擬中的應(yīng)用

楊 靖

(上海核工程研究設(shè)計(jì)院，上海 200233)

介紹了Fluidyn-Pollusol模型的基本理論，從輸入、運(yùn)行、輸出三方面，闡述了Fluidyn-Pollusol的主要處理模塊，并以某核電廠工程為例，分析了Fluidyn-Pollusol在地下水?dāng)?shù)值模擬中的應(yīng)用，以供參考借鑒。

Fluidyn-Pollusol，數(shù)值模型，地下水，核電廠

0 引言

目前我國核電處于快速增長的階段，隨著國內(nèi)外核電廠的選址建設(shè)，核電廠對(duì)地下水環(huán)境影響特別是內(nèi)陸核電廠址中放射性核素在地下水中遷移、吸附和沉積問題是一個(gè)需要特別關(guān)注的問題。尤其在2011年3月11日，日本福島核電站事件后，地震海嘯引起核電站故障并造成放射性物質(zhì)向地下水中釋放引起公眾對(duì)核電站的環(huán)境影響問題的關(guān)注空前加大。數(shù)值模擬方法是目前模擬放射性核素在地下水中遷移的主要技術(shù)方法之一[1]。這種方法以其有效性、靈活性和經(jīng)濟(jì)性成為地下水環(huán)境評(píng)價(jià)領(lǐng)域一種不可或缺的重要方法。在眾多的地下水模擬軟件中，F(xiàn)luidyn-Pollusol主要應(yīng)用于國外的核電行業(yè)地下水環(huán)境評(píng)價(jià)，例如WESTLAKES Scientific Consulting(英國WESTLAKES)，NUPEC (日本動(dòng)力工程設(shè)計(jì)試驗(yàn)中心)，CEA (法國原子能委員會(huì))，EDF (法國電力公司)，AREVA (法國阿海琺)，ANDRA (法國放射性廢棄物管理局)，IGCAR (印度甘地原子能研究中心)，NPCIL (印度國有核電公司)，BARC (印度巴巴原子研究中心)。

1 Fluidyn-Pollusol模型基本理論

Fluidyn-Pollusol地下水?dāng)?shù)值模型是基于Navier-Stokes方程[2]和達(dá)西定律[3]進(jìn)行的，Navier-Stokes方程詳見下式：

(1)

(2)

(3)

(4)

式中：ρ——流體密度；u，v，w——x，y，z方向上的流速；Kx，Ky，Kz——x，y，z方向有效水力傳導(dǎo)率；Sm——質(zhì)量源項(xiàng)；Su，Sv，Sw——x，y，z方向上的動(dòng)量源項(xiàng)。

在考慮地下水的對(duì)流、水動(dòng)力彌散、流體源匯項(xiàng)、平衡吸附以及一級(jí)不可逆反應(yīng)(如放射性衰變)的情況下，單一化學(xué)組分遷移的數(shù)學(xué)方程為：

2 Fluidyn-Pollusol主要處理模塊

2.1 輸入模塊

建立模型就是建立概念模型并將之轉(zhuǎn)化為數(shù)字形式，用于模擬地下水變化過程。將所評(píng)價(jià)區(qū)域劃分為若干層網(wǎng)狀單元格，確定邊界類型，然后根據(jù)含水層地質(zhì)巖性資料及其他水文特征確定傳導(dǎo)系數(shù)K，入滲系數(shù)R，儲(chǔ)水系數(shù)S及初始水位等各種參數(shù)。Fluidyn-Pollusol輸入支持DXF格式文件，可以實(shí)現(xiàn)快速建模，還能夠?qū)氲貓D(bmp格式)文件，可以輔助模型建立。

2.2 運(yùn)行模塊

Fluidyn-Pollusol可以用于模擬飽和流狀態(tài)的水流和溶質(zhì)運(yùn)移問題。可以模擬地下水中持續(xù)或偶發(fā)污染物(含放射性核素)組分的傳播，熱傳遞，可變密度流場(鹽水或海水入侵問題)以及非飽和帶流場以及物質(zhì)運(yùn)移問題。

2.3 輸出模塊

模型輸出允許用戶以三種不同方式展示其模擬結(jié)果。第一種方式就是在計(jì)算機(jī)屏幕上直接彩色立體顯示所有的模擬結(jié)果；第二種方式就是直接在各類打印機(jī)上輸出各種模擬評(píng)價(jià)的成果表格和成果圖鑒；最后一種方式就是將所有模擬結(jié)果以圖形或文本的文件格式輸出，輸出圖形包括可以標(biāo)記出滲流速度矢量大小的平面、剖面等值線和平面、剖面示蹤流線圖以及局部區(qū)域水均衡圖等一系列圖件。

3 模型應(yīng)用實(shí)例

3.1 水文地質(zhì)概況

基巖區(qū)巖性主要為中風(fēng)化板巖，回填區(qū)為人工填土，地下水流動(dòng)方向?yàn)橹鲝S房基巖區(qū)→回填區(qū)→廠外第四系地層→地表水。

計(jì)算時(shí)研究的是放射性同位素氚。這主要是因?yàn)樵诤穗姀S運(yùn)行時(shí)會(huì)產(chǎn)生大量的氚且不能通過常規(guī)的廢液處理系統(tǒng)進(jìn)行處理。假設(shè)污染源位置位于流出物廢液收集箱附近，假設(shè)其活度值為3.74×1010Bq/m3，釋放形式為瞬時(shí)釋放。

3.2 Fluidyn-Pollusol模型方法

在模擬放射性核素在地下水中遷移時(shí)，我們一般僅分析潛水層的情況。所以將模型中含水層劃分為一層，將模擬區(qū)域劃分為若干個(gè)單元網(wǎng)格。然后確定邊界條件、確定層頂部與底部高程、確定時(shí)間、初始水位、水平傳導(dǎo)率、儲(chǔ)水系數(shù)等相關(guān)參數(shù)，并將其與單元格相對(duì)應(yīng)。

一般做核電廠地下水環(huán)境模擬時(shí)，所建立的是一個(gè)長期模型，時(shí)間為90年。這個(gè)預(yù)測時(shí)間范圍主要考慮了我國核電廠的使用壽命以及核電廠退役所需時(shí)間，并對(duì)該時(shí)間適當(dāng)放寬，預(yù)測核電廠如果存在放射性核素對(duì)地下水環(huán)境的影響。

利用前期勘測資料中的地下水水位對(duì)模型進(jìn)行校準(zhǔn)，并對(duì)設(shè)置的參數(shù)進(jìn)行調(diào)整。模型校準(zhǔn)過程就是調(diào)整參數(shù)的過程，使得水位觀測值與模擬值最大程度的一致。

3.3 模型預(yù)測結(jié)果

對(duì)于土壤中不被吸附滯留的核素氚，會(huì)隨著地下水水流方向運(yùn)動(dòng)。如果污染源在核電廠放射性廢液收集箱附近，根據(jù)該廠址的廠址特性需要較長的時(shí)間內(nèi)才能產(chǎn)生明顯的遷移距離。本文實(shí)例中，氚大約需要在20年，濃度峰才能到達(dá)地表水，并被地表水稀釋。氚在地下水中遷移的速度主要由地下水巖土層特征所決定。模擬區(qū)邊界觀測點(diǎn)H-3濃度變化曲線見圖1。

4 結(jié)語

以上介紹了Fluidyn-Pollusol模型的原理、模塊、使用方法等等。Fluidyn-Pollusol的使用多見于國外的核電工程，具有獨(dú)特的優(yōu)越性?？梢灶A(yù)見Fluidyn-Pollusol地下水模型將在國內(nèi)核電工程地下水環(huán)境影響評(píng)價(jià)中發(fā)揮越來越大的作用。

[1] 郭曉東,田 輝，張梅佳，等.我國地下水?dāng)?shù)值模擬軟件應(yīng)用進(jìn)展[J].地下水,2010,32(4):5-7.

[2] 約翰D，安德森.計(jì)算流體力學(xué)基礎(chǔ)及其應(yīng)用[M].吳頌平，劉趙森，譯.北京:機(jī)械工業(yè)出版社,2007:170-172.

[3] 薛禹群，謝春紅.地下水?dāng)?shù)值模擬[M].北京:科學(xué)出版社,2007:12-13.

On application of Fluidyn-Pollusol in underground water numeric simulation in nuclear power plant

Yang Jing

(Shanghai Nuclear Engineering Research and Design Institute, Shanghai 200233, China)

The paper introduces the basic theory of Fluidyn-Pollusol model, illustrates its main processing model from the input, operation, and output, and analyzes its application of the underground water numeric simulation by taking some nuclear power plant as the example, so as to provide some reference.

Fluidyn-Pollusol, numeric model, underground water, nuclear power plant

1009-6825(2016)12-0193-02

2016-02-15

楊 靖(1982- )，女，工程師

TM623

A