• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.27-35
• /
• 2013

In this study, we developed the trial product on 3D magnetic inversion tentatively named 'KMag3D'. Also, we briefly introduced its own function and graphic user interface on which especially focused through the development in the form of user manual. KMag3D is consisted of two fundamental frame for the 3D magnetic inversion. First, algebraic reconstruction technique was selected as a 3D inversion algorithm instead of least square method conventionally used in various magnetic inversion. By comparison, it was turned out that algebraic reconstruction algorithm was more effective and economic than that of least squares in aspect of both computation time and memory. Second, for the effective determination of the 3D initial and a-priori information model required in the execution of our algorithm, we proposed the practical technique based on the assemblage of 2D forward modeling and inversion results for individual user-selected 2D profiles. And in succession, initial and a-priori information model were constructed by appropriate interpolation along the strke direction. From this, we concluded that our technique is both suitable and very practical for the application of 3D magentic inversion problem.

• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.101-109
• /
• 1998

Using airborne magnetic data, magnetic characteristics were studied at the Samryangjin caldera area developed in the volcanics of the Yuchon sub-basin, the south eastern part of the Gyeongsang basin. Residual magnetics, reduction to the pole, horizontal derivative, and vertical derivative maps are prepared. Using these maps, the magnetic lithofaces are zoned and the geological structures such as caldera and faults were qualitatively interpreted. In addition, the two quantitative interpretations were performed. Firstly, the forward modelling were done to the 14.5 line km crossing the caldera area to the northeast-southwest direction. Applying the 3-D Euler deconvolution method to the whole study area, the depth extent and the characteristics of the magnetic anomalous bodies were studied. According to the results, the magnetic lithofaces of the area are zoned by 4 units. In general, these are well matched with the geological distributions. But the biotite granites intruded in the northern boundary of the Samryangjin caldera show the high magnetic intensity, while the biotite granites of the other areas show the low magnetic intensity and the different magnetic lithofaces. Thus, we interpreted that the biotite granites are closely related with the volcanic activity of the Samryngjin caldera, and are intruded in the fracture zones developed along the caldera rim. The Samryangjin caldera and fault structures of the area can be easily recognized by the distinct magnetic structures from the various magnetic anomaly maps. Also the topographic characteristics well reflect these structures. The results of the forward modelling show that the magnetic basement depth of the Gyeongsang sedimentary basin is on the average about 6 km and in maximum 10 km. And the depth becomes shallower toward the caldera boundary due to the shallow intrusion of the volcanics. The results of the 3-D Euler method also show the caldera and fault structures. And the relatively shallow magnetic anomalous bodies which are related with the volcanics are generally developed to the east-west and northeast directions, while the deep magnetic anomalous bodies to the northwest direction.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.1
• /
• pp.30.3-30.3
• /
• 2021

In order to generate the 3D structure of the 21-cm signal during the reionization, semi-numerical simulations based on Excursion set formalism are broadly used. However, semi-numerical simulations in the realization of the 3D structure are known to be the ionizing photon non-conserving by the structure of the Excursion set approach. Recently, explicit photon conserving algorithms for semi-numerical simulations introduced, but they are still too slow when forward modelling the 21-cm signal with high-dimensional parameter spaces. Here, we introduce a new method for approximately correcting photon non-conservation, which can be applied on-the-fly. This method is tailored towards the efficient simulation and Bayesian inference with high-dimensional parameter space. Then, we investigate how large an impact that photon non-conservation has on astrophysical parameter inference by performing an MCMC analysis. We find that the ionizing escape parameter is deviated from the fiducial value by 2 sigma when we infer astrophysical parameters without this correction.

• Journal of the Korean Geophysical Society
• /
• v.6 no.3
• /
• pp.155-164
• /
• 2003

The validity and efficiency of the roll-along technique widely used in 2-D electrical resistivity survey are analyzed in case of the dipole-dipole and the Wenner-Schlumberger arrays by numerical modelling. The shallow anomalous resistivity bodies are successfully inverted both in the dipole-dipole and in the Wenner-Schlumberger arrays because the shallow data of pseudosection are not omitted by the roll-along technique. However, the deep anomalous resistivity bodies can not be well resolved due to the skip of observed data which is more significant in the Wenner-Schlumberger array having relatively poor horizontal coverage of obtaining data. Carrying out electrical survey adopting the dipole-dipole array, the skip of data is insignificant because it is unfeasible to expand the electrodes to the maximum electrode separation coefficient($n_max$) owing to low S/N ratio. In case of the Wenner-Schlumberger array, however, because it is generally feasible to expand the electrodes $n_max$ to the owing to high S/N ratio, it is highly possible that skip of data from the roll-along technique causes significant distortion of inversion results. Therefore, adopting the Wenner-Schlumberger array having deeper median depth(Edwards, 1977) than do the dipole-dipole array on condition of the same unit electrode spacing( ($a$) ) and $n_max$, it is recommended to determine $a$ based on not $n_max$but $n_prob$free from the skip of observing data and forward electrodes with keeping overlap interval 3/4 of the survey line length in order to reduce the distortion of resistivity structure and perform resistivity survey efficiently. These results are confirmed by numerical modelling.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.105-113
• /
• 2009

The interpretation of observed waveform characteristics identified in refraction and wide-angle reflection data increases confidence in the crustal structure model obtained. When calculating traveltimes and raypaths, wavefront methods on a regular grid based on graph theory are robust even with complicated structures, but basically compute only first arrivals. In this paper, we develop new algorithms to compute traveltimes and raypaths not only for first arrivals, but also for fast and later reflection arrivals, later refraction arrivals, and converted waves between P and S, using the modified wavefront method based on slowness network nodes mapped on a multi-layer model. Using the new algorithm, we can interpret reflected arrivals, Pg-later arrivals, strong arrivals appearing behind Pn, triplicated Moho reflected arrivals (PmP) to obtain the shape of the Moho, and phases involving conversion between P and S. Using two models of an ocean-continent transition zone and an oceanic ridge or seamount, we show the usefulness of this algorithm, which is confirmed by synthetic seismograms using the 2D Finite Difference Method (2D-FDM). Characteristics of arrivals and raypaths of the two models differ from each other in that using only first-arrival traveltime data for crustal structure analysis involves risk of erroneous interpretation in the ocean-continent transition zone, or the region around a ridge or seamount.