계산 방법 및 단계

(1) 수문지질학적 개념모델 확립

분석방법은 수리지질학적 조건에 대해 많은 제한을 두고 있으며 엄격한 이상화 요건을 가지고 있으나 실제 수문지질학적 조건은 매우 복잡한 경우가 많다. 분석적 방법을 사용하여 계산하기 위해서는 수문지질학적 조건을 합리적으로 단순화하고 일반화해야 하며, 단순화되고 일반화된 수문지질학적 조건을 수문지질학적 개념모델이라고 하며, 이는 지하수 시스템을 정성적으로 기술하는 것입니다.

1. 배수 유동장의 수리학적 특성을 분석합니다.

광상의 배수 유동장은 자연 유동장의 배경과 인간 채굴 요인의 중첩으로부터 진화되었습니다. 따라서 배수유동장을 분석하여 다양한 수리학적 특성을 구성할 때에는 모두 자연조건을 기초로 하고 채굴에 따른 영향을 충분히 고려해야 합니다.

(1) 비정상 흐름과 안정 흐름의 구별

일반적으로 배수 및 배수 시 광산 지역의 지하수는 대부분 불안정한 상태이지만, 배수 및 배수 시에는 지하수 함양량보다 적으면 안정된 상태가 될 수 있습니다.

채굴 초기 단계(채굴 단계)에서는 개발 터널이 계속 발전하고 변화하며 탈수 깔때기의 외부 경계가 계속 확장되고 광산에서 유입되는 물은 주로 대수층 저장고를 소비합니다. 이 단계에서 탈수 현장은 일반적으로 채굴의 후기 단계(추출 단계)에서 탈수 흐름은 주로 재충전이 부족한 채굴 지역의 외부 경계에서의 재충전 조건에 의해 제어됩니다. 조건에 따라 탈수 흐름장은 주로 대수층 저장소를 소비하며 이는 여전히 안정적인 흐름입니다. 충분한 재충전 조건이 있는 광산 지역이나 고정된 수두 재충전 경계가 있는 광산 지역에서는 광산수 유입(또는 탈수)이 발생합니다. 양)은 함양량에 의해 균형을 이루며 일반적으로 비교적 안정적인 흐름이 발생하므로 광산수 유입 예측은 우물 흐름 이론을 기본으로 안정화할 수 있습니다.

(2) 층류와 난류의 구별

건조 상태에서의 자연 이동 상태와 비교하면 광산 지역의 지하수는 여전히 넓은 지역에서 층류입니다. 건식 공학 근처에서 자주 발생하므로 Darcy의 법칙(직선 침투 법칙)은 결정론적 모델을 확립하는 기초로 남아 있습니다.

일반적으로 층류와 난류는 물 펌핑(방출) 테스트를 기반으로 한 단위 물 유입(qi) 방식으로 판단됩니다. 계산 공식은 다음과 같습니다.

압력 베어링 물.

특수 수문지질학

다이빙

특수 수문지질학

여기서: qi는 물 유입 단위입니다. Qi, Si는 안정된 물입니다. 펌핑 테스트의 각 프로세스에 대한 용량 및 안정적인 수위 감소 값.

매번 계산한 기가 일정하거나 상수에 가까우면 지하수는 층류이고, 그렇지 않으면 난류 또는 혼합류이다.

(3) 평면 흐름과 공간 흐름 구별

건류장에서 지하수의 이동 형태는 터널의 종류와 분포에 따라 제어되어 복잡한 흐름 패턴을 나타냅니다. , 이는 거시적으로는 두 가지 유형, 즉 완전한 터널로 흐르는 평면 흐름과 불완전한 터널로 흐르는 공간 흐름으로 일반화할 수 있습니다.

