• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.46-55
• /
• 2002

We investigated the distorting factors of velocity structure reconstructed by traveltime inversion. The set of models that fit the data in a numerical sense usually contains unrealistic models. Reconstructed velocity structure was enhanced because unreasonable models were eliminated by defining constraint of variable grid using a priori information. To correct time delay of source explosion, which distorts traveltime tomograms, terms for correction of time delay was formulated into equation of travel time tomography.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.140-145
• /
• 2011

In order to analyze the influence of the noise on a cross-well traveltime tomography to detect a small scale low velocity body in a homogeneous medium, the first arrival travel times were computed one a tunnel model by a finite-difference ray tracing scheme. Three different types and four different intensity levels of white noises were added to the computed first arrival travel times, and velocity tomograms were constructed using an iterative inversion method (SIRT). Tomograms with the noise intensity up to 10% of the maximum traveltime delay in the tunnel model, showed the exact location of the tunnel. However, the velocity shown at the tunnel location was not close to air velocity but only slightly less than the velocity of the background medium. The additive random noise showed significantly less degree of influence on the resulting tomogram than the source- and receiver consistent noise.

• Journal of the Korean Geophysical Society
• /
• v.8 no.2
• /
• pp.67-74
• /
• 2005

Two-dimensional S-wave velocity sections from SH-wave refraction tomography and surface wave dispersions were obtained by inverting traveltimes of first arrivals and surface wave dispersions, respectively. For the purpose of comparison, a P-wave velocity tomogram was also obtained from a P-wave refraction profiling. P and Rayleigh waves generated by vertical blows on a plate with a sledgehammer were received by 100- and 4.5-Hz geophones, respectively. SH-waves generated by horizontal blows on both sides of a 50 kg timber were received by 8 Hz horizontal geophones. The shear-wave signals were enhanced subtracting data of left-side blows from ones of the right-side blows. Shear-wave velocities from tomography inversion of first-arrival times were compared with ones from inverting dispersion curves of Rayleigh waves. Although the two velocity sections look similar to each other in general, the one from the surface waves tends to have lower velocities. First arrival picking of SH waves is troublesome since P and PS-converted waves arrive earlier than SH waves. Application of the surface wave method, on the other hand, is limited where lateral variation of subsurface tructures is not mild.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.11 s.254
• /
• pp.1051-1056
• /
• 2006

The mean and RMS velocity field of the respiratory gas flow in the human airway was studied experimentally by particle image velocimetry (PIV). Some researchers investigated the airflow for the mouth breathing case both experimentally and numerically. But it is very rare to investigate the airflow of nose breathing in a whole airway due to its geometric complexity. We established the procedure to create a transparent rectangular box containing a model of the human airway for PIV measurement by combination of the RP and the curing of clear silicone. We extend this to make a whole airway including nasal cavities, larynx, trachea, and 2 generations of bronchi. The CBC algorithm with window offset (64 $\times$ 64 to 32 $\times$ 32) is used for vector searching in PIV analysis. The phase averaged mean and RMS velocity distributions in Sagittal and coronal planes are obtained for 7 phases in a respiratory period. Some physiologic conjectures are obtained. The main stream went through the backside of larynx and trachea in inspiration and the frontal side in expiration. There exist vortical motions in inspiration, but no prominent one in expiration.

• Tunnel and Underground Space
• /
• v.18 no.4
• /
• pp.263-271
• /
• 2008

3D seismic tomography technique was assessed for applicability of developed 3D tomography algorithm based on Fresnel volume in the dam-planned area. Reconstructed 3D tomogram based on Fresnel volume and Fast Marching Method(FMM) reveals similar velocity structure to the other geotechnical survey results. With the correlation analysis between RMR data and seismic velocity information, it could provide reliable information of rock mass rate. The applicability of 3D seismic tomography was verified in this study. It would be expected to apply 3D tomography with new developed first arrival calculation and inversion algorithm to the engineering field economically.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.25-30
• /
• 1998

Time-lapse crosswell seismic data, recorded before and after the cavity filling, showed that the filling increased the velocity at a known cavity zone in an old mine site in Inchon area. The seismic response depicted on the tomogram and in conjunction with the geologic data from drillings imply that the size of the cavity may be either small or filled by debris. In this study, I attempted to evaluate the filling effect by analyzing velocity measured from the time-lapse tomograms. The data acquired by a downhole airgun and 24-channel hydrophone system revealed that there exists measurable amounts of source statics. I presented a methodology to estimate the source statics. The procedure for this method is: 1) examine the source firing-time for each source, and remove the effect of irregular firing time, and 2) estimate the residual statics caused by inaccurate source positioning. This proposed multi-step inversion may reduce high frequency numerical noise and enhance the resolution at the zone of interest. The multi-step inversion with different starting models successfully shows the subtle velocity changes at the small cavity zone. The inversion procedure is: 1) conduct an inversion using regular sized cells, and generate an image of gross velocity structure by applying a 2-D median filter on the resulting tomogram, and 2) construct the starting velocity model by modifying the final velocity model from the first phase. The model was modified so that the zone of interest consists of small-sized grids. The final velocity model developed from the baseline survey was as a starting velocity model on the monitor inversion. Since we expected a velocity change only in the cavity zone, in the monitor inversion, we can significantly reduce the number of model parameters by fixing the model out-side the cavity zone equal to the baseline model.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.31-42
• /
• 1998

