水環(huán)境數(shù)學(xué)模型

1、.,水環(huán)境數(shù)學(xué)模型 薛聯(lián)青 2011年10月,.,Simulation與仿真 建立水環(huán)境數(shù)學(xué)模型的目的 了解水環(huán)境系統(tǒng)內(nèi)部因子變化規(guī)律 對(duì)水環(huán)境系統(tǒng)變化進(jìn)行定性定量描述 規(guī)劃、管理、決策需要 水環(huán)境系統(tǒng)模擬任務(wù) 量化、優(yōu)化、決策、控制,水環(huán)境系統(tǒng)建模,.,水環(huán)境模擬涉及主要問(wèn)題,水流運(yùn)動(dòng) 污染物在水中的遷移轉(zhuǎn)化 水體的耗氧和復(fù)氧過(guò)程 河流水質(zhì)模型 湖泊與水庫(kù)水質(zhì)模型 面源污染分析 水污染控制系統(tǒng)規(guī)劃,.,水環(huán)境模擬模型,確定性模型模擬 不確定性模擬 隨機(jī)方法 概率統(tǒng)計(jì)方法 灰色模型 人工神經(jīng)網(wǎng)絡(luò)法 馬爾柯夫法 自組織法 多元回歸等,.,水環(huán)境系統(tǒng)模擬及污染控制發(fā)展,綜合水質(zhì)模型的完善 基于地理

2、信息系統(tǒng)平臺(tái)的研究 模擬預(yù)測(cè)的不確定性（敏感性）分析 基于可視化技術(shù)和技術(shù)的研究 水環(huán)境模擬及修復(fù)技術(shù),.,第一講 緒論,水環(huán)境分析 水文循環(huán)過(guò)程中水的污染和自?xún)?水體污染物及水體功能、水質(zhì)標(biāo)準(zhǔn) 水環(huán)境質(zhì)量的度量與評(píng)價(jià),.,定義：,水環(huán)境模擬及分析：是在研究河流、湖泊、水庫(kù)、海洋等水體的水質(zhì)變化機(jī)理和規(guī)律基礎(chǔ)上，建立水環(huán)境模擬預(yù)測(cè)等模型，根據(jù)將來(lái)的排污、水文氣象等條件，對(duì)未來(lái)水環(huán)境狀況進(jìn)行預(yù)測(cè)分析。,.,水環(huán)境分析內(nèi)容 水環(huán)境：是自然環(huán)境的一個(gè)重要組成部分，指自然界各類(lèi)水體，如河流、湖泊、水庫(kù)、海洋、地下水、空中水等的數(shù)量、質(zhì)量狀態(tài)的總和； 水量：降水、蒸發(fā)、下滲、徑流的變化； 水質(zhì)：泥沙、水

3、溫、溶解氧、有機(jī)物、無(wú)機(jī)物、重金屬、水生生物等； 水環(huán)境：水量與水質(zhì)的統(tǒng)一；,.,分析過(guò)程,針對(duì)要解決的水環(huán)境問(wèn)題，收集有關(guān)的水文、氣象、水質(zhì)觀測(cè)、實(shí)驗(yàn)資料和污染負(fù)荷情況； 根據(jù)被模擬水質(zhì)的物理、化學(xué)、生物變化規(guī)律，建立反映模擬物質(zhì)與其它因素間相互聯(lián)系的模型結(jié)構(gòu)； 率定模型參數(shù)； 模型檢驗(yàn)；,.,水文循環(huán)中水的污染與自?xún)?水循環(huán)：水的三態(tài)轉(zhuǎn)換； 自?xún)舻娜^(guò)程：物理、化學(xué)、生物過(guò)程； 物理凈化：污染物在水體中混合、稀釋、沉淀、吸附、凝聚、向大氣揮發(fā)和病菌死亡等物理作用過(guò)程； 化學(xué)凈化：污染物在水中由于分解化合、氧化還原、酸堿反應(yīng)等化學(xué)作用下濃度降低或喪失毒性等現(xiàn)象。 生物凈化：水體微生物群，在

4、分泌的酶作用下，使污染物分解和轉(zhuǎn)化為無(wú)害物質(zhì)的現(xiàn)象。,.,水體污染物,好氧有機(jī)物 可溶性鹽類(lèi)和酸、堿物質(zhì) 重金屬污染 有毒化學(xué)品 懸浮固體 油類(lèi)污染 熱污染 放射性污染 病源微生物污染,.,水功能區(qū)劃及納污能力計(jì)算,水域：國(guó)家自然保護(hù)區(qū)、生活飲用水、水源保護(hù)區(qū)、漁類(lèi)保護(hù)區(qū)、灌溉水源區(qū)等。 水體功能與水質(zhì)標(biāo)準(zhǔn)；如地面水環(huán)境質(zhì)量標(biāo)準(zhǔn) 其中: I 主要適用于源頭區(qū)，國(guó)家自然保護(hù)區(qū)； II 集中式生活飲用水水源地一級(jí)保護(hù)區(qū)， 珍貴魚(yú)類(lèi) 保護(hù)區(qū)等； III 集中式生活飲用水水源地二級(jí)保護(hù)區(qū)，一般魚(yú)類(lèi)保護(hù)區(qū)及游泳區(qū)； IV 一般工業(yè)用水區(qū)及人體非接觸的娛樂(lè)用水區(qū)； V 農(nóng)業(yè)用水區(qū)和一般景觀要求水域；,.,

5、第二講 數(shù)學(xué)模型概述,定義與分類(lèi) 數(shù)學(xué)模型的建立 模型的參數(shù)估值 模型的檢驗(yàn)與誤差分析 靈敏度分析,.,定義與分類(lèi),定義：根據(jù)觀察到的現(xiàn)象，歸結(jié)成一套反映其數(shù)量關(guān)系的數(shù)學(xué)公式與具體算法，用以描述對(duì)象的運(yùn)動(dòng)規(guī)律。 特征：抽象性與局限性 分類(lèi)：動(dòng)態(tài)模型和穩(wěn)態(tài)模型、線性與非線性模型、確定性模型與隨機(jī)模型、模擬模型和規(guī)律模型、參數(shù)模型和分布模型,.,數(shù)學(xué)模型的建立,對(duì)模型的要求 足夠的精度 可操作、實(shí)用 依據(jù)充分 存在可控變量 建模過(guò)程 數(shù)據(jù)收集與分析模型結(jié)構(gòu)選擇： 白箱模型、黑箱模型、灰箱模型 參數(shù)估值 模型檢驗(yàn)與修正 模型應(yīng)用于反饋,.,參數(shù)估值,圖解法：適用于線性關(guān)系 y= a+bx 一元線性回

6、歸 假設(shè)條件 自變量沒(méi)有誤差，因變量存在測(cè)量誤差 各測(cè)量點(diǎn)擬合最好的直線，為各點(diǎn)至直線的因變量偏差的平方和最小的直線，即,.,為了使偏差的平方和最小，必須滿足： 于是得到：,.,多元線性回歸：對(duì)于自變量的數(shù)目大于等于2的線性模型，可以采用多元線性回歸方法求解。,.,上式中：,.,最優(yōu)化方法：原理與線性回歸方法類(lèi)似 網(wǎng)格法：在可以預(yù)先估計(jì)參數(shù)區(qū)間的情況下，將各個(gè)參數(shù)的區(qū)間等分，在所有頂點(diǎn)處計(jì)算目標(biāo)值，并比較目標(biāo)值的大小，選優(yōu)。 經(jīng)驗(yàn)公式法,.,模型的檢驗(yàn),圖形表示法：如果測(cè)量值與計(jì)算值的交點(diǎn)位于 45o線附近一定范圍內(nèi)，則可以認(rèn)為模型的模擬結(jié)果是合格的該方法多用于模型計(jì)算誤差較大的場(chǎng)合。 相關(guān)系

7、數(shù)法：用相關(guān)系數(shù)來(lái)衡量曲線的擬合程度，適用于線性程度高的模型。,.,式中： 分別表示實(shí)測(cè)值和實(shí)測(cè)值的平均值； 分別表示計(jì)算值和計(jì)算值的平均值。 r 在 0 到 1 之間， r 值越大，擬合程度越高。 相對(duì)誤差法 式中， yi 為實(shí)測(cè)值， yi,為對(duì)應(yīng)的計(jì)算值,.,靈敏度分析,靈敏度分析的意義 估算模型計(jì)算結(jié)果的偏差 有利于根據(jù)需要探討建立高靈敏度或低靈敏度的模型 可以用來(lái)確定合理的設(shè)計(jì)裕量 環(huán)境系統(tǒng)的兩種靈敏度分析 狀態(tài)與目標(biāo)對(duì)參數(shù)的靈敏度，即研究參數(shù)變化對(duì)狀態(tài)變量和目標(biāo)產(chǎn)生的影響。 目標(biāo)對(duì)狀態(tài)的靈敏度，即研究狀態(tài)變量的變化對(duì)目標(biāo)值產(chǎn)生的影響。,.,狀態(tài)與目標(biāo)對(duì)參數(shù)的靈敏度 定義：在0 附近，

