基于GMS的數(shù)值模擬在某化工園地下水環(huán)境影響評(píng)價(jià)中的應(yīng)用

吳鵬飛, 彭 展, 陳小婷

(湖北省地質(zhì)局 第三地質(zhì)大隊(duì),湖北 黃岡 438000)

基于GMS的數(shù)值模擬在某化工園地下水環(huán)境影響評(píng)價(jià)中的應(yīng)用

吳鵬飛, 彭 展, 陳小婷

(湖北省地質(zhì)局 第三地質(zhì)大隊(duì),湖北 黃岡 438000)

在建設(shè)項(xiàng)目地下水環(huán)境影響評(píng)價(jià)中,通過對(duì)場地進(jìn)行水文地質(zhì)調(diào)查分析,運(yùn)用基于GMS的數(shù)值模擬分析建立三維數(shù)值模型,是評(píng)價(jià)地下水污染物運(yùn)移規(guī)律、污染范圍確定及濃度分布的重要手段。通過數(shù)值模擬分析,研究事故條件和非正常狀況條件下,地下水特征污染因子運(yùn)移規(guī)律,及在地下水環(huán)境中的濃度變化情況,對(duì)一定時(shí)間內(nèi)地下水可能受到污染的情況進(jìn)行分析和預(yù)測。

地下水環(huán)境影響評(píng)價(jià);GMS;數(shù)值模擬;化工園

中國地下水資源地域分布不均,地下水資源量總體上呈東南向西北逐漸降低的規(guī)律[1]。隨著社會(huì)發(fā)展,造成地下水的過度開采和污染等問題日益突出,地下水資源保護(hù)越來越重要。其中,開展建設(shè)項(xiàng)目地下水環(huán)境影響評(píng)價(jià)是項(xiàng)目建設(shè)中對(duì)地下水資源保護(hù)的重要預(yù)防手段之一。

針對(duì)地下水環(huán)境影響一級(jí)評(píng)價(jià),通常采用數(shù)值法進(jìn)行分析。本文以某在建化工園區(qū)為例,在研究園區(qū)水文地質(zhì)條件的基礎(chǔ)上,運(yùn)用GMS中的MODFLOW模塊建立三維數(shù)值模型,進(jìn)行地下水流場數(shù)值模擬。主要是分析事故條件和非正常狀況條件下,化工園區(qū)特征污染因子在地下水環(huán)境中的運(yùn)移特征及濃度變化情況,從而預(yù)測地下水可能受到污染的情況,為地下水保護(hù)提出針對(duì)性的預(yù)防措施。隨著計(jì)算機(jī)和數(shù)字化技術(shù)的飛速發(fā)展,采用數(shù)值法對(duì)項(xiàng)目建設(shè)中的地下水污染進(jìn)行分析評(píng)價(jià)、實(shí)施地下水資源保護(hù)已逐步成為常用的有效方法,值得推薦,并在實(shí)踐應(yīng)用中不斷總結(jié)提高和完善。

1 項(xiàng)目概況

1.1 基本情況

該化工園東西跨度6.4 km,南北跨度2.5 km,規(guī)劃面積20.32 km2,建有工業(yè)水廠、污水處理廠、消防中隊(duì)、?；穼＞€、?；坟泩?、?；反a頭等專用設(shè)施,同時(shí)規(guī)劃有基礎(chǔ)化工區(qū)、精細(xì)化工區(qū)、發(fā)展備用區(qū)和商住配套區(qū)四大功能區(qū),其中基礎(chǔ)化工區(qū)6 km2,精細(xì)化工區(qū)5 km2,商住及配套區(qū)5 km2,發(fā)展備用區(qū)4 km2。產(chǎn)業(yè)發(fā)展重點(diǎn)為液化天然氣、硫化工、80萬t乙烯配套項(xiàng)目、鈦化工、精細(xì)化工及醫(yī)藥化工等,目前建成投產(chǎn)、在建及待建各類企業(yè)50余家。本次模擬分析評(píng)價(jià)選取兩家代表企業(yè)作為模擬分析對(duì)象。

