• 한국해양환경ㆍ에너지학회지
• /
• 제11권4호
• /
• pp.191-198
• /
• 2008

$CO_2$ 해양지중저장 처리 시스템 설계를 수행하는데 있어 전산모사를 통한 공정 설계는 필수적이다. 즉, 수치 모델링을 이용하여 $CO_2$ 해양지중저장 처리 공정 중 발전소 및 제철소 등 대규모 발생원에서부터 $CO_2$를 포집하는 포집공정, 포집한 $CO_2$를 파이프라인이나 선박 등을 통해 이송하는 수송공정, 이를 깊이 800m 이상의 해저 지질구조내 대규모로 수백-수천년 이상 장기간 저장하는 저장공정 등 일련의 공정을 열역학 상태방정식 등을 이용하여 모사하는 것이다. 본 논문에서는 $CO_2$ 해양지중저장 처리를 위한 공정 설계에 사용되는 열역학 상태방정식들을 비교 분석하고 이들이 압축 및 수송공정에 미치는 영향 등을 평가하였다. 이와 같은 연구결과에 근거하여 $CO_2$ 해양지중저장 처리를 위한 압축 및 수송 공정 설계 시 유용한 열역학 상태방정식을 제안하였다. 상태방정식 계산결과 간의 비교 및 영향 평가를 위하여 순수 단일성분 $CO_2$와 순산소 연소 석탄 화력발전소에서 포집된 $CO_2$ 혼합물의 압축 및 수송거동을 분석하였다. 순수 단일성분 $CO_2$의 압축 및 수송공정 계산결과에는 이상기체 상태방정식을 제외한 상태방정식 간에 큰 차이가 존재하지 않으나, NO, Ar 그리고 $O_2$ 등이 포함된 순산소 연소 석탄 화력발전소에서 포집된 $CO_2$ 혼합물의 압축 및 수송공정에는 PR 계열의 상태방정식과 BWRS 상태방정식에서 커다란 차이가 나타났다. 즉, 8bar 이상의 압력, $30^{\circ}C$ 내외의 온도 영역에서 $CO_2$ 혼합물을 서로 다른 상태로 예측하였다. 이러한 불순물의 영향으로 인하여 $CO_2$ 혼합물의 압축 및 수송공정 설계에 PR계열의 PRBM 상태방정식의 적용이 BWRS식 적용보다 유용하다고 판단된다.

• 지구물리와물리탐사
• /
• 제9권1호
• /
• pp.60-66
• /
• 2006

세 개의 서로 다른 사암 샘플들 -두 개의 합성 샘플과 한 개의 현장 샘플- 에 대해 현장 저류층의 대표적인 구속압력과 공극압력하에서 초음파 시험이 수행되었다. 세가지 사암 샘플들은 (a) 캘사이트 시멘트(calcite intergranula. cement (CIPS))로 만든 합성사암, (b) 실리카 시멘트(silica intergranular cement)로 만든 합성 사암 (c) Otway Basin 의 Boggy Creek 1 시추공에서 시도되는 $CO_2$ 파일럿 프로젝트의 대상 암석층 중 Waarre 층으로부터 추출한 코아 샘플로 구성되어 있다. 공극률은 각각 32%, 33%, 26%이다. 초기시험은 실내건조(room-dried) 상태에 있는 코아들에 대해 구속응력을 5 MPa 씩 단계별로 65 MPa 까지 증가시키며 이루어졌다. 그리고 나서 모든 코아들에 처음에는 온도 $22^{\circ}C$에서 6 MPa 공극압력으로 기체상의 $CO_2$를, 그 다음에는 온도 $22^{\circ}C$ 에서 7 MPa 부터 17 MPa 까지 5 MPa 씩 증가시키면서 액체상의 $CO_2$를 주입하였다. 구속응력은 10MPa부터 65 MPa까지 달리 하였다. P와 S 초음파 파형들이 유효응력이 증가할 때마다 기록되었다. 속도-유효응력 반응들이 P 파와 S 파에 대해 실험 자료들로부터 계산되었으며, 감쇠(1/Qp)들은 스펙트럼 비 방법을 이용하여 파형들로부터 계산되었다. 각각의 사암들에 대한 이론적인 속도-유효응력 계산은 $CO_2$ 압력-밀도 와 $CO_2$ 체적계수-압력 상 다이어그램(phasediagram), Gassmann 유효 매질 이론(effective medium theory)을 이용하여 구하였다. 기체상의 $CO_2$ 주입은 속도-유효응력에서 건조상태(공기로 포화된 상태)에 비해 거의 무시할만한 변화를 가져왔다. 다양한 공극압력에서 액체상의 $CO_2$ 주입은 공기로 포화된 상태에 비해 속도-유효응력 반응을 평균 약 8% 정도 낮게 한다. 실험자료들은 높은 유효응력에서 Gassmann 계산들과 일치한다. 이러한 이론과 일치하는 "임계" 유효응력은 사암의 종류에 따라 달라진다. 이 차이는 각각의 사암 종류의 미세구조에서 미세 균열 수의 차이에 기인한 것이라 생각된다. 높은 유효응력에서의 이론과 의미있게 일치하였으며, $CO_2$ 주입 시 현장에서의 탄성파 거동을 예상하는데 있어서 어느 정도 확신을 준다.

