CalTOX 모델에 의한 휘발성유기화합물의 대기 중 예측 농도와 실측 농도간의 타당성 분석에 관한 연구 A Study on Analyzing the Validity between the Predicted and Measured Concentrations of VOCs in the Atmosphere Using the CalTOX Model원문보기
Objectives: This study calculated local residents exposures to VOCs (Volatile Organic Compounds) released into the atmosphere using the CalTOX model and carried out uncertainty analysis and sensitivity analysis. The model validity was analyzed by comparing the predicted and the actual atmospheric co...
Objectives: This study calculated local residents exposures to VOCs (Volatile Organic Compounds) released into the atmosphere using the CalTOX model and carried out uncertainty analysis and sensitivity analysis. The model validity was analyzed by comparing the predicted and the actual atmospheric concentrations. Methods: Uncertainty was parsed by conducting a Monte Carlo simulation. Sensitivity was dissected with the regression (coefficients) method. The model validity was analyzed by applying r2 (coefficient of determination), RMSE (root mean square error), and the Nash-Sutcliffe EI (efficiency index) formula. Results: Among the concentrations in the atmosphere in this study, benzene was the highest and the lifetime average daily dose of benzene and the average daily dose of xylene were high. In terms of the sensitivity analysis outcome, the source term to air, exposure time, indoors resting (ETri), exposure time, outdoors at home (ETao), yearly average wind speed (v_w), contaminated area in ㎡ (Area), active breathing rate (BRa), resting breathing rate (BRr), exposure time, and active indoors (ETai) were elicited as input variables having great influence upon this model. In consequence of inspecting the validity of the model, r2 appeared to be a value close to 1 and RMSE appeared to be a value close to 0, but EI indicated unacceptable model efficiency. To supplement this value, the regression formula was derived for benzene with y=0.002+15.48x, ethylbenzene with y ≡ 0.001+57.240x, styrene with y=0.000+42.249x, toluene with y=0.004+91.588x, and xylene with y=0.000+0.007x. Conclusions: In consequence of inspecting the validity of the model, r2 appeared to be a value close to 1 and RMSE appeared to be a value close to 0, but EI indicated unacceptable model efficiency. This will be able to be used as base data for securing the accuracy and reliability of the model.
Objectives: This study calculated local residents exposures to VOCs (Volatile Organic Compounds) released into the atmosphere using the CalTOX model and carried out uncertainty analysis and sensitivity analysis. The model validity was analyzed by comparing the predicted and the actual atmospheric concentrations. Methods: Uncertainty was parsed by conducting a Monte Carlo simulation. Sensitivity was dissected with the regression (coefficients) method. The model validity was analyzed by applying r2 (coefficient of determination), RMSE (root mean square error), and the Nash-Sutcliffe EI (efficiency index) formula. Results: Among the concentrations in the atmosphere in this study, benzene was the highest and the lifetime average daily dose of benzene and the average daily dose of xylene were high. In terms of the sensitivity analysis outcome, the source term to air, exposure time, indoors resting (ETri), exposure time, outdoors at home (ETao), yearly average wind speed (v_w), contaminated area in ㎡ (Area), active breathing rate (BRa), resting breathing rate (BRr), exposure time, and active indoors (ETai) were elicited as input variables having great influence upon this model. In consequence of inspecting the validity of the model, r2 appeared to be a value close to 1 and RMSE appeared to be a value close to 0, but EI indicated unacceptable model efficiency. To supplement this value, the regression formula was derived for benzene with y=0.002+15.48x, ethylbenzene with y ≡ 0.001+57.240x, styrene with y=0.000+42.249x, toluene with y=0.004+91.588x, and xylene with y=0.000+0.007x. Conclusions: In consequence of inspecting the validity of the model, r2 appeared to be a value close to 1 and RMSE appeared to be a value close to 0, but EI indicated unacceptable model efficiency. This will be able to be used as base data for securing the accuracy and reliability of the model.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
이는 대기오염으로 인한 국민건강이나 환경에 관한 위해를 예방하고, 대기환경을 적정하고 지속가능하게 관리·보전하여 모든 국민이 건강하고 쾌적한 환경에서 생활할 수 있도록 하는 것에 목적이 있다.
이러한 결과는 모델은 효율적이지 못하지만, 예측된 값의 분산이 관측된 값의 분산과 같다는 것을 의미 한다. 이러한 단점을 극복하기 위하여 본 연구는 R2의 기반이 되는 회귀의 기울기와 절편을 도출하였다.
그리고 모델의 대기 중 농도 예측값과 대기 중 실측된 농도와의 차이를 비교하여 모델의 타당성을 분석하였다. 이로써 CalTOX 모델에 의한 인체노출량 평가 결과의 정확성과 신뢰성을 확보하는 기초자료로 활용하고자 하였다.
가설 설정
37) 그러나 우리나라의 연구자들은 모델을 이용한 인체노출량 예측 시 입력변수 정보를 일일이 찾아서 수동으로 입력하거나, 서로 다른 출처로부터 입력 자료를 선정할 때는 공학적 판단까지 해야 하는 경우도 있어 어려움이 있다.37) 이러한 이유들로 인하여 인체노출량을 예측하는 모델 시스템 자체의 접근이 쉽지 않다. 그러나 최근 국가차원에서 노출평가 적용을 위한 한국인의 노출계수 핸드북48) 등의 개선뿐만 아니라, 국가 정보 공개 제도와 그 밖의 통계자료 시스템 등을 운영하며, 빅데이터들을 생산해 내고 있다.
제안 방법
그러나 분산만 정량화 된다는 사실은 R2을 항상 체계적으로 과대평가하거나 과소평가 하는 모델은 모델 예측이 잘 못되더라도 1에 가까운 R2 값을 예측한다는 단점이 있다.26,27) 따라서 본 연구에서는 R2의 기반이 되는 회귀의 기울기와 절편을 도출하였다. RMSE는 체계적 오류와 무작위 오류를 모두 측정하며, 모델에서 예측한 값과 실제로 관측된 값의 차이를 나타낸다.
39) 본 연구는 다매체ㆍ다경로 모델인 CalTOX 모델을 이용하여 산업체에서 대기로 배출된 VOCs의 대기 중 농도와 인체노출량을 예측하였다. 모델의 예측값이 실제 측정 농도와 일치할 수는 없다.
CalTOX 모델이 요구하는 배출량 자료와 물리·화학적 특성 자료, 지역 특성 자료, 인체 노출 특성 자료 중에서 본 연구는 지역 특성과 인체노출 특성 자료의 입력변수 값은 VOCs 각 물질마다 같은 값으로 입력하였고, 배출량과 물리·화학적 특성 자료만 다르게 입력하였다.
본 연구는 CalTOX 모델을 이용하여 산업체에서 대기로 배출된 VOCs에 대한 대기 중 농도를 예측하였고, 지역 주민들의 인체노출량을 산정하였으며, 불확실성 분석과 민감도 분석을 수행하였다. 그리고 모델의 대기 중 농도 예측값과 대기 중 실측된 농도와의 차이를 비교하여 모델의 타당성을 분석하였다. 이로써 CalTOX 모델에 의한 인체노출량 평가 결과의 정확성과 신뢰성을 확보하는 기초자료로 활용하고자 하였다.
CalTOX 모델이 예측한 연구대상 지역의 대기 중 농도에 대한 불확실성을 분석한 결과는 Table 1과 같다. 대기 중 농도는 가스상과 입자상으로 예측되는데, 본 연구에서 입자상 농도 예측값은 매우 미미한 것으로 분석되어 의미가 없는 것으로 판단하고, 가스상 농도만을 연구에 활용하였다. 벤젠의 농도는 95% 백분위 수에서 90% 신뢰구간이 1.
배출량은 Log-normal 확률분포 형태, VOCs의 물리·화학적 특성은 모두 Uniform 확률분포 형태로 적용하였다.
본 연구는 CalTOX 모델을 이용하여 산업체에서 대기로 배출된 VOCs에 대한 대기 중 농도를 예측하였고, 지역 주민들의 인체노출량을 산정하였으며, 불확실성 분석과 민감도 분석을 수행하였다. 그리고 모델의 대기 중 농도 예측값과 대기 중 실측된 농도와의 차이를 비교하여 모델의 타당성을 분석하였다.
그리고 이러한 입력변수들이 모델에 자동으로 입력되어지도록 네트워크화하고, 입력변수들이 정확하게 입력되었는지 자체 점검하는 과정도 포함시켜야 한다. 본 연구에서 입력변수들은 정부와 공공기관, 기존의 연구에서 신뢰할만하다고 판단되는 자료들을 최대한 이용하였다. 그러나 입력 변수 자료가 부재하거나, 존재한다고 하더라도 입력 변수로 사용이 어렵다고 판단되는 경우에는 미국 캘리포니아 지역의 입력변수를 사용하였음으로 이러한 측면에서 연구의 한계가 존재한다.
)와 Root mean-square error(RMSE)는 만족할 만한 결과가 수치로 제시되었으나, Nash-Sutcliffe 효율 지수(EI)를 이용한 결과는 음수로 나타나서 허용되지 않는 모델 성능으로 나타났다. 본 연구에서는 모델의 대기 중 농도 예측 결과와 대기 중 실측 농도와의 관계를 회귀식으로 제시하여 모델의 타당성 분석 결과를 보완하였다. 그 결과 본 연구 대상 물질 각각의 회귀식은 각각 다르게 나타났다.
지역 특성 입력변수 중 연평균 강수량은 Expon 확률분포 형태, 연평균 기온은 Kumaraswamy 확률분포 형태, 연평균 풍속은 Weibull 확률분포 형태, 지표수에서 물의 증발은 Expon 확률분포 형태를 적용하였다. 인체노출 특성 변수에서 몸무게와 호흡률(활동, 휴식)은 Triangular 확률분포 형태를 적용하였고, 노출기간과 노출량이 평균되어지는 기간은 Uniform 확률분포 형태를 적용하였다. 이외에는 Normal과 Log-normal 확률분포 형태를 적용하였다.
인체노출량 예측의 민감도 분석에서는 각각의 물질마다 7개 변수를 영향요인으로 도출하였다. 본 연구의 대상 물질들은 배출량(source term to air)이 결정계수 0.
대상 데이터
물리·화학적 특성 자료는 California Environmental Protection Agency (Cal/EPA)의 Department of Toxic Substances Control (DTSC)38)에서 제시한 값을 사용하였고, 지역 특성 자료 중 기상자료는 기상청46)과 국내 연구자료,24) 면적은 통계청47)에서 수집하였다.
연구 대상 물질의 배출량 자료는 화학물질배출·이동량정보(pollutant release and transfer register, PRTR)23)에서 수집하였다. 대상 지역은 화학물질관리법 시행규칙 제5조와 관련한 배출량 조사대상 화학물질 중, 80개 화학물질이 대기로 배출되고 있었으며, 배출량은 494,448 kg/년이었다. 본 연구의 대상물질은 총 98,580 kg/년이 배출되어 연구대상 지역 전체 배출량의 약 20%를 차지하는 것으로 나타났다.
실측한 대기 중 농도는 국립환경과학원의 대기환경연보43)의 자료를 이용하였다.
연구 대상 물질의 배출량 자료는 화학물질배출·이동량정보(pollutant release and transfer register, PRTR)23)에서 수집하였다.
인체노출 특성은 ‘한국인의 노출계수 핸드북’48) 등을 이용하였다. 이외에 국내 자료가 없는 경우에는 DTSC에서 제시한 캘리포니아 지역 자료를 이용하였다.
충청남도 서산시를 연구대상 지역으로 하였고, 대기 중으로 VOCs (벤젠, 톨루엔, 에틸벤젠, 자일렌, 스티렌)를 배출하는 산업체는 ‘국가 대기오염물질 배출량 산정방법 편람(III)’44)을 기준으로 하여, 2017년 코크스, 화합물 및 화학제품 제조업, 석유정제품 및 핵연료 산업체로 하였다.
데이터처리
민감도 분석도 @Risk 7.6 프로그램을 이용하였고, Regression (coefficients) 방법으로 하였다. 이 방법은 입력변수를 독립변수로 놓고 결과변수를 종속변수로 하는 다중회귀분석의 표준화 회귀계수(beta coefficients, β)로써, 결정계수(RSqr, r2)로 결과를 나타내는데, r2의 값은 “0”과 “1” 사이에 존재하며, 그 값이 클수록 독립변수들이 종속변수의 변이를 잘 설명함을 의미한다.
불확실성은 @Risk 7.6 (Palisade Co., USA) 프로그램을 이용하여, 확률분포를 이용한 각 1,000번의 몬테카를로 모의실험을 통하여 분석하였다. 입력 변수 각각에 대한 확률분포 모형은 김옥 등24)의 연구를 참조하였다.
이론/모형
모델에 의해 예측된 대기 중 농도의 타당성 분석을 위해 대기 중 실측된 농도와의 차이 비교를 위한 효율성 공식을 이용하였다. 결정계수(R2), Root mean-square error (RMSE), Nash-Sutcliffe 효율 지수(EI)25) 를 적용하였으며, Eq. (3), (4), (5)와 같다.
모델에 의해 예측된 대기 중 농도의 타당성 분석을 위해 대기 중 실측된 농도와의 차이 비교를 위한 효율성 공식을 이용하였다. 결정계수(R2), Root mean-square error (RMSE), Nash-Sutcliffe 효율 지수(EI)25) 를 적용하였으며, Eq.
, USA) 프로그램을 이용하여, 확률분포를 이용한 각 1,000번의 몬테카를로 모의실험을 통하여 분석하였다. 입력 변수 각각에 대한 확률분포 모형은 김옥 등24)의 연구를 참조하였다. 배출량은 Log-normal 확률분포 형태, VOCs의 물리·화학적 특성은 모두 Uniform 확률분포 형태로 적용하였다.
배출량은 Log-normal 확률분포 형태, VOCs의 물리·화학적 특성은 모두 Uniform 확률분포 형태로 적용하였다. 지역 특성 입력변수 중 연평균 강수량은 Expon 확률분포 형태, 연평균 기온은 Kumaraswamy 확률분포 형태, 연평균 풍속은 Weibull 확률분포 형태, 지표수에서 물의 증발은 Expon 확률분포 형태를 적용하였다. 인체노출 특성 변수에서 몸무게와 호흡률(활동, 휴식)은 Triangular 확률분포 형태를 적용하였고, 노출기간과 노출량이 평균되어지는 기간은 Uniform 확률분포 형태를 적용하였다.
향후 20년의 대기 중 농도와 인체노출량을 예측하기 위하여 CalTOX™4.0 모델(University of California)을 이용하였다.
성능/효과
CalTOX 모델을 이용하여 대기로 배출된 VOCs 물질에 대한 대기 중 농도와 인체노출량를 예측하고 불확실성 분석과 민감도 분석, 모델의 타당성을 분석한 결과, 본 연구 대상 물질인 VOCs의 대기 중 농도는 벤젠이 1.08×10−4 (95% CI; 9.88×10−5~1.30×10−4 ) mg/m3로 가장 높게 예측되었고, 대기 중 농도 분포 예측에 대한 민감도를 분석한 결과는 배출량 (source term to air), 연간 평균 풍속(yearly average wind speed), 면적(contaminated area in m2) 등으로 나타났다.
그럼에도 불구하고 모델의 대기 중 농도 분포 예측의 불확실도를 낮추기 위해서는 본 연구의 민감도 분석 결과에 기초하여, 정확성과 신뢰성 있는 입력 자료를 우선 확보하고 이를 데이터베이스화할 필요가 있다고 판단된다. 대기 중 가스상 농도 예측의 민감도 분석 결과는 배출량(source term to air), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area) 등으로 나타났으며, 입자상 농도는 대기 먼지 하중(atmospheric dust load; rhob_a), 유기 에어로졸의 분획(Aerosol organic fraction; foc_ap), 배출량(source term to air) 등으로 나타났다. 또한 인체노출량 예측에 적용된 다양한 입력 변수들의 민감도 분석을 수행한 결과는, 배출량 (source term to air), 실내 휴식시간(exposure time, indoors resting; ETri), 실외 노출시간(exposure time, outdoors at home; ETao), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area), 활동 시 호흡률(active breathing rate; BRa), 휴식 시 호흡률(resting breathing rate; BRr), 실내 활동시간(exposure time, active indoors; ETai) 등이 주요한 입력변수로 도출되었다.
또한 벤젠의 전생애평균노출량(LADDsim) 1.36×10−5 (95% CI;1.17×10−5~1.71×10−5 ) mg/kg/day과 자일렌의 일일평균노출량(ADDsim)의 1.34×10−5 (95% CI; 1.24×10−5~1.74×10−5 ) mg/kg/day는 다른 물질들보다 높은 수준으로 나타났다.
대기 중 가스상 농도 예측의 민감도 분석 결과는 배출량(source term to air), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area) 등으로 나타났으며, 입자상 농도는 대기 먼지 하중(atmospheric dust load; rhob_a), 유기 에어로졸의 분획(Aerosol organic fraction; foc_ap), 배출량(source term to air) 등으로 나타났다. 또한 인체노출량 예측에 적용된 다양한 입력 변수들의 민감도 분석을 수행한 결과는, 배출량 (source term to air), 실내 휴식시간(exposure time, indoors resting; ETri), 실외 노출시간(exposure time, outdoors at home; ETao), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area), 활동 시 호흡률(active breathing rate; BRa), 휴식 시 호흡률(resting breathing rate; BRr), 실내 활동시간(exposure time, active indoors; ETai) 등이 주요한 입력변수로 도출되었다. 그 중에서 지역 입력변수인 v_w와 Area 변수를 제외한 ETri, ETao, ETai는 지역 주민들의 시간활동 양상과 관련된 입력변수들이며, BRa와 BRr은 호흡률과 관련된 입력변수로 구분할 수 있다.
7 이상으로 가장 큰 영향을 미치는 것으로 나타났고, 그 다음으로는 실내 휴식시간(exposure time, indoors resting; ETri), 실외 노출시간(exposure time, outdoors at home; ETao)의 순으로 나타났다. 벤젠은 그 다음으로 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area), 활동 시 호흡률(active breathing rate; BRa), 휴식 시 호흡률(resting breathing rate; BRr)의 순으로 나타났으며, 에틸벤젠은 그 다음으로 연간 평균 풍속(yearly average wind speed; v_w), 활동 시 호흡률(active breathing rate; BRa), 면적(contaminated area in m2; Area), 휴식 시 호흡률(resting breathing rate; BRr)로 나타났고, 스티렌은 면적(contaminated area in m2; Area), 활동 시 호흡률(active breathing rate; BRa), 휴식 시 호흡률(resting breathing rate; BRr), 실내 활동시간(exposure time, active indoors; ETai)의 순으로 나타났으며, 톨루엔과 자일렌은 연간 평균 풍속(yearly average wind speed; v_w), 활동 시 호흡률(active breathing rate; BRa), 면적(contaminated area in m2; Area), 휴식 시 호흡률(resting breathing rate; BRr)로 나타났다.
벤젠의 회귀식은 y=0.002+15.485x로 96.8%의 설명력을 갖는 것으로 나타났으며, 에틸벤젠은 y ≡ 0.001+57.240x로 94.7%의 설명력을 갖는 것으로 나타났고, 스티렌은 y=0.000+42.249x로 94.8%의 설명력을 갖는 것으로 나타났다.
8) 이와 관련한 국내의 논문에는 이보람 등,42) 박진현 등,40) 류현수 등41)의 연구가 있다. 본 연구 결과와 함께 이러한 결과들을 통하여 VOCs의 인체노출량은, 대기 농도와 노출시간 또는 접촉률(강도)에 따라 달라질 수 있음을 다시 한 번 확인할 수 있었다.
인체노출량 예측의 민감도 분석에서는 각각의 물질마다 7개 변수를 영향요인으로 도출하였다. 본 연구의 대상 물질들은 배출량(source term to air)이 결정계수 0.7 이상으로 가장 큰 영향을 미치는 것으로 나타났고, 그 다음으로는 실내 휴식시간(exposure time, indoors resting; ETri), 실외 노출시간(exposure time, outdoors at home; ETao)의 순으로 나타났다. 벤젠은 그 다음으로 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area), 활동 시 호흡률(active breathing rate; BRa), 휴식 시 호흡률(resting breathing rate; BRr)의 순으로 나타났으며, 에틸벤젠은 그 다음으로 연간 평균 풍속(yearly average wind speed; v_w), 활동 시 호흡률(active breathing rate; BRa), 면적(contaminated area in m2; Area), 휴식 시 호흡률(resting breathing rate; BRr)로 나타났고, 스티렌은 면적(contaminated area in m2; Area), 활동 시 호흡률(active breathing rate; BRa), 휴식 시 호흡률(resting breathing rate; BRr), 실내 활동시간(exposure time, active indoors; ETai)의 순으로 나타났으며, 톨루엔과 자일렌은 연간 평균 풍속(yearly average wind speed; v_w), 활동 시 호흡률(active breathing rate; BRa), 면적(contaminated area in m2; Area), 휴식 시 호흡률(resting breathing rate; BRr)로 나타났다.
대상 지역은 화학물질관리법 시행규칙 제5조와 관련한 배출량 조사대상 화학물질 중, 80개 화학물질이 대기로 배출되고 있었으며, 배출량은 494,448 kg/년이었다. 본 연구의 대상물질은 총 98,580 kg/년이 배출되어 연구대상 지역 전체 배출량의 약 20%를 차지하는 것으로 나타났다.
모델에 의해 예측된 연구대상 물질들의 대기 중 농도 분포와 인체노출량에 유의한 영향을 미치는 변수들을 도출한 민감도 분석 결과는 Table 3, 4와 같다. 연구 대상 물질들은 대기 중에 가스상과 입자상 농도 분포로 예측되며, 가스상 농도의 민감도 분석 결과는 배출량(source term to air), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area)으로 나타났으며, 입자상 농도는 대기 먼지 하중(atmospheric dust load; rhob_a), 유기 에어로졸의 분획(Aerosol organic fraction; foc_ap), 배출량(source term to air), 대기 온도(Ambient environmental temperature; Temp), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area)으로 나타났다.
모델의 대기 중 농도 예측 결과의 타당성을 검증한 결과는 Table 5와 같다. 연구 대상 물질들은 모두 R2이 1에 가까운 값으로 나타났고, RMSE는 0에 가까운 값으로 나타났으나, EI는 음수로 나타났다. 이러한 결과는 모델은 효율적이지 못하지만, 예측된 값의 분산이 관측된 값의 분산과 같다는 것을 의미 한다.
74×10−5 ) mg/kg/day는 다른 물질들보다 높은 수준으로 나타났다. 인체노출량 예측에 대한 민감도를 분석한 결과는 배출량(source term to air), 실내 휴식시간(exposure time, indoors resting), 실외 노출시간(exposure time, outdoors at home), 연간 평균 풍속(yearly average wind speed), 면적(contaminated area in m2 ), 활동 시 호흡률(active breathing rate) 등으로 나타났다.
30×10−4 ) mg/m3로 가장 높게 예측되었고, 대기 중 농도 분포 예측에 대한 민감도를 분석한 결과는 배출량 (source term to air), 연간 평균 풍속(yearly average wind speed), 면적(contaminated area in m2) 등으로 나타났다. 타당성을 분석한 결과는 R2은 1, RMSE는 0에 가까운 값으로 나타났음으로 만족할 만한 검증 결과가 나타났으나, EI는 음수로 나타나서 허용되지 않는 모델 성능으로 판단하였다. 또한 벤젠의 전생애평균노출량(LADDsim) 1.
한편 모델 효율성 공식을 이용하여 타당성을 분석한 결과, 결정계수(R2)와 Root mean-square error(RMSE)는 만족할 만한 결과가 수치로 제시되었으나, Nash-Sutcliffe 효율 지수(EI)를 이용한 결과는 음수로 나타나서 허용되지 않는 모델 성능으로 나타났다. 본 연구에서는 모델의 대기 중 농도 예측 결과와 대기 중 실측 농도와의 관계를 회귀식으로 제시하여 모델의 타당성 분석 결과를 보완하였다.
후속연구
본 연구에서 입력변수들은 정부와 공공기관, 기존의 연구에서 신뢰할만하다고 판단되는 자료들을 최대한 이용하였다. 그러나 입력 변수 자료가 부재하거나, 존재한다고 하더라도 입력 변수로 사용이 어렵다고 판단되는 경우에는 미국 캘리포니아 지역의 입력변수를 사용하였음으로 이러한 측면에서 연구의 한계가 존재한다.
그 중에서 지역 입력변수인 v_w와 Area 변수를 제외한 ETri, ETao, ETai는 지역 주민들의 시간활동 양상과 관련된 입력변수들이며, BRa와 BRr은 호흡률과 관련된 입력변수로 구분할 수 있다. 그러므로 본 연구와 같이 모델을 이용한 인체노출량을 예측할 때 그 정확도를 향상시키기 위해서는 배출량과 지역주민들의 시간활동 양상, 호흡률에 대한 명확한 데이터 검증이 선행되어야 할 것으로 판단한다. WHO에서는 일반 및 대상 집단의 시간 활동 패턴을 정확하게 안정적으로 문서화하는 것은 인체 노출을 이해하고 완화하는 데 중요한 요소라고 지적하였다.
모델의 예측값이 실제 측정 농도와 일치할 수는 없다. 그럼에도 불구하고 모델의 대기 중 농도 분포 예측의 불확실도를 낮추기 위해서는 본 연구의 민감도 분석 결과에 기초하여, 정확성과 신뢰성 있는 입력 자료를 우선 확보하고 이를 데이터베이스화할 필요가 있다고 판단된다. 대기 중 가스상 농도 예측의 민감도 분석 결과는 배출량(source term to air), 연간 평균 풍속(yearly average wind speed; v_w), 면적(contaminated area in m2; Area) 등으로 나타났으며, 입자상 농도는 대기 먼지 하중(atmospheric dust load; rhob_a), 유기 에어로졸의 분획(Aerosol organic fraction; foc_ap), 배출량(source term to air) 등으로 나타났다.
본 연구의 이러한 결과들은 CalTOX 모델에 의한 인체노출량 평가 결과의 정확성과 신뢰성을 확보하는 기초자료로 활용될 것으로 예상한다.
따라서 배출량과 물리·화학적 특성의 차이가 회귀식의 도출에 영향을 미치는 것으로 추정할 수 있다. 앞으로 이러한 결과에 대한 연구가 진행되어야 할 것으로 생각한다.
그러나 최근 국가차원에서 노출평가 적용을 위한 한국인의 노출계수 핸드북48) 등의 개선뿐만 아니라, 국가 정보 공개 제도와 그 밖의 통계자료 시스템 등을 운영하며, 빅데이터들을 생산해 내고 있다. 앞으로도 국가 차원에서 모델의 입력 변수가 존재한다면 신뢰할 만한 자료로 보완하고, 존재하지 않는다면 이를 개발해 나가는 노력을 좀 더 기울여야 할 것이다. 그리고 이러한 입력변수들이 모델에 자동으로 입력되어지도록 네트워크화하고, 입력변수들이 정확하게 입력되었는지 자체 점검하는 과정도 포함시켜야 한다.
질의응답
핵심어
질문
논문에서 추출한 답변
휘발성유기화합물은 저농도로 장기간 인체에 노출되면 어떤 피해를 입는가?
그러나 이와 관련된 화학사고가 잇따르고, 이 물질들이 인체에 위험을 초래한다는 것이 확실해지면서1-4) 쾌적한 삶의 질에 대한 국민적 요구와, 위해(危害)지역 주민들의 건강에 대한 우려가 고조되었다.5) 특히 이 물질들 중 산업체에서 많이 사용하고 있는 휘발성유기화합물(volatile organic compounds; VOCs)은 저농도로 장기간 인체에 노출되면, 암을 유발하기도 하고 신경독성 질환 등이 나타날 수 있어 매우 위해하다.6)
국내에서는 VOCs 물질을 어떻게 지정고시하고 있는가?
우리나라는 대기환경보전법 제2조 제10호에 근거하여 배출시설(시행령 제45조제1항)의 관리대상 VOCs 물질 37종을 지정고시하였다. 그 중에서 톨루엔과 자일렌(o-, m-, p-포함)은 오존전구물질로 관리하고, 스티렌은 특정대기유해물질로 관리하며, 벤젠과 에틸벤젠은 오존전구물질과 특정대기유해물질로 관리하고 있다.
휘발성유기화합물의 인체노출평가에서 노출과 노출평가의 의미는 무엇인가?
이러한 VOCs의 인체노출평가에서 노출은 사람과 하나 이상의 VOCs 물질 사이의 시간과 공간에 대한 접촉을 의미하고,7,8) 노출평가는 이 물질에 대한 개인이나 집단의 환경에서 노출강도와 빈도, 지속 시간 등을 예측하거나 측정하는 과정을 의미한다.9) 화학물질의 저농도 만성노출로 인한 인체위해를 평가하여 정책 결정자의 의사 결정을 지원하고, 그 근거를 제공하기 위해 컴퓨터 시뮬레이션 모델을 많이 사용하고 있다.
참고문헌 (50)
Korea Industrial Complex Corporation. Available: http://www.kicox.or.kr [accessed 5 October 2018].
Lee BJ, Heo J, Jung DY, Kim SS, Ryu HS, Shuai JF, et al. Correlation relationship between personal exposure and biological monitoring for airborne toluene in an industrial complex and general environments. J Environ Health Sci. 2017; 43(4): 324-333.
Lopez-Aparicio S, Guevara M, Thunis P, Cuvelier K, Tarrason L. Assessment of discrepancies between bottom-up and regional emission inventories in norwegian urban areas. Atmos Environ. 2017; 154: 285-296.
Phetxumphou K, Dietrich AM, Shanaiah N, Smiley E, Gallagher DL. Subtleties of human exposure and response to chemical mixtures from spills. Environ Pollut. 2016; 214: 618-626.
Kim HJ, Im JY, Yun JH, Lee JH, Jeon JH, Lee CS. A Study on the characteristics of chemicals in major industrial complexes. J Environ Health Sci. 2018; 44(6): 515-523.
Park JH, Yoon CS. Model algorithms for estimates of inhalation exposure and comparison between exposure estimates from each model. J Korean Soc Occup Environ Hyg. 2019; 29(3): 358-367.
US NRC. Frontiers in assessing human exposures to environmental toxicants. Washington, D.C.: National Academy Press; 1991.
WHO. 2000. Available: http://www.inchem.org/documents/ehc/ehc/ehc214.htm [accessed 5 June 2018].
WHO. Harmonization project strategic plan: harmonization of approaches to the assessment of risk from exposure to chemicals. 2012.
Robert G. Sargent. Verification and validation of simulation models. Proceedings of the 30th Conference on Winter Simulation Conference. 1998.
Thomas E. McKone, Matthew MacLeod. Tracking multiple pathways of human exposure to persistent multimedia pollutants: regional, continental, and global-scale models. Annu. Rev. Environ. Resour. 2003; 28: 463-492.
US NRC. Models in environmental regulatory decision making. Washington, DC. The National Academies Press; 2007.
US EPA. Guidance on the development, evaluation, and application of environmental models. Office of the Science Advisor. Council for Regulatory Environmental Modeling. 2009.
US EPA. Guidelines for human exposure assessment. Washington, DC. Risk Assessment Forum. 1992.
Schlesinger S. Terminology for model credibility. Simulation. 1979; 32(3): 103-104.
US EPA. Framework for human health risk assessment to inform decision making. Office of the Science Advisor. Risk Assessment Forum. 2014.
Beven JK. Rainfall-runoff modelling-the primer. John Wiley & Sons Ltd., Chichester. 2001.
Yang JY. Estimation of human exposure to dioxins in korean urban residents by mulitimedia/multiroute model. [Seoul]: Yonsei University; 2001.
Moon JY, Yang JY, Lim YW, Park SE. Estimating human exposure to benzo(a)pyrene through multimedia/ multiroute exposure scenario. J. Environ. Toxicol. 2003; 18(4): 255-269.
Kim O. A multimedia & multiple pathway estimating human exposure of benzo(a)pyrene by using CalTOX model, focused on the fuel combustion in Seosan area. [Gongju]: Kongju National University; 2017.
Kim O, Lee MW, Song YH, Choi JH, Park SH, Park CY, et al. Uncertainty and sensitivity analyses of human aggregate risk assessment of benzene using the CalTOX Model. J Environ Health Sci. 2020; 46(2): 136-149.
Song YH. Assessment of integrated human risk of heavy metals emitted from coal-fired power plant using 3MRA model, [Gongju]: Kongju National University; 2018.
Ministry of environment. Available: http://icis.me.go.kr/prtr/prtrInfo/entrpsSearch.do [accessed 5 May 2020].
Kim O, Song YH, Choi JH, Park SH, Park CY, Lee MW, et al. Human exposure to BTEX and its risk assessment using the CalTOX model according to the probability density function in meteorological input data. J Environ Health Sci. 2019; 45(5): 497-510.
Nash JE, Sutcliffe JV. River flow forecasting through conceptual models part I-A discussion of principles. Journal of Hydrology. 1970; 10(3): 282-290.
Xiaohui Zhong, Utpal Dutta. Engaging Nash-Sutcliffe efficiency and model efficiency factor indicators in selecting and validating effective light rail system operation and maintenance cost models. Journal of Traffic and Transportation Engineering. 2015; 3: 255-265.
Krause P, Boyle DP, Base F. Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci. 2005; 5: 89-97.
Legates DR, McCabe GJ. Evaluating the use of ‘Goodness-of-Fit’ measures in hydrologic and hydroclimatic model validation. Water Resources Research. 1999; 35(1): 233-241.
Moriasi DN, Arnold JG, Van Liew MW, Binger RL, Harmel RD, Veith TL. Model evaluation guidelines for systematic qualification of accuracy in watershed simulations. American Society of Agricultural and Biological engineers. 2007; 50(3): 885-900.
Servat E, Dezetter A. Selection of calibration objective functions in the context of rainfall-runoff modeling in a sudanese savannah area. Hydrological Science Journal. 1991; 36(4): 307-330.
US EPA. Risk assessment guidance for superfund volume I: human health evaluation manual (Part F, supplemental guidance for inhalation risk assessment). Washingon, D.C., Office of Superfund Remediation and Technology Innovation. 2009.
Delmaar JE, Bremmer HJ. The ConsExpo spray model-Modelling and experimental validation of the inhalation exposure of consumers to aerosols from spray cans and trigger sprays. RIVM. 2009.
Delmaar J, Bokkers B, ter Burg W, Van Engelen JGM. First tier modeling of consumer dermal exposure to substances in consumer articles under reach: A quantitative evaluation of the ECETOC TRA for consumers tool. Regul Toxicol Pharm. 2013; 65(1): 79-86.
Riedmann R, Gasic B, Vernez D. Sensitivity analysis, dominant factors, and robustness of the ECETOC TRA v3, Stoffenmanager 4.5, and ART 1.5 occupational exposure models. Risk Anal. 2015; 35(2): 211-225.
Park JH, Yoon CH. Model algorithms for estimates of inhalation exposure and comparison between exposure estimates from each model. J Korean Soc Occup Environ Hyg. 2019; 29(3): 358-367.
McKone TE. CalTOX, A Multimedia Total Exposure Model for Hazardous-Waste Sites Part I: Executive Summary. California Environmental Protection Agency DTSC. 1993.
Davis PA, Olague NE, Goodrich MT. Approaches for the validation of models used for performance assessment of high-level nuclear waste repositories. 1991.
Park JH, Yang SY, Park YK, Ryu HS, Kim EC, Choe YT, et al. Exposure and risk assessment of benzene and PM10 for sub-populations using monte-carlo simulations. J Environ Health Sci. 2019; 45(3): 247-257.
Ryu HS, Yoon HJ, Eom IC, Park JH, Kim SS, Cho MS, et al. Time-activity pattern assessment for korean students. J Environ Health Sci. 2018; 44(2): 143-152.
Lee BR, Ban HK, Jang YL, Lee KY. Measurement of time/location for personal exposure assessment of air pollutants. J Environ Health Sci. 2016; 42(5): 314-323.
National Institute of Environmetal Research. Annual Report of Quality in Korea 2018. Incheon. National Institute of Environmetal Research. 2019.
National Institute of Environmental Research. National air pollutant emission calculation method manual (III). Incheon. National Institute of Environmental Research. 2013.
Hwang MY, Hong SY, Kwon YM, Jo HJ, Park CH. Differences of chemical exposure levels according to residential and personal life-style characteristics of korean adult population-from korean national environmental health survey. J Environ Health Sci. 2019; 45(2): 142-153.
Korea Meteorological Administration. ASOS. Available: https://data.kma.go.kr/cmmn/main.do[accessed 1 April 2019].
Korea Ststistical Information Service. Available: http://kosis.kr/index/index.do[accessed 1 December 2019].
National Institute of Environmental Research. Korean exposure factors handbook. Incheon. 2019.
McKone TE. CalTOX, A multimedia total-exposure model for hazardous-waste sites Part III: The multiple-pathway exposure model. California Environmental Protection Agency DTSC. 1992.
National Institute of Environmental Research. Korean exposure factors handbook for children. 2019.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.