• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.127-135
• /
• 1998

The Born approximation is widely used for solving the complex scattering problems in electromagnetics. Approximating total internal electric field by the background field is reasonable for small material contrasts as long as scatterer is not too large and the frequency is not too high. However in many geophysical applications, moderate and high conductivity contrasts cause both real and imaginary part of internal electric field to differ greatly from background. In the extended Born approximation, which can improve the accuracy of Born approximation dramatically, the total electric field in the integral over the scattering volume is approximated by the background electric field projected to a depolarization tensor. The finite difference and elements methods are usually used in EM scattering problems with a 2D model and a 3D source, due to their capability for simulating complex subsurface conductivity distributions. The price paid for a 3D source is that many wavenumber domain solutions and their inverse Fourier transform must be computed. In these differential equation methods, all the area including homogeneous region should be discretized, which increases the number of nodes and matrix size. Therefore, the differential equation methods need a lot of computing time and large memory. In this study, EM modeling program for a 2D model and a 3D source is developed, which is based on the extended Born approximation. The solution is very fast and stable. Using the program, crosshole EM responses with a vertical magnetic dipole source are obtained and the results are compared with those of 3D integral equation solutions. The agreement between the integral equation solution and extended Born approximation is remarkable within the entire frequency range, but degrades with the increase of conductivity contrast between anomalous body and background medium. The extended Born approximation is accurate in the case conductivity contrast is lower than 1:10. Therefore, the location and conductivity of the anomalous body can be estimated effectively by the extended Born approximation although the quantitative estimate of conductivity is difficult for the case conductivity contrast is too high.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.2
• /
• pp.231-243
• /
• 2016

This study presents the applicability of two-dimensional (2D) flood inundation model by applying to real irrigation reservoir failure with limited available data. The study area is Sandae Reservoir placed in Gyeongju and downstream area of it and the reservoir was failured by piping in 2013. The breach hydrograph was estimated from one-dimensional (1D) hydrodynamic model and the discharge was employed for upstream boundary of 2D flood inundation model. Topography of study area was generated by integrating digital contour map and satellite data, and Cartesian grids with 3m resolution to consider geometry of building, road and public stadium were used for 2D flood inundation analysis. The model validation was carried out by comparing predictions with field survey data including reservoir breach outflow, flood extent, flood height and arrival time, and identifying rational ranges with allowed error. In addition, the applicability of 2D model is examined using different simulation conditions involving grid size, building and roughness coefficient. This study is expected to contributed to analysis of irrigation reservoirs were at risk of a failure and setting up Emergency Action Plan (EAP) against irrigation reservoir failure.

• Korean Journal of Optics and Photonics
• /
• v.20 no.4
• /
• pp.211-216
• /
• 2009

Much research into two-dimensional (2-D) photonic crystal (PC) structures has been conducted for realization of ultrasmall optical integrated circuits. A 2-D silicon (Si) PC slab structure with air cladding (n=1) is one of the representative structures in 2-D PCs. While air-clad Si PC slab structures have good optical characteristics, their suspension in air can lead to mechanical weakness, making integration with some optical devices difficult. In this paper, we propose improving the mechanical robustness of PC structure by developing a 2-D Si PC structure with symmetric silica cladding (n=1.44) and comparing its optical properties to that of the air-clad structure. First, we investigate the optical properties of a 2-D Si PC slab structure with air cladding by using a 3-D finite difference time domain method. We determined that a photonic bandgap of 330 nm and a non-leaky propagating bandwidth of 100 nm in the optical communication range are possible. Next, we investigate the optical properties of 2-D Si PC slab structures with silica cladding. Even though the refractive index of the silica cladding is higher than that of air, we developed a silica-clad structure with good optical properties: a photonic band gap of approximately 230 nm and a non-leaky propagating bandwidth of 90 nm, comparable to that of the air-clad PC structures.

• Journal of Korea Water Resources Association
• /
• v.46 no.12
• /
• pp.1221-1234
• /
• 2013

The purpose of this study is to suggest the methodology for the computation of uplift pressure and discharge of the seepage flow under gravity dam. A 3-dimensional FDM model is developed for this purpose and this model can simulate the saturated Darcian flow in heterogeneous media. For the verification of the numeric model, test simulation has been executed and the mass balance has been checked. The error does not exceed 3%. Using the developed model, The uplift pressure and seepage flow discharge under gravity dam has been calculated. The uplift pressure shows the similar pattern, comparing with the result of flow-net method. As the length of grout curtain increases, the uplift pressure decreases linearly, but the seepage flow discharge shows the non-linear decreasing pattern. The coefficients of the formulas in the dam-design criteria have been analysed, and ${\alpha}=1/3$ corresponds to the value when the length of curtain grout is 70% of the aquifer height. The uplift pressure near the pressure relief drain has the big curvature vertically and horizontally. The developed model in this study can be used for the evaluation of the effects of seepage flow under gravity dam.

• Journal of The Geomorphological Association of Korea
• /
• v.19 no.3
• /
• pp.37-49
• /
• 2012