1) 완전한 터널까지의 평면 흐름. 터널을 완성하기 위해 흐르는 지하수의 평면유동은 평면복사유동과 단면유동의 두 가지 형태로 나눌 수 있으며, 터널의 종류에 따라 평면복사유동이 발생하며, 일반적으로 수직배수 시에는 평면복사유동이 발생하고, 수평배수 시에는 평면복사유동이 발생한다. 터널 배수 시에는 단면적 평면유동이 발생한다(실제로는 양쪽 끝에서 방사형 유동이 발생하는 경우가 많다). 평면 복사 흐름은 우물 흐름 공식으로 표현되고, 프로파일 평면 흐름은 단일 폭 흐름 공식으로 계산됩니다. 경사터널의 경우 수력유추법 연구에 따르면 터널의 경사도가 물유입량에 거의 영향을 미치지 않음이 입증되었으며, 터널의 경사도에 따라 터널의 종류에 따라 대략적인 계산이 가능하다. 수직 또는 수평 터널: 터널 경사도 α>45°이면 수직갱과 유사한 것으로 간주할 수 있으며, 복사정 흐름을 사용하여 계산할 수 있으며, 터널 경사도 α <45°인 경우 수평 터널과 유사한 것으로 간주할 수 있습니다. 프로파일 흐름의 단일 너비 흐름 공식을 사용하여 계산됩니다.

2) 불완전 터널로의 공간 흐름: 공간 흐름 영역은 불완전 터널 부근으로 제한되는 경우가 많으며 그 범위는 대수층 두께의 약 1.5~2.0배입니다. 완전한 수갱의 유입수량을 계산하기 위해서는 일반적으로 평면 분할법을 사용하고, 불완전한 터널의 유입수량을 근사화하기 위해서는 프로파일 분할법이나 경험식을 사용합니다.

(4) 지하수와 가압수 구별

탐사 시 지하수와 가압수는 쉽게 판별할 수 있지만, 광산을 감압하고 배수할 때 문제가 자주 발생합니다. 압력은 비가압수 또는 가압-비가압수로 변환됩니다. 또한 경사진 암석층이 있는 일부 광산 지역에서는 한쪽은 가압된 물의 원래 상태를 유지하고 다른 쪽은 가압되지 않은 물로 변환되는 일이 발생할 수 있습니다. 그림 13-1) 또는 압력이 가해지지 않은 물. 일반화할 때 위에서 언급한 변화하는 특성을 충분히 고려해야 하며 등가 근사 처리는 거시적 관점에서 수행되어야 합니다.

그림 13-1 광산 지역의 배수 깔때기 개략도

2. 경계 유형 결정

(1) 측면 경계 조건의 일반화

①경계 유형(또는 속성)의 구분: 경계 유형(또는 속성)은 물 격리 경계와 물 공급 경계의 두 가지 범주로 나눌 수 있습니다. 방수 경계는 대수층과 약투수층, 지하수 또는 물 차단 단층 사이의 경계선을 말하며, 물 공급 경계는 대수층과 같이 재충전 용량이 무제한인 일정한 수두 재충전 경계 윤곽을 의미합니다. (강한 수력학적 연결을 가짐) 지표수의 접촉 경계 또한 일부 강한 대수층은 약한 대수층의 일정한 수두 급수 경계가 될 수도 있습니다. 따라서 경계는 광산의 실제 상황에 따라 물격리경계 또는 물공급경계로 일반화되어야 한다. 일반화할 때 전반적인 일반화 효과를 기반으로 하고 마이닝 요소를 충분히 고려하고 경계 기하학에 주의를 기울이고 "등가" 원칙을 잘 활용해야 합니다.

그림 13-2 광산 지역의 복사 흐름 계산에 대한 개략도

1—1번과 2집수 지점의 수위 범위 2—유선을 나누는 블록 세그먼트; ; 3—관찰 구멍

②경계 형태의 단순화: 분석 계산 모델에서는 불규칙한 경계 형태를 무한 직선, 직교 경계, 경사 경계 및 평행 경계와 같은 이상적인 기하학적 패턴 경계로 단순화해야 합니다. . 기다리다. 또한 경우에 따라 경계를 단순화하기 위해 분할 방법을 사용할 수도 있습니다. 즉, 건류장의 특성에 따라 유동 표면을 따라 조건이 다른 여러 개의 부채꼴 분포 영역으로 분할되고 동일한 수위(압력) 표면(그림 13-2)이므로 분할 방법을 사용하여 물의 양을 계산할 수 있습니다.

(2) 수직 경계 유형의 일반화

