• 터널과지하공간
• /
• 제13권1호
• /
• pp.52-63
• /
• 2003

단층, 균열, 절리 등의 암반 내에 존재하는 수많은 불연속면은 암반의 역학적\ulcorner수리적 특성을 좌우하는 중요한 요소이다. 암반 내 지하수 유동에 큰 영향을 미치는 요소는 불연속면의 발생빈도와 기하학적 특성 그리고 불연속면 상호간의 연결성이라고 할 수 있다. 이 연구에서는 불연속면의 특성인 발생빈도, 크기, 방향, 간극의 크기 등의 분포함수를 가정하여 3차원 불연속면 연결망 내에서 지하수 유동 해석을 실시하는 프로그램을 작성하였다. 이 프로그램은 3차원 상에 불연속면을 발생시키고 불연속면간 연결성을 분석하여 수리해석을 실시한다. 이 프로그램을 통해 수리해석을 실시한 결과, 컴퓨터 메모리의 한계로 인해 대상지역의 대표요소체적을 정확하게 결정할 수는 없었지만, 대략 25$\times$25$\times$25 ㎥ 이상에서 결정될 것으로 추정할 수 있었다. 지하수 유동해석에 영향을 미치는 간극의 범위를 계산한 결과, 불연속면 평균 간극의 30% 이하의 간극을 갖는 불연속면은 지하수 유동에 미치는 영향이 미미한 것으로 나타났다. 또, 경계효과를 고려한 경우와 그렇지 않은 경우의 등가수리전도도의 차이는 거의 없었으며 이는 대상지역에서 간극이 큰 일부의 불연속면을 통한 유동이 전체적인 지하수 유동에 크게 영향을 주기 때문인 것으로 판단되었다. 입력자료 중에서 암반의 등가수리전도도에 영향을 미치는 요소의 중요도는 불연속면의 길이, 간극, 방향의 순서로 나타났다. 대상 암반에 단층면이 존재할 경우, 등가수리전도도는 단층면에 평행한 방향의 요소는 증가하며 이에 수직인 방향의 요소는 약간 증가하다 수렴하는 경향을 보였다.openyl methyl disulfide와 (E)-propenyl methyl disulfide 또한 동결건조 양파에서 증가되어 짐을 확인하였다. 올그루밀이 가장 많았다. Niacin함량은 탑동밀이 2.81 mg%로 가장 높은 함량을 보였고, 다음으로 알찬밀, 올그루밀 순이었다. 지방산 조성은 보리와 밀에서 Cl8:2>Cl6:0>Cl8:1 순으로 보리는 전체의 90%, 밀은 92%를 차지하였다. 단일 불포화지방산은 보리가 11∼13%, 밀이 21∼27%이며, 다중불포화지방산은 보리가 57∼59%, 밀은 36∼50%로 보리가 더 많은 것으로 나타났다. PUFA/SFA 비율은 보리의 경우 2.1로 품종별로 지방산 조성에 차이를 나타내지 않은 반면 밀은 1.0∼1.9 범위로 품종별로 지방산 조성에 차이를 나타내었다. 보리와 밀의 아미노산 함량은 glutamic acid를 가장 많이 함유하는 것으로 분석되었다. 보리 품종별 필수아미노산 함량을 살펴보면 Lysine, valine, tryptophan 함량은 두산8호에서, phenylalanine 함량은 서둔찰보리에 많이 함유되어 있는 것으로 분석되었다. 밀 품종별 Iysine, isoleucine 함량은 탑동밀에서 다소 낮았다.) 보였다. 체질량지수와 비만지수는 각각 HDL-콜레스테롤과는 부의 (각 P<0.05), 적혈구수와는 정의(각 p<0.05) 상관관계를 보였다. 허리엉덩이둘레비는 혈청 GPT, glucose, MCV와 각각 유의한 정의 상관관계를 보였다(각 p<0.05). 이상의 연구결과를 종합할 때 남녀 비만 중학생 모두 총 열량 섭취량 중 지질로 인한 열량 섭취비율이 높았고 비만도가 증가할수록 콜레스테롤의

• 자원환경지질
• /
• 제38권2호
• /
• pp.177-185
• /
• 2005