1.2 園區(qū)地質(zhì)環(huán)境特征

1.2.1 地形地貌

園區(qū)所在地屬丘陵崗地地貌,地勢較為平坦開闊,大氣擴(kuò)散條件較好,整個(gè)場地被第四紀(jì)沖積土覆蓋,呈黃色的亞粘土,土層地質(zhì)特性良好,地下水位低。園區(qū)地勢為西北高、東南低,最高海拔高程95.7 m,最低22.1 m。

1.2.2 氣象水文

園區(qū)屬亞熱帶大陸性季風(fēng)氣候,江淮小氣候區(qū),四季分明。全年太陽輻射量106.49～113.31千卡/cm2,年均日照時(shí)數(shù)1 959.4 h,西北部日照高于東南部,全市日照率為43%～49%;年平均氣溫為15.7～17.1 ℃;最高氣溫38.1 ℃,最低氣溫-5.3 ℃,年平均相對(duì)濕度77%,年平均氣壓1 010.6 hPa。全年無霜期在237—278 d,年平均降雨量1 223～1 493 mm,年降水總量222.37億m3,降雨日數(shù)(≥0.1 mm/d)在115—147 d。常年主導(dǎo)風(fēng)向N、頻率為19%,次主導(dǎo)風(fēng)向SE、頻率為14%,常年平均風(fēng)速1.8 m/s。

化工園區(qū)東側(cè)為巴河,與長江相連。園區(qū)距離長江約20 km。園區(qū)內(nèi)地表水以水塘為主,較大部分為近年來人工開挖魚塘、藕塘,深度約1～3 m。

1.2.3 地層巖性

1.2.4 水文地質(zhì)條件

區(qū)內(nèi)地下水主要為松散巖類孔隙水、基巖裂隙水兩種類型。

(1) 松散巖類孔隙水:主要賦存于第四系松散巖組的孔隙中,其賦水空間有限,一般受大氣降水及人工排水補(bǔ)給,水位水量隨季節(jié)而變化,水量較小。受地形、地貌控制,只有一定匯水面積的掌、杖形地及沖溝溝腦地帶有泉水出露。在地形切割較弱、海拔較高、降水量較大的地段,單泉流量0.1～1 L/s;地下水水化學(xué)類型屬HCO3—Ca·Mg型,礦化度0.5～1.0 g/L。該層地下水季節(jié)變化明顯,枯水期水量較小,雨季相對(duì)較大,但由于地表徑流排泄快,雨水滲透補(bǔ)給地下水的量有限,因此富水性較差,水量貧乏。

(2) 基巖裂隙水:主要賦存于變質(zhì)巖的風(fēng)化、構(gòu)造裂隙中,接受大氣降水及上層的孔隙水的滲入補(bǔ)給。區(qū)內(nèi)下伏基巖局部節(jié)理裂隙較發(fā)育,但多被次生礦物充填,連通性較差,水量甚微,為富水性微弱的裂隙含水層,上覆土層透水性弱,其儲(chǔ)水空間有限,地下水水量較貧乏。

1.3 園區(qū)地下水利用現(xiàn)狀及污染源調(diào)查

化工園區(qū)內(nèi)僅有少量分散式民井,居民飲用自來水,民井部分廢棄,少量作為生活用水,如洗衣、澆菜等。根據(jù)園區(qū)規(guī)劃,化工園建成后不設(shè)置居民區(qū)[3],現(xiàn)有村莊將按照規(guī)劃整體搬遷,現(xiàn)存民井將全部廢棄。總體來說,園區(qū)屬變質(zhì)巖貧水區(qū),地下水利用程度較低。