고갈 대수층의 지붕과 바닥의 특성과 해석 방법의 계산 요구 사항에 따라 일반적으로 수직 경계 유형으로 구분됩니다. 물 격리 경계와 수직 교차 흐름(교차층) 재충전 경계에는 두 가지 유형이 있습니다. 전자 대수층의 상부 및 하부 판은 불투수층이며, 후자 대수층의 상부 및 하부 판은 약투수층이며, 수직으로 인접한 대수층은 교차 흐름을 통해 약투수층을 통해 건조층을 재충전합니다( 교차층) 재충전 흐름(교차층) 재충전 경계는 고정 수두와 가변 수두의 두 가지 유형으로 나눌 수 있습니다. 현재 분석 방법은 전자의 문제를 해결할 수 있습니다.

(3) 내부 경계의 일반화

내부 경계는 채굴 시스템의 기하학적 형태이며, 원과 같은 이상적인 기하학적 패턴으로도 일반화되어야 한다. , 수평터널 등 일반적으로 , 대유정법을 계산에 사용하기 위해 터널 시스템을 일반화하거나 원형으로 변환하는 경우가 많다.

(2) 수학적 모델 수립(예측계산식)

수학적 모델은 수리지질학적 개념모델을 수학적으로 기술한 것으로 수학적 모델의 수립은 수문지질학적 개념에 기초한다. 개념모델과 배수정 및 차선의 특성을 파악하고, 계산방법을 정확하게 결정하고, 계산식을 선택하거나 설정합니다.

1. 큰 우물 방법

터널 시스템의 배수는 일반적으로 터널을 중심으로 하는 방사형 흐름을 형성합니다. , 반경 ro를 갖는 동등한 이상적인 "큰 우물", 전체 터널 시스템의 물 유입은 "큰 우물"의 물 유입과 동일하므로 우물 흐름 공식을 사용하여 광산의 물 유입을 예측할 수 있습니다. 무한 대수층에 대해서는 식 (13-21)과 식 (13-22)에 따라 광산수 유입량을 계산하고 예측할 수 있다. 다양한 특정 경계 조건에 대해 분석 공식은 다음과 같은 일반 형식으로 작성될 수 있습니다.

정상 흐름:

특수 수문지질학

불안정 흐름:< /p >

특수 수문지질학

수식에서: R?와 Rr은 안정 유동과 비정상 유동의 경계 유형 조건 계수이며, 다양한 이상화된 경계 유형 조건 계수는 표 13-1에 나와 있습니다. 이전과 같은 의미를 갖습니다. 표에 나열된 표현식은 분석 솔루션의 원형이 아니라 단순화된 근사 표현식입니다.

표 13-1 이상화된 경계 유형 조건 계수

계속 표

마찬가지로, 오버플로 재충전을 위해 수직 방향으로 일정한 수두가 있는 경우 오버플로가 발생할 수 있습니다. 광산으로 유입되는 물의 양을 계산하기 위한 우물 흐름 공식은 "Groundwater Dynamics"와 같은 책을 참조하시기 바랍니다. 여기서는 반복하지 않겠습니다.

2. 분할 방법

카민스키 복사 흐름 방법이라고도 합니다. 이 방법은 건유동장의 경계조건과 대수층의 이질적인 특성을 바탕으로 유동면과 등유수면을 따라 조건이 다른 여러 개의 부채꼴 분포지역으로 구분한다(그림 13-2). 형상분포면적 해당 지역의 지하수 흐름은 방사상이며, 유동면을 따라 분할된 부채꼴 영역의 경계는 물 차단 경계이고, 등수면을 따라 분할된 부채꼴 영역의 새로운 경계는 물 차단 경계입니다. 동등한 머리 경계. 현재 카민스키 평면 복사 흐름 공식은 각 섹터 구역의 물 유입 Qi를 계산하는 데 일반적으로 사용되며, 각 구역의 물 유입량(Q)의 합은 다음과 같이 계산됩니다.

제한 수역

특수 수문지질학

다이빙

특수 수문지질학

특수 수문지질학

여기: b1, b2는 방사형 유동의 상류 및 하류 구간의 폭(m), h1 및 h2는 구간 b1 및 b2의 방수 바닥판의 수두 높이(m)입니다. L은 구간 b1과 b2 사이의 거리입니다. 중). 각 구간의 하류 구간은 터널 배수 깔때기에 바로 인접한 등수선 부분을 기준으로 하고, 상류 구간은 터널의 급수경계에서 떨어진 등수두(위치)선 부분을 기준으로 합니다. .

3. 수평 터널 해석 방법