• Korean Chemical Engineering Research
• /
• 제50권6호
• /
• pp.980-987
• /
• 2012

비점오염원은 세정에 의해 하수관 내로 유입되어 수계까지 흘러 수계의 오염부하를 가중시킬 뿐만 아니라, 관의 상태에 따라 누수 및 월류를 야기하여 오염을 일으킨다. 이에 따라 하수관거의 유지관리 및 환경오염 방지를 위한 관거내 하수의 유량, 수질, 불명수 및 월류 등을 효율적으로 관리하기 위한 모니터링의 필요성이 증가하고 있다. 그러나 하수관거의 경우, 지하에 설계되며 그 구조 및 연결이 복잡한 특성으로 인해 실제 하수관거에 대한 모니터링은 쉽지 않다. 본 연구에서는 시스템해석이론에 기초한 하수관망 오염제거 해석과 관망 모니터링 및 이상진단방법을 제시하였다. 먼저 하수관망 공정모사 프로그램인 Stormwater & Wastewater Management Model for expert (XP-SWMM)을 이용하여 관망 내의 오염물질의 거동패턴을 해석하였다. 둘째, 다변량 통계 모니터링을 이용하여 하수관망 내의 수질 모니터링 및 하수관망 유출을 탐지하는 관망 이상 진단을 수행하였다. 정적/동적 상태 시스템에 기초한 하수관망 오염 매카니즘 해석결과, 강우시 총질소와 총인 부하량이 비강우시보다 급격하게 증가함을 확인하였으며, 이는 수계의 오염부하를 가중시킬 것으로 판단되다. 이에 따라 하수관망 내 유출은 강우로 생긴 유량 및 유입오염물질의 농도 증가로 인한 관망이상으로 사료된다. 제안된 하수관망 모니터링 및 이상진단 기법은 도시 유역에서의 비점오염원 관리와 지속적인 모니터링에 있어 효과적으로 사용될 수 있을 것으로 예상된다.

• 대한토목학회논문집
• /
• 제5권1호
• /
• pp.21-30
• /
• 1985

