巖溶管道網絡數值模擬方法研究進展

王益偉，渠光華，劉 埔，褚學偉，段先前

（1.貴州大學 資源與環境工程學院，貴陽550025；2.貴州大學 礦業工程學院，貴陽550025）

巖溶管道網絡控制著巖溶含水層中地下水流的水力特征及運移特征，對巖溶含水層中地下水流及溶質運移過程有著重要的影響［1-2］。要清晰刻畫地面以下巖溶管道網絡，需要獲得大量精確的管道網絡的空間信息［3］，然而，由于管道網絡的復雜性及非均質性，導致實際工作中是無法直接獲得充足的空間信息的。因此，采用數值模擬的方式獲得缺失的空間信息，是重現真實管道網絡的有效手段［4］。

目前，巖溶管道網絡特征的數值模擬研究，主要集中在：①巖溶管道網絡系統演化過程模擬；②巖溶管道網絡結構模擬（Structure-imitating）。

1 巖溶管道網絡系統演化過程模擬

巖溶管道網絡系統演化模擬，其目的是模擬控制碳酸鹽巖溶解的完整物理化學過程，包括模擬空隙、管道分布特征與水流及可溶巖溶解過程的關系，最終實現對整個巖溶管道系統演化的物理和化學過程的模擬［5-6］。這種方法對于清晰的理解巖溶演化過程是非常有用的，并能夠提供較為精確的關于空隙、管道的演化信息。但該種方法較難應用于實際中，很難反映出真實管道網絡系統的空間特征，原因在于重建形成當前管道系統的古氣候及地質條件是相當困難的。同時，構建這樣的模型，需要大量的數據，并且該類模型的解算過程通常是非線性的，對CPU的要求相當高。

2 巖溶管道網絡結構模擬

巖溶管道網絡結構模擬研究是一個極具挑戰的任務。原因在于巖溶管道網絡的復雜性，致使客觀上不可能存在有足夠詳細的信息去詳盡完整地反映出實際局部巖溶管道的特征［7-8］。針對這一問題，一種思路認為解決真實局部巖溶含水層的水文地質特征的預測，不如去解決理解整個巖溶含水系統的功能問題［9］。另一種思路，則采用統計學的觀點來重現管道網絡［10］，從而延伸彌補缺乏空間信息的區域。

2.1 全局式模型

基于第一種思路，全局式的模型被提出。全局式的模型，以黑箱理論為基礎，將巖溶管道網絡系統作為一個整體，利用系統的全局式反映，如利用泉的水動力及水化學特征，推測反演巖溶管道網絡特征。如，根據泉點流量系列、水化學系列及同位素數據，使用時間序列分析、主成分析等方法，建立起多種解譯模型［11-12］。Mayaud等［13］通過對泉點的時間序列分析，建立了有限元模擬模型，闡述了管道網絡對研究區域的水文地質條件的影響。長勇和吳吉春等［2］以泉流量系列為基礎，利用MODFLOW-CFP模塊構建了巖溶含水層的管道網絡模型，并分析了管道網絡對泉的流量過程的影響。盡管這些方法被廣泛的采用，但全局式的模型只揭示了系統全局的信息，大多數情況下并不能清晰地描述巖溶管道網絡的空間特征。這種模型能夠解決區域水量及流向問題，但對小尺度下局部水流的流量及流向問題是無法刻畫的，這樣的模型實際上是無法滿足需要明確知道水流運移路徑的研究，如巖溶含水層污染物運移及修復的研究。

2.2 分布式模型

基于第二種思想，一些分布式模型被提出。分布式模型能夠提供更真實的巖溶含水層中地下水的流量及水流方向的信息，但其難點及關鍵在于如何處理巖溶管道網絡的非均質性，既在空間中的變化。方法包括：地質統計的方法、分形的方法［14］、格子氣元胞自動機（lattice-gas cellular automaton，LGCA）方法［16］、偽并行遺傳算法［2］及反演模擬［3］。其中，基于統計學原理的模擬方法被廣泛使用。其余幾種方法通常會生成一個二維模型，但實際的管道網絡應該是一個三維的空間實體，同時，在這些模型中一些影響管道網絡形成的常規原則很難約束到所建的模型中。

