• 대한방사선치료학회지
• /
• 제29권2호
• /
• pp.9-17
• /
• 2017

목 적: 방사선치료에서 CT number에 의해 나타나는 조직의 묘사 및 전자밀도는 CT(Computed Tomography)기반의 전산화치료계획 정확성을 보장하는데 중요한 역할을 한다. 하지만 체내의 금속 이식물은 CT number의 정확성을 감소시킬 뿐 아니라 조직 묘사에 대한 불확실성을 나타내기 때문에 임상에서는 metal artifact를 감소시킬 수 있는 알고리즘이 개발되었다. 이에 본 연구에서는 본원에서 사용하고 있는 GE사의 SMART MAR의 CT number 정확도를 평가하고 방사선치료에서의 유용성에 대해 평가하고자 한다. 대상 및 방법: 영상평가를 위해 CIRS ED phantom을 이용하여 6개 조직의 rod를 삽입하여 동일한 체적의 ROI를 형성 후 original 영상을 획득하고 의료용 티타늄 rod를 삽입하여 non-SMART MAR 영상과 SMART MAR영상을 획득하여 CT number와 SD값을 비교하였다. Metal artifact로 인해 CT number 변화에 따른 선량변화 확인을 위해 전산화계획시스템 Epclipse를 사용하여 CIRS ED phantom CT 영상에 PTV를 형성하여 original 영상을 획득하고 편측 티타늄 삽입, 양측 티타늄 삽입 영상을 획득하여 non-SMART MAR와 SMART MAR영상에 동일한 치료계획을 수립하여 PTV가 받는 평균처방선량, HI(Homogeneity Index), CI(Conformity Index)를 비교, 분석하였다. 흡수선량 측정은 원통형 아크릴 팬텀과 0.125 cc ionchamber, electrometer를 이용하여 선량변환상수(cCy/nC)를 계산하고 CIRS phantom을 이용하여 편측, 양측 티타늄 rod를 삽입한 영상으로 non-SMART MAR와 SMART MAR 영상을 획득하여 동일한 지점에서의 흡수선량을 측정하여 전산화치료계획상의 point dose와 비교하였다. 결 과: 영상평가 결과 CT number는 non-SMART MAR영상보다 SMART MAR 영상이 original영상에 더 유사한 값이 나왔고 SD값은 SMART MAR영상에서 더 감소되었다. 선량평가 결과 평균처방선량과 HI 및 CI 값은 SMART MAR 영상보다 non-SMART MAR 영상이 original 영상에 더 근접한 결과가 나왔지만 통계적으로 유의하지 않았다. 흡수선량 측정결과 치료계획상의 point dose와 실제 흡수선량과의 차이가 non-SMART MAR의 영상에서는 각각 2.69, 3.63 %의 차이가 있었지만 SMART MAR영상에서는 0.56, 0.68 %로 감소하였다. 결 론: 금속 이식물을 삽입한 환자의 CT 영상에 SMART MAR를 적용했을 때 CT number 정확성 상승 및 SD 감소로 영상의 질이 향상되므로 종양과 정상조직의 윤곽도 생성 및 선량계산 시 유용할 것으로 사료된다.

• 대한환경공학회지
• /
• 제34권4호
• /
• pp.223-231
• /
• 2012