8、狀態(tài)變量 x （或目標(biāo) Z ）相對(duì)于原值 x*(或 Z*)的變化率和參數(shù)相對(duì)于0 的變化率的比值 狀態(tài)對(duì)參數(shù)的靈敏度： 目標(biāo)對(duì)參數(shù)的靈敏度,.,當(dāng) 0 時(shí)，可忽略高階微分項(xiàng)，得： 式中： 和 分別叫做狀態(tài)變量 和目標(biāo)函數(shù)的一階靈敏度系數(shù)，它反映了系統(tǒng)的靈敏度特征。,.,例：BOD降解規(guī)律為： ，若已知起點(diǎn) BOD5 濃度 L0 15mg/L ，BOD衰減速度常數(shù) kd=0.1d-1 ， kd 的變化幅度在 10 ，試求 t=2 d 處的 BOD5 值及其變化幅度。 解：t=2d處的BOD5為： BOD對(duì) kd 的靈敏度為：,.,BOD的變化幅度： 由 kd 的不確定性引起的BOD變化值：,.,

9、污染物在環(huán)境介質(zhì)中的運(yùn)動(dòng)特征 基本模型的建立 非穩(wěn)定源排放的解析求解 基本模型的穩(wěn)態(tài)解 環(huán)境質(zhì)量模型的數(shù)值求解,第三講 水環(huán)境模擬模型,.,水質(zhì)數(shù)學(xué)模型:是根據(jù)排入水體的污染物,分析預(yù)測(cè)未來(lái)水質(zhì)狀況的一種數(shù)學(xué)手段和工具.應(yīng)能全面準(zhǔn)確地反映污染物在水中的遷移轉(zhuǎn)化規(guī)律.（各種過(guò)程本身的特性是水質(zhì)分析和建模的基礎(chǔ)） 水環(huán)境系統(tǒng)數(shù)值模擬模型 確定性模型 隨機(jī)模型,水環(huán)境分析模型建立,.,污染物在水中的物理遷移過(guò)程：主要包括污染物隨水流的推移與混合，受泥沙顆粒和底岸的吸附與解吸、沉淀與再懸浮，底泥中污染物的輸送等作用過(guò)程。 水中有機(jī)污染物降解與轉(zhuǎn)化 污水生化反應(yīng)動(dòng)力學(xué),污染物在水中的遷移轉(zhuǎn)化特征,確定性

10、模擬模型：,.,遷移擴(kuò)散：污染物在水流作用下產(chǎn)生的轉(zhuǎn)移作用。包括：對(duì)流、分子擴(kuò)散、紊動(dòng)擴(kuò)散、離散 對(duì)流遷移通量的計(jì)算 式中： fx ， fy ， fz 分別為 x ， y ， z 方向上的污染物對(duì)流遷移通量； ux ， uy ， uz 環(huán)境介質(zhì)在 x ， y ， z 方向上的時(shí)均流速分量； C 是污染物在環(huán)境介質(zhì)中的時(shí)均濃度。,污染物在水中的遷移過(guò)程,.,過(guò)水?dāng)嗝嫖廴疚镙斠坡?斷面A上污染物輸移率為斷面平均流速和平均濃度及斷面面積乘積。,.,擴(kuò)散是由于物理量在空間上存在梯度使之在空間上趨于均化的物質(zhì)遷移現(xiàn)象。 分子擴(kuò)散：水中污染物由于分子的無(wú)規(guī)則運(yùn)動(dòng)，從高濃度區(qū)向低濃度區(qū)的運(yùn)動(dòng)過(guò)程。 Fick

11、 第一定律：分子擴(kuò)散質(zhì)量通量與擴(kuò)散物質(zhì)的濃度梯度成正比。 式中: I 分別表示 x ，y ，z方向上的污染物擴(kuò)散通量； Em 為分子擴(kuò)散系數(shù)m2/s，C是時(shí)均濃度。,分子擴(kuò)散作用輸移,.,湍流擴(kuò)散：湍流流場(chǎng)中質(zhì)點(diǎn)的瞬時(shí)值相對(duì)于平均值的隨機(jī)脈動(dòng)導(dǎo)致的分散現(xiàn)象。 式中: I分別表示 x ， y ， z 方向上由湍流擴(kuò)散引起的污染物擴(kuò)散通量； Ex ,Ey, Ez 為紊動(dòng)擴(kuò)散吸系數(shù)m2/s ；C為環(huán)境介質(zhì)中的污染物的時(shí)間平均濃度。,紊動(dòng)擴(kuò)散作用輸移,.,彌散輸移：為了補(bǔ)償由于采用狀態(tài)的空間平均值描述實(shí)際的空間分布不均所產(chǎn)生的輸移。 式中， I表示 x ， y ， z 方向上由湍流擴(kuò)散引起的污染物擴(kuò)散

12、通量； D 為離散系數(shù)； 為環(huán)境介質(zhì)中的污染物的時(shí)間平均濃度。,離散（彌散）作用輸移,.,湍流擴(kuò)散和彌散的引進(jìn)是為了彌補(bǔ)在實(shí)際計(jì)算中采用時(shí)間和空間平均值而引起的誤差。,.,廢水在河流中的混合,由于移流、擴(kuò)散、離散作用的存在，廢水排入河流后在河流中一般出現(xiàn)三種不同混合狀態(tài)的區(qū)段。 豎向混合河段：沿垂直方向達(dá)到混合均勻（三維） 橫向混合河段：從豎向均勻混合到下游污染物在整個(gè)橫斷面上均勻混合的區(qū)段（二維） 縱向混合河段：橫向混合均勻河段之后的河段（一維）,.,.,費(fèi)希爾（H.B.Fischer)公式,按有邊界限制水流中污染源對(duì)流擴(kuò)散公式； 斷面最小濃度和最大濃度之差在5%以?xún)?nèi)作為達(dá)到完全混合的標(biāo)準(zhǔn)；

13、 估算順直河流中達(dá)到斷面完全混合的距離的計(jì)算公式： 河流中心排污： 岸邊排污： L-排污口到斷面完全混合的距離 U-河流斷面平均流速；Ey-橫向擴(kuò)散系數(shù),.,吸附與解吸,吸附：水中溶解的污染物或膠狀物,當(dāng)與懸浮于水中的泥沙等固相物質(zhì)接觸時(shí),將被吸附在泥沙表面,并在適宜的條件下隨泥沙一起沉入水底,使水的污染物濃度降低,起到凈化作用; 解吸：被吸附的污染物質(zhì)當(dāng)水體條件(流速、濃度、PH）改變時(shí)，又溶于水中的過(guò)程。 吸附-解吸作用總的趨勢(shì)：水體污染濃度減少,.,吸附作用,一是弗勞德利希（Freundlich)吸附等溫式; 二是海納利(Henery)吸附等溫式；,Freundlich Henery,S

14、e:吸附達(dá)到平衡時(shí)水中泥沙的吸附濃度，等于泥沙吸附的污染物總量除以泥沙總量。 Ce：吸附平衡時(shí)水體的污染濃度，k,n為經(jīng)驗(yàn)常數(shù),.,沉淀與再懸浮計(jì)算,一、河流動(dòng)力學(xué)原理：先計(jì)算河段含沙量變化過(guò)程和沖淤過(guò)程，然后考慮泥沙對(duì)污染物的吸附解析作用，進(jìn)一步計(jì)算出污染物的沉淀與再懸浮。 二、采用一個(gè)系數(shù)直接對(duì)污染成分的減少和增加進(jìn)行估算。,C水中污染物在t時(shí)的濃度； Kc沉淀與再懸浮系數(shù)，沉淀取正，再懸浮取負(fù)；,.,有機(jī)污染物的衰減和轉(zhuǎn)化 降解：有機(jī)污染物在水中遷移擴(kuò)散的同時(shí)，還有微生物的生物化學(xué)作用下分解和轉(zhuǎn)化為其它物質(zhì)，從而使水體中有機(jī)污染濃度降低的現(xiàn)象。 根據(jù)溶解氧情況：分好氧和厭氧情況；并且在好