基于地質統計學原理的模擬方法根據地質統計學原理，建立管道網絡的空間幾何特征（方向、相互連接方式、連通性等）與可野外直接測量得出的參數（裂隙、基底面、水力梯度、水流出口位置）之間的相關關系，從而實現對管道網絡空間分布特征的模擬。其中，最常用的方法是序貫指示模擬，用指示隨機變量來描述管道的空間變化［16-17］。Jaquet及Jeannin［15］使用管道的密度作為相關空間變量，采用轉向帶法法模擬了在某一區域上的管道網絡，并根據模擬的密度值隨機的確定了管道的位置。但以上這些傳統模型都基于變差函數的兩點地質統計學方法。由于兩點統計測度只描述在已知一點的值的情況下，在特定位置上特定值的可能性，不受兩點之間結構及通過復雜的路徑是否連通的影響，因此傳統的地質統計模型不足以表達具有復雜空間結構的含水介質，是沒有能力描述含水介質的連通性的［4，18］然而，對復雜的巖溶管道網絡來說，連通特性無疑是一個非常重要的特性，是衡量巖溶含水層非均質性的一個重要參數。這種特性使得含水介質即便具有相同的網絡結構，當其連通性不同時，水流及溶質在管道網絡中的運移也會發生變化。為了解決模型中含水介質連通性問題，基于多點地質統計的高維度的非均質模型被提出。這些模型的基本思想是利用空間多個點的組合模式來描述含水介質結構信息，采用訓練圖像代替變差函數來描述地質體結構信息，克服兩點地質統計不能再現目標幾何形態的不足，這類型的模型更適合于進行具有復雜結構含水介質信息的模擬。采用多點地質統計方法的關鍵在于構建合理的訓練圖像，尤其是在模擬三維的地質體時［19-21］。訓練圖像構建的方法包括基于對象、 基于過程的方法及混合方法［20-22］。盡管，多點地質統計適用于模擬復雜的地質結構，同時巖溶含水層具有復雜的管道網絡，但目前多點地質統計學方法多用在模擬裂隙介質及多孔介質［23-25］，在油田工程方面應用較多，卻很少用于對巖溶管道網絡的模擬。

2.3 非均質模型

伴隨著復雜的非均質模型的出現，連通性的測度成為一個量化非均質模擬及標度非均質化模型的重要工具。連通性測度的研究主要基于兩種理論：逾滲理論［26］及圖論［27］。

以逾滲理論為基礎所建立的連通性，需要事先知道滲流的關鍵閾值，取一個簡單的滲透系數作為連通性的判斷。但Edery等［28］的研究結果卻表明低滲透地區也可能出現較大的流量，而一些低流速的區域并不連通。基于圖論的連通性測度研究，對連通性的定義更為清晰。Renard等［4］和Tyukhova等［29］基于圖論的原理提出了一個簡單的定量網絡連通性的方法，通過確定最小路阻的路徑，最終確定裂隙含水介質的連通性。同時，在Renard等［4］的研究中將含水介質的連通性分為靜態連通性及動態連通性兩種。靜態連通性取決于含水介質中存在高滲透性的通路，依賴于介質本身的滲透性能。而動態連通性則取決于水流及物質運移的物理化學過程，如水流及溶質的溝槽流現象。由于靜態連通性不易被直接測出，但動態連通性較易測出，通常可通過野外的示蹤實驗獲得，因此，揭示靜態連通性和動態連通性之間的關系，有助于利用野外現場實驗提高數值模型對真實巖溶管道網絡的反演及再現，是管道網絡模擬研究中的一個重要的內容。靜態連通性和動態連通性之間存在復雜的關系［20-32］。如，Ronayne等［33］提出，當在巖溶管道周圍出現高滲透性的圍巖時，靜態連通性不能保證存在高的動態連通性。Tyukhova and Willmann等［34］根據管道的路阻及管道寬度定義了靜態連通性的測度，并將其與動態連通性進行對比，最終將其測度值應用于對水流及溶質運移的模擬。盡管，目前對巖溶管道網絡連通性的研究不多，但該特性的研究將有助于提高模型對非均質含水層中水流及溶質運移模擬的預測能力；同時，量化巖溶管道網絡連通性的測度，分析動態連通性與靜態連通性的關系，是將衡量巖溶含水層非均質特性引入巖溶含水層水流及物質運移模擬研究關鍵方法。

3 結語

（1）巖溶管道網絡模擬研究歷經幾十年的努力，取得了許多成果。早期研究多集中建立全局式、分布式的模型，忽略含水層內部的復雜性，隨著計算機技術及物探技術的不斷提高，要求更高精度的非均質模型將成為解決巖溶含水層由于其復雜性所導致的工程問題的有用工具。

（2）在非均質模型中，目前多基于地質統計學方法構建管道網絡模型，可更真實地描述巖溶含水層中地下水的流量及水流方向的信息，但傳統的地質統計模型多基于兩點地質統計學方法，其所構建的巖溶管道模型不足以描述巖溶含水層中管道空間結構的復雜性，缺乏描述具有強烈非均質性含水介質連通性的能力。采用基于多點地質統計方法，結合空間中多點特征，綜合反應含水介質空間中的結構信息，構建能夠反映管道網絡連通特征的管道網絡模型，將是管道網絡數值模擬的新趨勢，其將對解決預測巖溶含水層中優勢水流路徑及溶質暈體分布特征等關鍵問題有著重要意義。