앞서 언급한 바와 같이 수평 터널 건조 및 배수는 일반적으로 단면 평면 흐름으로 단일 폭 흐름 공식으로 계산할 수 있습니다. "Groundwater Dynamics"” 및 관련 서적을 참조하십시오.

(3) 매개변수 결정

계산 매개변수에는 대수층 두께(M), 투수계수(K), 급수량(μ), 대형 우물 반경( ro), 영향 반경(R) 및 영향 구역 폭(L), 최대 수위 감소 값(Smax) 등 광산 내 유입수량을 예측하는 분석방법의 정확도는 앞서 언급한 조건뿐만 아니라, 사용된 다양한 매개변수의 정확도와 대표성, 광산 조건의 변화를 올바르게 추정했는지 여부에 따라 결정됩니다. 즉, 계산변수는 광산수 유입예측의 정확성에 직접적인 영향을 미치기 때문에 분석계산모델의 특성과 광산지역의 수문지질학적 조건, 향후 채굴계획 등을 고려하여 다양한 매개변수를 합리적으로 결정해야 한다. 여기서는 광산 유입수 특성을 바탕으로 특정 매개변수를 결정하는 방법만 소개합니다.

1. 대수층의 두께(M 또는 H)

분석 방법에서는 대수층의 두께를 동일한 두께로 간주하므로 전체 평균값이 필요합니다. 계산할 때에는 암석학 및 지질학적 특성을 기준으로 여러 계산 영역으로 나누어야 하며, 각 영역의 면적은 Fi로 하며, 각 영역의 평균 두께 Mi(또는 Hi)를 분할하여 계산한다. 광산 지역의 전체 평균 두께 Mcp(또는 Hcp)를 구하는 데 사용되는 방법:

특수 수문지질학

암층의 투과성 계수(K)

투과계수(K)는 분석계산 매개변수에서 매우 중요합니다. K 값을 정확하게 구하기 위해서는 먼저 대수층이 균질한지, 이질적인지 구별해야 합니다. 분석 방법은 주로 균질한 대수층에 적합합니다. 우리나라의 광물 매장지는 대부분 물이 채워진 이질적인 균열과 카르스트 대수층에서 생산됩니다. 일반적으로 이질성의 정도가 너무 크지 않은 경우에는 평균투과계수(Kcp)를 구하여 K 값을 결정할 수 있으며, 차이가 큰 경우에는 평균 일반화 또는 분할 계산을 수행해야 합니다. 다음과 같은 방법이 일반적으로 사용됩니다.

1) 가중평균법 : 대수층의 특성에 따라 두께평균법, 면적평균법, 방향평균법으로 구분할 수 있으며 그 중 앞의 두 가지 방법이 더 많이 사용되며, 계산 공식은 다음과 같습니다:

두께 평균 방법:

특수 수문지질학

면적 평균 방법:

특수 수문지질학

여기서: Kcp는 평균 투과성 계수(m/d)이고, Mi(Hi)는 제한된 수(또는 지하) 대수층의 각 수직 세그먼트의 두께(m)입니다. m/ d) Fi는 대수층의 각 구획 면적(m2)입니다.

2) 유동장 분석 방법: 이 방법은 물 펌핑(방출) 테스트 데이터를 기반으로 그려지거나 유동장의 전반적인 특성을 기반으로 여러 구체적인 특징으로 나누어 상대적으로 신뢰할 수 있는 등수소 지도가 필요합니다. 특성이 다른 단면에 대해서는 전자를 폐곡선법, 후자를 분할법이라고 합니다.

폐쇄형 등선법: 물 펌핑(방출) 테스트 데이터 또는 Darcy의 법칙에 따라 등수위를 그립니다.

특수 수리지질학

공식에서: Q는 광산으로부터의 물 유입이고, L1과 L2는 두 개의 닫힌 윤곽선의 길이이고, Δr은 두 개의 닫힌 윤곽선의 수위 차이입니다. 대수층의 평균 두께(m), 기타 기호의 의미는 이전과 동일합니다.

분할법 : 유동장의 특성에 따라 특성이 다른 여러 구간으로 나누어 분할법을 이용하여 계산한다. 식(13-27)과 식(13-28)에 따르면 분할법에 의한 K 값을 계산하는 공식을 유도할 수 있다. 예를 들어 다이빙인 경우 계산 공식은 다음과 같습니다.