15、氧或厭氧微生物的代謝作用下發(fā)生分解和轉(zhuǎn)化；,.,水中有機(jī)物的好氧降解轉(zhuǎn)化過(guò)程,.,水中有機(jī)物的降解轉(zhuǎn)化示意圖,.,降解轉(zhuǎn)化（生化反應(yīng)動(dòng)力學(xué)）,非守恒物質(zhì)進(jìn)入環(huán)境中以?xún)煞N方式發(fā)生降解： 由污染物自身的運(yùn)動(dòng)變化規(guī)律決定 在自然環(huán)境因素的作用下，由于化學(xué)或生物反應(yīng)而不斷衰減 式中， k 為降解速度常數(shù),.,生化反應(yīng)動(dòng)力學(xué)關(guān)系,生物降解反應(yīng)速度與有關(guān)因素，主要是與污染濃度、微生物變化關(guān)系： （1）水中微生物（菌、藻）增長(zhǎng)規(guī)律，直接影響污染物的降解； （2）水中有機(jī)污染物的降解規(guī)律，與水質(zhì)預(yù)測(cè)直接相關(guān)；,.,微生物增長(zhǎng)速度方程莫諾特方程,u微生物比增長(zhǎng)速度d-1 ，為微生物濃度增長(zhǎng)速度與當(dāng)時(shí)的微生物濃度

16、之比，即（Dx/dt)/X,X為微生物濃度 um基質(zhì)濃度較大情況時(shí)的最大比增長(zhǎng)速度 Ks半速常數(shù)，為u= um/2時(shí)的基質(zhì)濃度；,.,水體的耗氧過(guò)程和復(fù)氧過(guò)程,.,氧垂曲線,.,水體耗氧、復(fù)氧參數(shù)估值,K2，K1的處理與水動(dòng)力學(xué)因素，水文，PH值等因素有關(guān)；,.,第四講 水質(zhì)遷移轉(zhuǎn)化基本方程,基本假定 污染物與環(huán)境介質(zhì)相互溶合，污染物質(zhì)點(diǎn)與介質(zhì)質(zhì)點(diǎn)具有相同的流體力學(xué)特征。 污染物進(jìn)入介質(zhì)后能均勻地分散開(kāi)，不產(chǎn)生凝聚、沉淀和揮發(fā)，可將污染物質(zhì)點(diǎn)當(dāng)做介質(zhì)質(zhì)點(diǎn)進(jìn)行研究。,.,模型的推導(dǎo) 以函數(shù) C （ x ， y ， z ， t ）表示流體在點(diǎn) P （ x ， y，z ）處 t 時(shí)刻的污染物濃度。,

17、z,.,在 x 方向， t 時(shí)間內(nèi)，由于流體的推流遷移而造成微元體內(nèi)污染物質(zhì)的變化量為 在 y 方向和 z 方向的變化量分別為：,.,在 x 方向， t 時(shí)間內(nèi)，由于流體的擴(kuò)散作用而造成微元內(nèi)污染物質(zhì)的變化量為 在 y 方向和 z 方向的變化量分別為：,.,.,根據(jù)質(zhì)量守恒定律，可以得到： 將相應(yīng)公式代入上式，兩邊同時(shí)除以 x y z t ，并令 ，得到：,.,零維基本模型 一維基本模型 二維和三維基本模型,確定性水環(huán)境系統(tǒng)的基本模型,.,零維模型,所謂零維模型是描述在研究的空間范圍內(nèi)不產(chǎn)生環(huán)境質(zhì)量差異的模型這個(gè)空間范圍類(lèi)似于一個(gè)完全混合反應(yīng)器。,.,根據(jù)水量平衡方程可以寫(xiě)出： 由此得到零維水

18、質(zhì)遷移轉(zhuǎn)化基本方程： 式中， V 反應(yīng)器的容積； Q0 ， Q 流入與流出反應(yīng)器的物質(zhì)流量； C0 輸入反應(yīng)器的污染物濃度； C 輸出反應(yīng)器的污染物濃度，即反應(yīng)器中的污染物濃度； r 污染物的反應(yīng)速度； S 污染物的源與匯。,.,一維基本模型,一維基本模型是指描述在一個(gè)空間方向上存在環(huán)境質(zhì)量變化，即存在污染物濃度梯度的模型通過(guò)對(duì)一個(gè)微小的體積單元的質(zhì)量平衡過(guò)程的推導(dǎo)可以得到一維基本模型。,.,水流運(yùn)動(dòng)基本方程,1. 連續(xù)方程 2. 動(dòng)力方程 入流量-出流量+區(qū)間入流量=時(shí)段末蓄量-時(shí)段初蓄量（水量平衡基本方程） 根據(jù)水文氣象條件和河段地形資料，聯(lián)立求解上述方程，可得河段水位、流量、流速、水深沿

19、流程和時(shí)間的變化關(guān)系，從而作為求解水質(zhì)方程的條件給出。,.,一個(gè)微小體積元在 x 方向的污染物輸入,輸出關(guān)系：,.,根據(jù)流場(chǎng)中微小體積元（六面體）內(nèi)的輸入輸出關(guān)系，可以寫(xiě)出：,如果流場(chǎng)中的流速和彌散系數(shù)都是常數(shù)，則：,一維模型較多地應(yīng)用于比較長(zhǎng)而狹窄的河流水質(zhì)模擬。,.,穩(wěn)態(tài)一維遷移轉(zhuǎn)化方程,對(duì)于均勻河段，流量和排污穩(wěn)定時(shí)，各斷面污染濃度不隨時(shí)間變化，由此可得到具有源匯項(xiàng)的一維遷移轉(zhuǎn)化方程。,.,二維和三維基本模型,如果在 x 方向和 y 方向存在污染物的濃度梯度時(shí)，可以寫(xiě)出 x ， y 平面的二維基本模型，二維模型較多應(yīng)用于寬的河流,河口，也可用于空氣線源污染模擬。,.,二維模型的基本形式：

20、,如果在 x ， y,z 三個(gè)方向上都存在污染物濃度梯度，則可以寫(xiě)出三維空間的環(huán)境質(zhì)量基本模型：（海洋水質(zhì)模擬大多使用三維模型）,.,水質(zhì)方程的解析解,1.零維模型穩(wěn)態(tài)解 2.一維穩(wěn)態(tài)模型的解：二階線性偏微分方程,.,一維非穩(wěn)態(tài)水質(zhì)方程求解,瞬時(shí)排污情況的動(dòng)態(tài)解 拉氏變換后求解常微分方程得： 在起始斷面上，投放質(zhì)量為M的污染物質(zhì)瞬間排放于流量為Q的河水中，且污染物即刻與投放斷面的水相混合，初始時(shí)刻斷面濃度為M/Q.,.,非穩(wěn)定源排放的解析解,一維流場(chǎng)中的瞬時(shí)點(diǎn)排放源 忽略彌散，即 Dx 0,.,考慮彌散,.,瞬時(shí)點(diǎn)源排放的二維模型 其解析解如下，其中 M為污染物瞬時(shí)投放量，h為平均擴(kuò)散深度。,

21、.,瞬時(shí)點(diǎn)源排放的三維模型 其解析解為：,.,污染物在均勻流場(chǎng)中的分布特征,濃度場(chǎng)的正態(tài)分布 一維流場(chǎng)（瞬時(shí)點(diǎn)源）,上式可以寫(xiě)成：,令,.,作為典型的正態(tài)分布表達(dá)式，具有如下特征：,斷面處出現(xiàn)最大濃度的時(shí)間：,相應(yīng)的最大濃度值：,.,根據(jù)正態(tài)分布規(guī)律，在最大濃度發(fā)生點(diǎn)附近 2?t 的范圍內(nèi)，包含了大約 95 的污染物總量。,.,二維流場(chǎng)中的分布（穩(wěn)定源）,令,則有：,.,作為在 y 方向上存在正態(tài)分布的表達(dá)式，其最大濃度發(fā)生在 x 軸上，最大值為：,.,污染物到達(dá)岸邊所需的距離,定義：在中心排放的條件下，當(dāng)邊界處的污染物濃度達(dá)到斷面平均濃度的 5 ，則稱(chēng)污染物到達(dá)邊界由污染物排放點(diǎn)到污染物到達(dá)

22、邊界斷面的最小距離稱(chēng)為污染物到達(dá)岸邊所需的距離。,.,任意一個(gè)斷面的污染物平均濃度：,斷面上任意一點(diǎn)的濃度與平均濃度比值為：,.,中心排放時(shí)，y=B/2 ，可得：,根據(jù)定義,當(dāng)邊界濃度達(dá)到斷面平均濃度的 5% 時(shí)，被認(rèn)為污染物到達(dá)邊界，即：,于是可以求出： x 0.0137 相應(yīng)：,.,完成橫向混合所需的距離,定義：當(dāng)斷面上任意一點(diǎn)的污染物濃度與斷面平均濃度之比介于 0.95 至 1.05 之間時(shí)，則稱(chēng)該斷面已經(jīng)完成橫向混合由污染物排放點(diǎn)至完全混合斷面的最小距離稱(chēng)為完成橫向混合所需的距離。,.,根據(jù)斷面上任意一點(diǎn)的濃度與斷面平均濃度之間的關(guān)系，當(dāng) 時(shí)，求得 x=0.1 同時(shí)，斷面最大濃度發(fā)生在