본 연구는 유 무기성 폐기물 처리를 위한 다축스크류 난류접촉식 고효율 건조기의 최적 설계를 위한 연구의 일환으로 실험적 연구와 수치해석적 연구를 통하여 건조기 내부의 열유동 메카니즘을 규명하고 건조효율을 높일 수 있는 설계 기준을 제시 하고자 수행되었다. 사용된 건조기의 대표적인 특성은 건조용 가스의 바이패스 시스템으로 연소를 통해 얻어진 고온의 가스가 다축 스크류 내부 축을 따라 흐르면서 분공을 통해 고속으로 슬러지 내부로 분출되면서 열풍 건조를 하게 되는 것이다. 다양한 열원의 적용이 가능하고 고온 난류 분사식으로 높은 건조속도를 갖고 있으며 고점성 물질에도 적용이 가능한 것이 장점이다. 여기서 건조가스의 분배는 슬러지를 체적 가열하는 한편 슬러지가 뭉치지 않고 지속적으로 열과 물질 전달을 원활하게 하는 기공을 유지하도록 해야 한다. 실험결과 하수슬러지 200 kg/hr를 처리하는데 스크류를 1 rpm으로 회전시킬 때 적정 체류시간은 100분 정도로 나타났다. 또한 다양한 열량 공급 결과 150,000 kcal/hr로 공급한 경우 높은 건조효율을 유지하면서 잉여 열량공급으로 인한 과도한 열피로 및 열량의 낭비를 줄일 수 있는 것으로 나타났다. 유 무기성 슬러지의 설계기술 및 건조효율 향상을 위한 건조가스의 유동 및 온도분포를 수치해석적 연구를 통하여 계산하였으며 연소실과 건조실의 온도계산 결과는 실험 자료와 매우 유사하게 나타나서 성공적으로 비교검증을 수행하였으며, 향후 물질전달에 대한 세부적인 모델을 적용하여 연구를 지속적으로 수행할 예정이다.

• 대한토목학회논문집
• /
• 제3권2호
• /
• pp.1-15
• /
• 1983

