• Title/Summary/Keyword: Numerical inversion

Search Result 293, Processing Time 0.026 seconds

Three-dimensional Inversion of Resistivity Data (전기비저항 탐사자료의 3차원 역산)

  • Yi Myeong-Jong;Kim Jung-Ho;Cho Seong-Jun;Chung Seung-Hwan;Song Yoonho
    • Geophysics and Geophysical Exploration
    • /
    • v.2 no.4
    • /
    • pp.191-201
    • /
    • 1999
  • The interpretation of resistivity data has, so far, mainly been made under the assumption that the earth is of relatively simple structure and then using one or two-dimensional inversion scheme. Since real earth structure and topography are fully three-dimensional and very complicated In nature, however, such assumptions often lead to misinterpretation of the earth structures. In such situations, three-dimensional inversion is probably the only way to get correct image of the earth. In this study, we have developed a three-dimensional inversion code using the finite element solution for the forward problem. The forward modeling algorithm simulates the real field situation with irregular topography. The inverse problem is solved iteratively using the least-squares method with smoothness constraint. Our inversion scheme employs ACB (Active Constraint Balancing) to enhance the resolving power of the inversion. Including Irregular surface topography in the inversion, we can accurately define the earth structures without artifact in the numerical tests. We could get reasonable image of earth structure by Inverting the real field data sets taken over highway bridge construction site.

  • PDF

Analysis of cooling phenomenon of water with the supercooled (과냉각을 동반한 순수물의 냉각현상 해석)

  • Chu, Mi-Seon;Yun, Jeong-In;Kim, Jae-Dol;Kamata, Yoshinobu;Kato, Toyofumi
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.21 no.7
    • /
    • pp.862-872
    • /
    • 1997
  • Ice formation in a horizontal circular cylinder has been studied numerically. From the numerical analysis results, it was found that there were three types of freezing pattern and that freezing phenomenon was affected largely by density inversion and cooling rate. The type of freezing pattern largely depends on the secondary flow which is generated by density inversion. When supercooling energy is released before the development of the secondary flow, the annular ice layer grows. If the energy is released when the secondary flow is considerably developed and the supercooled region is removed to the upper half part of the cylinder, an asymmetric ice layer grows. And if the energy is released after perfect development of the secondary flow, instantaneous dendritic ice formation over the full region occurs. Furthermore, this secondary flow was found to have an effect on heat transfer characteristics. The heat transfer rate becomes small at the instant when the secondary flow is generated, but becomes large with the development of the flow. It's concluded that for the facilitation of heat transfer it is desirable to keep water in liquid phase until the secondary flow is perfectly developed. This study gave an instruction of performance improvement of capsule type ice storage tank.

Verification of Wavefront Inversion Scheme via Signal Subspace Comparison Between Physical and Synthesized Array Data in SAT Imaging (SAR Imaging에서 Physical Array와 합성 Array 신호의 Subspace 비교를 통한 Wavefront Inversion 기법 입증)

  • 최정희
    • Journal of the Korean Institute of Telematics and Electronics D
    • /
    • v.36D no.4
    • /
    • pp.34-41
    • /
    • 1999
  • Unlike the traditional radar system, Synthetic Aperture Radar(SAR) system is capable of imaging a target scene to ceertain degree of cross-range resolution. And this resolution is mainly depends on the size of aperture synthesized. Thus, a good system model and inversion scheme should be developed to actually give effect of synthesizing aperture size, which in turn gives better cross range resolution of reconstructed target scene. Among several inversion schemes for SAR imaging, we used an inversion scheme called wavefront reconstruction which has no approximation in wave propagation analysis, and tried to verify whether the collected data with synthesized aperture actually give the same support as that with physical aperture in the same size. To do this, we performed a signal subspace comparison of two imaging models with physical and synthesized arrays, respectively. Theoretical comparisons and numerical analysis using Gram-Schmidt procedures have been performed. The results showed that the synthesized array data fully span the physical array data with the same system geometry. This result strongly supports the previously proposed inversion scheme valuable in high resolution radar imaging.

  • PDF

Inversion of Acoustical Properties of Sedimentary Layers from Chirp Sonar Signals (Chirp 신호를 이용한 해저퇴적층의 음향학적 특성 역산)

  • 박철수;성우제
    • The Journal of the Acoustical Society of Korea
    • /
    • v.18 no.8
    • /
    • pp.32-41
    • /
    • 1999
  • In this paper, an inversion method using chirp signals and two near field receivers is proposed. Inversion problems can be formulated into the probabilistic models composed of signals, a forward model and noise. Forward model to simulate chirp signals is chosen to be the source-wavelet-convolution planewave modeling method. The solution of the inversion problem is defined by a posteriori pdf. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the ranges for a priori uniform distribution is based. The genetic algorithm can be applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L₂norm of the difference between measured and modeled signals. The observed signals can be separated into a set of two signals reflected from the upper and lower boundaries of a sediment. The separation of signals and successive applications of the genetic algorithm optimization process reduce the search space, therefore improving the inversion results. Not only the marginal pdf but also the statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm. The examples applied here show that, for synthetic data with noise, it is possible to carry out an inversion for sedimentary layers using the proposed inversion method.

  • PDF

