• Journal of the Korean Geotechnical Society
• /
• v.35 no.11
• /
• pp.111-119
• /
• 2019

The high-speed railway system is mainly composed of tunnel, bridge, and viaduct to meet the straightness needed for keeping the high speed up to 400 km/s. Seismic fragility for the high-speed railway infrastructure can be assessed as two ways: one way is studying each element of infrastructure analytically or numerically, but it requires lots of research efforts due to wide range of railway system. On the other hand, empirical method can be used to access the fragility of an entire system efficiently, which requires case history data. In this study, we collect the 2004 M_W 6.6 Niigata earthquake case history data to develop empirical seismic fragility function for a railway system. Five types of intensity measures (IMs) and damage levels are assigned to all segments of target system for which the unit length is 200 m. From statistical analysis, probability of exceedance for a certain damage level (DL) is calculated as a function of IM. For those probability data points, log-normal CDF is fitted using MLE method, which forms fragility function for each damage level of exceedance. Evaluating fragility functions calculated, we observe that T=3.0 spectral acceleration (SA_T3.0) is superior to other IMs, which has lower standard deviation of log-normal CDF and low error of the fit. This indicates that long-period ground motion has more impacts on railway infrastructure system such as tunnel and bridge. It is observed that when SA_T3.0 = 0.1 g, P(DL>1) = 2%, and SA_T3.0 = 0.2 g, P(DL>1) = 23.9%.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• Economic and Environmental Geology
• /
• v.15 no.2
• /
• pp.65-88
• /
• 1982

In the last ten years, marine geological and geophysical survey and research were conducted by Japanese, Russian and American scientists in the East Sea of Korea (Japan Sea). Many research results were published. However, regional research of the geology of the continental margin of the Korean Peninsula was not conducted. This study has made on attempt to classify submarine strata and stratigraphic boundaries. The study has revealed characters of submarine geology and structure. Isopach maps of each identified stratigraphic unit have been constructed as the results of this study. The study was conducted on the basis of analyses of marine seismic surveys carried out in the continental margin of the East Sea between Kangneung and Pohang. Three depositional basins were identified in the study area and they were named as, Mukho Basin, Hupo Basin and Pohang Basin. The Mukho Basin is developed in continental slope and shelf in the area between Kangneung and Samcheog. Quaternary and Pliocene sediments attain a maximum thickness of 900 m. Basement rocks are interpreted as granite and gneiss. They are correlated with granite-gneiss of the Taebaecksan Series of Pre-cambrian age and the Daebo granite of Jurassic age. The Hupo Basin is developed in the continental shelf between Uljin and Youngdeok. Quaternary and Pliocene sediments attain a maximum thickness of 600 m. Basement rocks were interpreted as granite and gneiss and they are correlated with metamorphic rocks of Pre-cambrian age and the Daebo granites, comprising the Ryongnam Massif. The Pohang Basin is developed in the area between Pohang and Gangu. This basin contains Miocene and older sediments. Basement rocks are not shown. Many faults are developed within the continental shelf and slope. These faults strike parallel with the coast line. A north-south direction is predominant in the southern study area. However, in the northern study area the faults strike north, and north-west. The faults are parallel to each other and are step faults down-thrown to the east or west, forming horst and graben structures which develop into sedimentary basins. Such faults caused the development of submarine banks along the boundary between the continental shelf and slope. This bank has acted as a barrier for deposition in the Hupo Basin. Paleozoic sedimentary rocks distributed widely in the adjacent land area are absent in the Mukho Basin. This suggests that the area of the basin was situated above the sea level until the Pliocene time. The study area contains Pliocene sediments in general. These sediments overlie the basement complex composed of metamorphic rocks, granites, Cretaceous (Kyongsang System) sedimentary rocks and Miocene sedimentary rocks. These facts lead to a conclusion that the continental shelf and slope of the study area were developed as a result of displacements along faults oriented parallel to the present coast line in the post Miocene time.

• The Journal of Engineering Geology
• /
• v.7 no.2
• /
• pp.139-149
• /
• 1997

Electrical resistivity monitoring methods were adopted to detect groundwater table change in alluvium. Numerical modelling test using finite element method(FEM) and field resisfivity monitoring were conducted in the study. The field monitoring data were acquired in the alluvium deposit site in Jeong-Dong Ri, Geum River where pumping test had been conducted continuously for 20 days to make artificial changes of groundwater table. The unit distance of the electrode array was 4m and 21 fixed electrodes were applied in numerical calculation and field data acquisition. "Modified Wenner" and dipole-dipole array configurations were used in the study. The models used in two-dimensional numerical test were designed on the basis of the simplifving geological model of the alluvium in Jeong Dong Ri, Geum River. Numerical test results show that the apparent resistivity pseudosections were changed in the vicinity of the pootion where groundwater table was changed. Furthermore, there are some apparent resistivity changes in the boundary between aquifer and crystalline basement rock which overlays the aquifer. The field monitoring data also give similar results which were observed in numerical tests. From the numerical test using FEM and field resistivity monitoring observations in alluvium site of Geum River, the electrical monitoring method is proved to be a useful tool for detecting groundwater behavior including groundwater table change. There are some limitations, however, in the application of the resistivity method only because the change of groundwater table does not give enough variations in the apparent resistivity pseudosections to estimate the amount of groundwater table change. For the improved detection of groundwater table changes, it is desirable to combine the resistivity method with other geophysical methods that reveal the underground image such as high-resolution seismic and/or ground penetrating radar surveys.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• The Journal of Engineering Geology
• /
• v.31 no.3
• /
• pp.381-393
• /
• 2021