Tension Leg Platform (TLP)이란 평행위치로부터 일정 범위내에서 움직임으로 인하여 외 력의 효과를 완화시키는 compliant 구조물인 동시에, 기인장력을 받고 있는 연직 anchor cable 이 있으므로 부력이 자중을 초과하게 되는 안정한 platform 이다. 일반적으로 부체는 해상조건이 험할수록, 그리고 수심이 깊어질수록 동요가 심해지는데 TLP는 기인장 cable로 인하여 심해에서도 비교적 동요가 작아서 최근 대수심구조물의 총아로 각광받고 있다. 일찌기 Paulling 등이 TLP 거동의 예측을 위하여 수정된 Morison 방정식을 사용하는 선형동유체력합성방법을 발표하였다. 그러나 만일 TLP의 각 부재가 Morison 방정식의 가정이 성립할 수 없을 정도로 크다면 새로운 해석이 필요하다 하겠다. 일본의 Tanaka는 이런 경우에 McCamy-Fuchs 이론의 결과치를 이용하였으나, 완전한 해석이라기 보다는 일종의 간편법이라 하겠다. 본고에서는 큰 배수용적을 가진 연직부체가 있고, 이론적 해석의 결과를 검토해 볼 수 있는 수리모형 실험 결과가 있는 Deep Oil Technology (DOT) 회사의 TLP를 대상으로 하였다. 이 TLP는 부력을 전담하고 있는 연직축대칭 원통과 이들을 연결하고 있는 세부재로 이루어져 있어 축대칭부분에는 축대칭 Green 함수를 사용하여 동유체력을 구하고 세부재는 종래의 수정된 Morison 방정식의 항력항을 선형화하여 동유체력을 구하였다. 그리하여 부재의 각 미소부분에서 구한 힘들을 TLP의 중심에 원점을 둔 좌표계로 옮겨 동적응답을 구한 것이다. 본 해석에서 부재 상호간의 작용은 무시하였으며 단지 부재간의 거리효과만 고려하였다. 따라서 사용된 좌표계는 전체 (Global) 좌표계, 지점 (Local) 좌표계 및 파랑 (Wave) 좌표계 등이었고 각 좌표계간의 변환식이 필요하였다. 전체적인 해석정도는 선형이론으므로 케이블의 강성도 역시 선형적으로 구하였으며, 앞서 언급했다시피 Morison 방정식의 비선형항인 항력항은 Fourier 해석으로 선형화 하였다. 이러한 Fourier 해석은 잘 알려져 있는 Lorentz 원리와 같다고 볼 수 있다. 세부재의 경우 접선력은 무시하였고 수입자의 운동에 의한 부채에 대한 수직력만 고려하였다. 여기서 파랑좌표계에서 지점좌표계로의 좌표변환이 주의를 요하고 있다. 이제 이렇게 구한 각 힘들을 전체좌표 계에서 6개의 자유도별로 운동방정식에 대입하면 각 자유도별 동적응답이 구하여지는 것이다. DOT TLP의 Surge mode에 대한 동적응답을 실험치와 비교하여 본 결과, 세부재에 대한 고려를 뺄 수 없음을 알 수 있었다. 이는 연직축대칭 부체의 크기가 그리 크지 않으므로 인한 것이며, TLP의 원형의 경우에는 보다 더 관성력이 지배적일 것으로 사료된다.

• 대한토목학회논문집
• /
• 제4권2호
• /
• pp.133-142
• /
• 1984

