河流與淺水湖泊水流數值模擬及污染物輸運規(guī)律研究.pdf

1、隨著城市化和工業(yè)化進程的加速，在我國當前面臨的各類水問題中，河流與湖泊水污染問題十分突出，已成為制約我國社會可持續(xù)發(fā)展的關鍵問題。利用數值模擬方法分析預測水流水質變化是水環(huán)境污染治理規(guī)劃中的主要方法之一，其結果可為相關決策者進行水環(huán)境管理規(guī)劃提供科學依據。雖然國內外學者在河流與湖泊水流水質模擬領域開展了大量的理論研究工作，一些理論與方法有了重要進展，但仍有許多需要深入研究的問題，以解決計算區(qū)域地形不規(guī)則、氣象條件多變、污染物指標眾多、不

2、同污染物之間關系復雜、水流-水質動力耦合密切、水流形態(tài)多樣、計算效率過低、計算失穩(wěn)等一系列難題。因此，亟待探索和研究河流與湖泊的水動力學驅動機理，并建立可靠、高效和穩(wěn)定的水動力水質耦合模型。
　　本文圍繞河流與湖泊中水流水質數值耦合模擬面臨的若干關鍵科學問題，基于迎風格式的有限體積方法，建立了適用于河流與淺水湖泊的高性能水流水質耦合數學模型，并以大東湖水系為研究對象，通過推求引水連通湖泊在自然循環(huán)與人工循環(huán)條件下各種污染物濃度對不

3、同流場的響應關系曲線，深入分析了連通湖泊群水循環(huán)對水污染遷移轉化過程的影響。相關成果可在湖泊水網連通工程實踐中推廣應用，具有良好的工程應用前景。本文的主要研究內容和創(chuàng)新包括：
　?。?）綜合考慮了不規(guī)則地形、氣象條件、水流形態(tài)、計算精度和計算效率等因素，建立了復雜計算條件下河流與湖泊二維水流數值模擬模型。為提高模型在復雜流體運動中的適應能力，引入局部Froude數，提出了一種水深-水位加權變量重構技術，可準確模擬二維水流中的緩流、

4、急流、混合流、間斷流等復雜流態(tài)。同時，深入分析了通量計算可能導致的負水深問題，創(chuàng)新性的提出了一種水量通量校正方法，有效解決了數值計算過程中產生負水深的難題，保證了計算過程水量守恒性，提高了水流數值模擬的有效性和穩(wěn)定性，為進一步揭示水體污染物的遷移規(guī)律提供準確的流場支持。
　?。?）針對傳統的一、二維水流數值聯解過程中面臨計算時空不同步的難題，提出了基于黎曼近似解的一二維水流方程的耦合聯解模式，建立了一維河道水流與二維平面水流數學耦

5、合模型。采用近似黎曼算子計算一、二維以及耦合區(qū)域的水流數值通量，不僅具有時間同步性，并且解決了兩種模型獨立求解可能遇到的空間分辨率和計算精度差異問題，實現了河湖連通條件下的水流數值模擬。
　?。?）針對流體運動中污染物輸運方程求解出現數值阻尼過大和數值振蕩等問題，建立了基于分布式水動力學的連通湖泊污染物遷移模型。模型采用能有效捕捉流態(tài)與污染物濃度間斷的HLLC型近似黎曼算子耦合求解水流和污染物對流通量，有效減少了模型計算過程中的數

6、值阻尼現象，為準確揭示不同水動力條件下污染物遷移規(guī)律提供了可靠的理論依據。
　　（4）針對大尺度計算區(qū)域的水流水質數值模擬面臨計算效率過低以及網格分辨率要求不同的難題，發(fā)展了一種網格自適應動態(tài)加密方法，該方法能根據計算區(qū)域的水位梯度特征和污染物梯度特征動態(tài)加密或稀疏網格數量，提高數值模擬空間精度的同時減少計算單元數量，實現了模型計算精度與計算效率的均衡。
　　（5）采用java編程語言，提出了一種基于CPU多線程并行計算的水

7、流水質模型求解模式。采用區(qū)域分解技術對計算區(qū)域進行分解，并引入虛擬單元對不同分解區(qū)域進行無縫銜接，保證了模型在計算過程中CPU資源的合理分配，進而通過多線程并行計算技術加速模型求解，極大幅度的提高了水流水質模型的計算效率，有效解決了實際工程中水流水質模擬面臨計算用時過長的難題。
　?。?）運用已建立的水流水質耦合數學模型，以大東湖水系為研究對象，開展水流水質模擬在水環(huán)境規(guī)劃中的應用研究。首先采用實測數據對水流水質模型進行參數率定，

8、然后研究了連通湖泊水動力學特性并對比分析大東湖在自然條件以及引水條件下流場分布特點，闡述了自然循環(huán)條件下水流變量與水質變量之間的定量關系，探求了污染物受水流運動影響的主要因子，揭示了污染物隨水流變化的遷移規(guī)律，最后以人工循環(huán)條件下引水路線、引水流量、引水時間等變量參數優(yōu)化為突破口，推求了引水連通湖泊在不同引水方案下各種污染物濃度對流場的響應關系曲線，闡明了水流運動、水體污染物遷移以及水質改善三者之間的定量關系，可為水環(huán)境工程規(guī)劃和水污染