본문 바로가기
정보를 알려주는 블로그 정보를 알려주는 블로그

Darcy 법칙으로 지반 내 지하수 침투속도와 유량 계산하기

읽는 시간 약 21분

지하수는 지반 내부의 공극과 암반의 균열을 따라 이동한다. 이때 물의 흐름을 정량적으로 평가하는 데 사용되는 대표적인 관계식이 Darcy 법칙(Darcy’s Law)이다.

Darcy 법칙은 지반의 투수계수와 수리경사, 물이 통과하는 단면적을 이용해 지하수 유량을 계산하는 방법이다.

흙막이 굴착 시 발생하는 지하수 유입량을 검토하거나, 제방 내부의 침투 흐름을 평가하고, 지하수의 평균 이동속도를 추정할 때 기본적인 개념으로 활용된다.

그러나 Darcy 법칙을 사용할 때에는 중요한 구분이 필요하다.

Darcy 유속은 지반 전체 단면적을 기준으로 계산한 유량속도이며, 실제 공극을 따라 이동하는 물의 평균 속도와 동일하지 않다.

이번 글에서는 Darcy 법칙의 원리와 공식, 수리경사 및 투수계수의 의미, 지하수 유량과 침투속도 계산 방법을 실제 수치 예제와 함께 알아본다.

1. Darcy 법칙(Darcy’s Law)이란 무엇인가?

Darcy 법칙은 다공성 매질을 통과하는 물의 흐름이 수리경사와 투수계수에 비례한다는 관계를 나타낸다.

1856년 프랑스의 기술자 앙리 다르시(Henri Darcy)가 모래층을 통과하는 물의 흐름을 실험적으로 연구하면서 제시한 법칙이다.

일반적인 1차원 포화 흐름에서 Darcy 법칙은 다음과 같이 표현한다.

Q = k × i × A

  • Q: 지하수 유량(m³/s)
  • k: 투수계수 또는 수리전도도(m/s)
  • i: 수리경사(무차원)
  • A: 흐름 방향에 수직인 전체 단면적(m²)

이 식은 지반의 투수계수가 크거나 수리경사가 커질수록, 동일한 단면적을 통과하는 유량이 증가한다는 것을 보여준다.

다만 Darcy 법칙은 기본적으로 층류 흐름을 가정한다.

또한 실제 지반에서는 투수계수가 위치와 방향에 따라 달라질 수 있으므로, 단순한 1차원 계산 결과를 복잡한 현장에 그대로 적용해서는 안 된다.

2. Darcy 법칙의 기본 공식과 부호의 의미

Darcy 법칙은 수두의 변화 방향까지 고려하면 다음과 같이 표현할 수 있다.

q = -k × (dh/dl)

  • q: Darcy 유속 또는 비유량(m/s)
  • k: 투수계수(m/s)
  • h: 전수두(m)
  • l: 흐름 방향의 거리(m)

식 앞의 음수 기호는 물이 일반적으로 전수두가 높은 곳에서 낮은 곳으로 이동한다는 것을 의미한다.

실무에서 흐름 방향이 이미 정해져 있고 유량의 크기만 계산하는 경우에는 다음과 같이 사용할 수 있다.

q = k × i

Q = q × A = k × i × A

여기서 i는 수두 차이를 흐름 길이로 나눈 값의 크기다.

따라서 계산할 때에는 수리경사에 부호를 포함했는지, 절댓값을 사용했는지 구분해야 한다.

3. 수리경사(Hydraulic Gradient)란 무엇인가?

수리경사는 물이 이동하는 구간에서 단위 길이당 발생하는 수두 차이를 의미한다.

일반적인 1차원 계산에서는 다음과 같이 구한다.

i = Δh / L

  • i: 수리경사
  • Δh: 두 지점 사이의 전수두 차이(m)
  • L: 흐름 경로의 길이(m)

예를 들어 두 지점 사이의 전수두 차이가 2m이고 흐름 길이가 10m라면 다음과 같다.

i = 2 / 10

i = 0.2

수리경사는 길이를 길이로 나눈 값이므로 단위가 없다.

