• Geophysics and Geophysical Exploration
• /
• v.23 no.3
• /
• pp.168-177
• /
• 2020

Seabed multiple reflections (seabed multiples) are the main cause of misinterpretations of primary reflections in both shot gathers and stack sections. Accordingly, seabed multiples need to be suppressed throughout data processing. Conventional model-driven methods, such as prediction-error deconvolution, Radon filtering, and data-driven methods, such as the surface-related multiple elimination technique, have been used to attenuate multiple reflections. However, the vast majority of processing workflows require time-consuming steps when testing and selecting the processing parameters in addition to computational power and skilled data-processing techniques. To attenuate seabed multiples in seismic reflection data, input gathers with seabed multiples and label gathers without seabed multiples were generated via numerical modeling using the Marmousi2 velocity structure. The training data consisted of normal-moveout-corrected common midpoint gathers fed into a U-Net neural network. The well-trained model was found to effectively attenuate the seabed multiples according to the image similarity between the prediction result and the target data, and demonstrated good applicability to field data.

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

Seismic reflection surveying is one of the most widely used and effective techniques for coal seam structure delineation and risk mitigation for underground longwall mining. However, the ability of the method can be compromised by the presence of volcanic cover. This problem arises within parts of the Bowen and Sydney Basins of Australia and seismic surveying can be unsuccessful. As a consequence, such areas are less attractive for coal mining. Techniques to improve the success of seismic surveying over basalt flows are needed. In this paper, we use elastic wave-equation-based forward modelling techniques to investigate the effects and characteristics of seismic wave propagation under different settings involving changes in basalt properties, its thickness, lateral extent, relative position to the shot position and various forms of inhomogeneity. The modelling results suggests that: 1) basalts with high impedance contrasts and multiple flows generate strong multiples and weak reflectors; 2) thin basalts have less effect than thick basalts; 3) partial basalt cover has less effect than full basalt cover; 4) low frequency seismic waves (especially at large offsets) have better penetration through the basalt than high frequency waves; and 5) the deeper the coal seams are below basalts of limited extent, the less influence the basalts will have on the wave propagation. In addition to providing insights into the issues that arise when seismic surveying under basalts, these observations suggest that careful management of seismic noise and the acquisition of long-offset seismic data with low-frequency geophones have the potential to improve the seismic results.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.09a
• /
• pp.83-92
• /
• 2005

Fluid conductivity logging has been applied in the boreholes to identify the permeable fi:actures and estimate the origin of saline groundwater in coast area. Fluid replacement technique measures the fluid electrical conductivity with depth at different times in a well after the borehole is first washed out with different water by passing a tube to the borehole bottom. Then formation water flows into the borehole through aquifer such as permeable fractures or porous formation during ambient or pumping condition. Measured conductivity profiles with times therefore indicate the locations of permeable zone or fractures within the open hole or the fully slotted casing hole. As a result of fluid conductivity logging for three boreholes in the study area, it is interpreted that saline groundwater is caused by seawater intrusion through fractured rock, although the effect by land reclamation partially remains. We are planning the quantitative analysis to estimate the hydraulic characteristics using fluid replacement technique, and this approach might be usefully utilized for assessing the characteristics of seawater intrusion, the design of optimal pumping, and estimating the hydraulic properties in coastal aquifer.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.147-154
• /
• 2009

Physical properties of rocks, such as velocity, are strongly dependant on detailed pore structures, and recently, pore micro-structures by X-ray tomography techniques have been used to simulate and understand the physical properties. However, the smoothing effect during the tomographic reconstruction procedure often causes an artifact - overestimating the contact areas between grains. The pore nodes near a grain contact are affected by neighboring grain nodes, and are classified into grain nodes. By this artifact, the pore structure has higher contact areas between grains and thus higher velocity estimation than the true one. To reduce this artifact, we tried two image processing techniques - sharpening filter and neural network classification. Both methods gave noticeable improvement on contact areas between grains visually; however, the estimated velocities showed only incremental improvement. We then tried to change the resolutions of tomogram and quantify its impact on velocity estimation. The estimated velocity from the tomogram with higher spatial resolution was improved significantly, and with around 2 micron spatial resolution, the calculated velocity was very close to the lab measurement. In conclusion, the resolution of pore micro-structure is the most important parameter for accurate estimation of velocity using pore-scale simulation techniques. Also the estimation can be incrementally improved if combined with image processing techniques during the pore-grain classification.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.207-212
• /
• 2004

This article reviews recent developments in three-dimensional (3-D) magntotelluric (MT) imaging. The inversion of MT data is fundamentally ill-posed, and therefore the resultant solution is non-unique. A regularizing scheme must be involved to reduce the non-uniqueness while retaining certain a priori information in the solution. The standard approach to nonlinear inversion in geophysis has been the Gauss-Newton method, which solves a sequence of linearized inverse problems. When running to convergence, the algorithm minimizes an objective function over the space of models and in the sense produces an optimal solution of the inverse problem. The general usefulness of iterative, linearized inversion algorithms, however is greatly limited in 3-D MT applications by the requirement of computing the Jacobian(partial derivative, sensitivity) matrix of the forward problem. The difficulty may be relaxed using conjugate gradients(CG) methods. A linear CG technique is used to solve each step of Gauss-Newton iterations incompletely, while the method of nonlinear CG is applied directly to the minimization of the objective function. These CG techniques replace computation of jacobian matrix and solution of a large linear system with computations equivalent to only three forward problems per inversion iteration. Consequently, the algorithms are efficient in computational speed and memory requirement, making 3-D inversion feasible.

