• The Journal of Engineering Geology
• /
• v.30 no.3
• /
• pp.289-302
• /
• 2020

To determine the shallow subsurface structure and sliding surface of land creeping in 2016 at Hadong-gun, Gyeongsangnam-do, geophysical surveys (electric resistivity, and refraction seismic methods, borehole televiewer) and slope stability analysis were conducted. The subsurface structure delineated with borehole lithologies and seismic velocity structures provided the information that the sediment layer on the top of the slope was rather as thick as 20 m and the underlying weathered rock (anorthosite) was thinner than 1 m. Based on the tension cracks observed during the geological mapping, televiewer scanning was performed at the borehole BH-2 and detected the intensive fracture zones at the ground-water level, associated with the slip weak zones mapped in dipole-dipole electrical resistivity section. Downslope sliding and slightly upward pushing at the apex of high resistive bedrock explains the curved slip plane of the land creeping. Such a convex structure might play a role of natural toe abutment for preventing the downward development of slip weak zones. In slope stability analysis, the safety factors of the slip weak zone are calculated with varying the groundwater levels for dry and rainy seasons and the downslope is founded to be unstable with safety factor of 0.89 due to fully saturated material in rainy season.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.369-376
• /
• 2007

Electrical resistivity tomography was carried out at seawater intrusion monitoring wells located at watershed in coastal areas. It is difficult to identify the characteristics of resistivity near monitoring well in case of using high signalto-noise ratio array due to the high conductivity condition in coastal aquifer although electrical resistivity survey is well adopted to delineate hydrogeological characteristics with the distribution of electrical resistivity. To improve the quality of electrical resistivity survey for two sites with seawater intrusion monitoring wells, inversion with the results of holeto-surface electrical resistivity tomography using single well was executed. The results of inversion for aquifer near wells were verified with the results of drilling log with the informations of fracture, electrical conductivity logging and normal resistivity logging. The inversion for aquifer near one of two wells was also performed at low and high tide with the same electrodes, respectively. From the inversion result, it is possible to obtain the resistivity images with high resolution and to identify the characteristics of aquifer related to seawater intrusion with tidal fluctuation. From this study, it was demonstrated that the hole-to-surface electrical resistivity tomography method accompanied with drilling log, electrical conductivity logging and normal resistivity logging would be useful to delineate the hydrogeological structures near monitoring wells in coastal areas.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.223-235
• /
• 2002

Since weak Bones or geological lineaments are likely to be eroded, there may develop weak Bones beneath rivers, and a careful evaluation of ground condition is important to construct structures passing through a river. DC resistivity method, however, has seldomly applied to the investigation of water-covered area, possibly because of difficulties in data aquisition and interpretation. The data aquisition having high quality may be the most important factor, and is more difficult than that in land survey, due to the water layer overlying the underground structure to be imaged. Through the numerical modeling and the analysis of a case history, we studied the method of resistivity survey at the water-covered area, starting from the characteristics of measured data, via data acquisition method, to the interpretation method. We unfolded our discussion according to the installed locations of electrodes, ie., floating them on the water surface, and installing them at the water bottom, because the methods of data acquisition and interpretation vary depending on the electrode location. Through this study, we could confirm that the DC resistivity method can provide fairly reasonable subsurface images. It was also shown that installing electrodes at the water bottom can give the subsurface image with much higher resolution than floating them on the water surface. Since the data acquired at the water-covered area have much lower sensitivity to the underground structure than those at the land, and can be contaminated by the higher noise, such as streaming potential, it would be very important to select the acquisition method and electrode array being able to provide the higher signal-to-noise ratio (S/N ratio) data as well as the high resolving power. Some of the modified electrode arrays can provide the data having reasonably high S/N ratio and need not to install remote electrode(s), and thus, they may be suitable to the resistivity survey at the water-covered area.

• Tunnel and Underground Space
• /
• v.15 no.4 s.57
• /
• pp.275-287
• /
• 2005

