• Journal of Astronomy and Space Sciences
• /
• 제22권4호
• /
• pp.431-440
• /
• 2005

K 지수는 특정 지역에서의 지자기 변동을 나타내는 지수로서 3시간마다 결정된다. K 지수를 결정하는 알고리즘은 태양복사의 규칙적인 일변화(Sq)를 결정하는 방법에 따라 다양하다. 그 중에서 FMI 알고리즘은 K 지수를 구하고자 하는 날과 그 전 후의 1일 자료를 포함한 총 3일 동안의 지자기 수평성분(H, D) 측정 자료를 사용하여 태양복사의 규칙적인 일변화(Sq)를 구한다. 그러나 이 알고리즘은 K 지수를 구하고자 하는 날로부터 1일 후까지의 측정 자료가 필요하기 때문에 실시간으로 K 지수를 구하는 경우에는 어려움이 있다. 본 연구에서는 준실시간으로 K 지수를 구하기 위해 미국 우주환경예보센터(NOAA/SEC)에서 제공하는 과거의 전지구적 지자기 K지수(Planetary K index: Kp)를 검색하여 최근의 Kp지수가 모두 3 이하인 지자기 활동이 조용한 날을 찾고 그 날의 지자기측정 자료를 FMI 알고리즘의 3일치 입력 자료로 사용함으로써 준실시간 K 지수를 구하는 방법을 고안하였다. 이러한 방법을 검증하기 위해 우리는 2003년 1년 동안 한국지질자원연구원에서 운용하는 경주 지자기 관측소의 FLARE+ 시스템으로 측정된 자료를 사용하고 우리의 준실시간 결정 방법을 적용하여 K 지수를 구하고 이를 일반적인 FMI 알고리즘으로 구한 K지수와 서로 비교하였다. 그 결과, 우리의 방법으로 결정된 K 지수는 일반적인 방법으로 구한 K 지수와 잘 일치하였다. 또한 준실시간 결정 방법으로 구한 K 지수를 일본 Kakioka 관측소에서 결정한 K 지수와 비교한 결과 그 상관지수(R)가 0.81 정도로 높음을 확인하였다. 본 연구의 결과는 향후 한반도 지역에서의 K 지수 변동을 준실시간으로 모니터링 하는데 활용될 수 있을 것으로 기대된다.

• 한국지구과학회지
• /
• 제36권7호
• /
• pp.649-660
• /
• 2015

2012년부터 2014년까지 한반도 및 그 주변에서 일어난 지진과 지자기장 변동성과의 상관관계를 연구하였다. 웨이블릿 기반의 셈블런스 분석기술을 지진과 관련된 지자기장 변동성 분석에 사용하였다. 총자기장을 이용한 분석 결과 관측소 반경 100 km 이내에서 발생한 지진의 경우 일관적인 셈블런스 변동양상이 나타남을 확인하였고, Z축 성분 지자기장 자료를 이용한 웨이블릿 기반의 셈블런스 분석에서도 비슷한 변동양상을 확인하였다. 자료 질이 좋은 청양 관측소 지자기장 자료를 중심적으로 분석하였고, 추가적으로 보현산 관측소의 Z축 성분 자료를 해석하여 상관성에 대한 신뢰도를 확인하고자 하였다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2007년도 공동학술대회 논문집
• /
• pp.353-360
• /
• 2007

KIGAM has installed a FLARE+ continuous geomagnetic observation system in 1998 in Daejeon of which the IAGA identification code is DZN. The coordinates of the absolute measurement plinth precisely determined by the PDGPS(Post-Processing Differential Global Positioning System) is (127-21-37.19E, 36-22-43.96N, 45.93 m) in WGS84 for the horizontal and from the geoid surface for the vertical. Periodically we have conducted the absolute geomagnetic measurement on the plinth above. We have processed the continuous time-variation 3-axis geomagnetic data measured on arbitrary sensors' coordinates in the observatory and absolute geomagnetic data together to get as the results the time-variation H(orizontal), D(eclination), Z(vertical down), F(scalar calculated from 3 components) and P(Proton Precession Magnetometer Data). We have compared our own data with those calculated from the 10th generation IGRF(International Geomagnetic Reference Field). All the measured data in the DZN Observatory can be acquired through the website http://geomag.kigam.re.kr.

