水文地質模型怎麼做

A. 水文地質概念模型的同位素解釋

一、補給來源與機制

松嫩平原第四系孔隙潛水補給源為降水、山區和高平原側向徑流補給（圖5—31）。西部砂礫石台地為大氣降水和山區少量側向徑流，傾斜平原為砂礫石，降水入滲和大興安嶺山地河谷潛水補給，以及汛期河水補給。低平原潛水以降水補給為主，同時接受東部高平原潛水補給；烏裕爾河、雙陽河變為散流補給潛水，地下水的補給時期為近50 a來的補給。

圖5—27 松遼分水嶺邊界附近第四系承壓水的δ¹⁸O（‰）分布

圖5—28 松遼分水嶺邊界附近第四系承壓水的¹⁴C模型年齡（a）分布

圖5—29 松遼分水嶺邊界附近泰康組承壓水的¹⁴C模型年齡（a）分布

圖5—30 松遼分水嶺邊界附近大安組承壓水的¹⁴C模型年齡分布

圖5—31 松嫩平原地下水補給機制示意圖

北部訥謨爾河－科洛河地下水系統補給來源為當地降水與代表山區降水的地表水，地下水以垂直入滲補給為主，河流側滲為輔，其中當地降水補給占總補給的60%～73%。山區地表水補給占總補給的27%～40%。中部烏裕爾河－雙陽河地下水系統主要接受降水和地表水補給，地下水以垂直入滲補給為主；該系統北部高平原烏裕爾河附近補給來源為來小興安嶺山區降水通過地表河流補給，西南部低平原補給來源為當地降水。東部高平原呼蘭河－通肯河地下水系統和拉林河－阿什河地下水系統補給來源以當地降水垂直入滲補給為主。南部高平原第二松花江地下水系統補給來源為遠距離的長白山區降水，通過河流補給。北部山前傾斜平原和低平原的雅魯河－阿倫河－諾敏河地下水系統補給來源為西部山區降水，西南部山前傾斜平原和低平原霍林河－洮兒河－綽爾河地下水系統補給來源為距離較近的西部降水，地下水以垂直入滲補給為主。

第四系承壓水、新近系泰康組和大安組補給來源為西部山前和東部高平原地下徑流。第四系承壓含水層地下水多為近50 a來的補給，在吉林中部具有古補給的特徵。新近系泰康組和大安組除西部和北部地區有近期水補給外，其他地區為古補給，補給時期大於一萬年。第四系承壓含水層與上層潛水和下伏泰康組含水層之間水力聯系密切，大安組相對封閉性較好。除周邊山前平原和東部高平原附近發生現代補給外，泰康組和大安組含水層中的地下水多為古補給。

第四系承壓水通過西部山前扇形平原及北部和東北部高平原的弱透水邊界，接受側向潛水補給。在盆地北部、西部邊緣地帶，承壓含水層頂板變薄，局部缺失，且孔隙承壓水水位多低於上覆潛水水位，潛水以「天窗」形式向下補給承壓水，在嫩江與第二松花江匯合地帶，承壓水以頂托或越流形式向上補給河谷潛水。

新近系泰康組與上覆第四系承壓水和潛水有較好的水力聯系，在邊緣地帶與孔隙潛水接觸，通過強透水邊界可得到自上而下的直接或越流形式補給。在西部及西北部邊緣得到基岩裂隙水側向徑流補給。大部分地區，下伏的大安組、依安組承壓水水位高於該層承壓水水位，下伏承壓水自下而上的頂托和越流補給。該層承壓水徑流速度緩慢，由南北向中部徑流，在嫩江與松花江匯合地帶，以向上頂托或越流形式補給河谷潛水。

大安組承壓水主要接受北部、東部高平原孔隙潛水及西部丘陵山區基岩裂隙水側向徑流補給。地下水徑流緩慢，徑流方向總體上為由北向南，在嫩江與松花江匯合處，以越流形式向上補給上覆含水層。

高平原白堊系承壓水主要接受高平原孔隙潛水及丘陵山區基岩裂隙水的側向徑流補給，在局部地表出露地帶可直接或間接接受大氣降水的補給。地下水徑流比較緩慢，總體流向為北西—南東向。

二、地下水流動與邊界特徵

圖5—32 松嫩平原地下水流動模式

松嫩平原地下水流動主要存在局部地下水流和區域地下水流（圖5—32），局部地下水流系統存在於整個平原的淺部，主要形成於潛水含水層，由降水、灌溉水、地表水補給；在垂向入滲至含水層後，向排泄區流動，在某些低窪地帶溢出地表，循環深度一般為50 m，山前區可達100 m以下。區域地下水流存在於深部承壓含水層，地下水年齡老，為遠距離補給，補給來源為周邊降水和地表水的入滲補給；以側向水平徑流方式流向盆地中心運動，通過盆地中心承壓含水層向上越流補給潛水。總體來說，松嫩平原地下水流的總趨勢是由東部、北部和西部三面山區流入盆地，在第二松花江和嫩江匯流處形成松嫩盆地區域地下水的排泄中心。盆地地下水通過嫩江、第二松花江河谷第四系孔隙潛水排向松花江。

山前傾斜平原的砂礫石扇形地台地及白堊系統碎屑岩與大興安嶺山地基岩銜接，為弱透水補給邊界。高平原黃土狀土覆蓋區，白堊系碎屑岩及其風化帶直接與長白山地基岩接觸，為弱透水補給邊界。南部松遼邊界具有相對的隔水性，在天然狀態下可視為零通量邊界。

三、補給強度和地下水更新

松嫩平原地下水氚含量分布深度說明：山前平原現代地下水循環深度可達100m，平均補給強度126mm/a；中部低平原小於80m，平均補給強度60mm/a；東部和南部高平原現代地下水循環深度小於50 m，平均補給強度59 mm/a。這些地下水平均滯留時間小於50 a，更新性較好。

第四系承壓現代年齡地下水分布在林甸—齊齊哈爾以北、齊齊哈爾—泰來—白城以西的山前台地和乾安—肇州以東的高平原。這一地區為主要補給區，平均補給強度為6.2 mm/a，具有一定的更新能力，其餘地區¹⁴C模型年齡5～10 ka（BP），地下水更新性較差。

新近系泰康組現代年齡地下水分布在林甸—齊齊哈爾以北和乾安—前郭以東，現代補給僅發生在山前和接近東部高平原的局部地區，其餘地區地下水為古補給，低平原中部年齡8～15 ka（BP）；新近系大安組現代年齡的地下水分布在北部訥河、沿山前一帶以及與東部高平原相接地帶，為地下水的補給區，低平原中部年齡11～23 ka（BP）。對於該兩組含水層來說，現代補給僅發生在周邊地區，平均補給強度2.8 mm/a；低平原中部為古補給，這兩組含水層地下水流動緩慢，更新性較差。