化工園地下水主要污染源有工業(yè)企業(yè)生產(chǎn)廢水和居民生活污水。工業(yè)企業(yè)生產(chǎn)廢水:整個(gè)園區(qū)現(xiàn)狀企業(yè)生產(chǎn)廢水排放量是83萬t/年,其中主要污染因子為COD、氨氮,COD為167.9 t/年;氨氮為12.9 t/年。生活污水:根據(jù)初步調(diào)查統(tǒng)計(jì),園區(qū)現(xiàn)有居民生活污水排放量為530.64 t/年。

2 地下水流數(shù)值模擬和污染物溶質(zhì)運(yùn)移數(shù)值模擬

2.1 水文地質(zhì)概念模型

水文地質(zhì)概念模型是把含水層或含水系統(tǒng)實(shí)際的邊界性質(zhì)、內(nèi)部結(jié)構(gòu)、滲透性能、水力特征和補(bǔ)給排泄等條件進(jìn)行合理的概化,以便可以進(jìn)行數(shù)學(xué)與物理模擬[4]?？茖W(xué)、準(zhǔn)確地建立水文地質(zhì)概念模型是地下水環(huán)境影響預(yù)測評(píng)價(jià)的關(guān)鍵。本次數(shù)值模擬分析的水文地質(zhì)概念模型設(shè)置如下。

(1) 含水層概化:園區(qū)含水巖組主要為第四系松散巖類孔隙水和變質(zhì)巖基巖裂隙水。

(2) 邊界條件:北部、東部邊界為給定水頭邊界;其他邊界為零通量邊界。上邊界為降水補(bǔ)給、蒸發(fā),下邊界取中風(fēng)化層底板,等效定義為零通量邊界(相對(duì)隔水邊界)。

(3) 源匯項(xiàng)設(shè)定:園區(qū)含水層主要接受大氣降雨補(bǔ)給,最終排泄至巴河。

(4) 初始條件確定:在模型中輸入收集的園區(qū)監(jiān)測孔的長期觀測水位,可生成初始等水位線,確定園區(qū)地下水的初始流場。

2.2 滲流模型

2.2.1 數(shù)學(xué)方程與求解平臺(tái)

通過對(duì)水文地質(zhì)概念模型的分析,依據(jù)滲流連續(xù)性方程和達(dá)西定律,建立模擬區(qū)地下水系統(tǒng)水文地質(zhì)概念模型相對(duì)應(yīng)的三維非穩(wěn)定流數(shù)學(xué)模型:

上述數(shù)學(xué)控制方程的求解平臺(tái)采用DHI-WASY公司開發(fā)的基于有限單元法的FEFLOW(Finite Element subsurface FLOW system)軟件。

在眾多模擬軟件中,由德國水資源規(guī)劃與系統(tǒng)研究所(WASY)開發(fā)出來的地下水流動(dòng)及物質(zhì)遷移模擬軟件系統(tǒng)FEFLOW具有獨(dú)到的特點(diǎn),它是迄今為止功能最為齊全的地下水模擬軟件包之一,可用于復(fù)雜三維非穩(wěn)定水流和污染物運(yùn)移的模擬。

2.2.2 初始網(wǎng)格與地質(zhì)模型

基于FEFLOW平臺(tái),輸入模擬區(qū)域矢量數(shù)據(jù),利用Advancing Front剖分方法,將區(qū)域離散為不規(guī)則三角剖分網(wǎng)格,剖分過程嚴(yán)格遵循Delaunay法則,使三角網(wǎng)格內(nèi)的三角形內(nèi)角角度為銳角,三邊長度盡量相等,三角形網(wǎng)中任一三角形的外接圓范圍內(nèi)不會(huì)有其他點(diǎn)存在,在散點(diǎn)集可能形成的三角剖分中,Delaunay三角剖分所形成的三角形的最小角最大。

在園區(qū)適當(dāng)加密三角剖分單元,最終得到模擬區(qū)初始二維剖分結(jié)果如圖1所示,其中結(jié)點(diǎn)數(shù)5 190個(gè),有限單元數(shù)10 033個(gè)。

根據(jù)水文地質(zhì)概念模型,地質(zhì)模型(含水系統(tǒng))由潛水含水層構(gòu)造,共分為三層(layer)四片(slice)。