여기서 주의할 점은 수리경사가 단순히 지표면의 경사와 같은 개념이 아니라는 것이다.

지하수는 지표면의 높낮이가 아니라 전수두의 차이에 따라 흐른다.

4. 전수두(Total Hydraulic Head)의 구성

지하수 흐름을 이해하려면 전수두의 개념을 알아야 한다.

일반적인 지반 내 침투 해석에서는 전수두를 다음과 같이 표현한다.

h = z + u/γw

  • h: 전수두(m)
  • z: 위치수두(m)
  • u: 간극수압(kPa 또는 kN/m²)
  • γw: 물의 단위중량(kN/m³)

위 식은 속도수두를 무시할 수 있는 일반적인 지반 침투 조건에서 사용한다.

위치수두는 기준면으로부터의 높이를 나타낸다.

압력수두는 간극수압을 물의 단위중량으로 나눈 값이다.

예를 들어 어떤 지점의 위치수두가 3m이고 압력수두가 5m라면 전수두는 다음과 같다.

h = 3 + 5

h = 8m

서로 다른 지점의 전수두를 비교하면 지하수의 흐름 방향을 판단할 수 있다.

5. 투수계수(k)가 의미하는 것은 무엇인가?

투수계수는 지반이 물을 얼마나 잘 통과시키는지를 나타내는 지표다.

일반적으로 공극이 크고 서로 잘 연결된 지반은 높은 투수계수를 나타낼 수 있다.

반면 작은 공극이 발달한 점성토는 낮은 투수계수를 나타내는 경우가 많다.

다만 투수계수는 흙 입자의 크기만으로 결정되지 않는다.

다음과 같은 조건의 영향을 받는다.

  • 입도분포와 세립분 함량
  • 간극비와 다짐 상태
  • 공극의 연결성
  • 흙의 구조와 층리
  • 물의 점성과 온도
  • 포화 상태

또한 같은 지반에서도 수평방향 투수계수와 수직방향 투수계수가 다를 수 있다.

따라서 지층의 방향성을 고려하지 않고 하나의 투수계수만 적용하면 계산 결과에 오차가 발생할 수 있다.

6. 토질별 투수계수의 일반적인 경향

토질에 따라 투수계수는 매우 넓은 범위에서 달라질 수 있다.

토질일반적인 투수 특성주요 특징
자갈매우 높은 편큰 공극을 통한 빠른 흐름 가능
모래높은 편입도와 세립분 함량에 따라 차이
실트낮거나 중간 정도입도와 구조에 따라 변동
점토매우 낮은 편작은 공극과 낮은 투수성
균열이 있는 암반균열 상태에 따라 크게 변동균열의 연결성이 중요

이 표는 일반적인 경향을 설명한 것이다.

실제 투수계수는 같은 토질에서도 수십 배 이상 차이가 날 수 있으므로, 설계에서는 시험자료와 현장 조건을 우선해야 한다.

특히 균열이 발달한 점토나 암반은 흙 입자 자체의 투수성이 낮더라도 균열을 따라 물이 빠르게 이동할 수 있다.

7. Darcy 법칙을 이용한 지하수 유량 계산 예제

다음은 Darcy 법칙을 이용해 지하수 유량을 계산하는 가상 사례다.

주어진 조건

항목적용값
투수계수(k)1.0 × 10⁻⁴m/s
수두 차이(Δh)2m
흐름 길이(L)10m
흐름 단면적(A)5m²

① 수리경사 계산

i = Δh / L

i = 2 / 10

i = 0.2

② Darcy 유속 계산

q = k × i

q = 1.0 × 10⁻⁴ × 0.2

q = 2.0 × 10⁻⁵m/s

③ 유량 계산

Q = q × A

Q = 2.0 × 10⁻⁵ × 5

Q = 1.0 × 10⁻⁴m³/s

④ 하루 유량으로 환산

하루는 86,400초이므로 다음과 같이 계산한다.

Qday = 1.0 × 10⁻⁴ × 86,400

Qday = 8.64m³/day

따라서 주어진 조건에서 계산한 지하수 유량은 하루 약 8.64m³다.

