• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.85-92
• /
• 2008

National Oceanographic Research Institute is carrying out an oceanographic survey for the entire sea areas around Korean Peninsula annually starting with the East Sea from 1996 by establishing a national oceanographic basic map survey plan for the sea areas under the jurisdiction of Korea, so this paper used the oceanographic geomagnetism data measured at the southern area of the Yellow Sea using 'Hae Yang 2000' in 1999, aiming at clarifying the cause of geomagnetic abnormality zone during the course of treating and analyzing the geomagnetic data. For treatment of magnetic data, we obtained electromagnetic force values and geomagnetic abnormality values around the investigated sea area through a process of searching and removal of bad data, correction of sensor positions, correction of magnetic field effects around the hull, correction of diurnal variation, normal correction, correction of cross point errors, etc. The electromagnetic force distribution around the investigated sea area was $49000\;{\sim}\;51600\;nT$, which is judged to be within the normal electromagnetic force intensity distribution range around the Yellow Sea. The isodynamic lines are distributed in Northeast-Southwest direction, and electromagnetic force values are increasing toward the northwest. The result of comparing the magnetic abnormality around the sea area among $124^{\circ}$ 49' 48" E, $35^{\circ}$ 10' 48" N $\sim$ $125^{\circ}$ 7' 48" E, and $35^{\circ}$ 33' 00" N sections with the elastic wave cross section and the result of modeling coincide well with the underground geological structure clarified from the existing elastic wave survey cross section. Therefore, it is judged that the distribution of magnetic force abnormality generally shows the effect pursuant to the distribution of the sedimentary basins in the Tertiary period and the bedrocks in the Cretaceous period which are well developed in the bottom of the sea.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.6
• /
• pp.440-447
• /
• 2013

In this paper, a higher harmonic control (HHC) methodology is applied to find the optimum input scenario for the vibratory hub loads reduction. A comprehensive aeroelastic analysis code, CAMRAD II, is used to model the HART (Higher-harmonic-control Aeroacoustic Rotor Test) II rotor, and parametric study is conducted for the best HHC inputs leading to a minimum vibration (MV) condition. The resulting outcomes are compared with the earlier HART II test results. It is indicated that the control input adopted in the MV condition showed less satisfactory results. The new MV condition obtained in the present investigation can achieve 45% lower vibration level than the baseline uncontrolled condition. The optimum HHC input results lead to 3/rev harmonic input having $0.8^{\circ}$ amplitude and $350^{\circ}$ phase angle. About 5% reduction in the required power is possible but accompanies with the increase of vibration level.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.131-135
• /
• 2016

Coarray is a parallel processing technique introduced in the Fortran 2008 standard. Coarray can implement parallel processing using simple syntax. In this research, we examined applicability of Coarray to seismic parallel processing by comparing performance of seismic data processing programs using Coarray and MPI. We compared calculation time using seismic wave propagation modeling and one to one communication time using domain decomposition technique. We also compared performance of parallel reverse-time migration programs using Coarray and MPI. Test results show that the computing speed of Coarray method is similar to that of MPI. On the other hand, MPI has superior communication speed to that of Coarray.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.173-178
• /
• 2008

In this paper, we will introduce rock physics modeling technique, which interrelate reservoir properties with seismic properties, and apply the technique to the Donghae-1 gas reservoir. From well-log data analysis, we obtained velocityporosity (Vp-$\phi$) relations for each formation. These relations can used to predict porosity from seismic data. In addition, we analyzed permeability data, which were obtained from core measurements and computational rock physics simulations. We then obtained permeability-porosity ($\kappa-\phi$) relations. Combining $\kappa-\phi$ with Vp-$\phi$ relations, we finally present quantitative Vp-$\kappa$ relations. As to Vp-$\phi$ modeling, we found that the degree of diagenesis and clay contents increase with depth. As to Vp-$\kappa$ relations, though \kappa-\phi relations are almost identical for all formations, we could obtain distinct Vp-$\kappa$ relations due to Vp-$\phi$ variations. In conclusion, the rock physics modeling, which bridges between seismic properties and reservoir properties, can be a very robust tool for quantitative reservoir characterization with less uncertainty.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.7-17
• /
• 2005

