• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.357-363
• /
• 2010

In this study, we investigated the geologic structure of the basement and overlying sediments of the construction site of the dinosaur egg fossil museum in Hwasung, Gyeonggi Province through refraction seismology, drilling, and downward seismic velocity measurements in the drill holes. The construction site ($350{\times}750\;m^2$) is located in the reclaimed area south of Sihwa Lake, Gojeong-ri. About 6,950 m of seismic refraction data consisting of 11 lines were acquired using a sledge hammer source. Drilling to the basement was performed at five sites. Sediment samples from drilling were analysed for grain-size distribution and age dating. At two drill holes, seismic velocity was measured with depth using a hammer as a seismic source. The geological structure of the study area consists of, from top to bottom, a tidal flat layer (5 ~ 12 m thick), a weathered soil layer (2 ~ 8 m thick), and the basement. The basement is interpreted as Cretaceous sedimentary rocks that tend to be shallow eastward. The volume of the tidal flat sediments and weathered soil in the study area is estimated as $1.4{\times}10^6\;m^3$, weighing $3.5{\times}10^6$ tons. The rate of sea level rise since 8,000 yrs BP is estimated to be 0.1 ~ 0.15 cm/yr.

• Korean Journal of Heritage: History & Science
• /
• v.46 no.4
• /
• pp.4-23
• /
• 2013

Songsalli Ancient Tombs of Gongju consists of seven tombs. King Muryeong's tomb, the seventh tomb, is a brick chamber tomb discovered during the drainage works for the fifth and the sixth tombs in 1971. The excavation at the time focused on topographic surveys of the tomb entrance and the inside of the burial chamber as well as collection of the remains. The burial mount survey confirmed the status of some stone slab remaining and lime-mixed soil layers, but the survey did not examine the exterior structure of the whole tomb as the mounds were removed even more deeply. The excavation revealed damages to the bricks and mural damages due to moisture and fungus in the sixth and the seventh tombs. Between 1996 and 1997, Gongju National University conducted a comprehensive detailed survey of Songsalli Ancient Tombs including a geophysical survey, with an aim to identify the root causes of such degradation. Based on the results, repair took place in 1999 and the fifth, sixth and seventh tombs were placed under permanent conservation to conserve the cultural assets. General public is currently denied access. The purpose of this study was to conduct a three-dimensional resistivity and GPR surveys on the ground surface of the fifth, sixth and seventh tombs of Songsalli Ancient Tombs in order to understand the underground status after repair. The study also aimed to understand the thickness of all the tomb walls and exterior structure based on GPR inside King Muryeong's tomb. The exploration on the ground surface found that the three tombs and soil adjacent to the tombs had resistivity as low as 5 to $90{\Omega}m$, which confirmed that the soil water content was still as high as that prior to the repair work. Additionally, GPR found that the wall construction of the burial chamber of King Muryeong's tomb was approximately 70cm in thickness, while the structure was of 2B with two bricks, about 35cm in length, put together longitudinally(2B brick masonry). The pathway to the burial chamber was of the 2B structure just like that of the burial chamber walls, while its thickness was 80cm with an eyebrow-type arch connected to it. Also, the ceiling exterior appears to have an arch structure, identical to the shape inside.

• Geophysics and Geophysical Exploration
• /
• v.26 no.1
• /
• pp.18-30
• /
• 2023

Submarine mud volcanos are topographic features that resemble volcanoes, and are formed due to eruptions of fluidized or gasified sediment material. They have gained attention as a source of subsurface heat, sediment, or hydrocarbons supplied to the surface. In the continental slope of the Canadian Beaufort Sea, mud volcano exists at various water depths. The MV420, is an active mud volcano erupting at a water depth of 420 meters, and it has been the subject of extensive study. The Korea Polar Research Institute(KOPRI) collected high-resolution seismic data and heat flow data around the caldera of the mud volcano. By analyzing the multi-channel seismic data, we confirmed the reverse-polarity reflector assumed by a gas hydrate-related bottom simulating reflector(BSR). To further elucidate the relationship between the BSR and gas hydrates, as well as the thermal structure of the mud volcano, a numerical geothermal model was developed based on the steady-state heat equation. Using this model, we estimated the base of the gas hydrate stability zone and found that the BSR depth estimated by multi-channel seismic data and the bottom of the gas hydrate stability zone were in good agreement., This suggests the presence of gas hydrates, and it was determined that the depth of the gas hydrate was likely up to 50 m, depending on the distance from the mud conduit. Thus, this depth estimate slightly differs from previous studies.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2010