횡방향(橫方向)과 종방향(縱方向), 대입열용접(大入熱鎔接)과 소입열용접(小入熱鎔接), 모재(母材)와 같은 강도(强度)의 용접재료(鎔接材料) 사용(使用)과 모재(母材)보다 낮은 강도(强度)의 용접재료(鎔接材料) 사용(使用) 등으로 서로 비교되는 용접(鎔接)이음의 소인장시험편(小引張試驗片)을 만들어서 피로시험(疲勞試驗)을 행하여 피로구열(疲勞龜裂) 진전속도(進展速度) ${\frac{da}{dN}}$와 구열선단(龜裂先端) 부근의 응력확대계수(應力擴大係數)의 변동범위(變動範圍) ${\Delta}K$와의 관계곡선(關係曲線)을 그려서 비교시험간(比較試驗間)의 모재(母材), 열영향부(熱影響部) 및 용착금속부(鎔着金屬部)로 구분, 혹은 지금까지의 연구자료(硏究資料) 등과 비교검토하였다. 이 결과, 다음과 같은 현상들을 알 수 있었다. 기본적(基本的)으로, 용접방향(鎔接方向), 용접입열량(鎔接入熱量), 용접재료(鎔接材料)의 강도(强度), 혹은 모재(母材), 열영향부(熱影響部) 및 용착금속부(鎔着金屬部)의 구분 등에 따라 ${\frac{da}{dN}}-{\Delta}K$관계에 큰 차이가 없었다. 다만, 첫째, 소재(素材)에 대한 경우에 비해 대개 처음은 같은 ${\Delta}K$에서 ${\frac{da}{dN}}$가 상당히 늦다가 점차 증가하여 중간쯤에서 같아진 후 끝 부분에서 같은 ${\Delta}K$에서${\frac{da}{dN}}$가 다소 빨라짐을 알 수 있었다. 둘째, 열영향부(熱影響部)에서 용접금속부(鎔接金屬部)로 진전(進展)하면서 ${\frac{da}{dN}}$가 다소 늦어지는 것을 알 수 있었다. 셋째, 용접방향(鎔接方向)과 구열방향(龜裂方向)이 평행인 경우가 직각인 경우에 비해, 모재(母材)와 같은 용접재료(鎔接材料)를 쓴 경우 소입열용접(小入熱鎔接)인 경우가 대입열용접(大入熱鎔接)인 경우에 비해, 모재(母材)보다 낮은 용접재료(鎔接材料)를 쓴 경우 대입열용접(大入熱鎔接)인 경우가 소입열용접(小入熱鎔接)인 경우에 비해, 대입열(大入熱)의 평행용접(鎔接)의 경우를 제외한 모든 용접(鎔接)에서 모재(母材)와 같은 강도(强度)의 용접재료(鎔接材料)를 쓴 경우가 모재(母材)보다 낮은 강도(强度)의 용접재료(鎔接材料)를 쓴 경우에 비해 ${\Delta}K$가 낮은 시기(時期)에 일찍 저속(低速)으로 구열(龜裂)이 시작되어 ${\Delta}K$의 큰 증가없이 ${\frac{da}{dN}}$가 빠르게 증가한 다음 다른 경우와 같은 성상(性狀)으로 진전됨을 알 수 있었다. 넷째, 소입열용접(小入熱鎔接)의 경우 평행용접(鎔接)이 직각용접(鎔接)에 비해, 소입열용접(小入熱鎔接)이 대입열용접(大入熱鎔接)에 비해 같은 ${\Delta}K$에서 ${\frac{da}{dN}}$가 다소 늦은 것을 알 수 있었다.

• 대한화학회지
• /
• 제37권11호
• /
• pp.951-960
• /
• 1993

수용액상에서 $Hg^{2+}$에 의한 trans-[Co(3,2,3-tet)$X_2]^+$ (3,2,3-tet = 4,7-diazadecane-1,10-diamine, $X_2\;=\;Cl_2,\;(NO_2)Cl,\;Br_2,\;(NO_2)Br,\;(NO_3)_2)$ 착물의 아쿠아 반응이 연구되었다. 생성물을 확인하고 반응메카니즘을 추정하기 위하여 크로마토그래피를 사용하였고 전자흡수 스펙트럼을 측정하였다. 그 결과 네자리 리간드인 3,2,3-tet가 배위된 여러가지 $trans-[Co(3,2,3-tet)X_2]^+$ 착물은 각각 아쿠아된 trans-[Co(3,2,3-tet)$(OH_2)_2]^{3+}$ 착물을 거쳐 cis-${\beta}$-[Co(3,2,3-tet)$(OH_2)_2]^{3+}$ 착물이 생성되었다. $Hg^{2+}$에 의한 trans-$[Co(3,2,3-tet)Cl_2]^+$ 착물과 trans-[Co(3,2,3-tet)$(NO_2)Cl]^+$ 착물의 아쿠아 반응에 대한 메카니즘을 추정하기 위하여 속도론적 조사를 하였따. 그 결과 $trans-[Co(3,2,3-tet)Cl_2]^+$ 착물은 D(dissociative)-메카니즘으로 진행되었고, trans-[Co(3,2,3-tet)$(NO2_)Cl]^+ $착물은 $I_d$(interchange dissociative)-메카니즘으로 진행되었다. 그리고 입체화학적인 거동을 조사하기 위하여 라세미(R,R:S,S)3,2,3-tet 대신에 키랄성이 R,R인 3,2,3,-tet를 배위시킨 trans-$[Co(R,R-3,2,3-tet)Cl_2]^+$ 착물에 $Hg^{2+}$를 용리시켰을 때 아쿠아 반응에 대한 원편광이색성(circular dichroism) 스펙트럼을 측정하여 그 절대구조를 확인한 결과 ${\Delta}-cis-{\beta}$-[Co(R,R-3,2,3-tet)$(OH_2)_2]^{3+}$ 착물이 생성되었다.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2007년도 추계학술대회 논문집
• /
• pp.551-557
• /
• 2007