단, 이 값은 일정한 투수계수와 수리경사를 갖는 1차원 정상류 조건을 가정한 결과다.

실제 굴착 현장의 지하수 유입량은 지층 구조와 경계 조건, 굴착 형상 등에 따라 달라질 수 있다.

8. Darcy 유속과 실제 침투속도의 차이

Darcy 법칙에서 계산하는 q는 일반적으로 Darcy 유속 또는 비유량(Specific Discharge)이라고 부른다.

이는 물이 통과하는 지반의 전체 단면적을 기준으로 계산한 값이다.

그러나 지반의 전체 단면적이 모두 물의 흐름에 사용되는 것은 아니다.

실제로 물은 연결된 공극을 통해 이동한다.

따라서 평균적인 공극 내 이동속도를 추정할 때에는 유효공극률을 고려한다.

vs = q / ne

또는

vs = (k × i) / ne

  • vs: 평균 선형유속(m/s)
  • q: Darcy 유속(m/s)
  • ne: 유효공극률(무차원)

여기서 유효공극률은 물의 이동에 실제로 기여하는 공극의 비율을 의미한다.

총공극률과 유효공극률은 동일하지 않을 수 있다.

9. 실제 평균 침투속도 계산 예제

앞선 계산에서 Darcy 유속은 다음과 같았다.

q = 2.0 × 10⁻⁵m/s

여기에 유효공극률이 0.25라고 가정해 보자.

vs = q / ne

vs = 2.0 × 10⁻⁵ / 0.25

vs = 8.0 × 10⁻⁵m/s

이를 하루 이동거리로 환산하면 다음과 같다.

vs,day = 8.0 × 10⁻⁵ × 86,400

vs,day = 6.912m/day

따라서 계산된 평균 선형유속은 하루 약 6.91m다.

구분계산 결과
Darcy 유속2.0 × 10⁻⁵m/s
평균 선형유속8.0 × 10⁻⁵m/s
Darcy 유속의 하루 환산값1.728m/day
평균 선형유속의 하루 환산값6.912m/day

두 값이 다른 이유는 전체 단면적과 실제 물이 이동하는 유효공극 공간이 다르기 때문이다.

다만 평균 선형유속 역시 이상화된 평균값이다.

개별 물 입자의 실제 이동 경로는 구불구불할 수 있으며, 지층의 불균질성과 분산 등의 영향을 받는다.

10. Darcy 법칙을 이용한 지하수 이동시간 계산

평균 선형유속을 이용하면 일정 거리의 지하수 이동시간을 개략적으로 계산할 수 있다.

t = L / vs

  • t: 이동시간
  • L: 이동거리
  • vs: 평균 선형유속

앞선 예제에서 평균 선형유속이 하루 6.912m이고 이동거리가 100m라고 가정하자.

t = 100 / 6.912

t ≈ 14.47일

따라서 단순한 1차원 정상류 조건에서는 약 14.5일의 이동시간이 계산된다.

하지만 실제 오염물질의 이동시간은 이 값과 다를 수 있다.

오염물질은 지반 내부에서 흡착, 확산, 분산 및 화학적 반응 등의 영향을 받을 수 있기 때문이다.

따라서 이 계산값을 오염물질이 실제로 도착하는 시간을 확정하는 근거로 사용해서는 안 된다.

11. 정수위 투수시험(Constant Head Test)

정수위 투수시험은 시험 중 일정한 수두 차이를 유지하면서 시료를 통과하는 유량을 측정하는 방법이다.

일반적으로 비교적 투수성이 높은 조립토에 활용된다.

정수위 조건에서는 다음 식을 이용해 투수계수를 계산할 수 있다.

k = (V × L) / (A × Δh × t)

  • V: 측정 시간 동안 통과한 물의 체적(m³)
  • L: 시료 길이(m)
  • A: 시료 단면적(m²)
  • Δh: 수두 차이(m)
  • t: 측정 시간(s)

계산 예제

항목적용값
통과한 물의 체적(V)0.001m³
시료 길이(L)0.20m
단면적(A)0.01m²
수두 차이(Δh)0.50m
측정 시간(t)400초