At the planning and design stages of a development of underground space or tunneling project, the information regarding ground conditions is very important to enhance economical efficiency and overall safety In general, the information can be expressed using RMR or Q-system and with the geophysical exploration image. RMR or Q-system can provide direct information of rock mass in a local scale for the design scheme. Oppositely, the image of geophysical exploration can provide an exthaustive but indirect information. These two types of the information have inherent uncertainties from various sources and are given in different scales and with their own physical meanings. Recently, RMR has been estimated in unsampled areas based on given data using geostatistical methods like Kriging and conditional simulation. In this study, simulated annealing(SA) is applied to overcome the shortcomings of Kriging methods or conditional simulations just using a primary variable. Using this technique, RMR and the image of geophysical exploration can be integrated to construct the spatial distribution of RM and to evaluate its uncertainty. The SA method was applied to solve an optimization problem with constraints. We have suggested the practical procedure of the SA technique for the uncertainty evaluation of RMR and also demonstrated this technique through an application, where it was used to identify the spatial distribution of RMR and quantify the uncertainty. For a geotechnical application, the objective functions of SA are defined using statistical models of RMR and the correlations between RMR and the reference image. The applicability and validity of this application are examined and then the result of uncertainty evaluation can be used to optimize the tunnel layout.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2002.09a
• /
• pp.101-125
• /
• 2002

A lake seismic survey was carried out to investigate possible geohazards for construction of the underground LPG storage at Namyang Lake. The proposed survey site has a land-lake combined geography and furthermore water depth of the lake is shallow. Therefore, various seismic methods such as marine single channel high resolution seismic reflection survey, sonobuoy refraction survey, land refraction survey and land-lake combined refraction survey were applied. Total survey amounts are 34 line-km of high resolution lake seismic survey, 14 lines of sonobuoy refraction survey, 890 m of land refraction survey and 8 lines of land-lake combined refraction survey. During the reflection survey, there were severe water reverberations from the lake bottom obscured subsurface profiling. These strong multiple events appeared in most of the survey area except the northern and southern area near the embankment where seems to be accumulated mainly mud dominated depositions. The sonobuoy refraction profiles also showed the same Phenomena as those of reflection survey. Meanwhile the results of the land-lake combined refraction survey showed relatively better qualities. However, the land refraction survey did not so due to low velocity soil layer and electrical noise. Summarized results from the lake seismic survey are that acoustic basement with relatively flat pattern appeared 30m below water level and showed three types of bedrock such as fresh, moderately weathered and weathered type. According to the results of the combined refraction survey, a velocity distribution pattern of the lake bottom shows three types of seismic velocity zone such as >4.5 km/s, 4.5-4.0km/s and <4.0km/s. The major fault lineament in the area showed NW-SE trend which was different from the Landsat image interpretation. A drilling was confirmed estimated faults by seismic survey.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.84-88
• /
• 2001

For the elastic migation, the velocity errors between the initial velocity model and true velocity model seriously affect the migrated images. The assumption of an initial velocity model, thus, is one of the critical factor for the successful migration. In case of applying the layered earth model as an initial velocity model, the layer boundary having large velocity contrast can not be defined well with conventional traveltime calculation algolithms and we have the difficulties for expressing the characteristics of the real subsurface. Smooth Background Model (SBM) we have applied as an initial velocity model in our study is characterized to be linearly varying the velocity with the depth, which can express the velocity variation in the subsurface properly. Thus it can properly be applied to traveltime calculation algolithms such as Vidale's method. In this study, Kirchhoff operator for prestack migration was used and the absolute amplitude obtained by modeling was applied as a weighted value to consider the true amplitude for initial model. Initial velocity model for migration was determined by using stacking velocity and we applied this model to real data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.44 no.1
• /
• pp.49-62
• /
• 2024