B. 放水試驗數值模擬及預測礦井疏水量

水文地質計算分析是水文地質條件評價由定性上升至定量的過程，是定量評價含水層和隔水層水文地質性質的重要工作，也是充分利用各種勘探、試驗、檢測、監測資料深化對礦井水害條件認識的重要工作。目前通用且成熟先進的礦井水文地質計算方法是水文地質數值模擬技術。

（一）數值模擬方法的基本概念

水文地質中的數值模擬方法就是利用刻畫地下水系統空間結構和水力特徵的數學模型作為工具，以數字模擬方法為手段來定量分析、評價、預測地下水系統的水文地質條件、參數結構、行為規律及其在擾動條件下的變化與響應。

數值模擬方法較之解析法乃至其他評價方法來說，它能夠比較全面充分地刻畫含水層的內部結構特點和模擬處理比較復雜含水層系統邊界及其他一般解析方法難以處理的水文地質問題。可以說，無論多麼復雜的水文地質問題，只要能歸結為利用一組數學方程刻畫的數學問題之後，藉助於大型計算機這個現代科技手段，總可以用數值模擬方法獲得對問題的定量化解答。所以，數值模擬方法是目前水文地質計算中一種強有力的數學工具，它的推廣應用標志著水文地質條件定量計算與分析進入了新的發展階段。

（二）數值模擬基本過程

採用數值模擬方法定量模擬評價礦井水文地質條件基本上可分為六大步驟：認真分析和研究礦區地質與水文地質條件，在礦井水文地質條件分析的基礎上建立模擬計算域的水文地質概念模型；根據水文地質概念模型及其礦井採掘條件建立計算域的數學模型；根據模擬計算區域的水文地質結構特點採用合理的方法離散化模擬計算區域；依據模擬計算區及其相鄰區域的水文地質試驗資料或水文地質長期觀測資料校正（識別）計算區域的水文地質參數，以獲得礦井水文地質條件的預測預報數學模型；利用未參與水文地質參數識別的水文地質試驗或其他觀測資料驗證（檢驗）所建立的礦井水文地質條件預測預報數學模型；運行所建立的礦井水文地質條件預測預報數學模型進行礦井涌水量及其他水文地質條件的預測預報（模型運轉）。現分別敘述如下。

1.建立模擬計算區的水文地質概念模型

在礦區水文地質調查和專門水文地質勘探的基礎上，根據對模擬計算區域內水文地質條件的認識和分析，綱要性地概化出研究計算區的水文地質概念模型。水文地質概念模型既取決於研究計算區的具體水文地質條件，但又不完全等於該區的實際水文地質條件。它是實際水文地質條件的概化和功能綱要，礦井水文地質概念模型要求明確和概化的主要內容有：

（1）概化確定模擬計算區的范圍及邊界條件

根據礦井水文地質勘探資料和礦井採掘要求，在明確了礦井主要充水含水層和模擬計算的含水層後，根據礦井對水文地質評價的要求，首先應圈出模擬計算區的范圍。一般情況下，模擬計算區最好是一個具有自身補給、徑流和排泄的獨立的天然水文地質系統，它具有自然邊界，便於較為准確地利用其客觀真實的邊界條件，避免人為劃定邊界時在資料提供上述的困難和誤差。但在實際工作中，我們所關心或劃定的模擬計算區域常常不能完全利用上述自然邊界。這時就需要充分利用水文地質調查、勘探和長期觀測資料等通過深入系統的水文地質條件分析建立人為的模擬計算邊界。

在利用含水層自然邊界有困難或在模擬計算區邊界因勘探試驗和觀測資料缺乏，不足以建立較為精確的人為邊界時，常常將已確定的計算范圍適當地向外延伸設置一層緩沖帶，緩沖帶的寬度視具體的水文地質條件和評價要求而定，一般為2～3層計算單元的寬度。緩沖帶的邊界一般以定水頭邊界或隔水邊界處理為宜。這種方法實際上就是對無限邊界的概化處理。

在計算范圍明確規定後，就要對所有邊界的水文地質性質進行詳細的研究和確定。一般情況下，只要含水層與常年有水的湖泊、河流、水庫等地表水體有直接的水力聯系時，不管是地表水排泄地下水，還是補給地下水，只要兩者之間存在密切的水力聯系，均可處理為第一類邊界條件。但是，對於自由入滲的地表水體，則必須作為第二類邊界條件處理。

（2）概化模擬計算區域內含水層的內部結構特徵

通過對含水層結構類型、埋藏條件、導儲水空隙結構及水力特徵的分析研究，確定模擬計算區內含水層類型，如要明確所研究的目標含水層是承壓含水層、潛水含水層、半承壓含水層，或是承壓潛水含水層並存，在此基礎上要對含水層的空間分布狀態進行概化。對於承壓含水層來說，主要明確含水層厚度的變化規律及其在模擬計算區內厚度的分布，對於潛水含水層來說，主要是要明確含水層底板標高的變化規律及其在模擬計算區內底板標高的分布。其結果最好通過含水層等厚線圖或含水層底板等高線圖反映出來。含水層的滲透性（導水性）概化是根據含水層的滲透系數（或導水系數）及其主滲透方向和儲水系數在空間上的變化規律，進行均質化分區。所謂含水層水文地質參數的均質化分區就是根據對所模擬研究的含水層區域內地質與水文地質條件的分析，將研究區劃分為若干個亞區域，而且認為在每個亞區內含水層水文地質參數是相等的（含水層是均勻的）。實際上，絕對均質或各向同性的岩層是不存在的，均質性劃分也只是相對的，只要含水層的水文地質參數變化不大，則可相對地在亞區內視為均質。一般情況下，鬆散岩層中的孔隙含水層多屬於非均質各向同性，基岩裂隙或岩溶裂隙含水層則多屬於非均質各向異性含水層。

（3）概化模擬計算目標含水層的水力特徵