본(本) 논문(論文)에서는 활성(活性)슬러지법(法)에 의(依)한 분뇨처리방법(糞尿處理方法) 대신 BFB를 이용(利用)한 저희석(低稀釋) 분뇨처리(糞尿處理)의 가능성(可能性)과 수학적(數學的) 해석(解析)을 시도(試圖)하였다. 실험(實驗)은 BFB의 메디아로써 폴리프로필렌계(系) 합성섬유(合成纖維)인 부직포(不織布)를 사용(使用)하였으며, 수온(水溫) $20^{\circ}C$에서 F/M비(比) 및 분뇨(糞尿)의 희석비(稀釋比)를 0.12~0.37/day와 2~10배(倍)로 변화(變化)시켜 연속식(連續式) 반응조(反應槽)에 의해 수행(修行)되었다. 여기서 얻어진 결과(結果)로는 저희석(低稀釋) 분뇨처리(糞尿處理)에 대(對)한 BFB의 적용(適用)은 다른 생물학적(生物學的) 처리법(處理法)에 비(比)해 더 좋은 효과(効果)를 얻을 수 있는 것으로 나타났는데 본연구(本硏究)에서 나타난 분뇨(糞尿)의 최적희석비(最適稀釋比)는 5배(倍)로 나타났고, 그 때의 최적(最適) 체류시간(滯留時間)은 17시간(時間)으로 나타났다. 따라서 반응조(反應槽)의 체적(體積)은 70%, 희석수(稀釋水)는 80%까지 절감(切減)할 수 있는 것으로 계산(計算)되었다. 또 완전혼합(完全混合) 활성(活性)슬러지법(法)에 적용(適用)되는 수학적(數學的) 해법(解法)으로도 BFB의 실험결과(實驗結果)를 해석(解析)할 수 있었다. $20^{\circ}C$에서 McKinney 식(式)($K_m$)과 Eckenfelder 식(式)($K_e$)을 이용(利用)하여 구(求)한 유기물(有機物) 제거속도(除去속도(速度))는 1.784/hr와 $2.0{\times}10^{-3}l/mg{\cdot}day$이었으며, 다른 방법(方法)에 비(比)해 유기물(有機物) 제거(除去)가 빠른 속도(速度)로 진행(進行)되었다. 미생물(微生物) 번식계수(繁殖係數)($a_5$), 내호흡(內呼吸) 계수(係數)(b), 미생물(微生物) 합성(合成)에 필요(必要)한 요소(酸素) 요구율(要求率)($a{_5}^{\prime}$) 및 내호흡(內呼吸)에 필요(必要)한 산소(酸素) 요구율(要求率)(b')은 각각(各各) 0.349, 0.0237/day, 0.495와 0.0336으로 계산(計算)되었다.

• 한국의학물리학회지:의학물리
• /
• 제24권3호
• /
• pp.154-161
• /
• 2013

본 연구는 총 10명의 전립선 암 환자를 대상으로 세기조절방사선치료(IMRT), 균일스캐닝양성자치료(USPT), 그리고 세기조절양성자치료(IMPT)기술을 이용한 치료계획의 결과를 비교, 평가 하였다. 각 치료 계획은 타깃 체적의 95%에 70 Gy가 28회 분할 조사되도록 하였으며 세기조절방사선치료(IMRT)에서는 step-and-shoot 기법을 이용하여 총 7개의 빔을 사용하여 방사선을 조사하였고, 균일세기양성자치료(USPT)와 세기조절양성자치료(IMPT)에서는 동일한 방사선 가중치의 측방향대향조사면(lateral opposing field)를 사용하여 타깃에 처방선량이 전달되도록 하였다. 한편, 세기조절양성자치료(IMPT)의 최적화를 위해 IMRT치료와 유사한 Inverse planning을 수행하였다. 결과 비교를 위해 타깃의 균질성지수(homogeneity index) 및 동형지수(conformity index)와 정상조직의 정상조직합병증확률(NTCP)을 계산하였다. 비록 치료기법간에 균질성지수(homogeneity index), 동형지수(conformity index)차이가 크지 않았지만, 직장의 경우 각 세기조절방사선치료(IMRT), 균일스캐닝 양성자치료(USPT) 및 세기조절양성자치료(IMPT)에서 2.233, 3.326 및 1.707로 계산되었다. 또한 방광의 정상조직합병증확률(NTCP)는 0.008, 0.003, 및 0.002를 나타내었다. 직장과 방광의 NTCP 값이 IMPT을 사용할 때 유의하게 낮은 값을 보이는 것을 확인하였다. 본 연구를 통해 전립선 암의 방사선 치료 시 세기조절방사선치료(IMRT)보다 양성자를 이용한 방사선 치료, 특히 최적화된 세기조절양성자치료(IMPT)가 치료 효과를 높일 수 있는 치료계획이 될 수 있음을 확인할 수 있었다.

• 한국의학물리학회지:의학물리
• /
• 제17권1호
• /
• pp.24-31
• /
• 2006

유방암 환자의 방사선치료에 있어 치료도중(intrafractional) 및 분할 치료 간(intefractional)에 발생되는 오차를 측정하는 자동분석소프트웨어를 개발하였다. 오차 분석 결과는 3차원 입체조형 방사선치료를 임상에 적용하기에 앞서 적절한 치료계획용적(Planning Target Volume, PTV)을 설정하는 데 있어 매우 중요하다. 본 연구에서는 전자포탈영상장치(Electrical Portal imaging Device, EPID)로써 Portal Vision LC250 액체 충전형 이온화 검출기를 사용하였다(fast frame-averaging 모드, 초당 1.4 프레임, 256X256 픽셀). 12명의 환자에 대해 최소 7일 이상씩 영상을 획득하였다. 매 치료마다 평균 8 내지 9개의 영상을 각 빔에 대해 얻었다(분당 400 MU 선량률). 총 2,931 (720 측정을 포함하는)개의 영상을 정량적으로 분석할 수 있는 자동화 영상 분석 소프트웨어를 개발하였다. 이를 통해 호흡으로 인해 발생하는 치료도중 오차와 분할 치료간 발생하는 분할치료오차의 표준편차($\sigma$)들을 계산하였다. 신뢰 구간 95%로 임상표적체적(Clinical Target Volume, CTV)을 포함할 수 있는 PTV 마진은 $2\;(1.96\;{\sigma})$으로 계산되었다. 주로 호흡으로 인해 유발되는 치료도중오차를 보상하기 위해 필요한 PTV 마진은 2 mm에서 4 mm이었다. 반면에 분할 치료간 오차를 보상하기위해 필요한 PTV 마진은 7 mm에서 31 mm이었다. 12명의 환자에 대한 전체 평균오차는 17 mm이었다. 분할치료 간 오차는 호흡에 의해 유발되는 치료도중 오차에 비해 2배에서 15배까지 더 크게 나타났다. 유방암 치료에 있어 3차원 입체정형조사나 세기조절방사선치료(Intensity Modulated Radiation Therapy, IMRT)를 적용하기에 앞서 반드시 셋업 오차의 크기를 측정하여 PTV에 적절히 반영되어야 한다. 유방에 대한 3차원 입체정형조사나 세기조절방사선치료를 위해 반드시 필요한 것은 아니지만, 큰 PTV 마진을 줄여주기 위해서는 영상유도방사선치료(Image Guided Radiation Therapy, IGRT)가 매우 유용하게 이용될 수 있다. 전자포탈영상장치 들은 본 보고서에서 기술한 바와 같은 자동분석소프트웨어를 반드시 포함하여야 한다. 이를 통해 수많은 EPID 영상들을 자동화 처리하고 오차분석을 시행함으로써 각 병원의 임상적용 방법 및 환경에 따라 상이하게 나타날 수 있는 오차의 크기를 감안한 적절한 PTV마진을 구하는데 도움을 얻을 수 있다. 이러한 장치들은 또한 최소의 노력으로 환자 치료를 관찰할 수 있는 귀중한 정보를 제공해 준다.

• 대한화학회지
• /
• 제29권3호
• /
• pp.239-246
• /
• 1985

[Cu(dl-trans-[14]-diene)]$^{2+}$ 착이온과 음이온(L$^{n-}$ = S$_2$O$_3^{2-}$, SCN$^-$, I$^-$, NO$_2^-$)사이의 화학평형 관계를 온도 15 ~ 35$^{\circ}C$, 압력 1 ~ 1500bar 범위에서 분광광도법으로 연구하였다. S$_2$O$_3^{2-}$, SCN$^-$, I$^-$ 및 NO$_2^-$일때 평형상수 K는 25$^{\circ}C$, 1500bar에서 각각 3.0, 1.9, 0.6 및 0.5이었다. K값은 온도와 압력이 증가함에 따라 감소하였다. K에 대한 온도 영향으로 부터 열역학적 파라메터(${\Delta}G^{\circ}$, ${\Delta}H^{\circ}$, ${\Delta}S^{\circ}$)를 계산하였으며, 이 결과 NO$_2^-$이온을 제외한 다른 이온(S$_2$O$_3^{2-}$, SCN$^-$, I$^-$)에서는 모두 발열반응으로 나타났다. K에 대한 압력 영향에서 얻은 반응체적변화(${\Delta}$V)는 각 음이온에서 모두 양의 값이었다. S$_2$O$_3^{2-}$ 이온일 때 1,500, 1,000 및 1,500bar에서의 ${\Delta}$V(cm$^3$/mole)는 각각 26,22,19 및 16이었고, 상압에서 S$_2$O$_3^{2-}$, SCN$^-$, I$^-$ 및 NO$_2^-$이온일 때 ${\Delta}$V(cm$^3$/mole)는 각각 26, 30, 64 및 45이었다. S$_2$O$_3^{2-}$이온인 경우 Fuoss식에 의하여 계산한 평형상수와 실험치를 비교하여 Cu(Ⅱ) 착이온과 L$^{n-}$간의 결합성을 고찰하였다.

• 대한방사선기술학회지:방사선기술과학
• /
• 제35권3호
• /
• pp.265-273
• /
• 2012

선형가속기를 기반으로 하는 세기조절방사선치료와 정위적방사선수술에서는 치료계획시스템의 소조사면에 대한 신뢰할만한 선량분포를 계산하기 위해서는 우선적으로 소조사면의 정확한 빔 자료 측정이 선행되어야 한다. 특히 소조사면의 빔 자료 측정에서 조사면 가장자리에서의 급격한 선량 변화, 측면 전자비평형, 그리고 검출기의 체적 영향으로 인한 적절한 검출기 선택이 중요하다. 따라서 본 연구에서는 선형가속기의 소조사면에 대한 빔 자료 측정에 있어서 검출기의 선량 특성을 알아보고자 하였다. 검출기는 0.01 cc 부피와 0.13 cc 부피의 이온전리함과 정위적다이오드를 사용하였으며, 빔 자료는 광자선(6 MV와 15 MV)에 대하여 조사면 크기를 $2{\times}2cm^2$에서 $5{\times}5cm^2$까지 변화시켜 각 검출기를 이용하여 깊이선량백분율, 선량출력계수, 그리고 빔측면도를 측정하였다. CC01 이온전리함과 정위적다이오드 검출기를 이용한 $PDD_{20}/PDD_{10}$은 $2{\times}2cm^2$ 조사면의 경우 6 MV와 15 MV에서 각각 1.02%와 0.12% 차이를 보였다. $3{\times}3cm^2$ 이상의 조사면에서는 각 검출기를 이용하여 얻어진 $PDD_{20}/PDD_{10}$의 차이가 6 MV와 15 MV에서 각각 평균 1.15%와 0.71% 이였다. CC01 이온전리함과 정위적다이오드 검출기를 이용한 선량출력계수 측정 결과, $2{\times}2cm^2$ 조사면의 경우 6 MV와 15 MV에서 0.5%와 1.5%이내에서 일치하였다. $3{\times}3cm^2$ 이상의 조사면에서는 각 검출기의 차이가 0.5% 이내이였다. 3개의 깊이에서 측정된 빔측면도의 반음영은 정위적다이오드 검출기의 경우 6 MV와 15 MV에서 각각 평균 2.7 mm와 3.5 mm, CC01 이온전리함의 경우 각각 평균 3.4 mm와 4.3 mm, CC13 이온전리함의 경우 각각 평균 5.2 mm와 6.1 mm이였다. 이를 통해 깊이선량백분율과 선량출력계수 측정 시 $2{\times}2cm^2$ 조사면에서는 CC01 이온전리함과 정위적다이오드 검출기를 $3{\times}3cm^2$에서 $5{\times}5cm^2$ 조사면에서는 각 검출기의 사용이 가능할 것으로 판단된다. 또한 소조사면에 대한 정확한 빔측면도의 반음영을 측정하기 위해서는 유효체적이 작은 CC01 이온전리함과 정위적다이오드 검출기 사용하는 것이 타당하겠다.

• 대한방사선치료학회지
• /
• 제32권
• /
• pp.17-29
• /
• 2020

목 적: 전립선 방사선 치료 시 회피 영역을 지정한 혼합 에너지(Mixed Energy Partial Arc, MEPA) VMAT 치료 계획의 선량평가를 통해 임상 적용 가능성을 평가한다. 대상 및 방법: 전립선암 환자 9명을 대상으로 치료 장비 True Beam Stx(Varian Medical System, USA), 치료 계획시스템(Eclipse, Varian Medical System, USA)을 이용하여 치료 계획용적(PTV)에 2.5Gy, 30회로 총 선량 70Gy를 처방하여, 회피 영역을 지정한 혼합 에너지 치료 계획과 6MV 1 ARC, 6MV, 10MV 그리고 15MV 2 ARC 치료 계획과 비교 분석하였다. 혼합 에너지 치료 계획은 첫 번째 ARC에서 갠트리 각도 180°~230°, 310°~50°은 6MV로 지정하고 230°~310°, 50°~130°는 회피 영역으로 지정하였고 두 번째 ARC에서 갠트리 각도 130°~50°, 310°~230°는 에너지 15MV로 설정 후 50°~310°는 회피 영역으로 지정하였다. 2가지 치료 계획을 각각 최적화 및 선량 계산 후 Sum plan을 생성하였다. 각각의 치료 계획은 선량체적히스토그램(Dose Volume Histogram, DVH)을 이용하여 PTV의 최대선량, 평균선량, D2%, 처방 선량 지수(Conformity Index, CI), 선량 균질 지수(Homogeneity Index, HI) 및 손상위험장기(Organ at Risk, OAR)의 선량을 분석하였다. 또한 monitor unit(MU)와 치료시간(beam on time)을 분석하였다. 결 과: MEPA 치료 계획은 PTV에서 최대선량, 평균선량, D2%값은 75.86±0.47, 71.95±0.51, 73.92±0.45로 6MV 1 ARC 계획보다 낮은 값(전체 p<0.0005)을 보였고 6MV, 10MV, 15MV 2 ARC 치료 계획과의 비교에서 최대선량은 평균 차이 각각 0.24%, 0.39%, 0.60%(p<0.450, 0.321, 0.139), D2% 값 0.42%, 0.49%, 0.59%(p<0.073, 0.087, 0.033)로 높게 분석되었다. 평균선량은 10MV 2 ARC 치료 계획보다 0.09%(p<0.615) 낮았지만 6MV 2 ARC와 15MV 2 ARC와의 비교에서 각각 0.27%, 0.12%(p<0.184, 0.521) 높게 나타났다. HI는 0.064±0.006로 가장 낮은 값으로 분석(p<0.005) 되어 가장 좋은 균질도를 보이며 CI는 MEPA에서 1.12±0.038, 6MV 1 ARC, 6MV, 10MV, 15MV 2 ARC 각각 1.12±0.036, 1.11±0.024, 1.11±0.030, 1.12±0.027로 모든 치료 계획에서 유사한 값으로 분석되었다. 직장 선량은 MEPA가 모든 선량 분석 파라미터에서 낮은 값으로 분석되었다. 특히 V70Gy, V50Gy, V30Gy, V20Gy은 3.40, 16.79, 37.86, 48.09으로 직장에 조사되는 저 선량 영역의 체적을 의미 있는 만큼 줄일 수 있었다. 방광선량은 V30Gy, V20Gy이 다른 치료 계획에 비해 낮은 값을 보였다. 하지만 양측 대퇴골두 비교에서 MEPA는 평균선량 9.69±2.93, 9.88±2.5로 비교된 치료 계획에 비해 2.8Gy~3.28Gy 높은 값을 보였다. 평균 MU값은 MEPA가 6MV 1 ARC보다 19.53%, 10MV 2 ARC에 비해 5.7% 감소하였다. 결 론: 저에너지 X선(6MV)과 고에너지 X선(15MV)을 사용 한 MEPA는 각각 에너지의 장점을 취해 비교한 치료 계획보다 개선된 HI와 동등한 CI값으로 분석되어 만족할 만한 PTV coverage를 보였으며 방광과 직장을 보호하는데 효과적이라고 생각된다.

• 대한방사선기술학회지:방사선기술과학
• /
• 제35권3호
• /
• pp.255-263
• /
• 2012

본 연구에서는 AAPM TG-119 보고서를 통해 본원에 도입된 세기조절방사선치료(IMRT)와 부피적세기조절회전치료(VMAT)의 치료 전 환자별 정도관리(patient specific QA)를 시행하고자 하였다. 본원의 치료계획장치를 이용해 각각의 치료계획을 수립하여 절대선량은 표적 및 위험장기에 전리함을 이용해 측정하였고 상대선량분포는 $DELTA^{4PT}bi$-planar diode array를 사용하여 측정하였다. 치료계획의 평가는 선량체적히스토그램을 이용하였고 선량검증은 측정값과 계산값을 비교하여 시행하였다. 치료계획평가에서 전립선의 경우에는 두 치료법 모두 표적과 위험장기의 목표선량에 도달하였으며 두경부와 Multi target의 경우, 세기조절방사선치료는 표적에서는 목표선량에 도달하지 못하였지만 부피적세기조절 회전치료는 표적과 위험장기 모두 목표선량에 도달하였다. C-shape(easy)은 두 치료법이 표적과 위험장기 모두 목표선량에 도달하였고, C-shape(hard)의 경우엔, 두 치료법이 표적에서는 목표선량에 도달하였으나 위험장기에서는 목표선량에 도달하지 못하였다. 절대선량평가에서는 세기조절방사선치료의 경우 평균오차율의 평균값이 표적에서 $1.24{\pm}2.06%$, 위험장기에서 $1.4{\pm}2.9%$였고 신뢰구간은 표적에서 3.65%, 위험장기에서 4.39%였다. 부피적세기조절회전치료는 평균오차율의 평균값이 표적에서 $2.06{\pm}0.64%$, 위험장기에서 $2.21{\pm}0.74%$ 였고 신뢰구간은 표적에서 4.09%, 위험장기에서 3.04%로 두 치료법 모두 제안된 허용기준인 표적에서 4.5%, 위험장기에서 4.7% 이내였다. 상대선량평가에서는 세기조절방사선치료의 경우 허용기준을 통과하는 감마인덱스의 평균값은 $98.3{\pm}1.5%$였고 신뢰구간은 3.78%였다. 부피적세기조절회전치료의 경우엔 평균값이 $98.2{\pm}1.1%$였고 신뢰구간은 3.95%로 두 치료법 모두 제안된 허용기준인 7.0% 이내였다. 따라서 이번 연구에서 실시한 시험을 통해 본원의 세기조절방사선치료와 부피적세기조절회전치료가 보고서에서 제안된 허용기준에 모두 부합하고 사용에 적합함을 확인할 수 있었다.

• 한국토양비료학회지
• /
• 제37권3호
• /
• pp.131-135
• /
• 2004

본 연구는 일반 기상 측정에서 에너지 수지를 구할 때 이용되는 지중열류 값과 토양 내 10 cm 깊이에서 측정한 토양온도를 활용하여 태양 복사에너지와 지표면 복사 에너지의 도달과 분배에 중요한 역할을 하는 토양표면의 온도를 예측하고자 수행하였다. 2003년 6월 10일부터 6월 24일까지 측정한 지중열류와 토양온도 그래프에서 토양온도나 지중열류는 주기성을 나타내며 일중 최저 지중열류와 최고 토양온도 사이에는 위상편차가 존재한다. 토심 5 cm에서 측정하여 시간별로 평균한 토양온도와 지중열류 사이에는 2시간의 시간지연이 존재하며 토양온도와 시 공간상에서의 지중열류는 정의 상관관계를 보였다. 이는 단위체적 당 열용량과 깊이에 따른 열량의 변화율이 태양에너지와 지표면 복사를 통한 지중열류에 비례한다는 것을 의미한다. 예측된 토양 표면온도는 시간별로 평균하였을 때 그 평균온도가 $20^{\circ}C$를 넘어 여름철의 기온을 반영하였으며 모양도 주기함수 형태를 보였다. 진폭은 $4.5^{\circ}C$로서 10 cm 깊이에서의 진폭인 $3.4^{\circ}C$보다 $1.1^{\circ}C$ 높았으며 최저온도가 나타난 시간은 토양표면의 경우는 오전 8시, 10 cm 깊이에서는 오전 9시였으며 최고온도가 나타난 시간은 토양표면은 16시, 10 cm 깊이는 19시이었다. 시간별로 평균하지 않았을 때의 최고와 최저온도는 각각 33.2, $16.5^{\circ}C$ 였으며, 토양표면온도 분포는 $15-20^{\circ}C$가 5.3%, $20-25^{\circ}C$가 65.6%, $25-30^{\circ}C$가 28.1%, $30-35^{\circ}C$가 1%로 대부분의 온도는 $20-25^{\circ}C$ 범위의 값을 나타냈다. 예측된 토양표면온도의 검증을 위해 토양표면 온도와 10 cm 깊이의 토양온도를 가지고 계산한 산술 평균과 토심 5 cm에서 측정한 토양온도를 비교하였다. 또한, 그 과정을 통해 얻어진 추정회귀모형은 P값이 0.001보다 작아 유의성이 인정 되었다. 회귀모형의 결정계수는 0.968이었고 표준오차는 0.38로 예측된 토양표면온도는 추정 회귀모형에 의해 실제 값에 가깝게 추정할 수 있을 것이다.