이 연구는 양수시험 해석해(Theis, 1935; Cooper-Jacob, 1946; Papadopulos-Cooper, 1967; Hantush, 1962a,b; Moench, 1985; Hantush-Jacob, 1955) 및 일반 방사상 유동 모델을 이용하여 우리나라의 균열암반 대수층(화강암, 화산암, 변성암, 백악기퇴적암, 제3기 퇴적암에 굴착된 100개 조사공)에서 수행되어진 양수시험으로부터 얻은 122개의 양수시험자료(수위강하 자료)를 분석하여 종합한 것이다. AQTESOLV 전산프로그램을 이용한 양수시험자료 분석에 의 하면, 122개 자료중 86개($71\%$)의 자료들이 이 연구에 사용된 해석해와 일치하며, 앙수시험자료 해석해 중에 누수(leaky) 및 경계조건(boundary condition)을 고려한 해석해들이 53개($43\%$)로 가장 많이 나타났다. 따라서, 양수시험자료의 해석은 균열암반 대수층의 수리지질학적 특성에 적합한 개념모델의 설정이 중요하다. 일반 방사상 유동(GRF)모델을 적용해보면, 122개의 자료 중 77개($63\%$)의 자료들이 Barker(1988)의 표준곡선에 의한 차원(1.1차원-2.9차원)을 보여준다. 이중 $44.2\%$에 해당하는 39개 자료가 1.1차원과 1.9차원 사이의 실수 유동차원을 보여주는 반면에 26개($6.5\%$)만이 Theis 이론에 맞는 2차원의 방사상 흐름을 보여주며, 38개($49.3\%$)는 2.1차원에서 2.9차원에 속한다. 따라서 우리나라 균열암반 대수층에서 지하수 유동은 대부분 실수차원의 유동을 보여주는 것으로 평가된다.

• 한국토양환경학회지
• /
• 제1권1호
• /
• pp.89-102
• /
• 1996

지하의 유체 유동 및 물질 변환을 해석하기 위하여 다중다상이론을 이용한 통합 모형을 개발하였다 종합적 지배식은 4개의 상내의 화합물들의 물질 및 힘평형 관계를 고려하여 유도되었다. 복합한 이동 및 변환 현상을 설명하고, 공간적 차원을 변동적으로 나타내기 위하여 관계된 모든 변수 및 식들을 함축적이면서 조직적으로 표현하였다. 도출된 비선형시스템은 다차원 유한요소프로_I램으로서 해를 구하였다. 본 개발된 프로그램은 역동적으로 메모리 용량을 조절하여 일이삼차원 문제를 PC부터 SP2슈퍼컴퓨터까지 여러 종류의 기종에서 해석할 수 있다. 계산시간과 저장용량을 줄이기 위하여 시스템식을 분리시키고, 슈퍼컴의 벡터 및 병렬처리를 이용하여 띠행렬의 해를 구하였다. 유속이 우세한 경우의 수치해석상의 불안정한 문제를 해결하기 위하여 상류가중, 질량묶음, 요소별 파라미터 평가법 등을 적용하였다. 일차원 이동문제에 대하여 유한요소법과 유한차분법의 수치해의 안정성 조건을 검토하였다. 구체적인 지하수 유동 및 오염문제에 대한 모델링 예는 본 논문집의 연계 논문에 수록하였다.

• 자원환경지질
• /
• 제30권5호
• /
• pp.433-442
• /
• 1997

균열암반 이중대수층 (1대수층과 2대수층)내 지하수흐름을 모식화한 프락탈모델이 개발되었다 이때 양수정내 지하수유입은 1대수층에서만 유래한다. 본 모델은 우물저장과 우물손실을 고려할 수 있어서 양수초기의 지하수위 변화를 정확하게 파악할 수 있는 잇점이 있다. 여러 가지 유동차원에 대해서, 1대수층의 수리전도도가 2대수층보다 더 크고, 비저유율이 2대수층보다 더 낮은 경우 (case 1)와 1대수층의 수리전도도가 2대수층보다 더 작고, 비저유율이 2대수층보다 더 큰 경우 (case 1)에 대해서 표준곡선을 작성하였으며, 이들 곡선의 특성을 분석하였다. 본 모델은 양수정과 관측정의 수위하강이 서로 다른 거동을 보여주는 양수시험자료를 해석하는데 유용하게 쓰일 수 있을 것이다.

• 터널과지하공간
• /
• 제16권4호
• /
• pp.281-287
• /
• 2006