水力條件是驅動地下水運動的力源條件，不同的水力條件會形成不同的地下水運動形式。含水層水力特徵的概化主要包括三方面內容：一是滲流是否符合達西地下水流規律；二是含水層中的地下水流呈一維運動、平面二維運動還是空間三維運動；三是地下水水流運動是穩定流還是非穩定流。一般情況下，在鬆散沉積的孔隙含水層、構造裂隙含水層以及溶洞不大，均勻發育的裂隙岩溶含水層中，地下水流在小梯度水力驅動下多符合達西地下水流規律。只有在大溶洞和寬裂隙中的地下水在大梯度水力條件的驅動下才不符合達西水流規律。嚴格地講，在開采狀態下，地下水的運動都存在著三維流特徵，特別是在礦井排水形成區域地下水位降落漏斗附近以及大降深的疏放水井孔附近地下水的三維流特徵更加明顯。但是，在實際工作中，由於三維滲流場的水位資料難以取得，因此目前在實際模擬計算過程中，多數情況下將三維流問題按二維流近似處理，所引起的計算誤差基本上也能滿足礦井水文地質計算的要求。

（4）概化計算區域的初始水文地質條件

根據模擬計算區礦井水文地質定量評價的要求，選定模擬計算的初始時刻，求出模擬計算的初始流場（也就是計算起始的地下水流場）。模擬計算的初始條件包括計算區內的水力場，初始水文地質參數場，一類邊界的水位值，二類邊界的水力梯度值以及計算區內自然存在的地下水源、匯項。其中最常見的確定計算區內的水力場的方法就是根據區內觀測孔的水位資料，作出計算區在選定的初始時刻的等水位線圖，再根據等水位線圖最後求出所有剖分節點的水位。此外，也可通過計算機來模擬初始流場，即利用所選定的初始時刻以前時段的水位資料，來模擬計算出所選定的初始時刻的水位，這種方法只適用於被校正後的數學模型。否則模擬出來的初始流場可靠性也不大。一類邊界的初始水位及其源、匯項可根據實際觀測資料直接給定，二類邊界的初始水力梯度可根據邊界內外的水位觀測值通過等水位線分析或水力計算確定。計算區內初始參數亞區的劃分及其初始參數值一般根據含水層水文地質結構分析及其解析法所獲得的水文地質參數確定。

2.建立計算區刻畫地下水運動規律的數學模型

通過對上述概化後的水文地質概念模型的分析，就可建立計算區描述地下水運動的數學模型。實際上數學模型就是把水文地質概念模型的數學化，是用一組數學關系式來刻畫模擬計算區內實際地下水流在數量上和空間上的一種結構關系，它具有復制和再現實際地下水流運動狀態的能力。我們所談的數學模型主要是指由線性和非線性偏微分方程所表示的數學模型。對於一個實際的地下水系統來說，這樣的數學模型一般應包括描述計算區內地下水運動和均衡關系的微分方程和定解條件組成，定解條件中包含有邊界條件和初始條件。這樣的數學模型一般情況下很難通過常規的解析方法而獲得其精確解，通常都需藉助於現代化計算機，用數值方法對其進行求解以獲得其近似解。這就是數值模擬方法的來源。

地下水系統的數學模型根據研究的出發點和具體方法的不同，可分為以下幾種：線性模型和非線性模型、靜態模型與動態模型、集中參數模型與分布參數模型、確定型模型與隨機模型，等等。目前在礦井水文地質條件模擬預測中最常用的、最容易被一般水文地質技術人員所掌握的是確定型的分布參數模型。

3.數學模型數值求解的一般過程

（1）從空間和時間上離散計算域

當建立了刻畫地下水流特徵的數學模型之後，需要利用數值方法對模型進行求解，用於求解地下水流數學模型的方法較多，最常見的有有限單元法和有限差分法。無論是採用有限單元法還是有限差分法，都需要對模擬計算區域進行離散化剖分，剖分網格的形狀多種多樣，最常見的平面二維水流剖分網格有三角形和矩形，空間三維水流剖分網格有四面體和六面體，不管採用何種剖分方法，其解的收斂性與穩定性在很大程度上都取決於單元剖分的大小，為了保證解的收斂與穩定，剖分的單元一般不宜過大，特別是在水力坡度變化大的地方，單元應變小加密。對於非穩定流問題，還需要對模擬計算的時間段進行離散化，在水頭變化較快的時段內，時間步長應取的小些。在時段劃分上，一般原則是：在水頭變化快的時期，例如在疏排水的初期，時段步長應取得小些，劃分的時段應多些；在水頭變化緩慢的時期，例如在疏排水的中後期，時段步長可取得大些。一般情況下，有限差分法對時段步長的要求不像有限單元法那麼嚴格。

（2）校正（識別）計算區的數學模型

數學模型應是實際含水層及其水流特徵的復製品。根據水文地質模型所建立的數學模型，必須反映實際流場的特點，因此，在進行模擬預報之前，必須對數學模型進行校正，即校正其方程、參數以及邊界條件等是否能夠確切地反映計算區的實際水文地質條件。由此可見，校正模型實際上就是通過擬合實際觀測到的水文地質現象而反過來求得反映含水層水文地質條件的有關參數的過程。在數學上常稱之為反演問題或逆問題。

目前常用的識別數學模型所採用的方法大體可分為直接解法和間接解法兩大類。直接解法就是從含有水頭、水量和參數的偏微分方程或從已離散的線性方程組出發，把實際觀測的水頭代入，從中直接解出水量或參數的方法，即直接解逆問題。這類方法有數學規劃法、擬線性化法等。由於直接解法所需結點的水頭均應是實際觀測值，這在實際上很難辦到，所以該法應用較少。間接解法就是先給定一組參數或水量，代入已離散的方程，求解正問題，將計算值與實測值比較是否接近。在這個過程中，要不斷地去解正問題，不斷地比較計算值與實測值，最後求得最佳解。目前採用間接解法較為廣泛。間接解法又可分為兩種形式：一是人工調試計算參數，二是機器自動優選計算參數。人工調試就是人為給定未知量（參數或水量）進行正演計算，求得目標函數，並不斷地修改未知量，重復進行正演計算，直至求得的目標函數滿足誤差要求為止，這時的未知量即是所要求的參數或水量。人工調試方便、簡單，特別是在掌握計算區水文地質條件的基礎上，容易盡快達到誤差要求。機器調試是給定未知量的約束條件和參數自動尋優的數學方法，讓機器自動尋優，不斷地解正問題，求得目標函數達到極小值時的未知量，即是所要求的參數或水量。常用方法有單因素優選法、最優控製法等。

（3）數學模型的校驗

當通過參數反演獲得了數學模型的有關定量水文地質參數後，我們就獲得了用於礦井水文地質條件模擬預測的唯一確定的數學模型。為了在運行模型之前進一步確認模型的可靠性，可利用已知的水文地質觀測資料與模型運行的計算結果進行比較分析，以確認模型的正確性。如果校驗結果較好，則可利用模型進行礦井水文地質條件的預測分析，否則，尚需重新考核和校正數學模型。

