A new finite difference model is developed for solute transport in a fractured medium that can consider advection, adsorption, first-order decay, and scale-dependent dispersivity of individual fractures. In the model, the dispersivity of individual fractures is employed as a variable increasing with travel distance from a source. The model is verified using an analytical solution for a single fracture. A solution from the new model is independent of the outlet boundary condition of fractures, and has little numerical dispersion error.
Khakimov, Elyorbek;Chung, Sang Yong;Senapathi, Venkatramanan;Elzain, Hussam Eldin;Son, JooHyeong
한국지하수토양환경학회지:지하수토양환경
/
제22권3호
/
pp.27-41
/
2017
This study dealt with the characteristics and the interrelations of hydrogeological parameters such as hydraulic conductivity, dispersivity and effective porosity of unconsolidated sediments for providing the basic data necessary for the planning of the management and preservation of groundwater quality in the Nakdong River Delta of Busan City, Korea. Groundwater quality in this area has been deteriorated due to seawater intrusion, agricultural fertilizer and pesticide, industrial wastewater, and contaminated river water. The physical properties (grain size distribution, sediment type, sorting) and aquifer parameters (hydraulic conductivity, effective porosity, longitudinal dispersivity) were determined from grain size analysis, laboratory permeability test and column tracer test. Among 36 samples, there were 18 Sand (S), 7 Gravelly Sand (gS), 5 Silty Sand (zS), 5 Muddy Sand (mS), and 1 Sandy Silt (sZ). Hydraulic conductivity was determined through a falling head test, and ranged from $9.2{\times}10^{-5}$ to $2.9{\times}10^{-2}cm/sec$ (0.08 to 25.6 m/day). From breakthrough curves, dispersivity was calculated to be 0.35~3.92 cm. Also, effective porosity and average linear velocity were obtained through the column tracer test, and their values were 0.04~0.46 and 1.06E-04~6.49E-02 cm/sec, respectively. Statistical methods were used to understand the interrelations among aquifer parameters of hydraulic conductivity, effective porosity and dispersivity. The relation between dispersivity and hydraulic conductivity or effective porosity considered the sample length, because dispersivity was affected by experimental scale. The relations between dispersivity and hydraulic conductivity or effective porosity were all in inverse proportion for all long and short samples. The reason was because dispersivity was in inverse proportion to the groundwater velocity in case of steady hydrodynamic dispersion coefficient, and groundwater velocity was in proportion to the hydraulic conductivity or effective porosity. This study also elucidated that longitudinal dispersivity was dependent on the scale of column tracer test, and all hydrogeological parameters were low to high values due to the sand quantity of sediments. It is expected that the hydrogeological parameter data of sediments will be very useful for the planning of groundwater management and preservation in the Nakdong River Delta of Busan City, Korea.
In this paper, solute transport in heterogeneous aquifers using a modified Fokker-Planck equation (MFPE) is investigated. This newly developed mathematical model is characterised with a time-, scale-dependent dispersivity. A two-dimensional finite volume quadrilateral mesh method (FVQMM) based on a quadrilateral background interpolation mesh is developed for analysing the model. The FVQMM transforms the coupled non-linear partial differential equations into a system of differential equations, which is solved using backward differentiation formulae of order one through five in order to advance the solution in time. Three examples are presented to demonstrate the model verification and utility. Henry's classic benchmark problem is used to show that the MFPE captures significant features of transport phenomena in heterogeneous porous media including enhanced transport of salt in the upper layer due to its parameters that represent the dependence of transport processes on scale and time. The time and scale effects are investigated. Numerical results are compared with published results on the some problems.
투수성이 높고 대수층의 저유량이 풍부한 강변여과수 개발 예정지역의 충적층에서 시험규모가 2m (IW-1공와 5m (IW-2공)인 경우의 수렴흐름 추적자시험이 수행되었다. 수렴흐름 추적자시험의 양수율은 $2,500m^3/day$ 이었으며, 주입정 (IW-1 및 IW-2공)에 염소이온 5kg을 순간 주입하였다. 염소이온의 농도이력곡선을 작성하여 초기도달시간과 최고농도의 차이를 분석하였으며, 누적질량회수곡선을 통해 양수 후 경과시간에 따른 염소이온의 질량회수율을 분석하였다. 그리고, 염소이온농도 대 누적질량회수율의 이력그래프를 작성하여 누적질량회수율에 따른 염소이온농도의 증가와 감소 변화를 분석하였다. 또한, 염소이온농도의 증가/감소 구간에 대한 선형회귀분석을 수행하여 농도증가율과 감소율의 변화를 파악하였다. 양수정에서 관측된 경과시간별 염소이온농도 자료를 CATTl 코드의 quot;Converging Radial Flow With Instantaneous Injectionquot; 해석법에 적용하여 종분산지수를 추정하였다. 추정된 종분산지수는 양수정과 IW-1공 사이의 충적층에서는 0.4152 m, 양수정과 IW-2공 사이의 충적층에서는 3.2665 m이었다. 양수정에서 이격거리가 멀수록 종분산지수는 증가하였으며, 이격거리에 대한 종분산지수의 비는 각각 0.21과 0.65 정도이었다.
본 연구에서는 복합토양층(조립질 모래, 세립질 모래, 실트 점토의 혼합토)에서 수행된 현장추적자시험과 3개 토양층에서 채취된 시료를 이용한 실내주상실험의 자료를 이용하여, 복합토양층의 수리분산특성을 분석하였다. 토양층별 투수성과 유동 분석에 의해, 유효공극률이 낮은 실트 점토의 혼합토와 세립질 모래에서 평균선형유속이 높고, 유효공극률이 큰 조립질 모래에서는 수리전도도가 높은 것을 알 수 있었다. 평균토양입경에 따른 수리전도도 함수는 Y=$3.49{\times}10^{-8}e^{15320x}$로 추정되었으며 결정계수는 0.90이었다. 평균토양입경에 따른 평균선형유속 함수는 Y=$1.88{\times}10^{-7}e^{11459x}$로 추정되었으며, 결정 계수는 0.81이었다. 그리고 평균토양입경에 따른 종분산지수 함수는 Y = 0.00256$e^{5971x}$이었으며, 결정계수는 0.98 정도로 나타났다. 평균선형유속과 종분산지수의 선형회귀분석 결과, 함수는 Y = 21.7527x+0.0063로 추정되었으며 결정계수는 0.9979로서 매우 높게 나타났다. 본 연구에서 산정된 현장규모/실내쥬모의 종분산지수비는 54.09로서 규모종속효과를 나타내었다. 현장추적자시험을 실시하여 산정한 종분산지수(1.39m)가 Xu와 Eckstein(1995)의 방법에 의해 산정된 종분산지수(0.186m)에 비해 7.47배 정도 크게 나타났다. 이는 시험 대수층 내 중 조립질 모래층에서의 우세한 확산 흐름에 의한 것이다.
여러 개의 층으로 구성된 토양층에서 연속주입 추적자시험을 수행하여, 불포화대와 포화대 구간에서의 수리분산특성 차이를 추적자의 농도이력곡선, 시간에 따른 농도변화 및 농도변화율, 그리고 거리에 따른 농도비 분석을 통해 비교 연구하였다. 연속주입 추적자시험에 의하면, 약 160시간이 경과된 후에 불포화대와 포화대에서 Rhodamine WT 최대농도의 차이는 약 $13{\sim}15$배 정도에 달하였고, 시험시간 대 추적자 농도증가율의 차이는 약 10배 정도로 나타났다. 또한 시간에 따른 추적자 농도이력곡선의 변화와 농도변화율이 불포화대에 비하여 포화대에서 크게 나타났다. 주입공에서의 이격거리에 따른 추적자의 농도비는 포화대 구간에서 더 빠르게 선형적으로 감소하였으며, 그 이후에 농도비가 2배 이하로 완만해 지는데 걸린 시간은 포화대 구간이 더 길었다. 이러한 차이들은 여러 개의 층서로 구성된 토양층의 포화대에서는 지하수가 존재하며, 또한 매질의 불균질성이 크고 투수성이 다양하기 때문에, 비교적 균질한 층서를 이루는 불포화대에 비하여 추적자 용액의 농도는 낮고, 추적자의 확산이 느리게 진행되어진데서 기인하였다. 그리고 포화대 구간에서의 유효공극율은 $10.19{\sim}10.50%$, 종분산지수는 $0.80{\sim}l.98m$, 횡분산지수는 $0.02{\sim}0.04m$의 범위로 산정되었으며, 실내주상시험의 종분산지수와 비교할 때 12배 이상의 규모종속효과가 나타났다.
대수층의 저유량이 풍부한 강변여과수 개발 예정지역의 충적층(지표면하 $25{\sim}35\;m$ 구간)에서 수리전도도와 종분산지수의 규모종속효과를 규명하기 위해 양수시험과 수렴흐름 추적자시험이 수행 되었다. 양수시험과 추적자시험의 규모는 2 m 와 5 m 이었으며 양수시험은 5개 공, 추적자시험은 3개 공을 이용하여 수행되었다. 양수시험은 일정한 양수율($2,500\;m^3/day$)로 수행되었으며, 양수 시작 후 경과시간에 따른 수위변화 자료를 AQTESOLV 3.5 프로그램에 입력하여 해석하였다. 시험 대수층의 수리전도도는 양수정에서 $1.745{\times}10^{-3}\;m/sec$, 양수정에서 이격거리가 2 m 구간에서는 $2.161{\times}10^{-3}\;m/sec$와 $2.270{\times}10^{-3}\;m/sec$ 이었으며, 이격거리가 5 m 구간에서는 $2.452{\times}10^{-3}\;m/sec$와 $2.591{\time}10^{-3}m/sec$로 산정되었다. 그리고, 양수정에서 회복시험 시 Theis(Recovery) 방법에 의해 해석된 수리전도도는 $1.603{\times}10^{-3}\;m/sec$이었다. 양수정에서 관측정의 이격거리(d)에 따른 수리전도도(K) 증가함수는 log K=0.0693 log d-2.671와 log K=0.0817 log d-2.655로 추정되었으며, 결정 계수는 각각 0.965와 0.979로서 매우 높게 나타났다. 따라서 양수정에서의 이격거리가 멀수록 수리전도도가 증가하는 규모종속을 확인하였으며, 또한 시험대수층의 수리전도도가 방사상으로 유사하게 분포하고 있음을 알 수 있었다. 수렴흐름 추적자시험의 양수율은 $2,500\;m^3/day$ 이었으며, 2개의 주입정에 염소이온 5 kg을 순간 주입하였다. 염소이온의 농도이력곡선을 작성하여 초기도달시간과 최고농도의 차이를 분석하였으며, 누적질량회수곡선을 통해 양수 후 경과시간에 따른 염소이온의 질량회수율을 분석하였다. 그리고, 염소이온농도 대 누적질량회수율의 이력그래프를 작성하여 누적질량회수율에 따른 염소이온농도의 증가와 감소 변화를 분석하였다. 또한, 염소이온농도의 증가/감소 구간에 대한 선형회귀분석을 수행하여 농도 증가율과 감소율의 변화를 파악하였다. 양수정에서 관측된 경과시간별 염소이온농도 자료를 CATTI 코드의 "Converging Radial Flow With Instantaneous Injection" 해석법에 적용하여 종분산지수를 추정하였다. 양수정에서 이격거리가 2 m인 경우의 종분산지수는 0.4152 m, 이격거리가 5 m인 경우의 종분산지수는 3.2665 m이었다. 따라서 양수정에서 이격거리가 멀수록 종분산지수가 증가하는 규모종속효과를 확인하였으며, 또한 이격거리에 대한 종분산지수의 비는 각각 0.21과 0.65 정도로서 증가하였다.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.