• Journal of the korean society of oceanography
• /
• v.38 no.4
• /
• pp.166-172
• /
• 2003

Multi­Sensor Core Logger (MSCL) is a useful system for logging the physical properties (compressional wave velocity, wet bulk density, fractional porosity, magnetic susceptibility and/or natural gamma radiation) of marine sediments through scanning of whole cores in a nondestructive fashion. But MSCL has a number of problems that can lead to spurious results depending on the various factors such as core slumping, gas expansion, mechanical stretching, and the thickness variation of core liner and sediment. For the verification of MSCL data, compressional wave velocity, wet bulk density, and porosity were measured on discrete samples by Hamilton Frame and Gravimetric method, respectively. Acoustic impedance was also calculated. Physical property data (velocity, wet bulk density, and impedance) logged by MSCL were slightly larger than those of discrete sample, and porosity is reverse. Average difference between MSCL and discrete sample at both sites is relatively small such as 22­24 m/s in velocity, $0.02­-0.08\;g/\textrm{cm}^3$ in wet bulk density, and 2.5­2.7% in porosity. The values also show systematic variation with sediment depth. A variety of factors are probably responsible for the differences including instrument error, various measurement method, sediment disturbance, and accuracy of calibration. Therefore, MSCL can be effectively used to collect physical property data with high resolution and quality, if the calibration is accurately completed.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.12
• /
• pp.127-139
• /
• 2017

Site investigation including boring and various in-situ borehole test (Pressuremeter test, Borehole shear test, Downhole test, Suspension PS logging, Density logging) and X-ray fluorescence analysis for rock core sample were performed to estimate geotechnical properties and weathering degree of weathered granite rock in Goyang. Deformation modulus, shear strength parameter and shear wave velocity estimated through in-situ borehole test had a tendency to increase with depth. And several chemical weathering indices evaluated by X-ray fluorescence analysis had a general tendency of reducing weathering degree in accordance with depth. Also, relationship between VR determined as a representative weathering index and the geotechnical properties was analyzed.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.09a
• /
• pp.446-449
• /
• 2004

High resolution geophysical imaging to delineate costal aquifer and seawater- freshwater interface has been applied in Baesu-eup, Yeonggwang-gun, Jeolla province Electrical resistivity information from vertical electrical sounding and 2-D electrical resistivity survey is key parameter to map equivalent Nacl concentration map over the survey area. Seismic velocity from refraction tomographic survey, on the other hand, gives more reliable information on the subsurface stratagraphy than electrical resistivity methods which frequently suffer from low resolution due to masking effect. We imaged high-resolution 3-D structure of costal aquifer by correlating the electrical resistivity with seismic velocity, and mapped equivalent NaCl concentration map using resistivity and hydro-geological information from well logging.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.19-29
• /
• 2006