（4）數學模型的運行與應用

經過識別和校驗後的數學模型，即可作為礦井水文地質條件和礦井涌水量預測預報的計算模型，可根據礦井開采條件、礦井水文地質要求進行多種問題的數值模擬計算。目前主要用於模擬預測不同條件下礦井疏降水量和疏降條件下的地下水流場。

4.數值方法的應用條件

雖然數值模擬方法在礦井水文地質條件定量分析和礦井涌水量預測方面有著明顯的優勢，但並不是在任何條件下都可得到很好的應用。數值模擬方法的成功應用必須建立在特定的條件之上。一般情況下，對一個礦區的礦井水文地質條件及其礦井涌水量進行數值模擬與預測時應具備下列基本條件：

1）必須有專門的地質與水文地質勘探資料嚴格控制礦井主要充水含水層（模擬的目標含水層）的空間賦存特徵，包括含水層的埋深、厚度、產狀、空間延展情況、結構類型（如含水層是單層的還是多層的）、頂底板岩層條件（有無天窗、缺失等），以及與主採煤層之間的位置關系。

2）要有專門的資料控制擬模擬的目標含水層的邊界條件。包括邊界的位置、物理結構、水文地質性質、可能出現的邊界隨時間變化（如分水嶺的移動、水位的動態變化、斷層受采礦擾動而發生活化等）、邊界外水體與邊界之間的關系等。

3）要有專門的水文地質試驗資料控制地下水的水動力學性質及其含水層的水文地質參數結構。包括地下水的流態（如層流還是紋流、一維流還是多維流、承壓水流還是無壓水流等）、含水層的滲透性能、越流條件、地下水水力梯度等。

4）要有大型群網觀測的抽放水試驗資料或具有區域性控製作用的地下水水力信息長期觀測資料。包括抽放水水量及其動態變化過程、抽放水過程中含水層水位及其變化過程、抽放水結束後地下水位回復程度及其回復過程。這些信息是進行水文地質條件反演和水文地質參數識別必不可少的信息。

5）其他影響含水層行為的相關信息。包括大氣降水及其時間分布、蒸發條件及其季節性變化、地表水系及其季節性變化、當地工農業用水及其開采情況、地表植被發育狀況等。這些因素會直接影響所建立的水文地質模型的准確性和真實性。

（三）超化礦水文地質計算的主要任務

1）分析處理L_1-3灰岩放水試驗的水量、觀測孔水位資料，建立礦井疏水量預測預報的水文地質概念模型。

2）通過水文地質參數的反演計算，形成礦井目標充水含水層的定量水文地質參數場和礦井疏水量預測預報的水文地質數學定解模型。

3）計算預測礦井不同開采水平L_1-3灰岩含水層的最大疏水量和最小疏水量，為建立礦井防排水系統提供依據。

4）計算預測礦井不同疏水條件下的地下水流場及其地下水位漏斗的擴展形態。

5）建立地下水疏降最優決策模型，提出最優疏水工程方案和疏水量時空分配方案。

（四）礦井水文地質計算方案

1）計算模型採用二維承壓水流數學模型。

2）計算方法採用有限元數值模擬技術。

3）計算所依據的基礎資料以井下放水試驗所獲得的所有可利用信息並結合歷史的勘探資料和礦井開采規劃資料。

4）計算程序為：水文地質概念模型的建立—水文地質條件模擬數學模型的建立—水文地質參數反演—礦井涌水量預測預報—疏水降壓孔的優化設計與計算。

C. 地下水是怎樣流動的呢拜託各位大神

1 地下水流系統動力學 地下水系統包括地下水、含水地質介質、水流邊界以及水源（如補給）和泄水口（如水泉、層際水流或水井）。水透過該系統流動，並保存在該系統中。在自然條件下，水從補給區到排泄區的流動時間范圍可從不足一天到100 多萬年。地下水流系統中儲存的水既有近期的降雨，也有沉積物沉積時被圈閉在沉積物中的水。 在發生水井抽水或補給率變化等水文變化之後，地下水系統中的水頭（水位）趨於平衡，可以用這種趨於平衡所需的時間來說明含水層反應時間的變化性。對於承壓地下水系統來說，反應時間可用下式計算： T * = S s L c2 / K （l ） 式中T *為盆地水力反應時間（T ）, S s為單位儲水量（L -1 ) , L c為盆地的某一特徵長度（L ）, K 為水力傳導率（L / T ）。水力傳導率是滲透率的一種量度，其范圍可達12 個數量級，而且地下水系統邊界之間的距離可從幾米到數百km 。利用（1 ）式計算出的兩種理想地下水系統的水力反應時間是：承壓型河流含水層系中的水平水流為0.1 天（144 分鍾），而厚的區域性低滲透率單元中的垂直水流為4.0 x l07 天（約11萬年）。 流經地下水系統的時間取決於水頭的時空梯度、水力傳導率和系統的孔隙度。流經系統的時間不同於趨於平衡的水力反應時間。例如，上面已計算出，承壓型河流含水層系中的水頭可在不到一天內對水力變化做出反應，而在自然條件下水穿越系統整個厚度所需的時間則大約為3 萬天（82 年）。基岩中破碎岩石系統的有效孔隙度小於砂、礫石等未固結多孔隙介質系統的有效孔隙度，而且穿過破碎岩石系統的流速可能較快。例如，對於在威斯康星州破碎白雲岩中所打的城市水井來說，水流經幾公里距離所需的時間估計不到一年，在這種情況下，補給量和抽水量的季節性變化會對流動時間的變化產生影響。在較緩慢的地下水系統中，如在泰國的曼谷盆地，在了解數萬年時間內地下水的流動時，要考慮到氣候和地質條件的長期變化。地下水的長期流動實際上也影響著所有地質作用過程，其中包括成岩作用、成礦作用以及石油的聚集。 在確定污染物在地下水系統內的遷移時，水的運移時間是重要的。全世界由於地表污染源所造成的大范圍地下水污染，反映出淺部地下水的年齡通常為幾十年或更短。由出水量大的水井引起的水力梯度會進一步縮短污染物向水井匯流的時間。 從地下水系統抽取的水最初來自於儲水區。隨著時間的流逝，水頭會在距抽水點更遠處逐漸下降，抽水造成的影響會通過地下水系統傳播。最終，抽水造成的影響會到達某一邊界（如河流），在該邊界上，要麼加大對地下水系統的補給量，要麼減少該系統的排水量。通常假設，如果對地下水的抽取速率不超過自然補給率，則這種抽水速率是「安全的」或「可持續的」，然而這種假設是錯誤的，因為它忽視了地下水系統排水量和補給量的這些變化。 多年來，水流和溶質遷移的計算機模型一直是評價地下水資源的綜合手段，它們被用來闡明從局部污染到由於大陸范圍的流體遷移而造成的大型礦體的一系列間題。這些模型的預測能力可用於對假說進行檢驗，這樣既加深了我們對現狀的了解，又加強了有關含水層對未來氣候或人為壓力的反應的預測。最近有人把地下水流動模型與地表-大氣模型和遷移模型與地球化學反應模型聯系在一起，從而擴大了可以闡述的問題種類的范圍。自動校準方案和對不確定性的分析增強了模型的實用性，而新的計算機可視化手段則提高了我們對含水層性質的變化性對地下水流動型式的影響的了解。 模型的預測精度是受模型的正確性（即對相關作用過程的正確顯示）和模型參數的不確定性制約的。後一種不確定性是由於參數值的可測定精度有限而造成的，但更重要的是由於含水層的特徵中所固有顯著非均質性造成的。地下水水文學中的一個根本問題，是無法充分描述和顯示這種非均質性，即便採用經過改進的模型，這個問題仍將繼續制約著模型預測的可靠性。空間非均質性與模型不確定性之間的聯系還取決於所要回答的問題的類型。例如，對水頭在含水層中的分布進行合理估計，只要求對空間非均質性進行有限的了解。而另一方面，某一特定位置的化學濃度預測的置信度，對水力性質空間分布的較小不確定性可能都是十分敏感的，甚至對於相對均勻的多孔隙介質也是如此。