• Journal of Astronomy and Space Sciences
• /
• 제22권4호
• /
• pp.419-430
• /
• 2005

이 연구에서는 지자기 폭풍의 주상 기간 동안 발생하는 태양풍 동압력 펄스에 대한 통계적 분석을 하였다. 이를 위해 먼저 1997년부터 2001년까지 5년간의 기간으로부터 지자기 폭풍 지수인 Dst 값이 -50nT 이하인 지자기 폭풍을 모두 111개 얻었다. 이러한 지자기 폭풍의 주상기간 동안에 발생한 태양풍 동압력 펄스를 정확히 조사하기 위해 태양풍 자료 뿐만 아니라 지구 저위도 여러 관측소에서 관측된 지자기 수평 성분 H값을 이용하였다. 즉 동압력 펄스가 자기권에 충돌하면 저위도 H 값이 전 지구적인 증가를 보여야 한다는 사실을 이용하였다. 이러한 과정을 통해 얻은 통계적 결과는 다음과 같다. 첫째, 자기 폭풍 중에 발생하는 H의 증가는 평균적으로 그 크기가 자기 폭풍의 강도와 비례하는 경향을 보인다. 이는 강한 자기폭풍일 수록 강한 태양풍 펄스를 동반한다는 것이다. 둘째로 자기폭풍 중에 발생하는 동압력 펄스의 발생 빈도 역시 자기 폭풍의 강도와 비례한다. 셋째, 동압력 펄스 발생 빈도가 0.4회/hr 이상인, 즉 2.5시간에 1회 이상의 동압력 펄스를 동반하는, 지자기 폭풍은 여기서 다루어진 전체 지자기 폭풍 중 약 $30\%$를 차지한다. 2.5시간은 서브스톰의 평균 지속 시간으로 볼 수 있으며, 따라서 자기 폭풍중에 서브스톰이 연속적으로 발생하는 것 만큼 자주 동압력 펄스가 나타나는 자기폭풍이 전체의 $30\%$라는 것이다. 한편 이러한 동압력 펄스의 기원을 이해하기 위해 먼저 지자기 폭풍 유도체에 대해 조사하였다. 그 결과 여기서 다루어진 지자기 폭풍의 약 $65\%$가 CME(Coronal Mass Ejection)에 의해 발생되었고 CIR(Corotating Interaction Regions)과 Type II bursts에 의해 발생한 것이 각각 6.3, $7.2\%$인 것으로 나타났다. 그런데 CME에 의해 발생된 지자기폭풍 중에서 $70\%$ 이상이 그 주상 기간이 CME와 충격파 사이의 공간인 sheath 영역 혹은 CME 앞부분에 해당되는 것으로 나타났다. 따라서 이들 지자기폭풍 주상기간에 빈번히 발생하는 동압력 펄스는 CME와 충격파 사이의 sheath 영역, 그리고 CME 앞부분 영역에서의 빈번한 태양풍 밀도 증가에 기인하는 것으로 보인다.

• 한국항해항만학회지
• /
• 제35권6호
• /
• pp.477-482
• /
• 2011

태양 흑점수의 증감주기 (약 11년)에 따른 태양폭발 (태양에서의 플레어 현상)은 태양 코로나 물질을 대방출하는 태양폭풍을 야기한다. 미국해양대기청 (NOAA: National Oceanic and Atmospheric Administration)은 태양 흑점활동이 2013년과 2014년 사이에 극대화 될 것이라고 예상했다. 강력한 태양폭풍의 영향이 지구에 미쳤을 경우 인공위성을 이용한 전 세계 측위시스템의 교란, 각종 통신수단 및 TV, 라디오 방송 등이 영향을 받을 것으로 예상된다. 실제로 1989년 태양폭풍은 캐나다에서 정전사태를 일으켜 9시간동안 약 600만 명이 정전으로 인한 피해를 입은 사례가 있다. 이와 같은 초강력 태양폭풍은 인공위성의 수명을 약 5~10년 정도 단축시키며 이로 인한 경제적 손실 및 파급효과를 고려하면 액수는 수십조 원에 달할 것으로 예상된다. 최근 2011년 2월 15일 10시 45분경 (01:30 - UTC)에 발생했던 X급 태양폭발에 의해 발생한 태양폭풍의 영향이 2011년 2월 18일 오전 10시 30분경 우리나라 (보현산 관측소)에서 관측되었다. 본 논문에서는 현재 흑점수가 증가하고 있는 시점에서 2월 18일의 태양폭발 일주일 전후 지자기 데이터를 비교하고, 또한 대전과 서울지역에서 관측한 RINEX 데이터를 이용하여 측위결과를 비교 분석하였다. 태양폭풍이 지구에 도달한 2011년 2월 18일의 지자기 관측값은 일주일 전후 데이터와 비교하여 양자(Proton) 자력계 관측결과가 요동하였고, 대전과 서울지역에서의 측위결과도 태양폭풍 일주일 전후와 비교하여 2월 18일에 가장 큰 측위오차를 보였다.