Recently crosshole seismic tomography has come to be widely used especially for the civil engineering, because it can provide more detail information than any other surface method, although the resolution of tomogram will be inevitably deteriorated to some extent due to the limited wavefield aperture on the nonuniqueness of traveltime inversion. In addition, our field sites often consist of a high-velocity bed rock overlain by low-velocity rock, sometimes with a contrast of more than 45 percent, and furthermore the bed rock is folded. The first arriving waves can be then the refracted ones that travel along the bed rock surface for some source/receiver distances. Thus, the desirable first arrivals can be easily misread that cause severe distortion of the resulting tomogram, if it is concerned with (straight ray) traveltime inversion procedure. In this case, comparision with synthetic data (forward modeling) is a valuable tool in the interpretation process. Besides, abundant information is contained in the crosshole data. For instance, examination of tube waves can be devoted to detecting discontinuities within the borehole such as breakouts, faults, fractures or shear zones as well as the end of the borehole. Specific frequency characteristics of marine silty mud will help discriminate from other soft rocks. The aim of this paper is to present several strategies to analyze and interpret the crosshole data in order to improve the ability at first to determine the spatial dimensions of interwell anomalies and furthermore to understand the underground structures. To this end, our field data are demonstrated. Possibility of misreading the first arrivals was illustrated. Tube waves were investigated in conjunction with the televiewer images. Use of shot- and receiver gathers was examined to benefit the detectabilities of discontinuities within the borehole.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.58-68
• /
• 2011

As part of basic studies of monitoring carbon dioxide ($CO_2$) storage using electrical and seismic surveys, laboratory experiments have been conducted to measure resistivity and P-wave velocity changes due to the injection of $CO_2$ into water-saturated sandstone. The rock sample used is a cylinder of Berea sandstone. $CO_2$ was injected under supercritical conditions (10 MPa, $40^{\circ}C$). The experimental results show that resistivity increases monotonously throughout the injection period, while P-wave velocity and amplitude decrease drastically due to the supercritical $CO_2$ injection. A reconstructed P-wave velocity tomogram clearly images $CO_2$ migration in the sandstone sample. Both resistivity and seismic velocity are useful for monitoring $CO_2$ behaviour. P-wave velocity, however, is less sensitive than resistivity when the $CO_2$ saturation is greater than ~20%. The result indicates that the saturation estimation from resistivity can effectively complement the difficulty of $CO_2$ saturation estimations from seismic velocity variations. By combining resistivity and seismic velocity we were able to estimate $CO_2$ saturation distribution and the injected $CO_2$ behaviour in our sample.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.1479-1483
• /
• 2008

The three main physiological functions of nose are air-conditioning, filtering and smelling. Knowledge of airflow characteristics in nasal cavities is essential to understand the physiological and pathological aspects of nasal breathing. Several studies have utilized physical models of the healthy nasal cavity to investigate the relationship between nasal anatomy and airflow. In our laboratory, there have been a series of experimental investigations on the nasal airflow in normal and deformed nasal cavity models by PIV under both constant and periodic flow conditions. In this time, airflow inside normal nasal cavity is investigated numerically by the FVM general purpose code. The comparisons with PIV measurement are appreciated. Heat and humidity transfer is dealt numerically. Dense CT data and careful treatment of model surface under the ENT doctor’s advice provide more sophisticated cavity models for both PIV experiment and numerical grid system. Average and RMS velocity distributions have been obtained for inspirational and expirational nasal. Temperature distribution, heat and humidity transfer through the mucosa are obtained.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.103-109
• /
• 2005

Airflow in the nasal cavity of a normal Korean adult is investigated experimentally by tomographic PIV measurement. Knowledge of airflow characteristics in nasal cavities is essential to understand the physiology and pathology aspects of nasal breathing. Several studies have utilized physical models of the healthy nasal cavity to investigate the relationship between nasal anatomy and airflow. All of these researches on nasal airflow are under the condition of constant flow-rate. In this study, nasal cavity flow with the physiological period is investigated by tomographic PIV, for the first time. A pumping system that can produce the periodic flow is created. Thanks to a new method for the model casting by a combination of the rapid prototyping and curing of clear silicone, a transparent rectangular box containing the complex nasal cavity can be made for PIV, The CBC PIV algorithm is used for analysis. Phase-averaged mean and RMS velocity distributions are obtained for inspirational and expiration nasal airflows. The comparison with the constant flow case is appreciated. There exist many flow patterns depending on each phase.