D. 計算方法、步驟

（一）建立水文地質概念模型

解析法對水文地質條件限制較多，有嚴格的理想化要求，而實際水文地質條件往往十分復雜，為了能夠用解析法計算，必須對水文地質條件進行合理的簡化和概化，經過簡化和概化後的水文地質條件稱水文地質概念模型，它是對地下水系統的定性描述。

1.分析疏幹流場的水力特徵

礦床的疏幹流場，是在天然流場背景下，疊加人為開采因素演變而成的，因此分析疏幹流場各種水力特徵時，均應以天然條件為基礎，充分考慮開採的影響。

（1）區分非穩定流與穩定流

一般，疏干排水時，礦區地下水多為非穩定狀態，但當疏干排水量小於地下水補給量時，可出現穩定狀態。

礦山開采初期（開拓階段），開拓井巷不斷發展變化，疏干漏斗的外邊界不斷擴展，礦坑涌水量以消耗含水層儲存量為主，該階段疏干場一般為非穩定流，礦山開采後期（回採階段），疏幹流量主要受流場外邊界的補給條件所控制，在補給條件不充分的礦區，疏幹流場以消耗含水層儲存量為主，仍為非穩定流，在補給條件充足的礦區，或具定水頭補給邊界的礦區，礦坑涌水量（或疏乾量）被補給量平衡，一般出現相對的穩定流，礦坑涌水量預測可以穩定井流理論為基礎。

（2）區分層流與紊流

礦區地下水在疏干條件下與天然運動狀態相比，在大面積內仍為層流，僅在疏干工程附近常出現紊流，故達西定律（直線滲透定律）仍然是建立確定性模型的基礎。

一般，常以抽（放）水試驗為依據，用單位涌水量（q_i）法對層流、紊流進行判別，計算式為：

承壓水

圖13-7 水位降深為S_k的Q-t曲線

E. 水文地質概念模型

一、基本模型

本區歷史上曾做過兩次地下水流數值模型（中國地質科學院水文地質環境地質研究所、甘肅省第二水文地質工程地質大隊，2000年6月；清華大學水利系水文水資源研究所與甘肅省水利水電勘測規劃設計研究院，2004年10月），兩次都是地下水二維流模型，且對研究區建立了三個模型分別模擬三個盆地。那麼，本研究如何設計模型?

如前所述，研究區諸盆地內，洪積扇頂部出山口河水滲漏補給後，地下水將以垂向為主地向下流動；洪積扇前緣地下水向泉口溢出也存在垂向分量；含水層之間弱透水層中的地下水也以垂向流動為主；本區所有的河流都是高度非完整的河流，非完整河底部的地下水也存在明顯的垂直流動；非完整抽水井附近的地下水流存在垂直流速分量，等等。因此，本區地下水流屬於三維流類型。

那麼，採用准三維流模型還是三維流模型呢?考慮到具體條件和建模目標——區域地下水資源評價與管理，實際含水系統的層數很多而模型分層又不宜過多，因此概化後的「弱透水層」中含有一定數量的砂層，即「弱透水層」是各向異性的。陳崇希等（2001）在蘇州市地面沉降模型研究中指出：「准三維流模型『誤差小於5%』的結論僅適用於弱透水層為各向同性介質；對於弱透水層的各向異性比為10時，誤差已超過27.7%，而且隨著弱透水層單位儲水系數的減小和模擬時間的延長，誤差還要增大，已存在超過30.5%的情況」。薛禹群（2003）也提出相似的定性見解。而本區概化的「弱透水層」滲透系數的各向異性比要遠遠地大於10，其誤差要明顯地增大。基於上述研究的成果，盡管三維流模型比准三維流模型的模擬工作量要大得多，本項目仍要採用三維流模型。這是本項研究與上述兩個模型設計上第一個重要差別。

另外，特別需要強調的是，本區的井孔中含一定數量的混合井孔，如果忽視它們的存在，將會導致模擬失真。因此，本項研究所建立的基本模型是含混合井孔的地下水三維不穩定流模型。

二、模型範圍及邊界條件

上面提到的兩個先期完成的模型都將北山山前戈壁前緣作為零流量邊界。本項研究考慮到北山仍得到降水入滲補給，且北山山前戈壁含水層組本身具有相當的地下水儲存功能，即有很高的調節儲量的能力，為此將模型範圍向北擴至北山地下水分水嶺處，即北山地下水分水嶺以南的北山山區均劃入模型內。此外，考慮到分隔玉門-踏實盆地和安西-敦煌盆地的北截山較為狹長，若作為隔水邊界處理則與實際條件不符；若作為弱透水邊界處理則難以給定邊界流量。為此，將北截山作為非均質的一個分區，將相鄰兩個盆地連成一個整體。玉門-踏實盆地與花海盆地之間原本存在數百米寬北石河沖積層的聯系帶，更不宜人為分隔成兩部分。如此，本項研究將三個盆地作為一個整體模型研究（圖5-1），其優點是顯而易見的，它既客觀地反映了諸盆地之間存在的水力聯系，又避免了給定人為邊界流量的困難。這是本項研究與上述兩個模型設計上第二個重要差別。

