• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.67-75
• /
• 1986

This paper introduces the nonlinear Stream Function Wave Theory for design waves efficiently to compute the wave energy and energy transport quantities and to analyze the effects of nonlinearities on them. The Stream Function Wave Theory was developed by Dean for case of the observed waves with assymmetric wave profiles and of the design waves with symmetric theoretical wave profiles. Dalrymple later improved the computational procedure by adding two Lagrangian constraints so that more efficient convergence of the iterative numerical method to a specified wave height and to a zero mean free surface displacement resulted. And the Stream Function coefficients are computed numerically by the improved Marquardt algorithm developed for this study. As the result of this study the effects of nonlinearities on the wave quantities of the average potential energy density, the average kinetic energy density result in overestimation by linear wave theory compared to the Stream Function Wave Theory and increase monotonically with decreasing $L^*/L_O$ and with increasing $H/H_B$. The effects of nonlinearities on the group velocity and the wavelength quantities result in underestimation by linear wave theory and increase monotonically with increasing $H/H_B$. Finally the effect of nonlinearity on the average total energy flux results in overestimation for shallow water waves and underestimation for deep water waves by linear wave theory.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.37B no.10
• /
• pp.889-900
• /
• 2012

In wireless sensor networks, the positioning scheme using received signal strength (RSS) has been widely considered. Appropriate estimation of path-loss exponent (PLE) between a sensor node and an anchor node plays a key role in reducing position error in this RSS-based positioning scheme. In the conventional researches, a sensor node directly uses the PLEs measured by its nearest anchor node to calculate its position. However, the actual PLE between a sensor node and the anchor node can be different from the PLE measured by its nearest anchor node. Thus, if a sensor node directly uses the PLEs measured by its nearest anchor node, the estimated position is different from the actual position of the sensor node with a high probability. In this paper, we describe the method how a sensor node estimates PLEs from the anchor nodes of interest by itself and calculates its position based on these self-estimated PLEs. Especially, our proposal suggests the mechanism to iteratively calculate the PLEs depending on the estimated distances between a sensor node and anchor nodes. Based on the recalculated PLEs, the sensor node reproduces its position. Through simulations, we show that our proposed positioning scheme outperforms the traditional scheme in terms of position error.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.31 no.4
• /
• pp.331-339
• /
• 2013

LiDAR has been one of the widely used and important technologies for 3D modeling of ground surface and objects because of its ability to provide dense and accurate range measurement. The objective of this research is to develop a method for automatic detection and modeling of railroad power lines using high density LiDAR data and RANSAC algorithms. For detecting railroad power lines, multi-echoes properties of laser data and shape knowledge of railroad power lines were employed. Cuboid analysis for detecting seed line segments, tracking lines, connecting and labeling are the main processes. For modeling railroad power lines, iterative RANSAC and least square adjustment were carried out to estimate the lines parameters. The validation of the result is very challenging due to the difficulties in determining the actual references on the ground surface. Standard deviations of 8cm and 5cm for x-y and z coordinates, respectively are satisfactory outcomes. In case of completeness, the result of visual inspection shows that all the lines are detected and modeled well as compare with the original point clouds. The overall processes are fully automated and the methods manage any state of railroad wires efficiently.

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

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

• Journal of Korean Society of Environmental Engineers
• /
• v.27 no.2
• /
• pp.151-157
• /
• 2005

A one-dimensional numerical model was developed to simulate vertical profiles of electron acceptors and their reduced species in benthic sediments. The model accounted for microbial degradation of organic matter and subsequent chemical reactions of interest using stoichiometric relationships. Depending on the dominant electron acceptors utilized by microorganisms, the benthic sediments were assumed to be vertically subdivided into six zones: (1) aerobic respiration, (2) denitrification, (3) manganese reduction, (4) iron reduction, (5) sulfate reduction, and (6) methanogenesis. The utilizations of electron acceptors in the biologically mediated oxidation of organic matter were represented by Monod-type expression. The mass balance equations formulated for the reactive transport of organic matter, electron acceptors, and their corresponding reduced species in the sediments were solved utilizing an iterative multistep numerical method. The ability of model to simulate a freshwater sediments system was tested by comparing simulation results against published data obtained from lake sediments. The simulation results reasonably agreed with field measurements for most species, except for ammonia. This result showed that the C/N ratio (106/16) in the sediments is lower than what the Redfield formula prescribes. Since accurate estimates of vertical profiles of electron acceptors and their reduced species are important to determine the mobility and bioavailability of trace metals in the sediments, the model has potential application to assess the stability of selected trace metals in the sediments.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.184-190
• /
• 1999

