摘要:本文在論述山區土壤水資源量計算模型的基礎上，以河北省太行山區為研究對象，初步計算太行山區耕地1980-2000年21年間最小產生的土壤水資源總量為312.6億m3,土壤水資源量平均深度為349mm。土壤水資源量、徑流量分別占評價區降水量的61.45%和33.05%。通過計算分析得出土壤厚度是決定土壤蓄水能力的主要因素。由于不同土質土壤儲水能力不同，因此在同一土層厚度范圍內各土質土壤水資源變化規律不同，表現為砂土>砂壤土>壤土>粘壤土>粘土。同時。土壤水資源量的年內分配規律表現為土壤水資源量與降水量年內分配同頻率。握了上述規律，有效調整農作物種植結構，降低了農業生產成本，開展農業節水措施是十分必要的。
　　關健詞：土壤水資源量,計算模型,太行山,分配規律,種植結構
　　
　　1前言
　　土壤水是地球水體的重要組成部分之一。陸地上植被和作物所利用的水分主要是從土壤中直接吸取的。它影響著農作物生長、生態環境建設以及水資源的合理分配與高效利用。因此提出一套適合不同分區、不同土質的土壤水資源的計算方法，了解河北山區土壤水資源分布時空規律與土壤蓄水量的評估對實現農業乃至整個山區經濟的可持續發展都具有重要意義。
　　2土壤水資源評價模型             

                                   
　　土壤水資源計算模型的目的是摸清現狀我省山區土壤分區和土地利用條件下土壤水資源的分布規律，為今后我省改進耕作方式、調整農村產業結構，高效利用土壤水提供基本依據。并以此為目的，研究天然條件下土壤水資源的時空變化規律，分析在天然條件下可以得到補充，并且對作物有效的那一部分水量是本模型建立的基礎。因而在評價土壤水資源時不考慮灌溉水。
　　物理模型基礎

                                                    
　　范如華先生認為：降雨補給土壤水，其過程如圖2-2示，是由土層上部達飽和含水量后逐漸下移，最后使土壤包氣帶達田間最大持水量后暫時平衡，這就是依靠土壤毛管力支持的懸著毛管水，可供作物根系吸收和棵間蒸發的雨水。
　　2.1模型建立的幾個基本假定
　　本次模型建立是研究大范圍內的土壤水資源轉化規律，考慮到降雨與下墊面條件的空間變異復雜，因此對模型進行假設概化。
　　2.1.1研究的狀態是一種假定的理想狀態。土壤水分的上行消耗為裸地的蒸發，而不考慮植物蒸發。得到的結果是一個潛在的最小土壤水資源量，從而根據河北省太行山區最少能產生土壤水資源的數量，達到合理利用有限的土壤水資源，指導土地利用結構調整和生產的目的。
　　2.1.2評價層內土壤為均質土壤
　　2.1.3采用的五種土壤質地類型（砂土、砂壤土、壤土、粘壤土、粘土），就某一種土壤質地而言，土壤調蓄能力（儲水、入滲等特性）完全相同；
　　2.1.4年降雨及年內分布決定了土壤水資源的狀況，因此采用降雨的頻率分析土壤水資源的頻率（程伍群）；
　　2.2.土壤水資源量模型確定

                                                          
　　本模型采用補給法計算土壤水資源，從土壤水資源可利用角度出發，利用含水量變化過程來推算土壤水資源量，不僅包含了本年降水的有效補給數量，由于年初土壤含水量是上一年土壤蓄水量的反映，因而使計算的土壤含水量，更接近本年的實際土壤水資源可利用量，較好地說明水分循環和水量平衡的全過程。
　　對于一個研究單元，確定的土壤水資源評價模型的數學表達式為：（1）
　　其中：
　　——一次降雨產生的土壤水資源量；
　　——日降水量
　　——凝結水量；
　　——為入滲到評價層以下的水量；
　　——無效降水量；
　　——雨期蒸發量；
　　——為潛水蒸發量；
　　——可以通過不同坡度、土層厚度的降雨徑流關系確定。
　　其中：
　　1、（凝結水量），經試驗證明，該值的總量較大，但平均每天的數量較小，在轉化為土壤水之前主要用于蒸發，對土壤水資源量的影響不大，因此忽略不計；
　　2、(無效降雨量)，當本日有降水時，根據降水歷時t和當日蒸散發能力計算無效降水量。即＝/24×t采用氣象學方法計算，或采用E601蒸發皿折算。）
　　3、（雨期蒸發量），由于雨期的濕度大，溫度偏低，蒸發非常小，因此忽略不計；
　　4、（潛水蒸發量），河北省太行山區地下水位較深（＞30m），并有巖石不透水層隔離，因此可以忽略不計，對于山區丘可以認為該值為零；
　　因此：
　　式（1）簡化為：
　　（2）
　　在山丘區一般情況下，與最終流出評價層，不參與土壤水資源評價，因此，把、統稱為徑流（）。
　　因此得出：
　　（3）
　　假定評價層深度內田間持水量為Wmax（mm），田間調萎含水量為W0（mm）其中
　�。�4）
　�。�5）
　　徑流量是區域降水特征和下墊面特征相互作用的產物。由于包氣帶的土壤組成、水分動態的變化等因素的差異，造成產流機制不同。主要包括兩種基本的產流機制，即超滲產流與蓄滿產流。由于河北省太行山區土壤砂性顆粒含量高,降水特點為雨強大,歷時短，主要以超滲產流為主；又因為該地土層薄,降水到一定程度又容易蓄滿。因此，為了提高土壤水資源量計算精度，本模型徑流量計算同時考慮超滲、蓄滿兩種產流機制。
　　對于每一種類型土質土壤，設置一個包氣帶日入滲能力FM，當（P－Pw（無效降水））大于FM時，超過部分產生徑流R超。即：
　　當P－Pw>FM時
　　R超=P－Pw－FM
　　當P－Pw≤FM時
　　R超=0
　　如果單次降雨量、無效降水量、計算超滲產流量之差大于土壤田間持水量與起始含水量（W’）的差，則，如果單次有效降雨量小于土壤田間持水量與起始含水量的差，則。
　　該值計算過程中需要確定降雨開始時的土壤實際平均含水量W/，需要通過的推求關系確定。計算公式為：
　�。�6）
　�。�7）
　　（8）
　　其中：
　　——第n+1日的土壤含水量；
　　——第n日的土壤含水量；
　　——第n日的有效降水量；
　　——第n日的蒸發量；
　　——為土壤濕度系數；
　　——為水面蒸發量；
　　——為土壤蒸發量與水面蒸發量的折減系數；
　　將各次降雨轉化而成的土壤水資源量累加，就可以得到一年內的土壤水資源量。對于一年，可以表示為：
　　（9）
　　該值的計算可以的按日進行全年累計，獲得全年的土壤水資源量，也可以對各作物的生長階段進行累加，獲得不同生長階段的土壤水資源量，這樣得到的結果包括兩個部分，即：年土壤水資源總量和土壤水資源的年內分配。
　　2.3．太行山土壤水資源計算分區及計算單元
　　土壤水資源分區是土壤水資源評價的基礎。根據土壤類型分布及其性狀等，同時采用有關專題分區結果確定河北省太行山土壤水資源分區。即太行山土壤水資源分區共劃分為一個一級區，兩個二級區，三個三級區。（見表2-1）。在三級分區的基礎上，以縣為單位，根據五種土壤質地（砂土、砂壤土、壤土、粘壤土、粘土）劃分計算單元。
　　太行山區土壤分區表
　　                                                                                                 表2-1

                                      
　　　　2.4太行山區土壤水資源評價層確定
　　太行山區土層薄且礫石含量多，因此評價層確定為整個土層。由于山區土層厚度差異非常大，因此我們把太行山區不同土質土層厚度劃分為四個范圍,分別為>100cm、55-100cm、30-55cm、<=30cm，并計算位于不同土層范圍內各計算單元土層厚度加權平均值，作為該范圍內的評價層厚度。
　　2.5典型區土壤水資源模型參數率定
　　選用各試驗站八十年代以來的實測及調查的降水、徑流、土壤含水量、灌溉水量資料，分別建立各典型區土壤水資源數量計算模型，進行模型參數率定和檢驗。
　　2.5.1入滲能力率定
　　包氣帶日入滲能力FM主要結合各典型區逐年的地表徑流量資料進行率定。率定得到的包氣帶日入滲能力FM隨土質不同變化較大，砂土、砂壤土、壤土、粘壤土、粘土的FM分別為50mm/d、45mm/d、40mm/d、35mm/d、30mm/d。
　　2.5.2蒸發損失系數K率定
　　蒸發損失系數K與日蒸發能力和土壤水含量有關，由于兩影響因素均具有年變化規律，故率定k值時，按月采用同值進行率定。
　　率定得到的各月蒸發損失系數k為多年平均值，綜合反映了氣象因素和土壤含水量的影響。比較各典型區率定得到的蒸發損失系數，代表不同土質各分區的率定k值，基本差別不大，因此，推廣到土壤水資源計算單元時，同土層不同質地可采用同一k值，作為模型初始試算值。見表2-2
　　
　　                                 土壤水資源計算分區不同土層蒸發衰減系k值表
　　                                    　　
　　2.6水分常數的確定
　　計算土壤水資源量必須首先確定田間持水量（θc(z)）、凋萎含水量θwp(z)兩個基本土壤水分常數。
　　經過收集整理河北省山區不同土壤質地代表點兩個基本土壤水分常數的有關資料，然后根據土壤計算單元的質地及鄰近代表點相應質地，確定該土壤水資源計算單元兩個基本土壤水分常數的采用值。見表2-3
　　太行區各分區土壤水分常數表
　　                                                                                  表2-3

                                       
　　　　3.河北省太行山區土壤水資源量計算
　　利用山區土壤水資源模型計算，在電子表格中分別輸入分區計算單元日降雨過程資料和日水面蒸發量，并根據各因子，利用建立的土壤水資源量（mm）、無效降水量、產生徑流量(mm)公式進行計算。
　　河北省太行山區耕地總面積為4473.15km2。經計算1980-2000年土壤水資源總量為312.360億m3；土壤水資源評價區降水總量為508.295億m3；無效降水總量為27.908億m3；徑流總量168.027億m3。土壤水資源量、徑流量分別占評價區降水量的61.45%和33.05%。
　　說明土壤水資源的總量比地表水與地下水之和還要多。正是因為有了土壤水資源的存在，才滿足了作物生長發育的大部分需水要求，占太行山區耕地80%以上的旱地農作物，完全靠土壤水生長發育，而且土壤水是農、林、牧業和陸地生態環境最重要的水資源。土壤水資源實際上在很大程度上對地表淡水和地下淡水的使用起到了很大的替代作用。
　　3.1土壤水空間變化規律
　　由圖3-1可看出，1980-2000年太行山耕地土壤水資源量與逕流量變化幅度與土層厚度關系密切,土層越薄土壤水資源量與逕流量變化幅度越小,產生逕流量越大,土層越厚土壤水資源量與逕流量變化幅度越大,產生逕流量越小，可見土層厚度是決定土壤蓄水能力的主要因素。由于不同土質土壤儲水能力不同，因此在同一土層厚度范圍內各土質土壤水資源變化規律不同，表現為砂土>砂壤土>壤土>粘壤土>粘土。
　　                                
　　
　　太行山行政分區多年平均土壤水資源量表
　　表3-1              

                   　　　　表3-1可以看出太行山區耕地多年平均土壤水比較豐富的地區為永年、井陘礦區，土壤水資源量達到400mm以上；而涉縣、贊皇、行唐、淶水土壤水資源量比較少年平均量在300mm以下。