23、 y=0 處，當(dāng) x 0.1 時(shí)，可以求得： 所以可以認(rèn)為，當(dāng) x=0.1 時(shí)，已經(jīng)完成橫向混合在中心排放時(shí)，完成橫向混合所需的距離為：,.,解析模型的應(yīng)用,環(huán)境質(zhì)量的模擬預(yù)測(cè) 解析模型的形式比較簡(jiǎn)單，應(yīng)用比較方便一維解析模型被廣泛應(yīng)用于各種中小型河流的水質(zhì)模擬，三維解析模型在空氣環(huán)境質(zhì)量預(yù)測(cè)中被普遍采用。 在流暢均勻穩(wěn)定的條件下，二維解析模型也可以用于模擬河流的水質(zhì) 在采用解析模型時(shí)一定要注意解析模型的定解條件,.,題例：,在流場(chǎng)均勻的河段中，河寬 B=500m，平均水深h=3m，流速Ux=0.5m/s ，橫向彌散系數(shù) Dy=1 m2/s 岸邊連續(xù)排放污染物，排放量 Q=1000 Kg/h

24、。試求下游 2 km 處的污染物最大濃度，污染物的橫向分布,擴(kuò)散域的寬度，以及完成橫向混合所需的時(shí)間。,.,解：已知污染物的源強(qiáng)Q=1000kg/h=277.78g/s 下游 2km處的污染物分布方差： 污染物的最大濃度發(fā)生在 y=0 處，計(jì)算如下式：,.,污染物的橫向分布可以通過(guò)計(jì)算不同的 y 值處的濃度值，然后作圖考察 完成橫向混合所需的距離 完成橫向混合所需的時(shí)間,.,估計(jì)彌散系數(shù),對(duì)于一維瞬時(shí)投放在投放點(diǎn)下游某處測(cè)得一組時(shí)間 ti 濃度 Ci 過(guò)程數(shù)據(jù)。將模型的解析解改寫(xiě)成： 對(duì)等式兩邊取對(duì)數(shù)：,.,在直角坐標(biāo)系上對(duì)下列兩個(gè)變量作圖，得到的直線斜率即為（ 1/Dx ）：,.,.,(2)

25、 矩法求解Dx、Dy,對(duì)于函數(shù)y=f (x)，可以寫(xiě)出： 零階矩（表示污染物的排放總量） 一階矩（表示污染物重心的位置）,.,二階矩（表示污染物分布的方差） 三階矩（表示分布曲線的對(duì)稱(chēng)程度）,.,二階矩M2表示分布的方差，對(duì)于一維瞬 時(shí)排放， 同時(shí)，由于： 可以得到： 對(duì)于二維穩(wěn)態(tài)模型：,.,例：在一維河流中瞬時(shí)投放若丹明染料若干，在下游 8km處測(cè)得若丹明的濃度過(guò)程線如下表所示試用矩法求河流的縱向彌散系數(shù)Dx。 ti(h)4.04.14.24.34.44.54.64.74.84.95.0 Ci 0.29 29 810 6690 18000 17000 6100 870 53 1.4 0.01

26、8,.,習(xí)題,例：在一河流岸邊排放口下游 1.5km處測(cè)量半江的欺騙橫向濃度分布，得到如下數(shù)據(jù)： Yi(m)10 20 30 40 50 70 100 150 200 300 Ci(毫克/L)35.0 31.2 28.3 20.5 14.5 7.6 1.05 0.02 00 已知河流平均流速u(mài)x=1.0m/s ，流場(chǎng)在觀察時(shí)間內(nèi)是穩(wěn)定的，降解可以忽略，試用圖解法求解河段的橫向彌散系數(shù)。,.,環(huán)境質(zhì)量模型的數(shù)值解,有限差分法：將一個(gè)空間和時(shí)間連續(xù)的系統(tǒng)變成一個(gè)離散系統(tǒng)，形成空間和時(shí)間的網(wǎng)格體系，然后計(jì)算各個(gè)網(wǎng)格節(jié)點(diǎn)處的系統(tǒng)狀態(tài)值，用以代表節(jié)點(diǎn)附近的值。其核心是用一個(gè)差分方程（差分商）來(lái)近似代替相

27、應(yīng)的微分方程（微分商）。,.,.,有限差分法的步驟 空間坐標(biāo)和時(shí)間坐標(biāo)的離散 以x 為步長(zhǎng)，把 x 方向坐標(biāo)劃分為 n 等分，每個(gè)節(jié)點(diǎn)坐標(biāo)為 xi i x （ i=1,2 n ） 第一以 t 為步長(zhǎng)，把 t 方向坐標(biāo)劃分為 m 等分，每個(gè)節(jié)點(diǎn)坐標(biāo)為 ti i t （ i=1,2 m ） 任一點(diǎn)的濃度 C( xi ， tj) Cij,.,用差分商代替導(dǎo)數(shù) 向前差分,.,向后差分,.,中心差分,.,用差分方程代替微分方程，得到一線性方程組解的穩(wěn)定性及收斂性檢驗(yàn): 收斂性：數(shù)值解是否收斂于原微分方程的真解; 穩(wěn)定性：解的過(guò)程中引入的舍入誤差是否會(huì)逐漸消失，即保持有界; 顯式差分方程 隱式差分方程,.

28、,設(shè)有方程 顯式差分：,.,.,得到 令 對(duì) i=1 對(duì) i=2 對(duì) i=3 到 n,.,顯式差分的穩(wěn)定條件 矩陣形式,.,隱式差分,.,最后得到隱式差分方程的一般格式,.,寫(xiě)成矩陣形式,.,.,對(duì)第 j+1 時(shí)刻的濃度空間分布，可由下式解出： 對(duì)第i=1個(gè)和i=n個(gè)方程， 是上下邊界的值，同時(shí)令,.,習(xí)題,1 、已知一組數(shù)據(jù)，試用： (1) 和 (2) 分別估計(jì) a1 ， b1 ， a2 ， b2 ，并作出模型檢驗(yàn)（相關(guān)系數(shù)法），說(shuō)明那一種模型結(jié)構(gòu)更適合上述數(shù)據(jù)。,.,2、一維穩(wěn)態(tài)河流，初始斷面污染物濃度C0=50mg/L，縱向彌散系數(shù)Dx=2.5mg/L，衰減系數(shù)k=0.2d-1 ，斷面平

29、均流速u(mài)x=0.5m/s ；試求下游 500m處在下述各種條件下的污染物濃度，并討論各種方法的計(jì)算結(jié)果的異同： （ 1 ）一般解析解； （ 2 ）忽略彌散作用時(shí)的解； （ 3 ）忽略推流作用時(shí)的解； （ 4 ）忽略衰減作用時(shí)的解。,.,第五講 河流水質(zhì)模型,基本水質(zhì)問(wèn)題 單一河段水質(zhì)模型 多河段水質(zhì)模型 河口水質(zhì)模型,.,河流的基本水質(zhì)問(wèn)題,污染物與河水的混合 豎向混合 橫向混合 縱向繼續(xù)混合過(guò)程 生物化學(xué)降解 碳BOD（ CBOD ）的降解 河流中的有機(jī)物的降解符合一級(jí)反應(yīng)動(dòng)力學(xué)規(guī)律,.,上式中，千周為含碳有機(jī)物降解速度常數(shù)，在其它條件不變的情況下，它是溫度的函數(shù)： T0 為參照溫度，通常取

30、 20 攝氏度是反應(yīng)活化能和水溫的函數(shù)，通常取 1.047 含氮BOD的降解 氮的降解動(dòng)力學(xué)：,.,該過(guò)程可用微分方程表達(dá),.,大氣復(fù)氧 復(fù)氧過(guò)程：一個(gè)流動(dòng)的水體從大氣中吸收氧氣的過(guò)程稱(chēng)為復(fù)氧過(guò)程（再曝氣過(guò)程）。 上式：一為氣體擴(kuò)散表面積， V 是水體體積。,.,溫度函數(shù) 光合作用：假定光合作用的速率與光照強(qiáng)度有關(guān)，光照強(qiáng)度可表示為時(shí)間的函數(shù)，所以有： T 為白天發(fā)生光合作用的持續(xù)時(shí)間； t 為光合作用開(kāi)始以后的時(shí)間； Pm 為光合作用產(chǎn)氧的最大速率。,.,藻類(lèi)的呼吸作用 藻類(lèi)的呼吸要消耗溶解氧，通常呼吸耗氧速度可以看做常數(shù) 底棲動(dòng)物和沉積物的耗氧,.,單一河段水質(zhì)模型,單一河段：研究河段內(nèi)的