산악 지형 내 터널 굴착에 따른 지하수 유입은 실제 암반 내 발달한 절리 및 파쇄구간을 따라 유동하는 지하수에 의해 발생한다. 이러한 현상을 보다 합리적으로 분석하기 위해 정밀 지반조사 물리탐사 및 시추 탐사 등의 분석 결과를 토대로 모델영역을 선정하고 대표 절리군 1, 2 및 3을 선택하였다. 3차원 절리망을 생성하기 위해 프로그램(FracMan)을 이용하였고 불확실성을 줄이기 위해 확률적 통계기법 중 하나인 Monte Carlo Simulation 을 적용하였다. 수리특성을 파악하기 위해 실시한 불연속면 수리특성 시험 양수시험 및 수압시험 결과 대표 절리군 1, 2 및 3의 투수량 계수($T_f$)는 각각 $9.60{\times}10^{-6},\;7.35{\times}10^{-6}$ 및 $1.01{\times}10^{-5}m^2/s$ 이였다. 모델 영역의 수리적 초기조건을 파악하기 위해 지하수 연속체 프로그램(Visual MODFLOW)를 병행하여 적용하였다. 지하수 불연속체 프로그램(MAFIC)을 이용하여 해석한 결과 터널내부의 단위 길이 당 유입량은 평균 $5.40{\times}10^{-1}m^3/min/km$, 표준편차 $3.04{\times}10^{-1}m^3/min/km$로 예측되었다. 이 결과를 토대로 터널 굴착에 따른 광역적인 지하수계 변화를 분석한 결과 $8.64{\times}10^{-2}{\sim}2.30{\times}10^{-1}m^3/min/km$로 기존 및 본 설계 적용 기준인 3.00 및 $2.00m^3/min/km$의 값보다 낮은 범위로 예측되어 터널 굴착 시 지하수 유입에 따른 공극수압의 변화에 따른 응력 변화가 적어서 터널 구조의 안정성에 미치는 영향이 적을 것으로 분석되었다.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제11권2호
• /
• pp.45-55
• /
• 2006

강변여과수 취수 방식의 도입이 고려되었던 경상북도 고령군 낙동강변의 충적대수층을 대상으로 2차원 지하수 유동모사를 수행하여, 취수정 및 주변 대수층에서의 지하수위와 하천수 유입량의 계절적 변화를 분석하였다. 개념모델의 하천 경계조건은 연평균 하천수위를 고정 수두로 설정하는 경우와 실제 관측된 시계열 자료를 이용하여 변동 수두로 설정하는 경우로 구분하여 모사를 수행하였다. 모사 결과 산정된 하천수 유입 비율은 고정하천수두 모사의 경우 연중 $74.3{\sim}87.0%$의 범위를 보이며 연평균 82.4%인 반면, 변동하천수두 모사의 경우 $52.7{\sim}98.1%$, 연평균 78.5%로 나타났다. 따라서 시기별로 변동하는 하천수위가 하천수 유입 비율에 큰 영향을 미치며, 하천수 유입 비율을 정확하게 산정하기 위해서는 변동하천수두 조건을 이용해야 하는 것으로 해석된다. 강변으로부터 취수정까지의 이격거리를 변화시키면서 모사를 반복 수행한 결과 취수정 이격거리가 커질수록 하천수 혼합 비율은 다소 감소하지만 그 영향이 수 % 정도로 매우 작게 나타났다. 강변여과 취수 시 지하수위 및 하천수 혼합비율에 영향을 미치는 인자인 대수층의 수리전도도 및 저유계수, 지하수 함양량, 취수량 등에 대한 민감도 분석을 수행하였다. 대수층의 평균 지하수위와 하천수 유입 비율은 취수량과 지하수 함양량 변화에 민감한 반면, 취수정의 지하수위는 수리전도도와 취수량에 민감한 것으로 나타났다.

• 지질공학
• /
• 제24권2호
• /
• pp.217-225
• /
• 2014

대규모 수막재배 시설은 동절기 수막재배 기간 동안 군집 취수정의 영향으로 지하수위가 하강하고 이로 인해 지하수 취수량 확보에 어려움을 겪고 있다. 본 연구에서는 동절기 수막재배지역에서의 집중취수에 따른 정밀한 지하수위 변화 분석을 수행하고자 한다. 대상지역인 청원군의 수막재배지역을 포함한 소유역의 지하수 유동을 모의하기 위해 3차원 유한차분모형인 MODFLOW를 구축하였으며, 중요한 입력자료인 지하수 함양량 시계열 자료는 유역수문모형 SWAT모형을 이용하여 구하였다. 여기에 수막재배지의 일최저기온과 취수가동시간과의 관계를 통해 지하수 이용량을 추정한 선행 연구(Moon et al., 2012)를 참고하여 취수계획을 수립, 지하수 모델링에 반영하였다. 부정류 모델링의 적합도를 평가하기 위해 동절기의 모의된 지하수위 강하 패턴을 관측지하수위 자료와 비교하였다. 모델링 기간은 수막재배가 본격적으로 이루어지는 2012년 11월부터 2013년 3월까지의 107일간으로 정하였다. 모의결과 동절기 지하수위 관측값과 비교적 유사한 지하수 하강을 재현할 수 있었고, 일최저기온을 이용한 취수계획방식이 다른 수막재배지역의 지하수 유동해석에도 적용가능할 것으로 판단되었다.