　　3.2壤水時間變化規律
　　經過不同質地典型單元豐、平、枯三個典型年進行逐月土壤儲水量初步分析，不同土質土壤典型年土壤儲水量年內變化有如下規律：土質為砂土時，各典型年1-5月份儲水量基本一致，6-12月份儲水量豐水年>平水年>枯水年;土質為其它類型時，1-5月份豐水年與平水年儲水量基本一致，枯水年則比豐、平水年偏小較多，6-12月份儲水量仍有豐水年>平水年>枯水年的規律。同時。土壤水資源量的年內分配規律表現為土壤水資源量與降水量年內分配同頻率。
　　4．土壤水巨大的開發利用潛力
　　土壤水資源被使用則對農作物生長和生態環境質量改善發揮作用。如果不被農作物、喬灌草所使用，則土壤水也同樣通過蒸發作用白白消耗掉。所以，科學使用土壤水資源是節約水資源，發揮土壤水對地表水和地下水資源替代作用的最有效措施。
　　太行山區多年平均降水量為556mm，平均土壤水資源量為349mm，大氣降水轉化為土壤水的比例達62.8%。說明土壤水資源十分豐富。但土壤水不能長期儲存，隨每次降水量和降水頻率而增大，降水后如得不到使用，又因土壤水的蒸發而逐漸減少，土壤水資源處于不斷的動態變化之中。冬春季降水總量少，降水頻率低，所以，土壤水也較少，掌握了上述規律，有效調整農作物種植結構，降低了農業生產成本，節約了水資源，大幅度提高了單位土地耕地面積的產量。
　　在年土壤水資源量達到350~400mm以上的區域，遵循降水和土壤水季節變化的規律，使各個區域的農業的產業結構和作物的產品結構和土壤水最優的周期相適應，即增加適水性作物的種植面積，減少非適水的作物的種植面積，充分發揮土壤水資源潛力，這樣，既可以減少灌溉成本，又可達到成倍地提高水分生產率和農產品產量的目的。針對我省太行山地區的實際情況說，增加夏種秋收的棉花和雜糧的種植面積，減少冬小麥的種植面積，就能達到充分有效地利用土壤水，使土壤水資源潛力得到最大限度的發揮，減少灌溉用水量，緩解水資源供水矛盾的目的。
　　參考文獻:
　　1.靳孟貴，張人權，等.土壤水資源評價的研究[J].水利學報，1999，8：30～34
　　2夏自強，李瓊芳.土壤水資源及其評價方法研究[J].水科學進展，2001，（4）：535～540
　　3由懋正，王會肖.農田土壤水資源評價[M].北京：氣象出版社，1996
　　4楊詩秀，雷志棟，等．農田尺度土壤水分監測[J]．水科學進展，1996，7（1）：14～19