양수시험

유체역학 시뮬레이션 수치해석 마지막 수정: 2026-09-05 04:34:09

1. 개요[편집]

양수시험 Pumping test / Aquifer test
목적	현장에서 투수량계수 T와 저류계수 S를 구한다
측정량	양수량 Q(일정) · 시간에 따른 수위강하 s(r, t)
기본 해	타이스 해(1935) — 우물함수 W(u), u = r²S/4Tt
실무 해석	쿠퍼-야콥 반로그 직선법(1946)
확장	한투시-야콥(누수) · 노이만(지연산출) · 한투시(부분관통)
부속 시험	단계양수시험(우물손실) · 회복시험 · 순간충격시험
귀착점	MODFLOW 보정의 사전정보 · PEST 역해석의 관측치

대수층의 물성을 재는 유일하게 정직한 방법은, 물을 뽑아 보고 땅이 어떻게 반응하는지 보는 것이다.

양수시험은 관정에서 일정한 유량으로 물을 뽑으면서 주변 관측공의 수위강하를 시간에 따라 기록하고, 그 곡선을 해석해를 맞춰 대수층의 투수량계수 $T$ 와 저류계수 $S$ 를 역산하는 현장 시험이다. 대수층시험(aquifer test)이라고도 한다. 지하수 유동 문서가 방정식과 해석해를 다룬다면, 이 문서는 그 해석해를 실제 현장 자료에 거꾸로 붙이는 절차와 그때 벌어지는 일들을 다룬다.

핵심 감각은 이렇다. 시험이 재는 것은 시료 하나의 물성이 아니라 수위강하 원추가 덮은 부피 전체의 유효 물성이다. 코어 시료의 투수시험이 손바닥만 한 영역을 재고 순간충격시험이 관정 주변 수 m 를 잰다면, 양수시험은 수백 m 반경을 잰다. 그래서 같은 지점에서 세 방법이 서로 다른 값을 주는 것은 오차가 아니라 각자 다른 부피를 평균한 결과다. 그리고 MODFLOW 격자 셀이 대변하는 부피에 가장 가까운 것이 양수시험이라는 사실이, 이 오래된 시험이 아직도 모델 보정의 1순위 입력인 이유다.¹

2. 시험 구성[편집]

양수정(pumping well) — 물을 뽑는 관정. 여기서 잰 수위강하에는 뒤에 나올 우물손실이 섞여 있어 $T$ 산정에 그대로 쓰면 안 된다.
관측공(observation well) — 최소 하나, 가능하면 서로 다른 거리·방향에 여러 개. 방향별로 다른 $T$ 가 나오면 이방성이나 경계의 신호다.
양수량 $Q$ — 시험 내내 일정해야 한다. 실무 허용 변동폭은 대개 ±5 % 이내이며, 이걸 못 지키면 해석해의 전제부터 무너진다. 수중펌프의 토출량은 수위가 내려가면 함께 변하므로 밸브로 잡아 준다.
기간 — 피압대수층 24~~72시간, 자유면대수층은 3~~7일. 자유면 쪽이 긴 이유는 뒤의 지연산출 때문이다. 며칠 동안 뽑아낸 물을 어디로 흘려보낼지도 시험 설계의 일부다.²
배경 보정 — 시험 전 며칠간 수위를 미리 기록해 지역 추세를 뺀다. 피압대수층은 기압 변화에 반응하므로(기압효율 0.2~0.75) 기압 자료로 보정해야 하고, 해안이나 큰 강 근처에서는 조석 성분까지 뺀다. 10 cm 짜리 수위강하를 해석하는 시험에서 기압 보정을 빼먹으면 신호보다 잡음이 크다.

3. 타이스 곡선 정합[편집]

무한·균질·등방·피압 대수층에 완전관통 우물을 놓고 순간부터 일정 $Q$ 로 뽑으면 수위강하는

s(r,t)=\frac{Q}{4\pi T}\,W(u),\qquad u=\frac{r^2S}{4Tt},\qquad W(u)=\int_u^{\infty}\frac{e^{-y}}{y}\,\mathrm{d}y