Reckless development of the underground by rapid urbanization causes inspection delay on replacement of existing structure and installation new facilities. However, frequent accidents occur due to deviation in construction design planned by inaccurate location information of underground structure. Meanwhile, the electrical resistivity survey, knowns as non-destructive method, is based on the difference in the electric potential of electrodes to measure the electrical resistance of ground. This method is significantly advanced with multi-electrode and deep learning for analyzing strata. However, there is no study to quantitatively assess change in electrical resistance according to geometric conditions of structures. This study evaluates changes in electrical resistance through geometric parameters of electrodes and structure. Firstly, electrical resistance numerical module is developed using generalized mesh occurring minimal errors between theoretical and numerical resistance values. Then, changes in resistances are quantitatively compared on geometric parameters including burial depth, diameter of structure, and distance electrode and structure under steady current condition. The results show that higher electrical resistance is measured for shallow depth, larger size, and proximity to the electrode. Additionally, electric potential and current density distributions are analyzed to discuss the measured electrical resistance around the terminal electrode and structure.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.199-206
• /
• 2003

Characteristics of the static shift are discussed by comparing the three-dimensional MT inversion with/without static shift parameterization. The galvanic distortion by small-scale shallow feature often leads severe distortion in inverted resistivity structures. The new inversion algorithm is applied to four numerical data sets contaminated by different amount of static shift. In real field data interpretations, we generally do not have any a-priori information about how much the data contains the static shift. In this study, we developed an algorithm for finding both Lagrangian multiplier for smoothness and the trade-off parameter for static shift, simultaneously in 3-D MT inversion. Applications of this inversion routine for the numerical data sets showed quite reasonable estimation of static shift parameters without any a-priori information. The inversion scheme is successfully applied to all the four data sets, even when the static shift does not obey the Gaussian distribution. Allowing the static shift parameters have non-zero degree of freedom to the inversion, we could get more accurate block resistivities as well as static shifts in the data. When inversion does not consider the static shift as inversion parameters (conventional MT inversion), the block resistivities on the surface are modified considerably to match possible static shift. The inhomogeneous blocks on the surface can generate the static shift at low frequencies. By those mechanisms, the conventional 3-D MT inversion can reconstruct the resistivity structures to some extent in the deeper parts even when moderate static shifts are in the data. As frequency increased, however, the galvanic distortion is not frequency independent any more, and thus the conventional inversion failed to fit the apparent resistivity and phase, especially when strong static shift is added. Even in such case, however, reasonable estimation of block resistivity as well as static shift parameters were obtained by 3-D MT inversion with static shift parameterization.

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

A pulse radar system has been developed recently to detect dormant underground tunnels that are deeply located at depths of hundreds of metres. To check the ability of the radar system to detect an obliquely oriented tunnel, five different borehole pairs in the tunnel test site were chosen so that the horizontal lines-of-sight cut the tunnel axis obliquely, in $15^{\circ}$ steps. The pulse radar signatures were measured over a depth range of 20 m around the centre of the air-filled tunnel. Three canonical parameters, consisting of the arrival time, attenuation, and dispersion time were extracted from the first and second peaks of the measured radar signatures. Using those parameters, the radar system can detect obliquely oriented tunnels at various angles up to 45 from the transmitter-receiver line of sight.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.34-44
• /
• 2001

Forward modeling by construction of synthetic data is usually practiced in a horizontal surface and a few subsurface structures. However, in-situ surveys often take place in such topographic changes that the corrupted field data always make it difficult to interpret the right signals. To examine the propagation characteristic of elastic waves on the irregular surface, a general mesh generation code for finite element method was modified to consider the topography. By implementing this algorithm, the time domain modeling was practiced in some models with surface topography such as mound, channel, etc. The synthetic data obtained by receivers placed on surface also agreed with the analytic solution. The snapshots showing the total wave-field revealed the propagation characteristic of the elastic waves through complex subsurface structures and helped to identify the signals on the time traces. The transmission of Rayleigh waves along the surface, compressive waves, and sheer waves was observed. Moreover, it turned out that the Rayleigh waves behave like a new source at the edge.