Seismic reflection surveying is one of the most widely used and effective techniques for coal seam structure delineation and risk mitigation for underground longwall mining. However, the ability of the method can be compromised by the presence of volcanic cover. This problem arises within parts of the Bowen and Sydney Basins of Australia and seismic surveying can be unsuccessful. As a consequence, such areas are less attractive for coal mining. Techniques to improve the success of seismic surveying over basalt flows are needed. In this paper, we use elastic wave-equation-based forward modelling techniques to investigate the effects and characteristics of seismic wave propagation under different settings involving changes in basalt properties, its thickness, lateral extent, relative position to the shot position and various forms of inhomogeneity. The modelling results suggests that: 1) basalts with high impedance contrasts and multiple flows generate strong multiples and weak reflectors; 2) thin basalts have less effect than thick basalts; 3) partial basalt cover has less effect than full basalt cover; 4) low frequency seismic waves (especially at large offsets) have better penetration through the basalt than high frequency waves; and 5) the deeper the coal seams are below basalts of limited extent, the less influence the basalts will have on the wave propagation. In addition to providing insights into the issues that arise when seismic surveying under basalts, these observations suggest that careful management of seismic noise and the acquisition of long-offset seismic data with low-frequency geophones have the potential to improve the seismic results.

• 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.67-72
• /
• 2006

Injection of $CO_2$ into underground saline formations, due to their large storage capacity, is probably the most promising approach for the reduction of $CO_2$ emissions into the atmosphere. $CO_2$ storage must be carefully planned and monitored to ensure that the $CO_2$ is safely retained in the formation for periods of at least thousands of years. Seismic methods, particularly for offshore reservoirs, are the primary tool for monitoring the injection process and distribution of $CO_2$ in the reservoir over time provided that reservoir properties are favourable. Seismic methods are equally essential for the characterisation of a potential trap, determining the reservoir properties, and estimating its capacity. Hence, an assessment of the change in seismic response to $CO_2$ storage needs to be carried out at a very early stage. This must be revisited at later stages, to assess potential changes in seismic response arising from changes in fluid properties or mineral composition that may arise from chemical interactions between the host rock and the $CO_2$. Thus, carefully structured modelling of the seismic response changes caused by injection of $CO_2$ into a reservoir over time helps in the design of a long-term monitoring program. For that purpose we have developed a Graphical User Interface (GUI) driven rock physics simulator, designed to model both short and long-term 4D seismic responses to injected $CO_2$. The application incorporates $CO_2$ phase changes, local pressure and temperature changes. chemical reactions and mineral precipitation. By incorporating anisotropic Gassmann equations into the simulator, the seismic response of faults and fractures reactivated by $CO_2$ can also be predicted. We show field examples (potential $CO_2$ sequestration sites offshore and onshore) where we have tested our rock physics simulator. 4D seismic responses are modelled to help design the monitoring program.

• 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.12 no.1
• /
• pp.33-40
• /
• 2009