An Inversion Package for Interpretation of Microgravity Data (고정밀 중력탐사 자료 역산 패키지)

  • Park, Yeong-Sue;Rim, Hyoungrea;Lim, Mutaek;Chung, Hojoon
    • Geophysics and Geophysical Exploration
    • /
    • v.20 no.4
    • /
    • pp.226-231
    • /
    • 2017
  • Since microgravity survey aims to delineate subsurface density structures in small scale, it requires inversion method, which is able to resolve small scale structures. It can be achieved by adopting a stabilizing functional which separates density boundary distinctly, which is different concept from general inversion routines. We composed Matlab-based interactive two-dimensional microgravity data inversion package containing several kinds of inversion routines with different stabilizing functional, for handling various geologic conditions and survey purposes. Different kinds of inversion routines in the package were verified and examined with representative synthetic data sets generated by numerical modeling. In addition, we applied the developed package to a real microgravity survey data.

Three-dimensional Resistivity Inversion Including Topographic Effect (지형효과를 포함한 3차원 전기비저항 역산)

  • 박종오;김희준;송무영
    • The Journal of Engineering Geology
    • /
    • v.14 no.1
    • /
    • pp.21-28
    • /
    • 2004
  • Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

Development of Automated Inversion Method for HWAW Method Using Genetic Algorithm (유전자 알고리즘을 이용한 HWAW 방법을 위한 자동화 역산 방법의 개발)

  • Park, Hyung-Choon;Hwang, Hea-Jin
    • Journal of the Korean Geotechnical Society
    • /
    • v.28 no.8
    • /
    • pp.55-63
    • /
    • 2012
  • The evaluation of shear modulus (or shear wave velocity) profile of the site is very important in various fields of geotechnical engineering and various surface wave methods have applied to determine the shear wave velocity profiles and showed good performance. Surface wave methods evaluate the dispersion curve in the field and determine the shear wave velocity profile through the inversion process. In this paper, the automated inversion process using the genetic algorithm is developed for HWAW method which is one of surface wave methods recently developed. The proposed method uses the error function based on the wavelength domain dispersion curve and can determine the reliable shear wave velocity profile not only in shallow depth but also in deep depth. To estimate the validity of the proposed method, numerical simulations and field test were performed and the proposed method was applied to determine the shear wave velocity profiles. Through the numerical simulations and field applications, the promising potential of the proposed method was verified.

Anewwaveform inversion methodto determine the rupture directivity of moderate earthquakes: numerical tests for rupture models (중간 규모 지진의 단층 파해 방향성 결정을 위한 새로운 주파수 영역 역산방법: 파쇄 전파 모델을 이용한 수치 시험)

  • Yoo, Seung-Hoon;Rhie, Jun-Kee
    • Geophysics and Geophysical Exploration
    • /
    • v.12 no.1
    • /
    • pp.114-120
    • /
    • 2009
  • Rupture directivity is the important parameter in estimating damage due to earthquakes. However, the traditional moment tensor inversion technique cannot resolve the real fault plane or the rupture directivity. To overcome these limitations, we have developed a new inversion algorithm to determine the moment tensor solution and the rupture directivity for moderate earthquakes, using the waveform inversion technique in the frequency domain. Numerical experiments for unilateral and bilateral rupture models with various rupture velocities confirm that the method can resolve the ambiguity of the fault planes and the rupture directivity successfully. To verify the feasibility of the technique, we tested the sensitivity to velocity models, which must be the most critical factor in practice. The results of the sensitivity tests show that the method can be applied even though the velocity model is not perfect. If this method is applied in regions where the velocity model is well verified, we can estimate the rupture directivity of a moderate earthquake. This method makes a significant contribution to understanding the characteristics of earthquakes in those regions.

Effects of Pressure Ratio on Population Inversion in a DF Chemical Laser with Concurrent Lasing

  • Park, Jun-Sung;Baek, Seung-Wook
    • Proceedings of the Korean Society of Propulsion Engineers Conference
    • /
    • 2004.03a
    • /
    • pp.287-293
    • /
    • 2004
  • A numerical simulation is presented for investigating the effects of pressure ratio of $D_2$ injector to supersonic nozzle on the population inversion in the DF chemical laser cavity, while a lasing concurrently takes place. The laser beam is generated between the mirrors in the cavity and it is important to obtain stronger population inversion and more uniform distribution of the excited molecules in the laser cavity in order to produce high power laser beam with good quality. In this study, these phenomena are investigated by means of analyzing the distributions of the DF excited molecules and the F atom used as an oxidant, while simultaneously estimating the maximum small signal and saturated gains and power in the DF chemical laser cavity. For the numerical solution, an 11-species (including DF molecules in various excited states of energies), 32-step chemistry model is adopted for the chemical reaction of the DF chemical laser system. The results are discussed by comparison with two $D_2$injector pressure cases; 192 torr and 388.64 torr. Major results reveal that in the resonator, stronger population inversions occur in the all transitions except DF(1)-DF(0), when the $D_2$injection pressure is lower. But, the higher $D_2$injection pressure provides a favorable condition for DF(1)-DF(0) transition to generate the higher power laser beam. In other words, as the pressure of $D_2$injector increases, the maximum small signal gain in the $V_{1-0}$ transition, which is in charge of generating most of laser power, becomes higher. Therefore, the total laser beam power becomes higher.r.

  • PDF