• Geophysics and Geophysical Exploration
• /
• v.13 no.3
• /
• pp.268-276
• /
• 2010

Earthquake Hazard Mitigation Law was activated into force on March 2009. By the law, the obligation to monitor the effect of earthquake on the facilities was extended to many organizations such as gas company and local governments. Based on the estimation of National Emergency Management Agency (NEMA), the number of free-surface acceleration stations would be expanded to more than 400. The advent of internet protocol and the more simplified operation have allowed the quick and easy installation of seismic stations. In addition, the dynamic range of seismic instruments has been continuously improved enough to evaluate damage intensity and to alert alarm directly for earthquake hazard mitigation. For direct visualization of damage intensity and area, Real Time Intensity COlor Mapping (RTICOM) is explained in detail. RTICOM would be used to retrieve the essential information for damage evaluation, Peak Ground Acceleration (PGA). Destructive earthquake damage is usually due to surface waves which just follow S wave. The peak amplitude of surface wave would be pre-estimated from the amplitude and frequency content of first arrival P wave. Earthquake Early Warning (EEW) system is conventionally defined to estimate local magnitude from P wave. The status of EEW is reviewed and the application of EEW to Odesan earthquake is exampled with ShakeMap in order to make clear its appearance. In the sense of rapidity, the earthquake announcement of Korea Meteorological Agency (KMA) might be dramatically improved by the adaption of EEW. In order to realize hazard mitigation, EEW should be applied to the local crucial facilities such as nuclear power plants and fragile semi-conduct plant. The distributed EEW is introduced with the application example of Uljin earthquake. Not only Nation-wide but also locally distributed EEW applications, all relevant information is needed to be shared in real time. The plan of extension of Korea Integrated Seismic System (KISS) is briefly explained in order to future cooperation of data sharing and utilization.

• 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.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.30-36
• /
• 2006

Japan's first pilot-scale $CO_2$ sequestration experiment has been conducted in Nagaoka, where 10400 t of $CO_2$ have been injected in an onshore aquifer at a depth of about 1100 m. Among various measurements conducted at the site for monitoring the injected $CO_2$, we conducted time-lapse crosswell seismic tomography between two observation wells to determine the distribution of $CO_2$ in the aquifer by the change of P-wave velocities. This paper reports the results of the crosswell seismic tomography conducted at the site. The crosswell seismic tomography measurements were carried out three times; once before the injection as a baseline survey, and twice during the injection as monitoring surveys. The velocity tomograms resulting from the monitoring surveys were compared to the baseline survey tomogram, and velocity difference tomograms were generated. The velocity difference tomograms showed that velocity had decreased in a part of the aquifer around the injection well, where the injected $CO_2$ was supposed to be distributed. We also found that the area in which velocity had decreased was expanding in the formation up-dip direction, as increasing amounts of $CO_2$ were injected. The maximum velocity reductions observed were 3.0% after 3200 t of $CO_2$ had been injected, and 3.5% after injection of 6200 t of $CO_2$. Although seismic tomography could map the area of velocity decrease due to $CO_2$ injection, we observed some contradictions with the results of time-lapse sonic logging, and with the geological condition of the cap rock. To investigate these contradictions, we conducted numerical experiments simulating the test site. As a result, we found that part of the velocity distribution displayed in the tomograms was affected by artefacts or ghosts caused by the source-receiver geometry for the crosswell tomography in this particular site. The maximum velocity decrease obtained by tomography (3.5%) was much smaller than that observed by sonic logging (more than 20%). The numerical experiment results showed that only 5.5% velocity reduction might be observed, although the model was given a 20% velocity reduction zone. Judging from this result, the actual velocity reduction can be more than 3.5%, the value we obtained from the field data reconstruction. Further studies are needed to obtain more accurate velocity values that are comparable to those obtained by sonic logging.

• Journal of the Korean Geographical Society
• /
• v.42 no.4
• /
• pp.488-505
• /
• 2007