• 천문학회보
• /
• 제40권1호
• /
• pp.50.3-50.3
• /
• 2015

서울전파천문대(Seoul Radio Astronomy Observatory; SRAO)는 서울대학교 물리 천문학부(천문 전공)이 교내에 설치 운영하고 있는 전파천문대로서, 지름 6미터의 안테나와 밀리미터파 수신기를 갖추고 있다. 전파망원경은 1999년 10월에 설치를 시작한 후 1년 반만인 2001년 3월 8일에 완공되어 '최초의 전파'를 검출하였다. 안테나 조립에서부터 수신기 개발까지 대부분의 일이 대학원생들에 의해 이루어졌다는 점에서 유일무이한 전파망원경이다. 2002년 4월 2일에는 관측소를 완공하고 SRAO의 개관식을 가졌다. 설치 초기에는 85-115GHz 수신용 SIS 수신기와 최대 대역폭 50MHz의 1024채널 자기상관 분광기 등을 갖추었으며, 이후 관련 기술 개발을 지속적으로 추진하여 2002년 홀로그래피 방법을 이용해서 경면 정밀도를 90마이크론 수준으로 향상시켰으며, 2003년에는 230GHz 대역용 수신기 개발에 착수하고 2008년 완성하여 우리나라에서 처음으로 230GHz대역 전파 창을 개척하였다. 운영 초기부터 공동 활용을 도모하여 서울대학교의 학생, 교수뿐만 아니라 국내외의 다른 대학, 연구소의 연구자들에게도 시간을 할애하였다. 지난 10여년간 학술지에 약 25편의 논문을 게재하였고, 관측과 기기개발을 주제로 약 12명의 석박사를 배출하였다. 서울대학교 전파망원경은 국내 대학이 보유한 연구용 장비로는 그동안 가장 경쟁력이 있는 장비로 자리매김을 하였고, 천문기기분야의 후진을 양성하는데 중요한 기여를 하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.1214-1218
• /
• 2007

