• 제목/요약/키워드: computational geophysics

검색결과 36건 처리시간 0.029초

CCQC modal combination rule using load-dependent Ritz vectors

  • Xiangxiu Li;Huating Chen
    • Structural Engineering and Mechanics
    • /
    • 제87권1호
    • /
    • pp.57-68
    • /
    • 2023
  • Response spectrum method is still an effective approach for the design of buildings with supplemental dampers. In practice, complex complete quadratic combination (CCQC) rule is always used in the response spectrum method to consider the effect of non-classical damping. The conventional CCQC rule is based on exact complex mode vectors. Sometimes the calculated complex mode vectors may be not excited by the external loading and errors in the structural responses always arise due to the mode truncation. Load-dependent Ritz (LDR) vectors are associated with the external loading and LDR vectors not excited can be automatically excluded. Also, contributions of higher modes are implicitly contained in the LDR vectors in terms of static responses. To improve the calculation efficiency and accuracy, LDR vectors are introduced in the CCQC rule in the present study. Firstly, the generation procedure of LDR vectors suitable for non-classical damping system is presented. Compared to the conventional LDR vectors, the LDR vectors herein are complex-valued and named as complex LDR (CLDR) vectors. Based on the CLDR vectors, the CCQC rule is then rederived and an improved response spectrum method is developed. Finally, the effectiveness of the proposed method in this paper is verified through three typical non-classical damping buildings. Numerical results show that the CLDR vector is superior to the complex mode with the same number in the calculation. Since the generation of CLDR vectors requires less computational cost and storage space, the method proposed in this paper offers an attractive alternative, especially for structures with a large number of degrees of freedom.

모델 제약조건이 적용된 MWNI (Minimum Weighted Norm Interpolation)를 이용한 트레이스 내삽 (Trace Interpolation using Model-constrained Minimum Weighted Norm Interpolation)

  • 최지현;송영석;최지훈;변중무;설순지;김기영;이정모
    • 지구물리와물리탐사
    • /
    • 제20권2호
    • /
    • pp.78-87
    • /
    • 2017
  • MWNI (Minimum Weighted Norm Interpolation)를 이용한 내삽 방법은 고차원으로 확장이 용이하고 상대적으로 계산 속도가 빠르다는 장점을 가지고 있으나 알리아스 효과가 존재하는 자료의 내삽에 취약하다. 이런 문제의 개선을 위해 제안된 방법이 모델제약(model-constrained) MWNI이다. 이 논문에서는 MWNI를 이용한 방법과 모델제약 MWNI 방법의 두가지 모듈을 개발한 후 알리아스 효과가 존재하는 자료의 내삽 결과를 비교하였다. 시공간 영역(t-x domain)과 주파수-파수 영역(f-k domain)의 결과 그림을 통해서 모델제약 MWNI를 적용했을 때의 결과가 더 효과적임을 확인할 수 있었다. 동해 울릉분지의 가스 하이드레이트 부존 지역의 현장 자료에 내삽을 적용한 결과, 가스침니 구간 전후로 진폭이 급격하게 변하는 자료에서도 내삽이 가능함을 확인할 수 있었다. 또한 매우 불규칙하고 넓은 구간에서 누락된 인공지진파 자료의 정규화를 통해 신호의 연결성 향상이 가능함을 보일 수 있었다. 결과적으로 이 논문에서 개발된 모듈은 현장의 다양한 여건에 의해 불규칙하거나 넓은 간격으로 얻어진 탄성파 자료의 정규화나 내삽에 효율적으로 활용될 수 있음을 확인하였다.