The purposes of this research are twofold; 1) to verify spatial differences of tectonic movement using the spatial distribution of earthquakes, and 2) to infer mechanisms that generate spatial accumulation patterns of earthquakes in the Korean Peninsula. The first part of this sequential paper (Park, 2007) argues that the Korean Peninsula consists of four geostructural regions in which tectonic deformation and consequent geomorphological development patterns are different from each other Since this conclusion has been made by terrain analyses alone, it is necessary to verify this suggestion using other independent geophysical data. Because earthquakes are results of movement and deformation of land masses moving in different directions, the distribution of earthquake epicenters may be used to identify the direction and rates of land mass movement. This paper first analysed the spatial distribution of earthquakes using spatial statistics, and then results were compared with the spatial arrangement of geostructural regions. The spatial distribution of earthquakes in the Korean Peninsula can be summarized as the followings; firstly, the intensity of earthquakes shows only weak spatial dependency, and shows large difference even at adjacent regions. Secondly, the epicenter distribution has a clear spatial accumulation pattern, even though the intensity of earthquake shows a random pattern. Thirdly, the high density area of earthquakes shows a clear 'L' shape, passing through Pyeongannam-do, centered at Pyeongyang, and Hwanghae-do, Seosan and Pohang. The correlation coefficient between the density of earthquakes and distance from geostructral region boundaries is much higher than those between the density of fault lines and distance from tectonic division boundaries. Since fault lines and tectonic divisions in the Korean Peninsula are the results of long-term geological development, there is an apparent scale discrepancy to find significant correlations with earthquakes. This result verifies the research hypothesis that the Korean Peninsula is divided into four geostructral regions in which each has its own moving direction and spatial deformation characteristics. The existence of geostructural regions is also supported by the movement parrerns of land masses estimated from the GPS measurements. This conclusion is expected to provide a new perspective to understand the geomorphological developments and the earthquake occurrences in the Korean Peninsula.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.164-173
• /
• 2004

To delineate geothermal water movement at the Pohang geothermal development site, Self-Potential (SP) survey and monitoring were carried out during pumping tests. Before drilling, background SP data have been gathered to figure out overall potential distribution of the site. The pumping test was performed in two separate periods: 24 hours in December 2003 and 72 hours in March 2004. SP monitoring started several days before the pumping tests with a 128-channel automatic recording system. The background SP survey showed a clear positive anomaly at the northern part of the boreholes, which may be interpreted as an up-flow Bone of the deep geothermal water due to electrokinetic potential generated by hydrothermal circulation. The first and second SP monitoring during the pumping tests performed to figure out the fluid flow in the geothermal reservoir but it was not easy to see clear variations of SP due to pumping and pumping stop. Since the area is covered by some 360 m-thick tertiary sediments with very low electrical resistivity (less than 10 ohm-m), the electrokinetic potential due to deep groundwater flow resulted in being seriously attenuated on the surface. However, when we compared the variation of SP with that of groundwater level and temperature of pumping water, we could identify some areas responsible to the pumping. Dominant SP changes are observed in the south-west part of the boreholes during both the preliminary and long-term pumping periods, where 3-D magnetotelluric survey showed low-resistivity anomaly at the depth of $600m\~1,000m$. Overall analysis suggests that there exist hydraulic connection through the southwestern part to the pumping well.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.34-40
• /
• 2008

Underwater seismic refraction with advanced interpretation approaches makes important contributions to shallow marine exploration and geotechnical investigations in Australia's coastal areas. A series of case studies are presented to demonstrate the recent applications of continuous and static USR methods to river crossing and port infrastructure projects at various sites around Australia. In Sydney, static underwater seismic refraction (USR) with bottom-placed receivers and borehole seismic imaging assisted the development of improved geotechnical models that reduced construction risk for a tunnel crossing of the Lane Cove River. In Melbourne, combining conventional boomer reflection and continuous USR with near-bottom sources and receivers improved the definition of a buried, variably weathered basalt flow and assisted dredging assessment for navigation channel upgrades at Geelong Ports. Sand quality assessment with continuous USR and widely spaced borehole information assisted commercial decisions on available sand resources for the reclamation phase of development at the Port of Brisbane. Buried reefs and indurated layers occur in Australian coastal sediments with the characteristics of laterally limited, high velocity, cap layers within lower velocity materials. If these features are not recognised then significant error in depth determination to deeper refractors can occur. Application of advanced refraction inversion using wavefront eikonal tomography to continuous USR data obtained along the route of a proposed offshore pipeline near Fremantle allowed these layers and the underlying bedrock refractor to be accurately imaged. Static USR and the same interpretation approach was used to image the drowned granitic regolith beneath sediments and indurated layers in the northern area of Western Australia at a proposed new berthing site where deep piling was required. This allowed preferred piling sites to be identified, reducing overall pile lengths. USR can be expected to find increased application to shallow marine exploration and geotechnical investigations in Australia's coastal areas as economic growth continues and improved interpretation methods are developed.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.52-59
• /
• 2008