특수 수문지질학

공식에서 기호의 의미는 이전과 동일합니다.

3. '큰 샤프트'(ro)의 기준 반경

갱의 모양은 매우 불규칙하여 복잡한 내부 경계를 형성하지만 터널 시스템이 배수되면 둘레는 점진적으로 통일된 낙하 깔때기를 형성할 수 있으므로 이론적으로 복잡한 모양의 터널 시스템은 동등한 이상적인 대형 우물로 간주될 수 있으며 광산의 유입량은 우물 흐름 공식에 따라 계산됩니다. 글쎄 방법. "큰 우물"의 기준 반경 ro는 일반적인 상황에서 다음 공식으로 계산할 수 있습니다:

특수 수문지질학

공식에서: F는 분포 범위로 구분되는 영역입니다. 터널 시스템의 (m2).

4. 건식 터널(시스템) 배수 시 영향 반경(Ro) 또는 영향 구역 폭(Ltp)

'큰 우물' 방법을 사용하는 경우 광산에 유입되는 물의 양을 예측합니다. 착지 깔때기의 영향 반경 Ro는 "큰 우물"의 중심에서 계산되어야 하며 "큰 우물"의 기준 반경(ro)에 다음을 더한 값과 같습니다. 배수 영향 반경(R), 즉:

Ro=ro R (13- 36)

건조 깔대기의 모양이 불규칙한 경우가 많으므로 사용하는 것이 더 합리적입니다. Ro를 나타내는 "영향력 반경"(Rcp). 길고 좁은 수평터널의 유입수량을 계산할 때, 배수 및 배수의 영향범위는 '영향대 폭(Lcp)'으로 표현됩니다.

기준 영향 반경(Rcp)과 영향 대역 폭(Lcp)을 결정하는 방법에는 여러 가지가 있습니다. 일반적으로 사용되는 방법은 다음과 같습니다.

1) 경험적 및 반경험적 공식:

다이빙 Kusarkin 공식:

특수 수리지질학

제한 수역 West Hatt 공식:

특수 수문지질학

2) Serovatko의 공식: 복잡한 터널 시스템의 경우 영향 반경은 터널 시스템 가장자리와 자연 수문지질 경계선(예: 재충전 경계선) 사이의 거리의 가중 평균을 기반으로 계산되어야 합니다. ):

특수 수문지질학

공식에서: ro는 터널 시스템의 "큰 유정"의 기준 반경(m)이고, bcp는 터널 외곽선 사이의 평균 거리입니다. 다양한 유형의 수문지질학 경계선(m) L은 다양한 유형의 수문지질학 경계선의 너비(m)입니다(그림 13-3).

3) 외삽 방법: 다방울 펌핑 테스트 데이터를 기반으로 R과 S 또는 R과 Q 사이의 선형 관계를 결정하고 특정 탈수 수준 또는 특정 탈수의 해당 탈수 영향을 외삽합니다. 금액. 예:

그림 13-3 터널 시스템 공급 경계의 개략도

그림 13-4 터널 배수 수위 지도

특수 수문지질학

또는

특수 수문지질학

공식에서: α는 비례 계수입니다.

4) 그리기 방법: 여러 낙하 스트로크의 펌핑 테스트 데이터를 기반으로 R=f(S) 및 R=f(Q) 곡선을 그리고 이를 외삽하여 R을 계산합니다.

5. 최대 수위 하강(Smax)

안전한 채굴을 위해 직접 물이 채워진 층의 경우 광산 배수에 대한 최대 수위 하강(Smax)은 다음과 같아야 합니다. 지하수위 터널 바닥까지 낮추거나(그림 13-4a), Smax=H(대수층 수위 또는 두께)를 취합니다(그림 13-4b). 이론적으로 말하면, 최대 수위하강의 최대 탈수량을 계산하기 위해 분석적 방법을 사용하는 것은 적절하지 않으며 종종 "왜곡"됩니다. 실제 문제를 해결할 때 일반적으로 광물 매장지의 탈수를 위한 최대 수위 강하를 계산하는 것이 필요합니다.

(4) 광산수 유입량 예측

매개변수가 결정된 후 수학적 모델(해석식)과 관련 생산 요구 사항.

1. 안정적인 흐름 분석 방법