댐 하류로 탁수를 선택적으로 배제하기 위해서는 방류 탁수가 하류에 미치는 영향을 정확히 예측할 수 있는 하천 탁도 예측 및 관리시스템 구축이 필요하다. 낙동강과 반변천의 합류부에서의 이차원적인 혼합에 관한 수치해석 결과는 완전혼합을 가정하는 일차원 수질모델링의 초기 입력자료에 사용됨으로써 낙동강 본류 전체구간의 탁도 모의결과의 정확성을 높이는 데 사용될 수 있다. 본 연구는 낙동강의 중상류에 위치한 반변천 합류부에 평면 이차원 비정상 수치모형인 KU-RLMS 모형을 적용하여 탁도 변화 특성을 규명할 목적으로 수행하였다. KU-RLMS 모형은 하천 및 저수지의 국부적인 수리, 수질, 유사이동 해석을 위해 개발된 평면 이차원 비정상 수치모형이다. 직사각형 격자를 사용하는 유한차분법의 단점을 보완하기 위해, 수심적분된 2차원 연속방정식, 운동량방정식, 이송확산방정식을 불규칙한 경계를 현실적으로 모사할 수 있는 직교곡선 좌표계로 변환한 방정식을 사용한다. 이 모형은 흐름, 농도, 지형변화를 조합하여 계산할 수 있는 모형으로서 점착성 및 비점착성 유사의 이동을 모의할 수 있다. 수치모형 적용을 위한 현황분석으로 안동 및 임하 조정지댐의 방류량, 안동 수위관측소의 수위, 법흥교 및 포진교 지점의 탁도 자료를 분석하였다. 이송확산모형의 보정을 위해, 안동대교 지점의 탁도 횡분포 측정 자료를 사용하여 확산계수에 대한 매개변수 추정 및 검증을 수행하였다. 또한, 안동조정지댐과 임하조정지댐의 방류량 및 방류탁도을 고려하여 수치모의조건을 결정하였으며, 각 조건에 대한 탁도 변화 특성을 분석하였다.된 주변국이 될 수밖에 없을 것이다. 21세기 문화산업에서 우리가 판단하게 될 디자인의 가치는 계몽의 원리에 대한 '역사성'과 '현재성'의 변증법에 달려있는 것이며, 새로운 철학, 새로운 문명, 새로운 세계를 열어가는 것이다.r$ (地理志)에는 추현리와 이미 외리를 언급하면서 상주의 자기제작의 위상을 짐작하는 기록이 언급되면서 전국의 상품의 절반을 담당하고 있었음을 알 수 있었다. $\ulcorner$경상도지리지$\lrcorner$(慶尙道地理志)에는 상주가 8곳으로 1/3의 자기 생산을 담당하고 있었다. $\ulcorner$경상도지리지$\lrcorner$(慶尙道地理志)에는 $\ulcorner$세종실록$\lrcorner$(世宗實錄) $\ulcorner$지리지$\lrcorner$(地理志)와 동년대에 동일한 목적으로 찬술되었음을 알 수 있다. $\ulcorner$경상도실록지리지$\lrcorner$(慶尙道實錄地理志)에는 $\ulcorner$세종실록$\lrcorner$(世宗實錄) $\ulcorner$지리지$\lrcorner$(地理志)와의 비교를 해보면 상 중 하품의 통합 9개소가 삭제되어 있고, $\ulcorner$동국여지승람$\lrcorner$(東國與地勝覽) 에서는 자기소와 도기

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2023년도 학술발표회
• /
• pp.316-316
• /
• 2023

2022년 8월 수도권 이상폭우로 인해 서울 도심지역의 지하시설, 도로, 주택 등에 침수가 발생하면서 인명 및 재산피해가 발생하였으며, 특히 동서로 가로지르는 정체전선으로 좁고 긴 비구름이 집중되면서 국지적으로 피해가 집중되었다. 서울시의 경우 도시화에 따른 불투수지역 증가 및 내수배제 불량에 따른 빗물 역류로 인한 피해가 지속적으로 발생하고 있으며, 최근에는 기후변화에 따른 방재성능목표 강우량을 초과하는 빈도의 이상폭우로 인해 하천범람과 내수배제 불량에 따른 복합적인 원인으로 침수피해 가중되고 있는 실정이다. 또한 서울시의 경우 전체 자연적, 사회적, 경제적, 환경적 요인 등의 지역적 편차가 매우 큰 도시로 지형적인 특성뿐만 아니라 취약시설(병원, 학교 등), 수방시설물(하천, 배수시설, 빗물펌프장 등) 및 방재시설(대피소, 구호소 등) 밀도 등에 따른 침수 취약성 및 위험성 등의 편차가 매우 크기 때문에 지역특성에 대한 피해사례가 다원화 되고 있는 실정이다. 본 연구에서는 30년 이상의 종관기상관측(ASOS)과 서울시 자치구별 20년 이상의 방재기상관측(AWS)자료를 기반으로 CMIP6 SSP(Shared Socioeconomic Pathways, 공통사회 경제경로)시나리오에 따른 극한기후 지수(강수강도, 호우일수, 지속기간, 1일 최대강수량, 95퍼센타일 강수일수 등)에 대한 재현성을 평가하고 공간자기상관분석 등 시공간적인 강우특성에 대한 변화를 전망하였다. 특히 여름철 강우의 경우 자치구별 편차가 크게 나타났고 이를 통해 대도시의 도심지역의 경우 세분화하여 지역의 정확한 강우특성을 파악하는 것이 필요하다는 것을 확인할 수 있었다. 본 연구의 결과는 도심지의 방재성능 초과강우 정의와 기준을 수립하고, 장기적인 수자원 및 도시계획 차원의 대책을 마련하는데 활용될 수 있을 것으로 판단된다. 기후위기에 따른 기록적인 호우(지역별 방재성능을 초과하는 강우)에 따른 재해는 구조적인 대책을 통해 모두 저감할 수 없는 한계가 있다. 하지만 인명피해를 최소화하는 것을 목표로 기후위기에 대한 적응단계로 인식하고 수리·수문학적, 사회경제학적 등 지역특성에 따른 방재성능목표 강우량에 대한 재검토와 더불어 법제도(풍수해보험, 저류조설치 의무화 등), 개인별 재해예방, 취약계층 안전망 확보, 반지하주택 침수안전대책, 재해지도 개선 등 구조적/비구조적인 대책을 통합 수립 및 보완하는 것이 필요한 시점이다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.1115-1119
• /
• 2007