We have developed a new deep-towed marine DC resistivity survey system. It was designed to detect the top boundary of the methane hydrate zone, which is not imaged well by seismic reflection surveys. Our system, with a transmitter and a 160-m-long tail with eight source electrodes and a receiver dipole, is towed from a research vessel near the seafloor. Numerical calculations show that our marine DC resistivity survey system can effectively image the top surface of the methane hydrate layer. A survey was carried out off Joetsu, in the Japan Sea, where outcrops of methane hydrate are observed. We successfully obtained DC resistivity data along a profile ${\sim}3.5\;km$ long, and detected relatively high apparent resistivity values. Particularly in areas with methane hydrate exposure, anomalously high apparent resistivity was observed, and we interpret these high apparent resistivities to be due to the methane hydrate zone below the seafloor. Marine DC resistivity surveys will be a new tool to image sub-seafloor structures within methane hydrate zones.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.1-14
• /
• 2008

We present a method for interpreting seismic records with arrivals and waveforms having characteristics which could be generated by extremely inhomogeneous velocity structures, such as non-typical oceanic crust, decollement at subduction zones, and seamounts in oceanic regions, by comparing them with synthetic waveforms. Recent extensive refraction and wide-angle reflection surveys in oceanic regions have provided us with a huge number of high-resolution and high-quality seismic records containing characteristic arrivals and waveforms, besides first arrivals and major reflected phases such as PmP. Some characteristic waveforms, with significant later reflected phases or anomalous amplitude decay with offset distance, are difficult to interpret using only a conventional interpretation method such as the traveltime tomographic inversion method. We find the best process for investigating such characteristic phases is to use an interactive interpretation method to compare observed data with synthetic waveforms, and calculate raypaths and traveltimes. This approach enables us to construct a reasonable structural model that includes all of the major characteristics of the observed waveforms. We present results here with some actual observed examples that might be of great help in the interpretation of such problematic phases. Our approach to the analysis of waveform characteristics is endorsed as an innovative method for constructing high-resolution and high-quality crustal structure models, not only in oceanic regions, but also in the continental regions.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.146-162
• /
• 2017

Multibeam and seismic data in the northwestern part of the Ulleung Basin were analyzed to study stratigraphy and evolutionary history of submarine canyon. A detailed analysis reveals that the sedimentary sequences in this area consist of four stratigraphic units separated by erosional unconformities. On the continental slope, these units are dominated by well-stratified facies with some slope failures, whereas these units show well-stratified and chaotic facies toward the basin floor. Generally, the sediment thickness is relatively thin on the slope, whereas thick sediment accumulation occurs on the base of slope and basin floor. Based on seismic characteristics and distribution, the deposition of each units are well correlated with the evolutionary history of the submarine canyon. Unit 1 directly overlying the acoustic basement has thin sediment layer on the slope, whereas its thickness gradually increase toward the basin floor. Compared to other units, Unit 2 is relatively thick accumulations on the slope and contains some slope failures related to faults systems. The mass transport sediments due to slope failures, mainly deposited on the base of slope as a submarine fan. The width and depth of submarine canyon increase due to dominant of the erosional process rather than the sediment deposition. Unit 3 is thin accumulation on the slope around the submarine canyon. Toward the basin floor, its thickness gradually increases. Unit 4 is characterized by thin layers including slides and slumps on the slope, whereas it formed thick accumulations at the base of slope as a submarine fan. The increase in the width and depth of submarine canyon results from the dominant of the erosional process and slope failures around the submarine canyon. Consequently, the formation of sedimentary units combined with the development of submarine canyon in this area is largely controlled by the amounts of sediment supply originated from slope failures, regional tectonic effects and sea-level fluctuations.

• Economic and Environmental Geology
• /
• v.49 no.5
• /
• pp.381-388
• /
• 2016

The purpose of this study is to integrate and visualize using mapping techniques based on precise seabed geomorphology, seafloor backscattering images and high-resolution underwater images of the nearshore area around the Dokdo, in the East Sea. We have been obtained the precise topography map using multibeam echosounder system around the nearshore area(~50 m) of the southern part of the Seodo. Side scan sonar survey for analysis seafloor backscattering images was carried out in the same area of topography data. High-resolution underwater images(zone(a), zone(b), zone(c)) were taken in significant habitat scope of the nearshore area of the southern part of the Seodo. Using the results of bathymetry, seafloor backscattering images, high-resolution underwater images, we performed an integrated visualization about the nearshore area of the Dokdo. The integrated visualizing techniques are possible to make the seabed characteristic mapping results of the nearshore area of the Dokdo. The integrated visualization results present more complex and reliable information than separate geological products for seabed environmental mapping study and it is useful to understand the relation between seafloor characteristics and topographic environments of the study area. The integrated visualizing techniques and mapping analysis need to study sustainably and periodically, for effective monitoring of the nearshore ecosystem of the Dokdo.