• The Journal of Engineering Geology
• /
• v.25 no.3
• /
• pp.349-356
• /
• 2015

Although geophysical methods are useful and generally provide valuable information about the subsurface, it is important to recognize their limitations. A common limitation is the lack of sufficient contrast in physical properties between different layers. Thus, multiple methods are commonly used to best constrain the physical properties of different layers and interpret each section individually. Ground penetrating radar (GPR) and shallow seismic reflection (SSR) methods, used for shallow and very shallow subsurface imaging, respond to dielectric and velocity contrasts between layers, respectively. In this study, we merged GPR and SSR data from a test site within the Cheongui granitic mass, where the water table is ~3 m deep all year. We interpreted the data in combination with field observations and existing data from drill cores and well logs. GPR and SSR reflections from the tops of the sand layer, water table, and weathered and soft rocks are successfully mapped in a single section, and they correlate well with electrical resistivity data and SPS (suspension PS) well-logging profiles. In addition, subsurface interfaces in the integrated section correlate well with S-wave velocity structures from multi-channel analysis shear wave (MASW) data, a method that was recently developed to enhance lateral resolution on the basis of CMP (common midpoint) cross-correlation (CMPCC) analysis.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.127-136
• /
• 2008

To identify the potential area of gas hydrate in the Ulleung Basin, 2-D and 3-D seismic surveys using R/V Tamhae II were conducted in 2005 and 2006. Seismic survey equipment consisted of navigation system, recording system, streamer cable and air-gun source. For reliable velocity analysis in a deep sea area where water depths are mostly greater than 1,000 m and the target depth is up to about 500 msec interval below the seafloor, 3-km-long streamer and 1,035 $in^3$ tuned air-gun array were used. During the survey, a suite of quality control operations including source signature analysis, 2-D brute stack, RMS noise analysis and FK analysis were performed. The source signature was calculated to verify its conformity to quality specification and the gun dropout test was carried out to examine signature changes due to a single air gun's failure. From the online quality analysis, we could conclude that the overall data quality was very good even though some seismic data were affected by swell noise, parity error, spike noise and current rip noise. Especially, by checking the result of data quality enhancement using FK filtering and missing trace restoration technique for the 3-D seismic data inevitably contaminated with current rip noises, the acquired data were accepted and the field survey could be conducted continuously. Even in survey areas where the acquired data would be unsuitable for quality specification, the marine seismic survey efficiency could be improved by showing the possibility of noise suppression through onboard data processing.

• The Bulletin of The Korean Astronomical Society
• /
• v.43 no.1
• /
• pp.72.4-73
• /
• 2018

북극 스발바르의 사면 지형에는 걸리가 발달되어 있다. 이러한 걸리는, 그 성인에는 여러 의견이 있으나, 화성에도 중고위도를 중심으로 다수 분포한다. 화성의 걸리는 2000년대에 들어 비로소 본격적으로 규명되고 있으나, 지형적 특성으로 인한 탐사의 한계로 지구에 분포하는 유사지형을 통한 비교 연구가 일반적이다(Costard, et al. 2007 등). 이 연구에서는 스발바르의 주도 롱이어비엔에서 UAV을 이용하여 획득한 DEM으로 스발바르 걸리를 측량하고, 이를 화성 중위도의 테라 사이메리아, 테라 시레넘, 노아 키스 테라에 분포하는 걸리와 비교하였다. Longyearbreen 빙하 전방에 위치한 사면을 UAV로 촬영하고, 이를 SfM-MVS(Structure from Motion & MultiView Stereo) 기법으로 3차원 점군 모델과 고해상도 DEM을 제작하여 분석하였다. 화성의 경우 MRO궤도 탐사선이 촬영한 HiRise DTM을 이용하여 분석하였다. 두 걸리는 기후와 지질 조건에 차이가 있음에도 불구하고 유사한 패턴을 보였다. 특히 테라 사이메리아에 위치한 걸리와 롱이어비엔 북사면의 걸리는 기준거리, 단면적, 폭, 경사, 제방 두께 등에서 상당한 정량적 유사관계가 있었다. 이는 두 행성의 걸리가 유사한 성인 및 형성 프로세스를 거쳤을 가능성을 시사한다. 측량 기법과 UAV 의 안정성을 개선시키면 지형 모델의 품질 향상 및 극지에서의 UAV 운용이 용이해질 것으로 기대된다. 또한 스발바르의 기후 요소 및 물리량 적용은 향후 화성 지형연구에도 응용할 수 있을 것으로 사료된다.

• Journal of the Korean Geophysical Society
• /
• v.7 no.4
• /
• pp.313-324
• /
• 2004

The finite element method and statistics of the convergence measurement are useful method of the stability analysis of the tunnel adjacent to the pier. It is the purpose of the this case study to certificate of validity of the application of those methods. The safety of the pilot tunnel method and LW pre-grouting has been evaluated from the FEM analysis. The three-dimensional finite element method is carried out for the decision of the level of stress redistribution at the two-dimensional numerical analysis. An analysis of the convergence is carried out by the estimation of preceding convergence at tunnel excavation. F-examination is applied for this estimation. As results of that analysis, The F-value is from 10.81 to 158.74 and the coefficient of determination is from 0.82 to 0.99. An analysis of convergence is carried out by using regression analysis. Consequently, it is shown that the convergence can be modeled as following function C(t) = a[1-exp(-bt)].