三層:第一層為第四系素填土,第二層為粉質(zhì)粘土,第三層為片麻巖層。

四片:地表高程、第四系素填土層底板、粉質(zhì)粘土層底板、片麻巖層底板。

其中地表高程數(shù)據(jù)采用ASTER GDEM數(shù)據(jù)(數(shù)據(jù)來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心科學(xué)數(shù)據(jù)中心),利用ESRI公司的ArcGIS軟件處理以上數(shù)據(jù),輸入FEFLOW后,即可建立模擬區(qū)三維地質(zhì)模型,如圖2所示,其中結(jié)點(diǎn)數(shù)10 380個(gè),有限單元數(shù)23 590個(gè)。

圖1 模擬區(qū)二維網(wǎng)格剖分Fig.1 Two dimensional grid dissection of simulation area

圖2 模擬區(qū)三維網(wǎng)格剖分Fig.2 3d grid profile of simulation area

2.2.3 水文地質(zhì)參數(shù)

本次模擬工作所用到的初始水文地質(zhì)參數(shù)主要依據(jù)園區(qū)布設(shè)的7個(gè)地下水監(jiān)測點(diǎn)資料及試驗(yàn)獲取的水文地質(zhì)參數(shù),同時(shí)根據(jù)模擬區(qū)水文地質(zhì)概念模型,對(duì)其滲透系數(shù)進(jìn)行了概化分區(qū),水文地質(zhì)初始參數(shù)取值如表1所示,第一層為第四系素填土,第二層為粉質(zhì)粘土,第三層為片麻巖層。

表1 評(píng)價(jià)區(qū)水文地質(zhì)初始參數(shù)取值表Table 1 Evaluation of the initial parameters of hydrogeological parameters

本次溶質(zhì)運(yùn)移模型中彌散度的確定主要依據(jù)是Geihar等(1992)對(duì)世界范圍內(nèi)所收集的59個(gè)大區(qū)域彌散資料進(jìn)行的整理分析。按照偏保守原則,最終確定的溶質(zhì)運(yùn)移模型參數(shù)見表2。

表2 溶質(zhì)運(yùn)移模型參數(shù)表Table 2 Parameter list of solute transport model

2.2.4 識(shí)別驗(yàn)證

利用正演試錯(cuò)法,將模擬區(qū)邊界條件、參數(shù)分區(qū)、參數(shù)取值等輸入模型,并執(zhí)行正演模擬,反復(fù)調(diào)整需要識(shí)別的參數(shù),直到模型結(jié)果與現(xiàn)狀調(diào)查中的水位觀測點(diǎn)擬合程度較好為止。

在參數(shù)識(shí)別基礎(chǔ)上,調(diào)整模型為非穩(wěn)定流模式,設(shè)置時(shí)間為30年,觀察水位觀測點(diǎn)的動(dòng)態(tài)特征,并記錄模型水均衡數(shù)據(jù)。對(duì)出現(xiàn)水動(dòng)態(tài)異常、水均衡失穩(wěn)等情況的識(shí)別結(jié)果,重新開展參數(shù)識(shí)別,直到識(shí)別結(jié)果能通過驗(yàn)證工作的檢驗(yàn)。

對(duì)調(diào)查評(píng)價(jià)區(qū)的地下水滲流場進(jìn)行參數(shù)識(shí)別,經(jīng)反復(fù)調(diào)整參數(shù)和均衡量,識(shí)別水文地質(zhì)條件,確定模型結(jié)構(gòu)、參數(shù)和均衡要素;當(dāng)所取參數(shù)如表3時(shí),模擬的地下水流場與實(shí)際地下水流場基本一致,識(shí)別的水文地質(zhì)參數(shù)符合實(shí)際水文地質(zhì)條件,基本反映了地下水系統(tǒng)的水力特征,可利用模型進(jìn)行地下水位預(yù)報(bào)。

