Abstract
We develop the weighted-averaging finite-element method which uses four kinds of element sets. By constructing global stiffness and mass matrices for four kinds of element sets and then averaging them with weighting coefficients, we obtain a new global stiffness and mass matrix. With the optimal weighting coefficients minimizing grid dispersion and grid anisotropy, we can reduce the number of grid points required per wavelength to 4 for a $1\%$ upper limit of error. We confirm the accuracy of our weighted-averaging finite-element method through accuracy analyses for a homogeneous and a horizontal-layer model. By synthetic data example, we reconfirm that our method is more efficient for simulating a geological model than previous finite-element methods.
일반 반사법 탐사 모델링에서 효율적인 주파수영역 수치모델링의 실용화를 위해 무엇보다 해결해야할 과제는 파장당 격자수를 줄이는 것이다. 본 연구에서는 이에 착안하여 수치분산 및 수치이방성을 최소화시키면서 한 파장당 필요한 격자수를 줄일 수 있는 가중평균 유한요소법을 개발하였다. 강성행렬과 질량행렬은 네 개의 사각형 유한요소로 확장하였으며, 모든 격자점이 포함되도록 유한요소를 배열하여 조합하였다. 확장된 네 개의 강성행렬과 질량행렬은 가중평균계수를 주어 선형결합하는 방법으로 가중평균하였다. 가중평균계수는 확장된 25점 평균차분법을 사용하여 가중평균계수를 결정하였다. 또한, 정확도 향상은 2차원 균질모델 과 수평층 모델에서 해석해와 한 파장당 4개의 격자점을 준 가중평균 유한요소법 수치해 비교를 통하여 검증하였다. 또한, 석유탐사에서 활용성이 높은 향사구조 모델을 선정하여 이의 반응을 관찰한 결과 지층경계면외에 네 개의 사각형 유한요소들의 구성으로 인한 인위적인 파의 도달이 인식되지 않았다. 따라서, 본 연구에서 고안된 가중평균 유한요소법은 주파수영역에서 폭 넓은 수치모델링연구을 가능하게 할 것이다.