The determination of seismic velocities in refractors for near-surface seismic refraction investigations is an ill-posed problem. Small variations in the computed time parameters can result in quite large lateral variations in the derived velocities, which are often artefacts of the inversion algorithms. Such artefacts are usually not recognized or corrected with forward modelling. Therefore, if detailed refractor models are sought with model based inversion, then detailed starting models are required. The usual source of artefacts in seismic velocities is irregular refractors. Under most circumstances, the variable migration of the generalized reciprocal method (GRM) is able to accommodate irregular interfaces and generate detailed starting models of the refractor. However, where the very-near-surface environment of the Earth is also irregular, the efficacy of the GRM is reduced, and weathering corrections can be necessary. Standard methods for correcting for surface irregularities are usually not practical where the very-near-surface irregularities are of limited lateral extent. In such circumstances, the GRM smoothing statics method (SSM) is a simple and robust approach, which can facilitate more-accurate estimates of refractor velocities. The GRM SSM generates a smoothing 'statics' correction by subtracting an average of the time-depths computed with a range of XY values from the time-depths computed with a zero XY value (where the XY value is the separation between the receivers used to compute the time-depth). The time-depths to the deeper target refractors do not vary greatly with varying XY values, and therefore an average is much the same as the optimum value. However, the time-depths for the very-near-surface irregularities migrate laterally with increasing XY values and they are substantially reduced with the averaging process. As a result, the time-depth profile averaged over a range of XY values is effectively corrected for the near-surface irregularities. In addition, the time-depths computed with a Bero XY value are the sum of both the near-surface effects and the time-depths to the target refractor. Therefore, their subtraction generates an approximate 'statics' correction, which in turn, is subtracted from the traveltimes The GRM SSM is essentially a smoothing procedure, rather than a deterministic weathering correction approach, and it is most effective with near-surface irregularities of quite limited lateral extent. Model and case studies demonstrate that the GRM SSM substantially improves the reliability in determining detailed seismic velocities in irregular refractors.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.21 no.3
• /
• pp.269-276
• /
• 2001

Ultrasonic testing of composite materials is much more difficult than that of isotropic materials, because of the beam skew phenomenon caused by their elastic anisotropy. An established analytical method exists for elastic wave propagation in anisotropic media as a result of previous research efforts. Yet, due to the complexity of the analytical method, solution of real problems must resort to the numerical method. In this work, analytical solutions have first been obtained for the wavefield due to a point source in a unidirectional fiber-reinforced composite, which may be modeled as transversely isotropic. Then, the corresponding numerical solutions have been obtained using the mass-spring lattice model(MSLM). The two solutions have agreed well with each other. Other problems such as reflection from free boundaries and scattering from cracks have also been solved numerically, and the results have been investigated from the viewpoint of wave mechanics. The numerical model whose validity has been confirmed by this work will be of great use in simulating ultrasonic testing of composite materials.

• 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.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Composites Research
• /
• v.18 no.5
• /
• pp.15-20
• /
• 2005

In this paper, damage induced acoustic emission in the composite plate in numerically simulated by using the three dimensional finite element method and explicit time integration. Acoustic source is modeled by equivalent volume source. To verify the proposed method, dynamic displacements due to the elastic wave are compared with the experiment when the fiber is broken in the single fiber embedded isotropic plate. For the laminated composite plates, the results are compared between homogenized model and DNS approach which models fibers and matrix separately. To capture high frequencies in the elastic wave, small time step size and a large number of meshes are required. The parallel computing technology is introduced to solve a large scale problem efficiently.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.203-208
• /
• 2007

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms for the travel time computation commonly used, however, may not give sufficiently precise results particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. We assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution.