Improving Water Saturation Estimation in a Low Resistivity Shaly Sand Reservoir Through Integrated Electrofacies, Rock Typing, Clay Mineral Analysis, and Saturation Height Modeling: NEAG-5 Field, Egypt

기술 요약: 저저항도 셰일 샌드 저류층에서의 수층 포화도 추정 개선 (이집트 NEAG-5 유전)

문제 정의
이집트 Abu Gharadig 분지의 Middle Bahariya Formation에 위치한 NEAG-5 유전의 경우, 저류층의 특성이 저저항도 셰일 샌드(low-resistivity shaly sand)라는 점 때문에 역사적으로 정확한 수층 포화도($S_w$) 추정이 어려웠습니다. 전도성 점토 광물의 존재와 강한 저류층 불균질성 및 복잡한 공극계가 결합되어, 기존의 Archie 기반 접근 방식은 수층 포화도를 체계적으로 과대평가하는 결과를 초도했습니다. 이러한 과대평가는 생산적인 구간을 수층 구간으로 오분류하게 만들어, 순 페이(net pay)의 과소 산정, 보수적인 탄화수소 체적 계산, 그리고 조기 수분 돌파(water breakthrough)를 포함한 최적화되지 않은 유전 개발 전략을 초래합니다.

방법론
본 연구는 점토 전도성과 저류층 불균질성의 영향을 포화도 계산으로부터 분리하기 위해 설계된 통합 물성학적 워크플로우를 제안합니다. 방법론은 다음과 같은 네 가지 주요 단계로 진행됩니다:

1. 저류층 암석 유형 결정 (Reservoir Rock Typing, RRT): 저자는 루틴 코어 분석(RCA) 데이터를 사용하여 저류층 품질 지수(RQI), 정규화된 공극률 지수(NPI), 유동 존 지표(FZI)를 계산했습니다. 이 유압 유동 단위(hydraulic flow unit) 분석을 통해 공극 목(pore throat) 구조와 유동 용량에 기초하여 네 가지 뚜렷한 저류층 암석 유형(RT-1 ~ RT-4)을 식별했습니다.
2. 전기적 상(Electrofacies) 모델링: 코어 데이터가 없는 시추공까지 암석 유형 분류를 확장하기 위해, 일반적인 와이어라인 로그(GR, RHOB, NPHI, Rt)와 코어 유래 암석 유형을 학습시킨 색인 확률적 자기 조직화 지도(Indexed Probabilistic Self-Organizing Map, IPSOM)를 사용했습니다. 이를 통해 유전 전역에 걸친 연속적인 전기적 상 예측치를 생성했습니다.
3. 광물학적 특성 규명: X선 회절법(XRD)과 Quanti Elan 역산 워크플로우를 사용하여 정량적 광물 구성을 확립했습니다. 이를 통해 특정 점토 광물(카올리나이트, 일라이트, 혼합층 일라이트-스멕타이트, 글루코나이트)의 분포와 저류층 품질 간의 상관관계를 정량화했습니다.
4. 포화도 모델링: 네 가지 독립적인 포화도 모델을 평가하고 비교했습니다:
   • 기존 Archie 모델: 코어 유래 지수($m=1.68, n=2.01$)와 형성수 비저항($R_w = 0.13 \Omega\cdot m$)으로 보정되었습니다.
   • 저저항도 페이 (Low Resistivity Pay, LRP): 중성자-밀도 및 감마선 크로스플롯 엔드포인트를 사용하여 도구 해상도 미만의 얇은 셰일 층상 구조에 대해 보정되었습니다.
   • Quanti Elan: 점토 분포 효과를 조정하기 위해 광물 부피를 포함했습니다.
   • 포화도 높이 함수 (Saturation Height Function, SHF): 네 가지 특정 저류층 암석 유형에 맞춰 보정된 Leverett J-함수 변환을 사용하여 개발되었습니다. 이 모델은 모세관 압력 데이터와 형성 압력 구배를 활용하여 암석 유형별 $S_w$ 방정식을 생성했습니다.

주요 결과

• 암석 유형 분포: 워크플로우는 뚜렷한 유동 특성을 가진 네 가지 암석 유형을 성공적으로 식별했습니다. 고품질인 RT-4는 투과율 >1000 mD 및 낮은 불변 수층 포화도를 나타냈으며, 저품질인 RT-1은 투과율 <10 mD 및 높은 결합수(bound water)를 나타냈습니다.
• 광물학적 영향: XRD 및 Quanti Elan 결과는 고품질 구간이 낮은 총 점토 함량을 가진 카올리나이트에 의해 지배되는 반면, 저품질 구간은 투과율을 저하시키고 전기 전도도를 높이는 일라이트, 혼합층 점토 및 글루코나이트의 농도가 더 높다는 것을 확인했습니다.
• 포화도 모델 성능:
  • 기존 Archie 모델은 점토가 풍부한 구역에서 $S_w$를 지속적으로 과대평가하여, 탄화수소가 존재하는 구간에서도 종종 75%에 가까운 값을 예측했습니다.
  • LRP 모델은 계산된 $S_w$를 약 50~55%로 낮추어 탄화수소 식별 능력을 개선했습니다.
  • 암석 유형 의존적 SHF 모델은 가장 지질학적으로 일관된 결과를 제공했으며, 모든 투과율 범위에서 생산 데이터와 더 잘 일치하는 모습을 보였습니다. 저품질 셰일 샌드(RT-2)에서 SHF 및 LsRP 모델은 Archie 모델 대비 계산된 수층 포화도를 약 50% 감소시켰습니다.
• 동적 검증: 수정된 특성 분석은 생산 이력과 강한 상관관계를 보였습니다. 고품질 유정인 NEAG-5-1(0.703 MMstb)은 유리한 전기적 상과 낮은 점토 함량에 의해 지배된 반면, 저조한 성능을 보인 NEAG-5-3(0.057 MMstb)은 파쇄 시도에도 불구하고 점토가 풍부한 저품질 상에 의해 지배되었습니다.

의의 및 현장 적용
본 연구는 불균질한 Middle Bahariya Formation에 대해 단일한 유전 전역 포화도 모델을 사용하는 것이 부적절함을 입증합니다. 전기적 상, 암석 유형 및 광물 분석을 통합함으로써 워크플로우는 평가 불확실성을 크게 줄였습니다.

이 연구의 주요 영향은 탄화수소 체적 추정치의 상당한 수정입니다. 통합 워크플로우를 통해 초기 저류층 내 석유량(STOIIP)이 약 5 MMstb에서 19 MMstb로 증가했습니다. 결과적으로 본 연구는 다음을 촉진했습니다:

• 이전에 과소 평가되었던 순 페이 구간의 식별.
• 유전 전역에 걸쳐 90미터 이상의 순 오일 샌드(net oil sand)를 추가한 완치(perforation) 작업 수행.
• 수분 함량(water cut)을 완화하기 위해 주입을 고품질 구간으로부터 격리함으로써 수층 주입 관리 최적화.
• 일라이트가 지배적인 저품질 상에 대해 유일하게 경제적인 개발 옵션으로서 수평정을 사용하는 전략적 권고.

결론적으로, 이 통합 프레임워크는 불균질한 저저항도 셰일 샌드 저류층에서의 매장량 추정 및 개발 계획을 위한 견고하고 반복 가능한 방법론을 제공합니다.