우리나라는 몬순기후대에 속하여 매년 홍수가 발생하고 이시기에 저수지로 유입된 탁수로 말미암아 저수지 인근 및 하류하천에 거주하는 주민들의 민원의 대상이 되고 있다. 이에 따라 유역, 저수지, 및 하천을 포함한 포괄적인 대책을 수립하기 위하여 유역과 저수지 그리고 저수지와 하천을 연계한 통합시스템을 구축함으로써 탁수의 거동을 고려한 저수지 운영에 대한 의사결정을 지원하고자 한다. 본 연구에서는 안동 및 임하댐의 연계운영을 통한 하류하천 탁도관리 방안을 수립하고자 동적수질모형인 Koriv1-win 모형을 구축하여 저수지 수질모의 시스템과 연계 가능하도록 하였다. 모형은 안동댐하류로부터 구미수위표 지점까지 134 km에 달하는 구간에 대하여 구축하였으며 대상구간 내에는 내성천을 비롯 5개의 주요지천이 존재하게 된다. 그림 1에는 수리모듈 및 수질모듈의 검증결과를 나타내고 있으며 탁도모의결과를 검증하기 위하여 임하댐 관리단에서 측정한 자료를 이용하였다. 수리해석결과 구미지점의 유량변동을 잘 반영하는 것으로 모의되었고 수질해석결과도 관측치와 상당히 비슷하게 나타나는 것을 알 수 있었다. 한편 안동댐방류수와 임하댐 방류수가 혼합되는 영역에서 2차원적 흐름양상과 각 댐으로부터 방류되는 방류수의 수온과 탁도의 혼합양상을 모의하였으며 실측을 통하여 모의결과와 비교해 보았다. 저수지 하류하천에 대하여 1차원 및 2차원적 수리해석과 수질해석결과를 통하여 하천에서의 탁수전파 거리 및 고탁도 탁수의 도달 시간을 산정할 수 있도록 함으로써 저수지 연계운영을 통한 하류단 탁도저감에 대한 계획수립이 가능하다.을 것으로 판단된다. 철거 선정 체계는 보 철거를 위해서 보 주변 지역 거주민을 설득하며, 하천의 환경 개선을 위한 합리적인 대안 제시를 목적으로 하고 있다. 선정 체계를 바탕으로 주민 협의 및 대안 제시를 통해 결정된 대상 보는 선정 체계 안에서 보 철거 영향 판단 절차에 따라서 보 철거로 인한 수문, 수리, 지형, 수질, 생태 영향을 판단하게 되며, 이와 더불어 사회 경제적인 영향을 평가하게 된다. 평가결과에 따라서 보를 완전히 철거하거나 다른 대안을 고려하여 보를 부분적으로 철거하거나 개량하게 된다.곳으로 1/3의 자기 생산을 담당하고 있었다. $\ulcorner$경상도지리지$\lrcorner$(慶尙道地理志)에는 $\ulcorner$세종실록$\lrcorner$(世宗實錄) $\ulcorner$지리지$\lrcorner$(地理志)와 동년대에 동일한 목적으로 찬술되었음을 알 수 있다. $\ulcorner$경상도실록지리지$\lrcorner$(慶尙道實錄地理志)에는 $\ulcorner$세종실록$\lrcorner$(世宗實錄) $\ulcorner$지리지$\lrcorner$(地理志)와의 비교를 해보면 상 중 하품의 통합 9개소가 삭제되어 있고, $\ulcorner$동국여지승람$\lrcorner$(東國與地勝覽) 에서는 자기소와 도기소의 위

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.1614-1618
• /
• 2007

