• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.203-208
• /
• 2007

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms for the travel time computation commonly used, however, may not give sufficiently precise results particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. We assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution.

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.327-335
• /
• 2006

Three processing strategies of seismic refraction data are tested in terms of velocity and depth profiles or structures for mapping of geological discontinuities: GRM(generalized reciprocal method), GLI(generalized linear inversion), Tomography. The test data used in this study are the shot gathers reconstructed by numerical modeling for the structures of 3 planar layers(horizontal, inclined), the buried vertical fracture zones, and vertical fault zones. Tomography is shown to be very efficient for mapping of more complicated tone such as vertical fault and buried fracture zones, whereas GRM and GLI can be useful for horizontal and/or inclined layers, probably on their bases of analysis of first arrivals in travel time curves.

• The Bulletin of The Korean Astronomical Society
• /
• v.43 no.1
• /
• pp.45.1-45.1
• /
• 2018

Formation of filaments and subsequent dense cores in ISM is one of the essential questions to address in star formation. To investigate this scenario in detail, we recently started a molecular line survey namely 'Filaments, the Universal Nursery of Stars (FUNS)' toward nearby filamentary clouds in Gould Belt using TRAO 14m single dish telescope equipped with a 16 multi-beam array. In the present work, we report the first look results of kinematics of a low mass star forming region L1478 of California molecular cloud. This region is found to be consisting of long filaments with a hub-filament structure. We performed On-The-Fly mapping observations covering ~1.1 square degree area of this region using C18O(1-0) as a low density tracer and 0.13 square degree area using N2H+(1-0) as a high density tracer, respectively. CS (2-1) and SO (32-21) were also used simultaneously to map ~290 square arcminute area of this region. We identified 10 filaments applying Dendrogram technique to C18O data-cube and 13 dense cores using FellWalker and N2H+ data set. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that filaments in L~1478 are velocity coherent and supercritical. Especially the filaments which are highly supercritical are found to have dense cores detected in N2H+. Non-thermal velocity dispersions derived from C18O and N2H+ suggest that most of the dense cores are subsonic or transonic while the surrounding filaments are transonic or supersonic. We concluded that filaments in L~1478 are gravitationally unstable which might collapse to form dense cores and stars. We also suggest that formation mechanism can be different in individual filament depending on its morphology and environment.

• Geomechanics and Engineering
• /
• v.13 no.6
• /
• pp.893-906
• /
• 2017

Automatic large scale soil model generation is very critical stage for earthquake hazard simulation of urban areas. Manual model development may cause some data losses and may not be effective when there are too many data from different soil observations in a wide area. Geographic information systems (GIS) for storing and analyzing spatial data help scientists to generate better models automatically. Although the original soil observations were limited to soil profile data, the recent developments in mapping technology, interpolation methods, and remote sensing have provided advanced soil model developments. Together with advanced computational technology, it is possible to handle much larger volumes of data. The scientists may solve difficult problems of describing the spatial variation of soil. In this study, an algorithm is proposed for automatic three dimensional soil and velocity model development of southern part of the European side of Istanbul next to Sea of Marmara based on GIS data. In the proposed algorithm, firstly bedrock surface is generated from integration of geological and geophysical measurements. Then, layer surface contacts are integrated with data gathered in vertical borings, and interpolations are interpreted on sections between the borings automatically. Three dimensional underground geology model is prepared using boring data, geologic cross sections and formation base contours drawn in the light of these data. During the preparation of the model, classification studies are made based on formation models. Then, 3D velocity models are developed by using geophysical measurements such as refraction-microtremor, array microtremor and PS logging. The soil and velocity models are integrated and final soil model is obtained. All stages of this algorithm are carried out automatically in the selected urban area. The system directly reads the GIS soil data in the selected part of urban area and 3D soil model is automatically developed for large scale earthquake hazard simulation studies.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.111-120
• /
• 2016

In order to enhance the connectivity of fracture network as fluid path in enhanced/engineered geothermal system (EGS), the exact locating of hydraulic fractured zone is very important. Hydraulic fractures can be tracked by locating of microseismic events which are occurred during hydraulic fracture stimulation at each stage. However, since the subsurface velocity is changed due to hydraulic fracturing at each stage, in order to find out the exact location of microseismic events, we have to consider the velocity change due to hydraulic fracturing at previous stage when we perform the mapping of microseimic events at the next stage. In this study, we have modified 3D locating algorithm of microseismic data which was developed by Kim et al. (2015) and have developed 3D velocity update algorithm using occurred microseismic data. Eikonal equation which can efficiently calculate traveltime for complex velocity model at anywhere without shadow zone is used as forward engine in our inversion. Computational cost is dramatically reduced by using Fresnel volume approach to construct Jacobian matrix in velocity inversion. Through the numerical test which simulates the geothermal survey geometry, we demonstrated that the initial velocity model was updated by using microseismic data. In addition, we confirmed that relocation results of microseismic events by using updated velocity model became closer to true locations.

• The Journal of Engineering Geology
• /
• v.31 no.4
• /
• pp.507-522
• /
• 2021

