• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.207-215
• /
• 2015

Electrical resistivity surveys were conducted at two subsidence areas near and at limestone mine sites, respectively, in order to estimate their causes of subsidence and the regions of potential occurrence. In addition, the linkages of mine development with these subsidences were investigated by the rock engineering analysis. Two study areas have different geological setting. One study site is the land subsidence area, which contains clay and sandy soil near the limestone mine, The other study site is the land subsidence area located just above the mine, which is expected to be relevant to the limestone mine. As results of two-dimensional (2D) electrical resistivity surveys at the sites 1 and 2, low resistivity zones, which are 70 ~ 120 ohm-m and 20 ~ 50 ohm-m, respectively, were found under the subsidence zones. For the study site 1, the possibility of subsidence was confirmed by using three-dimensional (3D) inversion performed with 2D resistivity profiles. For the study site 2, the cause of the subsidence and the possibility of subsidence occurrence were confirmed by rock engineering computation with regard to measurement line 7 in which low resistivity accompanied by subsidence area was observed.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.64-73
• /
• 2015

Recently, hydraulic fracturing is used in various fields and microseismic monitoring is one of the best methods for judging where hydraulic fractures exist and how they are developing. When locating microseismic events using single vertical well data, distances from the vertical array and depths from the surface are generally decided using time differences between compressional (P) wave and shear (S) wave arrivals and azimuths are calculated using P wave hodogram analysis. However, in field data, it is sometimes hard to acquire P wave data which has smaller amplitude than S wave because microseismic data often have very low signal to noise (S/N) ratio. To overcome this problem, in this study, we developed a grid search algorithm which can find event location using all combinations of arrival times recorded at receivers. In addition, we introduced and analyzed the method which calculates azimuths using S wave. The tests of synthetic data show the inversion method using all combinations of arrival times and receivers can locate events without considering the origin time even using only single phase. In addition, the method can locate events with higher accuracy and has lower sensitivity on first arrival picking errors than conventional method. The method which calculates azimuths using S wave can provide reliable results when the dip between event and receiver is relatively small. However, this method shows the limitation when dip is greater than about $20^{\circ}$ in our model test.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.33-38
• /
• 2012

A resistivity method has been applied to wide range of engineering and environmental problems with the help of automatic and precise data acquisition. Thus, more accurate modeling and inversion of time-lapse monitoring data are required since resistivity monitoring has been introduced to quantitatively find out subsurface changes With respect to time. Here, we used the finite element method (FEM) for 3D resistivity modeling since the method is easy to realize complex topography and arbitrary shaped anomalous bodies. In the FEM, the linear elements, also referred to as first order elements, have certain advantages of simple formulation and narrow bandwidth of system equation. However, the linear elements show the poor accuracy and slow convergence of the solution with respect to the number of elements or nodes. To achieve the higher accuracy of finite element solution, high order elements are generally used. In this study, we developed a 3D resistivity modeling program using high order Serendipity elements. Comparing the Serendipity element solutions for a cube model with the linear element solutions, we assured that the Serendipity element solutions are more accurate than the linear element solutions in the 3D resistivity modeling.

• Journal of Korea Water Resources Association
• /
• v.39 no.2 s.163
• /
• pp.99-108
• /
• 2006

To describe the behavior of suspended-sediment particles in turbulent open-channel flows, the advection-diffusion equation or its simplified form has been used. Though this equation was derived upon several assumptions, only a few studies tried to evaluate the limit of the assumptions. The reason is that it is very difficult to measure turbulence in open-channel flows and to discriminate the velocities of water and sediment particles. The present study aims to measure the velocity profiles of water and sediment particles in open-channel flows by using PTV (Particle Tracking Velocimetry), a kind of PIV (Particle Image Velocimetry). The measured results showed that sediment particles moved slower than water tracers did in the outer region. In the present study, the amount of velocity-lag reached about $5\%$ of the mom flow velocity and the position of the maximum velocity-lag was $g/h\approx0.05\;(g^{+}=30\~50)$ The main cause of the velocity-lag of sediment particles seems that the sediment particles have larger density than water has. On the other hand, in the viscous sublayer, sediment particle has a larger velocity than water tracers. The reason of the inversion of velocity-lag may be due to the no-sleep condition of water at the solid boundaries.

• Membrane Journal
• /
• v.19 no.4
• /
• pp.261-267
• /
• 2009

Polyimide which is the glassy polymer has high chemical resistance, thermal stability and high mechanical property. In this study, the polyimide hollow fiber membranes were prepared by the dry-jet wet phase inversion in order to investigate the permeation porperties of the $H_2S$ and $CH_4$. The morphology of prepared hollow fiber membranes and their permeation behaviors of $H_2S$ and $CH_4$ before and after silicon coating were evaluated. The permeance of $H_2S$ and $H_2S/CH_4$ selectivity increased due to plasticization with increasing the feed pressure. The permeance of KSM03b and selectivity of KSM03d were highest among the three type membranes used this experiments. The permeance decreased but the $H_2S/CH_4$ selectivity increased with increasing the air gap. The permeance reduced after silicon coating. However, the selectivity increased and the selectivity of KSM03d was 275 at 7 atm.

• Journal of the Korean earth science society
• /
• v.30 no.2
• /
• pp.153-164
• /
• 2009