表3 模型擬合水文地質(zhì)參數(shù)列表Table 3 List of hydrogeological parameters

依據(jù)識(shí)別后的參數(shù),水位擬合情況如圖3所示,擬合數(shù)據(jù)與實(shí)測水位數(shù)據(jù)基本一致,模擬結(jié)果可靠。

圖3 初始流場水位擬合折線圖Fig.3 Fitting line drawing of initial flow field water level

2.2.5 初始條件

經(jīng)參數(shù)識(shí)別驗(yàn)證后,運(yùn)用模擬所得水文地質(zhì)參數(shù),可得初始地下水流場圖如圖4所示。將其作為模擬的初始條件,使得其基本反映出實(shí)際地下水流場特征,可以此為基礎(chǔ)開展后續(xù)地下水環(huán)境影響預(yù)測評(píng)價(jià)工作。

圖4 初始流場圖Fig.4 Initial flow field diagram

2.3 模擬因子選擇和工況設(shè)定

根據(jù)擬建項(xiàng)目特點(diǎn),施工期及服役期滿后污染極小,主要產(chǎn)污時(shí)段為運(yùn)營期,故選取運(yùn)營期作為總模擬時(shí)間,假定時(shí)長為30年。計(jì)算時(shí)間步長為自適應(yīng)模式,記錄第100 d、1 000 d及每年的模擬預(yù)測結(jié)果,共計(jì)32個(gè)時(shí)間點(diǎn)的數(shù)據(jù),為污染物遷移規(guī)律的分析工作提供數(shù)據(jù)支撐。

2.3.1 模擬預(yù)測因子

本項(xiàng)目針對(duì)7個(gè)地下水監(jiān)測點(diǎn)采集7個(gè)樣品進(jìn)行水質(zhì)監(jiān)測分析,通過收集前期的資料及地下水水質(zhì)檢測報(bào)告得知,擬建項(xiàng)目生產(chǎn)過程中產(chǎn)生的廢水主要成分為COD、BOD5、SS、氨氮、重金屬等,結(jié)合本項(xiàng)目的工程特點(diǎn),選定COD作為擬建項(xiàng)目地下水污染特征因子。

2.3.2 工況設(shè)定

工況1:非正常狀況下。

模擬情景:根據(jù)《環(huán)境影響評(píng)價(jià)技術(shù)導(dǎo)則——地下水》(HJ610—2016),非正常排放情況下,預(yù)測源強(qiáng)可考慮設(shè)施老化或腐蝕情況。對(duì)于本項(xiàng)目地下水污染非正常排放源強(qiáng),模擬企業(yè)1、企業(yè)2兩家企業(yè)污染裝置、防滲裝置發(fā)生老化的情景,防滲等級(jí)降至10-5cm/s,污染物發(fā)生滲透。兩家企業(yè)的位置如圖5所示。

圖5 預(yù)測評(píng)價(jià)企業(yè)位置示意圖Fig.5 Position sketch of predictive evaluation enterprise

模擬污染物:COD。污染源概化:連續(xù)恒定排放,面源。泄漏點(diǎn):企業(yè)1、企業(yè)2。泄漏面積:分別為235 m2、145 m2。泄漏時(shí)間:持續(xù)性泄露,共30年。泄漏濃度:COD初始濃度分別為500 mg/L、800 mg/L。

工況2:事故狀況下。

模擬情景:事故排放源強(qiáng)主要考慮企業(yè)1、企業(yè)2兩家企業(yè)污染物直接滲入地下含水層中情景。

模擬污染物:COD。污染源概化:瞬時(shí)排放,面源。泄漏點(diǎn):污染裝置破裂導(dǎo)致泄漏。泄漏面積:分別為235 m2、145 m2。泄漏時(shí)間:泄漏10 d。泄漏濃度:COD初始濃度分別為500 mg/L、800 mg/L。

2.4 污染物濃度變化趨勢

2.4.1 非正常狀況評(píng)價(jià)結(jié)果

