The multiscale restriction smoothed basis method for fractured porous media (F-MsRSB)

详细信息    查看全文

• 作者：Swej Shah^a ; ^{S.Y.Shah@tudelft.nl" class="auth_mail" title="E-mail the corresponding author} ; Olav Mø ; yner^b ; ^{Olav.Moyner@sintef.no" class="auth_mail" title="E-mail the corresponding author} ; Matei Tene^a ; ^{M.Tene@tudelft.nl" class="auth_mail" title="E-mail the corresponding author} ; Knut-Andreas Lie^b ; ^{Knut-Andreas.Lie@sintef.no" class="auth_mail" title="E-mail the corresponding author} ; Hadi Hajibeygi^a ; ^{H.Hajibeygi@tudelft.nl" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Fractured porous media ; Embedded fractured modeling ; Multiscale finite-volume methods ; Restriction smoothed basis ; Algebraic iterative multiscale solver ; Scalable linear solvers ; Multiphase flow ; Reservoir simulation
• 刊名：Journal of Computational Physics
• 出版年：2016
• 出版时间：1 August 2016
• 年：2016
• 卷：318
• 期：Complete
• 页码：36-57
• 全文大小：12566 K

文摘

A novel multiscale method for multiphase flow in heterogeneous fractured porous media is devised. The discrete fine-scale system is described using an embedded fracture modeling approach, in which the heterogeneous rock (matrix) and highly-conductive fractures are represented on independent grids. Given this fine-scale discrete system, the method first partitions the fine-scale volumetric grid representing the matrix and the lower-dimensional grids representing fractures into independent coarse grids. Then, basis functions for matrix and fractures are constructed by restricted smoothing, which gives a flexible and robust treatment of complex geometrical features and heterogeneous coefficients. From the basis functions one constructs a prolongation operator that maps between the coarse- and fine-scale systems. The resulting method allows for general coupling of matrix and fracture basis functions, giving efficient treatment of a large variety of fracture conductivities. In addition, basis functions can be adaptively updated using efficient global smoothing strategies to account for multiphase flow effects. The method is conservative and because it is described and implemented in algebraic form, it is straightforward to employ it to both rectilinear and unstructured grids. Through a series of challenging test cases for single and multiphase flow, in which synthetic and realistic fracture maps are combined with heterogeneous petrophysical matrix properties, we validate the method and conclude that it is an efficient and accurate approach for simulating flow in complex, large-scale, fractured media.