The objective of this paper is to propose a methodology for surface sediments classification in tidal flats that can combine ground survey data with high-resolution remote sensing data by multivariate kriging. Unlike conventional methodologies that have classified remote sensing data by using pre-classified sediment components, a new classification methodology presented in this paper first generates sediment component fraction maps and then classifies the sediments on a final stage. For generating sediment component fractions, regression kriging, as one of multivariate kriging algorithms, is applied to integrate ground survey data and remote sensing data. First, trend components of sand, silt, and clay are derived through regression analysis of ground survey data and spectral information from remote sensing data. Then, residuals at sample locations are computed and interpolated to generate residual components in the study area. Finally, the sediment component fractions are computed by adding the residuals to the trend components and are classified on a final stage. A case study at the Baramarae tidal flats with KOMPSAT-2 imagery is carried out to evaluate the classification capability of the proposed classification methodology. Through the case study, the proposed methodology showed the best classification accuracy, compared with the conventional classification methodologies. Especially, much improvement of classification accuracy for fine-grained sediments were also obtained. Therefore, it is expected that the presented classification methodology would be an effective one for surface sediments classification in tidal flats.

• Korean Journal of Remote Sensing
• /
• v.35 no.6_1
• /
• pp.933-944
• /
• 2019

This study analyzes the tide deformation of land boundary regions on the east (Region A) and west (Region B) sides of the Ross Ice Shelf in Antarctica using Double-Differential Interferometric Synthetic Aperture Radar (DDInSAR). A total of seven Sentinel-1A SAR images acquired in 2015-2016 were used to estimate the accuracy of tide prediction model and Young's modulus of ice shelf. First, we compared the Ross Sea Height-based Tidal Inverse (Ross_Inv) model, which is a representative tide prediction model for the Antarctic Ross Sea, with the tide deformation of the ice shelf extracted from the DDInSAR image. The accuracy was analyzed as 3.86 cm in the east region of Ross Ice Shelf and it was confirmed that the inverse barometric pressure effect must be corrected in the tide model. However, in the east, it is confirmed that the tide model may be inaccurate because a large error occurs even after correction of the atmospheric effect. In addition, the Young's modulus of the ice was calculated on the basis of the one-dimensional elastic beam model showing the correlation between the width of the hinge zone where the tide strain occurs and the ice thickness. For this purpose, the grounding line is defined as the line where the displacement caused by the tide appears in the DDInSAR image, and the hinge line is defined as the line to have the local maximum/minimum deformation, and the hinge zone as the area between the two lines. According to the one-dimensional elastic beam model assuming a semi-infinite plane, the width of the hinge region is directly proportional to the 0.75 power of the ice thickness. The width of the hinge zone was measured in the area where the ground line and the hinge line were close to the straight line shown in DDInSAR. The linear regression analysis with the 0.75 power of BEDMAP2 ice thickness estimated the Young's modulus of 1.77±0.73 GPa in the east and west of the Ross Ice Shelf. In this way, more accurate Young's modulus can be estimated by accumulating Sentinel-1 images in the future.

• Journal of Radiation Protection and Research
• /
• v.41 no.1
• /
• pp.31-39
• /
• 2016

Background: It is necessary to consider the overall countermeasure for analysis of nuclear activities according to the increase of the nuclear facilities like nuclear power and reprocessing plants in the neighboring countries including China, Taiwan, North Korea, Japan and South Korea. South Korea and comprehensive nuclear-test-ban treaty organization (CTBTO) are now operating the monitoring instruments to detect radionuclides released into the air. It is important to estimate the origin of radionuclides measured using the detection technology as well as the monitoring analysis in aspects of investigation and security of the nuclear activities in neighboring countries. Materials and methods: A three-dimensional forward/backward trajectory model has been developed to estimate the origin of radionuclides for a covert nuclear activity. The developed trajectory model was composed of forward and backward modules to track the particle positions using finite difference method. Results and discussion: A three-dimensional trajectory model was validated using the measured data at Chernobyl accident. The calculated results showed a good agreement by using the high concentration measurements and the locations where was near a release point. The three-dimensional trajectory model had some uncertainty according to the release time, release height and time interval of the trajectory at each release points. An atmospheric dispersion model called long-range accident dose assessment system (LADAS), based on the fields of regards (FOR) technique, was applied to reduce the uncertainties of the trajectory model and to improve the detective technology for estimating the radioisotopes emission area. Conclusion: The detective technology developed in this study can evaluate in release area and origin for covert nuclear activities based on measured radioisotopes at monitoring stations, and it might play critical tool to improve the ability of the nuclear safety field.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.183-193
• /
• 2000

A series of studies such as geological logging data analysis, detailed geological survey, rock mass evaluation, in-situ and laboratory tests, rock strength and mechanical properties of the rock were concerned. The stability of the slope were carried out inorder to design the pit slope and individual benches using the stereographic projection analysis and numerical methods in Roto Pit of Pasir coal fetid. The bedding plane was one of the major discontinuities in the Roto Pit and the dip of which is about $60^{\circ}$in the northern part and $83^{\circ}$in the southern part. The dip of bedding becomes steeper from north to south. The plane and toppling failures are presented in many slopes. In laboratory test the average uniaxial compressive strength of mudstone was 9 MPa and that of weak sandstone was 10 MPa. In-situ test showed that the rocks of Roto north mining area are mostly weak enough to be classified in grade from R2(weak) to R3(medium strong weak) and the coal is classified in grades from R1(Very weak) to R2(Weak). The detailed stability analysis were carried out on 4 areas of Roto north(east, west, south and north), and 2 areas of Roto south(east and west). In this paper, the minimum factor of safety was set to 1.2 which is a general criterion for open pit mines. Using the stereographic projection analysis and the limit equilibrium method, slope angles were calculated as 30~$36^{\circ}$for a factor of safety greater than 1.2. Then these results were re-evaluated by numerical analysis using FLAC. The final slope angles were determined by rational described abode. A final slope of 34 degrees can guarantee the stability for the eastern part of the Roto north area, 33 degrees for the western part, 35 degrees for the northern part and 35 degrees for the southern part. For the Roto south area, 36 degrees was suggested for both sides of the pit. Once the pit slope is designed based on the stability analysis and the safety measures. the stability of 니ope should be checked periodically during the mining operations. Because the slope face will be exposed long time to the rain fall, a study such aspreventive measures against weathering and erosion is highly recommended to be implemented.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.430-440
• /
• 2000

A series of studies such as geological logging data analysis, detailed geological survey, rock mass evaluation, in-situ and laboratory tests, rock strength and mechanical properties of the rock were concerned. The stability of the slope were carried out inorder to design the pit slope and individual benches using the stereographic projection analysis and numerical methods in Roto Pit of Pasir coal field. The bedding plane was one of the major discontinuities in the Roto Pit and the dip of which is about 60$^{\circ}$ in the northern part and 83$^{\circ}$ in the southern part. The dip of bedding becomes steeper from north to south. The plane and toppling failures are presented in many slopes. In laboratory test the average uniaxial compressive strength of mudstone was 9MPa and that of weak sandstone was 10MPa. In-situ test showed that the rocks of Roto north mining area are mostly weak enough to be classified in grade from R2(weak) to R3(medium strong weak) and the coal is classified in grades from R1(Very weak) to R2(Weak). The detailed stability analysis were carried out on 4 areas of Roto north (east, west, south and north), and 2 areas of Roto south(east and west). In this paper, the minimum factor of safety was set to 1.2 which is a general criterion for open pit mines. Using the stereographic projection analysis and the limit equilibrium method, slope angles were calculated as 30∼36$^{\circ}$ for a factor of safety greater than 1.2. Then these results were re-evaluated by numerical analysis using FLAC. The final slope angles were determined by rational described above. A final slope of 34 degrees can guarantee the stability for the eastern part of the Roto north area, 33 degrees for the western part, 35 degrees for the northern part and 35 degrees for the southern part. For the Roto south area, 36 degrees was suggested for both sides of the pit. Once the pit slope is designed based on the stability analysis and the safety measures, the stability of slope should be checked periodically during the mining operations. Because the slope face will be exposed long time to the rain fall, a study such aspreventive measures against weathering and erosion is highly recommended to be implemented.

• Korean Journal of Remote Sensing
• /
• v.36 no.5_1
• /
• pp.847-860
• /
• 2020

The 2-pass DInSAR (Differential Interferometric SAR) processing steps for DEM generation consist of the co-registration of SAR image pair, interferogram generation, phase unwrapping, calculation of DEM errors, and geocoding, etc. It requires complicated steps, and the accuracy of data processing at each step affects the performance of the finally generated DEM. In this study, we developed an improved method for enhancing the performance of the DEM generation method based on the 2-pass DInSAR technique of TanDEM-X bistatic SAR images was developed. The developed DEM generation method is a method that can significantly reduce both the DEM error in the unwrapped phase image and that may occur during geocoding step. The performance analysis of the developed algorithm was performed by comparing the vertical accuracy (Root Mean Square Error, RMSE) between the existing method and the newly proposed method using the ground control point (GCP) generated from GPS survey. The vertical accuracy of the DInSAR-based DEM generated without correction for the unwrapped phase error and geocoding error is 39.617 m. However, the vertical accuracy of the DEM generated through the proposed method is 2.346 m. It was confirmed that the DEM accuracy was improved through the proposed correction method. Through the proposed 2-pass DInSAR-based DEM generation method, the SRTM DEM error observed by DInSAR was compensated for the SRTM 30 m DEM (vertical accuracy 5.567 m) used as a reference. Through this, it was possible to finally create a DEM with improved spatial resolution of about 5 times and vertical accuracy of about 2.4 times. In addition, the spatial resolution of the DEM generated through the proposed method was matched with the SRTM 30 m DEM and the TanDEM-X 90m DEM, and the vertical accuracy was compared. As a result, it was confirmed that the vertical accuracy was improved by about 1.7 and 1.6 times, respectively, and more accurate DEM generation was possible with the proposed method. If the method derived in this study is used to continuously update the DEM for regions with frequent morphological changes, it will be possible to update the DEM effectively in a short time at low cost.