本項研究所建模型的總面積約為46140km²。根據前面對水文地質條件的論述，模型邊界條件的設置如下（見圖5-1）：

（一）北邊界

取北山地下水分水嶺為模型的北邊界，該邊界為零通量邊界。由於缺少地下水位資料，假定地下水分水嶺與地表分水嶺一致，其誤差對模型地區不致產生大的影響。

圖5-1 模型範圍及邊界條件圖

（二）南邊界

取南部寬灘山、南截山、火焰山和卡拉塔什塔格山山前基岩與第四系鬆散沉積物的分界線為模型的南邊界，該邊界為弱透水邊界。它反映南部祁連山區地下水對本研究區的補給。

（三）東邊界

取花海盆地和金塔盆地的地下水分水嶺處為模型的東邊界，該邊界為零通量邊界。

（四）西邊界

西邊界的設置是比較困難的。歷史上疏勒河是出甘肅省西界流入新疆的羅布泊，盡管疏勒河河水早已斷流，但估計地下水流仍向下游徑流。考慮到下游缺乏資料，包括含水系統的結構岩性和地下水位等基礎數據，我們取安西-敦煌盆地最西的一個觀測孔處（庫穆塔格沙漠的東界）為模型的西邊界。安西-敦煌盆地有部分地下水向西流入沙漠。該邊界在模型識別中作為第一類邊界，在模型預測中作為第二類邊界。

（五）上邊界

取潛水面為模型的上邊界。地下水通過該邊界接受河流、渠系、大氣降水（凝結水）及田間灌溉回歸水的入滲補給，並以蒸發蒸騰、泉水溢出的方式排泄。南陽鎮-雙塔水庫區段的疏勒河地下水與地表水直接聯系，取為第一類邊界。

（六）下邊界

因為下更新統（Qp₁）為膠結砂礫岩，故取中更新統（Qp₂）的底界為模型的底邊界，該邊界為零通量邊界。

F. 三維水文地質建模的應用需求

地質體三維模擬與可視化軟體應能輔助地質人員進行地質數據分析、地質體模型構建、地質體可視化和各種地質應用分析。主要功能需求包括如下方面：

1.地質體三維模擬資料庫管理

能夠管理用於構建地質體三維模型的基礎地質數據（指由地質人員整理、分析形成的各種地質數據，如鑽孔數據、樣品數據、物探成果數據、地質圖紙數據、影像數據等）和地質模型數據（地質界面和地質體模型構建後所生成的結果數據，如DEM數據、地質體的幾何特徵數據和屬性特徵數據）。

要具備對基礎地質數據的錄入、查詢、坐標變換、統計分析、2D和3D顯示、導入和導出功能。要具備對三維地質模型數據的查詢、修改等功能。

2.基於表面的三維地質體建模

能夠為用戶提供一個進行基於表面的三維地質體建模的可視工作環境，根據經過預處理的地質模型基礎數據，具有自動生成各種地質界面的功能，並可進行人機交互編輯，依據這些經過多次修改的面能夠自動生成三維地質模型數據。

3.基於剖面的三維地質體建模

剖面圖是常見的地質數據形式。所以應能夠為用戶提供一個進行基於剖面的三維地質體建模的可視工作環境，根據經過預處理的地質模型基礎數據（鑽孔數據、剖面圖數據等），能互動式生成三維地質模型數據。

4.三維地質體屬性變化分析

礦體等地質體內部的屬性（如品位）在空間上是變化的，進行屬性變數的空間變化分析是地質工作者的一項重要任務。為了達到這一目的，可以對已建地質體幾何模型進行立方體、長方體等形式的剖分，進而構建屬性變數的空間變化模型，實現屬性變數的空間變化分析。

5.地質體的三維可視化和屬性查詢

地質模型是由地質體組成的。對所建立的一個地區的三維地質模型要能夠進行整體顯示和部分顯示。因為當整體顯示時，礦體的空間變化可能就難以觀察清楚，所以需要將一個礦體單獨顯示或將某幾個礦體組合顯示。

對於所顯示的地質體，要能夠通過滑鼠點擊後顯示出其有關相應的屬性信息，如斷層的產狀、長度等信息。

6.地形面和地質界面的可視化飛行

為清晰地展現不同面（如地表面）的空間分布狀態，採用飛行可視化技術，依據軟體實現飛行觀察。飛行的空間面對象主要包括地表面和地質界面，其中地表面飛行的對象是DEM和衛星影像合成的三維表面，地質界面飛行針對的是一個地下的空間界面，要求建立虛擬的地下環境。

7.選定范圍地質模型的生成功能

用戶可以在有效區域范圍內任意選定特定區域，系統能自動按垂向提取范圍內的地質體進行動態顯示。這個過程應根據人的視覺感受，通過針對數據對象的提取與運算，生成基於地質模型不同對象的三維顯示效果；通過對圖形觀察視角的不斷改變，顯示三維立體對象的不同側面，形成動態變化效果。

8.任意地質剖面生成

提供用戶觀察與X-Y平面垂直的任意剖面數據分布情況的功能，要求剖面將指定深度范圍內的剖面圖繪制出來。

9.等值線和等值面的生成

基於三維地質模型數據，繪制選定范圍內任一空間面的等值線（如煤層底板等高線）和等值面（如礦體內部等品位面），繪制結果可採用公共數據格式存儲。

10.地質體剖面面積、體積的計算和礦體儲量計算

能夠在所建立的三維地質模型的基礎上，方便地計算出地質體的剖面面積、地質體的體積和礦體的儲量。

11.基於地質模型的地下遨遊

為更好地顯示地下三維地質對象的空間關系，用戶可以以鑽洞的方式對地下進行遨遊。

12.地質制圖

在構建的三維地質模型的基礎上，可以派生出許多有用的地質分析數據，如等值線、等值面、任意方位和任意傾角的地質剖面，地質人員需要將這些分析數據以圖紙的形式表達出來，所以系統還需要具有制圖的功能。

G. 水文地質概念模型概化導則

地下水流模擬是一實用性很強的技術，解決現實問題是它的根本目的。因此，建專立的水文地屬質概念模型需與一定時期的科學技術水平以及研究區的水文地質調查研究程度相適應，能用於解決社會、經濟發展中所面臨的地下水模擬與管理問題。

水文地質概念模型概化導則按 GWI-C6 執行。

H. 　水文地質模型概化及其數值模擬模型

一、水文地質模型概化

（一）邊界條件