계산 및 예측 내용은 주로 두 가지 측면을 포함합니다.

1) 알려진 지점에서의 최대 수위 감소(Smax) 광산 수준), 광산의 총 물 유입량(Qtotal)을 예측합니다.

2) 주어진 배수능력(Q)을 전제로 해당 면적의 수위강하(또는 수압강하)값(S)을 계산한다.

계산 예 탄광의 탄층은 페름기 사암으로 둘러싸인 대수층 아래에 ​​묻혀 있습니다. 석탄층은 투과성이 있지만 수분이 많지 않은 단층에 의해 절단됩니다. 터널 시스템의 면적은 109.4×104m2이며, 그 윤곽은 불규칙한 원형이다. 사암 대수층의 평균 수두는 100m, 대수층 두께는 30m, 투수계수는 0.2m/d이다. 터널 시스템은 방수 셰일 바닥(그림 13-5)에 배치되었으며, 형성 경사각은 13°입니다. 광산에 유입되는 물의 양을 예측해 보세요.

해결책 알려진 조건에 따르면 단층의 투수성은 물격리경계가 아니라는 것을 의미하고, 물풍부도가 부족하다는 것은 물 공급경계가 아니라는 것을 의미하므로, 무한한 경계로 간주됩니다. 지층의 경사는 13°로 수평 대수층으로 간주될 수 있습니다. 대수층 두께와 투과성 계수가 일반화되었습니다. 그림 13-5에서 볼 수 있듯이 광산이 배수되면 수위가 방수지붕 아래로 떨어지게 되는데, 이는 압력이 가압되지 않은 물이므로 광산수를 계산하기 위해서는 무한경계압력-비압력 공식을 사용해야 한다. 유입.

그림 13-5 탄광의 도식적 단면

특수 수리지질학

r0과 R0은 다음 공식으로 계산됩니다.

특수 수문지질학

특수 수문지질학

h=0 (Smax=H)

관련 값을 공식에 ​​넣어서 구하세요.

Q=3666m3/d

실제 채굴 과정에서 광산 내 유입수량은 3600m3/d로 예상 결과에 근접했습니다.

2. 예측량 계산을 위한 비정상류 해석방법

1) 수위강하(S)와 탈수시간(t)을 고려하여 예상탈수량(Q)은 다음과 같다. 일정 기간(보통 두 우기 사이) 내에 특정 채굴 수준의 탈수 작업을 완료하고 합리적인 탈수량을 선택하거나, 탈수량이 채굴 깊이에 도달하면 수량을 예측하는 것이 필요합니다. 시간에 따른 광산 피트의 유입량(계절) 패턴을 변화시켜 장마철 최대 유입량과 발생 시간을 구합니다.

2) 배수량(Q)과 수위 저하(S)를 고려하여 탈수 시간(t)을 예측합니다. 즉, 배수 용량을 기준으로 특정 배수에 도달하는 데 필요한 시간을 계산합니다. 수준 또는 추가로 깔때기가 중요한 외부 경계로 확장되는 시기를 예측합니다.

3) 물 용량(Q)과 배수 시간(t)을 고려하여 배수 깔때기의 각 지점에서 수위 강하(S)를 예측합니다. 즉, 배수 용량에 따라 광산 지역의 지하수 배수 건식 감압 깔때기의 형성 및 확장, 깔때기 범위 내의 각 지점에서 시간에 따른 수위 변화 계산, 채굴 순서, 속도 및 기타 채굴 조치를 계획합니다.

4) 일정한 낙하 및 깊은 가변 흐름 계산. 탐사단계에서 광산수 유입량 예측을 위한 비정상유량 계산은 주로 탈수량 선정 및 최대 유입량 예측을 기반으로 한다. 탈수시간과 배수능력이 제한되지 않는 경우 사전배수는 일반적으로 탈수유량에 따른 수위저하 및 탈수시간을 먼저 계산한 후 탈수유량에 따른 S-t 곡선을 그린다(그림 13-6). )과 동일한 설계 드롭의 Q-t 곡선(그림 13-7)을 제공하여 설계 부서에서 합리적인 솔루션을 선택할 때 이를 참조할 수 있도록 합니다.

그림 13-6 다양한 배수 흐름에 대한 S-t 곡선

그림 13-7 수위 감소 Sk에 대한 Q-t 곡선