Specific survey objectives often cannot be met using only one geophysical method, as each method's results are influenced by the specific physical properties of subsurface materials. In particular, areas susceptible to geological hazards require investigation using more than one method in order to reduce risks to life and property. Instead of analyzing the results from each method separately, this work develops a four-quadrant criterion for classifying areas of levees as safe or weak. The assessment is based on statistically determined thresholds of seismic velocity (P-wave velocity from seismic refraction and S-wave velocity from multichannel analysis of surface waves) and electrical resistivity. Thresholds are determined by subtracting the standard deviation from the mean during performance testing of this correlation technique applied to model data of four horizontal and inclined fracture zones. Compared with results from the crossplot of resistivity and P-wave velocity, crossplot analysis using resistivity and S-wave velocity data provides more reliable information on the soil type, ground stiffness, and lithological characteristics of the levee system. A loose and sandy zone (represented by low S-wave velocity and high resistivity) falling within the second quadrant is interpreted to be a weak zone. This interpretation is well supported by the N values from standard penetrating test for the central core.

• Journal of KIISE:Computer Systems and Theory
• /
• v.36 no.5
• /
• pp.404-411
• /
• 2009

We introduce in this paper a new method for smooth foldover-free warping of images, based on the vector field deformation technique proposed by Von Funck et al. It allows users to specify the constraints in two different ways: positional constraints to constrain the position of a point in the image and gradient constraints to constrain the orientation and scaling of some parts of the image. From the user-specified constraints, it computes in the image domain a C1-continuous velocity vector field, along which each pixel progressively moves from its original position to the target. The target positions of the pixels are obtained by solving a set of partial derivative equations with the 4th order Runge-Kutta method. We show how our method can be useful for texture mapping with hard constraints. We start with an unconstrained planar embedding of a target mesh using a previously known method (Least Squares Conformal Map). Then, in order to obtain a texture map that satisfies the given constraints, we use the proposed warping method to align the features of the texture image with those on the unconstrained embedding. Compared to previous work, our method generates a smoother texture mapping, offers higher level of control for defining the constraints, and is simpler to implement.

• The Journal of the Convergence on Culture Technology
• /
• v.9 no.1
• /
• pp.539-544
• /
• 2023

The subject of this study is the Maae Buddha statue in granodiorite of the Mesozoic Cretaceous period, which is concerned about stability as a standing stone cultural property located in ◯◯-dong, Gyeongsangbuk-do. For stability analysis, three-dimensional face mapping, geological properties of joints, three-dimensional scanning, ultrasonic velocity, polarization microscopy, electron microscopy analysis and XRD analysis were performed. In addition, the safety factor of the Maaebul was calculated by analyzing the damage status investigation, stereographic projection analysis, rock classification, and limit equilibrium analysis. The types and scales of damage and possible collapse by section depend on the degree of weathering of the rock and the orientation and characteristics of the joints, but wedge-failure and toppling-failure are expected to be small-scale. The safety factor of Maaebul in dry and wet conditions is less than 1.2, so stability is concerned. The types of damage were mainly observed, such as exfoliation, cracking, granular decomposition, and vegetation growth. The Maaebul rock is granodiorite, and the surface discoloration materials are K, Fe, and Mg. The 4 sets of joints are developed, J1 is tensile joint and the others are shear joint. The uniaxial compressive strength estimated by ultrasonic exploration is 514kgf/cm2, which corresponds to most soft rocks and some weathered rocks. Rock classification(RMR) is estimated to be grade 5, very poor rock mass. These technique along with the existing methods of safety diagnosis of cultural properties are expected to be a reasonable tool for objective interpretation and stability review of stone cultural properties.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2003.03a
• /
• pp.359-364
• /
• 2003

For safety evaluation of a rock-fill dim, it is often necessary to investigate spatial distribution of weak zones such as fracture. Both DC-resistivity survey and seismic(SASW) method are usually used for the purpose. Recently, Multichannel analysis of surface waves(MASW) method which makes up for the weak point of SASW method is developed and the site examination which is simple came to be possible comparatively. In order to obtain 2-D shear-wave velocity(Vs) profile along the dam axis that can be associated with dynamic properties of filled materials, MASW method was adapted. Then, DC-resistivity survey and drilling survey were performed to compare with each results. We confirmed that the MASW method and DC-resistivity survey show complementary result that corresspond with drilling result. Therefore, MASW method is an efficient method for dynamic characterization of dam-filling materials and also the combination of related methods such as DC-resistivity can lead to an effective safety evaluation of rock-fill dam.

• Journal of Mechanical Science and Technology
• /
• v.16 no.12
• /
• pp.1652-1663
• /
• 2002

It is desirable to perform nondestructive evaluation to assess material properties and part homogeneity because manufacturing of carbon/carbon (C/C) composites requires complicated and costly processes. In this work several ultrasonic techniques were applied to carbon/carbon composites for the evaluation of spatial variations in material properties that are attributable to the manufacturing process. In a large carbon/carbon composite manufactured by chemical vapor infiltration (CVI) method, the spatial variation of ultrasonic velocity was measured and found to be consistent with the densification behavior in CVI process in order to increase the density of C/C composites. Ultrasonic velocity and attenuation depend on a density variation of materials. Low frequency through-transmission scans based on both amplitude and time-of-flight of the ultrasonic pulse were used for mapping out the material property inhomogeneity These results were compared with that obtained by dry-coupling ultrasonics. Pulse-echo C-scans was used to image near-surface material property anomalies such as the placement of spacers between disks during CVI. Also, optical micrograph had been examined on the surface of C/C composites using a destructive way.