본 논문은 도시철도공사 시 하저통과구간의 지반투수저하공법에 대한 연구결과이다. 본 연구에서는 하저지반의 투수저하공법들 중 특히 주입공법에 대하여 연구의 초점을 두었다. 이에 따라 하저지반 투수저하를 위하여 적용되고 있는 기존주입공법들의 이론적 배경과 시공방법들을 문헌과 시공사례자료 등을 고찰하여 비교 분석하였다. 이와 같이 다양한 고찰을 통하여 주입공법들을 하저통과 도시철도터널 안전시공 측면에서 평가하였다. 또한 분석 결과로부터 하저지반에서 장단기적으로 안전하게 지반수리거동을 도모하기위한 주입공법의 설계기법에 대하여 제시하였다. 이와 같이 하저지반 투수저하주입공법들의 평가결과를 근거로 하저통과 터널공사 시 안전시공과 완공 후 터널구조물의 장기적인 안전유지관리 측면에서 적용성이 양호하다고 판단되는 N.D.S(Natural Durable Stabilizer)공법과 S.M.I(Space-Multi- Injection Grouting) 공법을 제시하였다. 그러나 이 두공법의 실질적인 시공효과를 평가하기위하여 OO지역에 건설예정인 하저통과 도시철도터널공사구간에서 시험시공을 실시하였다. 시험시공기간 중에는 실내 및 현장에서 지반투수특성시험을 실시하여 시공전후의 지반의 투수성저하효과를 평가하였으며, 두공법의 시험시공결과는 이론적인 배경과 함께 비교 분석하였다. 결론적으로 S.M.I공법의 경우 N.D.S공법보다 투수성 저하효과가 큰 것으로 나타났다. 즉, 주입공법에 있어서 시멘트 무기질계 재료보다는 비 알카리성 실리카 졸계가 투수저하 효과를 더 크게 할 수 있을 것으로 판단되었다. 끝으로 본 연구결과로부터 특히 철도터널설계시 하저지반의 투수저하 증대를 위한 주입공법 설계기술발전에 공헌 할 것이다. 값을 가지며 표준 편차가 작았다 PR와 $SpO_2$는 두 그룹 모두 정상범위내의 값을 보였다. 점막하 주사를 이용한 Midazolam의 추가 투여는 안전성을 크게 침해하지 않는 범위 내에서 진정 깊이를 깊게 해주며 안정적인 진정을 이루게 함으로써 보다 바람직한 진정을 유도해 줄 수 있는 것으로 사료된다.보이지 않았다.tag의 표면에 형성된 뚜렷한 측부 가지들도 발견되었다.}60cm$에 위치시켜 표면선량의 over dose나 under dose의 발생 방지를 위한 주의가 반드시 필요하겠다.군의 전단결합강도가 가장 낮았다.본에서 III군에 비해 크기가 큰 기포가 더욱 많이 관찰되었다.e의 약동학은 cimetidine에 의해 유의한 차이를 보였으며 CYP1A2유전자형에 따른 영향은 관찰할 수 없었다. CYP1A2유전자형에 따른 생체내 대사능을 관찰하는 실험이 향후 이루어 져야 할 것으로 사료된다.san film보다 큰 수증기 투과도를 보였다.적으로 유의한 차이를 보이지 않았다.y tissue layer thinning은 3 군모두에서 관찰되었고 항암 3 일군이 가장 심하게 나타났다. 이상의 실험결과를 보면 술전 항암제투여가 초기에 시행한 경우에는 조직의 치유에 초기 5 일정도까지는 영향을 미치나 7 일이 지나면 정상범주로 회복함을 알수 있었고 실험결과 항암제 투여후 3 일째 피판 형성한 군에서 피판치유가 늦어진 것으로 관찰되어 인체에서 항암 투여후 수술시기는 인체면역계가 회복하는 시기를 3주이상 경과후 적어도 4주째 수술시기를 정하는 것이 유리하리라 생각되었다.한 복합레진은 개발의 초기단계이며, 물성의 증가를 위한 연구가 필요할 것으로 사료된다.또 다른 약물인 glyc