전자탐사 자료 해석을 위한 1차장 계산 (Primary Solution Evaluations for Interpreting Electromagnetic Data)

  • 김희준;최지향;한누리;송윤호;이기하
    • 지구물리와물리탐사
    • /
    • 제12권4호
    • /
    • pp.361-366
    • /
    • 2009
  • 전자탐사에서 수평 층서구조에 대한 Green 함수의 계산은 전자기 반응의 모델링에서 핵심적인 부분을 담당한다. 해석적으로 구해진 핵함수의 Hankel 변환으로 계산되는 Green 함수는 핵함수의 대수적 등가 표현방식에 의해 그 정확도가 결정된다. 특히 3차원 모델링의 경우 Green 함수 계산 횟수가 매우 많아서 Hankel 변환 계산이 전체 계산시간의 상당량을 차지하므로, 빠르고 정확한 Hankel 변환의 계산을 위해서 선형 수치필터를 이용한다. 최근 많이 시도되는 3차원 역산을 위한 모델링에서는 송신점에서의 특이성 문제를 피하기 위해 전기장을 1차장과 2차장으로 나누어 계산하는 것이 보통이다. 이 연구에서는 균질 반무한공간에 대해 지표면에 놓인 다섯 종류의 송신원에 대한 지하 매질에서의 전기장 세성분을 Hankel 변환을 이용하여 정리하고, 그 계산 방법에 대해 고찰하였다. 그리고 2중 반공간에서 EM1D를 이용하여 공기와 바다의 영향을 모두 고려한 전자기장을 계산할 때, 보다 정확한 해의 계산을 위해 TE 및 TM 모드에서의 반사계수를 유도하였다. 여기서 정리한 해를 이용하면 MT 문제는 물론, 해양 전자탐사의 경우에도 1차장을 정확히 계산할 수 있으므로 3차원 역산에서 보다 정확하고 효율적인 감도 계산이 가능할 것이다.

심층 신경망을 이용한 탄성파 속도 모델 구축 사례 분석 (Case Analysis of Seismic Velocity Model Building using Deep Neural Networks)

  • 조준현;하완수
    • 지구물리와물리탐사
    • /
    • 제24권2호
    • /
    • pp.53-66
    • /
    • 2021
  • 속도 모델 구축은 탄성파 탐사 자료처리에서 필수적인 절차이다. 주시 토모그래피나 속도 분석과 같은 기존 기법들은 하나의 속도 모델을 예측하는 데 계산 시간이 오래 걸리며 역산 결과의 품질이 전문가의 판단에 크게 의존한다. 전파형 역산 또한 초기 속도 모델에 크게 의존한다는 문제가 있다. 최근 심층 신경망 기법이 복잡하고 비선형적인 문제를 푸는데 적용되는 사례가 많아지면서 널리 보급되고 있다. 이 논문에서는 심층 신경망 기법을 이용한 탄성파 속도 모델 구축 사례들을 각 연구에 사용한 신경망에 따라 분류하며 조사하였다. 또한 훈련용 인공 속도 모델 생성 사례도 포함하였다. 심층 신경망은 대량의 데이터로부터 신경망을 훈련함으로써 모델 매개변수를 자동으로 최적화한다. 따라서 기존 기법들에 비해 역산 결과에 사람의 판단이 개입될 여지가 적으며 훈련을 마친 후 하나의 속도 모델을 예측하는 비용은 무시할 수 있다. 또한, 심층 신경망은 전파형 역산과 달리 초기 속도 모델이 필요하지 않다. 여러 연구에서 계산 비용뿐만 아니라 역산 결과에서도 심층 신경망 기법이 뛰어난 성과를 달성하는 것을 보여주었다. 연구 결과들을 바탕으로 속도 모델 구축에 사용된 심층 신경망 기법의 특징에 대해 분석하고 논의하였다.

다변량 통계 기법을 이용한 물리검층 자료로부터의 암석물리학상 결정 (Automatic Electrofacies Classification from Well Logs Using Multivariate Statistical Techniques)

  • 임종세;김정환;강주명
    • 지구물리와물리탐사
    • /
    • 제1권3호
    • /
    • pp.170-175
    • /
    • 1998
  • 이 연구는 다변량 통계 기법을 이용한 물리검층 자료로부터의 암석물리학상 결정으로 암상을 예측하는 것이다. 기술 통계 분석으로 물리검층 자료의 특성을 파악하고 주성분 분석에 의한 다변량 검층 자료들의 상관도 분석을 통해 변수들을 변환시켜 새로운 변수인 주성분을 구하고 변수들의 차원을 축소한다. 통계적 방법에 의한 주성분 검층 자료의 구획에 의한 효율적 자료 축소와 계산의 효율성을 높여 양질의 해석결과를 얻을 수 있다. 구획된 주성분 검층 자료로부터 계보적 군집 분석에 의해 암석물리학상을 결정한다. 최적 암석물리학상의 수는 전체 변동과 군집내의 변동사이의 비와 코어자료 등에 의해 비교 결정된다. 이 연구에서 개발된 암석물리학상 결정법을 국내대륙붕 물리검층자료에 적용한 결과 결정된 암석물리학상은 시추 코어 및 시추 암편 분석에 의한 암상 구분화와 잘 일치하였다. 이러한 연구는 저류층 특성인자의 신뢰성 있고 정량적인 평가로 유전 개발 및 생산 계획 시 유용한 도구로 활용될 수 있을 것이다.

  • PDF