In field surveys using the dipole-dipole electrical resistivity method, we often encounter negative apparent resistivity. The term 'negative apparent resistivity' refers to apparent resistivity values with the opposite sign to surrounding data in a pseudosection. Because these negative apparent resistivity values have been regarded as measurement errors, we have discarded the negative apparent resistivity data. Some people have even used negative apparent resistivity data in an inversion process, by taking absolute values of the data. Our field experiments lead us to believe that the main cause for negative apparent resistivity is neither measurement errors nor the influence of self potentials. Furthermore, we also believe that it is not caused by the effects of induced polarization. One possible cause for negative apparent resistivity is the subsurface geological structure. In this study, we provide some numerical examples showing that negative apparent resistivity can arise from geological structures. In numerical examples, we simulate field data using a 3D numerical modelling algorithm, and then extract 2D sections. Our numerical experiments demonstrate that the negative apparent resistivity can be caused by geological structures modelled by U-shaped and crescent-shaped conductive models. Negative apparent resistivity usually occurs when potentials increase with distance from the current electrodes. By plotting the voltage-electrode position curves, we could confirm that when the voltage curves intersect each other, negative apparent resistivity appears. These numerical examples suggest that when we observe negative apparent resistivity in field surveys, we should consider the possibility that the negative apparent resistivity has been caused by geological structure.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.250-262
• /
• 2008

In order to assess the feasibility of the dipole-dipole electric method to the investigation of metallic ore deposit, both field data simulation and inversion are carried out for several simplified ore deposit models. Our interest is in a vein-type model, because most of the ore deposits (more than 70%) exist in a vein type in Korea. Based on the fact that the width of the vein-type ore deposits ranges from tens of centimeters to 2m, we change the width and the material property of the vein, and we use 40m-electrode spacing for our test. For the vein-type model with too small width, the low resistivity zone is not detected, even though the resistivity of the vein amounts to 1/300 of that of the surrounding rock. Considering a wide electrode interval and cell size used in the inversion, it is natural that the size of the low resistivity zone is overestimated. We also perform field data simulation and inversion for a vein-type model with surrounding hydrothermal alteration zones, which is a typical structure in an epithermal ore deposits. In the model, the material properties are assumed on the basis of resistivity values directly observed in a mine originated from an epithermal ore deposits. From this simulation, we can also note that the high resistivity value of the vein does not affect the results when the width of the vein is narrow. This indicates that our main target should be surrounding hydrothermal alteration zones rather than veins in field survey. From these results, we can summarize that when the vein is placed at the deep part and the difference of resistivity values between the vein and the surrounding rock is not large enough, we cannot detect low resistivity zone and interpret the subsurface structures incorrectly using the electric method performed at the surface. Although this work is a little simple, it can be used as references for field survey design and field data Interpretation. If we perform field data simulation and inversion for a number of models and provide some references, they will be helpful in real field survey and interpretation.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.242-249
• /
• 2008

Recently, AVO analysis has been widely used in oil exploration with seismic subsurface section as a direct indicator of the existence of the gas. In the case of the deep reservoirs like the gas reservoirs in the East-sea, it is often difficult to observe AVO responses in CMP gathers even though the bright spots are shown in the stacked section. Because the reservoir becomes more consolidated as its depth deepens, P-wave velocity does not decrease significantly when the pore fluid is replaced by the gas. Thus the difference in Poisson's ratio, which is a key factor for AVO response, between the reservoir and the layer above it does not increase significantly. In this study, we analyzed the effects of Poisson's ratio difference on AVO response with a variety of Poisson's ratios for the upper and lower layers. The results show that, as the difference in Poisson's ratio between the upper and lower layers decreases, the change in the reflection amplitude with incidence angle decreases and AVO responses become insignificant. To consider the limitation of AVO responses shown in the gas reservoir in East-sea, the velocity model was made by simulation Gorae V structure with seismic data and well logs. The results of comparing AVO responses observed from the synthetic data with theoretical AVO responses calculated by using material properties show that the amount of the change in reflection amplitude with increasing incident angle is very small when the difference in Poisson's ratio between the upper and lower layers is small. In addition, the characteristics of AVO responses were concealed by noise or amplitude distortion arisen during preprocessing. To overcome such limitations of AVO analysis of the data from deep reservoirs, we need to acquire precisely reflection amplltudes In data acquisition stage and use processing tools which preserve reflection amplitude in data processing stage.