공식에 대입하면 다음과 같다.

k = (0.001 × 0.20) / (0.01 × 0.50 × 400)

k = 1.0 × 10⁻⁴m/s

따라서 계산된 투수계수는 1.0 × 10⁻⁴m/s다.

이 값은 설명용 가상 시험자료를 이용한 결과다.

실제 시험에서는 시료의 포화 상태와 측정 장비, 시험 절차 등을 관련 기준에 따라 관리해야 한다.

12. 변수위 투수시험(Falling Head Test)

변수위 투수시험은 시간에 따라 수두가 감소하는 과정을 측정하여 투수계수를 계산하는 방법이다.

일반적으로 정수위 시험보다 투수성이 낮은 흙을 평가할 때 활용할 수 있다.

기본적인 변수위 시험의 계산식은 다음과 같다.

k = [a × L / (A × t)] × ln(h1 / h2)

  • a: 스탠드파이프의 단면적(m²)
  • L: 시료 길이(m)
  • A: 시료 단면적(m²)
  • t: 측정 시간(s)
  • h1: 초기 수두 차이(m)
  • h2: 최종 수두 차이(m)

여기서 ln은 자연로그를 의미한다.

이 식은 일정한 시료 단면적과 적절한 경계 조건 등을 가정한 기본적인 변수위 시험 관계식이다.

실제 시험에서는 장비 구성과 적용 표준에 맞는 계산식을 사용해야 한다.

13. 정수위 시험과 변수위 시험의 차이

구분정수위 투수시험변수위 투수시험
수두 조건일정하게 유지시간에 따라 감소
주요 측정값통과 유량과 시간수두 변화와 시간
일반적인 활용비교적 투수성이 높은 흙비교적 투수성이 낮은 흙
계산 특징일정 수두 차이를 이용시간에 따른 수두 변화를 이용
주의사항안정적인 흐름 유지정확한 수두·시간 측정

두 시험은 투수계수를 평가하기 위한 방법이지만, 모든 흙에 동일하게 적용할 수 있는 것은 아니다.

특히 매우 낮은 투수계수를 가진 재료에서는 유연벽 투수시험과 같은 적절한 시험방법을 검토할 필요가 있다.

14. Darcy 법칙이 성립하기 위한 조건

Darcy 법칙은 다양한 지반 침투 문제에 활용되지만 모든 흐름 조건에서 정확하게 성립하는 것은 아니다.

① 층류 흐름

Darcy 법칙은 기본적으로 점성력이 지배적인 층류 흐름에 적용된다.

매우 큰 유속이 발생하는 조립질 지반 등에서는 비Darcy 흐름이 나타날 수 있다.

② 적절한 포화 조건

기본적인 포화 투수계수는 물로 포화된 지반에서의 흐름 특성을 나타낸다.

불포화 지반에서는 함수 상태와 흡수력에 따라 수리전도도가 달라질 수 있다.

③ 투수계수의 적용 범위

투수계수는 지반의 상태와 유체의 성질에 영향을 받는다.

따라서 다른 온도나 포화 상태에서 얻은 값을 그대로 적용하면 오차가 발생할 수 있다.

④ 흐름의 차원과 경계 조건

단순한 1차원 공식은 흐름 방향과 단면적을 일정하게 가정한다.

실제 굴착 현장처럼 2차원 또는 3차원 흐름이 발생하는 경우에는 추가적인 침투 해석이 필요할 수 있다.

15. Darcy 법칙과 불포화 지반의 수분 이동

Darcy 법칙은 포화 지반의 지하수 흐름을 설명하는 기본 관계식이지만, 불포화 지반의 수분 이동을 해석할 때에도 확장된 형태로 활용된다.

불포화 지반에서는 공극의 일부에 공기가 존재한다.

따라서 수분 함량과 흡수력에 따라 수리전도도가 변화할 수 있다.

일반적으로 불포화 상태에서는 포화 상태보다 수리전도도가 낮아지는 경향이 나타난다.

이러한 수분 이동을 해석할 때에는 토양수분특성곡선(SWCC)과 불포화 수리전도도 함수 등을 고려할 수 있다.

