[논문 리뷰] Pore Network and Medial Axis simultaneous extraction through Maximal Ball Algorithm

Mariane Barsi-Andreeta, Everton Lucas‐Oliveira외 3명|arXiv (Cornell University)|2019. 12. 10.

Enhanced Oil Recovery Techniques참고 문헌 18인용 수 5

한 줄 요약

이 논문은 X-ray 미세-CT 데이터로부터 극대 구의 알고리즘(MBA)을 개선하여 기공-목장망과 중추축 스켈레톤을 동시에 추출하는 새로운 방법을 제시한다. 극대 구를 가중치가 부여된 그래프로 모델링하고 다익스트라 알고리즘을 적용하여 최단 경로를 찾음으로써, 기하학적 위상 구조를 유지하는 중추축 네트워크를 생성한다. 이는 모래석 및 석회석 암석에서의 밀도 예측을 향상시키며, 하겐-포아젤 유량 모델 기반 가중치가 가장 정확한 단상류 유동 시뮬레이션을 제공한다.

ABSTRACT

Network Extraction algorithms from X-ray microcomputed tomography have become a routine method to obtain pore connectivity and pore morphology information from porous media. The main approaches for this extraction are either based on Max Ball Algorithm or Medial Axis Extraction. The first is a robust method to separate the pore space into pore and throats, which provides the pore and throat sizes distributions directly. The second simplifies the medium by thinning the volume into a one voxel wide centerline that preserves the volume's topological information. Since each method has drastically different implementations, it is not usual to use both to characterize the porous structure. This work presents a method to extract a simplified centerline of porous materials by adding few more steps to the well-established Max Ball algorithm. Results for sandstones and carbonate rock samples show that this Medial-Axis network can be used for single phase flow simulations, as well as preserves the mediums morphology.

연구 동기 및 목표

• 다공성 매체의 X-ray 미세-CT 데이터로부터 기공-목장망과 중추축 스켈레톤을 통합적으로 추출하는 방법을 개발하는 것.
• 복잡한 암석 구조에서 표면 거칠기에 민감한 기존 중추축 알고리즘의 한계를 해결하는 것.
• 기공 네트워크 모델에서 목장의 길이 추정을 정교화하여 단상류 유동 시뮬레이션의 정확도를 향상시키는 것.
• 지역 두께 측정값과 비교하여 추출된 중추축 네트워크의 형태적 충실도를 평가하는 것.
• 그래프 기반 경로 탐색 과정에서의 다양한 가중치 설정이 밀도 예측 정확도에 미치는 영향을 조사하는 것.

제안 방법

• X-ray 미세-CT 데이터로부터 기공 영역을 식별하고, 포함된 극대 구의 완전한 네트워크를 계산한다.
• 지역 최대값 구를 기공 중심으로 식별하고, 기공 간 연결 영역을 목장 영역으로 정의한다.
• 극대 구를 노드로, 연결성을 나타내는 간선을 갖는 가중치가 부여된 그래프를 구성하며, 가중치는 구 면적 또는 하겐-포아젤 유량 모델에 기반한다.
• 다익스트라 알고리즘을 적용하여 기공 중심 간 최단 경로를 계산하고, 1바voxel 두께의 중추축 네트워크를 형성한다.
• 기존 MBA 프레임워크를 사용하여 동일한 극대 구 구조에서 기공-목장망을 추출한다.
• 기하학적 보정 계수를 적용하여 기공 네트워크의 목장 길이를 보정함으로써 밀도 추정 정확도를 향상시킨다.

실험 결과

연구 질문

• RQ1 최대 구 알고리즘의 확장 기반으로 다공성 매체의 중추축을 기공-목장망과 동시에 추출할 수 있는가?
• RQ2 그래프 기반 경로 탐색 과정에서의 다양한 가중치 설정이 단상류 유동 시뮬레이션 정확도에 어떤 영향을 미치는가?
• RQ3 중추축 네트워크는 직접 두께 측정값과 비교할 때 원래 다공성 매체의 형태를 어느 정도 충실하게 유지하는가?
• RQ4 중추축에서 유도된 보정된 기공 네트워크는 원래 MBA 네트워크와 비교해 목장 길이 및 밀도 예측 측면에서 어떻게 다른가?
• RQ5 다양한 다공성 및 기공 크기 분포를 가진 실제 암석 샘플에서 제안된 방법의 계산 성능과 메모리 사용량은 어떠한가?

주요 결과

• 가중치가 부여된 극대 구 그래프 위에서 다익스트라 알고리즘을 적용하여 추출한 중추축 네트워크는 원래 다공성 매체의 위상 및 형태적 특성을 유지한다.
• 하겐-포아젤 기반 가중치 설정이 가장 정확한 밀도 예측을 제공하였으며, 테스트된 모든 샘플에서 보고된 값과 10% 이내의 예측 정확도를 보였다.
• 목장 길이를 중추축 기반으로 보정한 기공 네트워크(PNC) 모델은 원래 MBA 네트워크보다 밀도 과소평가를 최대 2배까지 감소시켰다.
• 벤스하이머 샘플의 경우, 보정된 PNC 네트워크가 2595 mD(kx)의 밀도를 예측하였고, 중추축 결과인 2916 mD와 유사하여 높은 일致성을 보였다.
• 샘플의 복잡성과 시드 수에 따라 처리 시간은 4~60분으로 변동하였으며, 가장 복잡한 샘플(Berea)의 경우 최고 50 GB의 메모리 사용량을 기록하였다.
• 중추축 네트워크의 지름 변화가 지역 두께 측정값과 양호한 일치를 보여 형태적 충실도를 확인하였다.