• 마이크로전자및패키징학회지
• /
• 제12권4호통권37호
• /
• pp.307-313
• /
• 2005

선형구조의 시편을 사용하여 Sn96.5Ag3.0Cu0.5의 electromigration 특성을 살펴보고 공정 조성의 SnPb의 electromigration특성과 비교, 분석하였다. SnAgCu의 electromigration에 관한 특성 중 시간에 따른 물질의 이동 양상과 여러 가지 electromigration 매개변수 (활성화 에너지, 임계 전류밀도, 확산계수와 유효전하수의 곱)을 살펴보았다. 임계 전류 밀도는 $140^{\circ}C$에서 $2.38{\times}10^4A/cm^2$ 이고 이 값은 $140^{\circ}C$에서 electromigration에 의한 물질 이동이 발생하지 않는 최대 전류 밀도를 나타낸다. 활성화 에너지는 $110-160^{\circ}C$ 온도 범위에서 0.56 eV가 측정되었다. DZ$\ast$의 값은 $110^{\circ}C$에서 $3.12{\times}10^{-10}\;cm^2/s$, $125^{\circ}C$에서 $4.66{\times}10^{-10}\;cm^2/s$, $140^{\circ}C$에서 $8.76{\times}10^{-10}\;cm^2/s$, $160^{\circ}C$에서 $2.14{\times}10^{-9}cm^2/s$ 이었다. 그리고 SnAgCu와 공정 조성의 SnPb 물질의 electromigration 거동 특징은 크게 다른데, SnPb의 경우 음극에서 보이드(void) 형성이 발생하기 전에 잠복 시간이 존재하고 SnAgCu의 경우 잠복 시간이 존재하지 않는다는 점이다. 이는 각 원소들의 확산 기구(diffusion mechanism)의 차이에 기인한 것이라 생각된다.

• 대한토목학회논문집
• /
• 제38권2호
• /
• pp.357-367
• /
• 2018

