1多点源滴条件下土壤水分运动的数学模型.docVIP

1 多点源滴条件下土壤水分运动的数学模型 1.1土壤水分运动方程 在一些实际的滴灌工程中，由于滴头间距较小，相邻两滴头之间的湿润锋经常会发生交汇现象，从而形成多孔介质，不考虑空气及温度对土壤水分运动的影响，也不考虑土壤水分运动的滞后效应，则土壤水分运动可以用Richard方程来描述 1 式中：—土壤体积含水率，cm3/cm3;h—土壤负压水头，cm;k h —土壤非饱和导水率，cm/min;t—时间，min;x、y、z—平面坐标，cm;规定z向下为正。 1.2定解条件 1）初始条件 H x、y、z h0 x、z 0≤x≤X,0≤z≤Z,T 0 2 式中：h0（x,z）—计算区域的初始土壤负压水头，cm;X—模拟计算区域的垂向最大距离，取40cm。 假定计算区域的初始土壤负压水头均匀分布，可以通过土壤水分特征曲线，将土壤的初始含水率转换为可用来作数值模拟时的土壤水分初始负压剖面。 2）边界条件 图1给出了多点源滴灌条件下土壤水分运动数值模拟计算区域示意图。图1中的A（0，0）及B（X，0）分别为2个滴头的位置。 AB边界为上边界，上边界是一个运动的边界，比较复杂，分2种情况考虑，一是滴灌流量较小，当滴头的出流速率小于土壤的入渗速率时，滴头所供应的水分能在瞬间渗入到土壤中，地表不会形成积水，即：地表无积水情况；二是，当滴头的出流速率大于土壤的入渗速率时，滴头所供应水分不能迅速渗入到土壤中，在地表形成积水，即:地表形成积水情况。 （1）地表无积水 在滴头处 ） （3） 式中：E一土壤表面蒸发强度，cm/min; q q0/A,q0滴头流量，mL/min;A一土壤表面湿润区面积，cm2。 在上边界非滴处 （0﹤x﹤X,z 0） 4 （2）地表积水 在地表饱和区 h 0 0≤x≤Rs,X-RS≤z≤X,z 0 5 在地表非饱和去 （Rs﹤x﹤X-Rs,z 0） 6 式中Rs一地表积水半径，cm。 在地表积水情况下，关于地表积水半径准确确定问题，目前仍然是个难点。本文借鉴文献【18】中的方法：先假定地表积水半径为一定值，用数值方法来模拟土壤水分运动（模拟的水量平衡误差为0.5%左右），将土体增加水量与时间的比值近似为滴头流量，然后通过一系列的模拟数据来拟合地表积水半径与滴头流量之间的关系。 CD边界为下边界，考虑地下水埋深的较大情况，可认为土壤水分和土壤压力水头保持不变，即 h x,z,t h0 0≤x≤X,z Z,T 0 7 AC边界和BD边界分别为左、右边界，由于其对称原因，左右边界上的水分通量为0