根據工作區區域地質、水文地質條件，以鬆散岩類孔隙水含水層為本次計算的目的層，將黑龍江、烏蘇里江、松花江、鴨蛋河、梧桐河、撓力河、七星河、別拉洪河及安邦河概化為一類水流邊界，即地下水與江河水位有較密切的水力聯系。平原與山地接觸地帶，地下水自山地向平原側向徑流的補給邊界，概化為二類水流邊界；而季節性積水的沼澤濕地為三類邊界。山地區（含區內的殘山丘陵區）及前第四系為隔水含水層，在計算模型中為不計算單元。邊界條件概化結果見圖4-1 。

圖4-1 三江平原邊界條件概化圖

（二）含水層參數分區

第四紀不同時期，工作區不同位置沉積了不同厚度與粒度的鬆散堆積物，使含水層在水平方向與垂直方向上均有分區性，根據工作區不同時期的岩相古地理特徵，將工作區的含水層參數分成6個區，見圖4-2～圖4-5。

二、數值模擬模型

將工作區地下水流系統概化為非均質各向同性、三維、非穩定流系統，用以下微分方程的定解問題描述：

三江平原地下水資源潛力與生態環境地質調查評價

式中：Ω為滲流區域，即工作區內有效計算單元所構成的區域，面積為39 415.68km²；Γ₀為滲流區上邊界，即地下水的自由表面；Γ₁為一類邊界，包括黑龍江、烏蘇里江、松花江、鴨蛋河、梧桐河、撓力河、外七星河、別拉洪河及安幫河等常年流水且與地下水有密切水力聯系的河流；Γ₂為二類邊界，即除黑龍江、烏蘇里江以外的工作區邊界；Γ₃為三類邊界，即工作區內季節性積水的沼澤濕地；h為含水層水位標高，m；h₀為初始水位，m；h₁為一類邊界江河水位，m；h_s為三類邊界沼澤濕地水位，m；K為含水層滲透系數，m/d；S為地下水自由面以下含水層的儲水系數，m^-1；μ為潛水含水層在潛水面上的重力給水度；ε為含水層源匯項，d^-1；p為潛水面的蒸發和降水補給等源匯項；q₂為二類邊界單位面積流量，m³/d；σ為沼澤濕地底部弱透水層的阻力系數，σ＝L/K_S，L為弱透水層厚度，m，K_S為弱透水層垂向滲透系數，m/d。

上述數學模擬模型應用三維模擬計算軟體進行求解計算。

I. 礦山地質環境多級模糊識別模型評價方法研究

——以小秦嶺礦區為例

邢永強¹鄭釗²宋鋒¹薛良偉¹王衛¹

（1.河南省國土資源科學研究院，鄭州 450016；2.天津大學建築工程學院，天津 300072）

《信陽師范學院學報》（自然科學版），文章編號：1003-0972-(2008)-02-0240-03

摘要 應用多級模糊識別模型進行礦山地質環境評價，克服最大隸屬度原則所不適用的地方，而且以相對隸屬度、隸屬函數為基礎，使隸屬度、隸屬函數的計算更加容易。建立了礦山地質環境評價模型，並應用於小秦嶺礦區地質環境質量綜合評價中，評價結果顯示，該方法對礦山地質環境評價可行、可靠。

關鍵詞多級模糊識別模型 地質環境評價 小秦嶺

礦山地質環境評價的數學方法有多種，如地理信息系統（GIS）、BP人工神經網路、灰色系統理論、模糊數學方法等都是較為常用的評價方法，這些方法各有其優缺點（田婧等，2007）。礦山地質環境質量評價是解決評價對象的多指標模式識別和排序問題，模糊數學理論中的最大隸屬度原則在此問題中得到了廣泛的應用。然而，最大隸屬度原則有時忽略某些信息，導致結果失真。陳守煜（1998）提出的多級模糊識別模型以相對級別特徵值作為判斷或識別的依據，克服了最大隸屬度原則所不適用的地方，使隸屬度、隸屬函數的計算更容易。因此，筆者採用多級模糊識別模型來評價某典型礦山的地質環境狀況。

1 模糊識別模型介紹

1.1 指標相對隸屬度公式

設某樣本需要對模糊子集A進行識別，樣本m個指標組成的特徵向量為X=(x₁x₂…x_ix_m)^T，其中：x_i為樣本指標i 的特徵值，i=1，2，…，m。

如樣本依據m個指標c個級別的已知指標標准特徵值進行識別，則有指標標准特徵值矩陣：

環境·生態·水文·岩土：理論探討與應用實踐

式中：y_ih為指標i級別h的標准特徵值，h=1，2，…，c，i=1，2，…，m。

通常將從1級至c級指標標准特徵值減小稱為指標遞減型；從1級至c級指標標准特徵值增加稱為指標遞增型，則遞減型指標對A的相對隸屬度（隸屬函數）為

環境·生態·水文·岩土：理論探討與應用實踐

h級指標標准特徵值對A的相對隸屬度為

環境·生態·水文·岩土：理論探討與應用實踐

遞增型指標對A的相對隸屬度（隸屬函數）為

環境·生態·水文·岩土：理論探討與應用實踐

h級指標標准特徵值對A的相對隸屬度為

環境·生態·水文·岩土：理論探討與應用實踐

1.2 模糊識別模型的建立

通過隸屬度分別把指標與指標標准特徵值矩陣變換為A的指標與指標標准特徵值的相對隸屬度矩陣：

R=(r₁r₂… r_m)^T=(r_i)

環境·生態·水文·岩土：理論探討與應用實踐

設樣本A的各個級別相對隸屬度矩陣為

U=(u₁u₂…u_c)^T=(u_h)

式中：u_h為樣本對A的級別h相對隸屬度，h=1，2，…，c。

由樣本在級別區間[a，b]范圍內，故矩陣U應滿足規一化約束條件：

u_h=1，同時應有u_h=0，且h＜a或h＞b。

為了求解樣本對模糊概念A的級別h的最優相對隸屬度，建立目標函數：

環境·生態·水文·岩土：理論探討與應用實踐

求解

環境·生態·水文·岩土：理論探討與應用實踐

求得模糊識別模型為

環境·生態·水文·岩土：理論探討與應用實踐

1.3 相對級別特徵值的求解

設級別變數h以對應的相對隸屬度為權重，其總和H_(u)=

u_h·h，稱為相對級別特徵值或相對狀態特徵值。

相對級別特徵值H_(u)是一個描述狀態或級別的無量綱數，且1≤H_(u)≤c，相對級別特徵值由於利用了全部相對隸屬度信息，其判斷結論比最大隸屬度原則更符合實際情況。據此，可應用相對級別特徵值H_(u)反饋得到相應的級別，對樣本作出歸屬何種級別的判定。

