馬少杰,李 茜,符云琳,程乾坤,方 皓,劉 舸
(1.廣東環境保護工程職業學院,廣東廣州 510655;2.佛山市水體新型污染物防治工程技術研究中心,廣東佛山 528216)
研究范圍主要為廣東黃茅海東部區域,上游河流邊界至銀洲湖官沖水位站,雞啼門水道上的雞啼門大橋處,虎跳門水道西炮臺水位站,外海東邊界至三灶潮位站,外海西邊界至大襟島以西鵝頭咀,外海向南至海水水質國控監測點位GD0411處,上游河道計算流量,海域范圍計算潮位。控制點見圖1。

圖1 研究范圍及水文計算點位分布圖
大范圍外海模型上邊界為西江高要、北江石角、東江博羅、流溪河及潭江上游,考慮到老鴉崗和石咀在潮流界內,故模擬過程將其流量邊界各上延35 km,外海下邊界設置在100 m 等深線處,計算范圍采用三角形非結構網格劃分。
水環境容量計算時往往選取不利水文條件,若水文條件太苛刻則容量計算結果會很小,在現有技術經濟條件下難以實施。河口一般綜合考慮上游河道來水和潮差的綜合影響,對于上游河道來水則為90%保證率最枯月流量(低水位)或近10年最枯月平均流量,對于潮汐影響則考慮保證率為10%、50%和90%潮差,近岸海域的潮位邊界一般選用不利潮差對應的歷史實測潮位過程或人工構造潮型作為設計水文條件。
綜合研究區域特點與管理需求,此研究設定四個計算情景,具體如下:情景一,上游入海河道典型豐水過程疊加海域10%頻率大潮差;情景二,上游入海河道典型平水過程疊加海域10%頻率大潮差;情景三,上游入海河道典型枯水過程疊加海域50%頻率中潮差;情景四,上游入海河道典型枯水過程疊加海域90%頻率小潮差。
3.1.1 計算方法
目前海域環境容量計算主要有整體海域水質控制法、水質控制點法和排污區允許范圍控制法三大類。此研究采用以排污區允許范圍為控制原則,利用數值模型進行試算,計算范圍采用非結構三角形網格劃分,網格邊長在700~1 000 m,靠近岸邊海域根據模型計算需要予以細化,局部細化區域網格邊長在100~300 m。
3.1.2 規劃排污區確定
排污區主要基于現狀直排口和規劃直排口設置,除1#排放口(位于高欄港口航運區)和2#排放口(位于雞啼門保留區)外,其他直排口和規劃排污口均位于上述排污口附近或依托排放。因此,按1#和2#兩個排污區進行計算。
3.1.3 排污混合區設置
排污混合區指污水排放口附近的允許超標區域,目標海域較為開闊,根據相關標準和區域管理實際,此研究排污混合區的范圍確定為Aa≤3 km2。
3.2.1 計算情景
此研究設計水文邊界條件共設置四種情景,上游入海河道考慮豐、平、枯水過程,外海邊界考慮10%、50%、90%潮差;設計水質邊界條件根據水質現狀監測和水質目標限值進行確定。各情景水文水質邊界條件見表1。

表1 水文水質邊界條件表
3.2.2 計算結果與討論
經計算,1#排污區各情景的鉛環境容量計算結果為1 088.51~1 224.51 kg/a,2#排污區各情景的計算結果為80.01~140.82 kg/a,情景一至情景四計算結果逐次變小。各情景鉛環境容量計算結果見表2。

表2 鉛環境容量計算結果表
從情景一和情景二計算結果對比來看,在10%大潮差工況下1#排污區平水過程比豐水過程環境容量計算結果小2.46%,2#排污區平水過程比豐水過程環境容量計算結果小10.94%,說明靠近河口的B 排污區環境容量計對受入海河流豐枯影響較大,遠離河口海域相對較開闊的1#排污區,受河海河流豐枯影響相對較小。從情景三和情景四計算結果對比來看,在入海河流枯水過程工況下,1#排污區小潮差比中潮差環境容量計算結果小8.31%,2#排污區小潮差比中潮差環境容量計算結果小31.51%,河口區受潮差情況影響大可能與局部環流或統計小潮差結果更小有關,后續研究中需要進一步分析可能的影響參數并做敏感性分析。為避免水文參數不確定性對水環境安全的影響,一般設置各排污區允許分配排污量時會保留10%~20%的富余度,即各排污口可利用允許排污量取計算環境容量的80%~90%。
總體來說目標區域潮差對計算結果影響大于入海河流流量變化,排污區附近河流徑流量越小、距河口越遠受入海河流豐枯影響越小受潮差影響越大。相比之下,入海河道枯水過程(90%頻率)疊加外海90%小潮差情景最為不利,較計算結果最大情景入海河道豐水過程(10%頻率)疊加10%大潮差小14.41%。因受水文地理特征等因素影響,不同水文情景容量結算結果往往差別較大,在日常工作中是否選用入海河流徑流過程和潮差雙不利情景的疊加進行計算值得商榷,應結合區域管理需要合理確定,可根據水域功能、敏感目標分布與區域發展定位適當取舍。因此,近岸海域水環境容量計算時可將不利潮型作為重點研究對象,結合規劃排污區實際和管理需要適當考慮入海河流水文參數,綜合確定設計水文條件。