[논문 리뷰] Fundamentals in Peaceman Model for Well-Block radius For Non-Linear Flows Near Well

Akif Ibragimov, Э.С. Закиров외 2명|arXiv (Cornell University)|2022. 03. 18.

Heat and Mass Transfer in Porous Media인용 수 4

한 줄 요약

이 논문은 유한 차분 해법과 해석적 반경 흐름 해법을 잘 고려한 엄밀한 수학적 프레임워크를 사용하여 Peaceman의 웰블록 반경 공식을 재해석한다. 비선형 Forchheimer 흐름을 위한 일반화된 웰블록 반경을 유도하여, 흐름 속도와 Forchheimer β 계수에 따라 달라지며, 널리 사용되는 Dake 공식을 수정함으로써 비-Darcy 흐름 영역에서 정확한 오일 샌드 시뮬레이션을 가능하게 한다.

ABSTRACT

We consider sewing machinery between finite difference and analytical solutions defined at different scale: far away and near source of the perturbation of the flow. One of the essences of the approach is that coarse problem and boundary value problem in the proxy of the source model two different flows. We are proposing method to glue solution via total fluxes, which is predefined on coarse grid. It is important to mention that the coarse solution "does not see" boundary. From industrial point of view our report provide mathematical tool for analytical interpretation of simulated data for fluid flow around a well in a porous medium. It can be considered as a mathematical "shirt" on famous Peaceman well-block radius formula for linear (Darcy) radial flow but can be applied in much more general scenario. As an important case, we consider nonlinear Forchheimer flow. In the article we rigorously obtained well-block radius, explicitly depending on $β-$Forchheimer factor and total rate of the flow on the well, and provide generalization of the Dake Formula and evaluation of the $D-$factor.

연구 동기 및 목표

• 비선형 흐름, 특히 Forchheimer 유형의 흐름에 대해 Peaceman의 웰블록 반경 개념을 수학적으로 재구성하여 오일 샌드 시뮬레이션의 정확도를 향상시키는 것.
• 기존 Darcy 기반 Peaceman 반경을 비선형 흐름 영역에 적용할 때 D 계수를 일정하다고 가정하는 데서 발생하는 모순을 해결하는 것.
• 통합된 물질 보존 방정식과 해석적 반경 흐름 해법을 연결하는 엄밀한 '접합' 방법을 수립하는 것.
• Forchheimer β 계수와 유량에 직접적으로 의존하는 새로운 웰블록 반경 공식을 유도하여 선형 및 비선형 흐름 영역 모두에 적용 가능하게 하는 것.
• 오일 샌드 시뮬레이터에서 사용되는 D 계수는 보편적이지 않으며, 격자 크기 Δ과 유량 q에 따라 달라지며, 비-Darcy 조건에서 특히 중요하다는 것을 입증하는 것.

제안 방법

• 내부 영역은 비선형 Forchheimer 흐름, 외부 영역은 선형 물질 보존을 갖는 굵은 격자 영역인 이중 영역 원환형 모델을 사용하며, 유량 연속성 조건으로 연결한다.
• 내부 원환 영역의 해석적 해를 표현하기 위해 그린 함수 방법을 적용하여 경계에서 정확한 유량 일치를 달성한다.
• 유한 차분 스텐실에서 유출된 총 유량과 해석적 해에서 유출된 총 유량을 동일시함으로써 웰블록 반경 R₀를 유도하며, 이는 굵은 격자 물질 보존과 일치시킨다.
• 비선형 흐름의 맥락에서 수학적 일관성과 해의 유일성을 확보하기 위해 수정된 Peaceman의 잘 정의된 조건을 도입한다.
• Forchheimer 법칙에 따르는 비선형 반경 흐름 방정식 −∂p/∂r = αv + βv²를 해결하고, 웰보어에서 격자 블록의 외부 경계까지 통합한다.
• 표준 Darcy 기반 R₀를 대체하여, β, q, 격자 크기 Δ에 따라 달라지는 보정 항을 포함하는 새로운 표현식을 유도한다. 이는 경계에서의 유량 일치에서 유도된다.

Figure 1 . Five-spot grid of size $\Delta$ with the well at center, and auxiliary $U(0,R_{w},R_{1})$ and $U(0,R_{w},R_{0})$

실험 결과

연구 질문

• RQ1 기존 Peaceman 웰블록 반경 공식을 비선형 Forchheimer 흐름으로 일반화할 수 있으며, 유한 차분 물질 보존과 일관성을 유지할 수 있는가?
• RQ2 비선형 영역에서 웰블록 반경 R₀는 Forchheimer β 계수와 유량 q에 어떻게 의존하는가?
• RQ3 비-Darcy 흐름이 존재할 경우 오일 샌드 시뮬레이터에서 사용되는 D 계수의 정확한 형태는 무엇이며, 격자 크기 Δ과 독립적인가?
• RQ4 Forchheimer 흐름 시뮬레이션에서 Darcy 기반 R₀를 사용할 경우 체계적인 오차가 발생하는가? 만약 그렇다면 어떤 조건에서 발생하는가?
• RQ5 유량이 다스러지는 웰 주변의 다중 척도 흐름 문제에서 유한 차분 해법과 해석적 해법을 연결하는 통합된 수학적 프레임워크를 구축할 수 있는가?

주요 결과

• Forchheimer 흐름의 웰블록 반경은 더 이상 일정하지 않으며, 기존 Peaceman 공식과 달리 흐름 속도 q와 Forchheimer β 계수에 직접적으로 의존한다.
• 유도된 R₀ 공식은 q²에 비례하는 보정 항을 포함하며, Δ가 작아지거나 q가 증가할수록 이 항이 두드러지게 되어 일정한 D 계수를 가정하는 것이 타당하지 않음을 보여준다.
• 압력 강하의 정확한 표현은 p₀ − p_w = (α/2π)q ln(R₀_Darcy / R_w) + (βq² / 4π²)(1/R_w − 1/Δ)이며, 이는 격자 크기 Δ에 명시적인 의존성을 보인다.
• 기존 시뮬레이터에서 사용되는 D 계수는 보편적이지 않으며, Δ과 q에 따라 조정되어야 하며, 비-Darcy 흐름에서는 격자 및 유량에 의존한다.
• 고전적 Peaceman 반경 R₀_Darcy = Δ·e^(-π/2)는 큰 Δ 또는 낮은 q에서만 좋은 근사값을 제공하며, 미세 격자 또는 고속도에서는 보정 항이 상당히 크다.
• 이 방법은 유한 격자 블록 압력 p₀와 실제 웰보어 압력 p_w 사이의 물리적으로 일관된 유량 일치 원리를 통해 시뮬레이터에서 계산된 압력을 해석하는 데 엄밀한 수학적 기반을 제공한다.

Figure 2 . Five-point stencil for the finite difference scheme, and domain of the flow with source at $0$