따라서 포화 투수계수 하나만으로 모든 불포화 지반의 수분 이동을 설명하는 것은 적절하지 않다.

16. Darcy 법칙과 모관상승의 관계

모관상승은 물의 표면장력과 흙의 공극 구조 등에 의해 발생하는 현상이다.

반면 Darcy 법칙은 수두 차이와 수리전도도를 이용해 물의 흐름을 설명하는 관계식이다.

두 개념은 서로 관련되지만 동일한 것은 아니다.

예를 들어 지하수면 위의 불포화 지반에서 수분이 이동할 때에는 모관흡수력과 수분 보유 특성이 중요한 역할을 할 수 있다.

이때 수분 이동량을 평가하려면 불포화 수리전도도와 수두 분포 등을 함께 고려해야 한다.

따라서 모관상승 높이만으로 지하수 유량을 계산하거나, 포화 Darcy 법칙만으로 모관상승 현상을 모두 설명해서는 안 된다.

17. 흙막이 굴착과 제방에서의 활용

Darcy 법칙은 다양한 지반공학적 문제에 활용된다.

① 흙막이 굴착

굴착면 주변의 지하수 흐름과 유입량을 검토할 때 활용할 수 있다.

다만 굴착 형상과 지층 분포가 복잡한 경우에는 2차원 또는 3차원 침투 해석이 필요할 수 있다.

② 제방 및 흙댐

제방 내부의 침투 흐름과 배수 조건을 검토하는 데 활용된다.

침투 유량뿐 아니라 동수경사와 간극수압 분포도 중요한 검토 대상이다.

③ 지하수 오염 이동

지하수의 평균 이동속도와 이동시간을 개략적으로 추정하는 데 활용할 수 있다.

다만 실제 오염물질의 이동은 흡착과 분산 등의 영향을 함께 받는다.

④ 배수시설 설계

지반의 투수성과 수두 조건을 이용해 배수 성능을 검토할 수 있다.

이때 실제 배수시설의 형상과 경계 조건을 반영해야 한다.

18. Darcy 법칙을 이용할 때 흔히 발생하는 계산 오류

오류 1. Darcy 유속을 실제 물의 이동속도로 해석한다.

Darcy 유속은 전체 단면적을 기준으로 계산한 값이다.

평균 선형유속을 추정할 때에는 유효공극률을 고려해야 한다.

오류 2. 수리경사를 지표면 경사와 동일하게 본다.

수리경사는 전수두 차이를 흐름 길이로 나눈 값이다.

오류 3. 투수계수 단위를 확인하지 않는다.

투수계수는 m/s 또는 cm/s 등으로 표시될 수 있다.

계산 전 반드시 단위를 통일해야 한다.

오류 4. 유량과 유속을 혼동한다.

유속은 길이/시간 단위를 사용하고, 유량은 체적/시간 단위를 사용한다.

오류 5. 모든 지반에 동일한 투수계수를 적용한다.

실제 지반에서는 깊이와 위치, 흐름 방향에 따라 투수계수가 달라질 수 있다.

오류 6. Darcy 법칙의 적용 한계를 무시한다.

조립질 지반에서 유속이 매우 높아지는 경우 등에는 비Darcy 흐름을 검토할 필요가 있다.

19. 지하수 유량 계산을 위한 체크리스트

검토항목확인 내용
투수계수시험값과 적용 단위
수리경사전수두 차이와 흐름 길이
단면적흐름 방향에 수직인 면적
포화 상태포화 또는 불포화 조건
흐름 방향수평·수직 및 복합 흐름
유효공극률평균 선형유속 계산 시 필요
지층 구조불균질성과 이방성
경계 조건지하수위와 배수 조건
흐름 형태층류 또는 비Darcy 흐름 가능성
계산 목적유량·유속·이동시간 중 무엇을 구하는지

20. 자주 묻는 질문

Q1. Darcy 법칙의 기본 공식은 무엇인가요?

일반적인 1차원 포화 흐름에서 유량의 크기는 Q = k × i × A로 계산할 수 있다.

Q2. Darcy 유속과 실제 침투속도는 같은 값인가요?