$W(u)$ 는 우물함수이고 수학적으로는 지수적분 $E_1$ 이다. 타이스(C. V. Theis, 1935)가 열전도 문헌에서 그대로 옮겨 온 해이며, 유도와 물리적 배경은 지하수 유동 문서에 있다.

원래의 해석법은 곡선 정합(type-curve matching)이다. $\log W$ 대 $\log(1/u)$ 의 표준곡선을 투명지에 그려 두고, 실측 $\log s$ 대 $\log t$ 그래프 위에 축을 평행하게 유지한 채 겹쳐 미끄러뜨려 가장 잘 맞는 자리를 찾은 뒤, 임의의 정합점에서 네 값 $(s, t, W, 1/u)$ 을 읽어

T=\frac{Q\,W(u)}{4\pi s},\qquad S=\frac{4T\,u\,t}{r^2}

로 계산한다. 로그-로그 축에서 곡선을 평행이동하는 행위가 곧 두 미지수 $T$ , $S$ 를 정하는 일이라는 것이 이 방법의 우아한 점이다 — 세로 이동이 $T$ , 가로 이동이 $S$ 에 대응한다.

4. 쿠퍼-야콥 직선법[편집]

$u$ 가 작아지면(먼 시간, 가까운 거리) $W(u)\approx-\gamma-\ln u$ 로 잘라 쓸 수 있고, 이것이 쿠퍼-야콥 근사(1946)다. 상용로그로 정리하면

s\approx\frac{2.303\,Q}{4\pi T}\log_{10}\frac{2.25\,T t}{r^{2}S}

$u<0.01$ 이면 오차 1 % 미만, $u<0.05$ 까지 늘려도 5 % 수준이다. 이 형태가 100년 가까이 실무를 지배한 이유는 계산이 쉬워서가 아니라 $s$ 대 $\log t$ 가 직선이기 때문이다. 반로그 모눈에 찍고 자를 대면 끝난다.

시간-강하 해석(관측공 하나, 시간축). 한 로그 사이클당 기울기 $\Delta s$ 에서 $T=2.303Q/(4\pi\Delta s)$ , 직선을 $s=0$ 까지 외삽한 절편 $t_0$ 에서 $S=2.25Tt_0/r^2$ .
거리-강하 해석(여러 관측공, 같은 시각). 기울기가 두 배가 되어 $T=2.303Q/(2\pi\Delta s)$ , 절편 $r_0$ 에서 $S=2.25Tt/r_0^{2}$ .

숫자를 하나 넣어 보자. $Q=1000\,\mathrm{m^3/일}$ , 로그 사이클당 $\Delta s=0.5\,\mathrm{m}$ 이면 $T=2.303\times1000/(4\pi\times0.5)\approx366\,\mathrm{m^2/일}$ . $r=50\,\mathrm{m}$ 관측공에서 $t_0=1\,$ 분 $=6.9\times10^{-4}\,$ 일이면 $S=2.25\times366\times6.9\times10^{-4}/50^2\approx2.3\times10^{-4}$ — 피압대수층의 전형적인 값이다.

여기서 두 계수의 신뢰도가 전혀 다르다는 점을 봐야 한다. $T$ 는 직선의 기울기에서 나오고, $S$ 는 절편을 외삽한 위치에서 나온다. 기울기를 몇 퍼센트만 잘못 잡아도 $t_0$ 는 배수로 움직이고, $S$ 는 그만큼 통째로 틀린다. 게다가 $S$ 는 $r^2$ 로 나뉘므로 관측공 거리 측량 오차가 제곱으로 들어간다. 그래서 실무 격언이 이렇다 — $T$ 는 믿고, $S$ 는 자릿수만 믿어라.

5. 조건을 하나씩 깨기[편집]

타이스 해의 전제는 여덟 개쯤 되고, 현장은 그중 최소 서너 개를 어긴다. 각 위반이 곡선 모양을 어떻게 바꾸는지가 곧 진단법이다.