31、流場(chǎng)保持均勻，且只有一個(gè)排放口，則該河段可被當(dāng)做單一河段。 S-P 模型 基本假定：河流中BOD的衰減復(fù)合一級(jí)反應(yīng)動(dòng)力學(xué)；反應(yīng)速度為常數(shù)；河流中的溶解氧的來(lái)源是大氣復(fù)氧。,.,基本形式,上式： L 河水中的BOD濃度，毫克/L D 河水中的氧虧值 k1 河水的BOD降解系數(shù)， 1/d k2 河水的復(fù)氧系數(shù)， 1/d,.,模型的解析解 L0 和 D0 分別為河流起始斷面的BOD和氧虧值。 臨界氧虧點(diǎn)：,.,托馬斯模型 在 S P 模型的基礎(chǔ)上，引進(jìn)沉淀作用對(duì)BOD去除的影響,.,其解析解為：,.,歐康奈爾模型 在托馬斯模型的基礎(chǔ)上，引進(jìn)含氮有機(jī)物對(duì)水質(zhì)的影響：,.,該模型的解析解為：,.,多河段

32、水質(zhì)模型,BOD-DO耦合矩陣模型,.,斷面劃分的原則： 斷面形狀變化處 支流或污水匯入處 取水口處 其它 多河段矩陣模型 根據(jù)連續(xù)性原理，寫(xiě)出每一個(gè)斷面的流量 Q 平衡關(guān)系：,.,根據(jù) S-P 模型寫(xiě)出由 i-1 斷面至 i 斷面之間的BOD衰減關(guān)系： 令 同時(shí)根據(jù)連續(xù)性方程，可得到：,.,令 可以得到：,.,該方程組可以用一個(gè)矩陣方程表達(dá)： 式中： 為 n 維向量， g1 0 L0,.,BOD-DO耦合矩陣模型 根據(jù) S-P 模型可以寫(xiě)出第 i 斷面的溶解氧濃度： 同樣，根據(jù)質(zhì)量平衡原理可得到：,.,即 令 再令,.,最后得到 該方程用矩陣形式表達(dá)： 式中 C 和 D 為 n 維矩陣。,.

33、,對(duì)于每個(gè)斷面的溶解氧，可表達(dá)為： 將BOD的表達(dá)式代入，得到耦合方程,.,令 可以得到多河段BOD-DO的耦合矩陣模型。 上式中的 U 是BOD的響應(yīng)矩陣， V 是DO的響應(yīng)矩陣。,.,已知一維河流的輸入,輸出數(shù)據(jù)如下圖所示設(shè)河流的飽和溶解氧值操作系統(tǒng) 10mg/L 。試用多河段模型模擬河流的BOD和DO。 單位： Q m3/s,L 毫克/L ， O 毫克/L ， k1 1/d ， k2 1/d ， t d,.,解： 第一步，計(jì)算 i=1 4 的 i ， i ， i, i 第二步，計(jì)算矩陣 A ， B ， C,D 的元素 第三步，計(jì)算逆矩陣，求出BOD和DO的響應(yīng)矩陣及向量 第四步，利用 U

34、,V 計(jì)算各斷面的BOD和DO濃度。,.,含支流的河流矩陣模型,假設(shè)主流含有 1,2 ， i n 個(gè)斷面，支流含有 1,2,第一 m 個(gè)斷面，在主流斷面匯入主流對(duì)支流寫(xiě)出矩陣方程，計(jì)算支流最下游斷面 m 的水質(zhì)，將支流作為污染源計(jì)入主流的矩陣方程。,.,其它水質(zhì)模型,QUALL-2 模型 模型描述的主要成分之間的作用 1、復(fù)氧作用； 2、河底生物 (包括底泥) 耗氧； 3、碳化階段降解耗氧； 4、光合作用產(chǎn)氧； 5、氨氮氧化耗氧； 6、亞硝酸鹽氮氧化耗氧： 7、碳化合物的沉淀；,.,8 、浮游植物對(duì)硝酸鹽氮的吸收； 9 、浮游植物對(duì)磷 (磷酸鹽) 的吸收； 10 、浮游植物呼吸產(chǎn)生磷（碳酸鹽）

35、； 11 、浮游植物的死亡、沉淀， 12 、浮游植物呼吸產(chǎn)生氨氮 13 、底泥釋放氨氮； 14 、氨氮轉(zhuǎn)化為亞硝酸鹽氛； 15 、亞硝酸鹽氮轉(zhuǎn)化為硝酸鹽氮 16 、底泥釋放磷。,.,.,模型的基本方程 式中， C 水質(zhì)組分的濃度；A河流斷面面積； Dx 縱向彌散系數(shù)； ux 縱向平均流速； S 水質(zhì)組分的來(lái)源（源）與消減（匯）項(xiàng)。 模型中各主要成分的基本方程,.,碳階段的生化降解 底泥耗氧 式中，底泥耗氧量； Kb 單位河段長(zhǎng)度上的底泥上浮速度。,.,氮循環(huán) 式中， N1 氨氮的濃度； N2 亞硝酸鹽氮的濃度； N3 硝酸鹽氮的濃度； 1 藻類(lèi)生物量中的氨氮分量； SN 單位河段長(zhǎng)度底泥中釋放

36、的氨氮速度； KN1 氨氮的衰減速度常數(shù)； KN2 亞硝酸鹽氮的衰減速度常數(shù)； Ab 藻類(lèi)的生物量； r 藻類(lèi)的呼吸速度常數(shù)； r 藻類(lèi)的比生長(zhǎng)速度常數(shù)。,.,藻類(lèi)生物量的增長(zhǎng) 式中， Sr 藻類(lèi)的沉淀速度； H 河流平均水深 藻類(lèi)的比增長(zhǎng)速度用下式計(jì)算： 第一式中， r ，最大最大的藻類(lèi)比增長(zhǎng)速度； P 正磷酸鹽的濃度；光照密度；河流的消光系數(shù)； KN 氮的半飽和濃度； KP 磷的半飽和濃度； KL 光線的半飽和系數(shù)。,.,磷循環(huán)： 式中， P 正磷酸鹽的濃度； 2 藻類(lèi)生物量中磷的分量； SP 單位長(zhǎng)度河底磷的懸浮速度。 溶解氧：,.,大腸菌的衰減： 式中： F 河流中的大腸菌濃度（個(gè)/升

37、） Kf 大腸菌死亡速度常數(shù) 其它可降解物質(zhì)： 式中，C 和 K 為任意可降解物質(zhì)的濃度與相應(yīng)的降解速度常數(shù)。,.,重金屬水質(zhì)模型 河流中重金屬的遷移轉(zhuǎn)化包括： 水體中重金屬污染物通過(guò)水面向空氣中的擴(kuò)散過(guò)程； 底部沉積物中已吸附的重金屬向水中釋放而重新進(jìn)入水體的過(guò)程； 吸附于懸浮物和沉積物后向固相遷移過(guò)程； 水體中懸浮物吸附重金屬污染物后向底層沉降過(guò)程； 水體中懸浮物吸附重金屬污染物后的絮凝過(guò)程； 重金屬污染物在水體中的擴(kuò)散遷移過(guò)程。,.,式中：C河流中溶解態(tài)的重金屬濃度；河流中的懸浮物濃度；CP河流中懸浮態(tài)的重金屬濃度；CS懸浮物中的重金屬含量；KS懸浮物吸附重金屬的速度常數(shù)；Kd重金屬在懸

38、浮物和水中的分配系數(shù)；b底泥懸浮物的懸浮速度。,.,河口水質(zhì)模型,河口是指入海河流受到潮汐作用的一段水體，它表現(xiàn)出明顯的時(shí)變特征。 一維穩(wěn)態(tài)模型 該模型在河口斷面面積定常，淡水流量穩(wěn)定的情況下，可以得到解析解：,.,排放口上游（ x0 ）：,.,C0 是在 x=0 處（排放口）的污染物濃度，可以用下式計(jì)算： 上式中的 W 為單位時(shí)間內(nèi)排放的污染物總量； Q 為淡水的平均流量； Dx 是縱向彌撒系數(shù)。,.,河口的有限段模型基本方程： 把空間坐標(biāo)離散化，每一個(gè)小的有限段則是一個(gè)完全混合反應(yīng)器,.,河口模型 第 i 個(gè)河段的推流遷移量為： 由于彌散作用導(dǎo)致的第 i 河段的質(zhì)量變化為 式中： 河段內(nèi)衰

