Validation of a Modified Normalized Spotted Hyena Optimization Fractal Dimension of Characterizing Pore Structure in Permo–Triassic Khuff Carbonate Reservoirs, Central Saudi Arabia

기술 요약: Permo–Triassic Khuff 탄산염암 저류층의 공극 구조 특성화를 위한 수정된 정규화된 스포티드 하이에나 최적화 프랙탈 차원 검증

문제 정의
중부 사우디아라비아의 Permo–Triassic Khuff 층과 같은 탄산염암 저류층은 복잡한 퇴적 및 속성 작용의 역사로 인해 극도로 높은 공극 네트워크 불균질성을 나타낸다. 실트질/사질 저류층과 달리, 이 저류층들은 유체 저장 및 흐름에 중요한 영향을 미치는 다중 스케일 공극 네트워크(입자 간 공극, 보그(vug), 균열)를 포함하고 있다. 프랙탈 이론이 공극 공간의 기하학적 구조를 특성화하고 모세관 압력과 투과율을 연결하는 틀로서 확립되어 있으나, 공극 네트워크의 불균질성을 기술하면서도 물리적 해석력을 유지할 수 있는 견고하고 대안적인 통계적 정식화가 필요하다. 특히, 프랙-탈 프레임워크 내에서 비모수적 통계 기술자(non-parametric statistical descriptors) 및 최적화 기반 기술의 적용 가능성에 대한 종합적인 평가와 검증이 필요하다.

방법론
본 연구는 Khuff Formation의 Khartam Member 표면 노출 샘플을 분석하여 세 가지 서로 다른 접근 방ا를 통해 공극 네트워크 불균질성을 평가하였다:

1. 전통적인 정규화된 공극 반경 모델: 이 기본 방법은 수은 침투 모세관 압력(MICP) 데이터를 활용하였다. 워시번(Washburn) 방정식을 사용하여 모세관 압력($P_c$)으로부터 공극 목 반경($R$)을 도출하였다. 정규화된 공극 반경($R/R_{max}$)과 수분 포화도($S_w$) 사이의 로그 관계를 설정하였다. 데이터에 2차 다항식 회귀($Y = aX^2 + bX + c$)를 적합시켰으며, 여기서 프랙탈 차원은 $D_f = 2 + b$로 계산되었다.
2. 수정된 정규화된 크루스칼-월리스(MNKW): 공극 목 분포에서 도출된 크루스칼-월리스 파라미터를 정규화하여 독립적인 통계적 기술자를 개발하였다. 전통적인 방법과 유사하게, $\log(MNKW)$와 $\log_{10}(S_w)$ 사이의 회귀 분석을 수행하여 계수 $b$를 추출하고 $D_f$를 계산하였다.
3. 수정된 정규화된 스포티드 하이에나 최적화(MNSHO): 이 최적화 기반 접근 방식은 공극 목 크기 확률 분포로부터 계산된 섀넌 엔트로피(Shannon entropy)를 활용하였다. 엔트로피 값을 정규화하여 MNSHO 파라미터를 생성하였다. $\log_{10}(S_w)$에 대한 $\log(MNSHO)$의 이차 회귀를 적합시켜 프랙탈 차원을 결정하였다.

모든 샘플은 헬륨 공극률 및 정상 상태 가스 투과도 측정을 포함한 표준 암석물리 테스트를 거쳤다. 본 연구는 프랙탈 차원과 공극률 및 투과율 간의 상관관계를 규명하기 위해 회귀 분석을 사용하였으며, MNKW와 MNSHO 방법 간의 일치도를 평가하기 위해 Bland–Altman 분석을 사용하였다.

주요 결과

• 프랙탈 차원 범위: 계산된 프랙탈 차원은 약 2.36에서 2.86 사이였다. 낮은 값은 연결성이 제한된 균질한 공극 시스템을 나타냈고, 높은 값은 개선된 유체 흐름 경로를 가진 복잡하고 불균질한 네트워크를 나타냈다.
• 암석물리적 상관관계: 공극률은 2.269%에서 11.253% 범위를 나타냈고, 투과도는 0.264에서 3.445 mD 범위를 나타냈다. 프랙탈 차원과 투과율 사이에 강한 양의 상관관계가 관찰되었다. 높은 프랙탈 차원을 가진 샘플은 일관되게 더 나은 투과도와 공극 네트워크 연결성을 보였으며, 이는 프랙탈 기하학이 공극률 단독보다 흐름 경로의 연속성을 더 효과적으로 포착함을 시사한다.
• 방법론적 검증:
  • 회귀 분석: MNKW 접근 방식은 전통적인 $\log(R/R_{max})$ 방법과 거의 완벽한 일치($R^2 = 0.9999986$)를 보였다. MNSHO 접근 방식 또한 전통적인 방법과 매우 강력한 양의 상관관계($R^2 \approx 0.9999$)를 입증하였다.
  • Bland–Altman 분석: 이 분석은 MNKW와 MNSHO 방법 간의 우수한 일치도를 확인하였으며, MNSHO 접근 방식에서 작은 체계적 편향이 관찰되었다. 좁은 일치 한계(limits of agreement)는 두 방법 모두 신뢰할 수 있고 재현 가능한 공극 네트워크 불균질성 추정치를 제공함을 나타낸다.

핵심 기여

• MNKW 및 MNSHO 도입: 본 연구는 탄산염암 저류층을 위한 견고한 프랙탈 기술자로서 수정된 정규화된 크루스칼-월리스(MNKW) 및 수정된 정규화된 스포티드 하이에나 최적화(MNSHO) 파라미터를 성공적으로 도입하고 검증하였다.
• 최적화 기반 기술의 검증: 본 연구는 최적화 기반 기술(MNSHO)이 모세관 압력 데이터로부터 도출된 전통적인 프랙탈 차원 추정치를 정확하게 재현할 수 있음을 입증하였다.
• 불균질성의 정량화: 본 연구는 프랙탈 차원이 공극 크기 분포, 연결성 및 구조적 복잡성을 통합하는 종합적인 정량적 기술자로서, 공극률 측정보다 더 미묘한 저류층 품질의 관점을 제공함을 확립하였다.

의의 및 주장
본 논문은 프랙탈 차원이 탄산염암 공극 시스템을 특성화하는 견고한 정량적 기술자이며, 불균질한 탄산염암 저류층의 특성을 효과적으로 설명한다고 주장한다. 결과적으로, 통계적(MNKW) 및 최적화 기반(MNSHO) 접근 방식 모두가 모세관 압력 데이터로부터 도출된 전통적인 방법의 유효한 대안임을 입증하였다.

저자들은 이러한 방법론이 다음 분야에 적용될 잠재력이 있다고 주장한다:

• 저류층 특성화 및 디지털 암석 분석.
• 불균질한 탄산염암 형성층의 투과도 예측.
• 공극 연결성이 유체 흐름의 핵심 제어 요소인 고도화된 저류층 품질 평가.

본 연구는 제안된 프레임워크가 Khuff Formation 및 기타 불균질한 탄산염 시스템을 평가하기 위한 신뢰할 수 있는 도구를 제공하며, 프랙탈 차원의 증가가 향상된 공극 네트워크 연결성과 복잡성에 의해 구동되는 저류층 성능 개선의 직접적인 지표임을 확인하며 결론을 맺는다.