污染裝置發(fā)生老化后,污染物下滲進(jìn)入地下水中,形成超標(biāo)污染暈,其遷移方向主要受水動(dòng)力場控制。本評(píng)價(jià)區(qū)域,污染物逐步向東南方向遷移擴(kuò)散,污染范圍持續(xù)擴(kuò)大,選取100 d、1 000 d、10 950 d三個(gè)時(shí)間段進(jìn)行統(tǒng)計(jì)。

針對(duì)三個(gè)典型時(shí)間段,統(tǒng)計(jì)污染暈的運(yùn)移距離、污染面積,見表4。

圖6展示了模型運(yùn)行100 d、1 000 d、10 950 d天三個(gè)時(shí)段下地下水中污染物的遷移擴(kuò)散情況。

圖6 非正常狀況下COD滲漏超標(biāo)污染暈平面圖Fig.6 COD leakage exceeds the level of pollution in abnormal condition

2.4.2 事故狀況下評(píng)價(jià)結(jié)果

泄露事故發(fā)生后,污染物下滲進(jìn)入地下水中,形成超標(biāo)污染暈,遷移方向主要受水動(dòng)力場控制,逐步向東南部擴(kuò)散,污染暈逐步擴(kuò)大,在地下水稀釋作用下,最高濃度逐步降低,在10 950 d(30年),超標(biāo)污染暈最大、濃度最小。

表4 非正常狀況下COD超標(biāo)污染暈預(yù)測結(jié)果Table 4 The result of abnormal COD pollution in abnormal conditions

針對(duì)三個(gè)典型時(shí)間段,統(tǒng)計(jì)了污染暈的運(yùn)移距離、污染面積,數(shù)據(jù)見表5。

表5 事故情景下COD超標(biāo)污染暈預(yù)測結(jié)果Table 5 The result of the prediction of COD contamination in the accident situation

模型運(yùn)行100 d、1 000 d、3 650 d三個(gè)時(shí)段下地下水中污染物的遷移擴(kuò)散情況,如圖7。

圖7 事故情景下COD滲漏超標(biāo)污染暈平面圖Fig.7 COD leakage exceeds the level of pollution in the accident situation

2.4.3 評(píng)價(jià)結(jié)論

本次評(píng)價(jià)采用三維非穩(wěn)定流水文地質(zhì)概念模型,對(duì)園區(qū)代表企業(yè)在非正常和風(fēng)險(xiǎn)事故情景下進(jìn)行了預(yù)測分析,模擬結(jié)果顯示:非正常狀況下,污染物下滲進(jìn)入地下水中,形成超標(biāo)污染暈,其遷移方向主要受水動(dòng)力場控制,逐步向東南部擴(kuò)散,污染暈面積持續(xù)擴(kuò)大,至第30年時(shí),超標(biāo)污染暈范圍最大;事故情景下,污染物下滲進(jìn)入地下水中,形成超標(biāo)污染暈,呈現(xiàn)面積逐漸擴(kuò)大、濃度逐漸降低的特點(diǎn),至第30年時(shí),污染物濃度最低。

在非正常狀況及事故情景下,超標(biāo)污染暈呈現(xiàn)逐步擴(kuò)大的趨勢,運(yùn)移出廠界,項(xiàng)目停止運(yùn)行后,超標(biāo)污染暈有可能繼續(xù)向東南擴(kuò)散,建議在污染裝置下布設(shè)防滲膜等防滲措施,并在其下游布設(shè)監(jiān)測井和應(yīng)急抽水井,防止地下水污染物對(duì)廠區(qū)外地下水環(huán)境造成的影響。

3 結(jié)論與建議