39、減量為,.,對(duì)每一個(gè)河段可以寫(xiě)出COD的質(zhì)量平衡關(guān)系方程： 式中的 WiL 為系統(tǒng)外輸入到第 i 河段的COD量。 如果以 Di 表示第 i 河段的氧虧濃度，同樣可以寫(xiě)出每一個(gè)河段的氧虧平衡關(guān)系的方程,.,式中的 WiL 為系統(tǒng)外輸入到第 i 河段的氧虧量。 對(duì)于潮周平均狀態(tài)，可以作為穩(wěn)態(tài)問(wèn)題處理，即：,.,對(duì)河口的BOD分布，即可以寫(xiě)出矩陣方程： 式中的 G 為 n 階系數(shù)矩陣，其中的第 i 行,第 j 列的矩陣元素 Gij 可按下式計(jì)算： 當(dāng)j=i 當(dāng)j=i-1 當(dāng)j=i+1,其余矩陣元素為零。,.,河口做模型 同樣對(duì)于河口的氧虧，也可以寫(xiě)出矩陣方程。 矩陣 H 為 n 階三對(duì)角矩陣，其元

40、素可通過(guò)下式計(jì)算： 當(dāng)j=i 當(dāng)j=i-1 當(dāng)j=i+1,.,F 為 n 階對(duì)角矩陣，其對(duì)角線上的元素為： 將BOD的矩陣表達(dá)式代入氧虧的矩陣表達(dá)式中，得： 邊界條件處理： i=1 時(shí)，出現(xiàn)了 Q0,L0 和 D0 ，在計(jì)算時(shí)可將這些已知項(xiàng)計(jì)入源項(xiàng)中。 i=n時(shí)，需要直到第n+1個(gè)河段的BOD和DO濃度及參數(shù)，對(duì)此可以將下游的濃度梯度視為 0 ，即令Ln+1= Ln ， Dn+1= Dn,.,第六講 湖泊與水庫(kù)水質(zhì)模型,湖泊水庫(kù)的水質(zhì)特征 營(yíng)養(yǎng)源與營(yíng)養(yǎng)負(fù)荷 箱式水質(zhì)模型 垂向溫度分布模型 綜合水質(zhì)模型,.,湖泊水庫(kù)的水質(zhì)特征,流速小，與河流相比湖泊和水庫(kù)中的水流處于相對(duì)靜止?fàn)顟B(tài)； 停留時(shí)間長(zhǎng)，

41、湖泊與水庫(kù)中的水流交換周期比較長(zhǎng)，屬于靜水環(huán)境； 水生生態(tài)系統(tǒng)相對(duì)比較封閉； 主要水質(zhì)問(wèn)題是富營(yíng)養(yǎng)化； 水質(zhì)的分層分布。,.,典型湖泊水溫垂向分層示意圖,.,營(yíng)養(yǎng)源和營(yíng)養(yǎng)負(fù)荷,主要營(yíng)養(yǎng)物,濕重下淡水中各種元素的含量,.,利貝希最小值定律 ( 最小量的 Liebig 法律 ) 植物生長(zhǎng)取決于外界提供給它的所需養(yǎng)料中數(shù)量最少的一種。 主要營(yíng)養(yǎng)源與營(yíng)養(yǎng)負(fù)荷計(jì)算 地表徑流的營(yíng)養(yǎng)負(fù)荷 Ijl 第 j 種營(yíng)養(yǎng)物質(zhì)的負(fù)荷， g/； Ai 第 i 種土地利用類(lèi)型的面積， m2 ； Eij 第 i 種土地利用類(lèi)型的單位面積上第 j 種污染物的流失量， g/m2 ； m 土地利用類(lèi)型的總數(shù)。,.,降水的營(yíng)養(yǎng)負(fù)荷

42、式中， Ijp 有降水輸入的第 j 種污染物的負(fù)荷， g/；如湖,庫(kù)的水面面積， m2 ； Cj 第 j 種營(yíng)養(yǎng)物在降水中的含量； P 年降水量， m/。 人為因素排放的營(yíng)養(yǎng)負(fù)荷 生活污水 式中， Ijs 流入湖泊或水庫(kù)的污水中含有的第 j 種營(yíng)養(yǎng)物的負(fù)荷， g/； S 產(chǎn)生污水的人數(shù)，人； Ejs 每人每年產(chǎn)生的第 j 種營(yíng)養(yǎng)物的量， g/人/。,.,工業(yè)污水 jk 第 k 種工業(yè)廢水中第 j 種營(yíng)養(yǎng)物的負(fù)荷， g/； Qk 第 k 種工業(yè)廢水的排放量， m3/； Ejk 第 k 種廢水中第 j 種營(yíng)養(yǎng)物的含量， g/m3 ； n 含第 j 種營(yíng)養(yǎng)物的污染源數(shù)。 內(nèi)部營(yíng)養(yǎng)負(fù)荷 ji 湖泊底泥

43、釋放的第 j 種營(yíng)養(yǎng)物的負(fù)荷， g/； Cji 底泥第 j 種營(yíng)養(yǎng)物的含量， g/m3 ； kj 底泥釋放第 j 種營(yíng)養(yǎng)物質(zhì)的速率常數(shù)。,.,湖泊、水庫(kù)的總營(yíng)養(yǎng)負(fù)荷： 式中， Ij 湖泊,水庫(kù)第 j 種污染物的總負(fù)荷。,.,湖泊水庫(kù)的箱式水質(zhì)模型,完全混合型 Vollenweider 模型 Vollenweider 模型把湖泊看作單一均勻的整體而不考慮湖泊的分層情況及湖水的對(duì)流擴(kuò)散作用。 式中： p 為湖水總磷濃度； J 為年輸入湖泊的總磷量； V 為湖泊容積；為磷沉積系數(shù)；為水力沖刷速率 (等于年輸出水量與湖泊容積 V 之比) 。,.,在假定,均不隨時(shí)間變化并給定初始條件：當(dāng) t=0 ，p=

44、p0時(shí)，得到該模型的解析解 在水體的入流、出流及營(yíng)養(yǎng)物質(zhì)的輸入穩(wěn)定的條件下， 令 ，得到穩(wěn)態(tài)方程為： 式中， L 為單位面積總磷負(fù)荷； Z 為湖水平均深度。,.,例：已知湖泊的容積 V 1.0 107m3，支流輸入水量Qin=0.5 108 m3/，河流中的COD濃度 3個(gè)毫克/L ；湖泊的COD本底濃度 C0 1.5個(gè)毫克/L ，COD在湖泊中的沉積速度常數(shù) =0.08/試求湖泊的COD平衡濃度，及達(dá)到平衡濃度的 99 所需的時(shí)間。,.,解： 代入各已知數(shù)據(jù)，得到 t 為 0.77,.,面源污染介紹,降雨徑流污染模型：以水文數(shù)學(xué)模型為基礎(chǔ) 20世紀(jì)70年代中期，是非點(diǎn)源模型發(fā)展時(shí)期； Hyd

45、rocomp公司的非點(diǎn)源系列模型： PTR、HSP、ARM、NPS 以及其它模型：SWMM、STORM、ACTMO、UTM等；,.,應(yīng)用研究,80年代后，模型應(yīng)用，開(kāi)發(fā)新的實(shí)用模型、非點(diǎn)源污染控制與管理措施階段： 特征污染單位線模型； 非點(diǎn)源污染負(fù)荷函數(shù)模型； 污染瞬時(shí)單位線模型；,.,降雨徑流污染形成過(guò)程及研究途徑,降雨徑流子過(guò)程水污染的載體； 產(chǎn)沙輸沙子過(guò)程； 污染物隨水流運(yùn)動(dòng)中的遷移轉(zhuǎn)化子過(guò)程； 最終表現(xiàn)為河流某一斷面的徑流過(guò)程和污染負(fù)荷過(guò)程； 受納水體污染子過(guò)程;,.,研究過(guò)程,根據(jù)降雨徑流污染過(guò)程的特點(diǎn)，步驟： 研究區(qū)域按地理土壤條件及土地利用類(lèi)型分類(lèi)； 在同類(lèi)型的小區(qū)域（流域）中選