분포형 수문모형은 유역을 격자단위로 세분화하여 매개변수를 부여하고, 증발산, 침투, 지표면유출, 중간유출, 지하수유출, 하도흐름 등 여러 가지 수문요소를 해석하는 종합적인 수문모형이다. 지표면에 내린 강우의 증발 및 침투, 유출은 토양수분의 함량에 크게 의존하게 된다. 따라서 토양수분에 대한 적절한 모의가 분포형 수문모형의 정확도를 좌우하는 핵심이라 할 수 있다. 토양수분은 식물의 생장 및 가용수자원 산정 등에 있어서 중요한 요소로서 토양층 상부의 수 미터 내에 존재하는 수분의 양을 일컫는다. 토양수분의 공간적 시간적 특성들은 증발, 침투, 지하수 재충전, 토양침식, 식생분포 등을 지배하는 매우 중요한 요소라 할 수 있다. 강우로 인한 지면과 지표하에서의 순간적인 포화공간의 형성 및 유출의 생성을 포함하는 과정과 증발산 등은 모두 비포화대(vadose zone) 혹은 토양층에서의 토양수분의 함량에 크게 의존하게 된다(이가영 등, 2005) 본 연구에서는 토양수분에 대한 정밀측정 자료가 있는 설마천 유역 범륜사 사면에 대하여 분포형 수문모형의 토양수분 해석 능력을 평가하고자 하였다. 토양수분 모의에 사용된 격자기반의 분포형 수문모형은 미공병단에서 개발한 GSSHA(USACE, 2006) 모형이다. 모형의 입력자료는 정밀토양도와 현장측정에 의한 토양매개변수를 반영하여 구축하였고, 강우 및 기상자료는 2003년 1월 1일 ${\sim}$ 2004년 12월 31일의 1시간 자료를 이용하였다. 모의기간 중 2003년은 초기 토양수분값 등 초기조건의 영향을 줄이기 위한 웜업 (Warm-up)기간으로 설정하였고, 2004년의 모의결과를 토양수분 관측값과 비교하였다.업지역, 상업지역 등과 같이 지형적 특성에 따른 유량측정망을 구축하는 것이다.의 의사결정 지원 도구가 될 것이다. 따라서, 본 연구에서는 도시유역의 물순환 해석을 위한 일련의 과정, 즉 자료의 조사 및 취득에서부터 물순환 해석 모형을 이용한 정량적 현황파악, 물순환 개선 기법 및 평가를 수행함에 있어 주요 착안점 및 실무에서의 기술적 가이드를 제공하고자 하였으며, 보다 세밀한 도시유역의 물순환 해석을 위하여 우리나라와 일본에서 적용이 활발한 물리적 기반의 분포형 모형(WEP, SHER, SWMM)의 적용사례를 통하여 국내 도시하천의 물순환 해석에 활용함에 있어서의 실질적인 적용절차 등을 제시하고자 하였다. 한다.호강유역의 급격한 수질개선을 알 수 있다.世宗實錄) $\ulcorner$지리지$\lrcorner$(地理志)와 동년대에 동일한 목적으로 찬술되었음을 알 수 있다. $\ulcorner$경상도실록지리지$\lrcorner$(慶尙道實錄地理志)에는 $\ulcorner$세종실록$\lrcorner$(世宗實錄) $\ulcorner$지리지$\lrcorner$(地理志)와의 비교를 해보면 상 중 하품의 통합 9개소가 삭제되어 있고, $\ulcorner$동국여지승람$\lrcorner$(東國與地勝覽) 에서는 자기소와 도기소의 위치가 완전히 삭제되어 있다. 이러한 현상은 첫째, 15세기 중엽 경제적 태평과 함께 백자의 수요 생산이 증가하자 군신의 변별(辨別)과 사치를 이유로 강력하게 규제하여 백자의 확대와 발전에