누수 대수층 — 한투시-야콥 해(1955). 위나 아래의 반투수층을 통해 물이 새어 들어오면 수위강하가 무한정 커지지 못하고 유한한 값으로 수렴한다. 해는 $s=\frac{Q}{4\pi T}W(u,\,r/B)$ 이고, 누수인자 $B=\sqrt{Tb'/K'}$ 가 반투수층 물성을 역산하는 통로가 된다. 반투수층 자체의 저류까지 넣은 것이 한투시(1960)의 수정해다.

자유면 대수층 — 지연산출. 자유면에서는 공극이 실제로 배수되는 데 시간이 걸리므로 곡선이 S자 3단이 된다. 초기에는 탄성 저류( $S\sim10^{-4}$ )로 빠르게 내려가고, 중간에 중력 배수가 따라오면서 평평해지며, 후기에 비산출률 $S_y\sim0.1\!-\!0.3$ 을 반영한 완만한 직선으로 돌아온다. 볼튼(1963)의 지연산출 이론과 노이만(1972·1974)의 해가 표준이며, 중간 평탄 구간을 후기 직선으로 착각해 $T$ 를 잡으면 수 배가 틀린다. 자유면 시험을 며칠씩 하는 이유가 이 구간을 지나가기 위해서다.

부분관통. 우물이 대수층 두께의 일부만 관통하면 흐름이 3차원이 되어 우물 근처 수위강하가 커진다. 한투시(1961)의 보정을 쓰거나, 더 간단하게는 관측공을 충분히 멀리 둔다. 통상 기준은

r>1.5\,b\sqrt{\frac{K_h}{K_v}}

이며, 이 거리를 넘으면 부분관통 효과가 실질적으로 사라진다. 이방성이 심한 층상 대수층( $K_h/K_v=100$ )이면 이 거리가 두께의 15배까지 늘어난다는 점이 함정이다.

경계 — 영상정 이론. 무한 대수층 가정이 깨지는 순간이 곡선에 나타난다. 불투수 경계가 있으면 같은 자리에 가상의 양수정을 하나 더 놓은 것과 같아 직선의 기울기가 두 배가 되고, 하천 같은 정수두 경계가 있으면 가상의 주입정 효과로 곡선이 평평해진다. 반로그 그래프에서 직선이 꺾이는 시점과 그 방향이 경계의 존재·성격·거리를 알려 준다. 자를 대고 그리던 시절의 기술이지만, 지금도 이 꺾임을 먼저 보고 나서 수치 해석에 들어간다.

우물 자체 — 정호저류와 스킨. 아주 초기 자료는 관정 케이싱 안의 물이 빠지는 것(wellbore storage)만 보여 주며, 로그-로그 그래프에서 기울기 1의 직선으로 나타난다. 이 구간은 대수층 정보가 없으므로 버려야 한다. 시추 이수나 목적층 손상으로 생긴 스킨은 우물 근처의 국소 저항이며, 석유공학의 정호시험에서 넘어온 미분 진단(수위강하의 $\log t$ 미분을 함께 그리는 부르데 도해)이 지하수 쪽에도 정착했다. 미분 곡선의 평탄부가 곧 $T$ 이고, 누수·경계·부분관통이 그 위에 각각 다른 모양을 얹는다 — 눈으로 모형을 고르는 데 원곡선보다 미분 곡선이 훨씬 유용하다.

6. 단계양수시험과 회복시험[편집]

단계양수시험(step-drawdown test). 양수정 자체의 수위강하는 대수층 손실과 우물손실의 합이다.

s_w=BQ+CQ^{n},\qquad n\approx2

$BQ$ 는 층류인 대수층 손실, $CQ^n$ 은 스크린 통과와 케이싱 내 난류에서 오는 우물손실이다(Jacob, 1947). 유량을 3~5단계로 계단식으로 올려 각 단계의 준안정 수위강하를 재고, $s_w/Q$ 를 $Q$ 에 대해 찍으면 직선이 되어 절편에서 $B$ , 기울기에서 $C$ 가 나온다. 이 시험의 산출물은 대수층 물성이 아니라 우물 효율 $BQ/(BQ+CQ^2)$ 과 적정 양수량이며, 관정 노후화(스크린 폐색)를 추적하는 지표이기도 하다. 그리고 이 시험을 해 봐야 “양수정 수위강하로 $T$ 를 구하면 왜 안 되는지”가 숫자로 보인다.