3차원 MT 역산에서 CG 법의 효율적 적용 (Conjugate Gradient Least-Squares Algorithm for Three-Dimensional Magnetotelluric Inversion)

  • 김희준;한누리;최지향;남명진;송윤호;서정희
    • 지구물리와물리탐사
    • /
    • 제10권2호
    • /
    • pp.147-153
    • /
    • 2007
  • CG (conjugate gradient) 법은 선형 연립방정식을 반복적으로 푸는 가장 효율적인 해법 중 하나이고, 또한 비선형 최소자승문제에도 적용할 수 있다. 자기지전류(MT) 역산 문제를 풀 때에는 최소자승문제의 목적함수 자체의 최소화에 직접 CG 법을 적용하거나, Gauss-Newton 법에 기초한 반복역산의 각 반복단계에서 모형의 변화량 계산에 CG 법을 이용할 수 있다. CG 법을 적용할 경우, 임의의 벡터에 대한 감도행렬의 영향 및 그 전치행렬의 전치행렬의 영향을 감도행렬을 직접 구하지 않고 계산할 수 있다는 장점이 있기 때문에 감도행렬의 계산 규모가 방대한 3차원 역산 문제에서 계산시간을 월등히 줄일 수 있다.

강인한 역산으로서의 하이브리드 $l^1/l^2$ norm IRLS 방법의 효율적 구현기법 (An Efficient Implementation of Hybrid $l^1/l^2$ Norm IRLS Method as a Robust Inversion)

  • 지준
    • 지구물리와물리탐사
    • /
    • 제10권2호
    • /
    • pp.124-130
    • /
    • 2007
  • 탄성파 역산에 있어서 가장 널리 사용되는 최소자승($l^2$ norm)해는 이상치(outlier)에 매우 민감하게 반응하는 경향이 있다. 이에 반해서 $l^1$ norm을 최소화하는 해는 이상치에 강인한 면을 보이나 일반적으로 좀 더 많은 계산이 필요하다. 반복적가중의 최소자승법(Iteratively reweighted least squares [IRLS] method)을 이용하면 이러한 $l^1$ norm 문제의 근사해(approximate solution)를 효율적으로 구할 수 있다. 본 논문에서는 작은 크기의 잔여분은 $l^2$ norm으로 처리하며, 큰 크기의 잔여분은 $l^1$ norm으로 처리하는 하이브리드 $l^1/l^2$ norm 최소화를 IRLS 방법에 쉽게 적용하는 구현 기법을 소개한다. 소개된 알고리즘은 특이치(singularity)처리를 위한 임계값의 결정에 민감하게 반응하는 기존의 $l^1$ norm IRLS 방법과는 달리 임계값 결정에 상관없이 늘 강인한 역산의 특성을 보여준다.

GRRI를 이용한 2차원 MT 탐사자료의 역산 (Generalized Rapid Relaxation Inversion of Two-Dimensional Magnetotelluric Survey Data)

  • 정용현;서정희;신창수
    • 지구물리와물리탐사
    • /
    • 제1권1호
    • /
    • pp.71-78
    • /
    • 1998
  • 2차원 MT탐사자료의 역산은 보편적으로 많은 계산시간과 기억용량이 소요되며 TM모드의 경우지표 근처의 불균질대에 야기될 수 있는 정적 효과에 대한 주의 깊은 고려도 필요하게 된다. 이에 본 연구에서는 GRRI를 이용하여 TM모드 MT자료의 2차원 역산 알고리듬을 구현하였다. 본 역산 알고리듬은 국부적인 2차원 해석을 통한 모형변수 증분 해석에 기초를 두는데 이는 모형변수의 증분을 구하는 과정에 있어서 탐사의 대상인 2차원 전 영역을 동시에 고려하지 않고 국부 영역들로 분할하여 순차적으로 고려함으로써 수평적인 변화가 고려되면서도 자코비안을 구하는 과정에서 사용되는 부가적인 모델링 과정을 배제하여 계산시간과 기억용량을 최소화하기 위한 것이다. 구현된 알고리듬의 타당성 검증은 간단한 이론모형에 대한 적용을 통하여 이루어졌으며, 이를 토대로 타 알고리듬과의 비교를 통하여 그 특징을 살펴보는 한편, 여러 가지 다양한 모형에 적용하여 그 현장 적용성을 고찰하였다. 본 알고리듬을 이용한 역산은 계산시간과 기억용량 면에서 매우 경제적이며 전류집중에 의해 탐사자료에 왜곡이 생기는 경우에는 특별한 전처리과정 없이도 타당한 분해능을 확인할 수 있었다. GRRI를 이용한 역산 알고리듬은 향후 현장에서 실시간 역산으로서 경제적이면서 타당한 자료 해석에 적용성이 기대된다.

  • PDF