運(yùn)用GMS軟件,建立某化工園區(qū)水文地質(zhì)概念模型,經(jīng)過模型的識(shí)別驗(yàn)證,模型模擬效果較好。預(yù)測結(jié)果表明:非正常狀況下,污染物下滲進(jìn)入地下水中,形成超標(biāo)污染暈,逐步向東南部擴(kuò)散,遷移距離最大達(dá)752 m,污染暈面積持續(xù)擴(kuò)大,至第30年時(shí),超標(biāo)污染暈范圍最大,污染面積達(dá)184 326 m2;事故情景下,污染物下滲進(jìn)入地下水中,形成超標(biāo)污染暈,呈現(xiàn)面積逐漸擴(kuò)大、濃度逐漸降低的特點(diǎn),至第30年時(shí),污染物濃度最低,污染物遷移距離最大為487 m,污染面積達(dá)137 583 m2。在非正常狀況及事故情景下,超標(biāo)污染暈呈現(xiàn)逐步擴(kuò)大的趨勢,運(yùn)移出廠界,項(xiàng)目停止運(yùn)行后,超標(biāo)污染暈有可能繼續(xù)向東南擴(kuò)散。

建議在污染裝置下布設(shè)防滲措施,并在其下游布設(shè)監(jiān)測井和應(yīng)急抽水井,防止地下水污染物對(duì)園區(qū)外地下水環(huán)境造成的影響。針對(duì)項(xiàng)目可能發(fā)生的地下水污染情況,建議園區(qū)進(jìn)行“可視化”處理,污水輸送管道盡可能地上敷設(shè),減少埋地管道;擬建項(xiàng)目以水平防滲為主,防滲設(shè)計(jì)根據(jù)項(xiàng)目場地天然包氣帶防污性能、污染物控制難易程度和污染物特性,對(duì)擬建項(xiàng)目采取整體分區(qū)防滲。同時(shí),按照規(guī)范要求做好地下水長期監(jiān)測。

[1] 環(huán)境保護(hù)部.關(guān)于印發(fā)《全國地下水污染防治規(guī)劃(2011—2020年)》的通知:環(huán)發(fā)[2011]128號(hào)[A/OL].(2005-08-20)[2017-09-20].http://www.360doc.com/content/15/0820/19/17103623_493700722.shtml.

[2] 陸浩.軍用鉛酸動(dòng)力電池建設(shè)項(xiàng)目水文地質(zhì)勘查報(bào)告[R].武漢:武漢地質(zhì)工程勘察院,2011.

[3] 肖威,黃瑞金,王銀國,等.黃岡化工園控制性詳細(xì)規(guī)劃[R].鄂州:鄂州市城市規(guī)劃勘測設(shè)計(jì)研究院,2013.

[4] 耿三方,徐月珍.何謂水文地質(zhì)概念模型[J].水文地質(zhì)工程地質(zhì),1985(4):38.

Based on the Numerical Simulation of GMS in a Chemical Industry GardenApplication of Water Environmental Impact Assessment

WU Pengfei, PENG Zhan, CHEN Xiaoting

(ThirdGeologicalBrigadeofHubeiGeologicalBureau,Huanggang,Hubei438000)

In groundwater environment impact assessment in construction project,by analyzing hydrogeology survey is carried out on the field,by using the numerical simulation analysis based on GMS three-dimensional numerical model,is the ranges of polluted groundwater pollutant migration rule,determine the evaluation and the concentration distribution of the important mean.Accident condition through numerical simulation analysis,research and under the condition of abnormal condition,the characteristics of groundwater pollutant migration rule,and its concentration in groundwater environmental change,for a certain time of inland water may be contaminated situation analysis and prediction.

the impact of groundwater environment; GMS; numerical simulation; chemical industry area

X523

A

1671-1211(2017)06-0728-07

10.16536/j.cnki.issn.1671-1211.2017.06.012

2017-08-21;改回日期2017-09-20

吳鵬飛(1988-),男,工程師,水文與水資源工程專業(yè),從事水工環(huán)及地質(zhì)環(huán)境調(diào)查評(píng)價(jià)工作。E-mail:260221964@qq.com

數(shù)字出版網(wǎng)址:http://www.cnki.net/kcms/detail/42.1736.P.20171026.0845.034.html數(shù)字出版日期2017-10-26 08:45

于繼紅)