2 應用實例

2.1 研究區概況

小秦嶺礦區位於河南省西部靈寶市境內的豫、陝、晉三省交界處，礦業經濟發展很快，礦區內有國家和地方黃金企業數十家，礦山坑口數千個，礦業已成為靈寶市的支柱產業。然而，由於當地礦業資源的管理、開發體制不健全，小秦嶺礦區礦業資源又作為淘金者的首選目標，該地區在淘金者的亂采濫挖下，不可再生的礦業資源和礦區地質環境遭到了嚴重的破壞。此外，人為的破壞還給該地區帶來了崩塌、滑坡、泥石流等地質災害。因此，運用適當的方法對小秦嶺礦區地質環境進行綜合評價，對以後礦山恢復治理工作有重要的指導意義。

2.2 評價單元劃分

結合對小秦嶺礦區的實地調研結果，並遵照客觀、公正、科學地反映礦區地質環境區域差異的原則，將評價區劃分了87個評價單元。採用先定性分析礦區的主要地質環境問題，並綜合考慮地形地貌特徵、水系發育特徵、人類活動強度等因素，對礦區內問題比較突出地區劃定評價單元網格；對於其餘地區，則按照3km×3km的正方形網格來劃分評價單元，在單元劃分的同時還注意與行政界域、水系界域的相互包容以及對邊緣單元、小單元的適當合並。此外，在遵循地質環境客觀特徵的基礎上，還考慮到兼顧局部特殊要求的情況，如在豫陝兩省交界的地區則按照行政分區邊界來劃定。劃分結果如下：單元dx1和dx2為大西峪區間；w1～w4為文峪區間；單元z1～z3為棗香峪區間；單元dh1～dh3為大湖峪區間；單元zy為藏馬峪和閻家峪區間；單元f為夫夫峪區間；單元g為觀音峪區間；單元j為荊山峪區間；單元i1～i10為蒼珠峪、白花峪、槍馬峪、楊砦峪、朱家峪相應的區間；單元1～61是按照正方形網格與各類界域邊界交匯並進行適當合並或裁減的評價單元。

2.3 礦區地質環境質量評價指標分級

在綜合比較成天翔等（2007）、徐友寧等（2003）、蔡斌等（2006）關於選取評價指標研究成果的基礎上，從小秦嶺礦區地質環境現狀條件出發，綜合考慮研究區自然條件、人類活動影響、資料收集情況等因素，選擇了地表坡度、岩土體抗侵蝕性、植被覆蓋率、年降水量、地質災害、水土流失、地表水污染、人類工程活動強度、礦渣堆積量9個評價指標。

對於選取的評價指標，按照地質環境質量「優」、「良」、「中」、「差」、「極差」劃分為5個級別，各級別相應指標的標准值如表1所示。

表1 小秦嶺礦區地質環境評價指標分級標准Table1 Index classification standard of geology environmental in Xiaoqinling field

2.4 確定權重系數

目前系數確定的方法很多，大致可分為德爾菲、層次分析等主觀賦權法和主成分分析、因子分析等客觀賦權法，運用主觀賦權法摻雜了決策者的主觀隨意性，而運用客觀賦權法卻缺乏決策者的意願，故本次研究採用主觀賦權與客觀賦權相結合的方法，先由主成分—因子分析賦權法計算出一組初始權重，再帶入評價模型進行計算，如果計算結果合理則直接採用該指標權重，如果計算結果差別較大，則在初始權重的基礎上再進行適當微調，最終求出一組合理的權重系數，如表2所示。

表2 評價指標的權重系數Table2 Weighing coefficient of evaluation index

2.5 小秦嶺礦區地質環境評價結果及分析

用VB語言編制相應的計算程序、開發評價模型，結合各單元評價指標的量值和小秦嶺礦區實地調研情況，給定地質環境質量評價的最終結果，見表3和圖1。

表3 小秦嶺礦區地質環境質量評價結果Table3 Evaluative result of geologic environmental quality in Xiaoqinling field

圖1 小秦嶺礦區地質環境質量評價效果圖

將該評價結果與已知物元分析法評價結果相比較，僅有39號、55號、58號、dx2號、i5號、dh3號單元評價結果有偏差，且偏差在允許范圍內。

3 結論

多級模糊識別模型分析較物元分析、模糊理論分析等方法計算簡便，又克服了最大隸屬度原則所不適用的地方。該方法用於礦山地質環境評價，對屬於相同級別的檢測單元間的差別亦可分區，可以直接看出各單元環境質量的優劣程度，從而為環境治理工作提供更直接的信息。實踐證明，用該方法對礦山地質環境進行評價可行、可靠。

參考文獻

蔡斌，胡卸文.2006.模糊綜合評判在綿陽市環境地質風險性分區評價中的應用.水文地質工程地質，（2）：67～74.

陳守煜.1998.工程模糊集理論與應用.北京：國防工業出版社.

成天翔，張駿，杜東菊.2007.天水地區斷裂活動性與地質災害的相關性研究.工程地質學報，15（1）：33～37.

田婧，韓秀麗.2007.工程地質環境評價方法.河北理工大學大學學報，29（1）：125～128.

徐友寧，袁漢春，何芳，陳社斌，張江華.2003.礦山環境地質問題綜合評價指標體系.地質通報，22（10）：829～832.

Research of Multi-Classification Fuzzy Pattern Recognition Model in Mine Geological Quality Evaluation

—The Case of the Mine in Xiaoqinling Field

Xing Yong-qiang¹，Zheng Zhao²，Song Feng¹，Xue Liang-wei¹，Wang Wei¹

（1.Henan Land and Resources Research Institute，Zhengzhou 450016；2.College of Civil Engineering，Tianjin University，Tianjin 300072）

Abstract：Based on the relative membership degree and function，the model of multi-classification fuzzy pattern recognition is applied in assessment of mine geo-environment quality，which counteracts the defects of maximum membership degree law and makes it easier to calculate the relative membership degree and function.The paper has established assessment model of geo-environment quality and uses it to assess the mine geo-environment quality of Xiaoqinling field.The result indicates that the method is feasible and reliable.

Key words：multi-classification fuzzy pattern recognition model；geological environment evaluation；Xiaoqinling field

J. 定義水文地質數學模型的定解條件有哪些

水文地質，地質學分支學科，指自然界中地下水的各種變化和運動的現象。回水文地質學是研究答地下水的科學。它主要是研究地下水的分布和形成規律，地下水的物理性質和化學成分，地下水資源及其合理利用，地下水對工程建設和礦山開採的不利影響及其防治等。