제온 파이 x200 프로세서를 이용한 3차원 음향 파동 전파 모델링 병렬 연산 성능 비교 (Comparison of Parallel Computation Performances for 3D Wave Propagation Modeling using a Xeon Phi x200 Processor)

  • 이종우;하완수
    • 지구물리와물리탐사
    • /
    • 제21권4호
    • /
    • pp.213-219
    • /
    • 2018
  • 본 연구에서는 제온 파이 x200 프로세서를 이용하여 3차원 파동 전파 모델링을 수행하고 기존의 제온 CPU를 사용한 경우와 병렬 연산 성능을 비교하였다. 제온 파이 1세대 프로세서인 제온 파이 나이츠 코너 보조프로세서와 달리 제온 파이 2세대 프로세서인 x200 프로세서는 직접 운영체제 실행이 가능하므로 내장 메모리와 주메모리 사이의 추가적인 통신이 필요 없다. 또한 제온 파이 x200 프로세서는 대용량 주메모리와 고대역폭 메모리를 이용하여 대규모 컴퓨팅을 독립적으로 실행할 수 있다. 병렬 연산 성능 비교를 위해 MPI (Message Passing Interface)와 OpenMP (Open Multi-Processing)를 이용해 모델링을 수행하였다. SEG/EAGE 암염돔 모델을 이용한 수치 실험 결과 제온 파이에서 다량의 연산 코어와 고대역폭 메모리를 이용해 12 코어 CPU 대비 2.69 ~ 3.24배 우수한 모델링 성능을 얻을 수 있었다.

BERT 기반 사전학습을 이용한 탄성파 자료처리: 송신원 모음 배열 비교 (Seismic Data Processing Using BERT-Based Pretraining: Comparison of Shotgather Arrays)

  • 신영재
    • 지구물리와물리탐사
    • /
    • 제27권3호
    • /
    • pp.171-180
    • /
    • 2024
  • 탄성파 자료처리는 탄성파 자료를 분석하여 지구 내부 구조와 특성을 파악하는 기술로, 높은 컴퓨터 연산력이 요구된다. 이러한 도전 과제를 해결하기 위해 머신러닝 기술이 도입되었으며, 잡음 제거, 속도 모델 구축 등 다양한 작업에서 활용되고 있다. 그러나, 대부분의 연구는 특정 탄성파 처리 작업에 집중되어 있어 자료에 내재된 유사한 특징과 구조를 충분히 활용하지 못하는 한계가 있다. 본 연구에서는 BERT (Bidirectional Encoder Representations from Transformers) 기반의 사전학습을 위해 단일 송신원 모음에서 수신기별 시계열 자료('수신기 배열')와 동일 시간에 기록된 수신기 신호('시간 배열')를 입력 자료로 활용하는 방법을 비교하였다. 이를 위해 단층을 포함한 속도 모델에서 생성한 합성 송신원 모음 자료를 이용하여 잡음 제거, 속도 추정, 그리고 단층 확인 작업을 수행하였다. 임의 잡음 제거 작업에서는 수신기 및 시간 배열에서 모두 좋은 성능을 보였으나, 공간적인 분포 파악이 요구되는 속도 추정 및 단층 확인 작업에서는 시간 배열의 결과가 상대적으로 더 우수함을 확인하였다.