A study of gas hydrate is a worldwide popular interesting subject as a potential energy source. A seismic survey for gas hydrate have performed over the East sea by the KIGAM since 1997. General indicators of natural submarine gas hydrates in seismic data is commonly inferred from the BSR (Bottom Simulating Reflection) that occurred parallel to the see floor, amplitude decrease at the top of the BSR, amplitude Blanking at the bottom of the BSR, decrease of the interval velocity, and the reflection phase reversal at the BSR. So the seismic data processing for detecting gas hydrates indicators is required the true amplitude recovery processing, a accurate velocity analysis and the AVO (Amplitude Variation with Offset) analysis. In this paper, we had processed the field data to detect the gas hydrate indicators, which had been acquired over the East sea in 1998. Applied processing modules are spherical divergence, band pass filtering, CDP sorting and accurate velocity analysis. The AVO analysis was excluded, since this field data had too short offset to apply the AVO analysis. The accurate velocity analysis was performed by XVA (X-window based Velocity Analysis). This is the method which calculate the velocity spectrum by iterative and interactive. With XVA, we could determine accurate stacking velocity. Geobit 2.9.5 developed by the KIGAM was used for processing data. Processing results say that the BSR occurred parallel to the sea floor were shown at $367\~477m$ depths (two way travel time about 1800 ms) from the sea floor through shot point 1650-1900, the interval velocity decrease around BSR and the reflection phase reversal corresponding to the reflection at the sea floor.

• Geophysics and Geophysical Exploration
• /
• v.18 no.1
• /
• pp.1-8
• /
• 2015

To estimate the shear-wave velocity (${\nu}_s$ beneath the OTVZ seismic station on Ngauruhoe volcano in New Zealand, we calculated receiver functions (RFs) using 127 teleseismic data ($Mw{\geq}5.5$) with high signal-to-noise ratios recorded during November 11, 2011 to September 11, 2013. The genetic inversion algorithms was applied to 21 RFs calculated by the iterative time-domain deconvolution method. In the 1-D ${\nu}_s$ model derived by the inversion, the Moho is observed at a 14 km depth, marked by a ${\nu}_s$ transition from 3.7 km/s to 4.7 km/s. The average ${\nu}_s$ of the overlying crust is 3.4 km/s, and the average ${\nu}_s$ of a greater than 9-km thick low-velocity layer (LVL) in the lower crust is 3.1 km/s. The LVL becomes thinner with increasing distance from the station. Another LVL thicker than 10 km with ${\nu}_s$ less than 4.3 km/s is found in the upper mantle. Those LVLs in the lower crust and the upper mantle and the relatively thin crust might be related to the magma activity caused by the subducting Pacific plate.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.19 no.5
• /
• pp.552-560
• /
• 2013

In this paper, an optimization algorithm for the block assembly accuracy control assessment is proposed with consideration for the current block assembly process and accuracy control procedure used in the shipbuilding site. The objective function of the proposed algorithm consists of root mean square error of the distances between design and measured data of the other control points with respect to a specific point of the whole control points. The control points are divided into two groups: points on the control line and the other points. The grouped data are used as criteria for determining the combination of 6 degrees of freedom in the registration process when constituting constraints and calculating objective function. The optimization algorithm is developed by using combination of the sampling method and the point to point relation based modified ICP algorithm which has an allowable error check procedure that makes sure that error between design and measured point is under allowable error. According to the results from the application of the proposed algorithm with the design and measured data of two blocks data which are verified and validated by an expert in the shipbuilding site, it implies that the choice of whole control points as target points for the accuracy calculation shows better results than that of the control points on the control line as target points for the accuracy of the calculation and the best optimized result can be acquired from the accuracy calculation with a fixed point on the control line as the reference point of the registration.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.66-73
• /
• 2014

To estimate the shear-velocity ($v_s$) structure beneath the WIZ station on White Island in New Zealand, we applied receiver function (RF) inversion and H-${\kappa}$ stacking methods to 362 teleseismic events (Mw > 5.5) recorded during April 20, 2007 to September 6, 2013. Using 71 RFs with errors less than 20% after 200 iterative computations, we determined that the depth to Moho of $v_s$ = 4.35 km/s is $24{\pm}1km$ within a 15 km radius of the station. In an 1-d $v_s$ model derived by RF inversions, a 4-km thick low-velocity layer (LVL) at depths of 18 ~ 22 km was identified in the lower crust. This LVL, which is 0.15 km/s slower than the rocks above and below it, may indicate the presence of a deep magma reservoir. The H-${\kappa}$ stacking method yielded an estimate of the depth to the Moho of 24.5 km, which agrees well with the depth determined by RF inversions. The low $v_p/v_s$ ratio of 1.64 may be due to the presence of gas-filled rock or hot crystallizing magma.

• Journal of the Korean earth science society
• /
• v.30 no.2
• /
• pp.153-164
• /
• 2009

We examined the regional 1-D deep resistivity structure of the Korean Peninsula using MT data acquired at seven sites located in the Kyongsang Basin and Kyonggi Massif. At the sites located in the Kyongsang Basin, surrounding sea distorts observed MT response and hence this distortion, so called "sea effect", is corrected using an iterative tensor stripping method. The 1-D layered inversion results for the seven MT sites reveal 4 layered structure, which is composed of 1) near surface layer, 2) upper crust, 3) lower crust and upper mantle, and 4) asthenosphere from the surface downward. Conrad interface, which is a boundary between upper and lower crust, is distinctly identified beneath all the MT sites. Conrad interface depth is estimated to about be 17km in the Kyongsang Basin and about 12km in the Kyonggi Massif, while the upper crust of the Kyongsang Basin is about 5 times more resistive than that of the Kyonggi Massif. Finally, asthenosphere is inferred to exist below a depth of approximately 100km with a resistivity of 200-300 ohm-m.