We examined the regional 1-D deep resistivity structure of the Korean Peninsula using MT data acquired at seven sites located in the Kyongsang Basin and Kyonggi Massif. At the sites located in the Kyongsang Basin, surrounding sea distorts observed MT response and hence this distortion, so called "sea effect", is corrected using an iterative tensor stripping method. The 1-D layered inversion results for the seven MT sites reveal 4 layered structure, which is composed of 1) near surface layer, 2) upper crust, 3) lower crust and upper mantle, and 4) asthenosphere from the surface downward. Conrad interface, which is a boundary between upper and lower crust, is distinctly identified beneath all the MT sites. Conrad interface depth is estimated to about be 17km in the Kyongsang Basin and about 12km in the Kyonggi Massif, while the upper crust of the Kyongsang Basin is about 5 times more resistive than that of the Kyonggi Massif. Finally, asthenosphere is inferred to exist below a depth of approximately 100km with a resistivity of 200-300 ohm-m.

• The Journal of Engineering Geology
• /
• v.11 no.2
• /
• pp.163-174
• /
• 2001

For quantitative evaluation of geotechnical engineering properties such as rippability and diggability, clear interpretation on the subsUJiace velocity structures should be preceded by figuring out top soil, weathered and soft rock layers, shape of basement, fracture zones, geologic boundary and etC. from the seismic refraction data. It is very important to set up suitable field parameters, which are the configuration of profile and its length, spacings of geophones and sources and topographic conditions, for increasing field data Quality. Geophone spacing of 3 to 5m is reconunended in the land slope area of house land development site. In refraction tomography technique, the number of source points should be more than a Cluarter of available channel number of instrument and the subsurface structure interpretation can be decreased the artifact of inversion by topographic effect. Compared with core logging data, it is shown that the velocity range of the soil is less than 700m/s, weathered rock 700~1,200m/s, soft rock 1,200~1,800m/s on the velocity tomogram section. And the upper limit of P-wave velocity for rippability is estimated 1,200 to 1,800m/s in land slope area of gneiss.

• Economic and Environmental Geology
• /
• v.35 no.3
• /
• pp.273-284
• /
• 2002

The Gravity-Geologic Method (GGM) was implemented for bathymetric determinations in the Drake Passage, Antarctica, using global marine Free-air Gravity Anomalies (FAGA) data sets by Sandwell and Smith (1997) and local echo sounding measurements. Of the 6548 bathymetric sounding measurements, two thirds of these points were used as control depths, while the remaining values were used as checkpoints. A density contrast of 9.0 gm/㎤ was selected based on the checkpoints predictions with changes in the density contrast assumed between the seawater and ocean bottom topographic mass. Control depths from the echo soundings were used to determine regional gravity components that were removed from FAGA to estimate the gravity effects of the bathymetry. These gravity effects were converted to bathymetry by inversion. In particular, a selective merging technique was developed to effectively combine the echo sounding depths with the GGM bathymetiy to enhance high frequency components along the shipborne sounding tracklines. For the rugged bathymetry of the research area, the GGM bathymetry shows correlation coefficients (CC) of 0.91, 0.92, and 0.85 with local shipborne sounding by KORDI, GEODAS, and a global ETOPO5 model, respectively. The enhanced GGM by selective merging shows imploved CCs of 0.948 and 0.954 with GEODAS and Smith & Sandwell (1997)'s predictions with RMS differences of 449.8 and 441.3 meters. The global marine FAGA data sets and other bathymetric models ensure that the GGM can be used in conjunction with shipborne bathymetry from echo sounding to extend the coverage into the unmapped regions, which should generate better results than simply gridding the sparse data or relying upon lower resolution global data sets such as ETOPO5.

• Economic and Environmental Geology
• /
• v.31 no.2
• /
• pp.127-139
• /
• 1998

A new technique of simultaneous inversion for 3-D seismic velocity structure by using direct, reflected, and refracted waves is applied to the southeast part of the Korean Peninsula including Pohang Basin, Kyongsang Basin and Ryongnam Massif. Pg, Sg, PmP, SmS, Pn, and Sn arrival times of 44 events with 554 seismic rays are inverted for locations and crustal structure. $6{\times}6$ with $0.5^{\circ}$ and 8 layers (4 km each layer) model was inverted. 3-D seismic crustal velocity tomography including eight sections from surface to Moho, ten profiles along latitude and longitude are analyzed. The results are as follows: 1) the average velocity and thickness of sediment are 5.04 km/s and 3-4 km, and the velocity of basement is 6.11 km/s. The shape of velocity in shallower layer is agreement with Bouguer gravity anomaly (Cho et al., 1997). 2) the velocities fluctuate strongly in the upper crust. The velocity distribution of the lower crust under Conrad appears basically horizontal. 3) the average depth of Moho is 30.4 km, and velocity is 8.01 km/s. 4) from the velocity and depth of the sediment, the thickness, velocity and form of the upper crust, and the depth and form of Moho, we can find the obvious differences among Ryongnam Massif, Kyongsang Basin and Pohang Basin. 5) the deep faults (a Ulsan series faults) near Kyongju and Pohang areas can be found to be normal and/or thrust faults with detachment extended to the bottom of the upper crust.