Physical, mechanical, hydraulic, and geophysical tests were applied to validate methods of inspecting the effectiveness of grouting on an agricultural reservoir dam. Data obtained from series of in situ and laboratory tests considered four stages: before grouting; during grouting; immediately after grouting; and after aging the grouting for 28 days. The results of SPT and triaxial tests, including the unit weight, compressive strength, friction angle, cohesion, and N-value, indicated the extent of ground improvement with respect to grout injection. However, they sometimes contained errors caused by ground heterogeneity. Hydraulic conductivity obtained from in situ variable head permeability testing is most suitable for identifying the effectiveness of grouting because the impermeability of the ground increased immediately after grouting. Electric resistivity surveying is useful for finding a saturated zone and a seepage pathway, and multichannel analysis of surface waves (MASW) is suitable for analyzing the effectiveness of grouting, as elastic velocity increases distinctly after grouting injection. MASW also allows calculation from the P- and S- wave velocities of dynamic properties (e.g., dynamic elastic modulus and dynamic Poisson's ratio), which can be used in the seismic design of dam structures.

• Geophysics and Geophysical Exploration
• /
• v.24 no.1
• /
• pp.6-19
• /
• 2021

Subsurface physical properties can be obtained and imaged by seismic exploration, and various algorithms have been developed for this purpose. In this regard, root mean square error (RMSE) has been widely used to quantitatively evaluate the accuracy of the developed algorithms. Although RMSE has the advantage of being numerically simple, it has limitations in assessing structural similarity. To supplement this, full-reference image quality assessment (FR-IQA) techniques, which reflect the human visual system, are being investigated. Therefore, we selected six FR-IQA techniques that could evaluate the obtained physical properties. In this paper, we used the full-waveform inversion, because the algorithm can provide the physical properties. The inversion results were applied to the six selected FR-IQA techniques using three benchmark models. Using salt models, it was confirmed that the inversion results were not satisfactory in some aspects, but the value of RMSE decreased. On the other hand, some FR-IQA techniques could definitely improve the evaluation.

• Economic and Environmental Geology
• /
• v.45 no.6
• /
• pp.661-672
• /
• 2012

Lab scale experiments to investigate the dissolution reaction among supercritical $CO_2$-sandstone-groundwater by using sandstones from Gyeongsang basin were performed. High pressurized cell system (100 bar and $50^{\circ}C$) was designed to create supercritical $CO_2$ in the cell, simulating the sub-surface $CO_2$ storage site. The first-order dissolution coefficient ($k_d$) of the sandstone was calculated by measuring the change of the weight of thin section or the concentration of ions dissolved in groundwater at the reaction time intervals. For 30 days of the supercritical $CO_2$-sandstone-groundwater reaction, physical properties of sandstone cores in Gyeongsang basin were measured to investigate the effect of supercritical $CO_2$ on the sandstone. The weight change of sandstone cores was also measured to calculate the dissolution coefficient and the dissolution time of 1 g per unit area (1 $cm^2$) of each sandstone was quantitatively predicted. For the experiment using thin sections, mass of $Ca^{2+}$ and $Na^+$ dissolved in groundwater increased, suggesting that plagioclase and calcite of the sandstone would be significantly dissolved when it contacts with supercritical $CO_2$ and groundwater at $CO_2$ sequestration sites. 0.66% of the original thin sec-tion mass for the sandstone were dissolved after 30 days reaction. The average porosity for C sandstones was 8.183% and it increased to 8.789% after 30 days of the reaction. The average dry density, seismic velocity, and 1-D compression strength of sandstones decreased and these results were dependent on the porosity increase by the dissolution during the reaction. By using the first-order dissolution coefficient, the average time to dissolve 1 g of B and C sandstones per unit area (1 $cm^2$) was calculated as 1,532 years and 329 years, respectively. From results, it was investigated that the physical property change of sandstones at Gyeongsang basin would rapidly occur when the supercritical $CO_2$ was injected into $CO_2$ sequestration sites.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.6 no.2
• /
• pp.81-92
• /
• 2001

According to analyses of high-resolution seismic profiles (air gun, sparker, and SBP) and a deep-drill core(YSDP 105) in the mid-eastern Yellow Sea, stratigraphic and geoacoustic models have been established and seismo-acoustic modeling has been fulfilled using ray tracing of finite element method. Stratigraphic model reflects seismo-, litho-, and chrono-stratigraphic sequences formed under a significant influence of Quaternary glacio-eustatic sea-level fluctuations. Each sequence consists of terrestrial to very-shallow-marine coarse-grained lowstand systems tract and tidal fine-grained transgressive to highstand systems tract. Based on mean grain-size data (121 samples) of the drill core, bulk density and P-wave velocity of depositional units have been inferred and extrapolated down to a depth of the recovery using the Hamilton's regression equations. As goo-acoustic parameters, the 121 pairs of bulk density and P-wave velocity have been averaged on each unit of the stratigraphic model. As a result of computer ray-tracing simulation of the subsurface strata, we have found that there are complex ray paths and many acoustic-shadow zones owing to the presence of irregular layer boundaries and low-velocity layers.