아니다. Darcy 유속은 전체 단면적을 기준으로 계산한 값이며, 평균 선형유속은 유효공극률을 고려해 추정할 수 있다.

Q3. 투수계수가 크면 지하수 유량도 증가하나요?

다른 조건이 동일하다면 투수계수가 클수록 계산되는 유량은 증가한다.

Q4. Darcy 법칙은 점토에도 적용할 수 있나요?

적절한 포화·층류 조건에서는 점토에도 적용할 수 있다. 다만 매우 낮은 투수성으로 인해 시험과 측정에 더 많은 시간과 정밀한 장비가 필요할 수 있다.

Q5. Darcy 법칙으로 지하수의 이동시간을 계산할 수 있나요?

평균 선형유속과 이동거리를 이용해 개략적인 이동시간을 계산할 수 있다. 그러나 실제 이동 경로와 지층의 불균질성 등에 따라 결과가 달라질 수 있다.

21. 결론: Darcy 법칙을 지반공학에서 어떻게 활용해야 할까?

Darcy 법칙은 지반 내부의 지하수 흐름을 정량적으로 평가하는 기본적인 관계식이다.

투수계수와 수리경사, 흐름 단면적을 이용하면 일정한 조건에서 지하수 유량을 계산할 수 있다.

또한 유효공극률을 고려하면 평균 선형유속과 이동시간을 개략적으로 추정할 수 있다.

그러나 Darcy 유속과 실제 공극 내 평균 이동속도는 동일하지 않으며, 지반의 불균질성과 이방성, 포화 상태 및 경계 조건에 따라 실제 흐름은 달라질 수 있다.

따라서 지반공학적 해석에서는 Darcy 법칙의 기본 원리를 이해하고, 투수계수와 수리경사 및 적용 조건을 정확하게 검토하는 것이 중요하다.


참고자료 및 출처

1. U.S. Geological Survey (USGS)

Ground-Water Flow and Darcy’s Law.

Darcy 법칙, 수리전도도, 수리경사 및 평균 선형유속에 관한 참고자료.

https://pubs.usgs.gov/circ/2002/circ1224/html/understanding.html

2. ASTM International — ASTM D2434-26

Standard Test Methods for Laboratory Measurement of Hydraulic Conductivity of Coarse-Grained Soils.

포화된 조립토의 수리전도도 측정에 관한 시험 표준.

https://store.astm.org/standards/d2434

3. ASTM International — ASTM D5084-24

Standard Test Methods for Measurement of Hydraulic Conductivity of Saturated Porous Materials Using a Flexible Wall Permeameter.

유연벽 투수시험 장치를 이용한 포화 다공성 재료의 수리전도도 측정에 관한 시험 표준.

https://store.astm.org/standards/d5084

4. Freeze, R. A. & Cherry, J. A. (1979)

Groundwater. Prentice-Hall.

지하수 흐름, Darcy 법칙 및 다공성 매질의 수리학적 특성에 관한 전문 참고문헌.

※ 본문의 계산 예제는 이해를 돕기 위한 가상 조건입니다. 실제 설계에는 현장 조사자료와 적용 기술기준을 우선해야 합니다.

지반인사이트
지반공학 및 토목공학 분야의 기술정보를 전문적으로 정리하고 있습니다. 지반조사, 토질시험, 연약지반 개량, 흙막이 공법, 기초공학, 지반계측 및 사면 안정 등 다양한 주제를 다룹니다. 국내외 건설기준과 기술문헌을 바탕으로 공법의 작동 원리, 시공방법, 계산 공식 및 현장 적용 시 주의사항을 이해하기 쉽게 설명하는 것을 목표로 합니다. 본 사이트는 지반공학 관련 기술자료를 체계적으로 제공하는 정보형 웹사이트입니다.
함께 보면 좋은 글

댓글 0

첫 댓글을 남겨보세요.

error: Content is protected !!

광고 차단 알림

광고 클릭 제한을 초과하여 광고가 차단되었습니다.

단시간에 반복적인 광고 클릭은 시스템에 의해 감지되며, IP가 수집되어 사이트 관리자가 확인 가능합니다.