46、擇典型試驗(yàn)區(qū)； 對(duì)典型流域進(jìn)行一定時(shí)間（一個(gè)水文年）的監(jiān)測(cè)； 建立代表流域的降雨徑流污染水質(zhì)模型； 應(yīng)用模型預(yù)測(cè)研究區(qū)域中要求地點(diǎn)的降雨徑流污染負(fù)荷過(guò)程；,.,降雨徑流污染監(jiān)測(cè),在整個(gè)降雨徑流過(guò)程中同步監(jiān)測(cè)降雨量、徑流量（地表、地下）和污染物濃度的連續(xù)時(shí)變過(guò)程。,.,降雨徑流污染負(fù)荷模型,降雨徑流污染形成的基礎(chǔ)降雨徑流過(guò)程； 降雨徑流污染負(fù)荷模型包括：模型結(jié)構(gòu)、模型參數(shù)率定、模型檢驗(yàn)、模擬預(yù)測(cè)； 一般包含3個(gè)子模型： 降雨徑流子模型； 流域侵蝕及泥沙子模型； 污染物遷移轉(zhuǎn)化子模型；,.,從水文學(xué)角度的模型分類(lèi),以經(jīng)驗(yàn)公式為主的污染負(fù)荷模型； 以單位線為特征的模型：流量過(guò)程和負(fù)荷過(guò)程的匯流計(jì)算均

47、采用單位線法，并分為時(shí)段單位線法和瞬時(shí)單位線法； 水量單位線流量過(guò)程； 污染單位線污染負(fù)荷過(guò)程； 以物理成因分析法為基礎(chǔ)的流域概念模型； 從物理成因的原理描述降雨徑流污染的水動(dòng)力學(xué)過(guò)程和污染物變化的物理、化學(xué)、生物過(guò)程。,.,降雨徑流污染負(fù)荷預(yù)測(cè),相關(guān)分析法： 將大的區(qū)域劃分為不同類(lèi)型的單元集水小區(qū)；對(duì)單元采用經(jīng)驗(yàn)公式計(jì)算徑流量、產(chǎn)沙量（土壤侵蝕量）、污染負(fù)荷量；集成單元結(jié)果，計(jì)算區(qū)域總污染負(fù)荷量。,.,面源污染負(fù)荷計(jì)算,劃分單元小區(qū)； 計(jì)算各單元的徑流量、土壤流失量、污染負(fù)荷量； 1、徑流量：徑流曲線數(shù)（CN）方程計(jì)算 Rs 地表徑流量；P降雨量； S 流域土壤蓄水能力； CN 徑流曲線數(shù)；

48、,.,2、土壤流失量,通用土壤流失方程（USLE） 單位面積土壤流失量，土壤侵蝕模數(shù)； K土壤可蝕性因子； 降雨能量因子，反映降雨能量對(duì)土壤侵蝕的作用，根據(jù)暴雨強(qiáng)度、雨量由綜合分析的 計(jì)算公式推求； Ls坡度長(zhǎng)度因子;C植被覆蓋因子;B侵蝕控制措施因子.,.,3、污染負(fù)荷量,某個(gè)小單元上第T天徑流輸出的污染物數(shù)量為： LDt0.1CDtRS,tTD LSt0.001CStMStTS LDt 、LSt單位面積上，某種溶解態(tài)污染物、固態(tài)物第t天的 流出量；CDt、CSt第t天的溶解態(tài)污染物、固態(tài)物濃度； RS,t第t天的地表徑流量；MSt第t天的土壤流失量；TD溶解態(tài)污染物沿地表向流域出口輸移的比

49、例系數(shù)；TS固態(tài)物沿地表向流域出口輸移的比例系數(shù)；,.,流域污染負(fù)荷量,將流域中各單元區(qū)第t天的某種溶解態(tài)的污染物相加，得到全流域第t天該污染物的負(fù)荷量。,.,降雨徑流污染負(fù)荷預(yù)測(cè)的單位線法,由流域產(chǎn)污過(guò)程（相當(dāng)于地面凈雨過(guò)程）和負(fù)荷過(guò)程（相當(dāng)于地面徑流過(guò)程）推求污染負(fù)荷單位線，根據(jù)預(yù)測(cè)的產(chǎn)污過(guò)程預(yù)報(bào)污染負(fù)荷過(guò)程。 時(shí)段特征污染單位線（CPG)； 地面（地下）徑流過(guò)程 地面（地下）徑流污染負(fù)荷過(guò)程,.,瞬時(shí)污染單位線法,基于納希（Nash)瞬時(shí)單位線IUH法，把流域?qū)Φ孛鎯粲旰臀廴疚锏膮R集過(guò)程，簡(jiǎn)化成一系列串聯(lián)的線性水庫(kù)對(duì)凈雨、產(chǎn)污過(guò)程的調(diào)蓄和混合轉(zhuǎn)化結(jié)果。,.,面源污染的流域數(shù)學(xué)物理模型,A

50、RM和SWMM模型（農(nóng)業(yè)徑流模型 ，城市暴雨水管理模型） ARM主要模擬流域的水文響應(yīng)、產(chǎn)沙、農(nóng)藥吸附與解析、農(nóng)藥降解及營(yíng)養(yǎng)物的轉(zhuǎn)化；土壤中分為四層（表土層、上層、下層和地下水層）模擬農(nóng)藥和營(yíng)養(yǎng)物在垂向上的遷移轉(zhuǎn)化； 由降雨蒸發(fā)得到植物截留、地面徑流、壤中流、下滲和地下水層的滲透。（融雪過(guò)程）,.,流域產(chǎn)沙、輸沙子模型：包括雨滴侵蝕和坡面流輸沙兩個(gè)過(guò)程。 在以上水文響應(yīng)和產(chǎn)沙基礎(chǔ)上，進(jìn)行農(nóng)藥流失量計(jì)算（包括在泥沙表面的吸附和解析過(guò)程）ARM模型采用下列過(guò)程模擬：X單位重土壤吸附的農(nóng)藥量；M永久固定態(tài)吸附的農(nóng)藥量；C溶液中農(nóng)藥的平衡濃度；N指數(shù)；K系數(shù)；,.,ARM模型,對(duì)土層中污染物隨水分的橫

51、向和豎向輸移做了模擬,但是對(duì)地下水的輸移作用沒(méi)有考慮； ARM模型可進(jìn)行單場(chǎng)暴雨和連續(xù)過(guò)程的模擬預(yù)測(cè)。,.,SWMM模型,是一個(gè)比較完善的城市暴雨水的徑流水質(zhì)預(yù)測(cè)和管理模型，根據(jù)降雨輸入（雨量過(guò)程線）和系統(tǒng)特征（流域、泄水、蓄水和處理）模擬一次暴雨事件的徑流水質(zhì)過(guò)程。,.,模塊構(gòu)成,徑流模塊：計(jì)算雨洪徑流、下滲及雨洪所挾帶的污染負(fù)荷； 輸送模塊和擴(kuò)充輸送模塊：把由徑流計(jì)算得到的 下水道進(jìn)水口的徑流過(guò)程和污染負(fù)荷過(guò)程作為輸入，經(jīng)過(guò)地下管網(wǎng)調(diào)蓄計(jì)算和水質(zhì)的遷移轉(zhuǎn)化計(jì)算，得到各個(gè)地點(diǎn)的水量、水質(zhì)的變化情況； 調(diào)蓄和處理模塊； 受納水體模塊；,.,第七講水污染控制系統(tǒng)規(guī)劃模型,系統(tǒng)的組成與分類(lèi)： 水污

52、染控制系統(tǒng)的組成 水污染控制系統(tǒng)的分類(lèi) 最優(yōu)規(guī)劃與方案選優(yōu),.,水污染控制系統(tǒng)的組成四部分,污染物發(fā)生系統(tǒng) 污染源是污水的發(fā)生源工業(yè)污染源和城鎮(zhèn)生活污染物是水污染的主要來(lái)源。隨著農(nóng)藥、化肥使用量的激增，農(nóng)業(yè)污染也變得日益突出。 污水收集輸送系統(tǒng) 污水收集、輸送系統(tǒng)是指將污水有污染源集中并輸送到污水處理廠的污水管道和污水提升泵站，亦指將污水有一個(gè)區(qū)域轉(zhuǎn)輸?shù)搅硗庖粋€(gè)區(qū)域的污水轉(zhuǎn)輸系統(tǒng)。,.,(3) 污水處理系統(tǒng) 污水處理系統(tǒng)是對(duì)污水進(jìn)行處理、改善水體是指的核心部分污水處理分：例行的污水一級(jí)、二級(jí)處理，氧化塘處理，土地處理等。在污水處理系統(tǒng)中，污染物的去除量是可控變量。通過(guò)通過(guò)調(diào)節(jié)污水處理程度調(diào)節(jié)污

