Kriging기법과 수치모형에 의한 이안지구 대수층의 투수량계수 추정

Aquifer Transmissivity Estimation with Kriging Techniques and Numerical Model in the LAN

대수층을 취급하는데 있어서 어려운 문제중에 하나는 수문학적인 변수의 공간적인 분포를 추정하는 것이다. 이안지구 대수층의 경우 면밀한 분석을 위한 관정시험 자료가 많지 않다. 이안지구 대수층의 투수량계수의 분포를 결정하기 위하여, 수위관측지점을 투수량계수 자료로 전환하기 위해서는 이 지점을 pilot point로 정하여 수문학적인 변수의 공간분포를 추정하는 3단계의 절차를 검토하였다. 구성된 절차는 기지점의 값을 기초로 하여 변수분포의 지질통계적인 평가, 수치모델에 의하여 산정된 값과 비교, 최소화 알고리즘에 의하여 pilot point값을 보정하는 것으로 이루어 진다. 이안대수층의 변수분포를 평가 하는데 총 258개의 소구역으로 분할하고 각 소구역의 대표 투수량계수를 대상으로 Variogram을 구해본 결과 w=0.623, a=2.743을 얻었다. 이 계수를 이용하여 Kriging을 수행하고 보간된 투수량계수 분포를 가지고 수치해를 얻고 이를 관측치와 비교하여 보정하는 절차를 통해 투수량계수의 등고선을 나타냈다.

One of the delicate problems in aquifer management is the identification of the spatial distribution of tile hydrological parameters. The observed data are insufficient to identify the distribution of transmissivities in LAN aquifer. To determine the distribution of the transmissivity in LAN aquifer, it would be required to transform the observed heads at the pilot points into transmissivities. Therefore, three procedures wire tackled for the identification of the spatial distribution of the hydrological parameters; geostatistical estimate of the parameter field on the basis of known well point, heads reconstructed by a numerical model, and modification of the values at pilot points by a minimization algorithm. The variogram of Kriging has been applied to a total of 258 transmissivity value in attempt to quantify their distribution of LAN aquifer. Variogram of the observed and optimized transmissivities at pilot points are adapted to the exponential form. So, it is fitted by theoretical one with coefficients of w=0.623, a=2.743. Values of head obtained through numerical analysis are adjusted to the observed values so that heads have been transformed completely into the transmissivities at the observation wells. The procedure represented contour map of the estimated transmissivities and the calculated head.