미소진동 자료를 이용한 EGS 사이트에서의 효율적인 모니터링 알고리듬 개발

Abstract

지열발전을 위해 심부에 인공적으로 균열대를 생성시키는 EGS (Enhanced/Engineered Geothermal System) 지열발전 기술에서는 유체의 이동통로가 되는 균열의 연결성 향상이 매우 중요하며, 다단계에 걸쳐 이루어지는 수압파쇄시 발생되는 균열의 정보는 미소진동 모니터링을 통해 확인이 가능하다. 하지만 각 단계별 수압파쇄시 발생되는 균열에 의해변화된 속도구조를 고려하지 않고 미소진동 모니터링을 수행하게 되면, 다음 단계의 수압파쇄시 발생된 균열의 위치정보는 실제 위치와는 차이를 보이게 된다. 이 연구에서는 Kim et al. (2015)에 의해 개발된 미소진동 위치역산 알고리듬을심부 수 km 하부를 대상으로 하는 EGS 지열발전에 적합하도록 개선시켰으며, 각 단계별 수압파쇄시 발생되는 균열에 의해 변화된 속도구조를 측정된 미소진동 자료를 이용하여 영상화할 수 있는 3차원 속도역산 알고리듬을 개발하였다. 아이코날 방정식(Eikonal equation)을 사용하여 단순 층서구조뿐만 아니라 복잡한 속도구조의 경우에도 적용가능하도록 하였고 그림자 영역(shadow zone)에 대해서도 어느 위치에서나 정확한 주시계산이 이루어지도록 하였으며, 프레넬 볼륨(Fresnel volume)을 이용한 자코비안(Jacobian) 계산을 통하여 속도역산의 계산시간을 효과적으로 단축시켰다. 또한, EGS 사이트를 모사한 속도모델에서 얻어진 미소진동 자료를 개발된 알고리듬에 적용시킨 결과, 전 단계에 이루어진 수압파쇄에 의해 변화된 속도를 반영하는 향상된 속도모델을 얻을 수 있었고 이를 이용하여 위치 재결정을 수행한 결과 실제 위치와거의 일치하는 결과를 얻었다. In order to enhance the connectivity of fracture network as fluid path in enhanced/engineered geothermal system (EGS), the exact locating of hydraulic fractured zone is very important. Hydraulic fractures can be tracked by locating of microseismic events which are occurred during hydraulic fracture stimulation at each stage. However, since the subsurface velocity is changed due to hydraulic fracturing at each stage, in order to find out the exact location of microseismic events, we have to consider the velocity change due to hydraulic fracturing at previous stage when we perform the mapping of microseimic events at the next stage. In this study, we have modified 3D locating algorithm of microseismic data which was developed by Kim et al. (2015) and have developed 3D velocity update algorithm using occurred microseismic data. Eikonal equation which can efficiently calculate traveltime for complex velocity model at anywhere without shadow zone is used as forward engine in our inversion. Computational cost is dramatically reduced by using Fresnel volume approach to construct Jacobian matrix in velocity inversion. Through the numerical test which simulates the geothermal survey geometry, we demonstrated that the initial velocity model was updated by using microseismic data. In addition, we confirmed that relocation results of microseismic events by using updated velocity model became closer to true locations.