53、染物的排放量，從而達(dá)到水污染控制目標(biāo)。 (4) 接受水體 水體是污水的最終出路，接受污水的水體包括河流、湖泊、海灣等。水體的水質(zhì)是一個(gè)地區(qū)環(huán)境質(zhì)量的一部分。,.,水污染控制方法,水污染控制方法很多早期的方法是針對(duì)每個(gè)小區(qū)的排水修建污水處理廠，控制污染物的排放量。 隨后，由于經(jīng)濟(jì)的發(fā)展和技術(shù)的進(jìn)步，有必要和有可能修建大型污水處理廠，區(qū)域性的污水處理廠日漸增多。,.,水污染控制系統(tǒng)的分類(lèi),按規(guī)劃層次分類(lèi)： 流域規(guī)劃：流域規(guī)劃的任務(wù)是在一個(gè)流域的范圍內(nèi)確定水污染控制的戰(zhàn)略目標(biāo)流域規(guī)劃的主要內(nèi)容是在流域范圍內(nèi)協(xié)調(diào)各個(gè)主要污染源（城市或區(qū)域）之間的關(guān)系，保證流域范圍內(nèi)的各個(gè)河段與支流滿足水質(zhì)要求。河流的

54、水質(zhì)要求主要取決于河流的功能。 流域規(guī)劃的結(jié)果可以作為污染源總量控制的依據(jù)，是區(qū)域規(guī)劃和流域規(guī)劃的基礎(chǔ),流域規(guī)劃是高層次規(guī)劃，需要高層次的主管部門(mén)主持和協(xié)調(diào)。,.,流域規(guī)劃示意圖,.,區(qū)域規(guī)劃 區(qū)域規(guī)劃是指流域范圍內(nèi)具有復(fù)雜污染源的城市或工業(yè)區(qū)水污染控制規(guī)劃區(qū)域規(guī)劃是在流域規(guī)劃指導(dǎo)下進(jìn)行的，其目的是將流域規(guī)劃的結(jié)果污染物排放總量分配給各個(gè)污染源，并為此制定具體的方案。 區(qū)域規(guī)劃既要滿足上層規(guī)劃流域規(guī)劃對(duì)該區(qū)域提出的限制，又要為下一層次的規(guī)劃設(shè)施規(guī)劃提供指導(dǎo)。,.,設(shè)施規(guī)劃 設(shè)施規(guī)劃的目的是按照區(qū)域規(guī)劃的結(jié)果，提出合理的污水處理設(shè)施方案，所選定的污水處理設(shè)施既要滿足污水處理效率的要求，又要使污水

55、處理的費(fèi)用最低。,.,按規(guī)劃方法分類(lèi),1. 排放口處理最優(yōu)規(guī)劃 排放口處理最優(yōu)規(guī)劃以每個(gè)小區(qū)的排放口為基礎(chǔ)，在水體水質(zhì)條件的約束下，求解各排放口的污水處理效率的最佳組合，目標(biāo)是各排放口的污水處理反映之和最低在進(jìn)行排放口處理最優(yōu)規(guī)劃時(shí)，各個(gè)污水處理廠的處理規(guī)模不變，它等于各小區(qū)收集的污水量。 （排放口處理最優(yōu)規(guī)劃又稱(chēng)水質(zhì)規(guī)劃）,.,2. 均勻處理最優(yōu)規(guī)劃 均勻處理最優(yōu)規(guī)劃的目的是在區(qū)域范圍內(nèi)尋求最佳的污水處理廠位置與規(guī)模的組合在同一的污水處理效率條件下，追求全區(qū)域的污水處理反映最低。 均勻處理最優(yōu)規(guī)劃也稱(chēng)污水處理廠群規(guī)劃問(wèn)題在某些國(guó)家或地區(qū)規(guī)定所有排入水體的污水都必須經(jīng)過(guò)二級(jí)處理（即機(jī)械處理生物

56、處理），盡管有的水體具有充裕的自?xún)裟芰?，也不允許降低污水處理程度。這就是污水均勻處理最優(yōu)規(guī)劃的基礎(chǔ)。,.,3. 區(qū)域處理最優(yōu)規(guī)劃 區(qū)域處理最優(yōu)規(guī)劃是排放口處理最優(yōu)規(guī)劃與均勻處理最優(yōu)規(guī)劃的綜合在區(qū)域處理最優(yōu)規(guī)劃中，既要尋求增加的污水處理廠位置與容量，又要尋求最佳的污水處理效率的組合。 采用區(qū)域處理最優(yōu)規(guī)劃方法既能廚房發(fā)揮污水處理習(xí)題的經(jīng)濟(jì)效能，又能合理利用水體的自?xún)裟芰^(qū)域處理最優(yōu)規(guī)劃問(wèn)題比較復(fù)雜，迄今尚未有成熟的求解方法。,.,最優(yōu)規(guī)劃與方案選優(yōu),最優(yōu)規(guī)劃的特點(diǎn)是根據(jù)污染源、水體、污水處理廠和輸水管線提供的信息，一次性求得水污染控制的最佳方案只有在資料詳盡、技術(shù)具備的情況下，才能順利求出最優(yōu)解

57、，最優(yōu)方案可以被視為理想方案。 與最優(yōu)規(guī)劃不同，方案選優(yōu)的工序是首先作出水污染控制的各種可能方案，然后對(duì)各個(gè)方案進(jìn)行水質(zhì)模擬檢驗(yàn)方案的可行性，并對(duì)方案的效益進(jìn)行分析，通過(guò)損益分析或多目標(biāo)規(guī)劃進(jìn)行方案選優(yōu)方案優(yōu)選是水污染控制規(guī)劃的實(shí)用方法。,.,規(guī)劃的依據(jù)有三點(diǎn),1. 污染控制系統(tǒng)費(fèi)用的構(gòu)成 包括：污水處理費(fèi)用與污水輸送費(fèi)用。 如果以一個(gè)地區(qū)的污水處理廠數(shù)目為變量，污水處理費(fèi)用和污水輸送費(fèi)用都可以表達(dá)為污水處理廠數(shù)量的函數(shù)。,.,隨著污水處理廠數(shù)量由大變小，即由分散處理逐步過(guò)渡到集中處理，系統(tǒng)的污水處理費(fèi)用將會(huì)由于規(guī)模經(jīng)濟(jì)效應(yīng)而明顯下降，但污水輸送的費(fèi)用將會(huì)迅速上升這種費(fèi)用的合成稱(chēng)為水污染控制系

58、統(tǒng)的全費(fèi)用。全費(fèi)用曲線上的最低點(diǎn)就是系統(tǒng)目標(biāo)的最優(yōu)點(diǎn)。,.,分析,對(duì)水污染控制費(fèi)用有著決定性影響的要素主要有下述三個(gè)方面：水體的自?xún)簦ㄍ┠芰Α⑽鬯幚砼c輸送的規(guī)模經(jīng)濟(jì)效應(yīng)和污水處理效率的經(jīng)濟(jì)效應(yīng)。,.,2. 水體的自?xún)裟芰?水體能夠同化污染物質(zhì)，保證水質(zhì)滿足某種既定功能要求的能力稱(chēng)為水體的自?xún)裟芰λw的自?xún)裟芰χ饕Q于它自身的物理、化學(xué)和生物學(xué)等方面的特性，也與對(duì)水質(zhì)的要求、與污水排放方式（如排放口的位置、集中排放還是分散排放）有關(guān)。 水體自?xún)裟芰梢员豢醋魇且环N特殊的自然資源，合理利用這一資源，可以降低污水處理的費(fèi)用，但水體自?xún)裟芰τ质且环N有限的資源，不能濫用在節(jié)省水污染控制費(fèi)用和防止水體污染之間應(yīng)該進(jìn)行協(xié)調(diào)。,.,3. 污水處理與輸送的規(guī)模經(jīng)濟(jì)效應(yīng),污水處理的費(fèi)用函數(shù)反應(yīng)了污水處理的規(guī)模、效率的經(jīng)濟(jì)特征目前，污水處理的費(fèi)用函數(shù)還只能作為經(jīng)驗(yàn)?zāi)Ｐ蛠?lái)處理。,式中， C 污水處理費(fèi)用； Q 污水處理規(guī)模：污水處理效率 K1 ， K2 ， K3,K4 費(fèi)用函數(shù)的參數(shù)。,.,在污水處理效率不變，即為常數(shù)時(shí)，上式可以寫(xiě)成： 式中， K K K4 為常數(shù)根據(jù)國(guó)內(nèi)外的研究成果，參數(shù) K2 的值在 0.70.8 之間。由于 K2 1 ，單位污水的處理費(fèi)用將隨著處理規(guī)模的增大而下降。費(fèi)用與規(guī)模的這種關(guān)系稱(chēng)為污水處理規(guī)模的經(jīng)濟(jì)效應(yīng)， K2 稱(chēng)為污水處理規(guī)模經(jīng)濟(jì)效應(yīng)指數(shù)。 污水處理規(guī)模經(jīng)濟(jì)效