• 지구물리와물리탐사
• /
• 제17권4호
• /
• pp.231-241
• /
• 2014

최근 엔지니어링 목적의 단일채널 탄성파 탐사가 많이 수행되고 있다. 단일채널 탄성파 탐사는 일반적으로 특별한 자료처리 없이 효율적으로 지하 지질구조를 파악할 수 있는 장점이 있지만, 복잡한 지질구조에 대한 정확한 영상화에는 한계를 가진다. 자원개발 목적의 다중채널 탄성파 탐사에서는 최근 파형역산에 대한 연구가 활발히 진행되고 있으며, 복잡한 지하구조에 대해서도 정확한 지하영상화 결과를 제시하고 있다. 이에 본 논문에서는 단일채널 탄성파 탐사 자료를 이용하여 지하 속도모델을 구하기 위한 탄성파 파형역산 알고리즘을 제안하였다. 제안한 알고리즘은 단일채널 탄성파 탐사를 고려하여 지하 매질을 1차원으로 가정하였으며, 벌림에 의한 지연시간을 제거하여 벌림에 의한 효과를 보정하였다. 파형역산은 안정적인 해의 계산이 가능한 가우스-뉴턴법을 이용하였다. 알고리즘은 수정된 Marmousi2 모델에 적용하여 검증하였으며, 부산항에서 얻은 현장자료에 적용해 보았다.

• 전기학회논문지
• /
• 제58권2호
• /
• pp.415-421
• /
• 2009

TOA(time-of-arrival) and TDOA(time-difference-of-arrival) positioning techniques are commonly used in many radio-navigation systems. From the literature, it is known that the position estimate and error covariance matrix of TDOA obtained by GN(Gauss-Newton) method is exactly the same as that of TOA when the error source of the range measurement is only an IID white Gaussian noise. In case of geo-location and indoor positioning, however, multi-path or NLOS(non-line-of-sight) error is frequently appeared in range measurements. Though its occurrence is random, the multipath acts like a bias for a stationary user if it occurs. This paper presents the comparisons of error characteristics between TOA and TDOA positioning in presence of multi-path or NLOS error. It is analytically shown that the position estimate of TDOA is exactly the same as that of TOA even when bias errors are included in range measurements with different magnitudes. By computer simulation, position estimation error and error distribution are analyzed in presence of range bias errors.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2003년도 추계학술대회 논문집 전력기술부문
• /
• pp.181-183
• /
• 2003

본 연구에서는 현재 국내-외 논문을 참조하여 배전계통에 적용된 조류계산 알고리즘의 적용사례를 조사하였다. 배전계통에 적용된 조류계산알고리즘은 Power Summation Method와 Current Summation Method가 있다. 배전계통은 선로의 구조가 거의 대부분 방사상으로 성되어 있기 때문에 Newton Raphson, Gauss Seidel, Fast Decoupled 방법을 그대로 사용하기는 문제가 있다. 그래서 DistFlow, Forward/Backward sweeping 법 등이 주로 사용되어 왔다. 그러나 미래의 계통은 현재와는 또 다른 Topology로 방사상과 루프가 혼합된 형태로 다르게 구성되어야 할 것이다. 이러한 상황에서 미래의 배전계통에 대하여 적용 가능한 알고리즘을 전망하고져 한다.

• 한국지반공학회지:지반
• /
• 제8권2호
• /
• pp.5-20
• /
• 1992

The physical-based and lumped-parameter hydrologic groundwater flow model for predicting the rainfall-triggered rise of groundwater levels in hillside slopes is developed in this paper to assess the risk of landslides. The developed model consists of a vertical infiltration model for unsaturated zone linked to a linear storage reservoir model(LSRM) for saturated zone. The groundwater flow model has uncertain constants like soil depttL slope angle, saturated permeability, and potential evapotranspiration and four free model parameters like a, b, c, and K. The free model parameters could be estimated from known input-output records. The BARD algorithm is uses as the parameter estimation technique which is based on a linearization of the proposed model by Gauss -Newton method and Taylor series expansion. The application to examine the capacity of prediction shows that the developed model has a potential of use in forecast systems of predicting landslides and that the optimal estimate of potential 'a' in infiltration model is the most important in the global optimum analysis because small variation of it results in the large change of the objective function, the sum of squares of deviations of the observed and computed groundwater levels. 본 논문에서는 가파른 산사면에서 산사태의 발생을 예측하기 위한 수문학적 인 지하수 흐름 모델을 개발하였다. 이 모델은 물리적인 개념에 기본하였으며, Lumped-parameter를 이용하였다. 개발된 지하수 흐름 모델은 두 모델을 조합하여 구성되어 있으며, 비포화대 흐름을 위해서는 수정된 abcd 모델을, 포화대 흐름에 대해서는 시간 지체 효과를 고려할 수 있는 선형 저수지 모델을 이용하였다. 지하수 흐름 모델은 토층의 두께, 산사면의 경사각, 포화투수계수, 잠재 증발산 량과 같은 불확실한 상수들과 a, b, c, 그리고 K와 같은 자유모델변수들을 가진다. 자유모델변수들은 유입-유출 자료들로부터 평가할 수 있으며, 이를 위해서 본 논문에서는 Gauss-Newton 방법을 이용한 Bard 알고리즘을 사용하였다. 서울 구로구 시흥동 산사태 발생 지역의 산사면에 대하여 개발된 모델을 적용하여 예제 해석을 수행함으로써, 지하수 흐름 모델이 산사태 발생 예측을 위하여 이용할 수 있음을 입증하였다. 또한, 매개변수분석 연구를 통하여, 변수 a값은 작은 변화에 대하여 목적함수값에 큰 변화를 일으키므로 a의 값에 대한 최적값을 구하는 것이 가장 중요한 요소라는 결론을 얻었다.