최근 경주와 포항에서 발생한 지진으로 국내에서 액상화 현상에 대한 관심이 커지고 있다. 지반의 액상화는 포화된 상태에서 지진과 같은 동하중을 받았을 때 과잉간극수압이 발생하여 흙이 강도를 상실하고 물과 같이 거동하는 현상이며 지반 침하와 상부구조물의 전도와 같은 심각한 문제를 야기한다. 따라서 액상화 발생 가능성을 미리 파악하고 대비할 필요가 있다. 액상화의 발생 가능성과 액상화 피해 정도는 일반적으로 액상화 가능 지수(Liquefaction Potential Index, LPI)에 의해 정량적으로 평가된다. LPI의 계산은 시추공 별로 이루어지며 지반응답해석이 필수적인 작업으로 선행되어 많은 시간과 노력이 필요하다. 본 연구에서는 다양한 지하수위 분포를 가지는 넓은 지역의 액상화 평가를 간단히 수행할 수 있도록 전단파 속도와 LPI의 상관관계를 이용한 액상화 평가 방법을 제안하였다. 제안된 방법은 액상화 가능 층의 평균 전단파 속도(${\bar{V}}s^{\prime}liquefiable$)와 액상화 가능 층의 두께로 나누어 정규화한 정규화 LPI의 상관관계를 분석하여 지하수위 별로 다양한 암반노두가속도에 대해 적용 가능한 상관관계식을 제시하고 이용한다. 상관관계를 이용한 액상화 평가 방법의 적용성을 확인하기 위해 서울특별시의 104개 시추조사자료를 이용하여 지하수위 0m, 1m, 2m, 3m에 대해 상관관계식을 제시하였으며 제시한 상관관계식을 이용하여 서울특별시와 경주시의 액상화 발생 가능성을 평가하였다. 지반응답해석을 이용해 계산한 LPI와 상관관계식을 이용해 계산한 LPI를 비교하였으며 제안된 액상화 평가 방법의 적용성을 확인하였다. 마지막으로 제안된 액상화 평가 방법에 따라 결정된 LPI의 분포를 지구통계학적 기법인 크리깅을 통해 지도로 나타내었다.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제9권1호
• /
• pp.79-86
• /
• 2004

암석 내 균열을 따른 수리전도도는 균열의 기하학적 요소, 즉 방향, 간극, 거칠기 그리고 상호 연결도에 주로 좌우된다. 따라서, 균열 내 투수계수를 정확하게 계산하기 위해서는 이와 같은 기하 요소들을 최대한 계산모델에 반영할 필요가 있다. 이 연구에서는 균열 기하양상을 최대한 정확히 반영한 균열모델에서 기존 수치해석과는 다른 새로운 방법인 균질화 해석법(homogenization analysis method)을 이용하여 균열을 따른 투수계수를 구하기 위해 수치해석을 수행하였다. 먼저, 공초점 레이저 스캔 현미경(Confocal Laser Scanning Microscope)을 이용하여 암석시료의 균열 조도와 균열에 가한 수직압축력의 변화에 따른 간극 변화량을 직접 측정하고, 이와 같이 획득한 자료는 균열모델 재현을 위한 입력자료로 사용되었다. 재현된 균열모델을 토대로 한 균질화 해석법은 미시규모(microscale) 매질특성과 거시규모(macroscale) 매질특성을 동시에 고려하여 투수계수를 계산할 수 있는 것이다. 즉, 균질화 해석법은 주기적 미세구조(microstructure)를 갖는 미소 불균질 물질의 거동특성을 구명하기 위해 개발된 새로운 형태의 섭동(perturbation) 이론이다. 이는 균질한 미시규모에서 미시 투수특성을 계산한 후, 거시규모에서의 균질화 투수계수를 계산하게 된다. 그러므로, 이 방법은 균열 기하양상의 국부적 영향을 고려한 투수특성을 정확히 해석할 수 있다. 균질화법을 이용한 투수계수 산정결과를 기존 연구에서 제안한 경험식과 비교하여 그 타당성을 검증하기 위해 전술한 2차원 균열모델을 이용한 투수계수 계산을 수행하였다. 균열모델은 거칠기(roughness)를 반영하고 동일한 간극을 할당한 평행판 모델을 가정하였다. 계산결과에 의하면, 균질화 해석법에 의해 계산한 C-투수계수는 실내투수시험에 의해 구한 투수계수와 같은 범위의 값을 가지거나 $10^1$ 정도의 차이를 보여, 그 계산결과는 타당하다고 볼 수 있다. 그러나, 균질화 해석법은 국부적으로 불균질한 균열 기하양상과 물질특성이 미시규모와 거시규모에서 모두 고려되므로, 이들 특성을 정확히 알고 있을 경우 기존에 제안된 경험식들에 의한 계산결과 보다 균질화 해석법의 결과가 훨씬 정확함을 주목하여야 한다.