• 물과 미래
• /
• 제29권5호
• /
• pp.215-222
• /
• 1996

본 논문에서는 다공성 매질의 특수율이 비정체형인 경우 대수투수계수-수두 교차공분산에 관한 식을 유도하였으며, 이 교차공분산은 수두분포로부터 특수장의 통계학적 특성을 유추하는데(inverse problem) 매우 중요한 역할을 담당한다. 비정체형 대수투수계수는 일정한 선형경향과 정체형인 미소 변동의 합으로 구성되었으며, 2차원 포화대수층에서 정상 유동문제를 추계학적으로 해석하여 수두분포를 얻었고 이로부터 교차공분산을 유도하였다. 투수계수의 상관함수가 가우스분포를 가지고 그 경향이 수두 경사와 평행이거나 직교하는 두 가지 경우에 대하여 교차공분산을 살펴 본 결과, 투수장의 경향이 주 흐름방향과 평행한 경우 흐름방향 쪽만 제외하고는 정체형임이 밝혀졌다. 또한, 흐름방향과 직교하는 쪽으로의 교차공분산은 정체형 모델 결과와 달리 영이 아님를 알 수 있었다. 따라서 지하수 유동이나 오염물질 확산문제를 다룰 경우, 투수계수장에 어떤 경향이 존재한다고 의심될 때에는 반드시 그 경향을 해석과정에 포함시켜야 한다.

• 한국수자원학회논문집
• /
• 제38권7호
• /
• pp.537-544
• /
• 2005

LPG등 고압가스를 지하에 저장하는 저장공동의 효과적인 관리를 위해서는 지하수위, 가스저장압, 수막공 주입압 등에 따른 공동주변의 유동장 해석, 공동내로의 지하수 유입량 해석을 실시해야 한다. 기존에는 공동의 정확한 형상을 반영하기 위해서 유한요소법이 보편적으로 사용되어 왔으나 한번 설정한 유한요소망으로부터 공동의 설계요소를 변경하는 작업은 수월하지 않아 설계전단계에서 공동 및 수막시설의 다양한 배치에 따른 모의를 수행하는데는 다소 무리가 있다. 이러한 불편함은 경계부의 형상과 조건만으로 내부점에서의 미지변수 계산을 효과적으로 수행할 수 있는 경계요소법을 도입함으로써 극복할 수 있다. 따라서 본 연구에서는 지하공동으로 배수되는 유입량 산정을 위해 경계요소법을 근간으로 한 2차원 지하수 흐름모형을 구성하였으며 지하공동으로 배수되는 정확한 유입량 산정을 위해 해석해 및 기존의 유한요소해석결과와 비교하여 수치해석의 타당성을 제시하였다. 이 모형은 공동의 설계단계에서 수막공의 간격과 배열형태, 공동의 최적형상을 결정하기 위한 모의실험에서 기존의 유한요소해석의 단점을 극복하는 대안으로 활용될 수 있을 것으로 판단된다.

• 터널과지하공간
• /
• 제21권1호
• /
• pp.66-81
• /
• 2011

암반 내 균열 생성, 진전, 파괴 등과 같은 지하 암반의 역학적 거동과 이들 균열을 통한 지하수 유동 및 온도 변화에 기인한 열응력이 역학적 거동에 미치는 상호작용을 모델링하기 위한 수치해석코드를 개발하였다. 개발된 수치해석코드에서는 기존의 2차원 FRACOD(Shen & Stephasson, 1993)에 열-역학 및 수리-역학 상호 거동을 모델링하기 위한 해석모듈을 개발하여 추가하였다. 열-역학 연계를 위해서는 가상열원법과 시간전진기법을 도입하였으며, 수리-역학 연계에서는 양해법에 의한 반복기법을 적용하였다. 수치해석결과와 해석해와의 비교를 통해 개발된 해석모듈의 유용성을 검증하고 해석사례를 통해 온도변화에 따른 열균열 발생 및 수압파쇄과정에서의 균열 진전 양상과 같은 절리암반 복합거동을 잘 표현할 수 있음을 확인하였다.