• 대한기계학회논문집
• /
• 제18권10호
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2007년도 공동학술대회 논문집
• /
• pp.329-336
• /
• 2007

For scaling of the gradient of misfit function, we develop a new pseudo-Hessian matrix constructed by combining amplitude field and pseudo-Hessian matrix. Since pseudo- Hessian matrix neglects the calculation of the zero-lag auto-correlation of impulse responses in the approximate Hessian matrix, the pseudo-Hessian matrix has a limitation to scale the gradient of misfit function compared to the approximate Hessian matrix. To validate the new pseudo- Hessian matrix, we perform frequency-domain elastic full waveform inversion using this Hessian matrix. By synthetic experiments, we show that the new pseudo-Hessian matrix can give better convergence to the true model than the old one does. Furthermore, since the amplitude fields are intrinsically obtained in forward modeling procedure, we do not have to pay any extra cost to compute the new pseudo-Hessian. We think that the new pseudo-Hessian matrix can be used as an alternative of the approximate Hessian matrix of the Gauss-Newton method.

• 콘크리트학회지
• /
• 제8권2호
• /
• pp.119-127
• /
• 1996

본 연구는 철근콘크리트 부재의 3차원 재료적 비선형해석을 하기 위한 것이다. 콘크리트는 3축 비선형 응력-변형률 거동, 균열, 파쇄 및 변형률완화를 포함하는 3차원 16절점 고체요소를 사용하고, 철근은 변형률경화를 갖는 3차원 3절점 트러스요소를 사용한다. 균열 후 골재의 맞물림을 고려하는 유효전단계수를 평가하기 위해서 균열의 진행여부에 따른 전단유지계수를 도입하였으며, 수치해를 얻기 위해 수정뉴턴방법을 사용하였다. 가우스점에서의 해석결과는 그래픽으로 확인된다. 수치예제로서 Krahl의 철근콘크리트 보와 Hedgren의 철근콘크리트 쉘을 채택하여 해석결과와 비교하였다.

• Nuclear Engineering and Technology
• /
• 제45권4호
• /
• pp.439-458
• /
• 2013

REX-10 is a fully-passive small modular reactor in which the coolant flow is driven by natural circulation, the RCS is pressurized by a steam-gas pressurizer, and the decay heat is removed by the PRHRS. To confirm design decisions and analyze the transient responses of an integral PWR such as REX-10, a thermal-hydraulic system code named TAPINS (Thermal-hydraulic Analysis Program for INtegral reactor System) is developed in this study. Based on a one-dimensional four-equation drift-flux model, TAPINS incorporates mathematical models for the core, the helical-coil steam generator, and the steam-gas pressurizer. The system of difference equations derived from the semi-implicit finite-difference scheme is numerically solved by the Newton Block Gauss Seidel (NBGS) method. TAPINS is characterized by applicability to transients with non-equilibrium effects, better prediction of the transient behavior of a pressurizer containing non-condensable gas, and code assessment by using the experimental data from the autonomous integral effect tests in the RTF (REX-10 Test Facility). Details on the hydrodynamic models as well as a part of validation results that reveal the features of TAPINS are presented in this paper.

• 한국지구과학회지
• /
• 제37권5호
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• 전자공학회논문지
• /
• 제54권8호
• /
• pp.99-106
• /
• 2017

전기 임피던스 단층촬영법은 주입 전류와 측정 전압을 기반으로 관심 도메인 내부의 도전율/저항률 분포를 복원하는 비파괴 영상 복원 기법이다. 본 논문에서는 역문제 계산시간을 줄이고 더불어 공간 해상도도 향상시키기 위해, 적응 문턱치 기반의 ROI(region of interest) 방법을 제안하였다. INTERMODES 방법에 의해 적응 문턱치가 계산이 되고 이 값을 기반으로 전체 도메인으로부터 ROI가 결정된다. 그리고 영상 복원의 계산 도메인을 ROI 내로 국한시켜 반복적 가우스-뉴턴 방법을 적용하여 저항률 분포를 추정하였다. 제안한 방법의 성능을 평가하기 위해 수치실험을 수행하고 그 결과를 비교분석하였다.