[发明专利]地下水在不同渗透介质边界渗流量的精确计算方法

权利要求书

1.一种地下水在不同渗透介质边界渗流量的精确计算方法，其特征在于，包括步骤： 

A、在已知均衡子区的前提下，以子区内部某边界渗流量为参数，参数化子区内部所有边界渗流量值； 

B、利用参数化的边界渗流量值，选择流速插值算法，参数化子区内所有单元流速场； 

C、以单元边界切向流速差最小为条件，最优化求解子区内所有边界渗流量。 

2.根据权利要求1所述的地下水在不同渗透介质边界渗流量的精确计算方法，其特征在于，所述步骤A包括： 

A1、根据收集到的水文地质资料，建立二维有限单元或控制体积有限单元模型并求解地下水预报问题，获得各结点的预报水位值； 

A2、建立基础数据的拓扑结构，包括对以下基础数据的存在和相互调用的结构：有限单元法的网格单元、地下水预报问题的预报水位值、边界条件、垂向补排、井参数； 

A3、利用Codes-Wolfgang或Durlofsky方法获得均衡子区数据，将由Durlofsky获得的子区数据转换为类似Codes-Wolfgang子区的形式； 

A4、假设子区内某一边界渗流量值为Q，并以此参数化所有未知渗流量边界。 

3.根据权利要求2所述的地下水在不同渗透介质边界渗流量的精确计算方法，其特征在于，所述步骤B还包括： 

利用参数化的边界插值单元内部流速场，参数化所有边界的切向流速差。 

4.根据权利要求3所述的地下水在不同渗透介质边界渗流量的精确计算方法，其特征在于，所述步骤C包括： 

C1、最优化边界切向流速差，计算出最优解或最优区间； 

C2、如果步骤C1中获得的是最优区间，则选取其中使得子区旋度最小的值为最优解； 

C3、将最优解代入参数化公式，得到所有未知渗流量边界的渗流量值； 

C4、针对其它任意子区重复进行上述步操A4之后的步骤。