회복시험(recovery test). 펌프를 끄고 수위가 돌아오는 과정을 재는 시험으로, 중첩 원리를 쓰면 잔류 수위강하는

s'=\frac{2.303\,Q}{4\pi T}\log_{10}\frac{t}{t'}

$t$ 는 양수 시작부터, $t'$ 는 정지부터의 시간이다. $s'$ 를 $\log(t/t')$ 에 대해 찍은 직선의 기울기에서 $T$ 가 나온다. 장점이 분명하다 — 펌프가 꺼져 있으므로 유량 변동도, 우물손실도, 펌프 진동도 없다. 그래서 회복시험의 $T$ 는 양수 구간의 $T$ 보다 대개 더 신뢰할 만하다. 대신 $S$ 는 나오지 않는다(정확히는 나오지만 믿을 수 없다). 양수시험을 설계할 때 회복 자료를 안 받고 철수하는 것은 공짜 검증 기회를 버리는 것이다.

7. 비유일성 — 이 시험의 진짜 난점[편집]

수위강하 곡선 하나에서 대수층을 복원하는 것은 대놓고 역문제이고, 그것도 잘 알려진 종류의 불량조건이다.

$T$ – $S$ 상관. 곡선 정합에서 $T$ 를 조금 키우고 $S$ 를 그에 맞춰 줄이면 초기 구간이 거의 그대로 겹친다. 목적함수의 골짜기가 길게 누워 있다는 뜻이고, 최적화로 풀면 상관계수 0.9 이상의 두 매개변수가 나온다.
모형 자체의 비유일성이 더 나쁘다. 후기 곡선이 평평해지는 현상은 (가) 누수, (나) 정수두 경계, (다) 자유면의 중력 배수 세 가지가 모두 만들 수 있다. 세 모형은 물리가 전혀 다른데 자료 적합도는 비슷하게 좋다. 어느 모형을 골랐느냐가 $T$ 를 두 배 이상 바꾸는데, 그 선택의 근거는 곡선이 아니라 지질 자료다. 곡선 정합의 결정계수를 자랑하는 보고서를 만나면 이 대목을 먼저 물어야 한다.
등방·균질 가정. 실제 얻는 $T$ 는 수위강하 원추가 덮은 영역의 어떤 평균이다. 2차원 로그정규 이질장에서 방사 흐름의 유효 투수량계수는 산술평균이 아니라 기하평균에 가깝다는 것이 알려져 있고, 이것이 코어 시험 평균보다 양수시험 값이 작게 나오는 흔한 이유다. 수리전도도 문서의 척도 의존성 논의와 같은 이야기.

8. 수치 역해석으로의 확장[편집]

해석해가 감당 못 하는 상황(다층, 이질, 여러 우물, 불규칙 경계, 시간에 따라 변하는 양수량)에서는 곡선 정합을 버리고 수치모형을 자료에 맞춘다. 국소 방사류 격자를 짠 MODFLOW 모형을 돌리고, 관측공의 수위강하 시계열을 관측치로 삼아 물성을 조정하는 것이다. 목적함수는 가중 잔차 제곱합

\Phi(\mathbf{m})=\sum_i w_i\left[s_i^{obs}-s_i^{sim}(\mathbf{m})\right]^2+\lambda\,\Phi_{reg}(\mathbf{m})

이고, 표준 도구가 PEST/PEST++ 다. 실무 절차는 대체로 이렇게 굳어져 있다.

매개변수를 파일럿 포인트로 두고 공간 보간은 크리깅으로 처리한다.
자코비안은 유한차분으로 채운다 — 매개변수 하나당 모형 실행 한 번이므로, 파일럿 포인트가 200개면 반복마다 201회 실행이다. 여기가 비용의 전부라 병렬 실행이 사실상 필수.
최적화는 레벤버그-마쿼트 방법, 불량조건은 티호노프 정규화나 SVD 절단으로 억누른다. 정규화항 $\Phi_{reg}$ 가 하는 일은 “자료가 말하지 않는 방향은 사전 추정값에 붙잡아 두는 것”이다.
불확실성은 선형 근사(자코비안 기반 사후 공분산)나 앙상블 계열(PESTPP-IES, 앙상블 칼만 필터의 반복형)로 낸다. 불확실성 정량화가 이 분야에서 유난히 일찍 자리 잡은 것은 취향이 아니라 필요였다.

주의할 것은 수치 역해석이 비유일성을 없애 주지 않는다는 점이다. 매개변수를 늘릴수록 적합은 좋아지고 예측 불확실성은 오히려 커진다. 그래서 성숙한 워크플로는 잔차를 0으로 만드는 대신, 양수시험의 해석해 결과를 사전 정보와 초기값으로 넣고, 정규화 강도를 자료 잡음 수준에 맞춰 잔차가 그 수준에서 멈추도록 조절한다. 자를 대고 그은 직선의 $T$ 가 아직도 쓸모 있는 이유가 여기 있다 — 그것이 최종 답이어서가 아니라, 수백 차원 최적화가 어디서 출발해야 하는지를 알려 주는 유일하게 저렴한 정보이기 때문이다.³

9. 자료를 볼 때의 점검표[편집]

로그-로그로 먼저 그린다. 초기 기울기 1 구간은 정호저류이므로 잘라낸다.
반로그로 다시 그려 직선 구간이 어디인지 본다. 꺾임이 있으면 경계·누수·지연산출 중 무엇인지 지질로 판정한다.
$u<0.05$ 를 만족하는 구간에서만 쿠퍼-야콥을 쓴다. 계산한 $T$ , $S$ 를 되넣어 $u$ 를 다시 계산해 자기모순이 없는지 확인한다.
여러 관측공의 $T$ 가 방향에 따라 계통적으로 다르면 이방성이나 선구조를 의심한다. 거리에 따라 계통적으로 다르면 이질성이다.
회복 자료로 $T$ 를 독립 산정해 양수 구간 값과 비교한다. 두 값이 2배 이상 벌어지면 유량이 일정하지 않았거나 자유면 배수가 진행 중이다.
마지막으로 물수지를 본다 — 뽑은 물의 총량과 수위강하 원추의 부피에 $S$ 를 곱한 값이 자릿수라도 맞는지. 검증 및 확인의 가장 싼 형태이고, 놀랄 만큼 자주 이 단계에서 오류가 잡힌다.

10. 관련 문서[편집]

이걸 “지지체적(support volume)이 맞는다”고 말한다. 격자 셀 한 변이 200 m 인 모형에 손바닥만 한 코어 시료의 $K$ 를 넣는 것은 단위가 아니라 척도가 안 맞는 일이며, 이질 매질에서 그 둘은 계통적으로 다른 값이다. 현장에서 “실험실 값이랑 왜 이렇게 다르냐”는 질문을 받으면 오차 이야기를 꺼내기 전에 이 이야기를 먼저 해야 한다. ↩
양수시험은 물리 시험 중에서도 유난히 사회적인 시험이다. 며칠간 지하수를 대량으로 뽑아 어딘가로 흘려보내야 하고, 이웃 관정의 수위가 실제로 내려가며, 그 사실이 곧 민원이 된다. 배출수를 시험 구역으로 되돌려 흘리면 함양 경계가 생겨 곡선이 평평해지는데, 이걸 누수로 해석해 반투수층 물성을 보고한 사례가 실제로 있다. 시험 설계에서 배수로 위치는 물리 문제다. ↩
그래서 이 분야 선배들이 신입에게 시키는 것이 아직도 모눈종이에 직선 긋기다. 자동 정합 소프트웨어가 뱉은 $T$ 가 손으로 그은 값과 자릿수부터 다르면 둘 중 하나는 틀렸고, 대개는 자동 쪽이 초기 정호저류 구간까지 통째로 적합에 넣은 경우다. 최적화는 사용자가 준 자료를 의심하지 않는다는 것이 이 도구의 유일한 단점이자 치명적인 단점이다. ↩