The first Japanese pilot-scale $CO_2$ sequestration project has been undertaken in an onshore saline aquifer, near Nagaoka in Niigata prefecture, and time-lapse well logs were carried out in observation wells to detect the arrival of injected $CO_2$ and to evaluate $CO_2$ saturation in the reservoir. $CO_2$ was injected into a thin permeable zone at the depth of 1110m at a rate of 20-40 tonnes per day. The total amount of injected $CO_2$ was 10400 tonnes, during the injection period from July 2003 to January 2005. The pilot-scale demonstration allowed an improved understanding of the $CO_2$ movement in a porous sandstone reservoir, by conducting time-lapse geophysical well logs at three observation wells. Comparison between neutron well logging before and after the insertion of fibreglass casing in observation well OB-2 showed good agreement within the target formation, and the higher concentration of shale volume in the reservoir results in a bigger difference between the two well logging results. $CO_2$ breakthrough was identified by induction, sonic, and neutron logs. By sonic logging, we confirmed P-wave velocity reduction that agreed fairly well with a laboratory measurement on drilled core samples from the Nagaoka site. We successfully matched the history changes of sonic P-wave velocity and estimated $CO_2$ saturation a(ter breakthrough in two observation wells out of three. The sonic-velocity history matching result suggested that the sweep efficiency was about 40%. Small effects of $CO_2$ saturation on resistivity resulted in small changes in induction logs when the reservoir was partially saturated. We also found that $CO_2$ saturation in the $CO_2$-bearing zone responded to suspension of $CO_2$ injection.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.457-465
• /
• 2013

We integrated and correlated datasets from surface and subsurface geophysics, drilling cores, and engineering geology to identify geological interfaces and characterize the joints and fracture zones within the rock mass. The regional geometry of a geologically weak zone was investigated via a fence projection of electrical resistivity data and a borehole image-processing system. Subsurface discontinuities and intensive fracture zones within the rock mass are delineated by cross-hole seismic tomography and analyses of dip directions in rose diagrams. The dynamic elastic modulus is studied in terms of the P-wave velocity and Poisson's ratio. Subsurface discontinuities, which are conventionally identified using the N value and from core samples, can now be identified from anomalous reflection coefficients (i.e., acoustic impedance contrast) calculated using a pair of well logs, comprising seismic velocity from suspension-PS logging and density from logging. Intensive fracture zones identified in the synthetic seismogram are matched to core loss zones in the drilling core data and to a high concentration of joints in the borehole imaging system. The upper boundaries of fracture zones are correlated to strongly negative amplitude in the synthetic trace, which is constructed by convolution of the optimal Ricker wavelet with a reflection coefficient. The standard deviations of dynamic elastic moduli are higher for fracture zones than for acompact rock mass, due to the wide range of velocities resulting from the large numbers of joints and fractures within the zone.

• Economic and Environmental Geology
• /
• v.43 no.6
• /
• pp.615-623
• /
• 2010

Seismic refraction survey is a geophysical method that delineates subsurface velocity structure using direct wave and critically refracted wave. The generalized reciprocal method(GRM) is an inversion technique which uses travel-time data from several forward and reverse shots and which can provide the geometry of irregular inclined refractors and structures underlain by hidden layer such as low velocity zone and thin layer. In this study, a simple Excel-GRM routine was tested for fast mapping of the interface between weathering layer and bedrock during the survey, with employing a pair of forward and reverse shots. This routine was proved to control the maximum dip of approximately $30^{\circ}C$ and maximum velocity contrast of 0.6, based on the panel tests in terms of dipping angle and velocity contrast for the two-layer inclined models. In contrast with conventional operation of five to seven shots with sufficient offset distance and indoor data analysis thereafter, this routine was performed in the field shortly after data acquisition. Depth to the bedrock provided by Excel-GRM, during the field survey for Cheongju granite area, correlates well with the elevation of the surface of soft rock from the drill core and SPS logging data. This cost-effective routine developed for quickly delineating the bedrock surface in the field survey will be readily applicable to mapping of weathering zone in narrow zone with small variation of elevation of bedrock.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.263-270
• /
• 2008

The Gas Hydrate Research and Development Organization (KGHDO) of Korea accomplished successfully geophysical logging (LWD: Logging While Drilling, MWD: Measurement While Drilling) for five sites in 2007, in order to investigate the presence of gas hydrate in the Ulleung Basin, the East Sea of Korea. The togging parameters acquired from LWD/MWD dre electrical resistivity, acoustic velocity, neutron density and porosity, and natural gamma. In addition, pressure, temperature, and diameter of borehole were measured. LWD/MWD data showed several evidences indicating the presence of gas hydrate. Based on LWD/MWD data, three coring sites were selected for sampling of gas hydrate. Subsequently, various gas hydrate samples were collected directly from three sites. Therefore. the presence of gas hydrates was verified by coring. LWD/MWD data will be significantly used to estimate the amount of gas hydrate. Also, they will provide important information to elucidate about sedimentologic characteristics of gas-hydrate bearing formation and sedimentary environment of the Ulleung Basin.

• Geophysics and Geophysical Exploration
• /
• v.8 no.4
• /
• pp.237-242
• /
• 2005

A seismic refraction study in estimation of depth to the bedrock demonstrates that 1) the average velocity in the medium is about 250 m/s in the surface layer (< 4 m), 2,500 m/s in the weathered formation, and greater than 3,000 m/s in the bedrock, 2) the depth to the deepest reflector assumed to be the bedrock is about 17 m; however, according to the cores collected in a borehole in study area, the bedrock (granite) occurred at depth 25 m, 3) according to the density and velocity logging, at depth 17 m, a measurable velocity and density increase are observed, and 4) the velocity of the weathered formation is relatively high and therefore, the acquisition offsets ($70{\sim}80m$) are turned out not to be long enough to record the refracted signal from the bedrock at depth 25 m as first arrivals.

• 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.22 no.4
• /
• pp.225-238
• /
• 2019

Migration velocity analysis (MVA) for creating optimum depth-domain velocities in seismic imaging was applied to marine long-offset multi-channel data, and the effectiveness of the MVA approach was demonstrated by the combinations of conventional data processing procedures. The time-domain images generated by conventional time-processing scheme has been considered to be sufficient so far for the seismic stratigraphic interpretation. However, when the purpose of the seismic imaging moves to the hydrocarbon exploration, especially in the geologic modeling of the oil and gas play or lead area, drilling prognosis, in-place hydrocarbon volume estimation, the seismic images should be converted into depth domain or depth processing should be applied in the processing phase. CMP-based velocity analysis, which is mainly based on several approximations in the data domain, inherently contains errors and thus has high uncertainties. On the other hand, the MVA provides efficient and somewhat real-scale (in depth) images even if there are no logging data available. In this study, marine long-offset multi-channel seismic data were optimally processed in time domain to establish the most qualified dataset for the usage of the iterative MVA. Then, the depth-domain velocity profile was updated several times and the final velocity-in-depth was used for generating depth images (CRP gather and stack) and compared with the images obtained from the velocity-in-time. From the results, we were able to confirm the depth-domain results are more reasonable than the time-domain results. The spurious local minima, which can be occurred during the implementation of full waveform inversion, can be reduced when the result of MVA is used as an initial velocity model.